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

    小地老虎不同組織表達(dá)差異揭示其翅發(fā)育機(jī)制

    2019-10-24 09:11:12鄭偉剛靳明輝SwapanChakrabarty肖彬蕭玉濤袁海濱
    生物技術(shù)通報(bào) 2019年10期
    關(guān)鍵詞:差異基因胸部成蟲

    鄭偉剛 靳明輝 Swapan Chakrabarty 肖彬 蕭玉濤 袁海濱

    (1. 吉林農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,長春 130118;2. 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)基因組研究所,深圳 518120)

    小地老虎(Agrotis ipsilon)、棉鈴蟲(Helicoverpa armigera)、帝王蝶、黏蟲(Mythimna separata)、稻飛虱(Nilaparvata lugens)、玉米螟(Ostriniafurnacalis、Ostrinia nubilalis)及東亞飛蝗(Locusta migratoria manilensis(Meyen))等都屬于遷飛性昆蟲[1-6]。遷飛是動物種群對棲息地資源發(fā)生季節(jié)性或偶然性變化時(shí)進(jìn)化出的一種適應(yīng)性對策,昆蟲的遷飛行為嚴(yán)重影響人類的生產(chǎn)和生活,如危害農(nóng)作物,造成農(nóng)作物大量減產(chǎn)、傳播人和牲畜的疾病等。影響昆蟲遷飛的因素眾多,包括生存環(huán)境、激素調(diào)控、遺傳因素等,而翅發(fā)育是昆蟲遷飛的一個(gè)重要影響因素。昆蟲翅膀的進(jìn)化對昆蟲覓食、求偶、避敵、抵御外界不良環(huán)境等起到重要作用[7-8]。隨著果蠅等模式昆蟲的深入研究,一系列與翅發(fā)育有關(guān)的信號通路逐漸被揭示出來,如Wnt、Hedgehog、Decapentaplegic(Dpp)、Hippo和 c-Jun N-terminal Kinase(JNK)等信號通路,其中Wnt通路是昆蟲翅發(fā)育的關(guān)鍵信號通路。

    Wnt信號通路包含3種Wnt信號傳導(dǎo)通路:經(jīng)典Wnt通路,非經(jīng)典Wnt /鈣通路以及非經(jīng)典planar cell polarity通路[9-11],這 3種途徑都是通過Wnt蛋白配體與Frizzled(Fz)家族受體結(jié)合而激活的。經(jīng)典Wnt通路引起基因轉(zhuǎn)錄的調(diào)節(jié),并且被認(rèn)為部分地被SPATS1基因負(fù)調(diào)節(jié)[11]。非經(jīng)典planar cell polarity通路調(diào)控細(xì)胞骨架,骨架決定細(xì)胞的形狀。非經(jīng)典Wnt /鈣通路,調(diào)控細(xì)胞內(nèi)的鈣。Decapentaplegic(Dpp)信號通路中,Dpp是參與果蠅發(fā)育的關(guān)鍵形態(tài)發(fā)生素,也是第一個(gè)經(jīng)過驗(yàn)證的分泌形態(tài)發(fā)生素[12-13]。早期果蠅胚胎和15個(gè)成蟲器的正確形成和發(fā)育需要Dpp調(diào)控,這些組織將生長發(fā)育為成年果蠅中的四肢和其他器官和結(jié)構(gòu)。DPP在調(diào)節(jié)組織的生長和器官組織大小方面起著重要的作用[12,14]。Hedgehog(Hh)信號通路在控制細(xì)胞增殖、組織模式形成、干細(xì)胞維持和發(fā)育方面發(fā)揮著多種作用[15-18]。Hedgehog信號通路是一種向胚胎細(xì)胞傳遞信息的信號通路,這些信息是胚胎細(xì)胞分化所必需的。胚胎的不同部位具有不同濃度的Hedgehog基因信號蛋白。該通路在成蟲個(gè)體中也有作用。Hippo信號通路也被稱為Salvador或者Warts通路,通過調(diào)控細(xì)胞增殖和凋亡來控制動物的器官大小。該通路的名字來源于它的一個(gè)關(guān)鍵信號元件——蛋白激酶Hippo(Hpo)[19-20]。該基因的突變導(dǎo)致組織過度生長,或類似“河馬”表型。

    小地老虎(Agrotis ipsilon)又名切根蟲、土蠶等,屬于一種世界性的鱗翅目夜蛾科農(nóng)業(yè)害蟲[5,21],也是典型的遷飛性農(nóng)業(yè)害蟲。由于翅是在成蟲階段發(fā)育形成,且屬于昆蟲胸部的組成部分,本試驗(yàn)對小地老虎幼蟲,成蟲頭部、胸部、腹部和進(jìn)行轉(zhuǎn)錄組測序,探究可能決定小地老虎翅發(fā)育的關(guān)鍵通路及相關(guān)基因[22-23]。

    1 材料與方法

    1.1 材料

    小地老虎幼蟲主要來源于中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)基因組研究所農(nóng)業(yè)生態(tài)基因組學(xué)中心實(shí)驗(yàn)室,小地老虎成蟲主要來源于山東長島野外采集。

    1.2 方法

    1.2.1 小地老虎總RNA的提取、文庫構(gòu)建和測序 小地老虎幼蟲6頭作為一個(gè)重復(fù),共取18頭做3組幼蟲整蟲轉(zhuǎn)錄組Ya、Yb、Yc的組織材料,為避免腸道微生物污染,幼蟲整蟲解剖去除圍食膜,再用液氮研磨。小地老虎成蟲每6頭作為一個(gè)重復(fù),分別取頭、胸、腹,做頭(Ta、Tb、Tc)、胸(Xa、Xb、Xc)、腹(Fa、Fb、Fc)轉(zhuǎn)錄組材料。

    取烘干研缽,加入少量乙醇,灼燒防止研缽內(nèi)殘留其他物質(zhì)。待研缽烘干后,加入液氮,將組織樣品加入研缽中研磨至粉狀,取適量放入裝有1 mL TRIzol的RNAase free 1.5 mL離心管中。通常情況下,1 mL TRIzol Reagent 能夠裂解50-100 mg組織樣品。

    向樣品中加入0.2 mL氯仿(每1 mL TRIzol加0.2 mL氯仿)后,用手輕輕搖動15 s左右,使其混合并于常溫孵育2-3 min,然后將樣品于12 000 r/min、4℃條件下進(jìn)行離心15 min。將離心管傾斜為45°并用移液器將水相吸出,然后移入新的離心管,在此過程中應(yīng)盡量避免吸入混合相或者有機(jī)層。取0.5mL 異丙醇(每1 mL TRIzol對應(yīng)0.5 mL異丙醇)加入之前移出的水相樣品中,常溫孵育10 min后在4℃、12 000 r/min條件下,離心10 min,倒掉上清液并留下RNA沉淀。加入1 mL 75%酒精(每1 mL的 TRIzol至少加1 mL的75%乙醇),輕微震蕩混勻后在4℃、7 500 r/min條件下離心5 min,倒掉上清液,只留底部沉淀。風(fēng)干RNA沉淀5-10 min,禁止真空離心RNA沉淀。加入0.02-0.05 mL無核酸酶水,使用移液器反復(fù)吸打,然后在55-60℃條件下水浴10-15 min,使RNA沉淀溶解,將得到的RNA溶液于-80℃條件下長期保存。

    利用 Nanodrop2000 對根據(jù)上述步驟提取出的total RNA,檢測其濃度和純度,利用瓊脂糖凝膠電泳檢測 RNA 完整性,使用Agilent 2100 測定total RNA RIN 值。單次建庫要求RNA 總量 5 μg,濃度≥ 200 μg/mL,OD260/OD2280介于 1.8-2.2 之間。利用Illumina公司的TruSeqTMRNA sample preparation kit完成cDNA文庫的構(gòu)建,并通過Illumina Hiseq 測序平臺的150 paired-end模式進(jìn)行轉(zhuǎn)錄組測序。

    1.2.2 轉(zhuǎn)錄組組裝和ORF預(yù)測 對于無參考基因組的轉(zhuǎn)錄組測序和相關(guān)研究,獲得RNA-seq高質(zhì)量測序數(shù)據(jù)后,需要將測序數(shù)據(jù)進(jìn)行組裝,才能進(jìn)行后續(xù)注釋及其他生物學(xué)功能分析。Trinity(http://trinityrnaseq.sourceforge.net/,版本號 :v2.6.6)是目前從頭組裝Illumina 短片段序列比較常用的軟件,本試驗(yàn)在對原始數(shù)據(jù)進(jìn)行過濾接頭和低質(zhì)量數(shù)據(jù)后,使用Trinity進(jìn)行從頭組裝。使用TransDecoder(https://github.com/TransDecoder/TransDecoder,版本號:5.5.0)軟件識別轉(zhuǎn)錄本序列中候選的編碼區(qū)域,結(jié)合 cd-hit(http://weizhongli-lab.org/cd-hit/,版本號:v4.8.1)進(jìn)行去冗余,得到獨(dú)立基因(Unigene)序列。

    1.2.3 轉(zhuǎn)錄組功能注釋 使用Trinity軟件組裝得到轉(zhuǎn)錄本序列后,要對生物學(xué)相關(guān)問題進(jìn)行探究,就要對這些轉(zhuǎn)錄本進(jìn)行功能注釋。將拼接所得所有(unigene)序列,使用Interperscan(http://www.ebi.ac.uk/interpro/,版本號:v5)進(jìn)行結(jié)構(gòu)域預(yù)測,使用 Blastp(https://www.ncbi.nlm.nih.gov/)分別與KEGG、eggNOG、NR、Trembl、Swissprot、數(shù)據(jù)庫進(jìn)行比對獲得相應(yīng)的功能注釋信息,然后進(jìn)行KEGG通路富集分析。KEGG(Kyoto Encyclopedia of Genes and Genomes,京都基因和基因組百科全書,http://www.genome.jp/kegg/)是基因組注釋方面的公共數(shù)據(jù)庫。KEGG數(shù)據(jù)庫中豐富的通路信息將有助于我們從系統(tǒng)水平去了解基因的生物學(xué)功能,例如代謝途徑、遺傳信息傳遞以及細(xì)胞過程等一些復(fù)雜的生物功能,這大大提高了該數(shù)據(jù)庫在實(shí)際生產(chǎn)和應(yīng)用中的價(jià)值。根據(jù)與KEGG數(shù)據(jù)庫比對結(jié)果,獲得轉(zhuǎn)錄本對應(yīng)的KO編號,根據(jù)KO編號可以獲得某轉(zhuǎn)錄本可能參與的具體生物學(xué)通路。eggNOG(evolutionary genealogy of genes:Non-supervised Orthologous Groups)數(shù)據(jù)庫包含了豐富的注釋信息,包含了KEGG等數(shù)據(jù)庫信息。Swiss-Prot是一個(gè)手工注釋的非冗余蛋白序列數(shù)據(jù)庫。它結(jié)合了從科學(xué)文獻(xiàn)中提取的信息和生物尿酸評估計(jì)算分析,TrEMBL(翻譯的EMBL核苷酸序列數(shù)據(jù)庫),為那些不在Swiss-Prot中的蛋白質(zhì)提供自動注釋。NR(NCBI 非冗余蛋白庫)為綜合數(shù)據(jù)庫,包括了SwissProt、PIR(Protein information resource)、PRF(Protein research foundation)、PDB(Protein data bank)蛋白質(zhì)數(shù)據(jù)庫中非冗余的數(shù)據(jù)以及從GenBank和RefSeq的CDS數(shù)據(jù)庫中翻譯所得的蛋白質(zhì)數(shù)據(jù)。

    1.2.4 差異表達(dá)基因分析 使用組裝軟件將轉(zhuǎn)錄組組裝好后,利用RSEM(http://deweylab.biostat.wisc.edu/rsem/)進(jìn)行表達(dá)量分析。RSEM利用bowtie的比對結(jié)果進(jìn)行表達(dá)量統(tǒng)計(jì),首先得到每個(gè)樣品比對到每個(gè)基因上的read count數(shù)目,然后對其進(jìn)行TPM數(shù)值轉(zhuǎn)換,TPM數(shù)值可以反映基因的表達(dá)水平。使用R包(edgeR和DESeq)分別對小地老虎轉(zhuǎn)錄組進(jìn)行基因表達(dá)差異分析,取兩者交集,得到最終差異基因。

    2 結(jié)果

    2.1 原始測序數(shù)據(jù)質(zhì)控

    測序小地老虎幼蟲,成蟲頭、胸和腹部,各3個(gè)樣本,共計(jì)12個(gè)樣本。小地老虎幼蟲共測數(shù)據(jù)19.59 G,小地老虎成蟲頭部20.73 G,小地老虎成蟲胸部21.36 G,小地老虎成蟲腹部21.32 G,共計(jì)數(shù)據(jù)量82.99 G,質(zhì)控去掉Adapter和低于Q20的數(shù)據(jù)后最終得到78.67 G(表1)。

    2.2 轉(zhuǎn)錄組組裝和ORF預(yù)測

    使用 Trinity(http://trinityrnaseq.sourceforge.net/)軟件,對小地老虎轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行組裝。最終組裝得到轉(zhuǎn)錄組243.4 M,最長序列27.4 K,平均長度0.6 K,N50 0.7 K,denove轉(zhuǎn)錄組組裝結(jié)果GC含量為39.66%(表2)。使用TransDecoder對轉(zhuǎn)錄組進(jìn)行ORF 預(yù)測,并去冗余后,得到37 744個(gè)unigenes。轉(zhuǎn)錄組組裝長度頻數(shù)分布(圖1),圖1-A為轉(zhuǎn)錄本(Transcripts)組裝長度頻數(shù)分布,在200-500范圍內(nèi)的轉(zhuǎn)錄本比例最多,為67.59%,其次為在500-1 000范圍內(nèi)的轉(zhuǎn)錄本比例為22.03%。圖1-B為Unigenes組裝長度頻數(shù)分布,在500-1 000范圍內(nèi)的轉(zhuǎn)錄本比例最多,為33.83%,其次為在200-500范圍內(nèi)的比例為30.29%。

    表1 小地老虎轉(zhuǎn)錄組測序原始數(shù)據(jù)及質(zhì)控?cái)?shù)據(jù)統(tǒng)計(jì)

    表2 轉(zhuǎn)錄組denovo組裝結(jié)果統(tǒng)計(jì)

    2.3 小地老虎轉(zhuǎn)錄組功能注釋

    使用Interperscan進(jìn)行結(jié)構(gòu)域預(yù)測,共有20 522個(gè)unigenes預(yù)測出結(jié)構(gòu)域信息。將拼接所得獨(dú)立基因,使用Blastp分別與NR、Sport、Trembl、KEGG數(shù)據(jù)庫進(jìn)行比對,比對后取besthit獲得相應(yīng)的注釋信息。

    KEGG 數(shù)據(jù)庫注釋功能基因23 334個(gè),注釋通路373條,KO 6 161個(gè);eggNOG 數(shù)據(jù)路注釋功能基因31 556個(gè),Orthologous Groups(OG,同源群)8 211個(gè);Nr 數(shù)據(jù)庫共注釋功能基因29 971個(gè);Swissprot 數(shù)據(jù)庫注釋功能基因本21 538個(gè);Trembl數(shù)據(jù)庫注釋功能基因32 354個(gè)。各個(gè)數(shù)據(jù)庫功能注釋Venn圖(圖2),其中5個(gè)數(shù)據(jù)庫能夠共同注釋到的基因有18 345個(gè)。5個(gè)數(shù)據(jù)庫分別注釋出了本數(shù)據(jù)庫特有基因,其中Trembl數(shù)據(jù)庫注釋出的特有基因最多,為426個(gè),KEGG數(shù)據(jù)庫注釋出的特有基因最少,有15個(gè)(圖2)。

    圖1 轉(zhuǎn)錄組組裝長度頻數(shù)分布

    2.4 小地老虎轉(zhuǎn)錄組表達(dá)差異分析

    分別使用EdgeR.R和DEseq2.R 進(jìn)行差異分析,并求出差異基因。再將兩種方法得到的差異基因求交集,得到最終的差異基因(圖3)。最終得到小地老虎成蟲頭部樣本與小地老虎胸部樣本轉(zhuǎn)錄組差異基因4 069個(gè),小地老虎成蟲胸部樣本與小地老虎腹部樣本轉(zhuǎn)錄組差異基因6 905個(gè),小地老虎成蟲腹部樣本與小地老虎頭部樣本轉(zhuǎn)錄組差異基因7 777個(gè),小地老虎幼蟲與小地老虎成蟲頭部樣轉(zhuǎn)錄組差異基因5 768個(gè),小地老虎幼蟲與小地老虎成蟲胸部樣轉(zhuǎn)錄組差異基因4 836個(gè),小地老虎幼蟲與小地老虎成蟲腹部樣轉(zhuǎn)錄組差異基因8 917個(gè)。小地老虎成蟲頭部與腹部差異基因最多,小地老虎成蟲頭部與胸部差異基因最少。小地老虎幼蟲與成蟲頭、胸、腹三部分相比,小地老虎幼蟲與成蟲腹部的差異基因最多,與胸部的差異基因最少(圖3)。

    圖2 轉(zhuǎn)錄組不同數(shù)據(jù)庫注釋結(jié)果Venn圖

    圖3 各組樣本DESeq2與EdegR基因差異表達(dá)火山圖和差異基因Venn圖

    通過對小地老虎轉(zhuǎn)錄組KEGG通路富集分析,結(jié)合Interperscan結(jié)構(gòu)域預(yù)測和其他注釋信息得到小地老虎翅發(fā)育機(jī)制相關(guān)通路基因,并進(jìn)行表達(dá)量分析與表達(dá)量差異比較分析,對差異表達(dá)基因繪制小地老虎Wnt信號通路基因表達(dá)熱圖(圖4);并根據(jù)KEGG通路差異基因富集分析繪制小地老虎Wnt信號通路基因圖(圖5)。由小地老虎Wnt相關(guān)基因KEGG通路基因富集分析和表達(dá)量差異比較分析可知,Notum、Frizzled、CaN等53個(gè)元件均被注釋出,但是PRP、GBp、ICAT等17個(gè)元件沒有在小地老虎轉(zhuǎn)錄組中發(fā)現(xiàn)(圖4-5)。Wnt信號通路包含經(jīng)典Wnt通路,非經(jīng)典Wnt /鈣通路以及非經(jīng)典planar cell polarity通路,其中非經(jīng)典Wnt /鈣通路以及非經(jīng)典planar cell polarity通路總基因數(shù)目較少,但是在小地老虎轉(zhuǎn)錄組中注釋出的基因較多,非經(jīng)典Wnt/鈣通路只有一個(gè)元件未注釋出,非經(jīng)典planar cell polarity通路只有兩個(gè)元件未注釋出;非經(jīng)典Wnt /鈣通路,調(diào)控細(xì)胞內(nèi)的鈣,非經(jīng)典planar cell polarity通路調(diào)控細(xì)胞骨架,骨架決定細(xì)胞的形狀。

    圖4 小地老虎Wnt信號通路基因表達(dá)熱圖

    小地老虎胸部組織與頭部組織相比,F(xiàn)rizzled(Fz)、BMP and activin membrane-bound inhibitor(BAMBI)、(CKIα)、casein kinase 1(CKIε)、Axin、Myc proto-oncogene protein(c-myc)、E3 ubiquitinprotein ligase(Siah-1)、glypican 4(Knypek)、serine/threonine-protein phosphatase 2B catalytic subunit(CaN)、cullin 1(Cul1)等在小地老虎胸部組織中都上調(diào)表達(dá),Notum、calcyclin binding protein(SIP)、classical protein kinase C alpha type(PKC)在胸部組織中下調(diào)表達(dá)。小地老虎胸部組織與腹部組織相比,Notum、Dally、Siah-1、CaN等在小地老虎胸部組織中都上調(diào)表達(dá),wingless-type MMTV integration site family、Frizzled、Rac等在小地老虎胸部組織中下調(diào)表達(dá)。小地老虎胸部組織與小地老虎幼蟲組織相比,F(xiàn)rizzled、Dally、SIP、Cul1、Siah-1等上調(diào)表達(dá),Rac、phosphatidylinositol phospholipase C beta(PLC)、Frizzled等下調(diào)表達(dá)(圖5)。

    小地老虎胸部組織相對于頭部、腹部及幼蟲組織,只有Siah-1和CaN全部上調(diào)表達(dá),Notum雖然在小地老虎胸部組織相對于頭部、幼蟲對比中上調(diào)表達(dá),但在小地老虎胸部組織相對腹部組織中下調(diào)表達(dá)。SIP同時(shí)在小地老虎胸部相對于頭部和腹部比較中下調(diào)表達(dá)。

    非經(jīng)典planar cell polarity通路中,小地老虎胸部相對于頭部有表達(dá)差異的基因全部為表達(dá)上調(diào),小地老虎胸部相對于腹部有表達(dá)差異的基因全部為表達(dá)下調(diào)。非經(jīng)典planar cell polarity通路中,PLC在小地老虎胸部相對于頭部和幼蟲的基因表達(dá)中均為下調(diào),下游元件CaN表現(xiàn)為表達(dá)上調(diào)。

    圖5 小地老虎Wnt信號通路基因圖

    3 討 論

    Wnt信號通路相關(guān)基因在小地老虎幼蟲、小地老虎成蟲頭、胸、腹部均有表達(dá)。Wnt信號通路包含3種Wnt信號傳導(dǎo)通路:經(jīng)典Wnt通路,非經(jīng)典Wnt /鈣通路以及非經(jīng)典planar cell polarity通路。非經(jīng)典planar cell polarity通路是通過Wnt與Fz及其共受體結(jié)合而激活的。然后受體結(jié)合Dishevelled(Dsh),Dsh利用其PDZ和DIX結(jié)構(gòu)域與morphogenesis 1(Daam1)的Dishevelled-associated激活因子1形成復(fù)合物。然后Daam1通過鳥嘌呤交換因子激活G-蛋白Rho[24]。Rho激活Rho相關(guān)激酶Rho-associated kinase(ROCK),這是細(xì)胞骨架的主要調(diào)節(jié)因子之一。Dsh還與Ras-related C3 botulinum toxin substrate 1(Rac1)形成復(fù)合物,并介導(dǎo)profilin與肌動蛋白的結(jié)合[25]。Rac1激活JNK,也可導(dǎo)致肌動蛋白聚合。丙氨酸與肌動蛋白結(jié)合可導(dǎo)致細(xì)胞骨架和原腸胚的重組。經(jīng)典Wnt通路相關(guān)基因有一些在小地老虎轉(zhuǎn)錄組中未被注釋到,這可能與小地老虎的發(fā)育有關(guān),有些基因沒有表達(dá),也有可能與物種有關(guān),小地老虎可能不含有某些基因。同時(shí)這些在小地老虎腹部表達(dá)的基因可能參與一些基因轉(zhuǎn)錄的調(diào)控,這也是此信號通路的主要功能。

    Wnt信號通路相關(guān)基因在小地老虎幼蟲和成蟲不同組織表達(dá)差異揭示小地老虎可能的翅發(fā)育機(jī)制,非經(jīng)典planar cell polarity通路調(diào)控細(xì)胞骨架,決定骨架細(xì)胞的形狀,Rac在小地老虎腹部相對小地老虎胸部上調(diào)表達(dá),Rac也能間接影響JNK,從而影響一些基因的轉(zhuǎn)錄。非經(jīng)典Wnt /鈣通路中,小地老虎胸部相對于頭部、腹部和小地老虎幼蟲,CaN、Siah-1都表現(xiàn)為上調(diào)表達(dá),CaN在肌肉形成和功能的基本過程中起著重要的作用[26],因此小地老虎成蟲胸部CaN的表達(dá)可能與其肌肉功能有關(guān)。本次試驗(yàn)的小地老虎成蟲為處于遷飛狀態(tài)的小地老虎成蟲,胸部肌肉相關(guān)基因表達(dá)較高與其生物性狀相適應(yīng)。

    Wnt信號通路相關(guān)基因在小地老虎幼蟲、成蟲頭、胸、腹部的表達(dá)存在差異,小地老虎胸部組織與頭部組織、腹部組織和小地老虎幼蟲組織相比,CaN、Siah-1相關(guān)基因都顯著上調(diào)表達(dá)。這些基因有可能在其他昆蟲中也會有類似現(xiàn)象。對于昆蟲翅發(fā)育相關(guān)基因的研究主要集中于昆蟲羽化之前,本試驗(yàn)表明,在昆蟲羽化后一些翅發(fā)育相關(guān)基因會繼續(xù)在成蟲胸部表達(dá),也可能與翅部肌肉等組織的發(fā)育和代謝有關(guān)。這些對于昆蟲或者對于遷飛昆蟲的研究可能提供一些研究思路。近年來,隨著illumina、pacbio、nonapore等高通量測序技術(shù)的迅猛發(fā)展,以及CRISPR-Cas9技術(shù)的不斷發(fā)展,越來越多的科學(xué)問題通過基因組測序成功解決,也為研究昆蟲遷飛相關(guān)基因功能解析提供了新的技術(shù)方法[27]。

    另外,影響昆蟲翅發(fā)育相關(guān)通路還有Decapentaplegic(Dpp)、Hedgehog、Hippo信號通路等。Dpp是參與果蠅發(fā)育的關(guān)鍵形態(tài)發(fā)生素,Dpp也是第一個(gè)經(jīng)過驗(yàn)證的分泌形態(tài)發(fā)生素[12-13]。Hedgehog(Hh)信號通路在控制細(xì)胞增殖、組織模式形成、干細(xì)胞維持和發(fā)育方面發(fā)揮著多種作用[15-18]。Hippo信號通路也被稱為Salvador或者Warts通路,通過調(diào)控細(xì)胞增殖和凋亡來控制動物的器官大小。該通路的名字來源于它的一個(gè)關(guān)鍵信號元件——蛋白激酶Hippo(Hpo)[19-20]。該基因的突變導(dǎo)致組織過度生長,或類似“河馬”表型。

    昆蟲翅的進(jìn)化對昆蟲的意義重大,探究影響昆蟲翅生長發(fā)育的相關(guān)通路,有助于了解昆蟲翅的進(jìn)化及發(fā)育的過程。在眾多的昆蟲中,有相當(dāng)一部分為農(nóng)業(yè)害蟲,有些如棉鈴蟲、小地老虎等還是可進(jìn)行遷飛的害蟲,探究昆蟲翅的發(fā)育機(jī)制,可以為揭開昆蟲遷飛面紗并研發(fā)新型防治農(nóng)業(yè)害蟲方法提供幫助。本研究通過對小地老虎轉(zhuǎn)錄組的組裝及注釋,結(jié)合影響果蠅翅發(fā)育相關(guān)通路基因,找到許多潛在與小地老虎翅發(fā)育相關(guān)通路基因,為后續(xù)研究小地老虎翅發(fā)育相關(guān)基因提供了良好的基礎(chǔ)。

    4 結(jié)論

    通過對小地老虎轉(zhuǎn)錄組的組裝及注釋,最終組裝得到轉(zhuǎn)錄組243.4 Mb,組裝出轉(zhuǎn)錄本429 288個(gè),獨(dú)立基因37 744個(gè)。小地老虎Wnt相關(guān)通路中,Notum、Frizzled、CaN等53個(gè)相關(guān)元件被注釋出。Wnt信號通路中的非經(jīng)典planar cell polarity通路中,PlC在小地老虎胸部相對于頭部和幼蟲的基因表達(dá)中均為下調(diào),下游元件CaN、Siah-1在小地老虎胸部相對于頭部、腹部及幼蟲組織的基因表達(dá)中均表現(xiàn)為表達(dá)上調(diào),可能與其肌肉功能有關(guān)。

    猜你喜歡
    差異基因胸部成蟲
    ICR鼠肝和腎毒性損傷生物標(biāo)志物的篩選
    胸部腫瘤放療后椎體對99Tcm-MDP的攝取表現(xiàn)及分析
    栗山天牛成蟲取食、交配和產(chǎn)卵行為
    基于RNA 測序研究人參二醇對大鼠心血管內(nèi)皮細(xì)胞基因表達(dá)的影響 (正文見第26 頁)
    柿樹害蟲枯葉夜蛾防治效果試驗(yàn)研究
    天牛,就在樹上趴著
    人體胸部
    新誘劑對柑橘大實(shí)蠅成蟲的誘殺效果
    胸部Castleman病1例報(bào)道
    SSH技術(shù)在絲狀真菌功能基因篩選中的應(yīng)用
    天堂av国产一区二区熟女人妻| 精品乱码久久久久久99久播| 琪琪午夜伦伦电影理论片6080| 午夜亚洲福利在线播放| 男人和女人高潮做爰伦理| 国产精品美女特级片免费视频播放器| 人妻久久中文字幕网| 国产在线精品亚洲第一网站| 亚洲成人免费电影在线观看| 精品久久久久久,| 久久精品综合一区二区三区| 日韩 亚洲 欧美在线| 99热这里只有是精品50| 国内精品美女久久久久久| 九九久久精品国产亚洲av麻豆| 我要搜黄色片| 婷婷丁香在线五月| 免费在线观看日本一区| 一区二区三区高清视频在线| 国产私拍福利视频在线观看| 国内精品宾馆在线| 亚洲欧美日韩高清在线视频| 两个人的视频大全免费| a级毛片免费高清观看在线播放| 在线观看午夜福利视频| 一级黄片播放器| 日韩精品有码人妻一区| 国产毛片a区久久久久| 国产高清三级在线| 91麻豆精品激情在线观看国产| 无遮挡黄片免费观看| 国产午夜精品久久久久久一区二区三区 | 99热只有精品国产| 中文字幕熟女人妻在线| 成人二区视频| 看黄色毛片网站| 国产精品98久久久久久宅男小说| 国产男靠女视频免费网站| 国产高清视频在线观看网站| 国内毛片毛片毛片毛片毛片| 久久久久国内视频| 一个人观看的视频www高清免费观看| 久久久久久久久久久丰满 | 搞女人的毛片| 精品人妻偷拍中文字幕| 在线免费观看的www视频| 一级毛片久久久久久久久女| 好男人在线观看高清免费视频| 精品日产1卡2卡| 欧美高清成人免费视频www| 变态另类成人亚洲欧美熟女| 日韩人妻高清精品专区| 婷婷色综合大香蕉| 天堂动漫精品| 欧美在线一区亚洲| 国产视频内射| 在线看三级毛片| 中出人妻视频一区二区| 看片在线看免费视频| 欧美日韩黄片免| 99国产极品粉嫩在线观看| 国产精品人妻久久久影院| 少妇猛男粗大的猛烈进出视频 | 欧美日韩瑟瑟在线播放| 午夜福利在线在线| 色综合亚洲欧美另类图片| 国产精品久久久久久av不卡| 国产精品伦人一区二区| 看黄色毛片网站| 成人毛片a级毛片在线播放| 免费电影在线观看免费观看| 欧美zozozo另类| 老司机午夜福利在线观看视频| 免费av观看视频| 午夜福利高清视频| 国产精品一及| 国产精品精品国产色婷婷| 少妇猛男粗大的猛烈进出视频 | 国产精品久久久久久久电影| 村上凉子中文字幕在线| 国产精品,欧美在线| 最近视频中文字幕2019在线8| 毛片女人毛片| 日韩欧美在线二视频| netflix在线观看网站| 日韩中文字幕欧美一区二区| 午夜福利欧美成人| 99九九线精品视频在线观看视频| 国产精品一区www在线观看 | 88av欧美| 两个人视频免费观看高清| 熟女人妻精品中文字幕| 亚洲avbb在线观看| 波多野结衣巨乳人妻| 国产极品精品免费视频能看的| 国产男靠女视频免费网站| 99久久中文字幕三级久久日本| 69av精品久久久久久| 国产视频内射| 联通29元200g的流量卡| 久久久久久久久久成人| 小说图片视频综合网站| 91av网一区二区| 国产成年人精品一区二区| 中国美女看黄片| 成年女人毛片免费观看观看9| 亚洲午夜理论影院| 欧美+亚洲+日韩+国产| 精品国产三级普通话版| 香蕉av资源在线| 狂野欧美激情性xxxx在线观看| 国产亚洲精品久久久com| 国产亚洲av嫩草精品影院| x7x7x7水蜜桃| 欧美日韩亚洲国产一区二区在线观看| 亚洲人成网站高清观看| 国产黄片美女视频| 亚洲欧美日韩无卡精品| 国产精品一及| 日韩欧美三级三区| 哪里可以看免费的av片| а√天堂www在线а√下载| 亚洲欧美日韩高清专用| 精品人妻一区二区三区麻豆 | 欧美高清成人免费视频www| 国产视频一区二区在线看| 一区二区三区激情视频| 免费观看在线日韩| 午夜福利在线在线| 美女大奶头视频| 少妇的逼好多水| 中亚洲国语对白在线视频| av在线老鸭窝| 亚洲精品影视一区二区三区av| 国产一区二区在线av高清观看| 深爱激情五月婷婷| 午夜老司机福利剧场| 亚洲人与动物交配视频| av黄色大香蕉| 人妻制服诱惑在线中文字幕| 在线免费观看的www视频| 午夜精品久久久久久毛片777| 联通29元200g的流量卡| 亚洲av五月六月丁香网| 身体一侧抽搐| 欧美又色又爽又黄视频| 少妇被粗大猛烈的视频| АⅤ资源中文在线天堂| 国内毛片毛片毛片毛片毛片| 午夜精品在线福利| 亚洲精品乱码久久久v下载方式| 少妇熟女aⅴ在线视频| 久久婷婷人人爽人人干人人爱| 天美传媒精品一区二区| 人妻少妇偷人精品九色| 国产欧美日韩一区二区精品| 一级av片app| 成人美女网站在线观看视频| 国内精品一区二区在线观看| 久久久久久九九精品二区国产| 日本精品一区二区三区蜜桃| 免费看av在线观看网站| 亚洲一区二区三区色噜噜| 噜噜噜噜噜久久久久久91| 午夜激情福利司机影院| 欧美国产日韩亚洲一区| 麻豆一二三区av精品| 免费无遮挡裸体视频| 日韩欧美国产一区二区入口| 在线免费观看的www视频| 少妇裸体淫交视频免费看高清| 午夜激情欧美在线| av福利片在线观看| 国产黄a三级三级三级人| 97超级碰碰碰精品色视频在线观看| 国产精品不卡视频一区二区| 99视频精品全部免费 在线| 精品人妻偷拍中文字幕| 婷婷精品国产亚洲av在线| 亚洲av美国av| 91久久精品国产一区二区成人| 美女 人体艺术 gogo| 国产色婷婷99| 少妇被粗大猛烈的视频| 18禁黄网站禁片午夜丰满| 欧美日韩瑟瑟在线播放| 亚洲av免费在线观看| 日日干狠狠操夜夜爽| 99久久精品国产国产毛片| 国产69精品久久久久777片| 国产又黄又爽又无遮挡在线| 久久精品国产亚洲av涩爱 | 夜夜夜夜夜久久久久| 变态另类成人亚洲欧美熟女| 欧美性感艳星| 成人一区二区视频在线观看| 免费观看在线日韩| eeuss影院久久| 我的女老师完整版在线观看| 少妇裸体淫交视频免费看高清| 一进一出抽搐动态| 久久人妻av系列| 久久精品国产亚洲av涩爱 | 不卡一级毛片| 天天躁日日操中文字幕| 中文资源天堂在线| 亚洲精品在线观看二区| 国产欧美日韩精品亚洲av| 制服丝袜大香蕉在线| 亚洲va在线va天堂va国产| 99热这里只有精品一区| 窝窝影院91人妻| 日本免费一区二区三区高清不卡| 99久久中文字幕三级久久日本| 久久久色成人| 在线观看66精品国产| 九九热线精品视视频播放| 成年女人毛片免费观看观看9| 国产精品亚洲美女久久久| 国产午夜福利久久久久久| 亚洲成人精品中文字幕电影| 久久久久久久久大av| 一夜夜www| 简卡轻食公司| 免费观看精品视频网站| 久久久成人免费电影| 很黄的视频免费| 变态另类成人亚洲欧美熟女| а√天堂www在线а√下载| 国产 一区 欧美 日韩| 波野结衣二区三区在线| 亚洲美女黄片视频| 亚洲无线在线观看| 国内精品宾馆在线| 中国美白少妇内射xxxbb| www.色视频.com| 国产中年淑女户外野战色| 麻豆成人av在线观看| 极品教师在线免费播放| 老女人水多毛片| 国产精品福利在线免费观看| 亚洲欧美日韩高清专用| 国产探花极品一区二区| 国产 一区精品| 国产高清视频在线观看网站| 岛国在线免费视频观看| 欧美xxxx性猛交bbbb| 黄色日韩在线| 欧美一区二区亚洲| 少妇的逼水好多| 亚洲av.av天堂| 色视频www国产| 国产精品一区二区性色av| 免费人成视频x8x8入口观看| 免费av观看视频| 亚洲美女视频黄频| 嫁个100分男人电影在线观看| 校园春色视频在线观看| 波多野结衣巨乳人妻| 国产极品精品免费视频能看的| 黄色女人牲交| videossex国产| 精品99又大又爽又粗少妇毛片 | 国产伦精品一区二区三区视频9| 男女之事视频高清在线观看| 97超级碰碰碰精品色视频在线观看| 欧美日韩瑟瑟在线播放| 欧美日韩中文字幕国产精品一区二区三区| 精品一区二区三区av网在线观看| 国产精品一及| 别揉我奶头 嗯啊视频| 国产精品福利在线免费观看| 国产人妻一区二区三区在| 又紧又爽又黄一区二区| 超碰av人人做人人爽久久| 亚洲欧美精品综合久久99| 欧美精品啪啪一区二区三区| 午夜福利18| 毛片一级片免费看久久久久 | 色综合色国产| 亚洲av日韩精品久久久久久密| 一本精品99久久精品77| 九九久久精品国产亚洲av麻豆| 久久久国产成人精品二区| 欧美一区二区精品小视频在线| 一个人观看的视频www高清免费观看| 中文资源天堂在线| 国产精品av视频在线免费观看| 日日干狠狠操夜夜爽| 91麻豆av在线| 久久精品国产亚洲av涩爱 | 老熟妇乱子伦视频在线观看| а√天堂www在线а√下载| 午夜日韩欧美国产| 亚洲无线在线观看| 联通29元200g的流量卡| 亚洲国产欧洲综合997久久,| 亚洲人成网站高清观看| 国产午夜精品久久久久久一区二区三区 | 国产亚洲av嫩草精品影院| 国产色婷婷99| 99久久久亚洲精品蜜臀av| 五月玫瑰六月丁香| 人妻久久中文字幕网| 亚洲精品一区av在线观看| 午夜精品久久久久久毛片777| 全区人妻精品视频| 午夜老司机福利剧场| 少妇高潮的动态图| 亚洲国产精品成人综合色| 熟妇人妻久久中文字幕3abv| 超碰av人人做人人爽久久| 国产欧美日韩精品一区二区| 精品久久久久久,| 99久久无色码亚洲精品果冻| 国产视频内射| 中文字幕高清在线视频| 自拍偷自拍亚洲精品老妇| 麻豆成人午夜福利视频| 亚洲一区二区三区色噜噜| a在线观看视频网站| 看十八女毛片水多多多| 欧美成人a在线观看| 亚洲欧美日韩东京热| 免费看美女性在线毛片视频| 欧美人与善性xxx| 国产精品爽爽va在线观看网站| 久久国产精品人妻蜜桃| 中文字幕av成人在线电影| 免费观看人在逋| 久久精品国产亚洲网站| 久久精品国产亚洲av涩爱 | 国产精品电影一区二区三区| 天堂√8在线中文| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av香蕉五月| 色av中文字幕| 日本 av在线| 人妻久久中文字幕网| 美女大奶头视频| 成人性生交大片免费视频hd| 国产亚洲精品综合一区在线观看| av国产免费在线观看| 99riav亚洲国产免费| 天堂网av新在线| 亚洲美女黄片视频| 春色校园在线视频观看| 波多野结衣巨乳人妻| 久久亚洲精品不卡| 久久精品国产亚洲av天美| 五月伊人婷婷丁香| 亚洲avbb在线观看| 九九久久精品国产亚洲av麻豆| 亚洲男人的天堂狠狠| 欧美bdsm另类| 亚洲精品乱码久久久v下载方式| 国产三级在线视频| www.www免费av| 午夜日韩欧美国产| 一区福利在线观看| 日本欧美国产在线视频| 欧美一级a爱片免费观看看| 五月伊人婷婷丁香| 蜜桃久久精品国产亚洲av| 午夜a级毛片| 亚洲男人的天堂狠狠| 午夜a级毛片| 中文字幕免费在线视频6| 午夜福利欧美成人| 亚洲第一电影网av| 精品乱码久久久久久99久播| 亚洲在线观看片| 天堂动漫精品| 久久精品国产亚洲av天美| 嫩草影视91久久| 亚洲性夜色夜夜综合| 欧美中文日本在线观看视频| 久久久久精品国产欧美久久久| 大型黄色视频在线免费观看| 在线看三级毛片| 22中文网久久字幕| 成人性生交大片免费视频hd| 12—13女人毛片做爰片一| 亚洲欧美日韩卡通动漫| 中文资源天堂在线| 日本爱情动作片www.在线观看 | 国产高清三级在线| 村上凉子中文字幕在线| 国产精品久久视频播放| 欧美一级a爱片免费观看看| 亚洲一区二区三区色噜噜| 长腿黑丝高跟| 亚洲av不卡在线观看| 婷婷色综合大香蕉| 色尼玛亚洲综合影院| 国产免费av片在线观看野外av| 久久精品国产亚洲av天美| 久久婷婷人人爽人人干人人爱| 久久精品国产鲁丝片午夜精品 | 国产激情偷乱视频一区二区| 变态另类丝袜制服| av黄色大香蕉| 精品久久久久久,| 春色校园在线视频观看| 免费人成在线观看视频色| 97人妻精品一区二区三区麻豆| 成人无遮挡网站| 尾随美女入室| 国产一区二区三区av在线 | 亚洲七黄色美女视频| 久久精品国产鲁丝片午夜精品 | x7x7x7水蜜桃| 日韩欧美 国产精品| 国模一区二区三区四区视频| 一级毛片久久久久久久久女| 成年免费大片在线观看| 99久久九九国产精品国产免费| 久久久久久久久大av| 亚洲美女视频黄频| 少妇猛男粗大的猛烈进出视频 | 亚洲va日本ⅴa欧美va伊人久久| 美女xxoo啪啪120秒动态图| 联通29元200g的流量卡| 午夜福利在线在线| 国内少妇人妻偷人精品xxx网站| 黄色女人牲交| 亚洲欧美日韩高清专用| 不卡视频在线观看欧美| АⅤ资源中文在线天堂| 欧美国产日韩亚洲一区| 深夜精品福利| 18+在线观看网站| 免费观看人在逋| 白带黄色成豆腐渣| 99久久精品国产国产毛片| 成人av一区二区三区在线看| xxxwww97欧美| 少妇丰满av| 国产在视频线在精品| 99久久精品热视频| 男人狂女人下面高潮的视频| 精品人妻偷拍中文字幕| 亚洲黑人精品在线| 久久久久久久亚洲中文字幕| 午夜福利18| 国产成人影院久久av| 国产白丝娇喘喷水9色精品| 久久天躁狠狠躁夜夜2o2o| 亚洲第一区二区三区不卡| 久99久视频精品免费| av中文乱码字幕在线| 国产v大片淫在线免费观看| 久久婷婷人人爽人人干人人爱| 又黄又爽又刺激的免费视频.| 少妇的逼好多水| 在线观看舔阴道视频| 国产精品女同一区二区软件 | 婷婷六月久久综合丁香| 中出人妻视频一区二区| 嫩草影视91久久| 亚洲成人久久性| 国产亚洲91精品色在线| 免费看美女性在线毛片视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲乱码一区二区免费版| 精品一区二区三区视频在线| 精品不卡国产一区二区三区| 亚洲五月天丁香| 伦理电影大哥的女人| 观看美女的网站| 午夜日韩欧美国产| 精品久久国产蜜桃| 亚洲三级黄色毛片| 中文字幕高清在线视频| 少妇被粗大猛烈的视频| 日本撒尿小便嘘嘘汇集6| 色哟哟哟哟哟哟| 成人无遮挡网站| 欧美激情久久久久久爽电影| 99九九线精品视频在线观看视频| 国产在线男女| a级毛片免费高清观看在线播放| 搞女人的毛片| 嫁个100分男人电影在线观看| 大型黄色视频在线免费观看| 成人亚洲精品av一区二区| 久久久久久久久久黄片| 18+在线观看网站| 亚洲aⅴ乱码一区二区在线播放| av.在线天堂| 一级黄片播放器| 午夜福利高清视频| 国产精品一区www在线观看 | 亚洲av.av天堂| 日韩欧美在线二视频| 亚洲精品日韩av片在线观看| 在线看三级毛片| 91午夜精品亚洲一区二区三区 | 中文字幕av成人在线电影| 可以在线观看的亚洲视频| 黄色一级大片看看| 日本熟妇午夜| 俄罗斯特黄特色一大片| 欧美一级a爱片免费观看看| 少妇人妻一区二区三区视频| 国产成人影院久久av| 欧美高清成人免费视频www| 亚洲精品粉嫩美女一区| 精品久久久噜噜| 他把我摸到了高潮在线观看| 亚洲久久久久久中文字幕| 久久精品国产自在天天线| 色av中文字幕| 天堂√8在线中文| 伊人久久精品亚洲午夜| 在现免费观看毛片| 日本色播在线视频| 黄色欧美视频在线观看| 欧美日韩乱码在线| 动漫黄色视频在线观看| 成人鲁丝片一二三区免费| 99久久中文字幕三级久久日本| 日韩在线高清观看一区二区三区 | 九色国产91popny在线| 综合色av麻豆| av福利片在线观看| 很黄的视频免费| 色综合亚洲欧美另类图片| 亚洲中文日韩欧美视频| 久久天躁狠狠躁夜夜2o2o| 又粗又爽又猛毛片免费看| 亚洲美女搞黄在线观看 | 亚洲人成伊人成综合网2020| 丰满的人妻完整版| 国产精品不卡视频一区二区| a在线观看视频网站| av在线观看视频网站免费| 午夜久久久久精精品| 噜噜噜噜噜久久久久久91| 小说图片视频综合网站| 桃色一区二区三区在线观看| 免费不卡的大黄色大毛片视频在线观看 | 一个人看视频在线观看www免费| 禁无遮挡网站| 搡老岳熟女国产| 亚洲三级黄色毛片| 久久久久久久精品吃奶| 男女啪啪激烈高潮av片| 亚洲精品乱码久久久v下载方式| 色综合站精品国产| h日本视频在线播放| 国产中年淑女户外野战色| 99热网站在线观看| 国产午夜精品论理片| 毛片女人毛片| 3wmmmm亚洲av在线观看| 国产欧美日韩精品一区二区| 亚洲精品影视一区二区三区av| 此物有八面人人有两片| 久久久久久久久久久丰满 | 99热精品在线国产| 亚洲黑人精品在线| 神马国产精品三级电影在线观看| 老师上课跳d突然被开到最大视频| 国内毛片毛片毛片毛片毛片| 亚洲中文字幕一区二区三区有码在线看| 免费看日本二区| 亚洲熟妇熟女久久| 国产精品1区2区在线观看.| 久久久久性生活片| 国产精品,欧美在线| 欧美绝顶高潮抽搐喷水| 简卡轻食公司| 精品一区二区免费观看| 亚洲无线观看免费| 亚洲欧美日韩高清专用| 国内揄拍国产精品人妻在线| 国产亚洲精品综合一区在线观看| 男人的好看免费观看在线视频| 深夜a级毛片| 啦啦啦韩国在线观看视频| 性插视频无遮挡在线免费观看| av天堂在线播放| 精品久久久噜噜| 亚洲国产精品成人综合色| 国产三级中文精品| 真实男女啪啪啪动态图| 国产视频一区二区在线看| 成人性生交大片免费视频hd| 国产免费av片在线观看野外av| 国产午夜精品论理片| 久9热在线精品视频| 国内揄拍国产精品人妻在线| 91在线观看av| 麻豆国产av国片精品| 男女做爰动态图高潮gif福利片| 18+在线观看网站| 久久九九热精品免费| 日韩中字成人| 中文字幕av在线有码专区| 美女cb高潮喷水在线观看| 五月伊人婷婷丁香| 看黄色毛片网站| 久久国内精品自在自线图片| 男人舔女人下体高潮全视频| 在线播放无遮挡| 最近在线观看免费完整版| 婷婷六月久久综合丁香| 国产精品亚洲美女久久久| 国产人妻一区二区三区在| 91麻豆精品激情在线观看国产| 天天一区二区日本电影三级| 少妇人妻精品综合一区二区 | 成人毛片a级毛片在线播放| 久久久久精品国产欧美久久久|