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

    大興安嶺草蓯蓉全長(zhǎng)轉(zhuǎn)錄組測(cè)序及生物信息學(xué)分析

    2024-12-31 00:00:00楊秋雪沃佳美雪徐曉敏郭忠錄劉樹民高宏偉

    摘 要 草蓯蓉為東北道地藥材,具有補(bǔ)腎壯陽,潤(rùn)腸通便,止血的功效,草蓯蓉中含有草蓯蓉多糖、草蓯蓉環(huán)烯醚萜苷、草蓯蓉酸等活性成分,具有豐富的藥用價(jià)值。目前對(duì)草蓯蓉的研究多集中于有效成分,藥理性等方面,在轉(zhuǎn)錄組及功能基因方面的內(nèi)容較少。為進(jìn)一步了解大興安嶺草蓯蓉的轉(zhuǎn)錄組,豐富其遺傳信息利用高通量測(cè)序技術(shù)獲得大興安嶺草蓯蓉種子的轉(zhuǎn)錄組信息。利用轉(zhuǎn)錄組測(cè)序獲得Unigene序列,將Unigene與七大功能數(shù)據(jù)庫(kù)注釋進(jìn)行比對(duì)注釋、SSR分析、CDS預(yù)測(cè),得到草蓯蓉轉(zhuǎn)錄組的信息。結(jié)果顯示,共獲得18.87 G的CleanData,各樣本的有效數(shù)據(jù)量分布在5.4~7.06 Gb,Q30堿基分布在95.66%~97.45%,平均GC含量為49.35%。拼接出Unigene 57 799條,總長(zhǎng)度為48 308 661 bp,平均長(zhǎng)度為835.8 bp。大興安嶺草蓯蓉與油菜(Brassica napus)序列相似度最高,47 035個(gè)(81.38%)基因注釋到非冗余蛋白序列數(shù)據(jù)庫(kù)(NR), "33 653個(gè)(58.22%)基因注釋到蛋白序列數(shù)據(jù)庫(kù)(Swissprot),13 174個(gè)(22.79%)基因注釋到京都基因與基因組百科全書數(shù)據(jù)庫(kù)(KEGG),27 886個(gè)(48.25%)基因注釋到真核生物蛋白相鄰類的聚簇?cái)?shù)據(jù)庫(kù)(KOG), "42 506個(gè)(73.54%)基因注釋到直系同源蛋白分組比對(duì)數(shù)據(jù)庫(kù)(eggNOG),30 928個(gè)(53.51%)基因注釋到基因本體論數(shù)據(jù)庫(kù)(GO),29 642個(gè)(51.28%)基因注釋到蛋白家族數(shù)據(jù)庫(kù)(Pfam)。預(yù)測(cè)出SSRs 11 031個(gè),包含SSRs的Unigene有57 799條,包含大于1個(gè)SSRs的Unigene有1 886條,復(fù)合型SSRs943個(gè)。預(yù)測(cè)出 "53 033條CDS序列,其中數(shù)據(jù)庫(kù)比對(duì)方法預(yù)測(cè)出47 140條,ESTScan預(yù)測(cè)出5 893條。研究結(jié)果豐富了大興安嶺草蓯蓉的轉(zhuǎn)錄組數(shù)據(jù),為深入研究大興安嶺[JP3]草蓯蓉生物學(xué)特性及分子機(jī)制等提供參考,也有利于其資源開發(fā)和保護(hù)利用。[JP]

    關(guān)鍵詞 草蓯蓉;全長(zhǎng)轉(zhuǎn)錄組測(cè)序;大興安嶺;代謝通路;生物信息分析

    草蓯蓉(Bosehniakia rossica)為列當(dāng)科草蓯蓉屬植物,多分布于中國(guó)黑龍江、吉林、內(nèi)蒙古等地區(qū)。全草入藥,味甘咸,性溫,功效補(bǔ)腎壯陽,潤(rùn)腸通便,止血,用于腎虛陽萎、腰關(guān)節(jié)冷痛、便秘等,具有豐富的藥用價(jià)值[1-5]。草蓯蓉種子中含有草蓯蓉多糖、草蓯蓉環(huán)烯醚萜苷,草蓯蓉苷、草蓯蓉酸等活性成分[6-8],具有抗腫瘤、抗炎、抗氧化、保肝等作用[9-10]。

    董學(xué)花等[11]研究發(fā)現(xiàn),草蓯蓉環(huán)烯醚萜苷能抑制人肝癌HepG2細(xì)胞活力,降低細(xì)胞黏附率,還能顯著下調(diào)TGF-β1誘導(dǎo)后HepG2中MMP2、MMP7、MMP9、Snail、Slug蛋白表達(dá),說明草蓯蓉環(huán)烯醚萜苷具有抗腫瘤作用;李香丹等[12]對(duì)采用四氯化碳植物油混懸液建立肝纖維化模型大鼠,與模型組相比,草蓯蓉環(huán)烯醚萜苷給藥組能明顯降低肝組織中的ALT、AST、TBIL水平,減輕肝纖維化程度,說明草蓯蓉環(huán)烯醚萜苷對(duì)肝損傷有一定的保護(hù)作用;鄭峰等[13]研究發(fā)現(xiàn),經(jīng)草蓯蓉環(huán)烯醚萜苷灌胃后可使肝癌大鼠血清中的MDA、GST水平降低,增加SOD活性,說明草蓯蓉化學(xué)醚萜苷能通過抗氧化作用進(jìn)而發(fā)揮抗腫瘤作用;此外,草蓯蓉多糖能抑制炎癥因子的表達(dá),如IL-β、IL-6、TNF-α等,還能降低LPS誘導(dǎo)的RAW264.7巨噬細(xì)胞中COX-2、MyD88蛋白含量[14]。

    草蓯蓉最早記載于唐《新修本草》[15]。大興安嶺以淺山丘陵為主,屬中溫帶大陸性季風(fēng)氣候。草蓯蓉性耐寒冷,喜愛冷涼氣候,大興安嶺的主要生態(tài)因子都適宜草蓯蓉生長(zhǎng),故草蓯蓉成為大興安嶺道地藥材之一[16-18]。有關(guān)草蓯蓉分子生物學(xué)、轉(zhuǎn)錄組和功能基因的相關(guān)信息較少,在一定程度上制約了草蓯蓉研究的進(jìn)展。[JP3]因此利用高通量測(cè)序技術(shù)對(duì)草蓯蓉種子進(jìn)行分析,為擴(kuò)大草蓯蓉資源,進(jìn)一步探索其功能基因信息提供理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    大興安嶺草蓯蓉轉(zhuǎn)錄組測(cè)序材料于2022年9月在大興安嶺地區(qū)采集,采集多個(gè)植株生長(zhǎng)成熟的草蓯蓉種子,用錫紙包裹并快速放入液氮處理,置于-80 ℃保存?zhèn)溆谩?/p>

    1.2 大興安嶺草蓯蓉種子RNA抽提和文庫(kù)構(gòu)建

    用TRIzo(ThermoScientific,貨號(hào)A33251)提取總RNA,采用NanoDrop2000分光光度計(jì)(ThermoScientific, USA)檢測(cè)RNA純度和定量,并用Agilent2100生物分析儀(AgilentTechnologies, SantaClara, CA, USA)評(píng)估RNA完整性。然后根據(jù)說明書使[JP3]用VAHTS Universal V6RNA-seqLibraryPrep試劑盒(Vazyme-innovation in enzyme technology,USA)構(gòu)建轉(zhuǎn)錄組文庫(kù)。[JP]

    1.3 轉(zhuǎn)錄組測(cè)序與組裝

    文庫(kù)質(zhì)量經(jīng)Agilent 2100生物分析儀質(zhì)檢合格后,使用IlluminaNovaseq6000測(cè)序平臺(tái)進(jìn)行測(cè)序,生成150 bp的雙端序列。建庫(kù)測(cè)序先經(jīng)過種子總RNA提取與RNA的質(zhì)量檢測(cè),再對(duì)RNA進(jìn)行富集、片段化逆轉(zhuǎn)錄生成cDNA,最后PCR擴(kuò)增完成上級(jí)測(cè)序,此步驟及部分?jǐn)?shù)據(jù)分析由上海歐易生物醫(yī)學(xué)科技有限公司完成。高通量測(cè)序得到的原始圖像數(shù)據(jù)文件經(jīng)堿基識(shí)別(BaseCalling)分析轉(zhuǎn)化為原始測(cè)序序列,稱之為Raw Data或Raw Reads,結(jié)果以FASTQ文件格式存儲(chǔ),其中包含測(cè)序序列(reads)的序列信息以及其對(duì)應(yīng)的測(cè)序質(zhì)量信息。通過測(cè)序,得到大量的樣本雙端測(cè)序數(shù)據(jù)。采用fastp軟件對(duì)原始數(shù)據(jù)進(jìn)行質(zhì)量預(yù)處理,首先對(duì)原始數(shù)據(jù)進(jìn)行質(zhì)控并去除接頭,再?gòu)?′端及5′端以不同方式去除低質(zhì)量堿基,最終得到高質(zhì)量的cleanreads。用Trinity軟件paired-end的拼接方法得到Transcript序列,根據(jù)序列相似性以及長(zhǎng)度,挑選出最長(zhǎng)的序列作為Unigene;之后再利用CD-HIT軟件聚類去冗余得到最終的Unigene,以此作為后續(xù)分析的參考序列。

    1.4 Unigene功能注釋

    對(duì)NR、COG、KOG、GO、Swiss-Prot、eggNOG、KEGG數(shù)據(jù)庫(kù)注釋,采用diamond軟件進(jìn)行比對(duì),取elt;1e-5的注釋,篩選具有最高序列相似性的蛋白,從而得到功能注釋信息。對(duì)于Pfam數(shù)據(jù)庫(kù)注釋,采用HMMER軟件和蛋白家族模型比對(duì),從而篩選出得分最高的家族。

    1.5 CDS預(yù)測(cè)

    將Unigene序列與NR、SwissProt和KOG數(shù)據(jù)庫(kù)作對(duì)比,直到與所有蛋白數(shù)據(jù)庫(kù)對(duì)比完,剩余對(duì)比不上的Unigene的編碼區(qū)采用ESTScan(3.0.3) 軟件預(yù)測(cè),獲得其編碼區(qū)的核酸序列(序列方向5′-gt;3′)和氨基酸序列。

    1.6 SSR分析

    利用MISA進(jìn)行SSR預(yù)測(cè),然后利用Primer3將獲得的大興安嶺草蓯蓉種子SSR位點(diǎn)的序列以退火溫度60 ℃±2 ℃、引物長(zhǎng)度20 bp左右、PCR產(chǎn)物長(zhǎng)度100~300 bp 為條件進(jìn)行引物設(shè)計(jì)。

    2 結(jié)果與分析

    2.1 轉(zhuǎn)錄組測(cè)序和denovo組裝

    對(duì)大興安嶺草蓯蓉種子樣本進(jìn)行轉(zhuǎn)錄組測(cè)序,每個(gè)樣本重復(fù)3次,共獲得18.87 Gb的CleanData,各樣本的有效數(shù)據(jù)量分布在5.4~7.06 Gb,Q30堿基分布在95.66%~97.45%,平均GC含量為49.35%。上述數(shù)據(jù)說明獲得的轉(zhuǎn)錄組數(shù)據(jù)質(zhì)量較好,可用于后續(xù)的分析處理(表1)。

    使用Trinity將cleanreads進(jìn)行denovo組裝成轉(zhuǎn)錄本,去冗余后共獲得57 799條Unigene,總長(zhǎng)48 308 661 bp,平均長(zhǎng)度為835.8 bp,N50為 "1 075 bp;其中長(zhǎng)度大于等于500 bp的有35 255條序列,占比61.00%;15 752條序列長(zhǎng)度大于等于1 000 bp,占比27.25%,見圖1。

    2.2 Unigene功能注釋

    大興安嶺草蓯蓉組裝后的Unigene序列與7個(gè)常見數(shù)據(jù)庫(kù)比對(duì),見表2。其中注釋到Unigene最多的數(shù)據(jù)庫(kù)為NR數(shù)據(jù)庫(kù),有47 035條,占 "81.38%,KEGG數(shù)據(jù)庫(kù)被注釋到的Unigene最少,有13 174條,占22.79%。

    為了解大興安嶺草蓯蓉種子的相近物種情況,將其轉(zhuǎn)錄組序列與NR庫(kù)數(shù)據(jù)比對(duì),匹配較多的物種主要有油菜(Brassica napus)、芝麻(Sesamum indicum)、蕓苔屬植物(Brassica rapa)、獨(dú)腳金(Striga asiatica)、紫花風(fēng)鈴木(Handroanthus impetiginosus),分別占12.09%、12.00%、8.83%、6.49%、5.84%。結(jié)果顯示大興安嶺草蓯蓉種子與油菜(Brassica napus)具有的同源序列匹配度比例最高, 見圖2。

    2.3 KEGG功能注釋

    KEGG數(shù)據(jù)庫(kù)注釋結(jié)果顯示有關(guān)大興安嶺草蓯蓉糖類的合成與代謝的Unigene被注釋335條,有關(guān)萜類和多酮化合物代謝的Unigene有263條,這3類Unigene可能與大興安嶺草蓯蓉的活性成分的合成與代謝有關(guān),見圖3。根據(jù)KEGG數(shù)據(jù)庫(kù)比對(duì),篩選出參與大興安嶺草蓯蓉種子活性成分合成通路的 Unigene,包括132條Unigene涉及萜類骨架生物合成,129條Unigene涉及精氨酸生物合成,128條Unigene涉及苯丙氨酸、酪氨酸和色氨酸生物合成,90條Unigene涉及類固醇生物合成,87條Unigene涉及泛醌和萜類醌生物合成,85條Unigene涉及葉酸生物合成,64條Unigene涉及賴氨酸生物合成,60條Unigene涉及三萜、哌啶和吡啶生物合成,56條Unigene涉及異喹啉生物堿生物合成,40條Unigene涉及類胡蘿卜素生物合成,38條Unigene涉及單環(huán)類生物合成,30條Unigene涉及類黃酮生物合成,13條Unigene涉及二萜類生物合成,9條Unigene涉及單萜類生物合成。

    2.4 GO功能注釋

    對(duì)大興安嶺草蓯蓉轉(zhuǎn)錄組Unigene進(jìn)行GO功能注釋,共分類到BP、CC和MF" 3個(gè)GO類別的44個(gè)條目中。其中,在第一大類BP中,共包含21個(gè)條目,其中Cellularprocess(細(xì)胞過程)和Metabolicprocess(代謝過程)基因數(shù)量最多;在第二大類Cellularcomponent(細(xì)胞成分)中,共包含有13個(gè)條目,其中Cell(細(xì)胞)、Cellpart(細(xì)胞組分)的基因數(shù)量最多;第三大類MF中,共包含有10個(gè)條目,Binding(結(jié)合因子)和Catalyticactivity(催化因子活性)的基因數(shù)量最多,見圖4。

    2.5 KOG功能注釋

    KOG數(shù)據(jù)庫(kù)注釋了60 500個(gè)轉(zhuǎn)錄本,按照KOG功能分類可分為25個(gè)組,見圖5。其中,“一般功能預(yù)測(cè)”的轉(zhuǎn)錄本最多,共4 778個(gè);“蛋白質(zhì)翻譯后修飾、蛋白周轉(zhuǎn)、分子伴侶”的轉(zhuǎn)錄本有3 380個(gè);“翻譯、核糖體結(jié)構(gòu)與生物發(fā)生”的轉(zhuǎn)錄本有2 998個(gè)。此外,發(fā)現(xiàn)“次生代謝物生物合成、轉(zhuǎn)運(yùn)與代謝”等相關(guān)的轉(zhuǎn)錄本有806個(gè),這些轉(zhuǎn)錄本可能與大興安嶺草蓯蓉中有效成分的合成途徑有關(guān)。

    2.6 編碼序列(CDS)分析

    對(duì)大興安嶺草蓯蓉轉(zhuǎn)錄組所有Unigene的CDS序列進(jìn)行預(yù)測(cè),共預(yù)測(cè)出53 033條CDS序列,其中數(shù)據(jù)庫(kù)比對(duì)方法預(yù)測(cè)出47 140條。對(duì)預(yù)測(cè)的CDS序列長(zhǎng)度進(jìn)行整合,其中長(zhǎng)度≥500 nt的序列占總序列數(shù)的比例為48.32%,長(zhǎng)度≥ "1 000 nt的序列占總序列數(shù)的比例為18.69%。ESTScan預(yù)測(cè)出5 893條,對(duì)預(yù)測(cè)的CDS序列長(zhǎng)度進(jìn)行整合,其中長(zhǎng)度≥500 nt的序列占總序列數(shù)的比例為16.14%,長(zhǎng)度≥1 000 nt的序列占總序列數(shù)的比例為1.49%,見圖6和圖7。

    2.7 簡(jiǎn)單重復(fù)序列(SSR)分析

    根據(jù)組裝結(jié)果,使用MISA軟件對(duì)Unigene的SSR進(jìn)行檢測(cè),將SSR統(tǒng)計(jì)結(jié)果按重復(fù)類型分為6類,共11 031條SSR。其中最豐富的重復(fù)類型為單堿基重復(fù),共有5 281條(占47.87%);其次為雙堿基重復(fù)和三堿基重復(fù)類型,分別為 "2 998條(占27.18%)和2 497條(22.64%);其余為四、五和六堿基重復(fù)類型,共255條(占 "2.31%),見表3。根據(jù)引物設(shè)計(jì)條件獲得了 "10 088條符合條件的引物(表4),這將為后續(xù)草蓯蓉分子標(biāo)記的開發(fā)提供參考。

    3 討論與結(jié)論

    轉(zhuǎn)錄組分析能獲得植物的功能基因從而了解植物在生長(zhǎng)發(fā)育階段的調(diào)節(jié)機(jī)制,草蓯蓉轉(zhuǎn)錄組測(cè)序結(jié)果顯示,大興安嶺草蓯蓉種子共獲得 "18.87" Gb CleanData,經(jīng)過建庫(kù)測(cè)序與生物信息分析得知Q30達(dá)到97.45%,GC含量為 "49.35%,表明大興安嶺草蓯蓉種子轉(zhuǎn)錄組測(cè)序數(shù)據(jù)質(zhì)量?jī)?yōu)良。通過denovo組裝,共獲得57 799條轉(zhuǎn)錄本,其中序列長(zhǎng)度大于等于500 bp的有 "35 255條,平均長(zhǎng)度為835.80 bp,N50為1 075 bp,表明大興安嶺草蓯蓉轉(zhuǎn)錄組序列組裝質(zhì)量較高。序列比對(duì)顯示大興安嶺草蓯蓉種子與油菜(Brassica napus)和芝麻(Sesamumin dicum)親緣關(guān)系最為接近。通過KEGG功能注釋,其中有關(guān)大興安嶺草蓯蓉糖類的合成與代謝的Unigene被注釋335條,有關(guān)萜類和多酮化合物代謝的Unigene有263條,有關(guān)其他次生代謝物的生物合成的Unigene有342條。草蓯蓉發(fā)揮藥效的主要化學(xué)成分為草蓯蓉多糖等,糖類作為物質(zhì)代謝的碳骨架能為植物提供能量,研究表明多糖類化合物參與了細(xì)胞各種生命現(xiàn)象的調(diào)節(jié),具有多種生理活性[19],因此對(duì)該類Unigene的挖掘和分析具有重要的研究意義。萜類化合物廣泛分布于自然界,具有保護(hù)植物細(xì)胞膜,刺激植物產(chǎn)生內(nèi)源激素并且保護(hù)植物免受強(qiáng)光的傷害等作用[20],其中三萜皂苷不僅能為植物抵御病原體還具有多種藥理活性[21],轉(zhuǎn)錄組測(cè)序結(jié)果顯示有關(guān)萜類和多酮化合物代謝的Unigene有263條,這些Unigene可能與大興安嶺草蓯蓉種子中萜類成分的合成與代謝有關(guān)。目前,SSR標(biāo)記已在蘋果等不同抗性植物的分子遺傳學(xué)檢測(cè)中得到很好的應(yīng)用,本研究共檢測(cè)到11 031條SSR,其中最豐富的重復(fù)類型為單堿基重復(fù),共有5 281條,為資源利用與藥材鑒定提供參考。大興安嶺草蓯蓉種子在調(diào)節(jié)碳代謝、氨基酸代謝、能量代謝和脂質(zhì)代謝等方面的基因可能是在生長(zhǎng)發(fā)育時(shí)期應(yīng)對(duì)低溫脅迫的基礎(chǔ)。本研究對(duì)大興安嶺草蓯蓉種子進(jìn)行了轉(zhuǎn)錄組測(cè)序分析,獲得大量轉(zhuǎn)錄水平數(shù)據(jù),豐富了對(duì)草蓯蓉道地產(chǎn)區(qū)基因的認(rèn)識(shí), 為今后草蓯蓉功能基因的進(jìn)一步探索提供了一定的數(shù)據(jù)參考, 同時(shí)也為日后基因工程技術(shù)對(duì)草蓯蓉進(jìn)行遺傳改良奠定了理論基礎(chǔ)。

    參考文獻(xiàn) Reference:

    [1] 許長(zhǎng)青,許順貴,劉春光,等.草蓯蓉的藥理作用研究進(jìn)展[J].武警醫(yī)學(xué),2023,34(1):75-78.

    XU CH Q,XU SH G,LIU CH G,et al.Research progress on pharmacological effects of Bosehniakia rossica[J].Medical Journal of the Chinese People’s Armed Police Forces,2023,34(1):75-78.

    [2]王曉琴,黃 慧,李彩峰,等.列當(dāng)屬、肉蓯蓉屬與草蓯蓉屬藥用植物親緣學(xué)初探[J].中國(guó)中藥雜志,2018,43(23):4608-4621.

    WANG X Q,HUANG H,LI C F,et al.Pharmacophylogenetic study of genus Orobanche,Cistanche and Boschniakia[J].China Journal of Chinese Materia Medica,2018,43(23):4608-4621.

    [3]李彩峰,王曉琴,劉 勇,等.草蓯蓉化學(xué)成分及藥理活性研究進(jìn)展[J].中草藥,2014,45(7):1016-1023.

    LI C F,WANG X Q,LIU Y,et al.Advances in studies on chemical constituents of" Boschniakia rossica" and their pharmacological activities[J].Chinese Traditional and Herbal Drugs,2014,45(7):1016-1023.

    [4]馬茜茜,元海丹,葉 利,等.不同草蓯蓉提取物中化學(xué)成分及藥理作用的研究進(jìn)展[J].安徽農(nóng)業(yè)科學(xué),2015,43(7):11-14.

    MA Q Q,YUAN H D,YE L,et al.Research advances in chemical constituents and" pharmacological effects of several Boschniakia rossica extracts[J].Journal of Anhui Agricultural Sciences,2015,43(7):11-14.

    [5]劉 青.草蓯蓉多糖調(diào)控miR-138-5p對(duì)H2O2誘導(dǎo)H9c2心肌細(xì)胞凋亡和氧化損傷的影響[J].中西醫(yī)結(jié)合心腦血管病雜志,2022,20(16):2914-2919.

    LIU Q.Effects of Boschniakia rossica" polysaccharides on H2O2-induced apoptosis and oxidative stress damage in H9c2 cardiomyocyte by regulating miR-138-5p[J].Chinese Journal of Integrative Medicine on Cardio-Cerebrovascular Disease,2022,20(16):2914-2919.

    [6]陳 虹,黃明偉,程振東,等.草蓯蓉多糖下調(diào)miR-302a-3p表達(dá)對(duì)缺氧/復(fù)氧誘導(dǎo)的心肌細(xì)胞損傷的影響[J].中國(guó)藥理學(xué)通報(bào),2023,39(7):1241-1247.

    CHEN H,HUANG M W,CHENG ZH D,et al.Effect of" Boschniakia rossica polysaccharides on hypoxia/reoxygenation-induced cardiomyocyte damage by down-regulating miR-302a-3p expression[J].Chinese Pharmacological Bulletin,2023,39(7):1241-1247.

    [7]崔 靜,王帥帥,曹萌萌,等.草蓯蓉多糖提取物對(duì)喉癌細(xì)胞增殖、遷移及程序性死亡分子1/程序性死亡配體-1水平的研究分析[J].中國(guó)耳鼻咽喉頭頸外科,2021,28(8):492-496.

    CUI J,WANG SH SH,CAO M M,et al.Effect of Cistanche deserticola extract on proliferation,migration and PD-1P/D-L1 level of laryngeal carcinoma cells[J].Chinese Archives of Otolaryngology-Head and Neck Surgery,2021,28(8):492-496.

    [8]華 龍,王晨宇,姚坤厚,等.草蓯蓉多糖對(duì)人肝癌細(xì)胞侵襲遷移的影響及機(jī)制[J].解放軍醫(yī)學(xué)雜志,2021,46(8):763-770.

    HUA L,WANG CH Y,YAO K H,et al.Effect of Boschniakia rossica polysaccharides on the invasion and migration of human hepatocarcinoma cells and its mechanism[J].Medical Journal of Chinese People’s Liberation Army,2021,46(8):763-770.

    [9]張 天.草蓯蓉多糖對(duì)脂多糖誘導(dǎo)的小鼠J774A.1巨噬細(xì)胞炎癥反應(yīng)的抑制作用[D].吉林延邊:延邊大學(xué),2022.

    ZHANG T.Inhibitory effect of" Boschniakia rossica polysaccharides on lipopolysaccharide-lnduced inflammatory response in mouse J774A.1 macrophages[D].Yanbian Jilin:Yanbian University,2022.[ZK)]

    [10] 金愛花,常 慧,張 釗,等.草蓯蓉多糖對(duì)過氧化氫致HL-7702細(xì)胞損傷的保護(hù)作用[J].食品與發(fā)酵工業(yè),2021,47(21):132-139.

    JIN A H,CHANG H,ZHANG ZH,et al.Protective effect of Boschniakia rossica polysaccharides on hydrogen peroxide induced damage to hl-7702 cells[J].Food and Fermentation Industries,2021,47(21):132-139.

    [11]董學(xué)花,延光海,朱潔波,等.草蓯蓉環(huán)烯醚萜苷聯(lián)合5-氟尿嘧啶對(duì)人肝癌HepG2、SK-Hep1細(xì)胞上皮間質(zhì)轉(zhuǎn)化的抑制效果研究[J].中草藥,2020,51(17):4498-4505.

    DONG X H,YAN G H,ZHU J B,et al.Inhibitory effects of iridoid glycosides from Boschniakia rossica combined with 5-fluorouracil on human hepatoma HepG2 and SK-Hep1 cells epithelial-mesenchymal transition[J].Chinese Traditional and Herbal Drugs,2020,51(17):4498-4505.

    [12]李香丹,金永日,金海燕,等.草蓯蓉環(huán)烯醚萜苷對(duì)大鼠肝纖維化的保護(hù)作用及其機(jī)制[J].延邊大學(xué)醫(yī)學(xué)學(xué)報(bào),2016,39(3):175-178.

    LI X D,JIN Y R,JIN H Y,et al.Protective effects and the mechanisms of Boschniakia rossica iridoid glycoside on liver fibrosis in rats[J].Journal of Medical Science Yanbian University,2016,39(3):175-178.

    [13]鄭 峰,崔香丹,金雪峰,等.草蓯蓉環(huán)烯醚萜苷對(duì)肝癌大鼠抗氧化活性的影響[J].時(shí)珍國(guó)醫(yī)國(guó)藥,2017,28(4):787-789.

    ZHENG F,CUI X D,JIN X F,et al.Effects of the antioxidant activity for" iridoid glucosides from Boschniakia rossica on rat" liver cancer[J].Lishizhen Medicine and Materia Medica Research,2017,28(4):787-789.

    [14]劉莉園,張 釗,葛乃嘉,等.草蓯蓉多糖對(duì)脂多糖誘導(dǎo)的RAW264.7巨噬細(xì)胞炎癥反應(yīng)的影響[J].中國(guó)藥學(xué)雜志,2021,56(18):1479-1485.

    LIU L Y,ZHANG ZH,GE N J,et al.Effect of Boschniakia rossica" polysaccharides on LPS-induced inflammation of RAW264.7 macrophage[J].Chinese Pharmaceutical Journal,2021,56(18):1479-1485.

    [15]張 航,王 輝,楊樹東.草蓯蓉的本草考證及藥理研究[J].長(zhǎng)春中醫(yī)藥大學(xué)學(xué)報(bào),2016,32(1):40-42.

    ZHANG H,WANG H,YANG SH D.Herbalogical study and pharmacological effect of Bosehniakia rossic[J].Journal of Changchun University of Chinese Medicine,2016,32(1):40-42.

    [16]張峰源,蘇遠(yuǎn)航,周婉瑩,等.晝夜不對(duì)稱增溫對(duì)大興安嶺森林凈初級(jí)生產(chǎn)力的影響[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2023,43(6):137-147.

    ZHANG F Y,SU Y H,ZHOU W Y,et al.Effects of different changes in daytime and nighttime temperatures on forest net primary productivity in the Greater Khingan Mountains[J].Journal of Central South University of Forestry amp; Technology,2023,43(6):137-147.

    [17]姜 玉,金立成,黃巖飛.黑龍江省大興安嶺林區(qū)中醫(yī)藥產(chǎn)業(yè)發(fā)展研究[J].中國(guó)產(chǎn)經(jīng),2020(14):101-102.

    JIANG Y,JIN L CH,HUANG Y F.Research on the development of traditional Chinese medicine industry in the" Daxing’an mountains of Heilongjiang province[J].Chinese Industry amp; Economy,2020(14):101-102.

    [18]薩如拉,白學(xué)良.大興安嶺南部山地藥用苔蘚植物多樣性[J].分子植物育種,2019,17(18):6186-6197.

    SA R L,BAI X L.Diversity of medicinal bryophytes in the southern Greater Khingan Mountains[J].Molecular Plant Breeding,2019,17(18):6186-6197.

    [19]師仁增,鄧 霞,田 琳,等.植物根系分泌物對(duì)解磷微生物的影響研究進(jìn)展[J].山東農(nóng)業(yè)科學(xué),2023,55(8):174-180.

    SHI R Z,DENG X,TIAN L,et al.Research progress of" influences of" plant root exudates on" phosphorus solubilizing microorganisms[J].Shandong Agricultural Sciences,2023,55(8):174-180.

    [20]戴瑤瑤,閆濱濱,顏雨豪,等.玄參中環(huán)烯醚萜類化合物的研究進(jìn)展[J].中草藥,2023,54(9):2993-3003.

    DAI Y Y,YAN B B,YAN Y H,et al.Research advances on iridoids of Scrophularia ningpoensi[J].Chinese Traditional and Herbal Drugs,2023,54(9):2993-3003.

    [21]韋筱媚,廖振鋒,胡 營(yíng),等.崗梅全長(zhǎng)轉(zhuǎn)錄組測(cè)序及其三萜皂苷生物合成相關(guān)基因的挖掘[J].中國(guó)現(xiàn)代中藥,2023,25(7):1515-1528.

    WEI X M,LIAO ZH F,HU Y,et al.Full-length transcriptome sequencing of Ilex asprella (Hook.et Arn.) Champ.ex Benth.and exploration of genes involved in its triterpenoid saponins biosynthesis[J].Modern Chinese Medicine,2023,25(7):1515-1528.

    Full-length Transcriptome Sequencing and Bioinformatics Analysis of Boschniakia" rossica in Daxing’an Mountains

    Abstract Boschniakia rossica is a valuable medicinal materials native to northeastern China,known for its therapeutic effect in strengthening kidney function,combating constipation, and promoting hemostasis.It contains active compounds such as Boschniakia rossica polysaccharides, Iridoid glucosides,Boschniakia rossica acid, all contributing to its significant medicinal value.To further explore the transcriptome of stanche from Daxing’an Mountains and enhance its genetic database, we conducted a full-length transcriptomic analysis of Boschniakia rossica seeds using high-throughput sequencing technology.Through transcriptomes squencing we obtained Unigene sequence,followed by functional annotation, SSR (Simple Sequence Repeat) analysis and CDS (Coding Sequence) prediction, and an in-depth analysis of genetic information.The results showed that the a total of 18.87 Gb" of clean data was obtained, with the effective data per sample ranging from 5.4 to 7.06 Gb.The Q30 base quality score ranged from 95.66%" to 97.45%, and the average GC content was 49.35%.We identified 57 799 Unigen sequnces, with a total length of 48 308 661 bp and an average length of" "835.8 bp.The Boschniakia rossica seeds transcriptome showed the highest sequence similarity to Brassica napus.Functional annotations revealed that 47 035(81.38%) genes were annotated in the NR database," "33 653(58.22%) in Swissprot, 13 174(22.79%) in KEGG, 27 886(48.25%) in KOG, 42 506 "(73.54%) in eggNOG, 30 928(53.51%) in GO, 29 642(51.28%) in Pfam.11 031 SSRs were predicted, with" "57 799 Unigene containing SSRs, 1 886 Unigene containing more than 1SSRs, and compound SSRs 943 predicted.A total of 53 033 CDS sequences were predicted, with 47 140 were predicted by the database alignment and 5 893 by ESTScan.These findings enrich the transcriptome data of Boschniakia rossica, and provide a reference for further studying the biological characteristics and molecular mechanism of Boschniakia rossica.This study also lays a foundation for the sustainable exploitation, protection, and utilization of Boschniakia rossica resources.

    Key words Boschniakia rossica; Full-length transcriptome sequencing;Daxing’an Mountains;Metabolic pathway; Bioinformatics analysis

    netflix在线观看网站| 黄色毛片三级朝国网站| 久久精品aⅴ一区二区三区四区| 国产一级毛片七仙女欲春2 | 多毛熟女@视频| 色av中文字幕| 极品人妻少妇av视频| 免费在线观看亚洲国产| 精品欧美一区二区三区在线| www.熟女人妻精品国产| 黄色丝袜av网址大全| 亚洲人成伊人成综合网2020| www国产在线视频色| 亚洲国产中文字幕在线视频| √禁漫天堂资源中文www| 91大片在线观看| 久久人妻av系列| 一级毛片高清免费大全| 又黄又粗又硬又大视频| 亚洲精品一卡2卡三卡4卡5卡| 女同久久另类99精品国产91| 12—13女人毛片做爰片一| 无遮挡黄片免费观看| 国产成人欧美| 亚洲国产欧美日韩在线播放| 久久久久亚洲av毛片大全| 亚洲国产毛片av蜜桃av| 亚洲片人在线观看| 老汉色av国产亚洲站长工具| 亚洲性夜色夜夜综合| 大型av网站在线播放| 桃红色精品国产亚洲av| 最新在线观看一区二区三区| 亚洲avbb在线观看| 亚洲在线自拍视频| 窝窝影院91人妻| 91在线观看av| 一边摸一边抽搐一进一出视频| 国产亚洲精品久久久久久毛片| 男人舔女人的私密视频| 成人av一区二区三区在线看| 黑人巨大精品欧美一区二区蜜桃| 免费少妇av软件| 欧美激情 高清一区二区三区| 免费看美女性在线毛片视频| av超薄肉色丝袜交足视频| 亚洲色图 男人天堂 中文字幕| 免费在线观看亚洲国产| 熟女少妇亚洲综合色aaa.| 欧美久久黑人一区二区| 色综合婷婷激情| 9热在线视频观看99| 12—13女人毛片做爰片一| 波多野结衣av一区二区av| 亚洲第一青青草原| 男女床上黄色一级片免费看| 久久久水蜜桃国产精品网| cao死你这个sao货| 亚洲一卡2卡3卡4卡5卡精品中文| 十分钟在线观看高清视频www| 国内久久婷婷六月综合欲色啪| 可以在线观看毛片的网站| 在线十欧美十亚洲十日本专区| 激情视频va一区二区三区| 成人国产一区最新在线观看| 一级a爱片免费观看的视频| 国产91精品成人一区二区三区| 亚洲男人天堂网一区| 国产精品自产拍在线观看55亚洲| 免费搜索国产男女视频| 欧美人与性动交α欧美精品济南到| 亚洲成国产人片在线观看| 久久热在线av| 一本综合久久免费| 亚洲三区欧美一区| 久久人人精品亚洲av| 亚洲视频免费观看视频| av在线天堂中文字幕| 国产av在哪里看| 亚洲全国av大片| 久久国产精品男人的天堂亚洲| 国产精品一区二区三区四区久久 | 久久婷婷人人爽人人干人人爱 | 国产人伦9x9x在线观看| 精品久久久精品久久久| 成人特级黄色片久久久久久久| 亚洲av成人不卡在线观看播放网| 国产高清激情床上av| 午夜福利免费观看在线| 国产精华一区二区三区| 中出人妻视频一区二区| 精品第一国产精品| 亚洲va日本ⅴa欧美va伊人久久| 身体一侧抽搐| 久久久国产欧美日韩av| 岛国视频午夜一区免费看| 久久久久久国产a免费观看| 啪啪无遮挡十八禁网站| 日韩一卡2卡3卡4卡2021年| 亚洲第一电影网av| 亚洲午夜理论影院| 侵犯人妻中文字幕一二三四区| 身体一侧抽搐| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲男人天堂网一区| 国产高清有码在线观看视频 | 精品国产超薄肉色丝袜足j| 97人妻精品一区二区三区麻豆 | 男女下面进入的视频免费午夜 | 日韩欧美一区视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 1024香蕉在线观看| 日本免费a在线| 99在线视频只有这里精品首页| 法律面前人人平等表现在哪些方面| 亚洲欧美日韩高清在线视频| 岛国在线观看网站| 亚洲专区中文字幕在线| 国产精品久久久人人做人人爽| 亚洲人成77777在线视频| 在线观看免费日韩欧美大片| 欧美日韩亚洲国产一区二区在线观看| 一级毛片精品| 一区二区三区激情视频| 一级a爱片免费观看的视频| 久久久久久国产a免费观看| 久久久久九九精品影院| 可以免费在线观看a视频的电影网站| 国产99久久九九免费精品| 国语自产精品视频在线第100页| av天堂在线播放| 久久精品国产清高在天天线| 国产精品久久久久久人妻精品电影| 久久中文字幕人妻熟女| 后天国语完整版免费观看| 熟女少妇亚洲综合色aaa.| 在线观看舔阴道视频| 老司机靠b影院| 精品国产美女av久久久久小说| 亚洲九九香蕉| 国产精品美女特级片免费视频播放器 | 国产在线精品亚洲第一网站| 少妇熟女aⅴ在线视频| 91字幕亚洲| 亚洲第一青青草原| 久久久久久久久久久久大奶| 亚洲欧美激情综合另类| 窝窝影院91人妻| 在线观看www视频免费| 在线观看一区二区三区| 极品教师在线免费播放| 国产欧美日韩综合在线一区二区| 欧美久久黑人一区二区| 777久久人妻少妇嫩草av网站| av中文乱码字幕在线| 欧美黑人欧美精品刺激| 最近最新免费中文字幕在线| 首页视频小说图片口味搜索| 一级毛片高清免费大全| 丝袜在线中文字幕| 欧美老熟妇乱子伦牲交| 国产精品久久久久久人妻精品电影| 亚洲一区高清亚洲精品| 黑人欧美特级aaaaaa片| 少妇 在线观看| 日韩精品免费视频一区二区三区| 精品国产乱码久久久久久男人| 无人区码免费观看不卡| 看免费av毛片| 亚洲精品中文字幕在线视频| 一区福利在线观看| 狠狠狠狠99中文字幕| 最新在线观看一区二区三区| 99国产精品99久久久久| 少妇粗大呻吟视频| 亚洲,欧美精品.| 母亲3免费完整高清在线观看| 久久久国产成人精品二区| 久久 成人 亚洲| 身体一侧抽搐| 中文字幕人妻丝袜一区二区| 嫁个100分男人电影在线观看| 多毛熟女@视频| 国产精品国产高清国产av| 国产野战对白在线观看| 免费高清在线观看日韩| 无人区码免费观看不卡| 免费一级毛片在线播放高清视频 | 久99久视频精品免费| 可以在线观看的亚洲视频| 人妻丰满熟妇av一区二区三区| 欧美成人午夜精品| 一进一出好大好爽视频| 欧美成人一区二区免费高清观看 | 精品福利观看| 无遮挡黄片免费观看| 亚洲性夜色夜夜综合| 99久久国产精品久久久| 啦啦啦观看免费观看视频高清 | 日韩三级视频一区二区三区| 亚洲成av人片免费观看| 久久久久精品国产欧美久久久| 女性被躁到高潮视频| 国产一区二区三区视频了| 少妇被粗大的猛进出69影院| 一级,二级,三级黄色视频| 亚洲avbb在线观看| 亚洲欧美精品综合一区二区三区| 一区二区日韩欧美中文字幕| av网站免费在线观看视频| 人妻久久中文字幕网| 99久久99久久久精品蜜桃| 亚洲精品中文字幕在线视频| 国产又色又爽无遮挡免费看| 俄罗斯特黄特色一大片| 国产97色在线日韩免费| 欧美乱色亚洲激情| 美女大奶头视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩一级在线毛片| 啦啦啦韩国在线观看视频| 成人av一区二区三区在线看| 天堂√8在线中文| 精品国产一区二区久久| 99久久99久久久精品蜜桃| 天天躁夜夜躁狠狠躁躁| 少妇被粗大的猛进出69影院| 免费观看人在逋| 国产国语露脸激情在线看| 免费在线观看日本一区| 少妇粗大呻吟视频| 久久久久久久午夜电影| 精品一区二区三区四区五区乱码| 女人爽到高潮嗷嗷叫在线视频| 精品人妻1区二区| 俄罗斯特黄特色一大片| 少妇的丰满在线观看| 又黄又爽又免费观看的视频| 国产高清激情床上av| 视频在线观看一区二区三区| 黄色a级毛片大全视频| 99久久久亚洲精品蜜臀av| 日本免费a在线| 级片在线观看| 欧美不卡视频在线免费观看 | 国产精品久久久人人做人人爽| 视频在线观看一区二区三区| 亚洲av电影不卡..在线观看| 国产熟女xx| 国产av在哪里看| 搡老妇女老女人老熟妇| 精品欧美一区二区三区在线| 欧美日韩精品网址| 久久天堂一区二区三区四区| 啪啪无遮挡十八禁网站| 天堂影院成人在线观看| 精品一区二区三区视频在线观看免费| 欧美日本中文国产一区发布| 亚洲人成电影免费在线| 国产亚洲精品综合一区在线观看 | 亚洲色图综合在线观看| 涩涩av久久男人的天堂| 高清黄色对白视频在线免费看| 99精品在免费线老司机午夜| 国产成年人精品一区二区| 校园春色视频在线观看| 国产三级在线视频| 久久这里只有精品19| 亚洲人成网站在线播放欧美日韩| 亚洲五月天丁香| 国产成+人综合+亚洲专区| 19禁男女啪啪无遮挡网站| 日韩欧美国产在线观看| 久热这里只有精品99| 国产精品亚洲av一区麻豆| 午夜精品国产一区二区电影| 久久久国产成人免费| 亚洲五月色婷婷综合| 亚洲七黄色美女视频| 黄色毛片三级朝国网站| 免费在线观看完整版高清| 欧美另类亚洲清纯唯美| 人人妻,人人澡人人爽秒播| tocl精华| 日本黄色视频三级网站网址| 这个男人来自地球电影免费观看| 九色亚洲精品在线播放| 电影成人av| 搡老熟女国产l中国老女人| 啦啦啦 在线观看视频| 久久久久精品国产欧美久久久| 国产精品久久久av美女十八| 操美女的视频在线观看| 一区二区三区高清视频在线| 亚洲五月婷婷丁香| 亚洲人成网站在线播放欧美日韩| 亚洲精品久久国产高清桃花| 黄色成人免费大全| av超薄肉色丝袜交足视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品99久久99久久久不卡| 制服丝袜大香蕉在线| 中文字幕高清在线视频| videosex国产| 久久久久久久久中文| 欧美黑人欧美精品刺激| 久久伊人香网站| 亚洲国产高清在线一区二区三 | 亚洲 欧美 日韩 在线 免费| 长腿黑丝高跟| 欧美丝袜亚洲另类 | 国产视频一区二区在线看| 女性生殖器流出的白浆| www国产在线视频色| 又紧又爽又黄一区二区| 长腿黑丝高跟| 性欧美人与动物交配| 久久国产精品影院| 亚洲中文av在线| 欧美成人免费av一区二区三区| 又紧又爽又黄一区二区| 精品久久久久久久毛片微露脸| 人妻久久中文字幕网| 黄色女人牲交| 国产色视频综合| 热re99久久国产66热| 丁香欧美五月| 桃红色精品国产亚洲av| 窝窝影院91人妻| 久久国产精品男人的天堂亚洲| 亚洲五月婷婷丁香| 黑人操中国人逼视频| 久久婷婷人人爽人人干人人爱 | 日本免费一区二区三区高清不卡 | 色婷婷久久久亚洲欧美| 精品国产美女av久久久久小说| 国产极品粉嫩免费观看在线| 不卡av一区二区三区| 88av欧美| 久久久国产成人精品二区| 激情视频va一区二区三区| 久久久精品欧美日韩精品| 欧洲精品卡2卡3卡4卡5卡区| 亚洲一区二区三区色噜噜| 精品国产乱码久久久久久男人| 欧美大码av| 国产亚洲精品第一综合不卡| 免费看十八禁软件| 国产成人免费无遮挡视频| 日韩精品免费视频一区二区三区| 久9热在线精品视频| 亚洲va日本ⅴa欧美va伊人久久| 88av欧美| 精品久久蜜臀av无| 欧美日韩黄片免| 亚洲七黄色美女视频| 国产99久久九九免费精品| 男人舔女人下体高潮全视频| 在线观看免费午夜福利视频| 亚洲午夜精品一区,二区,三区| 亚洲熟妇熟女久久| 老司机福利观看| 99久久久亚洲精品蜜臀av| 国产成人精品无人区| 中文字幕av电影在线播放| 日韩欧美一区视频在线观看| 黄色片一级片一级黄色片| 亚洲一区中文字幕在线| 欧美激情极品国产一区二区三区| 老司机福利观看| www.自偷自拍.com| cao死你这个sao货| 成人欧美大片| 男人舔女人下体高潮全视频| 99香蕉大伊视频| 曰老女人黄片| 淫秽高清视频在线观看| 久久午夜综合久久蜜桃| 亚洲国产精品999在线| 欧美人与性动交α欧美精品济南到| 欧美激情久久久久久爽电影 | 免费观看精品视频网站| 午夜a级毛片| 热99re8久久精品国产| 黄网站色视频无遮挡免费观看| 免费女性裸体啪啪无遮挡网站| 看免费av毛片| 国产av一区二区精品久久| 亚洲免费av在线视频| 免费在线观看影片大全网站| 中亚洲国语对白在线视频| 一级毛片高清免费大全| 免费在线观看完整版高清| 一卡2卡三卡四卡精品乱码亚洲| 日本免费一区二区三区高清不卡 | 757午夜福利合集在线观看| 国产精品日韩av在线免费观看 | 黑人巨大精品欧美一区二区mp4| 亚洲无线在线观看| 级片在线观看| 美女免费视频网站| 男人的好看免费观看在线视频 | 麻豆av在线久日| 搞女人的毛片| 黄色视频不卡| 午夜久久久久精精品| 国产三级在线视频| 男人操女人黄网站| 亚洲色图av天堂| 成在线人永久免费视频| 精品第一国产精品| 在线观看一区二区三区| 日本vs欧美在线观看视频| 老汉色∧v一级毛片| 最近最新中文字幕大全免费视频| 成人欧美大片| 国产成年人精品一区二区| 高清在线国产一区| 嫩草影视91久久| 成年人黄色毛片网站| 在线十欧美十亚洲十日本专区| 久久久久久大精品| 久久热在线av| 久久精品aⅴ一区二区三区四区| 日本vs欧美在线观看视频| 99久久久亚洲精品蜜臀av| 久久天躁狠狠躁夜夜2o2o| 欧美大码av| 精品国产超薄肉色丝袜足j| 亚洲黑人精品在线| 亚洲一区二区三区不卡视频| 欧美日本中文国产一区发布| 午夜久久久久精精品| 日韩精品免费视频一区二区三区| 性欧美人与动物交配| 此物有八面人人有两片| 午夜精品在线福利| 亚洲熟女毛片儿| 日韩av在线大香蕉| 午夜福利影视在线免费观看| 欧美日本中文国产一区发布| 制服人妻中文乱码| 午夜a级毛片| 狂野欧美激情性xxxx| 国产精品影院久久| 欧美日韩福利视频一区二区| 亚洲电影在线观看av| 国产精品久久视频播放| 少妇熟女aⅴ在线视频| 亚洲九九香蕉| 在线国产一区二区在线| 亚洲精品美女久久av网站| 九色国产91popny在线| 久久亚洲精品不卡| www国产在线视频色| 神马国产精品三级电影在线观看 | 国产欧美日韩精品亚洲av| 亚洲av片天天在线观看| 99国产精品一区二区蜜桃av| 国产99白浆流出| 很黄的视频免费| 啦啦啦免费观看视频1| 91大片在线观看| 国产亚洲欧美在线一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产三级在线视频| 国产成人欧美在线观看| 此物有八面人人有两片| 亚洲一卡2卡3卡4卡5卡精品中文| 黄网站色视频无遮挡免费观看| 久久天躁狠狠躁夜夜2o2o| 亚洲av五月六月丁香网| 99精品久久久久人妻精品| 不卡一级毛片| www.www免费av| 亚洲欧美精品综合一区二区三区| 99国产精品一区二区蜜桃av| 日韩国内少妇激情av| 一个人免费在线观看的高清视频| 淫妇啪啪啪对白视频| 午夜福利,免费看| xxx96com| 国产成年人精品一区二区| 久久久久精品国产欧美久久久| 国产又色又爽无遮挡免费看| 国产aⅴ精品一区二区三区波| 岛国视频午夜一区免费看| 热re99久久国产66热| 国产成人精品无人区| 妹子高潮喷水视频| 最新在线观看一区二区三区| 国产91精品成人一区二区三区| 看片在线看免费视频| 最近最新中文字幕大全免费视频| 国产一区二区三区视频了| 最新在线观看一区二区三区| 亚洲成人国产一区在线观看| 国产亚洲精品av在线| 中亚洲国语对白在线视频| 国产欧美日韩一区二区三区在线| 色av中文字幕| 黄色视频不卡| 曰老女人黄片| 久99久视频精品免费| 国产xxxxx性猛交| 亚洲欧洲精品一区二区精品久久久| 国产精品日韩av在线免费观看 | 一个人观看的视频www高清免费观看 | 男女做爰动态图高潮gif福利片 | 中文字幕人妻丝袜一区二区| 国产成+人综合+亚洲专区| 国产成人影院久久av| 亚洲全国av大片| 熟女少妇亚洲综合色aaa.| 国产一区二区三区综合在线观看| 日韩国内少妇激情av| 给我免费播放毛片高清在线观看| 久久久国产欧美日韩av| 9色porny在线观看| 欧美日韩乱码在线| 国产一区二区三区视频了| 久久久国产成人精品二区| 免费在线观看完整版高清| 女人爽到高潮嗷嗷叫在线视频| 日韩高清综合在线| 神马国产精品三级电影在线观看 | 亚洲欧美精品综合久久99| 国产精品一区二区三区四区久久 | 动漫黄色视频在线观看| 成人国语在线视频| 亚洲成av片中文字幕在线观看| 亚洲一区二区三区不卡视频| 成人手机av| 波多野结衣高清无吗| 国产片内射在线| 国产一区二区三区在线臀色熟女| 在线观看免费视频日本深夜| 亚洲成a人片在线一区二区| 国产精品,欧美在线| 午夜日韩欧美国产| 国产精品久久视频播放| 99国产极品粉嫩在线观看| 国产人伦9x9x在线观看| av在线天堂中文字幕| 少妇 在线观看| 久久这里只有精品19| 一进一出好大好爽视频| 一本综合久久免费| 亚洲色图 男人天堂 中文字幕| 国产又爽黄色视频| 啦啦啦韩国在线观看视频| 免费人成视频x8x8入口观看| 国产精品综合久久久久久久免费 | 在线视频色国产色| 啦啦啦观看免费观看视频高清 | 亚洲国产欧美日韩在线播放| 欧美另类亚洲清纯唯美| 91成年电影在线观看| 国产精品亚洲一级av第二区| 欧美黄色淫秽网站| 久久精品影院6| 一边摸一边抽搐一进一出视频| bbb黄色大片| 精品国产美女av久久久久小说| 咕卡用的链子| 亚洲成av人片免费观看| 一a级毛片在线观看| 日日爽夜夜爽网站| 日本 欧美在线| 日韩欧美国产一区二区入口| 久久人妻av系列| 亚洲欧美一区二区三区黑人| 日本免费一区二区三区高清不卡 | 久久久国产精品麻豆| 国产视频一区二区在线看| 国产成人欧美在线观看| 黄色视频不卡| 侵犯人妻中文字幕一二三四区| 桃色一区二区三区在线观看| 十八禁网站免费在线| 欧美日韩精品网址| 操美女的视频在线观看| 99国产精品免费福利视频| 手机成人av网站| 欧美色欧美亚洲另类二区 | 国产精品亚洲美女久久久| 亚洲色图 男人天堂 中文字幕| 精品一品国产午夜福利视频| 国产亚洲精品av在线| 后天国语完整版免费观看| 国产主播在线观看一区二区| 91九色精品人成在线观看| 好男人在线观看高清免费视频 | 深夜精品福利| 波多野结衣巨乳人妻| 人成视频在线观看免费观看| 动漫黄色视频在线观看| 亚洲欧美激情综合另类| 亚洲片人在线观看| 黄色毛片三级朝国网站| 99riav亚洲国产免费| 国产又爽黄色视频| 久久久久久久久中文| 久久性视频一级片| 日本 av在线| 午夜福利影视在线免费观看| av福利片在线| 亚洲激情在线av| 欧美日韩瑟瑟在线播放| 国产国语露脸激情在线看| 1024香蕉在线观看| 日韩国内少妇激情av| 麻豆国产av国片精品| 久久亚洲精品不卡| 久久精品国产综合久久久|