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

    杧果RAV基因家族的全基因組分析

    2021-09-17 23:49:30孫宇劉志鑫葉子羅睿雄劉曉妹蒲金基張賀
    江蘇農(nóng)業(yè)學(xué)報 2021年4期
    關(guān)鍵詞:表達分析

    孫宇 劉志鑫 葉子 羅睿雄 劉曉妹 蒲金基 張賀

    摘要: 為了揭示RAV家族基因在杧果中的序列特征及其表達特性,采用生物信息學(xué)方法對杧果RAV家族基因進行序列分析,并通過qRT-PCR技術(shù)研究杧果膠孢炭疽菌和細菌性黑斑病菌侵染過程中該家族基因的相對表達量。結(jié)果表明,從杧果基因組中鑒定出6個RAV基因家族成員,其編碼的蛋白質(zhì)均具有AP2和B3超家族保守域結(jié)構(gòu),命名為MiRAV1~MiRAV6。MiRAV1和MiRAV5為不穩(wěn)定親水堿性蛋白質(zhì),MiRAV2和MiRAV3為不穩(wěn)定親水酸性蛋白質(zhì),MiRAV4為穩(wěn)定親水酸性蛋白質(zhì),MiRAV6為穩(wěn)定親水堿性蛋白質(zhì)。杧果RAV基因編碼的蛋白質(zhì)主要定位于細胞核中,無規(guī)則卷曲和α-螺旋是6個MiRAVs蛋白二級結(jié)構(gòu)的主要元件。系統(tǒng)進化樹結(jié)合基序分析結(jié)果表明,在杧果與擬南芥、煙草、蘋果、菠蘿和毛果楊中,RAV家族基因具有多樣性,在進化上結(jié)構(gòu)具有保守性。qRT-PCR結(jié)果顯示,膠孢炭疽菌侵染過程中,杧果RAV基因的相對表達量均顯著下調(diào);細菌性黑斑病菌侵染過程中,MiRAV1~ MiRAV6的相對表達量在3 h時均顯著下調(diào),MiRAV1、MiRAV2、MiRAV3和MiRAV6的相對表達量在6 h時均顯著上調(diào),MiRAV1、MiRAV5和MiRAV6的相對表達量分別在12 h、24 h、48 h、72 h時均顯著上調(diào)。以上發(fā)現(xiàn)將為研究杧果RAV基因家族成員的功能和作用機制奠定基礎(chǔ)。

    關(guān)鍵詞: 杧果;RAV基因家族;全基因組鑒定;表達分析

    中圖分類號: S667.7?? 文獻標識碼: A?? 文章編號: 1000-4440(2021)04-0957-11

    Genome-wide analysis of the RAV gene family in mango

    SUN Yu1,2, LIU Zhi-xin2, YE Zi1,2, LUO Rui-xiong3, LIU Xiao-mei1,2, PU Jin-ji1,2, ZHANG He1,2

    (1.College of Plant Protection, Hainan University, Haikou 570228, China;2.Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences/Key Laboratory of Integrated Pest Management on Tropical Grops, Ministry of Agriculture and Rural Affairs, Haikou 571101, China;3.Tropical Crops Genetic Resources Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101, China)

    Abstract: In order to reveal the sequence and expression characteristics of RAV gene in Mangifera indica, the sequence of the RAV gene family was analyzed using bioinformatics methods, and the relative expression level of RAV gene family was studied by qRT-PCR during the infection of Colletotrichum gloeosporioides and Xanthomonas campestris pv. mangiferaeindicae. The results showed that six RAV gene family members were identified from the mango genome, and the encoded proteins had the conserved domains of the AP2 and B3 superfamily, named MiRAV1-MiRAV6. MiRAV1 and MiRAV5 were unstable hydrophilic basic proteins, MiRAV2 and MiRAV3 were unstable hydrophilic acidic proteins, MiRAV4 was a stable hydrophilic acidic protein, and MiRAV6 was a stable hydrophilic basic protein. The proteins encoded by RAV gene of mango were mainly located in the nucleus, and random coil and α-helix were the main elements of secondary structure of six MiRAVs proteins. Phylogenetic tree construction combined with motif analysis showed that, among mango and Arabidopsis thaliana, tobacco, apple, pineapple and Populus tomentosa, the RAV family genes were diverse and their structure was conservative in evolution. The results of qRT-PCR indicated that the relative expression of RAV gene in mango was significantly down regulated during the infection of Colletotrichum gloeosporioides. During the infection of Xanthomonas campestris pv. mangiferaeindicae, the relative expression levels of MiRAV1-MiRAV6 were significantly down regulated at 3 h, and the relative expression levels of MiRAV1, MiRAV2, MiRAV3 and MiRAV6 were significantly up-regulated at 6 h, and the relative expression levels of MiRAV1, MiRAV5 and MiRAV6 were significantly up-regulated at 12 h, 24 h, 48 h and 72 h, respectively. These results will lay a foundation for further study on the function and mechanism of RAV gene family members in mango.

    Key words: mango;RAV gene family;genomic identification;expression analysis

    轉(zhuǎn)錄因子,又稱反式作用因子,其具有特殊的結(jié)構(gòu)可調(diào)控基因表達,是植物中重要的調(diào)控因子。植物轉(zhuǎn)錄因子在脅迫反應(yīng)中,對調(diào)控基因表達發(fā)揮重要作用[1-3],當植物受到逆境脅迫時,轉(zhuǎn)錄因子與相應(yīng)的順式作用元件結(jié)合,可以調(diào)控下游相關(guān)基因的表達 [4]。RAV(Related to AB13/VP1)轉(zhuǎn)錄因子廣泛存在于植物中,該家族同時具有2個保守結(jié)構(gòu)域[5-8]:位于N端,可與GCC-box區(qū)域特異性結(jié)合的AP2結(jié)構(gòu)域;位于C端,可與CACCTG序列特異性結(jié)合的B3結(jié)構(gòu)域[9]。VP1/ABI3以及ARFs轉(zhuǎn)錄因子家族都具有該結(jié)構(gòu)域,主要調(diào)控植物對生長素和脫落酸的反應(yīng)[10]。RAV屬于AP2/ERF家族的一個亞家族,該家族還包含ERF、AP2、DREB亞家族[11]。

    迄今為止,RAV基因家族成員已在擬南芥(Arabidopsis thaliana)[12-13]、黃瓜(Cucumis sativus)[14]、棉花(Gossypium)[15]、東方山羊豆(Galega orientalis)[16]、煙草(Nicotiana tabacum)[17]、茶樹(Camellia sinensis)[18]等植物中被陸續(xù)克隆和鑒定。例如:在擬南芥中發(fā)現(xiàn)7個RAV基因家族成員AtRAV1~AtRAV7,AtRAV1和AtRAV2在擬南芥中參與應(yīng)對外界脅迫和葉片衰老 [12-13];黃瓜中CsRAV在各個組織中均有表達[14];棉花中GhRAV表達受黃萎病菌的誘導(dǎo)[15];東方山羊豆中GoRAV表達不僅受低溫和聚乙二醇的誘導(dǎo),還受茉莉酸甲酯、水楊酸和脫落酸的誘導(dǎo)[16];煙草中GmRAV過量表達會導(dǎo)致葉片失綠,同時抑制植株的生長 [17];茶樹中CsRAV表達受NaCl、乙烯和低溫的誘導(dǎo)[18]。RAV基因廣泛的分布表明它們在植物的生長發(fā)育和抗逆反應(yīng)調(diào)控過程中起著非常重要的作用。

    杧果(Mangifera indica),被稱為熱帶水果之王,是熱帶地區(qū)重要的果樹,也是農(nóng)民脫貧致富的“金杧果”和“致富樹”。隨著2020年杧果全基因組測序的完成 [19],杧果基因功能的研究進程明顯加快,然而目前國內(nèi)外對杧果RAV基因家族的研究還未見報道。為此,本研究對杧果RAV基因家族成員的序列特征及其表達特性進行分析,為杧果基因功能的分析以及杧果的遺傳育種等方面提供參考。

    1 材料與方法

    1.1 杧果生物脅迫處理

    杧果嫁接苗由農(nóng)業(yè)農(nóng)村部儋州杧果種質(zhì)資源圃提供,品種為貴妃杧。將株高約60 cm,均勻一致,樹齡1年的杧果嫁接苗移栽到花盆內(nèi),以室溫25 ℃,相對濕度70%~90%,光周期12 h(光)/12 h(暗)的條件下進行生物脅迫處理。設(shè)置膠孢炭疽菌和細菌性黑斑病菌 2個處理,每個處理3株苗,處理0 h作為對照。以分生孢子含量為1 ml 2×106 個的膠孢炭疽菌分生孢子懸浮液、含量為2×107CFU 的細菌性黑斑病菌懸浮液分別對杧果嫁接苗葉片均勻噴霧,取樣時間分別為0 h、3 h、6 h、12 h、24 h、48 h、72 h,將杧果葉片剪碎,液氮速凍,-80 ℃保存待用。

    1.2 杧果全基因組數(shù)據(jù)的來源

    杧果全基因組數(shù)據(jù)來源于NCBI(PRJNA487154),擬南芥、煙草、蘋果(Malus domestica)、菠蘿(Ananas comosus)、毛果楊(Populus trichocarpa)5個物種的基因組和RAV基因信息分別來源于NCBI(https://www.ncbi.nlm.nih.gov/)和PlantTFDB(http://planttfdb.gao-lab.org/)數(shù)據(jù)庫。

    1.3 杧果RAV基因家族成員的鑒定

    從擬南芥數(shù)據(jù)庫獲取AtRAV1~AtRAV7氨基酸序列,并將其作為查詢對象在杧果基因組數(shù)據(jù)庫中執(zhí)行Blastp比較(E≤10-10)。使用Pfam工具(http://pfam.xfam.org/)檢測得到的杧果蛋白質(zhì)氨基酸序列結(jié)構(gòu)域,進一步去除沒有典型AP2和B3結(jié)構(gòu)域的蛋白質(zhì)氨基酸序列,最終得到所有杧果RAV基因家族成員[20]。其他5個物種中的RAV基因家族成員采用類似的方法進行篩選。

    1.4 杧果RAV基因家族的生物信息學(xué)分析

    通過在線工具ProtParam預(yù)測RAV相對分子質(zhì)量、等電點、外顯子數(shù)、脂溶指數(shù)、氨基酸數(shù)等理化性質(zhì)[21];利用PSORT在線軟件進行亞細胞定位預(yù)測;通過SOPMA在線軟件預(yù)測蛋白質(zhì)的二級結(jié)構(gòu);利用NCBI Conserved Domain Search預(yù)測RAV的保守結(jié)構(gòu)域,再通過SWISS-MODEL在線工具對蛋白質(zhì)保守域的三級結(jié)構(gòu)進行預(yù)測;利用MEME在線軟件對杧果RAV家族的蛋白質(zhì)保守基序進行分析[22],再通過DNAMAN 6.0軟件進行高同源蛋白的氨基酸序列比對。

    1.5 系統(tǒng)進化樹分析

    利用ClustalW程序?qū)x果和擬南芥、煙草、蘋果、菠蘿、毛果楊的RAV家族成員的氨基酸序列進行多序列比對。獲得RAV蛋白氨基酸序列后使用MEGA 7.0軟件,通過鄰近法構(gòu)建進化樹,其中校驗參數(shù)Bootstrap值設(shè)置為1 000次重復(fù)[23]。

    1.6 RNA提取及cDNA合成

    杧果葉片總RNA提取參照RNAprep Pure多糖多酚植物總RNA提取試劑盒(天根生化科技有限公司產(chǎn)品)說明書中的方法。利用RevertAid First Strand cDNA Synthesis Kit試劑盒反轉(zhuǎn)錄成第1鏈cDNA [24]。

    1.7 杧果RAV基因家族表達分析

    利用 QuantStudio 6 Flex的實時熒光定量PCR檢測系統(tǒng)檢測膠孢炭疽菌和細菌性黑斑病菌侵染下杧果RAV家族基因的表達量,引物見表1。反應(yīng)體系為18.0 μl,設(shè)置2個復(fù)孔。反應(yīng)程序參照UltraSYBR Mixture試劑盒說明書(北京康為世紀生物科技有限公司產(chǎn)品)進行設(shè)置。以杧果MiActin作為內(nèi)參基因,0 h的表達量為對照,運用2-△△Ct法進行數(shù)據(jù)統(tǒng)計,確定MiRAV1~MiRAV6的相對表達量[25],并通過HemI 1.0.3.7軟件繪制熱圖,其中設(shè)置膠孢炭疽菌相對表達量熱圖參數(shù)(條形圖的數(shù)量)為11,小數(shù)點位數(shù)值為2,最小值為-20,階躍值為2;設(shè)置細菌性黑斑病菌相對表達量熱圖參數(shù)(條形圖的數(shù)量)為11,小數(shù)點位數(shù)值為2,最小值為-2,階躍值為2。整個表達矩陣每個值進行對數(shù)轉(zhuǎn)化,底數(shù)為2。

    2 結(jié)果與分析

    2.1 全基因組水平鑒定杧果RAV基因家族成員

    在杧果基因組數(shù)據(jù)庫中進行Blastp檢索,發(fā)現(xiàn)杧果RAV基因家族有6個成員,分別命名為MiRAV1~MiRAV6。通過生物學(xué)在線分析工具ProtParam對杧果RAV基因家族6個成員MiRAV1~MiRAV6編碼的蛋白質(zhì)氨基酸序列的理化性質(zhì)進行了分析,結(jié)果(表2)表明,6個MiRAVs蛋白相對分子質(zhì)量為32 523 340~44 122 280,蛋白質(zhì)等電點為5.83~9.34,外顯子數(shù)均為2,脂溶指數(shù)為58.93~74.31,氨基酸數(shù)為281~392個,不穩(wěn)定指數(shù)為34.14~46.61,親水性平均數(shù)為-0.506~-0.863。因此,預(yù)測這6個MiRAVs蛋白中MiRAV1和MiRAV5為不穩(wěn)定親水堿性蛋白質(zhì),MiRAV2~MiRAV3為不穩(wěn)定親水酸性蛋白質(zhì),MiRAV4為穩(wěn)定親水酸性蛋白質(zhì),MiRAV6為穩(wěn)定親水堿性蛋白質(zhì)。

    對杧果RAV家族成員進行亞細胞定位預(yù)測,結(jié)果(表3)顯示,它們主要定位于細胞核中,其次定位于葉綠體和細胞質(zhì)中,推斷該家族成員主要存在于進行光合作用的器官中。其中,只有MiRAV1和MiRAV4沒有定位于葉綠體,MiRAV4和MiRAV6沒有定位于細胞質(zhì),線粒體、過氧化氫酶體,質(zhì)膜定位的MiRAVs成員最少。

    通過SOPMA在線軟件對MiRAV1~MiRAV6蛋白進行二級結(jié)構(gòu)預(yù)測。如圖1所示,6個MiRAVs蛋白無規(guī)則卷曲結(jié)構(gòu)占比為44.58%~54.85%,α-螺旋結(jié)構(gòu)占比為21.88%~31.67%,延伸鏈結(jié)構(gòu)占比為16.86%~19.13%,β-轉(zhuǎn)角結(jié)構(gòu)占比為4.85%~7.56%。因此,預(yù)測無規(guī)則卷曲和α-螺旋是6個MiRAVs蛋白二級結(jié)構(gòu)的主要元件,其次是延伸鏈和β-轉(zhuǎn)角。

    通過NCBI Conserved Domain Search對杧果RAV家族蛋白質(zhì)氨基酸序列進行保守結(jié)構(gòu)域分析,發(fā)現(xiàn)MiRAV1~MiRAV6蛋白均含有AP2和B3超家族保守域(圖2)。6個MiRAVs蛋白AP2超家族保守域位置介于18~124;B3超家族保守域位置介于137~291。基于同源建模原理,通過SWISS-MODEL軟件對杧果RAV蛋白的2個保守結(jié)構(gòu)域(AP2、B3)三級結(jié)構(gòu)進行預(yù)測,發(fā)現(xiàn)在AP2保守結(jié)構(gòu)域中,含有1個 α-螺旋和3個反向平行的β-轉(zhuǎn)角;B3保守結(jié)構(gòu)域中,含有1個α-螺旋、2個反向平行的β-折疊和延伸鏈。

    2.2 杧果RAV蛋白系統(tǒng)進化樹分析

    為了更好地了解杧果RAV基因家族的系統(tǒng)發(fā)育關(guān)系,將杧果(6個)、擬南芥(7個)、煙草(5個)、蘋果(11個)、菠蘿(2個)、毛果楊(5個)RAV基因編碼的蛋白質(zhì)氨基酸序列通過Blast工具進行蛋白質(zhì)氨基酸序列分析,利用ClustalW程序和MEGA 7.0軟件構(gòu)建6個杧果RAV蛋白與其他物種(共5種)植物RAV蛋白系統(tǒng)進化樹(圖3),將36個RAV成員分成3個分支(Class Ⅰ、Class Ⅱ和Class Ⅲ),結(jié)果顯示,6個物種的RAV蛋白家族成員相互嵌合在一起,表明植物RAV在進化過程中相對保守。每個分支中均存在2個杧果RAV成員,其中,Class Ⅰ包含杧果2個、擬南芥2個、毛果楊3個、蘋果5個和煙草1個;Class Ⅱ包含杧果2個、蘋果4個、菠蘿2個和擬南芥5個;Class Ⅲ包含杧果2個、毛果楊2個、煙草4個和蘋果2個。在36個RAV成員中,MiRAV1、MiRAV6分別與MDP0000128924、MDP0000939633親緣關(guān)系最近,MiRAV2、MiRAV3分別與AtRAV1、AtRAV7親緣關(guān)系最近,MiRAV4與MDP0000207722、MDP0000485280親緣關(guān)系最近,MiRAV5與Potri.018G109200.1、Potri.018G109200.2親緣關(guān)系最近。因此,推測RAV家族蛋白質(zhì)不僅具有多樣性,對于親緣關(guān)系較近的蛋白質(zhì),它們還具有相似或相近的生物學(xué)功能。

    2.3 杧果RAV家族蛋白質(zhì)motif和保守結(jié)構(gòu)域分析

    基于在線軟件MEME對杧果、擬南芥、煙草、蘋果、菠蘿、毛果楊RAV轉(zhuǎn)錄因子家族的蛋白質(zhì)保守基序進行分析(圖4),在RAV轉(zhuǎn)錄因子家族中鑒定出5個保守基序,主要的保守基序有4個(motif1~motif4),N端均存在保守的motif1,C端均存在保守的motif2。其中,MiRAV4~MiRAV5、AtRAV3~AtRAV4、MDP0000526584、MDP0000534780、MDP0000165802、MDP0000207722、MDP0000485280、Potri.018G109200.1、Potri.018G109200.2、Potri.003G212800.1、XP_016492322.1不含有motif5。通過基序分析發(fā)現(xiàn),RAV基因家族在進化上結(jié)構(gòu)具有保守性。

    通過DNAMAN 6.0軟件對杧果6個RAV蛋白和擬南芥7個、煙草5個、蘋果11個、菠蘿2個、毛果楊5個RAV蛋白進行多序列比對。結(jié)果(圖5)發(fā)現(xiàn),MiRAV1~MiRAV6與5個物種的RAV轉(zhuǎn)錄因子家族存在2個相似的保守域A和保守域B,保守域A與motif1 基本重疊,保守域B與 motif2、motif3 基本重疊。因此,可以推測motif1、motif2和motif3是構(gòu)成AP2、B3結(jié)構(gòu)域的主要結(jié)構(gòu)。

    2.4 杧果RAV基因家族的脅迫響應(yīng)

    以杧果MiActin作為內(nèi)參基因,0 h處理作為對照,利用qRT-PCR檢測并分析了2種病原菌侵染過程中杧果RAV家族基因的表達程度。qRT-PCR結(jié)果顯示,MiRAV1~MiRAV6和內(nèi)參基因MiActin的熔解曲線為單峰,熔解溫度分別為81.69 ℃、74.15 ℃、71.79 ℃、73.75 ℃、77.79 ℃、80.52 ℃、83.73 ℃,將qRT-PCR產(chǎn)物經(jīng)2%的瓊脂糖凝膠進行電泳,結(jié)果(圖6)顯示,目標條帶單一,清晰,與預(yù)計大小一致,說明所設(shè)計的引物具有特異性;杧果膠孢炭疽菌侵染過程中(圖7A),MiRAV1~MiRAV6的相對表達量均顯著下調(diào),其中MiRAV1~MiRAV6相對表達量在6 h時達到最低值,MiRAV1~MiRAV5的相對表達量在72 h時達到最高值,MiRAV6的相對表達量在48 h時達到最高值;杧果細菌性黑斑病菌侵染過程中(圖7B),MiRAV1~ MiRAV6的相對表達量在3 h時均顯著下調(diào),MiRAV1、MiRAV2、MiRAV3和 MiRAV6的相對表達量在6 h時顯著上調(diào),MiRAV1、MiRAV5和 MiRAV6的相對表達量分別在12 h、24 h、48 h、72 h時均顯著上調(diào)。結(jié)果表明,在不同病原菌侵染過程中,MiRAV1~ MiRAV6的相對表達量與對照差異顯著。

    3 討論

    基因組測序技術(shù)的提高和分子生物學(xué)研究的深入,為植物基因家族的鑒定、功能基因的挖掘提供了基礎(chǔ),大量植物全基因組已經(jīng)被成功測序。RAV家族是一類植物特異的轉(zhuǎn)錄因子,廣泛存在于植物細胞內(nèi)。在擬南芥[12-13]、黃瓜[14]、棉花[15]、東方山羊豆[16]、煙草[17]、茶樹[18]等中關(guān)于RAV已開展相應(yīng)研究,主要體現(xiàn)在參與應(yīng)對外界脅迫和葉片衰老等方面功能,但關(guān)于其杧果RAV轉(zhuǎn)錄因子家族的研究至今仍未見報道。

    本研究在杧果全基因組水平鑒定到6個RAV基因,理化性質(zhì)分析結(jié)果表明,MiRAV1和MiRAV5為不穩(wěn)定親水堿性蛋白質(zhì),MiRAV2和MiRAV3為不穩(wěn)定親水酸性蛋白質(zhì),MiRAV4為穩(wěn)定親水酸性蛋白質(zhì),MiRAV6為穩(wěn)定親水堿性蛋白質(zhì)。杧果RAV蛋白主要存在于細胞核中,無規(guī)則卷曲和α-螺旋是6個MiRAVs蛋白二級結(jié)構(gòu)的主要元件,均具有AP2和B3超家族保守域結(jié)構(gòu),AP2保守域三級結(jié)構(gòu)含有1個 α-螺旋和3個反向平行的β-轉(zhuǎn)角,B3保守域三級結(jié)構(gòu)含有1個 α-螺旋、2個反向平行的 β-折疊和延伸鏈,據(jù)此推斷它們在調(diào)控生物性狀上存在協(xié)同作用。杧果與擬南芥、煙草、蘋果、菠蘿、毛果楊RAV蛋白家族構(gòu)建系統(tǒng)進化樹結(jié)合基序分析結(jié)果顯示,6個物種的RAV蛋白家族成員相互嵌合在一起,表明植物RAV基因在進化過程中相對保守;在模式植物擬南芥中參與應(yīng)對外界脅迫和葉片衰老的基因AtRAV1和AtRAV2被聚類到ClassⅡ中,因此,推斷該ClassⅡ分支下的杧果MiRAV2和MiRAV3基因可能與應(yīng)對外界脅迫和葉片衰老有關(guān)[12-13];對于親緣關(guān)系較近的蛋白質(zhì),推斷它們還具有相似或相近的生物學(xué)功能;杧果RAV蛋白主要的保守基序有4個(motif1~ motif 4),N端均存在保守的motif1,C端均存在保守的motif2,motif1~ motif 3與盧合均等[15]研究棉花時所鑒定到的motif1~motif3基本一致,motif1與韓林賀等[26]研究普通煙草時所鑒定到的motif1基本一致;杧果與其他5個物種RAV蛋白進行多序列比對結(jié)果顯示,存在2個明顯的保守域A和保守域B,保守域A與motif1基本重疊,保守域B與motif2、motif3基本重疊,因此,可以推測motif1是構(gòu)成AP2的主要結(jié)構(gòu),motif2和motif3是構(gòu)成B3的主要結(jié)構(gòu)。

    中國杧果病害危害嚴重,常見的病害有炭疽病、細菌性黑斑病、白粉病、蒂腐病、瘡痂病等[27-28],其中,杧果膠孢炭疽病與細菌性黑斑病為害杧果最為嚴重[29-30]。為了解杧果RAV基因在膠孢炭疽菌與細菌性黑斑病菌侵染下相對表達量的變化,對杧果葉片進行膠孢炭疽菌與細菌性黑斑病菌生物脅迫處理,并通過qRT-PCR技術(shù)研究膠孢炭疽菌與細菌性黑斑病菌侵染過程中杧果RAV家族基因的相對表達量。結(jié)果顯示,在不同病原菌侵染過程中,MiRAV1~MiRAV6的相對表達量與對照差異顯著,表明MiRAV1~MiRAV6對生物脅迫有響應(yīng)。RAV基因家族的表達分析研究,更多集中在MeJA[16]、NaCl[18]、ABA和低溫[26]等非生物脅迫上,也有前人在研究中提到RAV基因?qū)γ藁S萎病菌[15]、番茄青枯病菌[31]等生物脅迫的響應(yīng),本研究通過對杧果接種膠孢炭疽菌與細菌性黑斑病菌后的分析,發(fā)現(xiàn)了此家族的部分成員在膠孢炭疽菌與細菌性黑斑病菌侵染植株之后出現(xiàn)上調(diào)表達的現(xiàn)象, 表明RAV基因可能在杧果提高抗炭疽病和細菌性黑斑病的能力方面起著積極作用,如MiRAV1、MiRAV5、MiRAV6在細菌性黑斑病菌侵染后的12~72 h顯著上調(diào)表達,可以作為后續(xù)研究杧果抗細菌性黑斑病機制的候選基因,這也在分子水平上拓寬了杧果響應(yīng)生物脅迫的機制。

    參考文獻:

    [1] OHAMA N, SATO H, SHINOZAKI K, et al. Transcriptional regulatory network of plant heat stress response[J]. Trends in Plant Science, 2017, 22(1): 53-65.

    [2] KIM Y S, AN C, PARK S, et al. CAMTA-mediated regulation of salicylic acid immunity pathway genes in arabidopsis exposed to low temperature and pathogen infection[J]. The Plant Cell, 2017, 29(10): 2465-2477.

    [3] PRASAD K V S K, ABDRL-HAMEED A A E, XING D, et al. Global gene expression analysis using RNA-seq uncovered a new role for SR1/CAMTA3 transcription factor in salt stress[J]. Scientific Reports, 2016, 6(1): 443-448.

    [4] SCHWECHHEIMER C, BEVAN M. The regulation of transcription factor activity in plants[J]. Trend in Plant Science, 1998, 3(10): 278-283.

    [5] GIRAUDAT J, HAUGE B M, VALON C, et al. Isolation of the Arabidopsis ABI3 gene by positional cloning[J]. The Plant Cell, 1992, 4(10): 1251-1261.

    [6] SUZUKI M, KAO C Y, MCCARTY D R. The conserved B3 domain of VIVIPAROUS1 has a cooperative DNA binding activity[J]. The Plant Cell Online, 1997, 9(5): 799-807.

    [7] RIECHMANN J L, LAIN S, GARCIA J A. Highlights and prospects of potyvirus molecular biology[J]. Journal of General Virology, 1992, 73(1): 1-16.

    [8] KAGAYA Y, OHMIYA K, HATTORI T. RAV1, a novel DNA-binding protein, binds to bipartiterecognition sequence through two distinct DNA-binding domains uniquely found in higher plants[J]. Nucleic Acids Res, 1999, 27(2): 470.

    [9] RIECHMANN J L, HEARD J, MARTIN G, et al. Arabidopsis transcription factors:genome-widecomparative analysis among eukaryotes[J]. Science, 2000, 290(10): 2105.

    [10]FINKELSTEIN R R, WANG M L, LYNCH T J, et al. The arabidopsis abscisic acid response locusABI4 encodes an APETALA2 domain protein[J]. Plant Cell, 1998, 10(6): 1043.

    [11]孫 堯,孫 鑫,王 雷. 毛果楊RAV/PLC基因生物信息學(xué)分析[J].廣東農(nóng)業(yè)科學(xué), 2019, 46(12): 36-41,153.

    [12]HU Y X, WANG Y X, LIU X F, et al. Arabidopsis RAV1 is down-regulated by brassinosteroid and may act as a negative regulator during plant development[J].Cell Research, 2004, 14(1): 8-15.

    [13]SOHN K H, LEE S C, JUNG H W, et al. Expression and functional roles of the pepper pathogen-induced transcription factor RAV1 in bacterial disease resistance,and drought and salt stress tolerance[J].Plant Molecular Biology, 2006, 61(6): 897-915.

    [14]蔣倫偉,胡麗芳,袁永成,等. 黃瓜RAV基因家族的全基因組分析[J].中國蔬菜, 2012, 1(16): 15-21.

    [15]盧合均,周忠麗,陳浩東,等. 棉花RAV基因家族的全基因組分析[J].棉花學(xué)報, 2014, 26(6): 471-482.

    [16] CHEN X, WANG Z, WANG X, et al. Isolationand characterization of GoRAV, a novel geneencoding a RAV-type protein in Galegae orientalis[J].Genes Genet Syst, 2009, 84(2): 101.

    [17]SRINIVASAN C, LIU Z R, HEIDMANN I, et al. Heterologous expression of the BABY BOOM AP2/ERF transcription factor enhances the regeneration capacity of tobacco (Nicotiana tabacum L.)[J].Planta, 2007, 225: 341-351.

    [18]陳林波,李葉云,王 琴,等. 茶樹冷誘導(dǎo)基因RAV的克隆與表達特性分析[J].植物生理學(xué)通訊, 2010, 46(4): 354.

    [19]PENG W, YING F L, JIAN F H, et al. The genome evolution and domestication of tropical fruit mango[J]. Genome Biology, 2020, 21(1). Doi:10.1186/s13059-020-01959-8.

    [20]FINN R D, COGGILL P, EBERHARDT R Y, et al. The pfam protein families database: towards a more sustainable future[J].Nucleic Acids Research, 2016, 44:279-285.

    [21]GASTEIGER E, GATTIKER A, HOOGLAND C, et al. ExPASy: the proteomics server for in-depth protein knowledge and analysis[J].Nucleic Acids Research, 2003, 31(13): 3784-3788.

    [22]BAILEY T L, BODEN M, BUSKE F A, et al. MEME SUITE: tools for motif discovery and searching[J].Nucleic Acids Research, 2009, 37: 202-208.

    [23]KUMAR S, STECHER G, TAMUR K. MEGA7: molecular evolutionary genetics analysis Version 7.0 for bigger datasets [J].Molecular Biology and Evolution, 2016, 33(7): 1870-1874.

    [24]夏 楊,蘇初連,晁 駿,等. 菠蘿VOZ轉(zhuǎn)錄因子序列特征及其對非生物脅迫的響應(yīng)[J].西北植物學(xué)報, 2018, 38(7): 1228-1234.

    [25]LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2 [-Delta Delta C(T)] method [J].Methods-A Companion to Methods in Enzymology, 2001, 25(4): 402-408.

    [26]韓林賀,徐同慶,任昂彥,等.普通煙草RAV家族基因結(jié)構(gòu)與功能分析[J].基因組學(xué)與應(yīng)用生物學(xué), 2019, 38(4): 1658-1665.

    [27]蒲金基,張 賀,周文忠. 芒果病害綜合防治技術(shù)[J].中國熱帶農(nóng)業(yè), 2015(3): 38-42.

    [28]楊鄭州,黃柳芳,謝曉娜,等. 葉疫病菌侵染芒果后葉片細胞壁降解酶活性測定[J]. 江蘇農(nóng)業(yè)科學(xué),2019,47(10):138-141.

    [29]喻群芳,漆艷香,張輝強,等. 4種生防菌對芒果細菌性黑斑病的田間防效[J].中國熱帶農(nóng)業(yè), 2019(3): 23-25.

    [30]于海英,蘭建強,王曉燕,等. 芒果膠孢炭疽菌致病性的初步研究[C]//中國植物病理學(xué)會. 中國植物病理學(xué)會2012年學(xué)術(shù)年會論文集. 北京:中國農(nóng)業(yè)科學(xué)技術(shù)出版社, 2012: 210-213.

    [31]LI C W, SU R C, CHENG C P, et al. Tomato RAV transcription factor is a pivotal modulator involved in the AP2/EREBP-mediated defense pathway [J].Plant Physiology, 2011, 156(1):213-227.

    (責任編輯:陳海霞)

    收稿日期:2020-11-25

    基金項目:國家重點研發(fā)專項(2019YFD1000504);中國熱帶農(nóng)業(yè)科學(xué)院基本科研業(yè)務(wù)費專項(1630042017019)

    作者簡介:孫 宇(1995-),男,海南三亞人,碩士研究生,研究方向為植物保護。(E-mail)sunyuqp@qq.com

    通訊作者:張 賀,(E-mail)atzzhef@163.com;蒲金基,(E-mail)cataspjj@163.com

    猜你喜歡
    表達分析
    巴西橡膠樹CCaMK基因家族成員的鑒定和表達分析
    太子參玉米黃質(zhì)環(huán)氧化酶基因的克隆與表達分析
    木薯PHOR1基因的序列分析及其乙烯和茉莉酸甲酯誘導(dǎo)表達特性
    雷公藤貝殼杉烯酸氧化酶基因的全長cDNA克隆與表達分析
    紅花生育酚環(huán)化酶基因的克隆及表達分析
    三色堇DFR基因的克隆及表達分析
    香蕉抗壞血酸過氧化物酶基因的克隆及表達分析
    膠孢炭疽菌漆酶基因Lac2的序列特征與表達分析
    玉米紋枯病病菌y—谷氨酰轉(zhuǎn)肽酶基因克隆與表達分析
    信號分子與葉銹菌誘導(dǎo)下小麥病程相關(guān)蛋白1基因的表達分析
    91精品一卡2卡3卡4卡| 超色免费av| 免费看光身美女| 伊人久久精品亚洲午夜| 人妻制服诱惑在线中文字幕| 国产精品不卡视频一区二区| 精品人妻熟女毛片av久久网站| 伦精品一区二区三区| 99re6热这里在线精品视频| 99视频精品全部免费 在线| 欧美日韩视频精品一区| 欧美性感艳星| 丰满饥渴人妻一区二区三| 免费黄频网站在线观看国产| 在线天堂最新版资源| 欧美日韩在线观看h| 色哟哟·www| 日韩免费高清中文字幕av| 午夜精品国产一区二区电影| www.av在线官网国产| 亚洲精品,欧美精品| 日本猛色少妇xxxxx猛交久久| 大片免费播放器 马上看| 欧美亚洲日本最大视频资源| 看免费成人av毛片| 亚洲色图 男人天堂 中文字幕 | 亚洲国产精品专区欧美| 精品卡一卡二卡四卡免费| 美女视频免费永久观看网站| 国产成人freesex在线| 亚洲精品自拍成人| 国产视频首页在线观看| 99热6这里只有精品| 天天影视国产精品| 久久99一区二区三区| 亚洲av.av天堂| 国产成人免费无遮挡视频| 亚洲欧美日韩另类电影网站| 国产午夜精品一二区理论片| 2022亚洲国产成人精品| 一级毛片电影观看| 久久久久国产精品人妻一区二区| 国产不卡av网站在线观看| av网站免费在线观看视频| 亚洲一级一片aⅴ在线观看| 黄片播放在线免费| 国产一级毛片在线| 日韩熟女老妇一区二区性免费视频| 国产精品无大码| 亚洲精品国产色婷婷电影| 国产精品久久久久久精品电影小说| 久久久久久久久久久丰满| 国产欧美亚洲国产| 视频区图区小说| 午夜免费鲁丝| 22中文网久久字幕| 日本猛色少妇xxxxx猛交久久| 国产在视频线精品| 欧美一级a爱片免费观看看| 久久女婷五月综合色啪小说| 97超视频在线观看视频| 亚洲精品国产av成人精品| 人妻系列 视频| 欧美变态另类bdsm刘玥| 汤姆久久久久久久影院中文字幕| 国产精品女同一区二区软件| 在线亚洲精品国产二区图片欧美 | 制服人妻中文乱码| 秋霞在线观看毛片| 麻豆成人av视频| 大香蕉久久成人网| 校园人妻丝袜中文字幕| 久久久久久人妻| 国产成人freesex在线| 免费人妻精品一区二区三区视频| 黄色怎么调成土黄色| 麻豆乱淫一区二区| 人人妻人人爽人人添夜夜欢视频| 日韩中字成人| 亚洲高清免费不卡视频| 一个人免费看片子| 伊人亚洲综合成人网| 欧美精品国产亚洲| 九色亚洲精品在线播放| 亚洲国产毛片av蜜桃av| 天天躁夜夜躁狠狠久久av| 九色成人免费人妻av| 久久国产亚洲av麻豆专区| 国产无遮挡羞羞视频在线观看| 伦精品一区二区三区| 亚洲美女黄色视频免费看| 亚洲国产精品一区三区| 久久久久久人妻| 日韩熟女老妇一区二区性免费视频| 国产国语露脸激情在线看| 黄色毛片三级朝国网站| 久久精品国产亚洲av天美| 99国产综合亚洲精品| 桃花免费在线播放| 水蜜桃什么品种好| 国产探花极品一区二区| 日本wwww免费看| 日韩欧美精品免费久久| 国产免费又黄又爽又色| av网站免费在线观看视频| 久久久精品94久久精品| 热re99久久精品国产66热6| 亚洲性久久影院| 国产 一区精品| 国精品久久久久久国模美| 亚洲av不卡在线观看| 国产欧美亚洲国产| 亚洲美女视频黄频| 久久久久网色| 亚洲,欧美,日韩| 69精品国产乱码久久久| 久久人人爽av亚洲精品天堂| 国产免费视频播放在线视频| 人人妻人人添人人爽欧美一区卜| 午夜福利视频精品| 热99国产精品久久久久久7| 中文字幕制服av| 日韩精品免费视频一区二区三区 | 国产色爽女视频免费观看| 精品少妇黑人巨大在线播放| 亚洲av欧美aⅴ国产| 国产一区二区在线观看av| 99九九线精品视频在线观看视频| 精品一区二区三区视频在线| 精品人妻一区二区三区麻豆| 国产在线免费精品| 一级a做视频免费观看| 国产日韩欧美视频二区| 99精国产麻豆久久婷婷| 中文字幕最新亚洲高清| videos熟女内射| 纵有疾风起免费观看全集完整版| 校园人妻丝袜中文字幕| 精品一区二区三卡| 一区二区三区精品91| 日日摸夜夜添夜夜爱| 老女人水多毛片| av.在线天堂| 多毛熟女@视频| 免费看av在线观看网站| 老司机影院毛片| 大片电影免费在线观看免费| 在线观看www视频免费| 女人久久www免费人成看片| 七月丁香在线播放| 久久久精品区二区三区| 看十八女毛片水多多多| 免费av中文字幕在线| 亚洲av综合色区一区| 亚洲经典国产精华液单| 激情五月婷婷亚洲| 亚洲精品一区蜜桃| 日产精品乱码卡一卡2卡三| 一区二区三区四区激情视频| 秋霞伦理黄片| 亚洲国产av新网站| 一边摸一边做爽爽视频免费| 国产精品久久久久久久电影| 国产综合精华液| av线在线观看网站| 国产免费一级a男人的天堂| 99九九线精品视频在线观看视频| 亚洲欧美色中文字幕在线| 这个男人来自地球电影免费观看 | 五月开心婷婷网| 国产伦精品一区二区三区视频9| 国产成人免费无遮挡视频| 内地一区二区视频在线| 国产成人一区二区在线| 视频区图区小说| 国产老妇伦熟女老妇高清| 久久久久久久亚洲中文字幕| 最新的欧美精品一区二区| 欧美另类一区| 久久精品国产a三级三级三级| 人人妻人人澡人人爽人人夜夜| 国产国拍精品亚洲av在线观看| 成年人午夜在线观看视频| 亚洲不卡免费看| 国产高清不卡午夜福利| 亚洲欧洲国产日韩| 少妇人妻精品综合一区二区| 免费看光身美女| 精品久久蜜臀av无| 免费观看无遮挡的男女| 日本欧美视频一区| 亚洲av成人精品一二三区| 久久久久精品性色| 91精品一卡2卡3卡4卡| 亚洲欧洲国产日韩| 国产精品久久久久成人av| 精品一品国产午夜福利视频| 亚洲精品自拍成人| 日韩成人av中文字幕在线观看| 高清在线视频一区二区三区| 夜夜看夜夜爽夜夜摸| 久久影院123| 亚洲精品第二区| av视频免费观看在线观看| 欧美xxⅹ黑人| 一级毛片 在线播放| 99热全是精品| 成年人免费黄色播放视频| 大码成人一级视频| 成年人午夜在线观看视频| 亚洲av综合色区一区| 精品国产露脸久久av麻豆| 最近最新中文字幕免费大全7| 欧美 亚洲 国产 日韩一| 亚洲第一av免费看| av免费观看日本| 国产免费视频播放在线视频| av在线app专区| 免费高清在线观看视频在线观看| 国产黄片视频在线免费观看| 久热这里只有精品99| 天天操日日干夜夜撸| 黄色毛片三级朝国网站| 少妇精品久久久久久久| 亚洲欧美精品自产自拍| 国产高清有码在线观看视频| 亚洲美女搞黄在线观看| 亚洲精品乱码久久久久久按摩| 激情五月婷婷亚洲| 日日啪夜夜爽| 两个人的视频大全免费| 亚洲欧洲日产国产| 在线亚洲精品国产二区图片欧美 | 久久久国产一区二区| 女性被躁到高潮视频| 男男h啪啪无遮挡| 中文字幕人妻熟人妻熟丝袜美| 免费高清在线观看日韩| 久久97久久精品| 中文字幕亚洲精品专区| 亚洲精品亚洲一区二区| 久久青草综合色| 成年女人在线观看亚洲视频| 日产精品乱码卡一卡2卡三| 性色avwww在线观看| 丝袜喷水一区| 日韩,欧美,国产一区二区三区| 日韩伦理黄色片| 视频中文字幕在线观看| 精品久久久噜噜| 国产免费一级a男人的天堂| 国产av国产精品国产| 免费看av在线观看网站| 老司机影院毛片| 色视频在线一区二区三区| 水蜜桃什么品种好| 日韩在线高清观看一区二区三区| 国产男女超爽视频在线观看| 日本黄色日本黄色录像| 亚洲精品日韩av片在线观看| 国产永久视频网站| av免费在线看不卡| 你懂的网址亚洲精品在线观看| 看十八女毛片水多多多| 午夜免费鲁丝| 久久国产亚洲av麻豆专区| 性色av一级| 亚洲精品一二三| 69精品国产乱码久久久| a级毛片免费高清观看在线播放| 欧美人与性动交α欧美精品济南到 | 最后的刺客免费高清国语| 美女视频免费永久观看网站| 尾随美女入室| 国产高清不卡午夜福利| 国产极品粉嫩免费观看在线 | 一级毛片电影观看| 国产日韩欧美视频二区| 极品少妇高潮喷水抽搐| 国产毛片在线视频| 免费观看在线日韩| 七月丁香在线播放| 亚洲欧美日韩卡通动漫| 看十八女毛片水多多多| 久久久午夜欧美精品| 久久99精品国语久久久| 午夜福利网站1000一区二区三区| 久久久久久久久大av| 大香蕉97超碰在线| 久久久久久久国产电影| 国产成人一区二区在线| 国产极品天堂在线| 免费久久久久久久精品成人欧美视频 | av专区在线播放| 中国美白少妇内射xxxbb| 欧美日韩在线观看h| 青春草亚洲视频在线观看| 久久精品久久精品一区二区三区| 久久久久人妻精品一区果冻| 中文字幕精品免费在线观看视频 | 国产精品久久久久成人av| 飞空精品影院首页| 考比视频在线观看| 免费高清在线观看日韩| 欧美精品人与动牲交sv欧美| av播播在线观看一区| 啦啦啦在线观看免费高清www| 免费观看无遮挡的男女| 国产精品蜜桃在线观看| 肉色欧美久久久久久久蜜桃| 99热网站在线观看| 亚洲精品乱久久久久久| 婷婷色综合www| 亚洲第一av免费看| 九草在线视频观看| 一区二区三区精品91| 午夜激情久久久久久久| 日韩大片免费观看网站| 精品国产国语对白av| 亚洲av在线观看美女高潮| 寂寞人妻少妇视频99o| 女人精品久久久久毛片| 一区二区三区乱码不卡18| 国产乱来视频区| 日韩电影二区| 全区人妻精品视频| 免费播放大片免费观看视频在线观看| 国语对白做爰xxxⅹ性视频网站| 一级爰片在线观看| 91国产中文字幕| av在线播放精品| 国产免费又黄又爽又色| 久久精品久久久久久噜噜老黄| 一本一本综合久久| 国产成人精品久久久久久| 亚洲综合色网址| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品久久久久久婷婷小说| 欧美xxxx性猛交bbbb| 国产在线免费精品| 2018国产大陆天天弄谢| 80岁老熟妇乱子伦牲交| 日韩制服骚丝袜av| 老女人水多毛片| 欧美亚洲 丝袜 人妻 在线| 最近2019中文字幕mv第一页| 国产伦理片在线播放av一区| 超色免费av| 国产乱来视频区| 国产av精品麻豆| 国模一区二区三区四区视频| 简卡轻食公司| 成年av动漫网址| 黄片播放在线免费| 免费日韩欧美在线观看| 久久精品国产亚洲网站| 亚洲精品乱码久久久v下载方式| 欧美精品人与动牲交sv欧美| 日韩av免费高清视频| av在线观看视频网站免费| 99re6热这里在线精品视频| 边亲边吃奶的免费视频| 91午夜精品亚洲一区二区三区| 2018国产大陆天天弄谢| 日本av免费视频播放| 精品国产乱码久久久久久小说| 久久人人爽人人爽人人片va| 一本久久精品| 亚洲国产色片| av黄色大香蕉| 亚洲精品乱码久久久久久按摩| 亚洲综合色网址| 日韩亚洲欧美综合| 免费人成在线观看视频色| 亚洲av中文av极速乱| 十分钟在线观看高清视频www| 精品少妇久久久久久888优播| 欧美日韩亚洲高清精品| 天堂中文最新版在线下载| 少妇人妻久久综合中文| 午夜老司机福利剧场| 亚洲av不卡在线观看| 日本黄大片高清| 日韩强制内射视频| 国内精品宾馆在线| 99久久人妻综合| 日韩av不卡免费在线播放| 在线观看国产h片| 啦啦啦啦在线视频资源| 国产一区二区三区av在线| av卡一久久| 五月开心婷婷网| 欧美少妇被猛烈插入视频| 特大巨黑吊av在线直播| 最新的欧美精品一区二区| 91久久精品国产一区二区成人| 女人久久www免费人成看片| 国产亚洲最大av| 97超碰精品成人国产| 一区二区三区精品91| 校园人妻丝袜中文字幕| 免费av不卡在线播放| 97在线视频观看| 国产精品国产av在线观看| xxx大片免费视频| 美女中出高潮动态图| 亚洲国产av影院在线观看| 亚洲精品一区蜜桃| 欧美亚洲 丝袜 人妻 在线| 成人二区视频| 国产成人精品久久久久久| 国产老妇伦熟女老妇高清| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美三级亚洲精品| 天堂8中文在线网| 涩涩av久久男人的天堂| 国产一区亚洲一区在线观看| 亚洲精品乱码久久久v下载方式| 91国产中文字幕| 日日爽夜夜爽网站| 大又大粗又爽又黄少妇毛片口| 欧美日韩一区二区视频在线观看视频在线| 久久久久国产精品人妻一区二区| 亚洲av在线观看美女高潮| 亚洲欧美色中文字幕在线| 国产午夜精品久久久久久一区二区三区| 嫩草影院入口| 夜夜骑夜夜射夜夜干| 久久国产精品大桥未久av| 蜜臀久久99精品久久宅男| 中文字幕av电影在线播放| 欧美日韩亚洲高清精品| 免费观看av网站的网址| 一区二区三区四区激情视频| videos熟女内射| 日本黄大片高清| av国产久精品久网站免费入址| 久久久精品免费免费高清| 各种免费的搞黄视频| 国产亚洲一区二区精品| 性高湖久久久久久久久免费观看| 中文精品一卡2卡3卡4更新| 亚洲高清免费不卡视频| 久久99热6这里只有精品| 日韩欧美一区视频在线观看| 狂野欧美激情性bbbbbb| 亚洲av二区三区四区| 亚洲综合色网址| 夜夜爽夜夜爽视频| 国产精品国产三级国产av玫瑰| 亚洲丝袜综合中文字幕| 欧美精品国产亚洲| 亚洲精品国产av蜜桃| 街头女战士在线观看网站| 一边摸一边做爽爽视频免费| 人妻一区二区av| 亚洲欧美色中文字幕在线| 激情五月婷婷亚洲| 大又大粗又爽又黄少妇毛片口| 久久久久人妻精品一区果冻| 人妻 亚洲 视频| 人人妻人人添人人爽欧美一区卜| 秋霞伦理黄片| 啦啦啦中文免费视频观看日本| 青春草视频在线免费观看| 激情五月婷婷亚洲| 考比视频在线观看| 一二三四中文在线观看免费高清| 午夜视频国产福利| 国产黄片视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 国产极品粉嫩免费观看在线 | 夜夜爽夜夜爽视频| 亚洲欧美精品自产自拍| 免费观看性生交大片5| 日韩欧美一区视频在线观看| 久久久久久久久久久久大奶| 高清av免费在线| 十八禁高潮呻吟视频| 伦理电影大哥的女人| 亚洲美女视频黄频| 亚洲综合精品二区| 少妇的逼水好多| 午夜av观看不卡| 女性被躁到高潮视频| 午夜福利,免费看| 日韩大片免费观看网站| 久久久久国产精品人妻一区二区| 免费av不卡在线播放| 少妇被粗大猛烈的视频| 国产精品熟女久久久久浪| 日韩 亚洲 欧美在线| 伦精品一区二区三区| 亚洲av男天堂| 韩国av在线不卡| 亚洲性久久影院| 日产精品乱码卡一卡2卡三| av又黄又爽大尺度在线免费看| 在线亚洲精品国产二区图片欧美 | 精品一区二区免费观看| 桃花免费在线播放| 亚洲高清免费不卡视频| 五月开心婷婷网| 国产亚洲欧美精品永久| 国产精品一二三区在线看| 国产成人精品无人区| 黄片播放在线免费| av黄色大香蕉| 边亲边吃奶的免费视频| 大香蕉97超碰在线| 国产又色又爽无遮挡免| 亚洲天堂av无毛| a级毛色黄片| 伊人久久精品亚洲午夜| 日本午夜av视频| 亚洲精品乱久久久久久| 欧美日韩av久久| 国产成人午夜福利电影在线观看| 国产一区二区在线观看av| 韩国av在线不卡| 免费高清在线观看视频在线观看| 成人无遮挡网站| 亚洲国产毛片av蜜桃av| 一本—道久久a久久精品蜜桃钙片| 99热全是精品| 精品少妇黑人巨大在线播放| 久久ye,这里只有精品| 亚洲国产成人一精品久久久| 精品人妻偷拍中文字幕| a级毛片黄视频| 国产在线视频一区二区| 99久久中文字幕三级久久日本| 少妇丰满av| 99热网站在线观看| av国产久精品久网站免费入址| 亚洲人成网站在线观看播放| 大片电影免费在线观看免费| 亚洲欧美日韩另类电影网站| 成人亚洲精品一区在线观看| 综合色丁香网| 亚洲国产精品国产精品| 亚洲av中文av极速乱| 国内精品宾馆在线| 午夜福利,免费看| 久久人人爽人人爽人人片va| 一级爰片在线观看| 美女主播在线视频| 97精品久久久久久久久久精品| 久久久久国产网址| 欧美另类一区| 精品久久久噜噜| 免费观看性生交大片5| 女性生殖器流出的白浆| 97精品久久久久久久久久精品| 亚洲精品自拍成人| 丁香六月天网| 一级毛片电影观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产毛片av蜜桃av| 亚洲精品久久午夜乱码| 国产在线视频一区二区| 午夜激情av网站| 久久99精品国语久久久| 国产一区有黄有色的免费视频| 国产男女内射视频| 亚洲人成77777在线视频| 91久久精品国产一区二区成人| 日韩av在线免费看完整版不卡| 99热6这里只有精品| 国产精品国产三级国产av玫瑰| 日韩电影二区| 好男人视频免费观看在线| 国产欧美亚洲国产| 久久久久国产精品人妻一区二区| 日韩av不卡免费在线播放| 国产一级毛片在线| 亚洲美女视频黄频| 国产一区有黄有色的免费视频| 七月丁香在线播放| 久热这里只有精品99| 一级毛片我不卡| 插逼视频在线观看| 999精品在线视频| av播播在线观看一区| 男人添女人高潮全过程视频| 在线观看三级黄色| 久久久国产精品麻豆| 久久国产亚洲av麻豆专区| 亚洲精品久久久久久婷婷小说| 男人爽女人下面视频在线观看| 九九久久精品国产亚洲av麻豆| 亚洲精品一区蜜桃| 天美传媒精品一区二区| 亚洲欧美色中文字幕在线| 91久久精品国产一区二区成人| 伦理电影大哥的女人| 麻豆精品久久久久久蜜桃| 欧美精品国产亚洲| 91精品一卡2卡3卡4卡| 老司机影院成人| 日韩人妻高清精品专区| 高清在线视频一区二区三区| 午夜免费鲁丝| 大陆偷拍与自拍| 免费看av在线观看网站| 九九久久精品国产亚洲av麻豆| 亚洲,欧美,日韩| 亚洲精品久久午夜乱码| av.在线天堂| 九九爱精品视频在线观看| 各种免费的搞黄视频| 久久久久久久久久久免费av| 女性被躁到高潮视频| 欧美日韩精品成人综合77777| 国产亚洲精品第一综合不卡 |