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

    埃及伊蚊不同組織的基因共表達模式分析

    2020-12-11 08:05:00陳曉潔
    昆蟲學報 2020年10期
    關鍵詞:伊蚊共表達觸角

    劉 杰,李 勃,陳曉潔,陳 斌

    (重慶師范大學昆蟲與分子生物學研究所,媒介昆蟲重慶市重點實驗室,重慶 401331)

    截至2017年,全球已知蚊蟲約3 573種,其中我國約有419種(付文博和陳斌,2018)。它們分布廣,繁殖快。絕大多數(shù)種類的雌蚊需要吸食血液作為補充營養(yǎng)以使卵正常發(fā)育。因其吸血行為,蚊蟲傳播許多危險的人類疾病的病原體,是多種惡性傳染病的傳播媒介(Tolle,2009)。蚊媒病是人類健康的主要威脅,每年全球有7億人受感染,造成70多萬人死亡。近年來,登革熱發(fā)病率大幅度增長,全球約有一半人口面臨登革熱的威脅(WHO,2019)(https:∥www.who.int/zh/news-room/fact-sheets/detail/dengue-and-severe-dengue)。埃及伊蚊Aedesaegypti(Linnaeus,1762)是登革熱病毒的主要傳播媒介,同時也攜帶基孔肯雅熱、黃熱病以及寨卡病毒病等傳染病病原。這些疾病每年在全世界造成數(shù)億人感染,導致數(shù)萬人死亡,并呈逐年上升趨勢(Akbarietal.,2013)。媒介蚊蟲的控制是蚊媒病控制的關鍵,已成為當今社會所面臨的亟待解決的重大問題之一。蚊蟲的不同組織承擔不同的生理功能,掌握它們的基因表達差異有利于認知它們的生物學功能。蚊蟲的遺傳控制是蚊蟲潛在的重要控制措施之一,弄清蚊蟲不同組織在行使功能過程中的共表達基因以及起重要作用的核心(hub)基因有利于尋找遺傳控制的靶標基因。

    轉(zhuǎn)錄組是特定細胞或組織在特定時間或狀態(tài)下所有轉(zhuǎn)錄產(chǎn)物的集合(Wangetal.,2010),通過對轉(zhuǎn)錄組的研究可以揭示蚊蟲組織分化過程中的基因表達模式和結構變異,并發(fā)現(xiàn)新基因等。近年來,隨著microarray和RNA-Seq等高通量技術的飛速發(fā)展與成本的降低,產(chǎn)生了大規(guī)模蚊蟲組織轉(zhuǎn)錄組學數(shù)據(jù)(Martínez-Barnetcheetal.,2012;Matthewsetal.,2016;Chenetal.,2017)。傳統(tǒng)的轉(zhuǎn)錄組學研究主要是在兩類樣本相互比較的基礎上,鑒定差異表達基因,進而對差異基因進行注釋與富集分析,并挖掘這些基因的功能表現(xiàn)和參與的通路(pathway),從而解析相關的生物學問題或分子機制。但這樣的研究策略主要局限于對組間差異表達的基因進行研究,容易忽略掉基因之間的協(xié)同表達或表達互作。目前,基于全局的權重基因共表達網(wǎng)絡分析(weighted gene co-expression network analysis,WGCNA)正好為這個問題的解決帶來了契機。WGCNA分析基于大樣本轉(zhuǎn)錄組數(shù)據(jù),根據(jù)基因表達水平的不同將基因聚類,具有相似表達模式的基因會被劃分到同一個模塊,同一模塊中的基因常具有相似的生物學功能(楊宇昕等,2019)。作為一種系統(tǒng)生物學分析方法,WGCNA可以用來鑒定高度協(xié)同變化的基因集、挖掘與目標組織或性狀相關的特異性模塊,并可進一步研究模塊內(nèi)基因的互作關系,以及篩選hub基因(Langfelder and Horvath,2008),且在描述解析復雜表型下的分子作用機制和網(wǎng)絡關系等方面具有獨特優(yōu)勢。自WGCNA方法和軟件包開發(fā)以來,它已廣泛應用于多個研究領域,并取得了一系列的研究成果(Greenhametal.,2017;Liuetal.,2018;Dengetal.,2019)。

    在蚊蟲轉(zhuǎn)錄組研究領域,Hickner等(2015)用WGCNA方法研究了尖音庫蚊Culexpipiens滯育誘導的潛在調(diào)控機制,構建了包含5個模塊的基因共表達網(wǎng)絡,通過對網(wǎng)絡模塊中基因的富集分析,鑒定出了與滯育誘導相關的多個基因集;通過對埃及伊蚊環(huán)境驅(qū)動變異介導的模塊化基因網(wǎng)絡研究,Kang等(2018)揭示了先天免疫反應的關鍵轉(zhuǎn)錄調(diào)控因子和關鍵基因。這些研究表明,將轉(zhuǎn)錄組學數(shù)據(jù)與WGCNA方法相結合,可為探索蚊蟲生物學過程相關分子機制提供強有力的工具。然而,尚未有基于WGCNA方法對蚊蟲不同組織基因共表達模式的研究。

    蚊蟲的不同組織在其傳播疾病過程中起著重要的作用,可作為防控的靶標。本研究利用WGCNA方法,結合RNA-Seq轉(zhuǎn)錄組數(shù)據(jù),構建了埃及伊蚊9種組織加權基因共表達網(wǎng)絡,挖掘出與組織特異性相關的共表達模塊,研究模塊基因功能,并篩選出在組織特異性模塊內(nèi)起重要作用的hub基因,揭示基因間的相互作用關系,旨在為探究蚊蟲不同組織特有的基因資源信息,進而解析基因功能與組織功能的聯(lián)系,篩選遺傳控制的靶標基因提供新的思路和途徑。

    1 材料與方法

    1.1 轉(zhuǎn)錄組數(shù)據(jù)來源

    埃及伊蚊不同組織的轉(zhuǎn)錄組數(shù)據(jù)來源于NCBI(National Center for Biotechnology Information)的SRA(Sequence Read Archive)數(shù)據(jù)庫,獲取碼為PRJNA236239(Matthewsetal.,2016)。選擇其中未吸血狀態(tài)下雌蚊和雄蚊9種組織的雙端測序數(shù)據(jù)進行分析,包括雌蚊的觸角、腦、喙、下顎須、卵巢,以及雄蚊的觸角、腦、前足、中足、后足和腹部末端(被定義為雄蚊包括生殖器在內(nèi)的3個末端腹節(jié))等。不同組織數(shù)據(jù)包含的生物學重復次數(shù)為:雌蚊喙5次生物學重復,觸角、下顎須4次重復,卵巢3次重復,腦2次重復,雄蚊各組織均3次重復,總計36個樣本。WGCNA分析根據(jù)基因表達水平的不同將基因聚類。在轉(zhuǎn)錄組研究中,可通過比對到基因組區(qū)域或基因外顯子區(qū)的測序序列(reads)計數(shù)(read counts)來估計基因的表達水平。但reads計數(shù)除了與基因的真實表達水平成正比外,還與基因的長度和測序深度成正相關,因此使得不同基因的表達水平不具有可比性。為了解決這個問題,人們引入了FPKM(fragments per kilobase million)、RPKM(reads per kilobase million)和TPM(transcripts per million)的概念。其中TPM是每百萬條reads的轉(zhuǎn)錄本,其同時考慮了測序深度和基因長度對reads計數(shù)的影響,是目前常用的衡量基因表達水平的方法之一。本研究中埃及伊蚊不同組織基因的TPM值采用Matthews等(2016)文獻中Additional File 6中的結果。

    1.2 基因共表達網(wǎng)絡構建

    采用TPM值作為標準化的基因表達水平定量,本研究共納入36個樣本。在應用WGCNA前,基于公式log2(TPM+1)對數(shù)據(jù)進行轉(zhuǎn)換(趙志洪,2017)。經(jīng)過缺失值移除以及方差計算后,篩選出方差最大的5 000個基因,使用R(v.3.6.1)軟件中(https:∥cran.r-project.org/bin/windows/base/old/)的WGCNA(v.1.67)包來構建加權基因共表達網(wǎng)絡并劃分模塊。首先利用函數(shù)pickSoftThreshold(Langfelder and Horvath,2008)計算權重值,通過選擇合適的權重值,使網(wǎng)絡符合無尺度網(wǎng)絡分布。

    根據(jù)上一步分析結果,選擇最優(yōu)的軟閾值。使用自動網(wǎng)絡構建函數(shù)blockwiseModules(Langfelder and Horvath,2008)構建共表達網(wǎng)絡,除參數(shù)minModuleSize=50外,其余參數(shù)按照默認設置。

    1.3 基因表達組織特異性模塊鑒定

    一組表達模式高度相似的基因常聚集在同一模塊,而具有相似表達模式的基因被認為具有類似的生物學功能(楊宇昕等,2019)。一般而言,如果某模塊與樣品的某種特定表型的相關性顯著高于其他模塊,說明這一個模塊可能與該表型存在強的關聯(lián),為特異性模塊。本研究由于只有埃及伊蚊不同組織的轉(zhuǎn)錄組數(shù)據(jù),并沒有其他表型數(shù)據(jù),在此直接將不同組織作為研究的表型性狀,以分析模塊與組織的聯(lián)系。首先建立一個行為36個組織樣本,列為11個不同取樣組織名稱的矩陣,若行與列為同一組織,即為1,不是同一組織即為0。計算每個模塊特征向量基因(module eigengene,ME)和這個矩陣之間的相關系數(shù),該數(shù)值越接近1,表示模塊與樣品表型的正相關性越強;越接近-1,表示模塊與樣品表型負相關性越強。本研究中模塊與組織高度正相關的篩選標準為:模塊與組織相關系數(shù)大于0.65,且統(tǒng)計顯著性P<0.05。選擇與組織高度正相關的模塊,作為組織特異性模塊進行后續(xù)分析(Downsetal.,2013)。

    1.4 組織特異性模塊基因的GO和KEGG富集分析

    首先,使用Bioconductor中的AnnotationHub(v.2.18.0)包獲取埃及伊蚊的注釋信息,下載埃及伊蚊OrgDb數(shù)據(jù)庫(AH76482,v.2.1)中的注釋數(shù)據(jù),再利用Bioconductor中的clusterProfiler(v.3.14.3)包對模塊內(nèi)基因進行GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)富集分析(Yuetal.,2012)。一般情況下,P<0.01,q<0.05,則為富集功能項或通路。

    1.5 模塊基因互作網(wǎng)絡構建篩選hub基因

    導出WGCNA共表達網(wǎng)絡分析結果到Cytoscape(v.3.7.1)軟件(Suetal.,2014)中,對基因調(diào)控網(wǎng)絡進行可視化。通過CytoHubba插件的12種算法,得到連接度在前10的基因,選擇在不同算法中重復出現(xiàn)次數(shù)大于5的基因,在本研究中被定義為hub基因。

    2 結果

    2.1 權重基因共表達網(wǎng)絡構建

    根據(jù)基因表達水平對樣本進行聚類。如圖1所示,基于平均(average)距離和層次聚類算法聚類結果表明,基因表達模式除了在雄成蚊前足和中足中出現(xiàn)了一定的散亂之外,測序所用的埃及伊蚊成蚊組織大多具有很高的重復性。雌/雄成蚊的觸角、腦中基因表達模式相近,聚類在一起;雄成蚊前足、中足及后足中基因表達模式相近,聚類在一起;雌成蚊喙與下顎須中基因表達模式相近,聚類在一起。對軟閾值(β)進行篩選,最終選取β值為12。通過WGCNA方法成功構建埃及伊蚊成蚊不同組織加權基因共表達網(wǎng)絡,并將網(wǎng)絡劃分為11個模塊(圖2:A)。圖2(A)上半部分的聚類樹中每個枝代表一個模塊,每個枝上的葉代表一個基因,圖2(A)下半部分的模塊顏色中每個顏色代表一個模塊。共表達模塊中基因的具體數(shù)量如圖2(B)所示。模塊中基因個數(shù)123~1 423不等,共含有5 000個基因。其中turquoise模塊包含的基因數(shù)目最多,達1 423;greenyellow模塊內(nèi)的基因數(shù)量最少,為123;grey表示未劃分到模塊中的基因,圖2(A)中grey模塊基本不可見,因未劃分到模塊的基因僅有7個。

    圖1 基于基因表達水平的埃及伊蚊成蚊不同組織聚類Fig.1 Clustering of different adult tissues of Aedes aegypti based on gene expression levelsFe_Ov:雌蚊卵巢Female ovary;Ma_Br:雄蚊腦Male brain;Fe_Br:雌蚊腦Female brain;Ma_An:雄蚊觸角Male antenna;Fe_An:雌蚊觸角Female antenna;Ma_At:雄蚊腹部末端Male abdominal tip;Ma_HL:雄蚊后足Male hindleg;Ma_FL:雄蚊前足Male foreleg;Ma_ML:雄蚊中足Male midleg;Fe_Os:雌蚊喙Female proboscis;Fe_Pa:雌蚊下顎須Female maxillary palp。根據(jù)基因表達水平,采用平均距離和層次聚類算法對樣本進行聚類。橫軸表示埃及伊蚊成蚊不同組織,代號末位的數(shù)字表示該組織的重復樣本;縱軸代表基因間的聚類高度。According to gene expression level,average distance and hierarchical clustering algorithm were used to cluster the samples.The abscissa represents the different adult tissues of A.aegypti,the last digit of the code represents a duplicate sample of the tissue,and the ordinate represents the cluster height of genes.

    圖2 埃及伊蚊成蚊不同組織基因共表達模塊Fig.2 Gene co-expression modules in different adult tissues of Aedes aegyptiA:基因聚類數(shù)和模塊分割。上半部分表示基因聚類所得聚類樹,縱坐標代表各基因間的聚類距離;下半部分表示按樹的分枝切割得到的模塊,相同模塊用同一顏色表示。Clustering dendrogram of genes and module splitting.The upper part represents dendrogram tree obtained by gene hierarchical clustering,and the ordinate represents the clustering distance between genes.The lower half represents the module cut by the branches of the tree,and the same module is represented by the same color.B:共表達模塊中基因數(shù)量分布。橫坐標代表模塊,縱坐標代表模塊中的基因數(shù)目,不同的顏色與模塊名稱對應,每個柱子上的數(shù)字表示該模塊內(nèi)的基因數(shù)目。Distribution of gene number.The abscissa represents the module,and the ordinate represents the number of genes in the module.Different colors correspond to the module name.The number on each column represents the number of genes in the module.

    2.2 組織特異性模塊分析結果

    通過篩選發(fā)現(xiàn)在11個模塊中有6個為高度的組織特異性模塊(r>0.65,P<0.05)(圖3)。在雌成蚊觸角、喙、卵巢、下顎須以及雄成蚊腦和腹部末端組織中各鑒定出1個特異性模塊。green模塊與雌成蚊觸角存在高度相關性(r=0.69,P=4E-06),purple模塊與雌成蚊喙存在高度相關性(r=0.8,P=4E-09),turquoise模塊與雌成蚊卵巢存在高度正相關性(r=0.96,P=4E-20),red模塊與雌成蚊下顎須存在高度正相關性(r=0.69,P=4E-06)。與雄成蚊腦高度正相關的模塊為blue模塊(r=0.7,P=2E-06),與雄成蚊腹部末端高度正相關的模塊為magenta模塊(r=0.96,P=2E-37)。6個組織特異性模塊中沒有任何一個模塊和雄成蚊前足、中足、后足存在較高相關性。推測可能是這3個組織中的基因主要進行基本的生命活動,因此沒有特異性的模塊與之對應。

    圖3 埃及伊蚊成蚊不同組織基因共表達模塊與性狀關聯(lián)熱圖Fig.3 Heat map of module-trait relationship of different adult tissues in Aedes aegyptiFe_An:雌蚊觸角Female antenna;Fe_Br:雌蚊腦Female brain;Fe_Os:雌蚊喙Female proboscis;Fe_Ov:雌蚊卵巢Female ovary;Fe_Pa:雌蚊下顎須Female maxillary palp;Ma_An:雄蚊觸角Male antenna;Ma_At:雄蚊腹部末端Male abdominal tip;Ma_Br:雄蚊腦Male brain;Ma_FL:雄蚊前足Male foreleg;Ma_ML:雄蚊中足Male midleg;Ma_HL:雄蚊后足Male hindleg.MEgreenyellow:黃綠色模塊特征向量基因Module eigengene of green yellow module;MEturquoise:藍綠色模塊特征向量基因Module eigengene of turquoise module;MEmagenta:品紅色模塊特征向量基因Module eigengene of magenta module;MEpurple:紫色模塊特征向量基因Module eigengene of purple module;MEred:紅色模塊特征向量基因Module eigengene of red module;MEblack:黑色模塊特征向量基因Module eigengene of black module;MEbrown:棕色模塊特征向量基因Module eigengene of brown module;MEyellow:黃色模塊特征向量基因Module eigengene of yellow module;MEblue:藍色模塊特征向量基因Module eigengene of blue module;MEgreen:綠色模塊特征向量基因Module eigengene of green module;MEpink:粉紅色模塊特征向量基因Module eigengene of pink module;MEgrey:灰色模塊特征向量基因Module eigengene of grey module.該圖基于R語言所得,橫軸表示埃及伊蚊不同組織,縱軸表示每一個模塊的特征向量。括號內(nèi)數(shù)值表示模塊與組織相關系數(shù)和統(tǒng)計顯著性。相關系數(shù)越接近1,表示模塊與樣品正相關性越強;越接近-1,表示模塊與樣品負相關性越強。Based on the R language,the horizontal axis represents the different tissues of A.aegypti,and the vertical axis represents the module eigengene of each module.The values in parentheses represent the correlation coefficient and statistical significance between the module and the tissues.The closer the correlation coefficient is to 1,the stronger the positive correlation between the module and the sample;the closer to -1,the stronger the negative correlation between the module and the sample.

    2.3 組織特異性模塊GO和KEGG富集分析

    組織特異性模塊可能參與組織特有的生物學過程,為進一步解析組織特異性模塊內(nèi)基因的生物學功能,并與相應組織功能聯(lián)系,同時也驗證共表達網(wǎng)絡構建和組織特異性模塊的可靠性,對這6個特異性模塊進行GO和KEGG富集分析。GO富集分析結果表明(表1),與雌成蚊觸角高度正相關的green模塊在纖毛組裝、基于微管運動和細胞對刺激的反應等生物學過程中發(fā)揮主要作用,其主要的細胞組分為微管相關復合物、細胞外圍、質(zhì)膜等,具有氣味結合、嗅覺受體活性、跨膜信號受體活性等生物學功能;與雄成蚊腦高度正相關的blue模塊在生物學過程調(diào)節(jié)、信號轉(zhuǎn)導、神經(jīng)系統(tǒng)過程等生物學過程中發(fā)揮主要作用,其主要的細胞組分為質(zhì)膜、突觸、突觸后膜等,具有轉(zhuǎn)運活性、神經(jīng)遞質(zhì)受體活性、鐵離子結合等功能;與雌成蚊喙高度正相關的purple模塊具有絲氨酸型肽鏈內(nèi)切酶活性、絲氨酸水解酶活性等功能;與雌成蚊卵巢高度正相關的turquoise模塊在DNA復制、細胞分裂、細胞周期等生物學過程中發(fā)揮主要作用,其主要的細胞組分為核小體、染色體等,具有DNA結合、鋅離子結合、微管結合、解旋酶活性等功能;與雌蚊下顎須高度正相關的red模塊在有機酸代謝、羧酸代謝、細胞氨基酸代謝等生物學過程中發(fā)揮作用,其細胞組分為細胞外區(qū)域,具有絲氨酸型肽鏈內(nèi)切酶活性、碳水化合物結合等功能;與雄成蚊腹部末端高度正相關的magenta模塊在羧酸代謝、草酸代謝、有機酸代謝等生物學過程中發(fā)揮主要作用。

    表1 埃及伊蚊成蚊組織特異性基因共表達模塊的GO富集分析Table 1 GO enrichment analysis of tissue specific gene co-expression module in Aedes aegypti adults

    KEGG富集分析結果表明(表2),與雌成蚊觸角高度正相關的green模塊顯著富集在藥物代謝-細胞色素P450、不飽和脂肪酸的生物合成以及谷胱甘肽代謝等通路;與雄成蚊腦高度正相關的blue模塊顯著富集在刺激神經(jīng)組織中的交互、光傳導以及其他多糖降解通路;與雌成蚊卵巢高度正相關的turquoise模塊顯著富集在DNA復制、錯配修復、戊糖和葡萄糖醛酸鹽的相互轉(zhuǎn)化等通路;與雌成蚊下顎須高度正相關的red模塊顯著富集在甘氨酸、絲氨酸和蘇氨酸的代謝通路以及精氨酸和脯氨酸代謝通路;magenta和purple模塊未發(fā)現(xiàn)顯著富集通路。

    表2 埃及伊蚊成蚊組織特異性基因共表達模塊的KEGG富集分析Table 2 KEGG enrichment analysis of tissue specific gene co-expression module in Aedes aegypti adults

    2.4 組織特異性模塊hub基因識別

    整合CytoHubba插件12種不同算法的結果,取重復出現(xiàn)次數(shù)大于5次的基因,共得到47個hub基因(green,blue,purple,turquoise,red和magenta模塊分別有9,5,6,9,9和9個)(表3)。green,blue,purple,turquoise,red和magenta模塊中連通性最高的基因分別為AAEL010426,AAEL002896,AAEL002600,AAEL000961,AAEL007784和AAEL006429。

    表3 埃及伊蚊成蚊組織特異性基因共表達模塊內(nèi)的hub基因Table 3 Hub genes within a tissue-specific module in Aedes aegypti adults

    基于WGCNA共表達網(wǎng)絡分析結果,對組織特異性模塊內(nèi)連通性最高的核心基因進行可視化展示,因與核心基因互作的基因較多,僅對權重值前50個有關聯(lián)的基因進行可視化(圖4)。結果發(fā)現(xiàn),green模塊的基因在雌/雄成蚊觸角中高表達,與核心基因AAEL010426具有較高互作網(wǎng)絡關系的基因為AAEL017505,AAEL017129,AAEL017104,AAEL007147和AAEL013563;blue模塊的基因在雌/雄成蚊腦中高表達,與核心基因AAEL002896具有較高互作網(wǎng)絡關系的基因為AAEL012541,AAEL004246,AAEL007559,AAEL003472和AAEL002220;purple模塊的基因在雌成蚊喙中表達最高,與核心基因AAEL002600具有較高互作網(wǎng)絡關系的基因為AAEL014382,AAEL017265,AAEL011610,AAEL013032和AAEL014390;turquoise模塊的基因在雌成蚊卵巢中高表達,與核心基因AAEL000961具有較高互作網(wǎng)絡關系的基因為AAEL000898,AAEL013338,AAEL000923,AAEL007879和AAEL011516;red模塊的基因在雌成蚊下顎須、雄成蚊腹部末端較高表達,與核心基因AAEL007784具有較高互作網(wǎng)絡關系的基因為AAEL000044,AAEL003114,AAEL007033,AAEL002624和AAEL012726;magenta模塊的基因在雄成蚊腹部末端高表達,與核心基因AAEL006429具有較高互作網(wǎng)絡關系的基因為AAEL010040,AAEL006414,AAEL002000,AAEL001550和AAEL005107。

    圖4 埃及伊蚊成蚊組織特異性基因共表達模塊內(nèi)連通性最高的核心基因的基因網(wǎng)絡Fig.4 Gene networks of the most interconnected hub gene in the tissue-specific gene co-expression module of Aedes aegypti adults圖中綠色、藍色、紫色、藍綠色、紅色和品紅色線條分別表示green,blue,purple,turquoise,red和magenta模塊。每個網(wǎng)絡圖中位于中心的為該組織特異性基因共表達模塊內(nèi)連通性最高的核心基因,其余帶有顏色標記的基因表示與模塊內(nèi)連通性最高的核心基因具有較高互作網(wǎng)絡關系的前5個基因,顏色越深表示互作關系越強。網(wǎng)絡中只列出與核心基因互作權重值較高的基因。The green,blue,purple,turquoise,red,and magenta lines in the figure represent the green,blue,purple,turquoise,red,and magenta modules,respectively.At the center of each network diagram is the hub gene with the highest connectivity within the tissue-specific gene co-expression module.The remaining genes with color marks represent the top 5 genes with a higher interactive network relationship with the hub genes with the highest connectivity in the module.The darker the color,the stronger the interactive relationship.In the network,only the genes with high weights interacted with the hub genes were listed.

    已知AAEL007373和AAEL008046可能參與埃及伊蚊尋找宿主過程中氨的檢測(Matthewsetal.,2016)。基于WGCNA分析結果,發(fā)現(xiàn)AAEL007373被劃分到green模塊,該模塊基因在雌/雄成蚊觸角中高表達;AAEL008046被劃分到blue模塊,該模塊基因在雌/雄成蚊腦中高表達。以這兩個基因為核心構建局部網(wǎng)絡(圖5),發(fā)現(xiàn)green模塊內(nèi)的基因AAEL007373與離子通道型受體基因AAEL000047、氣味受體基因AAEL015147和AAEL005999等具有較高互作網(wǎng)絡關系;green模塊內(nèi)的核心基因AAEL014089,AAEL007669,AAEL012424,AAEL010426,AAEL001510和AAEL013197也均與其存在互作。blue模塊內(nèi)的基因AAEL008046與AAEL004246,AAEL006825和AAEL000504等存在關聯(lián);blue模塊內(nèi)的核心基因AAEL002896,AAEL010918和AAEL001696也與其存在互作網(wǎng)絡關系。

    圖5 埃及伊蚊AAEL007373和AAEL008046的基因網(wǎng)絡Fig.5 Gene networks of AAEL007373 and AAEL008046 of Aedes aegypti圖中綠色和藍色分別表示green和blue模塊。網(wǎng)絡中只列出與核心基因互作權重值較高的基因。Green and blue colors in the figure represent green and blue modules,respectively.In the network,only the genes with high weights interacted with the hub genes were listed.

    3 討論與結論

    基因的表達存在相關性,功能相關的一組基因常協(xié)調(diào)一致,共同表達。對蚊蟲不同組織的基因共表達模式進行研究,可獲得組織特有的基因資源信息,得到組織在行使功能過程中的共表達基因以及起重要作用的核心基因,幫助研究者深入理解組織中基因間的相互作用模式,進而解析基因功能與組織生物學功能的聯(lián)系,篩選遺傳控制的靶標基因。在大樣本轉(zhuǎn)錄組數(shù)據(jù)中,傳統(tǒng)的差異分析或趨勢分析,僅羅列了單一的基因列表,無法得到基因間的相關性,并對基因進行有效的分類,亦不能推測新的基因相互作用關系。WGCNA分析從系統(tǒng)角度出發(fā),將協(xié)同變化的基因進行聚類,以網(wǎng)絡模塊的形式呈現(xiàn)基因的共表達特性。其優(yōu)勢在于,將海量的數(shù)據(jù)進行歸納和整理,高效研究基因整體表達規(guī)律,并系統(tǒng)地反饋樣本中基因間的相互作用模式,同時可特異性地篩選出與樣本性狀相關的基因模塊和hub基因。作為一種高效的數(shù)據(jù)挖掘手段,其在多樣本RNA-Seq數(shù)據(jù)的應用中越來越流行(Zhaoetal.,2010;巨飛燕等,2019)。WGCNA分析可為蚊蟲不同組織的基因共表達模式研究提供新思路。

    本研究利用WGCNA方法對埃及伊蚊成蚊不同組織的基因共表達模式進行研究,構建了埃及伊蚊9種組織加權基因共表達網(wǎng)絡,結果顯示埃及伊蚊這9種組織涉及11個基因共表達模塊(圖2);將轉(zhuǎn)錄組數(shù)據(jù)對應的不同組織(雌蚊:腦、觸角、喙、下顎須、卵巢;雄蚊:腦、觸角、前足、中足、后足、腹部末端)作為研究的性狀,成功鑒定到與雌蚊觸角、喙、下顎須、卵巢以及雄蚊腦、腹部末端相關的6個組織特異性模塊,分別為green,purple,red,turquoise,blue和magenta模塊(圖2)。蚊蟲分化成熟的多個組織共同支撐其追尋寄主、吸取血液、交配產(chǎn)卵、傳播疾病等生命活動過程。研究表明,蚊蟲可利用多種化學和物理線索,如宿主釋放的氣味、CO2、周圍環(huán)境的化學信號以及雌蚊翅產(chǎn)生的特定頻率等,并據(jù)此定位宿主的位置、選擇配偶和找尋產(chǎn)卵地點(Hartberg,1971;Catoretal.,2009;Matthewsetal.,2016)。觸角、喙和下顎須在此過程中起著極其重要的作用,因此它們被認為是蚊蟲防治控制的重要靶標之一(張晶晶,2019)。其中觸角是蚊蟲最主要的嗅覺器官,在成蚊觸角上有許多被稱為感受器的感覺結構,這是進行化學感知的物理位置(張靜等,2019);喙是蚊蟲在進食過程中處理味覺代碼的重要頭部附屬物,其作為探針到達宿主皮膚下的血管(Maekawaetal.,2011)。下顎須是蚊蟲主要的化學感覺器官,可感知CO2和揮發(fā)性氣味物質(zhì)等(Grantetal.,1995)。腦是蚊蟲神經(jīng)調(diào)節(jié)的重要場所,也是目前蚊蟲生物防控研究的熱點之一,其內(nèi)含有多種神經(jīng)遞質(zhì),如乙酰膽堿、GABA以及谷氨酸等(Matthewsetal.,2016)。在腦的高級處理中心可整合嗅覺和其他感覺模式,最終將這些信號轉(zhuǎn)化為行為,其可能是蚊蟲防控的一個新靶標(Riabininaetal.,2016)。而卵巢、腹部末端是與蚊蟲生殖細胞發(fā)育、繁殖相關的重要場所,有研究認為可通過將雌蚊轉(zhuǎn)化為無害的雄蚊的策略而實現(xiàn)蚊蟲防控,提示調(diào)節(jié)蚊蟲性別分化過程可能是當前防控蚊蟲的一個重要突破口(Halletal.,2015)。對6個組織特異性模塊內(nèi)基因的注釋分析發(fā)現(xiàn),與雌成蚊觸角特異性相關的green模塊具有氣味結合和嗅覺受體活性等功能;與雌成蚊喙特異性相關的purple模塊具有絲氨酸型肽鏈內(nèi)切酶活性和絲氨酸水解酶活性等功能;與雌成蚊下顎須正相關的red模塊具有絲氨酸型肽鏈內(nèi)切酶活性和碳水化合物結合等功能;與雄成蚊腦特異性相關的blue模塊在生物學過程調(diào)節(jié)、信號轉(zhuǎn)導、神經(jīng)系統(tǒng)過程等生物學過程中發(fā)揮主要作用;與雌成蚊卵巢特異性相關的turquoise模塊在DNA復制、細胞分裂、細胞周期等生物學過程中發(fā)揮主要作用;與雄成蚊腹部末端特異性相關的magenta模塊在羧酸代謝、草酸代謝、有機酸代謝等生物學過程中發(fā)揮主要作用(表1和2)。其均能與對應組織功能相聯(lián)系,說明組織特異性模塊的識別可靠,可用于后續(xù)相關信息的挖掘。

    核心基因在基因表達網(wǎng)絡中具有高的連接度,其往往是重要的作用靶點和研究熱點,在生物學過程中起著關鍵作用。但通過Cytoscape構建基因互作網(wǎng)絡后,較多的node和edge常會掩蓋真正的核心基因,影響篩選過程。因此本研究采用Cytoscape軟件中的CytoHubba插件,進一步縮小了6個組織特異性模塊中核心基因的搜索范圍,共篩選出47個在模塊內(nèi)可能起重要作用的hub基因(表3)。這些基因大多可在黑腹果蠅Drosophilamelanogaster中找到對應直系同源物。green模塊內(nèi)AAEL005810,AAEL007669,AAEL010214和AAEL013197在黑腹果蠅中的直系同源基因分別為CG14085,Antdh,CCDC151和CG4362;blue模塊內(nèi)hub基因在黑腹果蠅中均有直系同源物;turquoise模塊內(nèi)AAEL008595,AAEL000961和AAEL003226在黑腹果蠅中有直系同源基因,分別為Mael,Clos和Cpr49Aa-Cpr65Az;red模塊內(nèi)AAEL007033,AAEL000311和AAEL013417在黑腹果蠅中直系同源基因分別為P5cr-2,Fon和CG30280;magenta模塊內(nèi)AAEL006427,AAEL006429,AAEL017567和AAEL004944在黑腹果蠅中的直系同源基因分別為CG17571-thetaTry,thetaTry-CG17571,Nep5和Muc12Ea(Matthewsetal.,2016)。47個hub基因中,29個基因為已知功能基因。如與雌成蚊觸角特異性相關的green模塊中,共篩選出9個hub基因(表3),其中AAEL010426和AAEL001510被描述為氣味受體家族基因(Heetal.,2015;Matthewsetal.,2016);AAEL014089屬于離子受體家族基因;AAEL010214編碼卷曲螺旋結構域蛋白,AAEL012424編碼微管蛋白。與雌成蚊喙特異性相關的purple模塊中,共篩選出6個hub基因(表3),其中AAEL000011和AAEL000039屬于離子受體家族基因(Matthewsetal.,2016);AAEL002600,AAEL013033和AAEL013032被描述為絲氨酸蛋白酶家族基因(Bartholomayetal.,2004;Matthewsetal.,2016)。與雌成蚊下顎須特異性相關的red模塊中,共篩選出9個hub基因(表3),其中AAEL011624屬于絲氨酸蛋白酶家族基因,AAEL000044屬于鳥氨酸脫羧酶家族基因,AAEL013417編碼血管生成素相關蛋白(Matthewsetal.,2016);AAEL007033編碼A-凝集素錨定亞基,AAEL003405為微管動態(tài)蛋白的調(diào)節(jié)因子,AAEL006920被描述為G蛋白偶聯(lián)受體家族基因(Bartholomayetal.,2004;Matthewsetal.,2016)。篩選到的hub基因功能與組織功能相關,表明篩選到的基因較為可靠,但有的hub基因的報道研究甚少,提示后續(xù)可進一步通過blast,RT-PCR以及RNA干擾等手段對篩選到的hub基因功能進行預測和檢驗。Matthews等(2016)認為2個與氨轉(zhuǎn)運有關的基因AAEL007373和AAEL008046可能參與埃及伊蚊尋找宿主過程中氨的檢測,本研究WGCNA分析結果發(fā)現(xiàn)AAEL007373被劃分到green模塊,其在雌/雄蚊觸角中高表達,AAEL008046被劃分到blue模塊,其在雌/雄蚊腦中高表達,這與Matthews等(2016)報道相符,此外,本研究構建了以這兩個基因為核心的局部網(wǎng)絡,揭示其基因網(wǎng)絡關系(圖5)。AAEL000252被報道為雌蚊富集表達基因(Matthewsetal.,2016),在本研究中被劃分到turquoise模塊,其在雌成蚊卵巢中高表達;AAEL010631被報道為雄蚊富集表達基因(Matthewsetal.,2016),被劃分到magenta模塊,其在雄成蚊腹部末端中高表達。

    本研究是基于WGCNA方法和轉(zhuǎn)錄組數(shù)據(jù)對蚊蟲不同組織基因共表達模式的首次探索。構建了埃及伊蚊成蚊9種組織加權基因共表達網(wǎng)絡,獲得了11個基因共表達模塊,其中,6個模塊與組織特異性相關。完成了組織特異性模塊內(nèi)基因的功能富集分析,發(fā)現(xiàn)這些模塊的生物學功能均與對應組織功能密切相關。此外,鑒定了6個組織特異性模塊內(nèi)共47個hub基因,揭示了其基因互作網(wǎng)絡關系。本研究結果可為蚊蟲基因共表達模式分析提供方法基礎,對探究蚊蟲不同組織特有的基因資源信息以及功能基因生物信息研究有參考價值。

    猜你喜歡
    伊蚊共表達觸角
    長觸角小怪獸
    羥哌酯驅(qū)蚊液和檸檬桉醇驅(qū)蚊液驅(qū)避白紋伊蚊的效果對比
    化工管理(2021年7期)2021-05-13 00:45:12
    鄭宇婷等發(fā)表于本刊2021年37卷第2期《利用ssr標記淺析云南登革熱重點地區(qū)埃及伊蚊種群遺傳特征》一文的勘誤
    云南省邊境地區(qū)埃及伊蚊和白紋伊蚊孳生特征調(diào)查*
    侵襲性垂體腺瘤中l(wèi)ncRNA-mRNA的共表達網(wǎng)絡
    夯實為僑服務機制 多方觸角深入基層
    華人時刊(2020年19期)2021-01-14 01:17:20
    試探著向硅谷伸出觸角
    能源(2018年5期)2018-06-15 08:56:20
    膀胱癌相關lncRNA及其共表達mRNA的初步篩選與功能預測
    中國流行株HIV-1gag-gp120與IL-2/IL-6共表達核酸疫苗質(zhì)粒的構建和實驗免疫研究
    胃癌患者癌組織HIF-1α、TGF-β共表達及其臨床意義
    久久鲁丝午夜福利片| 久久国产精品人妻蜜桃| 国产精品一国产av| 肉色欧美久久久久久久蜜桃| 国产伦理片在线播放av一区| 热99久久久久精品小说推荐| 青春草亚洲视频在线观看| 午夜av观看不卡| 在线亚洲精品国产二区图片欧美| 性色av乱码一区二区三区2| 亚洲专区中文字幕在线| 国产精品熟女久久久久浪| 免费在线观看黄色视频的| 51午夜福利影视在线观看| 亚洲人成电影免费在线| 欧美老熟妇乱子伦牲交| 五月天丁香电影| 午夜免费鲁丝| 黄片播放在线免费| 国产视频首页在线观看| 亚洲av男天堂| 日本一区二区免费在线视频| 国产成人精品无人区| 国产精品.久久久| 国产一卡二卡三卡精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产一级毛片在线| 国语对白做爰xxxⅹ性视频网站| 国产视频一区二区在线看| 男人添女人高潮全过程视频| 精品福利永久在线观看| 欧美中文综合在线视频| 水蜜桃什么品种好| 国产主播在线观看一区二区 | 欧美日韩福利视频一区二区| 黄色毛片三级朝国网站| 天天躁夜夜躁狠狠躁躁| 人人妻人人澡人人爽人人夜夜| 午夜av观看不卡| 制服人妻中文乱码| 激情视频va一区二区三区| 国产av一区二区精品久久| 天堂中文最新版在线下载| 日本wwww免费看| 91精品国产国语对白视频| 另类精品久久| 女性生殖器流出的白浆| 久久久精品国产亚洲av高清涩受| 亚洲成人手机| 精品人妻一区二区三区麻豆| 久久人妻熟女aⅴ| 尾随美女入室| 国产精品一区二区免费欧美 | 久久久久国产一级毛片高清牌| 亚洲国产av影院在线观看| 最新在线观看一区二区三区 | 五月天丁香电影| 国产精品一二三区在线看| 777久久人妻少妇嫩草av网站| 国产视频一区二区在线看| 韩国精品一区二区三区| 999久久久国产精品视频| 久久久久久久精品精品| 日日夜夜操网爽| 青青草视频在线视频观看| 我的亚洲天堂| 午夜日韩欧美国产| 亚洲欧洲精品一区二区精品久久久| 亚洲国产中文字幕在线视频| 黄片播放在线免费| av不卡在线播放| 国产精品99久久99久久久不卡| 欧美亚洲日本最大视频资源| 午夜精品国产一区二区电影| 妹子高潮喷水视频| 啦啦啦在线免费观看视频4| 国产精品欧美亚洲77777| 一区在线观看完整版| 91国产中文字幕| 国产熟女午夜一区二区三区| 后天国语完整版免费观看| 欧美人与善性xxx| 亚洲av成人不卡在线观看播放网 | 亚洲视频免费观看视频| 久久99精品国语久久久| 国产免费福利视频在线观看| 啦啦啦在线观看免费高清www| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区久久| av线在线观看网站| 中文字幕人妻丝袜制服| 午夜久久久在线观看| 男人爽女人下面视频在线观看| 美女中出高潮动态图| 亚洲av欧美aⅴ国产| 国产男女内射视频| 国产无遮挡羞羞视频在线观看| 99国产精品免费福利视频| 亚洲熟女精品中文字幕| 丰满少妇做爰视频| 欧美av亚洲av综合av国产av| 十八禁人妻一区二区| 伦理电影免费视频| 亚洲精品在线美女| 国产三级黄色录像| 亚洲成国产人片在线观看| 99国产精品99久久久久| 国产亚洲av片在线观看秒播厂| 叶爱在线成人免费视频播放| 伦理电影免费视频| 成人国产一区最新在线观看 | 美国免费a级毛片| 午夜福利视频在线观看免费| 亚洲av国产av综合av卡| av国产精品久久久久影院| 亚洲专区中文字幕在线| 亚洲欧美一区二区三区黑人| 精品国产一区二区三区四区第35| 极品少妇高潮喷水抽搐| 七月丁香在线播放| 久久影院123| 在线 av 中文字幕| 国产野战对白在线观看| 免费在线观看影片大全网站 | 巨乳人妻的诱惑在线观看| 考比视频在线观看| 夫妻性生交免费视频一级片| 亚洲精品自拍成人| 人人妻人人澡人人看| 激情五月婷婷亚洲| 亚洲av片天天在线观看| 久久国产精品影院| 丁香六月天网| 亚洲久久久国产精品| av网站免费在线观看视频| 国产精品秋霞免费鲁丝片| 制服人妻中文乱码| 少妇 在线观看| 国产精品麻豆人妻色哟哟久久| 大香蕉久久网| 亚洲人成电影观看| 9191精品国产免费久久| 欧美黑人精品巨大| 精品国产一区二区久久| 一二三四在线观看免费中文在| 深夜精品福利| 久久久精品国产亚洲av高清涩受| 97精品久久久久久久久久精品| 人人妻人人爽人人添夜夜欢视频| 日本vs欧美在线观看视频| 国产亚洲av高清不卡| 亚洲精品美女久久久久99蜜臀 | 免费少妇av软件| 最新的欧美精品一区二区| 校园人妻丝袜中文字幕| 精品亚洲成国产av| 又大又爽又粗| 精品卡一卡二卡四卡免费| 亚洲欧美色中文字幕在线| 国产免费福利视频在线观看| 国产99久久九九免费精品| 男女午夜视频在线观看| 欧美精品人与动牲交sv欧美| 国产成人精品久久二区二区免费| 久热这里只有精品99| 国产熟女午夜一区二区三区| 欧美日韩av久久| 精品少妇黑人巨大在线播放| 极品人妻少妇av视频| 免费日韩欧美在线观看| 欧美日本中文国产一区发布| 桃花免费在线播放| 黄色视频在线播放观看不卡| 日韩,欧美,国产一区二区三区| 性少妇av在线| 一二三四在线观看免费中文在| 亚洲国产日韩一区二区| 精品欧美一区二区三区在线| 欧美成人精品欧美一级黄| 狠狠婷婷综合久久久久久88av| 免费av中文字幕在线| 午夜视频精品福利| 18在线观看网站| 99国产精品99久久久久| a级片在线免费高清观看视频| 在线观看免费午夜福利视频| 手机成人av网站| 国产精品二区激情视频| 中文精品一卡2卡3卡4更新| 国产无遮挡羞羞视频在线观看| 叶爱在线成人免费视频播放| av国产精品久久久久影院| 亚洲欧美色中文字幕在线| 大香蕉久久成人网| 国语对白做爰xxxⅹ性视频网站| 国产欧美日韩一区二区三 | 婷婷丁香在线五月| 久久久久久人人人人人| 1024视频免费在线观看| 2021少妇久久久久久久久久久| 欧美老熟妇乱子伦牲交| 七月丁香在线播放| 中文字幕制服av| 大型av网站在线播放| 青青草视频在线视频观看| 波多野结衣av一区二区av| 欧美亚洲 丝袜 人妻 在线| 精品高清国产在线一区| 成人午夜精彩视频在线观看| 亚洲免费av在线视频| 国产日韩欧美视频二区| 国产熟女欧美一区二区| 一区二区三区四区激情视频| 另类亚洲欧美激情| 亚洲国产av影院在线观看| 亚洲精品成人av观看孕妇| 人人妻人人澡人人爽人人夜夜| 欧美少妇被猛烈插入视频| 国产欧美日韩一区二区三区在线| 亚洲精品久久久久久婷婷小说| 久久精品国产综合久久久| 热99久久久久精品小说推荐| 色综合欧美亚洲国产小说| 飞空精品影院首页| 国产熟女欧美一区二区| 91精品伊人久久大香线蕉| 亚洲人成网站在线观看播放| 精品国产乱码久久久久久小说| 老司机在亚洲福利影院| 欧美激情高清一区二区三区| 国产1区2区3区精品| av电影中文网址| 欧美变态另类bdsm刘玥| 中文精品一卡2卡3卡4更新| 久久久久久久国产电影| 老司机靠b影院| 一区福利在线观看| 国产成人免费观看mmmm| 在线亚洲精品国产二区图片欧美| 亚洲激情五月婷婷啪啪| 成年人午夜在线观看视频| 老汉色av国产亚洲站长工具| 亚洲av电影在线观看一区二区三区| 人人妻,人人澡人人爽秒播 | 中国国产av一级| 国产成人系列免费观看| 亚洲国产av影院在线观看| 亚洲中文字幕日韩| 国产精品 国内视频| 91九色精品人成在线观看| 99re6热这里在线精品视频| 高清av免费在线| 天天躁夜夜躁狠狠躁躁| 十八禁人妻一区二区| av不卡在线播放| 女人久久www免费人成看片| 一级毛片电影观看| 精品人妻一区二区三区麻豆| 国产成人一区二区三区免费视频网站 | 女人被躁到高潮嗷嗷叫费观| 超色免费av| 欧美人与性动交α欧美软件| 国产av精品麻豆| 黄色a级毛片大全视频| av欧美777| 国产一区二区在线观看av| 精品国产超薄肉色丝袜足j| 欧美在线一区亚洲| 国产精品久久久久成人av| 国产在线视频一区二区| 我的亚洲天堂| 国产成人一区二区在线| 精品一区在线观看国产| 久久久精品免费免费高清| 一级毛片 在线播放| 每晚都被弄得嗷嗷叫到高潮| 只有这里有精品99| 国产国语露脸激情在线看| 一级黄片播放器| 久久精品国产亚洲av高清一级| 亚洲中文日韩欧美视频| 国产欧美日韩精品亚洲av| 国产精品 欧美亚洲| 亚洲国产毛片av蜜桃av| 色婷婷av一区二区三区视频| 91字幕亚洲| 久久精品亚洲av国产电影网| 麻豆乱淫一区二区| 久久免费观看电影| a 毛片基地| 美女脱内裤让男人舔精品视频| 脱女人内裤的视频| 日韩,欧美,国产一区二区三区| 久久久久精品国产欧美久久久 | 国产一卡二卡三卡精品| 国产精品 欧美亚洲| 亚洲精品一区蜜桃| 免费观看a级毛片全部| 少妇 在线观看| 在线看a的网站| 亚洲av片天天在线观看| 精品久久久久久久毛片微露脸 | 日本欧美国产在线视频| 中文字幕人妻丝袜制服| 国产精品久久久久成人av| 午夜福利在线免费观看网站| 丝袜美足系列| 国产免费现黄频在线看| 91精品国产国语对白视频| 91精品三级在线观看| 制服人妻中文乱码| 国产伦理片在线播放av一区| 两个人看的免费小视频| 久久毛片免费看一区二区三区| 成人三级做爰电影| 日韩电影二区| av又黄又爽大尺度在线免费看| 亚洲国产精品国产精品| 老司机午夜十八禁免费视频| 男人舔女人的私密视频| 欧美日韩福利视频一区二区| 亚洲精品国产av蜜桃| 亚洲天堂av无毛| 日韩一区二区三区影片| 久久精品亚洲av国产电影网| 国产精品国产三级国产专区5o| 国产精品久久久久久人妻精品电影 | 欧美日韩福利视频一区二区| 狠狠精品人妻久久久久久综合| 国产黄频视频在线观看| 亚洲中文字幕日韩| 婷婷色麻豆天堂久久| 欧美精品亚洲一区二区| 日韩人妻精品一区2区三区| 国产老妇伦熟女老妇高清| 久久鲁丝午夜福利片| 国产男女内射视频| 欧美国产精品va在线观看不卡| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利在线免费观看网站| 国产一区二区三区av在线| 亚洲色图综合在线观看| 亚洲国产看品久久| 午夜老司机福利片| 中文乱码字字幕精品一区二区三区| 国产成人精品在线电影| 91老司机精品| 搡老岳熟女国产| 麻豆国产av国片精品| 手机成人av网站| 久久精品人人爽人人爽视色| 国产高清videossex| 男人舔女人的私密视频| 日韩av在线免费看完整版不卡| 啦啦啦在线免费观看视频4| 一区福利在线观看| 欧美黄色淫秽网站| 亚洲国产最新在线播放| 久热这里只有精品99| 亚洲av成人精品一二三区| 美女国产高潮福利片在线看| 日韩欧美一区视频在线观看| 99热国产这里只有精品6| 欧美性长视频在线观看| 精品少妇黑人巨大在线播放| 亚洲成国产人片在线观看| 国产成人欧美| 精品免费久久久久久久清纯 | 国产日韩一区二区三区精品不卡| 一区福利在线观看| 国产97色在线日韩免费| 国产日韩欧美视频二区| 在线看a的网站| 日日夜夜操网爽| 国产av一区二区精品久久| 亚洲免费av在线视频| 亚洲国产成人一精品久久久| 在线观看免费日韩欧美大片| 亚洲精品国产av成人精品| 亚洲人成网站在线观看播放| 性少妇av在线| 成人国产一区最新在线观看 | 一级黄片播放器| 中文字幕人妻丝袜制服| 久久青草综合色| 国产精品香港三级国产av潘金莲 | 麻豆乱淫一区二区| 99国产综合亚洲精品| 成人18禁高潮啪啪吃奶动态图| 女人被躁到高潮嗷嗷叫费观| 国产女主播在线喷水免费视频网站| xxxhd国产人妻xxx| 午夜福利视频在线观看免费| 日本av免费视频播放| 久久久久网色| 久久女婷五月综合色啪小说| svipshipincom国产片| 精品亚洲成a人片在线观看| 色精品久久人妻99蜜桃| 男男h啪啪无遮挡| 欧美激情 高清一区二区三区| 777米奇影视久久| 飞空精品影院首页| 电影成人av| 一本久久精品| 老汉色av国产亚洲站长工具| 免费久久久久久久精品成人欧美视频| 亚洲 国产 在线| 国产成人影院久久av| 夜夜骑夜夜射夜夜干| 国产精品.久久久| 国产一区二区 视频在线| 国产亚洲精品第一综合不卡| 精品国产乱码久久久久久小说| 久久久久国产精品人妻一区二区| 国产成人精品在线电影| 国产99久久九九免费精品| 1024香蕉在线观看| 国产男女内射视频| 十八禁高潮呻吟视频| 黄色片一级片一级黄色片| 少妇粗大呻吟视频| 欧美亚洲 丝袜 人妻 在线| 亚洲欧洲日产国产| 欧美国产精品va在线观看不卡| 亚洲成人国产一区在线观看 | 每晚都被弄得嗷嗷叫到高潮| 久久99精品国语久久久| 老司机亚洲免费影院| 90打野战视频偷拍视频| 精品一区二区三区av网在线观看 | 另类精品久久| 母亲3免费完整高清在线观看| 一本综合久久免费| 黄色视频不卡| 国产极品粉嫩免费观看在线| 777久久人妻少妇嫩草av网站| 中文字幕人妻丝袜制服| 久久精品国产综合久久久| 国产精品九九99| 免费在线观看日本一区| 国产成人一区二区在线| 久久ye,这里只有精品| 国产又爽黄色视频| 在线观看免费日韩欧美大片| 狂野欧美激情性xxxx| 日韩av免费高清视频| 亚洲少妇的诱惑av| 午夜老司机福利片| 亚洲欧美日韩高清在线视频 | 男女高潮啪啪啪动态图| 日日摸夜夜添夜夜爱| 欧美精品高潮呻吟av久久| 国产精品亚洲av一区麻豆| 一边摸一边抽搐一进一出视频| 九草在线视频观看| 看免费av毛片| 亚洲av美国av| 精品一区二区三区av网在线观看 | 人人妻人人爽人人添夜夜欢视频| 青春草亚洲视频在线观看| 一区福利在线观看| 国产精品三级大全| 少妇精品久久久久久久| 中文字幕人妻熟女乱码| 久久精品久久久久久久性| 亚洲精品日韩在线中文字幕| 人妻人人澡人人爽人人| av有码第一页| 校园人妻丝袜中文字幕| 亚洲国产精品999| 久9热在线精品视频| 欧美日韩av久久| 中文字幕人妻丝袜制服| 99国产精品99久久久久| a级片在线免费高清观看视频| 午夜激情久久久久久久| 欧美+亚洲+日韩+国产| 大香蕉久久成人网| a 毛片基地| 欧美日韩亚洲综合一区二区三区_| 久久国产亚洲av麻豆专区| 国产三级黄色录像| 老司机亚洲免费影院| 国产精品 国内视频| 亚洲欧美激情在线| 精品国产乱码久久久久久小说| 精品人妻1区二区| xxx大片免费视频| av视频免费观看在线观看| 亚洲国产精品一区三区| 女人久久www免费人成看片| 亚洲自偷自拍图片 自拍| 99国产精品一区二区蜜桃av | 久久狼人影院| 国产精品一区二区在线不卡| 中文字幕人妻丝袜一区二区| 国产免费视频播放在线视频| 夫妻性生交免费视频一级片| 国产免费一区二区三区四区乱码| 日韩,欧美,国产一区二区三区| 亚洲精品国产色婷婷电影| 国产亚洲精品第一综合不卡| 亚洲熟女精品中文字幕| 黄色视频在线播放观看不卡| 色网站视频免费| 男女高潮啪啪啪动态图| 久久亚洲精品不卡| 可以免费在线观看a视频的电影网站| 久久精品国产a三级三级三级| 后天国语完整版免费观看| 久久久久国产精品人妻一区二区| 亚洲欧美一区二区三区黑人| 久久国产亚洲av麻豆专区| 熟女av电影| 久久久精品94久久精品| 欧美97在线视频| 日韩大片免费观看网站| 午夜精品国产一区二区电影| 嫁个100分男人电影在线观看 | 成人免费观看视频高清| 亚洲精品av麻豆狂野| 91精品三级在线观看| 丰满迷人的少妇在线观看| h视频一区二区三区| 成人三级做爰电影| 午夜影院在线不卡| 成人国产av品久久久| 午夜福利免费观看在线| 校园人妻丝袜中文字幕| 亚洲av美国av| 9色porny在线观看| 尾随美女入室| 亚洲av电影在线观看一区二区三区| 日韩电影二区| 亚洲情色 制服丝袜| 90打野战视频偷拍视频| 亚洲精品一卡2卡三卡4卡5卡 | av网站免费在线观看视频| 久久久久久免费高清国产稀缺| 菩萨蛮人人尽说江南好唐韦庄| 巨乳人妻的诱惑在线观看| 久久国产精品人妻蜜桃| 亚洲人成77777在线视频| 在现免费观看毛片| 中文欧美无线码| 黑人巨大精品欧美一区二区蜜桃| 啦啦啦在线观看免费高清www| 黄片小视频在线播放| 王馨瑶露胸无遮挡在线观看| 国产成人精品久久二区二区91| 老汉色∧v一级毛片| 成人国语在线视频| 熟女av电影| 9191精品国产免费久久| 麻豆av在线久日| 18禁观看日本| 美女扒开内裤让男人捅视频| 国产欧美亚洲国产| 在现免费观看毛片| 国产一区二区三区av在线| 丁香六月天网| 欧美大码av| 美女脱内裤让男人舔精品视频| 亚洲专区国产一区二区| 国产国语露脸激情在线看| 人妻 亚洲 视频| 免费一级毛片在线播放高清视频 | 国产在线免费精品| 欧美日韩av久久| 亚洲国产av影院在线观看| 亚洲成人国产一区在线观看 | 国产精品偷伦视频观看了| 国产在线观看jvid| 丁香六月欧美| 日韩大码丰满熟妇| 日本欧美视频一区| 日日爽夜夜爽网站| 国产黄色免费在线视频| 少妇精品久久久久久久| 日日夜夜操网爽| 1024视频免费在线观看| av在线老鸭窝| 黄色毛片三级朝国网站| 又大又爽又粗| 亚洲第一av免费看| 国产成人a∨麻豆精品| 这个男人来自地球电影免费观看| av在线老鸭窝| 日日爽夜夜爽网站| 国产成人影院久久av| 国产日韩欧美在线精品| 久久国产亚洲av麻豆专区| 国产成人影院久久av| 国产又爽黄色视频| 中文字幕高清在线视频| 久久性视频一级片| 精品一品国产午夜福利视频| 黄色毛片三级朝国网站| 麻豆av在线久日| 午夜影院在线不卡| 亚洲中文日韩欧美视频| 国产精品成人在线| 精品一品国产午夜福利视频| 欧美乱码精品一区二区三区| 国产熟女午夜一区二区三区| 国产午夜精品一二区理论片| 我要看黄色一级片免费的| 免费在线观看黄色视频的| 亚洲第一青青草原| 一边摸一边抽搐一进一出视频| 性少妇av在线| 日韩,欧美,国产一区二区三区| 亚洲精品第二区| 亚洲综合色网址|