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

    不同品種綿羊肌內(nèi)脂肪沉積相關(guān)miRNA的篩選與功能預(yù)測(cè)

    2023-12-13 15:22:35孟科趙薇郭晨浩聶偉陶毛孩袁曉春孫昊然馮登偵
    關(guān)鍵詞:轉(zhuǎn)錄組綿羊

    孟科 趙薇 郭晨浩 聶偉 陶毛孩 袁曉春 孫昊然 馮登偵

    摘要: 為明確不同品種綿羊(Ovis aries)肌肉組織miRNA的表達(dá)情況,篩選影響綿羊肌內(nèi)脂肪沉積的關(guān)鍵miRNA,以灘羊、杜泊羊和小尾寒羊?yàn)檠芯繉?duì)象,測(cè)定背最長(zhǎng)肌肌內(nèi)脂肪含量并對(duì)其進(jìn)行轉(zhuǎn)錄組測(cè)序(RNA-Seq),篩選綿羊脂肪沉積相關(guān)的差異表達(dá)miRNA。結(jié)果表明,灘羊肌內(nèi)脂肪含量顯著高于杜泊羊,極顯著高于小尾寒羊;杜泊羊肌內(nèi)脂肪含量顯著高于小尾寒羊。轉(zhuǎn)錄組測(cè)序質(zhì)量良好。3個(gè)品種綿羊共鑒定到134個(gè)已知miRNA和153個(gè)預(yù)測(cè)miRNA,這些 miRNA的長(zhǎng)度主要介于21~23 nt,且首位堿基對(duì)U堿基具有明顯偏好性。以P ≤0.05和|log2FC|>1(FC為差異倍數(shù))為篩選條件,3個(gè)比較組共獲得42個(gè)差異表達(dá)miRNA;杜泊羊與灘羊比較組中篩選出22個(gè)差異表達(dá)miRNA;小尾寒羊與杜泊羊比較組中篩選出21個(gè)差異表達(dá) miRNA;小尾寒羊與灘羊比較組中篩選出12個(gè)差異表達(dá)miRNA。篩選出的42個(gè)miRNA靶向到326個(gè)靶基因,顯著富集到減數(shù)分裂I、肌肽-寡糖1,2-α-甘露糖苷酶活性和甘露寡糖甘露糖苷酶活性等303個(gè)基因本體(GO)條目,以及富集到α-亞麻酸代謝、亞油酸代謝、多巴胺能突觸、胰島素信號(hào)傳導(dǎo)途徑及N-聚糖生物合成等320個(gè)京都基因與基因組百科全書(shū)(KEGG)代謝通路。根據(jù)GO和KEGG的功能描述進(jìn)一步篩選出23個(gè)miRNA及其靶向的119個(gè)靶基因可能參與脂肪沉積。進(jìn)一步構(gòu)建miRNA靶向調(diào)控網(wǎng)絡(luò),最終篩選到oar-miR-133、oar-miR-485-5p、novel_6、novel_85和novel_103等綿羊肌內(nèi)脂肪沉積的關(guān)鍵miRNA。隨機(jī)選擇的8個(gè)差異表達(dá)miRNA的實(shí)時(shí)熒光定量PCR結(jié)果與轉(zhuǎn)錄組測(cè)序結(jié)果一致,表明轉(zhuǎn)錄組測(cè)序結(jié)果可靠。本研究為深入開(kāi)展綿羊肌內(nèi)脂肪沉積的分子調(diào)控機(jī)理、高品質(zhì)綿羊品種的選育提供理論依據(jù)。

    關(guān)鍵詞: 綿羊;肌內(nèi)脂肪沉積;轉(zhuǎn)錄組;miRNA

    中圖分類號(hào): ?S826?? 文獻(xiàn)標(biāo)識(shí)碼: A?? 文章編號(hào): 1000-4440(2023)07-1554-13

    Screening and functional prediction of miRNA related to intramuscular fat deposition in different sheep breeds

    MENG Ke1,2, ZHAO Wei1, GUO Chen-hao1, NIE Wei1, TAO Mao-hai1, YUAN Xiao-chun1, SUN Hao-ran1,F(xiàn)ENG Deng-zhen1

    (1.School of Agriculture, Ningxia University, Yinchuan 750021, China;2.Animal Husbandry and Fishery Technology Extension and Service Center, Shapotou District, Zhongwei City, Zhongwei 755000, China)

    Abstract: In order to clarify the expression of miRNAs in muscle tissues of different breeds of sheep (Ovis aries) and screen the key miRNAs affecting intramuscular fat deposition in sheep, Tan sheep, Dorper sheep and Small Tail Han sheep were used as experimental materials. The intramuscular fat content of longissimus dorsi muscle was measured and transcriptome sequencing (RNA-Seq) was performed to screen the differentially expressed miRNAs related to fat deposition in sheep. The results showed that the intramuscular fat content of Tan sheep was significantly higher than that of Dorper sheep, and extremely significantly higher than that of Small Tail Han sheep, and the intramuscular fat content of Dorper sheep was significantly higher than that of Small Tail Han sheep. The quality of transcriptome sequencing was good. A total of 134 known miRNAs and 153 predicted miRNAs were identified in three breeds of sheep. The length of these miRNAs mainly ranged from 21 nt to 23 nt, and the first base had obvious preference for U base. With P≤0.05 and |log2FC|>1 (FC was the fold change)as the screening conditions, a total of 42 differentially expressed miRNAs were obtained from the three comparison groups. Twenty-two differentially expressed miRNAs were screened in the comparison group of Dorper sheep and Tan sheep, 21 differentially expressed miRNAs were screened in the comparison group of Small Tail Han sheep and Dorper sheep, and 12 differentially expressed miRNAs were screened in the comparison group of Small Tail Han sheep and Tan sheep. The screened 42 miRNAs targeted to 326 target genes, which were significantly enriched into 303 gene ontology(GO) entries such as meiosis I, mannosyl-oligosaccharide 1,2-α-mannosidase activity and mannosyl-oligosaccharide mannosidase activity, and 320 Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways such as α-linolenic acid metabolism, linoleic acid metabolism, dopaminergic synapses, insulin signaling pathway and N-glycan biosynthesis. According to the functional description of GO and KEGG, 23 miRNAs and 119 target genes might be involved in fat deposition. The miRNA targeted regulatory network was further constructed, and finally the key miRNAs for fat deposition in sheep such as oar-miR-133, oar-miR-485-5p, novel_6, novel_85 and novel_103 were screened. The qRT-PCR results of the eight randomly selected differentially expressed miRNAs were consistent with transcriptome sequencing results, indicating that the transcriptome sequencing results were reliable. This study provides a theoretical basis for further research on the molecular regulation mechanism of intramuscular fat deposition in sheep and the breeding of high-quality sheep breeds.

    Key words: sheep;intramuscular fat deposition;transcriptome;miRNA

    隨著生活質(zhì)量的提升,人們對(duì)羊肉的需求不僅僅滿足于數(shù)量,羊肉的品質(zhì)更加成為人們關(guān)注的焦點(diǎn)。羊肉的品質(zhì)主要由肉質(zhì)性狀決定,而肌內(nèi)脂肪(IMF)含量是評(píng)價(jià)肉品質(zhì)的重要指標(biāo)。IMF含量的高低影響到羊肉的風(fēng)味、嫩度、適口性及多汁性[1-2]。肌內(nèi)脂肪的形成以及含量受到遺傳、環(huán)境、營(yíng)養(yǎng)等多種因素的共同影響[3]。已有研究結(jié)果表明肌內(nèi)脂肪含量具有較高的遺傳力;在相同的環(huán)境和營(yíng)養(yǎng)條件下,遺傳背景成為影響羊肉肌內(nèi)脂肪的重要因素[4]。探究不同品種綿羊IMF含量差異及其分子調(diào)控機(jī)制對(duì)綿羊肉質(zhì)性狀改善具有重要意義。

    MicroRNA(miRNA) 作為一類不編碼的小分子RNA,主要通過(guò)與mRNA的3′UTR區(qū)結(jié)合來(lái)抑制靶基因的表達(dá)而發(fā)揮作用[5]。近年來(lái),國(guó)內(nèi)外許多學(xué)者研究報(bào)道m(xù)iRNA在動(dòng)物脂肪的形成過(guò)程中發(fā)揮重要的調(diào)控作用。邵靜等[6]研究發(fā)現(xiàn),miR-17-3p通過(guò)靶向調(diào)控KCTD15基因的表達(dá)從而影響延邊黃牛前體脂肪細(xì)胞分化;Ma等[7]研究結(jié)果表明, miR-26b-5p通過(guò)靶向調(diào)控FGF21基因促進(jìn)脂肪生成,參與調(diào)控山羊肌內(nèi)脂肪細(xì)胞的分化。Chen等[8]研究發(fā)現(xiàn)miR-376a靶向調(diào)控KLF15基因,參與秦川牛的脂肪細(xì)胞增殖分化過(guò)程;Deng等[9]發(fā)現(xiàn)miR-27a 通過(guò)特異性結(jié)合視黃醇X受體α(RXRα) 抑制綿羊前脂肪細(xì)胞的分化。隨著測(cè)序技術(shù)不斷發(fā)展,應(yīng)用轉(zhuǎn)錄組測(cè)序(RNA-Seq)技術(shù)探究影響IMF沉積的研究也越來(lái)越多。Cheng等[10]通過(guò)RNA-Seq技術(shù)比較了不同肌內(nèi)脂肪含量豬背最長(zhǎng)肌的miRNA表達(dá)情況,篩選到18個(gè)差異miRNA可能參與調(diào)控脂肪酸合成和脂質(zhì)代謝;喻世剛等[11]通過(guò)比較高和低肌內(nèi)脂肪含量雞胸中miRNA的表達(dá)差異,篩選到5個(gè)miRNA參與到脂肪生成過(guò)程;Huang等[12]通過(guò)比較中國(guó)沼澤水牛肌肉和脂肪組織之間的差異miRNA,篩選到8個(gè)miRNA在水牛肌內(nèi)脂肪沉積過(guò)程中發(fā)揮重要作用。上述研究結(jié)果表明,應(yīng)用RNA-Seq技術(shù)探究IMF沉積相關(guān)miRNA表達(dá)譜的研究較為成熟且已取得一定成果。但目前有關(guān)不同品種綿羊肌內(nèi)脂肪沉積相關(guān)miRNA的研究較少。本研究基于杜泊羊、灘羊和小尾寒羊3個(gè)品種綿羊的肌內(nèi)脂肪含量差異,通過(guò)RNA-Seq技術(shù)鑒定并篩選綿羊 IMF 沉積過(guò)程中差異表達(dá)的miRNA,探索miRNA調(diào)節(jié)綿羊肌內(nèi)脂肪沉積的分子機(jī)制。旨在為探究分子水平調(diào)控綿羊肉質(zhì)性狀提供參考,為高品質(zhì)綿羊品種的選育提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    從寧夏宇泊科技有限公司石嘴山威震肉羊繁育場(chǎng)選擇飼養(yǎng)條件一致的8月齡健康綿羊品種杜泊羊、灘羊和小尾寒羊各8只(公母各4只)為試驗(yàn)材料。對(duì)綿羊進(jìn)行屠宰后,取第3肋骨到第12肋骨之間的背最長(zhǎng)肌組織作為樣品,放置于液氮中冷藏備用。

    1.2 綿羊肌內(nèi)脂肪含量的測(cè)定

    采用索氏脂肪抽提法對(duì)試驗(yàn)綿羊的肌內(nèi)脂肪含量進(jìn)行測(cè)定。具體測(cè)定方法是:(1)取適量冷藏樣品,剔除肉樣筋膜后,絞為肉糜,65 ℃烘干,測(cè)定100 g樣品中水分含量記為A;(2)取測(cè)試樣稱質(zhì)量記為W(g),用濾紙包好;(3)在102 ℃烘箱中烘干至恒質(zhì)量,記為W1(g);(4) 將濾紙包好的樣品放入索氏抽提器的抽提筒內(nèi),倒入2/3體積的乙醚后接上冷凝裝置,在55 ℃恒溫水浴中抽提6~8 h。抽提結(jié)束后,在102 ℃烘箱中烘干至恒質(zhì)量,記為W2(g)。每品種3個(gè)重復(fù)。肌內(nèi)脂肪含量(IMFC)計(jì)算公式為:

    IMFC=(W1-W2)/[(100-A)×W]。

    1.3 總RNA提取與質(zhì)量控制

    選擇上述3個(gè)綿羊品種母羊各3只,采用Trizol法提取樣品背最長(zhǎng)肌組織的總RNA,利用瓊脂糖凝膠電泳、Nanodrop分光光度計(jì)(Therom公司產(chǎn)品,美國(guó))和Agilent 2100 生物分析儀(Agilent公司產(chǎn)品,美國(guó))檢測(cè)RNA樣品的濃度、純度以及完整性。

    1.4 cDNA文庫(kù)的構(gòu)建與測(cè)序

    RNA樣品檢測(cè)合格后,使用小RNA樣品制備試劑盒(Small RNA Sample Pre Kit)進(jìn)行cDNA文庫(kù)的構(gòu)建。構(gòu)建方法為:首先利用小RNA的3′和 5′端的特殊結(jié)構(gòu),以總RNA為起始樣品在小RNA兩端加上接頭,進(jìn)行反轉(zhuǎn)錄合成cDNA,隨后經(jīng)過(guò)PCR擴(kuò)增和PAGE膠電泳對(duì)目標(biāo)DNA片段進(jìn)行分離,切膠回收后得到cDNA文庫(kù)。對(duì)構(gòu)建好的文庫(kù)進(jìn)行檢驗(yàn),檢驗(yàn)合格后在Illumina平臺(tái)進(jìn)行測(cè)序。

    1.5 miRNA鑒定與表達(dá)分析

    獲得原始序列后經(jīng)過(guò)數(shù)據(jù)質(zhì)量控制得到有效序列,對(duì)有效序列進(jìn)行小RNA的種類和長(zhǎng)度篩選。使用Bowtie_0.12.9軟件將篩選得到的長(zhǎng)度為18~35 nt的小RNA對(duì)比到綿羊基因組(Oar_rambouillet_v1.0)參考序列上。將對(duì)比到參考序列的數(shù)據(jù)與miRBase數(shù)據(jù)庫(kù)進(jìn)行比對(duì),注釋分析已知的miRNA;在miRBase數(shù)據(jù)庫(kù)未比對(duì)到的數(shù)據(jù),通過(guò)miREvo_v1.1軟件預(yù)測(cè)新的miRNA,并利用mirdeep2_0_0_5軟件進(jìn)行分析。對(duì)鑒定到的miRNA利用Stringtie v1.3.3軟件進(jìn)行表達(dá)量計(jì)算,采用TPM(transcript per million)值作歸一化處理。利用DESeq_1.2.10軟件分析樣本間miRNA的表達(dá)差異,以P≤0.05且|log2 FC|≥1(FC為差異位數(shù))作為標(biāo)準(zhǔn)篩選表達(dá)差異miRNA。

    1.6 差異miRNA靶基因預(yù)測(cè)及功能注釋和富集分析

    使用miRanda-3.3a和RNAhybrid v2.0軟件預(yù)測(cè)表達(dá)差異miRNA的靶基因,取交集作為預(yù)測(cè)結(jié)果。miRanda-3.3a軟件參數(shù)設(shè)置為:-sc 140 -en -10 -scale 4-strict -out;RNAhybrid v2.0軟件參數(shù)設(shè)置為:-e 10 -p 0.05。以P≤0.05為顯著閾值,分別對(duì)表達(dá)差異miRNA的靶基因進(jìn)行基因本體(GO)功能注釋及京都基因與基因組百科全書(shū)(KEGG)富集分析。

    1.7 候選miRNA與其靶基因的互作網(wǎng)絡(luò)分析

    利用Cytoscape 3.9.1軟件對(duì)篩選到的靶基因繪制互作網(wǎng)絡(luò)圖。

    1.8 測(cè)序結(jié)果可靠性分析

    為驗(yàn)證測(cè)序數(shù)據(jù)的準(zhǔn)確性,利用實(shí)時(shí)熒光定量聚合酶鏈?zhǔn)椒磻?yīng)(qRT-PCR)技術(shù),測(cè)定不同綿羊品種8個(gè)差異表達(dá)miRNA的相對(duì)表達(dá)量,與RNA-seq測(cè)序表達(dá)量進(jìn)行一致性比較。利用Primer Premier 6.0軟件并采用頸環(huán)法設(shè)計(jì)miRNA引物,以U6作為內(nèi)參基因(表1)。熒光染料為TB Green Premix Ex Taq TMⅡ(TaKaBa,大連),所用的引物由生工生物工程(上海)股份有限公司合成。PCR反應(yīng)體系為:2×TB Green Premix Ex Taq II 10.0 μl,正向引物0.8 μl,反向引物0.8 μl,cDNA模板2.0 μl,RNase-free雙蒸水6.4 μl。qPCR擴(kuò)增程序?yàn)椋?5 ℃預(yù)變性2 min;95 ℃變性5 s,58 ℃退火30 s,39個(gè)循環(huán);95 ℃延伸5 s,60 ℃維持1 min。每個(gè)樣品設(shè)置3次重復(fù)。

    1.9 統(tǒng)計(jì)分析

    以2-△△Ct方法進(jìn)行miRNA相對(duì)表達(dá)量的計(jì)算,利用SPSS 26.0軟件對(duì)3個(gè)品種綿羊的肌內(nèi)脂肪含量和miRNA的表達(dá)量進(jìn)行單因素方差分析。以P<0.05表示差異顯著,P<0.01表示差異極顯著。實(shí)時(shí)熒光定量結(jié)果圖片由GraphPad Prism V8.0軟件繪制。

    2 結(jié)果與分析

    2.1 綿羊肌內(nèi)脂肪含量差異分析

    3個(gè)綿羊品種肌內(nèi)脂肪含量如表2所示。灘羊肌內(nèi)脂肪含量最高,達(dá)5.42%,而小尾寒羊肌內(nèi)脂肪含量最低,2.10%。灘羊肌內(nèi)脂肪含量顯著高于杜泊羊,且極顯著高于小尾寒羊;杜泊羊肌內(nèi)脂肪含量顯著高于小尾寒羊。

    2.2 RNA-seq數(shù)據(jù)分析

    本研究建立3個(gè)品種綿羊的9個(gè)背最長(zhǎng)肌組織樣品文庫(kù),通過(guò)llumina平臺(tái)進(jìn)行深度測(cè)序后,獲得原始序列共108 813 537條,質(zhì)控后獲得有效序列共107 166 390條。每個(gè)樣本的Q20在99%以上,Q30在96%以上,錯(cuò)誤率均為0.01%,堿基組成穩(wěn)定均衡(表3)。上述結(jié)果表明測(cè)序質(zhì)量良好,可用于進(jìn)一步分析。

    2.3 miRNA的注釋與鑒定

    在18~35 nt范圍內(nèi),篩選各樣品中的有效序列[13]。篩選到的有效序列長(zhǎng)度在22 nt的序列所占比例最高,灘羊和小尾寒羊22 nt的read所占比例高于杜泊羊(圖1)。每個(gè)樣本至少有87%以上的小RNA序列比對(duì)到了綿羊基因參考序列上(表4)。在3種綿羊肌肉組織中共鑒定到134個(gè)已知miRNA和153個(gè)預(yù)測(cè)miRNA。對(duì)鑒定到的miRNA進(jìn)行首位堿基偏好性分析發(fā)現(xiàn),已知miRNA和預(yù)測(cè)miRNA的首位堿基均對(duì)U堿基具有明顯偏好性(圖2)。

    2.4 差異表達(dá)miRNA的篩選與分析

    對(duì)3個(gè)綿羊品種miRNA 的表達(dá)量進(jìn)行兩兩比較分析。 杜泊羊與灘羊比較組中篩選出差異表達(dá)miRNA 22個(gè),其中上調(diào)表達(dá)13個(gè),下調(diào)表達(dá)9個(gè);小尾寒羊與杜泊羊比較組中篩選出21個(gè)差異表達(dá)miRNA,其中上調(diào)表達(dá)13個(gè),下調(diào)表達(dá)8個(gè);小尾寒羊與灘羊比較組中篩選出12個(gè)差異表達(dá) miRNA,其中上調(diào)表達(dá)9個(gè),下調(diào)表達(dá)3個(gè)(圖3)。3個(gè)比較組中沒(méi)有共表達(dá)的差異表達(dá)miRNA,小尾寒羊與灘羊比較組與杜泊羊與灘羊比較組,小尾寒羊與杜泊羊比較組之間均有4個(gè)共表達(dá)的差異表達(dá)miRNA,杜泊羊與灘羊比較組與小尾寒羊與杜泊羊比較組之間有5個(gè)共表達(dá)的差異表達(dá)miRNA。3個(gè)比較組中共篩選出42個(gè)差異表達(dá)miRNA(圖4)。篩選出的miRNA包括23個(gè)已知miRNA和19個(gè)新預(yù)測(cè)的miRNA。篩選出的已知miRNA的差異表達(dá)信息如表5所示。

    2.5 差異表達(dá)miRNA靶基因預(yù)測(cè)及功能分析

    篩選出的42個(gè)差異表達(dá) miRNA預(yù)測(cè)到326個(gè)靶基因,形成411個(gè)靶向關(guān)系對(duì)。其中,novel_85、oar-miR-133和oar-miR-485-5p預(yù)測(cè)到的靶基因數(shù)量較多。novel_85預(yù)測(cè)到99個(gè)靶基因,包括糖原合酶基因(GYS1)、組蛋白脫乙酰酶3基因(HDAC3)、胰島素樣生長(zhǎng)因子II基因(IGF2)等;oar-miR-133預(yù)測(cè)到55個(gè)靶基因,包括肌球蛋白 IA基因(MYO1A)、鳥(niǎo)氨酸脫羧酶1基因(ODC1)、磷脂酶 A2基因(PLA2)等;oar-miR-485-5p預(yù)測(cè)到29個(gè)靶基因,包括甘露糖-6-磷酸異構(gòu)酶基因(MPI)、羥酰基谷胱甘肽水解酶基因(HAGH)、跨膜絲氨酸蛋白酶2基因(TMPRSS2)等。對(duì)每組預(yù)測(cè)到的靶基因進(jìn)行功能分析。 GO功能富集分析發(fā)現(xiàn),杜泊羊與灘羊比較組顯著富集到111個(gè)GO條目,其中生物學(xué)過(guò)程51個(gè),細(xì)胞組分24個(gè),分子功能36個(gè);小尾寒羊與杜泊羊比較組顯著富集到97個(gè)GO條目,其中生物學(xué)過(guò)程51個(gè),細(xì)胞組分11個(gè),分子功能35個(gè);小尾寒羊與灘羊比較組顯著富集到95個(gè)GO條目,其中生物學(xué)過(guò)程41個(gè),細(xì)胞組分6個(gè),分子功能48個(gè)。根據(jù)P值大小每組前20的GO條目如圖5所示,顯著富集在減數(shù)分裂、減數(shù)分裂前期、清道夫受體活性、貨物受體活性、肌肽-寡糖1,2-α-甘露糖苷酶活性和甘露寡糖甘露糖苷酶活性等GO條目。KEGG信號(hào)通路富集分析發(fā)現(xiàn),杜泊羊與灘羊比較組富集到106個(gè)通路,7個(gè)通路顯著富集;小尾寒羊與杜泊羊比較組富集到168個(gè)通路,4個(gè)通路顯著富集;小尾寒羊與灘羊比較組富集到46個(gè)通路,5個(gè)通路顯著富集。根據(jù)P值大小每組前20個(gè)信號(hào)富集通路如圖6所示,顯著富集的通路主要有α-亞麻酸代謝、亞油酸代謝、多巴胺能突觸、胰島素信號(hào)傳導(dǎo)途徑、N-聚糖的生物合成和苯丙氨酸、酪氨酸和色氨酸的生物合成等代謝通路。根據(jù)GO和KEGG的功能描述篩選出23個(gè)差異表達(dá)miRNA及其靶向的119個(gè)基因可能參與脂肪沉積。

    2.6 miRNA靶向調(diào)控網(wǎng)絡(luò)

    篩選出的119個(gè)可能參與脂肪沉積的靶基因與23個(gè)差異表達(dá)miRNA構(gòu)成129個(gè)靶向關(guān)系對(duì),其調(diào)控網(wǎng)絡(luò)如圖7所示。通過(guò)調(diào)控網(wǎng)絡(luò)圖發(fā)現(xiàn),oar-miR-133、oar-miR-485-5p、novel_6、novel_85和novel_103靶向到較多的靶基因且構(gòu)成了完整的調(diào)控網(wǎng)絡(luò),故推測(cè)其可能是調(diào)控綿羊脂肪沉積的關(guān)鍵miRNA。

    2.7 miRNA表達(dá)驗(yàn)證

    利用qRT-PCR得到3個(gè)綿羊品種8個(gè)差異表達(dá)miRNA(oar-miRNA-99a、oar-miR-26a、oar-miR-374b、oar-miR-19b、novel_9、novel_202、 novel_82、novel_62)的相對(duì)表達(dá)量與RNA-seq測(cè)序得到的這8個(gè)miRNA表達(dá)趨勢(shì)基本保持一致(圖8),說(shuō)明測(cè)序結(jié)果真實(shí)可靠。

    3 討論

    3.1 3個(gè)品種綿羊肌內(nèi)脂肪含量差異

    肉質(zhì)性狀是綿羊關(guān)鍵的經(jīng)濟(jì)性狀,也是育種者關(guān)注的重要問(wèn)題。IMF含量的高低影響羊肉的食用價(jià)值。灘羊作為中國(guó)特有的裘皮綿羊品種,其肉質(zhì)細(xì)嫩、無(wú)膻味;杜泊羊原產(chǎn)自南非國(guó)家,產(chǎn)肉性能高,是育肥生產(chǎn)中優(yōu)質(zhì)品種;而小尾寒羊具有繁殖率高的特點(diǎn),是一產(chǎn)多羔的代表性品種[14]。3個(gè)品種中,灘羊的肌內(nèi)脂肪含量最高,而小尾寒羊最低。灘羊和杜泊羊作為優(yōu)質(zhì)的肉羊品種,肌內(nèi)脂肪含量均高于小尾寒羊。

    3.2 3個(gè)品種綿羊背最長(zhǎng)肌miRNA的鑒定及特征

    微小核糖核酸(miRNA)是由20~24個(gè)核苷酸堿基組成的非編碼內(nèi)源性小分子RNA[15]。本研究從3個(gè)品種綿羊背最長(zhǎng)肌中共鑒定到134個(gè)已知miRNA和153個(gè)預(yù)測(cè)miRNA,其長(zhǎng)度主要分布在21~23 nt,且長(zhǎng)度為22 nt的miRNA占比最高,符合miRNA的長(zhǎng)度分布特征。對(duì)鑒定到的miRNA進(jìn)行首位堿基偏向性分析,發(fā)現(xiàn)首位堿基偏向于U堿基。這與Liu等[16]在綿羊骨骼肌中鑒定到miRNA的首堿基偏向性和付雪峰等[17]在藏絨山羊皮膚組織中鑒定出的miRNA的首位堿基偏向性一致。miRNA的首位堿基對(duì)U堿基的偏好性可能與miRNA的形成以及其對(duì)靶基因的沉默機(jī)制有關(guān)[18]。

    3.3 綿羊肌內(nèi)脂肪含量相關(guān)miRNA的篩選

    本研究在杜泊羊、灘羊和小尾寒羊的背最長(zhǎng)肌中篩選獲得42個(gè)差異表達(dá)miRNA,包括23個(gè)已知的miRNA和19個(gè)預(yù)測(cè)的miRNA。進(jìn)一步通過(guò)表達(dá)量差異篩選出23個(gè)可能與綿羊背最長(zhǎng)肌脂肪沉積相關(guān)的miRNA。篩選出的這些已知miRNA已有部分被證實(shí)參與調(diào)節(jié)脂肪沉積和脂質(zhì)代謝。例如,miR-433-3p對(duì)支鏈酮酸脫氫酶β(BCKDHB)有負(fù)調(diào)控作用,影響綿羊前體脂肪細(xì)胞分化從而參與調(diào)節(jié)綿羊的脂肪代謝[19]; miR-27a 通過(guò)調(diào)節(jié)靶基因RXRα的表達(dá)來(lái)增強(qiáng)綿羊前脂肪細(xì)胞增殖并抑制綿羊前脂肪細(xì)胞分化,在調(diào)節(jié)綿羊肌肉脂肪和皮下脂肪組織之間的差異脂肪積累中起重要作用[9];miR-143是與脂質(zhì)代謝密切相關(guān)的miRNA之一,通過(guò)調(diào)節(jié)靶基因的表達(dá)來(lái)調(diào)節(jié)奶牛脂肪細(xì)胞分化和甘油三酯合成[20]。Wang等[21]研究發(fā)現(xiàn)miR-26a 和miR-143在成年肉牛肌內(nèi)脂肪中的表達(dá)量顯著高于皮下脂肪; 汪海洋等[22]研究發(fā)現(xiàn)miR-143在西門塔爾牛肌內(nèi)脂肪組織中表達(dá)量高于皮下脂肪。本研究中,miR-26a作為杜泊羊與灘羊比較組和小尾寒羊與灘羊比較組中篩選出共有的差異表達(dá)miRNA,在灘羊肌肉組織中表達(dá)量顯著高于杜泊羊和小尾寒羊,且杜泊羊的表達(dá)量略高于小尾寒羊;同樣,miR-143在杜泊羊與灘羊比較組中被篩選出來(lái),該miRNA在灘羊肌肉組織中表達(dá)量高于杜泊羊。在3個(gè)品種綿羊肌內(nèi)脂肪含量分析中,灘羊IMF含量最高而小尾寒羊最低,上述2個(gè)miRNA表達(dá)量的研究結(jié)果也符合該特征。同樣,let-7b是小尾寒羊與杜泊羊比較組和小尾寒羊與灘羊比較組中篩選出共有的差異表達(dá)miRNA,且與miR-26a在3個(gè)品種綿羊肌肉組織中的表達(dá)趨勢(shì)一致,即高肌內(nèi)脂肪含量的組織中l(wèi)et-7b表達(dá)量亦高。這個(gè)結(jié)果與喻世剛等[11]發(fā)現(xiàn)的雞胸肌內(nèi)脂肪含量與let-7b表達(dá)量的關(guān)系一致。上述結(jié)果也進(jìn)一步說(shuō)明本研究篩選出的miRNA作為與綿羊肌內(nèi)脂肪含量相關(guān)的候選miRNA的可靠性。但這些miRNA對(duì)綿羊肌肉脂肪含量的影響機(jī)制還需進(jìn)一步研究。

    3.4 差異miRNA的靶基因預(yù)測(cè)及功能分析

    篩選出的42個(gè)差異表達(dá)miRNA預(yù)測(cè)到326個(gè)靶基因,形成411個(gè)靶向關(guān)系對(duì)。值得注意的是novel_85預(yù)測(cè)到99個(gè)靶基因,為差異表達(dá)miRNA中靶向到最多候選基因的miRNA。這些候選的靶向基因可以參與脂肪沉積。GYS1基因與糖原代謝相關(guān),其突變影響到肌糖原的沉積[23];HDAC3蛋白的復(fù)合物使過(guò)氧化物酶體增殖物激活受體γ2(Peroxisome proliferator-activated receptor γ2)基因(PPARG2)轉(zhuǎn)錄失活,從而參與到脂肪生成過(guò)程[24];IGF2基因的突變?cè)斐韶i肌內(nèi)脂肪和脂肪酸含量存在差異,可以參與到脂肪組織中脂肪酸的形成。上述研究結(jié)果表明,novel_85可能是參與綿羊脂肪沉積及脂肪代謝的關(guān)鍵miRNA。另外也有一些差異表達(dá)miRNA靶向的基因與脂肪沉積和脂肪代謝相關(guān)的研究。如oar-miR-133和oar-miR-485-5p同時(shí)靶向的谷胱甘肽硫轉(zhuǎn)移酶基因(P1GSTP1)可以參與到脂質(zhì)代謝過(guò)程[25];novel_103靶向的丙酮酸脫氫酶激酶4基因(PDK4),在中國(guó)草原紅牛的研究中發(fā)現(xiàn)該基因的多態(tài)性與其肌內(nèi)脂肪含量相關(guān)[26];oar-miR-433-3p靶向的色氨酸3基因(SESN3)參與調(diào)節(jié)前脂肪細(xì)胞脂肪生成[27]。為了進(jìn)一步了解篩選出的差異miRNA的功能,對(duì)其預(yù)測(cè)的靶基因進(jìn)行功能分析。GO功能富集分析發(fā)現(xiàn),富集顯著的條目包括糖原代謝過(guò)程、油酸酯水合酶活性、三羧酸循環(huán)過(guò)程、低密度脂蛋白顆粒受體結(jié)合等與糖代謝及脂質(zhì)代謝相關(guān)的條目。KEGG信號(hào)通路富集分析發(fā)現(xiàn),可以富集到α-亞麻酸代謝、亞油酸代謝、多巴胺能突觸、胰島素信號(hào)傳導(dǎo)途徑以及脂肪的消化和吸收通路。已有研究結(jié)果表明,α-亞麻酸會(huì)影響奶牛乳腺中的脂肪酸代謝[28],亞油酸可以增加呼吸鏈解偶聯(lián)蛋白1(UCP1)和肉堿-棕櫚油轉(zhuǎn)移酶1 (CPT1)的表達(dá),促進(jìn)脂肪降解[29],并且不同共軛亞油酸對(duì)小鼠的脂肪沉積具有顯著差異 [30]。上述研究結(jié)果表明,篩選出差異表達(dá)miRNA的靶基因可能參與脂肪沉積。根據(jù)GO和KEGG的功能描述篩選出可能參與綿羊脂肪沉積的119個(gè)靶基因,并對(duì)靶向關(guān)系進(jìn)行調(diào)控網(wǎng)絡(luò)構(gòu)建。研究認(rèn)為,oar-miR-133、oar-miR-485-5p、novel_6、novel_85和novel_103可能是參與綿羊肌內(nèi)脂肪沉積的關(guān)鍵miRNA。但miRNA 對(duì)靶基因和 IMF 沉積的影響和調(diào)控機(jī)制尚需進(jìn)行進(jìn)一步闡明。

    4 結(jié)論

    杜泊羊,灘羊和小尾寒羊這3個(gè)品種綿羊IMF含量存在差異,灘羊肌內(nèi)脂肪含量最高而小尾寒羊最低。篩選出23個(gè)已知的和19個(gè)未知的差異表達(dá)miRNA可能參與綿羊肌內(nèi)脂肪沉積,oar-miR-133、oar-miR-485-5p、novel_6、novel_85和novel_103可能是綿羊肌內(nèi)脂肪沉積的關(guān)鍵調(diào)控miRNA。

    參考文獻(xiàn):

    [1] FRANK D, WATKINS P, BALL A, et al. Impact of brassica and lucerne finishing feeds and intramuscular fat on lamb eating quality and flavor. a cross-cultural study using Chinese and non-Chinese Australian consumers [J]. J Agric Food Chem,2016, 64(36): 6856-6868.

    [2] 閆 偉,王玉濤,張永浩,等. 綿羊肌內(nèi)前體脂肪細(xì)胞CNR1和FABP4基因表達(dá)研究[J]. 中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào), 2022, 24(3): 95-102.

    [3] LIANG C, QIAO L, HAN Y, et al. Regulatory roles of SREBF1 and SREBF2 in lipid metabolism and deposition in two Chinese representative fat-tailed sheep breeds[J]. Animals, 2020, 10(8): 1317.

    [4] MOHSEN A, HOSEIN M S, MOHAMMAD M S, et al. Association of the calpastatin genotypes, haplotypes, and SNPs with meat quality and fatty acid composition in two Iranian fat- and thin-tailed sheep breeds[J]. Small Ruminant Research,2016, 149: 40-51.

    [5] HORAK M, NOVAK J, BIENERTOVA-VASKU J. Muscle-specific microRNAs in skeletal muscle development[J]. Dev Biol, 2016, 410(1): 1-13.

    [6] 邵 靜,張珈溯,尹寶珍,等. miR-17-3p靶向KCTD15調(diào)控延邊黃牛前體脂肪細(xì)胞分化[J].畜牧獸醫(yī)學(xué)報(bào), 2020, 51(11): 2689-2698.

    [7] MA J Q, LIN Y Q, ZHU J J, et al. MiR-26b-5p regulates the preadipocyte differentiation by targeting FGF21 in goats[J]. In Vitro Cellular & Developmental Biology-Animal, 2021,57: 257-263.

    [8] CHEN X Y, SAYED H A R, CHENG G, et al. Bta-miR-376a targeting KLF15 interferes with adipogenesis signaling pathway to promote differentiation of Qinchuan beef cattle preadipocytes[J]. Animals, 2020, 10(12): 2362.

    [9] DENG K, REN C, FAN Y, et al. miR-27a is an important adipogenesis regulator associated with differential lipid accumulation between intramuscular and subcutaneous adipose tissues of sheep[J]. Domest Anim Endocrinol,2020,71: 106393.

    [10]CHENG F, LIANG J, YANG L, et al. Systematic identification and comparison of the expressed profiles of lncRNAs, miRNAs, circRNAs, and mRNAs with associated co-expression networks in pigs with low and high intramuscular fat[J]. Animals, 2021,11(11):3212.

    [11]喻世剛,王 鋼,廖 娟,等. 沐川烏骨黑雞肌內(nèi)脂肪沉積相關(guān)microRNAs的篩查與鑒定[J].中國(guó)畜牧獸醫(yī), 2021, 48(8): 2713-2726.

    [12]HUANG J P, WANG S Z, FENG X, et al. miRNA transcriptome comparison between muscle and adipose tissues indicates potential miRNAs associated with intramuscular fat in Chinese swamp buffalo[J]. Genome, 2019, 62(11): 729-738.

    [13]BARTEL D P. Metazoan MicroRNAs[J]. Cell,2018, 173(1): 20-51.

    [14]孟 科,張?zhí)炻?,?鵬,等. 綿羊MYF5基因多態(tài)性及其與生長(zhǎng)性狀的關(guān)聯(lián)分析[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2022, 30(3): 496-505.

    [15]BUSHATI N, COHEN S M. microRNA functions[J]. Annu Rev Cell Dev Biol, 2007, 23: 175-205.

    [16]LIU Z J, LI C Y, LI X Y, et al. Expression profiles of microRNAs in skeletal muscle of sheep by deep sequencing[J]. Asian-Australasian Journal of Animal Sciences, 2019, 32(6): 757-766.

    [17]付雪峰,趙冰茹,索朗達(dá),等. 不同絨細(xì)度的西藏絨山羊皮膚組織miRNA分析與鑒定[J].農(nóng)業(yè)生物技術(shù)學(xué)報(bào),2021, 29(11): 2118-2128.

    [18]袁鈺潔,周婧雯,殷 實(shí),等.不同發(fā)育階段牦牛睪丸組織miRNA的分析及鑒定[J].中國(guó)獸醫(yī)學(xué)報(bào),2022,42(1):165-174.

    [19]閆曉茹,師 濤,潘洋洋,等. MiR-433-3p靶向調(diào)節(jié)綿羊BCKDHB表達(dá)[J].中國(guó)農(nóng)業(yè)科學(xué), 2017, 50(22): 4389-4397.

    [20]ZHANG L, WU Z Q, WANG Y J, et al. MiR-143 regulates milk fat synthesis by targeting smad3 in bovine mammary epithelial cells[J]. Animals, 2020, 10(9): 1453.

    [21]WANG H Y, ZHENG Y, WANG G L, et al. Identification of microRNA and bioinformatics target gene analysis in beef cattle intramuscular fat and subcutaneous fat[J]. Molecular Biosystems,2013, 9(8): 2154-2162.

    [22]汪海洋,鄭 月,李惠俠,等. 西門塔爾牛肌內(nèi)和皮下脂肪miRNA表達(dá)譜及miR-27b靶基因分析[J]. 中國(guó)農(nóng)業(yè)科學(xué),2013, 46(18): 3894-3900.

    [23]FANG L, SHEN B, IRWIN D M, et al. Parallel evolution of the glycogen synthase 1 (muscle) gene Gys1 between old world and new world fruit bats (order: chiroptera)[J]. Biochem Genet, 2014, 52: 443-458.

    [24]LAI P H, WANG W L, KO C Y, et al. HDAC1/HDAC3 modulates PPARG2 transcription through the sumoylated CEBPD in hepatic lipogenesis[J]. Biochimica et Biophysica Acta,2008, 1783(10): 1803-1814.

    [25]王 晶.延黃牛GSTP1基因遺傳變異分析及其對(duì)脂代謝的影響[D].延吉:延邊大學(xué),2021.

    [26]馬彥茹,于永生,曹 陽(yáng),等. 中國(guó)草原紅牛PDK4基因多態(tài)性及其與肉質(zhì)性狀的關(guān)聯(lián)分析[J]. 中國(guó)畜牧獸醫(yī), 2022, 49(6): 2186-2194.

    [27]LIN W, ZHAO J, YAN M, et al. SESN3 inhibited SMAD3 to relieve its suppression for MiR-124, thus regulating pre-adipocyte adipogenesis[J]. Genes,2021,12(12):1852.

    [28]YAMG G, BU D P, WAMG J Q, et al. Duodenal infusion of α-linolenic acid affects fatty acid metabolism in the mammary gland of lactating dairy cows[J]. Journal of Dairy Science,2012,95(10):5821-5830.

    [29]伏春燕,張 燕,魏祥法,等. 共軛亞油酸降低脂肪沉積的分子機(jī)制研究進(jìn)展[J]. 動(dòng)物營(yíng)養(yǎng)學(xué)報(bào),2019,31(8): 3456-3462.

    [30]陳 林,趙偉杰,張楓琳,等. 不同共軛亞油酸異構(gòu)體對(duì)小鼠脂肪沉積、能量代謝和腸道微生物的影響[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2022, 43(3): 1-8.

    (責(zé)任編輯:石春林)

    猜你喜歡
    轉(zhuǎn)錄組綿羊
    在雨中“躲雨”的綿羊
    山羊叭叭和綿羊咩咩
    數(shù)綿羊
    數(shù)綿羊
    奔跑的綿羊
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    基于轉(zhuǎn)錄組測(cè)序的山茱萸次生代謝生物合成相關(guān)基因的挖掘
    金釵石斛轉(zhuǎn)錄組SSR位點(diǎn)信息分析
    人參屬藥用植物轉(zhuǎn)錄組研究進(jìn)展
    巧計(jì)得綿羊
    无遮挡黄片免费观看| 久久久久国内视频| 桃色一区二区三区在线观看| 国产成人免费无遮挡视频| 亚洲一区中文字幕在线| 亚洲第一av免费看| 国产av精品麻豆| 在线观看日韩欧美| 亚洲色图综合在线观看| 亚洲五月色婷婷综合| 欧美一区二区精品小视频在线| 热99国产精品久久久久久7| 亚洲一区中文字幕在线| 自线自在国产av| 宅男免费午夜| 一级,二级,三级黄色视频| 首页视频小说图片口味搜索| 国产精品久久电影中文字幕| 母亲3免费完整高清在线观看| 嫩草影院精品99| 男女下面插进去视频免费观看| 18禁黄网站禁片午夜丰满| 天堂影院成人在线观看| 欧美成人午夜精品| 国产91精品成人一区二区三区| 91成年电影在线观看| 波多野结衣一区麻豆| 国产片内射在线| 色婷婷久久久亚洲欧美| 美女午夜性视频免费| 99精国产麻豆久久婷婷| 国产97色在线日韩免费| 两性午夜刺激爽爽歪歪视频在线观看 | 1024香蕉在线观看| 超色免费av| 国产激情久久老熟女| 亚洲一区二区三区不卡视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av第一区精品v没综合| 人人妻人人澡人人看| 精品国产乱子伦一区二区三区| 亚洲欧美日韩另类电影网站| 久久精品国产清高在天天线| 黄片大片在线免费观看| 国产精品一区二区免费欧美| 亚洲精品中文字幕在线视频| 日本精品一区二区三区蜜桃| 日韩大码丰满熟妇| 人人澡人人妻人| 久久香蕉精品热| 日韩欧美一区视频在线观看| 51午夜福利影视在线观看| 99久久精品国产亚洲精品| 脱女人内裤的视频| 欧美+亚洲+日韩+国产| 亚洲精品美女久久久久99蜜臀| 国产黄色免费在线视频| 天堂影院成人在线观看| 天堂俺去俺来也www色官网| 男女做爰动态图高潮gif福利片 | 久久久久久久精品吃奶| 久久伊人香网站| 国产高清国产精品国产三级| 欧美成人免费av一区二区三区| 免费观看精品视频网站| 久久精品国产综合久久久| 精品国产一区二区久久| 热99re8久久精品国产| 黄片播放在线免费| 老熟妇仑乱视频hdxx| 人人妻人人爽人人添夜夜欢视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产亚洲av麻豆专区| 悠悠久久av| 男女床上黄色一级片免费看| 99热国产这里只有精品6| 免费在线观看日本一区| 国产男靠女视频免费网站| 精品久久久久久久毛片微露脸| 一级作爱视频免费观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲 欧美 日韩 在线 免费| 欧美日本中文国产一区发布| 久久精品亚洲av国产电影网| av天堂久久9| 一级黄色大片毛片| 欧美黑人欧美精品刺激| 波多野结衣一区麻豆| 人人澡人人妻人| 久久久久久免费高清国产稀缺| 国产91精品成人一区二区三区| 在线观看免费日韩欧美大片| 免费一级毛片在线播放高清视频 | 美女高潮到喷水免费观看| 99精品久久久久人妻精品| 亚洲 欧美一区二区三区| 国产无遮挡羞羞视频在线观看| 国产高清视频在线播放一区| 久久中文字幕人妻熟女| 日韩有码中文字幕| 久久久久久久久免费视频了| 国产精品国产高清国产av| 午夜影院日韩av| 91av网站免费观看| 丰满人妻熟妇乱又伦精品不卡| 美女国产高潮福利片在线看| 精品国产美女av久久久久小说| cao死你这个sao货| 国内久久婷婷六月综合欲色啪| 又紧又爽又黄一区二区| 日韩精品免费视频一区二区三区| 在线观看免费视频网站a站| 搡老岳熟女国产| 亚洲av第一区精品v没综合| 亚洲avbb在线观看| 精品一区二区三区视频在线观看免费 | 亚洲国产精品一区二区三区在线| 国产成人av激情在线播放| 悠悠久久av| 午夜福利在线免费观看网站| 欧美中文综合在线视频| 不卡一级毛片| 国产有黄有色有爽视频| 亚洲欧洲精品一区二区精品久久久| 亚洲成人免费电影在线观看| av在线天堂中文字幕 | 日韩欧美国产一区二区入口| 精品国产乱码久久久久久男人| 亚洲欧美精品综合久久99| 午夜激情av网站| 免费av毛片视频| 亚洲av电影在线进入| 亚洲欧美日韩另类电影网站| 在线观看www视频免费| 两个人免费观看高清视频| 午夜福利在线免费观看网站| 久久久久久久精品吃奶| 成年人免费黄色播放视频| 一个人免费在线观看的高清视频| 精品免费久久久久久久清纯| 国产99白浆流出| 亚洲av五月六月丁香网| 一本综合久久免费| 夜夜夜夜夜久久久久| 老汉色∧v一级毛片| 免费在线观看日本一区| 99久久国产精品久久久| 日本vs欧美在线观看视频| 亚洲欧美日韩高清在线视频| 久久香蕉精品热| avwww免费| 久久人妻熟女aⅴ| 91国产中文字幕| 18禁黄网站禁片午夜丰满| 国产精品一区二区在线不卡| 91在线观看av| 久久久国产欧美日韩av| 日日干狠狠操夜夜爽| a在线观看视频网站| 婷婷丁香在线五月| 麻豆久久精品国产亚洲av | 一边摸一边抽搐一进一小说| 久久久久久久久免费视频了| 欧美一区二区精品小视频在线| 亚洲全国av大片| 99精国产麻豆久久婷婷| 精品午夜福利视频在线观看一区| 欧美日韩av久久| 巨乳人妻的诱惑在线观看| 精品免费久久久久久久清纯| 久久久国产一区二区| 91精品国产国语对白视频| 国产精品亚洲一级av第二区| 亚洲人成77777在线视频| 亚洲人成电影免费在线| netflix在线观看网站| 久久欧美精品欧美久久欧美| 免费高清在线观看日韩| 精品国产一区二区久久| 麻豆一二三区av精品| 美女大奶头视频| 俄罗斯特黄特色一大片| 欧美一级毛片孕妇| 女人爽到高潮嗷嗷叫在线视频| 久久久久国产一级毛片高清牌| 香蕉久久夜色| 天堂动漫精品| 麻豆国产av国片精品| xxx96com| 18禁观看日本| 免费女性裸体啪啪无遮挡网站| 国产av一区在线观看免费| 欧美日韩av久久| 亚洲国产看品久久| 亚洲五月天丁香| 久久久国产一区二区| 97人妻天天添夜夜摸| 美女午夜性视频免费| 欧美日韩视频精品一区| 久久国产亚洲av麻豆专区| 欧美亚洲日本最大视频资源| 国产欧美日韩精品亚洲av| 真人一进一出gif抽搐免费| 久久久久久久久免费视频了| 97碰自拍视频| 嫩草影院精品99| 别揉我奶头~嗯~啊~动态视频| 国产成人欧美在线观看| a级毛片黄视频| 国产xxxxx性猛交| 在线看a的网站| 性欧美人与动物交配| 亚洲一区二区三区不卡视频| 精品一品国产午夜福利视频| 91国产中文字幕| 精品久久蜜臀av无| 大型黄色视频在线免费观看| videosex国产| 性少妇av在线| 亚洲国产精品合色在线| 日本 av在线| 国产精品美女特级片免费视频播放器 | 国产精品 国内视频| 级片在线观看| 精品国产亚洲在线| 波多野结衣av一区二区av| 久久99一区二区三区| 中出人妻视频一区二区| 一级黄色大片毛片| 美女高潮到喷水免费观看| 久久伊人香网站| 国产亚洲精品久久久久久毛片| 国产91精品成人一区二区三区| 成人手机av| 国产精品一区二区免费欧美| 亚洲av成人不卡在线观看播放网| 久久中文字幕人妻熟女| 美女福利国产在线| 国产精品影院久久| 中文字幕另类日韩欧美亚洲嫩草| 国内久久婷婷六月综合欲色啪| 夜夜看夜夜爽夜夜摸 | 成年女人毛片免费观看观看9| 美国免费a级毛片| 亚洲国产毛片av蜜桃av| 亚洲国产欧美一区二区综合| 人妻丰满熟妇av一区二区三区| 午夜精品国产一区二区电影| 中文欧美无线码| 黄色怎么调成土黄色| 少妇粗大呻吟视频| 亚洲欧美一区二区三区黑人| 久久国产精品人妻蜜桃| 999久久久国产精品视频| 99精国产麻豆久久婷婷| 在线播放国产精品三级| 久久精品91蜜桃| www.999成人在线观看| 国产av又大| 色婷婷久久久亚洲欧美| 少妇被粗大的猛进出69影院| 精品日产1卡2卡| 纯流量卡能插随身wifi吗| 色精品久久人妻99蜜桃| 午夜福利影视在线免费观看| 露出奶头的视频| 91大片在线观看| 亚洲美女黄片视频| 久久这里只有精品19| 久久精品成人免费网站| 99国产精品99久久久久| 女人高潮潮喷娇喘18禁视频| 久久中文字幕人妻熟女| www.www免费av| 男女之事视频高清在线观看| 国产精品98久久久久久宅男小说| 动漫黄色视频在线观看| 啦啦啦 在线观看视频| 精品熟女少妇八av免费久了| 国产精品野战在线观看 | 国产一卡二卡三卡精品| 久久性视频一级片| 国产1区2区3区精品| 自线自在国产av| 男女之事视频高清在线观看| 天堂影院成人在线观看| 久久人人97超碰香蕉20202| 亚洲国产精品一区二区三区在线| av天堂久久9| 国产深夜福利视频在线观看| 亚洲黑人精品在线| 国产在线精品亚洲第一网站| 久久亚洲精品不卡| 久久九九热精品免费| aaaaa片日本免费| 桃色一区二区三区在线观看| 久久久久国产一级毛片高清牌| 精品日产1卡2卡| 看片在线看免费视频| 黄色成人免费大全| 日本五十路高清| av中文乱码字幕在线| 91精品三级在线观看| 大型av网站在线播放| 热99re8久久精品国产| 黄色a级毛片大全视频| 亚洲欧美精品综合一区二区三区| 日韩成人在线观看一区二区三区| 国产一区二区在线av高清观看| 久久久久九九精品影院| 亚洲精品在线美女| 涩涩av久久男人的天堂| 一边摸一边抽搐一进一小说| 色播在线永久视频| 日本五十路高清| 久久九九热精品免费| 咕卡用的链子| 久久精品国产亚洲av香蕉五月| tocl精华| 欧美 亚洲 国产 日韩一| 午夜a级毛片| 国产成人av教育| 欧美性长视频在线观看| a级片在线免费高清观看视频| 亚洲人成电影免费在线| 一边摸一边抽搐一进一出视频| 亚洲av美国av| www国产在线视频色| 国产亚洲精品久久久久久毛片| 露出奶头的视频| 国产亚洲欧美精品永久| 亚洲五月色婷婷综合| 女生性感内裤真人,穿戴方法视频| 免费不卡黄色视频| 国产精品99久久99久久久不卡| 亚洲七黄色美女视频| 法律面前人人平等表现在哪些方面| 亚洲成人久久性| 中文字幕另类日韩欧美亚洲嫩草| 看免费av毛片| av有码第一页| 国产成人啪精品午夜网站| 欧美一区二区精品小视频在线| 久久中文字幕人妻熟女| 久久热在线av| 黄色女人牲交| 精品无人区乱码1区二区| 欧美激情久久久久久爽电影 | 国产亚洲欧美在线一区二区| 欧美最黄视频在线播放免费 | 午夜免费观看网址| x7x7x7水蜜桃| 亚洲欧美激情在线| av网站免费在线观看视频| 欧美大码av| 免费看a级黄色片| 一边摸一边抽搐一进一小说| 啦啦啦在线免费观看视频4| 国产欧美日韩一区二区三区在线| 亚洲人成电影免费在线| 亚洲av五月六月丁香网| 在线天堂中文资源库| 很黄的视频免费| 欧美性长视频在线观看| 国产成人精品在线电影| 在线播放国产精品三级| 国产深夜福利视频在线观看| 啪啪无遮挡十八禁网站| 欧美成人午夜精品| 久久久水蜜桃国产精品网| 免费女性裸体啪啪无遮挡网站| 99国产精品免费福利视频| 一区福利在线观看| 在线观看免费日韩欧美大片| 午夜福利影视在线免费观看| 国产欧美日韩精品亚洲av| 久久精品国产99精品国产亚洲性色 | 啦啦啦在线免费观看视频4| 正在播放国产对白刺激| 欧美日韩一级在线毛片| 久久伊人香网站| 久久久久国产精品人妻aⅴ院| 自线自在国产av| 亚洲久久久国产精品| 精品福利永久在线观看| 最新在线观看一区二区三区| av国产精品久久久久影院| 午夜两性在线视频| 女人高潮潮喷娇喘18禁视频| www.999成人在线观看| a级毛片在线看网站| 日本vs欧美在线观看视频| 久久国产亚洲av麻豆专区| 天天躁狠狠躁夜夜躁狠狠躁| 又大又爽又粗| 久久热在线av| 一边摸一边做爽爽视频免费| 国产免费av片在线观看野外av| 免费搜索国产男女视频| 少妇的丰满在线观看| 亚洲精品久久午夜乱码| 久久久久久免费高清国产稀缺| 男女做爰动态图高潮gif福利片 | 中文字幕最新亚洲高清| 亚洲国产精品合色在线| 国产日韩一区二区三区精品不卡| 欧美黄色淫秽网站| 91九色精品人成在线观看| 亚洲成人免费电影在线观看| 日日干狠狠操夜夜爽| 757午夜福利合集在线观看| 人人妻人人爽人人添夜夜欢视频| 国产成人免费无遮挡视频| 国产三级黄色录像| 高清在线国产一区| 色在线成人网| 亚洲男人的天堂狠狠| 999久久久精品免费观看国产| 国产在线观看jvid| 无人区码免费观看不卡| 999久久久国产精品视频| 国产成人影院久久av| 黄网站色视频无遮挡免费观看| 精品久久久久久久毛片微露脸| 国产亚洲精品一区二区www| 国产精品一区二区三区四区久久 | 男女做爰动态图高潮gif福利片 | 人人妻人人澡人人看| 一级黄色大片毛片| 久久九九热精品免费| 精品福利观看| 午夜免费观看网址| 一级毛片女人18水好多| 欧美老熟妇乱子伦牲交| 久久久久久久久免费视频了| 亚洲一区中文字幕在线| 欧美激情极品国产一区二区三区| 亚洲精品av麻豆狂野| 成人亚洲精品一区在线观看| 99国产精品99久久久久| 亚洲激情在线av| 少妇的丰满在线观看| 欧美日韩亚洲高清精品| 99riav亚洲国产免费| 免费日韩欧美在线观看| 国产精品久久视频播放| 午夜激情av网站| 久久久久久久午夜电影 | 精品国内亚洲2022精品成人| 国产99白浆流出| 啦啦啦 在线观看视频| 国产亚洲欧美98| 一本大道久久a久久精品| 久久国产乱子伦精品免费另类| 高清黄色对白视频在线免费看| 国产片内射在线| 亚洲成人久久性| 叶爱在线成人免费视频播放| 操美女的视频在线观看| 欧美国产精品va在线观看不卡| 一a级毛片在线观看| 一级,二级,三级黄色视频| 欧美中文综合在线视频| 超碰成人久久| 叶爱在线成人免费视频播放| 国产av一区在线观看免费| 妹子高潮喷水视频| av有码第一页| 成年女人毛片免费观看观看9| 成年版毛片免费区| 宅男免费午夜| 国产av又大| 9191精品国产免费久久| 欧美在线黄色| 国产成人影院久久av| 日日爽夜夜爽网站| 日本欧美视频一区| 嫁个100分男人电影在线观看| 啦啦啦在线免费观看视频4| 精品人妻1区二区| 国产亚洲欧美在线一区二区| 久久人妻av系列| 欧美性长视频在线观看| 国产高清激情床上av| 男人操女人黄网站| 亚洲专区字幕在线| 99热只有精品国产| 最新在线观看一区二区三区| 女人被躁到高潮嗷嗷叫费观| www.www免费av| 一级黄色大片毛片| 午夜福利在线观看吧| 亚洲自偷自拍图片 自拍| 国产精品久久视频播放| 九色亚洲精品在线播放| av网站免费在线观看视频| 咕卡用的链子| 日本a在线网址| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩 欧美 亚洲 中文字幕| 国产精品av久久久久免费| 国内久久婷婷六月综合欲色啪| 国产高清视频在线播放一区| 欧美在线一区亚洲| 成人亚洲精品av一区二区 | 黄色怎么调成土黄色| 波多野结衣高清无吗| 亚洲av第一区精品v没综合| 国产精品久久电影中文字幕| 欧美日韩黄片免| 成在线人永久免费视频| 欧美中文日本在线观看视频| 麻豆成人av在线观看| 99国产综合亚洲精品| 自拍欧美九色日韩亚洲蝌蚪91| 久久香蕉精品热| 90打野战视频偷拍视频| 国产激情久久老熟女| 天天躁夜夜躁狠狠躁躁| 99在线视频只有这里精品首页| 成人特级黄色片久久久久久久| 一边摸一边抽搐一进一小说| 亚洲一区二区三区色噜噜 | 91精品三级在线观看| 淫秽高清视频在线观看| 亚洲人成电影免费在线| 可以免费在线观看a视频的电影网站| 国产乱人伦免费视频| 亚洲成人国产一区在线观看| 中文字幕人妻熟女乱码| 高清欧美精品videossex| 国产免费av片在线观看野外av| 操出白浆在线播放| 伊人久久大香线蕉亚洲五| 别揉我奶头~嗯~啊~动态视频| 久久婷婷成人综合色麻豆| 久久久国产成人免费| 黄网站色视频无遮挡免费观看| 久久天躁狠狠躁夜夜2o2o| 亚洲午夜精品一区,二区,三区| 免费搜索国产男女视频| 手机成人av网站| 久久久久久久精品吃奶| 午夜福利影视在线免费观看| 日本免费a在线| 搡老熟女国产l中国老女人| 亚洲精品久久成人aⅴ小说| 国产精品久久电影中文字幕| 国产蜜桃级精品一区二区三区| 国产精品永久免费网站| 黄色视频不卡| 国产亚洲精品综合一区在线观看 | 黄色视频,在线免费观看| 看黄色毛片网站| 国产精品日韩av在线免费观看 | 日日爽夜夜爽网站| 啦啦啦 在线观看视频| 看免费av毛片| 99国产精品一区二区蜜桃av| 9热在线视频观看99| 久久亚洲真实| 国产亚洲精品综合一区在线观看 | 日本精品一区二区三区蜜桃| 男男h啪啪无遮挡| 亚洲av第一区精品v没综合| 老司机在亚洲福利影院| 麻豆av在线久日| 在线观看66精品国产| 欧美日韩一级在线毛片| 男人舔女人的私密视频| 久久久精品国产亚洲av高清涩受| 日韩精品青青久久久久久| 精品久久久久久成人av| 精品一品国产午夜福利视频| 亚洲欧洲精品一区二区精品久久久| 美女午夜性视频免费| 极品人妻少妇av视频| 999久久久精品免费观看国产| 巨乳人妻的诱惑在线观看| 国产一区二区激情短视频| 久久草成人影院| 一区二区三区激情视频| а√天堂www在线а√下载| 黄片播放在线免费| 久久人人97超碰香蕉20202| 成人免费观看视频高清| 黄色丝袜av网址大全| 亚洲精品一区av在线观看| 一级,二级,三级黄色视频| 免费av中文字幕在线| 麻豆成人av在线观看| 日韩欧美免费精品| 欧美人与性动交α欧美软件| 亚洲七黄色美女视频| 国产精品乱码一区二三区的特点 | 三上悠亚av全集在线观看| 日本vs欧美在线观看视频| 欧美在线黄色| 高清在线国产一区| 午夜激情av网站| 亚洲精品美女久久久久99蜜臀| 国产成人av教育| 成年版毛片免费区| 亚洲精品美女久久久久99蜜臀| 久久久精品国产亚洲av高清涩受| 变态另类成人亚洲欧美熟女 | 国产片内射在线| 三上悠亚av全集在线观看| 国产一区二区激情短视频| 亚洲性夜色夜夜综合| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美激情在线| 亚洲精品中文字幕一二三四区| 看黄色毛片网站| 久热这里只有精品99|