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

    MYOD1 inhibits avian adipocyte differentiation via miRNA-206/KLF4 axis

    2021-12-17 11:52:38ZhengWangQiangsenZhaoXiaoqinLiZhongtaoYinSiruiChenSenWuNingYangandZhuochengHou

    Zheng Wang,Qiangsen Zhao,Xiaoqin Li,Zhongtao Yin,Sirui Chen,Sen Wu,Ning Yang and Zhuocheng Hou*

    Abstract

    Background:A considerable number of muscle development-related genes were differentially expressed in the early stage of avian adipocyte differentiation.However,the functions of them in adipocyte differentiation remain largely known.In this study,the myoblast determination protein 1(MYOD1)was selected as a representative of muscle development.We investigated its expression,function,and regulation in avian adipocyte differentiation.

    Results:The expression of MYOD1 decreased significantly in the early stage of avian adipocyte differentiation.CRIS PR/Cas9-mediated deletion of MYOD1 induced adipocyte differentiation,whereas over-expression of MYOD1 inhibited adipogenesis.The mRNA-seq data showed that MYOD1 could perturb the lipid biosynthetic process during differentiation.Our results showed that MYOD1 directly up-regulates the miR-206 expression by binding the upstream 1200 bp region of miR-206.Then,over-expression of miR-206 can inhibit the adipogenesis.Furthermore,MYOD1 affected the expression of endogenous miR-206 and its target gene Kruppel-like factor 4(KLF4),which is an important activator of adipogenesis.Accordingly,the inhibition of miR-206 or over-expression of KLF4 could counteract the inhibitory effect of MYOD1 on adipocyte differentiation.

    Conclusions:Our results establish that MYOD1 inhibits adipocyte differentiation by up-regulating miR-206 to suppress the KLF4 expression.These findings identify a novel function of MYOD1 in adipocyte differentiation,suggesting a potential role in body-fat distribution regulation.

    Keywords:Adipocyte differentiation,Avian,CRISPR/Cas9,miR-206/KLF4 axis,MYOD1

    Background

    Adipocytes are unique in the quantity of lipids that they can store,the rapid release of these calories,and protein for use by other organs,which can profoundly affect our health[1].Avian has been utilized as a good animal model for studying basic adipogenesis mechanisms,which can be directly used for the genetic improvement of avian fat deposition[2].Adipogenesis is driven by an increase in adipocyte cell size(hypertrophy)or number(hyperplasia)[3].Adipocyte differentiation is regulated by an elaborate network of transcription factors,and understanding the underlying transcriptional networks is relevant and timely both from a basic and medical research perspective[4].

    Several studies and observations showed that fat deposits in the muscles will cause the loss of muscle quality and are more likely to induce metabolic diseases in various animals[5–8].Adipose and muscle tissues originate from mesenchymal stem cells(MSCs)[9–12],which are the multipotent and relevant targets for therapies aiming to enhance tissue regeneration[13].In response to lineage-specific inducers,MSCs from different depots can differentiate into many different,mutually exclusive lineages,including the adipocyte and myoblast lineages[14–16].It has been shown that peroxisome proliferator-activated receptor γ (PPARγ)and myoblast determination protein 1(MYOD1)are the master regulators of adipogenesis and myogenesis,respectively[17,18].At the same time,the MYOD1-driven and PPARγdriven differentiation programs are mutually exclusive[19–21].However,the functional contributions of muscle development genes to adipocyte differentiation remain largely unexplored.

    The microRNA-206(miR-206)is a vertebrate-specific miRNA,which belongs to the skeletal muscle-specific myomiR family(myomiRs)along with miR-1,miR-133a,miR-133b.It has been confirmed to be involved in the pathogenesis of many diseases,including heart failure,chronic obstructive pulmonary disease,and various types of cancers[22,23].In chicken,miR-206 is significantly associated with broiler birthweight[24].Various TFs essential for skeletal muscle development have been shown to regulate miR-206 expression during myogenic differentiation,such as MYOD1,myogenic factor 5(MYF5),myogenin(MYOG),and myocyte enhancer factor 2C(MEF2C)[25–27].In mammals,miR-206 promotes apoptosis,induces cell cycle arrest,and inhibits cell migration and adipocyte differentiation by targeting c-MET and its downstream PI3K/AKT pathway[28,29].However,the role of miR-206 in avian adipocyte differentiation still remains unclear.

    Kruppel-like factor 4(KLF4)is an evolutionarily conserved zinc finger-containing transcription factor which is also known for being one of four key factors required for inducing pluripotent stem cells.KLF4 has been shown to be induced very early following the induction of adipogenesis,and acted as an activator of adipogenesis by inducing CCAAT enhancer binding proteins β(C/EBPβ)expression[30],which then trans-activate C/EBPα and PPARγ[31,32].Studies have shown that knocking down KLF4 suppresses adipocyte differentiation[30].More recent studies have shown that the ability of A2b adenosine receptor(ADORA2B)to block adipocyte development has been shown to be dependent on its ability to modulate KLF4 expression[33].

    Our previous work found many muscle development genes involved in the early regulation of adipocyte differentiation[34].In this study,MYOD1,a representative of muscle development,was significantly down-regulated in the early stage of avian adipocyte differentiation.Studies of loss-of-function and gain-of-function demonstrated that MYOD1 is a key repressor of adipocyte differentiation by interacting with miR-206/Kruppel-like factor 4(KLF4)axis.The mRNA-seq showed that overexpression of MYOD1 in adipocytes inhibits the expression of most lipid biosynthesis genes and also promotes the expression of some myogenic genes.Our findings imply that a novel function of MYOD1 in adipocyte differentiation.Thus,affecting MYOD1 might represent a viable strategy to improve fat ratio to muscle in avian.

    Materials and methods

    Plasmid construction

    MYOD1-knock-out and knock-in plasmids:Using the MYOD1 sequence obtained from the NCBI database(Accession:NC_006092.5),we designed gRNA sequences targeting exon1 of MYOD1,known as sgRNA1: CGACCCGTGCTTCAACACGT and sgRNA2:GCGGCTCAGCAAGGTCAACG.We synthesized the oligo-DNAs corresponding to these gRNAs.We annealed them to a T7 promoter-driven Cas9 and to a U6 promoter-driven gRNA vector in order to obtain two gRNA-expressing plasmids.In order to construct the MYOD1-over-expression vector,the full-length coding sequence of MYOD1(NCBI Reference Sequence:NM_204214.2)was amplified from chicken subcutaneous adipose cDNA by PCR,and cloned into the CMV promoter-driven piggyBac and an EF1α promoter-driven GFP plasmid by replacing GFP using EcoRI and SalI(New England Biolabs,Ipswich,MA,USA).

    pmirGLO dual-luciferase reporters:The 3′-UTR fragment of KLF4 (NCBI Reference Sequence:XM_004949369.3)containing the binding sites were amplified by PCR from chicken subcutaneous adipose cDNA and then cloned into pmirGLO vector.The mutant vectors were constructed by PCR mutagenesis.Six seed sequences were successfully mutated from CATTCC to GTGAAG for the KLF4-3′-UTR vector.

    Gene over-expression vector:The MYOD1 and KLF4 over-expression vector was constructed according to the user manual of the Easy Ligation Kit(Sidansai,Shanghai,China).MYOD1 and KLF4 coding sequence(NCBI Reference Sequence:NM_204214.2 and XM_004949369.3)were amplified from chicken subcutaneous adipose cDNA by PCR.The PCR product was cloned into the pcDNA3.1 vector.The successful MYOD1 and KLF4,over-expression vector,was confirmed by DNA sequencing.

    miR-206 promoter-reporter plasmid:A 1876 bp fragment of the miR-206 promoter was isolated by PCR using the primers listed in Table S1.The PCR product was digested with KpnI and SmaI restriction enzymes(Takara,Otsu,Japan),then the insertion was ligated into the pGL4.10 vector(Promega,Madison,WI,USA)to create the expression vector pGL4.10_?1876bp.The pGL4.10_?1876 bp was subjected to be sequenced,and this construct was used as a template,and pGL4.10_-1234 bp was isolated by PCR.All cloning plasmids were confirmed by sequencing.

    RNA extraction,cDNA synthesis,and quantitative realtime PCR

    According to the manufacturer’s instruction,the total RNA was isolated from the cells using RNAiso reagent(Takara,Otsu,Japan).According to the manufacturer’s manual,the reverse transcription reaction for mRNA was performed with PrimeScript RT reagent Kit(Perfect Real-Time)(Takara,Otsu,Japan).The reverse transcription reaction for miRNA was using miRNA First-Strand cDNA Synthesis SuperMix(Transgen,Beijing,China).The specific qRT-PCR Primer of mRNA and miRNA were designed using Primer 3 software(version 0.4.0,Howard Hughes Medical Institute).Primer sets are listed in Table S1.The qPCR reaction,which contained KAPA SYBR FAST qPCR Kit(KAPA Biosystems,MA,USA),was carried out in ABI-7500 PCR machine(Applied Biosystems,MA,USA)following a standard thermal protocol as previously described[35].All reactions were run in triplicate.

    Cell culture

    A cell line of immortalized chicken preadipocytes(ICPs)[36]was cultured in DMEM/F12(Gibco,Gaithersburg,MD,USA)supplemented with 10% fetal bovine serum(Hyclone,Logan,UT,USA),and 0.2% penicillin/streptomycin(Invitrogen,Carlsbad,CA,USA).To induce ICPs differentiation,we added 160μmol/L sodium oleate(Sigma Life Science,St.Louis,MO,USA)to the medium[37].

    For MYOD1OEand MYOD1KOcell selection,ICPs were seeded in 6-well plates for further transfection using Lipofectamine 3000(Invitrogen,Carlsbad,CA,USA).After a 48-h recovery period,the cells were supplemented with 3μg/mL of puromycin(Sigma-Aldrich,MO,USA)in the culture medium for 12days until clone formation.Cells were harvested using 0.25% trypsin/EDTA(Gibco,Gaithersburg,MD,USA),and the cell density was calculated using a handheld automated cell counter(Millipore,Darmstadt,Germany).Single cells were plated in each well of a 96-well plate by limiting dilution and then cultured for 10 d in the cell culture medium.The medium was replaced every 4 d.Confluent cell colonies were propagated and genotyped by PCR and sequencing.

    Transfections

    Transfections were performed with Lipofectamine 3000 reagent(Invitrogen,Carlsbad,CA,USA)according to the manufacturer’s direction.Nucleic acids were diluted in OPTI-MEM Medium(Gibco,Gaithersburg,MD,USA).All experiments were carried out at least three times independently.

    Oil red O staining and quantification

    The cells were washed with PBS and fixed in 4% formaldehyde for 10 min.The cells were then stained with Oil-Red-O working solution(Solarbio,Beijing,China)according to the manufacturer’s manual.After another wash with PBS,the cell nuclei were counterstained with Hoechst 33342(Solarbio,Beijing,China).Images were acquired by inverted fluorescent microscope(Nikon)with 20× objective lens(200×magnification),and quantified by Image J software(National Institutes of Health,Bethesda,MD,USA).Particle number was quantified with the analysis particles functions in threshold single sections with size(pixel2)setting from 0.1 to 10,and circularity from 0.1 to 1.The Oil-Red-O dyes were then extracted in isopropanol solution containing 4% Nonidet P-40 and quantified by NanoDrop 2000C spectrophotometers(Thermo Fisher Scientific,San Jose,CA,USA)at 510nm.

    RNA oligonucleotides

    The miR-206 mimics,negative control(NC)mimic,miR-206 inhibitors,and NC inhibitor were all purchased from GenePharma(GenePharma,Shanghai,China).

    Dual-luciferase reporter assay

    For the promoter activity assays,ICPs were cotransfected with reporter plasmid and MYOD1 over-expression vector or control vector.The TK-Renilla reporter was also cotransfected to each sample as an internal control using the Lipofectamine 3000 reagent(Invitrogen,Carlsbad,CA,USA)in 48-well plates.The miRNA target verification assay was also performed in ICPs.Wild-type or mutant KLF4-3′-UTR dual-luciferase reporter(200ng)and miR-206 mimic or NC mimic(50nmol/L)were cotransfected into ICPs.After 48h transfection,cells were washed by PBS twice,and the activities of Firefly and Renilla luciferase were measured according to the manual of Luc-pair Duo-Luciferase Assay Kit 2.0(GeneCopoeia,Rockville,MD,USA).All the data were acquired by averaging the results from three independent repeats.

    Western blot

    Cultured cells were washed with PBS and homogenized with RIPA buffer(Beyotime,Jiangsu,China)containing protease inhibitor cocktail(Beyotime,Jiangsu,China).Protein concentrations were determined using the BCA Protein Assay Kit(Beyotime,Jiangsu,China).Proteins were denatured and subjected to 10% polyacrylamide gel and transferred to methanol-activated PVDF membranes.Blots were probed using the primary antibodies:mouse anti-MYOD1(1:500;Santa Cruz Biotechnology,USA,Cat#sc-377460),rabbit anti-KLF4(1:500;Bioss,Beijing,China,Cat#52850R),mouse anti-GAPDH,(1:5000;Bioworld,St Louis Park,MN,USA,Cat#MB001),overnight at 4°C.After 1 h incubation with anti-mouse or antirabbit HRP-conjugated second antibody(1:5000,Bioss,Beijing,China,Cat#40296G,40295G).Immunodetection was performed using enhanced chemiluminescence(ECL)Western blotting substrate(Beyotime,Jiangsu,China),and was detected with FluoChem R imaging system(ProteinSimple,CA,USA).

    RNA-seq analysis

    Raw reads were trimmed to remove adapters and lowquality reads,with Trimmomatic(version 0.39)[38].Trimmed reads were mapped to the chickenreference genome(Ensembl release 100:ftp://ftp.ensembl.org/pub/release100/fasta/gallus_gallus/dna/Gallus_gallus.GRCg6a.dna.toplevel.fa.gz) using HISAT2[39].Read counts for each gene were calculated using Stringtie(v.2.1.2)and normalized by library sequencing depth using the R package DESeq2(v.1.28.1)after filtering the gene with low expression[40].We used the DEGSeq2(v.1.28.1)package to identify DEGs between MYOD1NC,MYOD1KO,and MYOD1OEcells at different days(day 0 and day 5).Therefore,samples were excluded from further analysis due to their low global Pearson correlation with the other repeat samples(R2<0.95).Genes with|log2FC|≥0.585(or≥1)and the Benjamini&Hochberg(BH)adjusted P-value(adjusted-P value)<0.05 were considered as differentially expressed genes.

    Functional enrichment and prediction of miRNA target genes

    Genes were annotated with gene Symbols from the Uniprot database for functional annotation.Gene Ontology(GO)analysis of the enriched genes was performed using the web-based Metascape[41](a gene annotation&analysis resource,http://metascape.org/gp/index.html#/main). Putative miRNA targets for miR-206 were predicted by online software,TargetScan (version 7.2,http://www.targetscan.org/),miRDB(http://mirdb.org/) as well as miRTarBase(http://mirtarbase.cuhk.edu.cn/) to choose target genes for validation.

    Statistical analysis

    Each experiment was repeated three times,and all results are represented as the mean±SD.Three independent sample t-test was used to perform the statistically significant difference between groups.The level of significance was presented as?(P<0.05).

    Results

    MYOD1 is a repressor of adipocyte differentiation

    In our previous study[34],we found that the expression of MYOD1 is significantly down-regulated after the beginning of Pekin duck adipocyte differentiation(Fig.1a).To further understand the relationship between MYOD1 and adipocyte differentiation,we detected its expression in adipogenic differentiation of ICPs by qPCR(Fig.1b)and Western blot(Fig.1c).MYOD1 was also significantly down-regulated its expression after adipocyte differentiation in chicken.These results indicate that MYOD1 is involved in the avian adipocyte differentiation process.

    Fig.1 MYOD1 is down-regulated in the early stage of avian adipocyte differentiation.a,b MYOD1 mRNA levels in different time points during the differentiation of Pekin duck subcutaneous preadipocytes(mRNA-seq)and ICPs(qPCR).c MYOD1 protein levels at different time points during ICPs differentiation were determined by Western blotting.Data are shown as mean±SD of three biological replicates.The t-test was used to analyze the statistical differences between groups.*,P<0.05

    To further determine the roles of MYOD1 in avian adipocyte differentiation,we first performed gain-offunction experiments by using piggyBac delivery[42]of MYOD1 into ICPs.As shown in Fig.2a and Fig.2b,transduced MYOD1 clone can significantly induce MYOD1 expression relative to cells transduced with a control vector.Strikingly,over-expression of MYOD1 in ICPs(MYOD1OE)significantly blocked adipogenesis(Fig.2c,d),as shown by oil red O staining of neutral lipids(Fig.S1a,b).Adipocyte markers,such as PPARγ,AFABP,C/EBPα,and C/EBPβ (Fig.2e),were significantly decreased in MYOD1OEcells.

    Fig.2 MYOD1 is a repressor of avian adipocyte differentiation.a,b MYOD1OE cell significantly promoted MYOD1 mRNA and protein expression in ICPs.c Representative images of MYOD1OE cells reduced the lipid droplet formation by oil red O staining on day 3.d Comparison of the lipid droplet content of MYOD1OE and MYOD1NC cells obtained by oil red O staining and extraction methods.e mRNA levels of adipocyte genes PPARγ,A-FABP,C/EBPα,and C/EBPβ were analyzed with qPCR.f Upper part:schematic diagram of MYOD1 exon1 region and the two targeting loci of MYOD1 sgRNA(red).Lower part:DNA sequence map around the targeting locus of the cleaved band amplified from ICPs transfected with both sgRNAs.g,h MYOD1KO significantly reduced MYOD1 mRNA and protein expression in ICPs.i Representative images of MYOD1KOcells promoted the lipid droplet formation by oil red O staining on day 3.j Comparison of the lipid droplet content of MYOD1OE and MYOD1NC cells obtained by oil red O staining and extraction methods.k mRNA levels of adipocyte marker genes were analyzed with qPCR.l Compare the fold change of known lipid biosynthesis genes in MYOD1OE,MYOD1NC,and MYOD1KOcells on day 5 compared to day 0.Data are shown as mean±SD of three biological replicates.The t-test was used to analyze the statistical differences between groups.*,P<0.05

    We further sought confirmation of MYOD1 antiadipogenic activity through loss-of-function studies,in which we would predict enhanced adipose conversion.We transfected ICPs with two single guide RNAs(sgRNAs)targeting the exon 1 of MYOD1 and a Cas9 vector.We established MYOD1-knock-out ICP Clones(MYOD1KO)with 260-bp frame-shifting deletions in MYOD1(Fig.2f),indicating that both gRNAs work efficiently with Cas9 to edit MYOD1.The mutation of MYOD1 also dramatically down-regulated its expression relative to cells transfected with an only Cas9 vector(Fig.2g,h).MYOD1KOcells demonstrated enhanced adipogenic potential,including greater lipid accumulation(Fig.2i,j;Fig.S2a,b),and increased the expression of adipocyte marker genes(Fig.2k).

    MYOD1 perturb the lipid biosynthetic process

    To further understand the effect of MYOD1 on adipogenesis,we performed mRNA-Seq experiments in MYOD1OE,MYOD1NC,MYOD1KOcells prior to differentiation(day 0)and day 5 after differentiation.In total,886 and 1335 differentially expressed genes(DEGs)were up-regulated and down-regulated in MYOD1OEcells compared to MYOD1NCcells on day0,respectively(Table S2,fold change>2;adjusted-P value<0.05).GO analyses for DEGs up-regulated in MYOD1OEcells are enriched for development processes such as blood vessel development,muscle structure development,and positive regulation of muscle tissue development.In contrast,DEGs downregulated in MYOD1OEcells are enriched for extracellular matrix organization,metabolism process,and cell morphogenesis involved in differentiation (Table 1,Table S3).Also,256 and 971 DEGs were up-regulated and down-regulated in MYOD1KOcells compared to MYOD1NCcells on day 0,respectively(Table S2,fold change>2;adjusted-P value<0.05).As expected,GO analysis for genes down-regulated in MYOD1KOcells also are enriched for striated muscle cell differentiation and muscle system process.In contrast,genes up-regulated in MYOD1KOcells are enriched inthe adipogenesis related pathway,including cell morphogenesis in differentiation,regulation of MAPK cascade,and positive regulation of lipid metabolic process(Table 1,Table S3).These results indicate that although differentiation has not yet begun,the knock-in or knockout of MYOD1 alone seems to have affected preadipocyte characteristics.

    Table 1 GO enrichment analysis of DEGs in MYOD1OE,MYOD1NC,and MYOD1KO cells

    To better understand the effect of MYOD1 on adipogenesis,we re-analyzed the DEGs before and after adipogenic differentiation in MYOD1KO,MYOD1NC,and MYOD1OEcells.A total of 2457,781,and 1636 DEGs were significantly up-regulated on day 5 of differentiation of MYOD1OE,MYOD1NC,and MYOD1KOcells compared to day 0,respectively(Table S4,fold change>1.5;adjusted-P value<0.05).The up-regulated DEGs of the three adipocyte lines on day 5 were significantly enriched in cell division and DNA replication,which is an essential step in adipocyte differentiation.At the same time,down-regulated DEGs are also significantly enriched in muscle structure development(Table 1,Table S5),indicating that muscle development-specific genes are also considerably suppressed during chicken adipocyte differentiation.Notably,we found that the lipid biosynthetic process was only enriched in the DEGs up-regulated on day 5 in MYOD1KOand MYOD1NCcells,but not in MYOD1OEcells(Table 1,Table S5).Furthermore,to test whether MYOD1 over-expression impacts the lipid biosynthetic process in a statistical threshold-independent manner,we directly compared expression fold-changes of known lipid biosynthesis genes(from Genecards database:Pathcards:Fatty Acyl-CoA Biosynthesis&Triglyceride Biosynthesis)during adipocyte differentiation and found that most lipid biosynthesis genes tended to be up-regulated on day 5 of MYOD1KOcells,but down-regulated in MYOD1OEcells(Fig.2l,Table S6),which consistent with observations from our phenotype-profiling experiments and indicating that MYOD1 is a repressor for the adipogenesis.

    The expression of miR-206 can be significantly induced in MYOD1OEcells

    MYOD1 has been confirmed to induce highly conserved myomiRs expression by binding to miRNA upstream regions,including miR-1,miR-133a,miR-133b,and miR-206,and these myomiRs are widely directly involved in the inhibition of other pathways[43–47].qPCR assay found that MYOD1OEcells expressed significantly higher levels of miR-1,miR-133a,miR-133b,and miR-206(Fig.3a),and miR-206 has the highest expression fold change(about 40 fold).MYOD1KOcells also significantly reduce miR-206 expression(Fig.3b).Also,we found multiple binding sites of MYOD1 in the 2000 bp upstream region of gga-miR-206,especially within 1200 bp(Table S7).Thus,we hypothesize that miR-206 may be an essential mediator for MYOD1 to inhibit adipocyte differentiation.To validate the regulatory relationship between MYOD1 and miR-206 transcription,two distinct lengths of upstream regions(1876 and 1234bp)of the gga-miR-206 transcription start site were amplified and cloned into pGL4.10 vector to detect the promoter activity.After cotransfecting with pGL-TK and pcDNA3.1-MYOD1 into ICPs,both of pGL4.10_?1876bp and pGL4.10_?1234bp showed a significantly increasing promoter activity compared to cotransfected with pGL-TK and pcDNA3.1(P<0.05)(Fig.3c).Our results showed that the chicken MYOD1 can bind to 1200bp upstream region of gga-miR-206 and promote the transcription activity of miR-206 in ICPs,as previously reported[24].

    miR-206 is a downstream gene of MYOD1 that inhibits adipocyte differentiation

    We then examined the expression and function of miR-206 during adipocyte differentiation.The expression of miR-206,similar to the protein level of MYOD1,was significantly down-regulated after differentiation(Fig.3d),suggesting a synergistic relationship between them,and both of them were involved in this process.Therefore,we transfected miR-206 mimic and inhibitor into ICPs,respectively.Transfection-mediated gene transfer resulted in up to 840-fold elevation in the expression of miR-206(Fig.3e).Over-expression of miR-206 also significantly inhibits lipid droplet accumulation and reduced adipocyte marker genes expression(Fig.3f,g;Fig.S3a,b),whereas the inhibition of miR-206 promotes adipocyte differentiation(Fig.3h-j;Fig.S4a,b).Together,these results demonstrated that miR-206 could inhibit adipocyte differentiation.

    To examine whether MYOD1 targets miR-206 to regulate adipocyte differentiation,we transfected MYOD1OEcells with miR-206 inhibitor and knock-down of miR-206 could counteract the inhibition effect of MYOD1 over-expression on adipocyte differentiation(Fig.3k,l;Fig.S5a,b).Similarly,over-expression of miR-206 also significantly inhibited lipid accumulation in MYOD1KOcells(Fig.S6a-e).Together,these results suggest that the inhibitory effect of MYOD1 on adipocyte differentiation was achieved by its downstream gene miR-206.

    Fig.3 The inhibitory effect of MYOD1 on adipocyte differentiation was achieved by its downstream gene miR-206.a Over-expression of MYOD1 up-regulated miR-1,miR-133a,miR-133b,and miR-206 expression in ICPs.b MYOD1 knockout down-regulated miR-1,miR-133a,miR-133b,and miR-206 expression in ICPs.c MYOD1 over-expression promoted the relative luciferase activity of the pGL4.10_-1234bp and pGL4.10_-1876bp reporter in ICPs.d Relative miR-206 expression during ICPs differentiation.e Transfected with miR-206 mimic significantly promoted miR-206 expression in ICPs.f Over-expression of miR-206 reduced the lipid droplet formation by oil red O staining on day 3.g mRNA levels of adipocyte marker genes were analyzed with qPCR.h Transfected with miR-206 inhibitor significantly reduced miR-206 expression in ICPs.i Inhibition of miR-206 promoted the lipid droplet formation by oil red O staining on day 3.j mRNA levels of adipocyte marker genes were analyzed with qPCR.k Inhibition of miR-206 could counteract the inhibition effect of MYOD1 over-expression on adipocyte differentiation.l mRNA levels of adipocyte marker genes were analyzed with qPCR.Data are shown as mean±SD of three biological replicates.The t-test was used to analyze the statistical differences between groups.*,P<0.05

    KLF4 is a miR-206 target gene,functioning as an activator of adipocyte differentiation

    In order to explore the potential mechanism of miR-206 in regulating adipocyte differentiation,we performed bioinformatic analysis.Three bioinformatic tools(TargetS-can,miRDB,and miRTarBase)were employed to identify the candidate targets of miR-206(Table S8).All three programs predicted a total of 16 miR-206 target genes,such as KLF4,CCND2,and UTRN(Fig.4a).Among them,we noticed that KLF4,a key activator of adipocyte differentiation,contains the miR-206 binding site in its 3′-UTR(Fig.4b).This binding site is conserved among avian(Fig.4c),suggesting the biological relevance of miR-206 in regulating KLF4 expression in avian.mRNA-seq also showed that the fold change(day 5 vs.day 0)of KLF4 expression was the highest in the differentiated MYOD1KOcells,followed by the MYOD1NC,and the lowest in the MYOD1OEcells(Fig.4d).To validate whether KLF4 is the target gene of miR-206,the 3′-UTR of chicken KLF4 containing the wild-type or mutated miR-206-binding sites were cloned into pmirGLO vector,and the luciferase activity was found to be significantly decreased in ICPs cotransfected with the vector carrying the wild-typebinding site in the presence of miR-206 mimics,but not in ICPs carrying the mutated-binding site(Fig.4e).On the other hand,the luciferase activity was found to be significantly increased in ICPs cotransfected with the vector carrying the wild-type-binding site in the presence of miR-206 inhibitor,but not in ICPs carrying the mutatedbinding site(Fig.4f).In addition,over-expression of miR-206 inhibited the mRNA and protein level of KLF4(Fig.4g,h),and the inhibition of miR-206 promoted the mRNA and protein level of KLF4(Fig.4i,j).Therefore,the above results indicated that KLF4 is the miR-206 target gene.

    Fig.4 KLF4 is the miR-206 target gene,functioning as an activator of adipocyte differentiation.a Overlap of three miRNA target bioinformatic prediction algorithms.b The potential binding site of miR-206 in the KLF4 mRNA 3′-UTR.c The potential binding site(red)of miR-206 in the KLF4 mRNA 3′-UTR is highly conserved among vertebrates.d Compare the fold change of KLF4 in MYOD1OE,MYOD1NC,and MYOD1KO cells at day 5 compared to day 0.e,f Dual-luciferase reporter assay indicated that miR-206 could bind to the predicted binding site of the KLF4 mRNA 3′-UTR.g,h miR-206 over-expression inhibited KLF4 mRNA and protein expression in ICPs.i,j miR-206 inhibition promoted KLF4 mRNA and protein expression in ICPs.k,l KLF4 mRNA and protein levels at different time points during ICPs differentiation.m Over-expression of KLF4 could counteract the inhibition effect of miR-206 over-expression on adipocyte differentiation.n mRNA levels of adipocyte marker genes were analyzed with qPCR.Data are shown as mean±SD of three biological replicates.The t-test was used to analyze the statistical differences between groups.*,P<0.05

    During adipocyte differentiation,KLF4 mRNA has gradually up-regulated its expression until day 3(Fig.4k).However,the protein level for KLF4 reached the highest on the day 1 after differentiation and disappeared on day 3(Fig.4l).Our resultsalso showed that KLF4 was induced very early following the induction of adipogenesis,and function as an early regulator of adipocyte development has been attributed to its capacity to induce C/EBPβ expression as the previous report[30].Furthermore,we cotransfected miR-206 and pcDNA3.1-KLF4 into ICPs,and KLF4 over-expression can counteract the inhibition effect of miR-206 on adipocyte differentiation(Fig.4m,n;Fig.S7a,b).Therefore,KLF4 is the miR-206 target gene,which can function as an activator of adipocyte differentiation.

    Based on the above results,we can deduce MYOD1 inhibits adipocyte differentiation by inhibiting the expression of KLF4.Western blot show MYOD1KOcells expressed significantly higher protein levels of KLF4,while MYOD1OEcells expressed lower protein levels of KLF4(Fig.5a).Finally,over-expression of KLF4 also significantly promoted adipocyte differentiation in MYOD1OEcells(Fig.5b,c;Fig.S8a,b).Also,over-expression of MYOD1 could significantly inhibit the expression of KLF4’s target genes C/EBPβ,while inhibiting miR-206 or over-expression of KLF4 in MYOD1OEcells will significantly increase C/EBPβ expression.

    Fig.5 MYOD1 affected KLF4 expression.a Western blot analysis of the protein levels of KLF4 in MYOD1OE,MYOD1NC,and MYOD1KO cells.b mRNA levels of adipocyte marker genes were analyzed with qPCR.c Representative images of KLF4 over-expression promoted the lipid droplet formation by oil red O staining(red).d Model of the MYOD1-mediated regulatory pathway for adipocyte differentiation.Data are shown as mean±SD of three biological replicates.The t-test was used to analyze the statistical differences between groups.*,P<0.05

    Discussion

    Avian are important farm animals throughout the world,producing eggs and high-quality meat for humans.Muscle mass and fat content are both important traits for meat-producing chickens but attaining an appropriate muscle and fat ratio is an excellent challenge for the broiler industry[48].Discovering functional genes that can simultaneously regulate muscle and adipocyte development is a key step to improve these traits[49].Previous and current transcriptomic data have shown that many muscle development genes were differentially expressed during adipocyte differentiation.Among these,our study investigated MYOD1 inhibited avian adipocyte differentiation via the miRNA-206/KLF4 axis(Fig.5d).We have also identified a negative regulation between MYOD1 and lipid biosynthesis genes during adipogenesis.Previous studies reported that MYOD1 expression in brown fat is significantly reduced,and it inhibits brown fat development through PRDM16[50].A lineagetracing study reveals that the MYOD1 lineage does not give rise to brown adipocytes,indicating a role of MYOD1 in myogenic cell fate switch in the common progenitors that give rise to both myoblasts and brown preadipocytes[20].Moreover,the loss of MYOD1 also facilitates the adipogenic trans-differentiation of C2C12 myoblasts by the miR-133/IGF1R/PI3K/AKT signaling pathway[47].These findings revealed a novel function of MYOD1 in adipogenesis.

    MYOD1 is considered as a master regulator of myogenesis as its expression can induce myogenic differentiation in myoblasts,fibroblasts,and a variety of other cell types[18,51,52].In this study,muscle structure development was significantly enriched in the upregulated DEGs of MYOD1OEcells,and also enriched in down-regulated DEGs of MYOD1KOon day 0,which all demonstrate the strong transcriptional regulatory activity of MYOD1 on downstream genes.Whether the ectopic expression of MYOD1 in preadipocytes promoted the myogenesis has not been examined in this study.Still,this evidence makes us believe that the overexpression of MYOD1 can promote the transdifferentiation of preadipocytes into muscle cells.In addition,several co-culture experiments have revealed that adipogenesis is strongly inhibited by the presence of satellite cell-derived myofibres[8,53].Supportively,the transcriptomic analysis also suggested that MYOD1 would be a repressor of adipogenesis by inhibiting the expression of lipid biosynthesis genes.Similarly,myogenesis was wholly blocked in both the MYOD1/MYF5 and MYOD1/IGF2 double knockout mice.Both results showed potential functions of MYOD1 on adipogenesis[54,55].On the other hand,myoblasts transdifferentiate into mature adipocytes by ectopic expression of adipogenic transcription factors under conditions permissive for adipogenesis[14].The mutual exclusion of the two lineage-specific transcription factors balances and determines the developmental separation of fat and muscle tissue.

    The current fast-growing and high-energy diet makes the chicken more likely to show myopathy,such as white striping,characterized by more fat in the breast muscles,resulting in meat with higher fat content and lower protein content[56,57].Combined with previously reported that MYOD1 can promote muscle development,increasing the expression of MYOD1 may be an effective strategy to treat these diseases.

    Finally,we show that miR-206 is an important mediator of MYOD1 induced inhibition of adipogenesis.Previous study has proven that miR-206 inhibits adipocyte differentiation by targeting c-MET and its downstream PI3K/AKT pathway in mammals[28,29].Through TargetScan and miRDB analysis,however,c-MET was not the target gene of gga-miR-206(Table S8),and the potential binding site of gga-miR-206 is also not in the 3′-UTR of c-MET(Table S9),which may be due to a divergence in evolution.Studies have shown that KLF4 expression is required for C/EBPβ expression observed in the early stages of adipogenesis.Indeed,in this study,we found that C/EBPβ and its downstream genes C/EBPα and PPARγ were significantly inhibited in MYOD1OEcells than in MYOD1NC,but they were increased after over-expression of KLF4 or inhibition of miR-206,suggesting a possible involvement of KLF4 and C/EBPβ in the inhibition of MYOD1 on C/EBPα and PPARγ transcription.However,given the over-expression of MYOD1 also promotes other myomiRs(miR-1,miR-133a,miR-133b)expression,we are implying that miRNA-206 may not be the sole target of MYOD1 to inhibit adipocyte differentiation.The inhibitory effect of MYOD1 on adipocyte differentiation should probably be accounted for by multiple targets.Therefore,identifying MYOD1’s other targets and understanding how it may affect these targets will warrant further investigation.

    Conclusions

    In summary,this study provides new insight into that MYOD1 also works as the repressor of adipocyte differentiation via miR-206/KLF4 axis in avian adipocyte model.Considering the significant role of adipocyte differentiation in the formation and function of adipose,clarification of the mechanism of MYOD1-mediated regulation of adipocyte differentiation is essential for exploring strategies for the treatment of metabolic disorders,including white striping.In combination with previous findings of the beneficial role of MYOD1 in muscle differentiation,we proposed that MYOD1 may be a crucial target for improving the ratio of muscle to fat in avian.

    Abbreviations

    ADORA2B:A2b adenosine receptor;A-FABP:Adipocyte-type fatty acid-binding protein;CCND2:Cyclin D2;CDS:Coding sequence;C/EBPs:CCAAT enhancer binding proteins;c-MET:MET proto-oncogene,receptor tyrosine kinase;CRIS PR:Clustered regularly interspaced short palindromic repeats;DEGs:Differentially expressed unigenes;FBS:Foetal bovine serum;GO:Gene Ontology;IGF1R:Insulin-like growth factor 1 receptor;IGF2:Insulin-like growth gactor 2;KLF4:Kruppel-like factor 4;MEF2C:Myocyte enhancer factor 2C;MYF5:Myogenic factor 5;MYOD1:Myogenic differentiation 1;MYOG:Myogenin;PPARγ:Peroxisome proliferator activated receptor gamma;PRDM16:PR/SET domain 16;RNA-seq:High-throughput sequencing of RNA;RT-qPCR:Quantitative real-time PCR;UTR:Untranslated region;UTRN:Utrophin

    Supplementary Information

    The online version contains supplementary material available at https://doi.org/10.1186/s40104-021-00579-x.

    Additional file 1.Table S1.Primer sequences used for the PCR and qPCR analysis.

    Additional file 2:Fig.S1.Statistics of the number and size of lipid droplets per cell.

    Additional file 3:Table S2.List of all DEGs in MYOD1OEcells vs.MYOD1NCcells and MYOD1KOcells vs.MYOD1NCon day 0.

    Additional file 4:Table S3.GO analysis of DEGs in MYOD1OEcells vs.MYOD1NCcells and MYOD1KOcells vs.MYOD1NCon day 0.

    Additional file 5:Table S4.List of all DEGs before(day 0)and after adipogenic differentiation(day 5)in MYOD1KO,MYOD1NCand MYOD1OEcells.

    Additional file 6:Table S5.GO analysis of DEGs before(day 0)and after adipogenic differentiation(day 5)in MYOD1KO,MYOD1NCand MYOD1OEcells.

    Additional file 7:Table S6.A list of the expression fold changes of known lipid biosynthesis genes in MYOD1KO,MYOD1NCand MYOD1OEcells.(Genes are selected from genecards database:Pathcards:Fatty Acyl-CoA Biosynthesis&Triglyceride Biosynthesis).

    Additional file 8:Table S7.Binding sites of MYOD1 in the 2000 bp upstream region of gga-miR-206.

    Additional file 9:Table S8.Predicted target genes of gga-miR-206.

    Additional file 1.:Table S9.Predicted binding miRNAs on 3′-UTR of c-Met gene.

    Acknowledgments

    We thank the Poultry Breeding Group of the College of Animal Science and Technology,Northeast Agricultural University,for providing the ICPs line.

    Authors’contributions

    ZCH conceived and designed the experimental plan.ZW,XQL and SRC collected samples and performed the experiments.ZW,QSZ and ZTY participated in bioinformatics analyses.ZW,QSZ,NY and ZCH drafted and revised this manuscript.All authors read and approved the final manuscript.

    Funding

    The work was supported by the National Waterfowl-Industry Technology Research System(CARS-42),National Nature Science Foundation of China(31972525,31572388),Beijing Municipal Science&Technology Commission(Z181100002418008),Key-Area Research and Development Program of Guangdong Province(2020B020222003).

    Availability of data and materials

    The data analyzed during the current study are available from the corresponding author on reasonable request.

    Declarations

    Ethics approval and consent to participate

    Collection of adipose samples for use in the described experiments were conducted following methods approved by the Animal Care and Use Committee of China Agricultural University(permit number:SYXK 2007-0023).

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Author details

    1National Engineering Laboratory for Animal Breeding and Key Laboratory of Animal Genetics,Breeding and Reproduction,MARA,College of Animal Science and Technology,China Agricultural University,Yuanmingyuan West Road No.2,Beijing 100193,China.2State Key Laboratory of Agrobiotechnology,College of Biological Sciences,China Agricultural University,Yuanmingyuan West Road No.2,Beijing 100193,China.

    Received:6 November 2020 Accepted:1 March 2021

    99在线视频只有这里精品首页| 亚洲av成人av| 免费少妇av软件| av有码第一页| 91国产中文字幕| 女人被狂操c到高潮| 黄色a级毛片大全视频| 一本大道久久a久久精品| 成人av一区二区三区在线看| 久久精品影院6| 国产高清有码在线观看视频 | 国产亚洲精品久久久久久毛片| 亚洲伊人色综图| 国产欧美日韩精品亚洲av| 中国美女看黄片| www国产在线视频色| 91精品国产国语对白视频| 日本黄色视频三级网站网址| 18禁国产床啪视频网站| 一边摸一边抽搐一进一小说| 可以免费在线观看a视频的电影网站| 一级毛片高清免费大全| 999精品在线视频| 免费av毛片视频| а√天堂www在线а√下载| 很黄的视频免费| 亚洲 欧美一区二区三区| 国产亚洲av高清不卡| 搡老岳熟女国产| 国产精品爽爽va在线观看网站 | 日韩有码中文字幕| 久久 成人 亚洲| 精品国产亚洲在线| 国产精品一区二区三区四区久久 | 如日韩欧美国产精品一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 精品久久久久久久人妻蜜臀av | 一级,二级,三级黄色视频| av片东京热男人的天堂| 精品高清国产在线一区| 自线自在国产av| 国产精品二区激情视频| 成人手机av| 久久伊人香网站| 少妇熟女aⅴ在线视频| 国产精品亚洲一级av第二区| 黄色成人免费大全| 电影成人av| 两个人视频免费观看高清| 性色av乱码一区二区三区2| 一夜夜www| 欧美日韩亚洲综合一区二区三区_| 宅男免费午夜| 欧美黑人精品巨大| 亚洲第一av免费看| 一进一出好大好爽视频| 长腿黑丝高跟| 国产亚洲欧美在线一区二区| 午夜a级毛片| 丝袜在线中文字幕| 久久人妻熟女aⅴ| 亚洲熟妇熟女久久| 天堂√8在线中文| 无遮挡黄片免费观看| 久久精品人人爽人人爽视色| 黄色片一级片一级黄色片| 日日爽夜夜爽网站| 激情视频va一区二区三区| 亚洲少妇的诱惑av| 悠悠久久av| 欧美 亚洲 国产 日韩一| 韩国av一区二区三区四区| 看免费av毛片| 高清在线国产一区| 日韩欧美国产在线观看| 欧美精品亚洲一区二区| 国产精品,欧美在线| 麻豆av在线久日| 国产激情久久老熟女| 亚洲午夜精品一区,二区,三区| 免费少妇av软件| а√天堂www在线а√下载| 欧美一级毛片孕妇| 久久国产精品影院| 国产精品一区二区三区四区久久 | 母亲3免费完整高清在线观看| 精品日产1卡2卡| 超碰成人久久| 久久精品aⅴ一区二区三区四区| 久久久久久久精品吃奶| 久久午夜综合久久蜜桃| 一进一出抽搐gif免费好疼| 变态另类丝袜制服| 日韩一卡2卡3卡4卡2021年| 丝袜美足系列| 巨乳人妻的诱惑在线观看| 最近最新中文字幕大全免费视频| 国产精品 欧美亚洲| 国产精品1区2区在线观看.| 久久精品aⅴ一区二区三区四区| 麻豆av在线久日| 欧美大码av| 夜夜看夜夜爽夜夜摸| 日韩欧美三级三区| 人人妻人人澡人人看| 成人精品一区二区免费| 国产成人系列免费观看| 日韩高清综合在线| 人妻丰满熟妇av一区二区三区| 亚洲精品国产区一区二| 成人亚洲精品一区在线观看| 日本撒尿小便嘘嘘汇集6| 好看av亚洲va欧美ⅴa在| 国产亚洲精品第一综合不卡| 国产精品久久视频播放| 欧美国产日韩亚洲一区| 看黄色毛片网站| 成人国产一区最新在线观看| 精品免费久久久久久久清纯| 多毛熟女@视频| 啦啦啦观看免费观看视频高清 | 一级作爱视频免费观看| www.精华液| 亚洲国产精品合色在线| 欧美黄色片欧美黄色片| 很黄的视频免费| 欧美乱码精品一区二区三区| 久久亚洲精品不卡| 视频在线观看一区二区三区| 老司机午夜福利在线观看视频| 欧美日本视频| 搞女人的毛片| 999久久久精品免费观看国产| 女人高潮潮喷娇喘18禁视频| 女性生殖器流出的白浆| 亚洲人成77777在线视频| 亚洲欧美激情综合另类| 国产1区2区3区精品| 久久婷婷成人综合色麻豆| 91麻豆精品激情在线观看国产| 亚洲av第一区精品v没综合| www国产在线视频色| 丝袜美足系列| 纯流量卡能插随身wifi吗| 人人妻人人爽人人添夜夜欢视频| a级毛片在线看网站| 女人被狂操c到高潮| 人妻丰满熟妇av一区二区三区| 亚洲五月色婷婷综合| av天堂在线播放| 午夜免费观看网址| 国产精品av久久久久免费| 91av网站免费观看| 淫秽高清视频在线观看| 免费看美女性在线毛片视频| 丝袜人妻中文字幕| 欧美一级毛片孕妇| 黑人操中国人逼视频| 中文字幕精品免费在线观看视频| 非洲黑人性xxxx精品又粗又长| 制服诱惑二区| 99久久精品国产亚洲精品| 两个人看的免费小视频| 成人特级黄色片久久久久久久| 18禁观看日本| 婷婷六月久久综合丁香| 中文字幕另类日韩欧美亚洲嫩草| 欧美色视频一区免费| 精品午夜福利视频在线观看一区| 午夜免费鲁丝| 一本大道久久a久久精品| 一进一出抽搐gif免费好疼| 欧美中文综合在线视频| 色av中文字幕| 中文字幕最新亚洲高清| 日日摸夜夜添夜夜添小说| 国产成人啪精品午夜网站| 欧美成人免费av一区二区三区| 亚洲伊人色综图| www日本在线高清视频| 国产1区2区3区精品| 一级毛片精品| 午夜久久久在线观看| 后天国语完整版免费观看| 国产精品精品国产色婷婷| 啦啦啦免费观看视频1| 日韩国内少妇激情av| 久久午夜综合久久蜜桃| 日韩欧美在线二视频| 高潮久久久久久久久久久不卡| 免费在线观看视频国产中文字幕亚洲| 久久久久亚洲av毛片大全| 91麻豆精品激情在线观看国产| 少妇裸体淫交视频免费看高清 | 午夜久久久在线观看| 国产亚洲精品久久久久久毛片| 国产一区二区在线av高清观看| 国产区一区二久久| 老司机在亚洲福利影院| 99国产精品免费福利视频| 热99re8久久精品国产| 国产免费男女视频| 免费高清视频大片| 精品一区二区三区四区五区乱码| 女同久久另类99精品国产91| 国内毛片毛片毛片毛片毛片| 久久国产精品人妻蜜桃| 韩国精品一区二区三区| 狠狠狠狠99中文字幕| 制服人妻中文乱码| 亚洲精品中文字幕在线视频| 亚洲九九香蕉| svipshipincom国产片| 国产亚洲av高清不卡| av有码第一页| 性欧美人与动物交配| 激情在线观看视频在线高清| 国产午夜精品久久久久久| 两人在一起打扑克的视频| 欧美不卡视频在线免费观看 | 91字幕亚洲| 制服丝袜大香蕉在线| 国产av一区在线观看免费| ponron亚洲| 欧美日韩亚洲综合一区二区三区_| 日韩欧美在线二视频| 老熟妇仑乱视频hdxx| 99国产综合亚洲精品| 欧美中文综合在线视频| 国产精品久久久久久亚洲av鲁大| av超薄肉色丝袜交足视频| 精品国内亚洲2022精品成人| 免费少妇av软件| 后天国语完整版免费观看| 一本大道久久a久久精品| 1024香蕉在线观看| 国产精品秋霞免费鲁丝片| 9热在线视频观看99| 日韩欧美三级三区| 精品日产1卡2卡| 每晚都被弄得嗷嗷叫到高潮| 久久亚洲真实| 咕卡用的链子| 免费高清视频大片| 丰满人妻熟妇乱又伦精品不卡| 日韩免费av在线播放| 亚洲一码二码三码区别大吗| 精品国内亚洲2022精品成人| 久久精品国产亚洲av香蕉五月| 日本一区二区免费在线视频| 美女国产高潮福利片在线看| 午夜福利免费观看在线| 国产成人啪精品午夜网站| 日韩免费av在线播放| 日韩av在线大香蕉| 亚洲午夜理论影院| 亚洲成人久久性| 国产精品香港三级国产av潘金莲| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品美女特级片免费视频播放器 | 叶爱在线成人免费视频播放| 国产一区二区激情短视频| 午夜福利一区二区在线看| 涩涩av久久男人的天堂| 亚洲成人久久性| 丝袜在线中文字幕| 变态另类成人亚洲欧美熟女 | 国产成人精品在线电影| 夜夜夜夜夜久久久久| 国产亚洲精品av在线| 亚洲专区中文字幕在线| 叶爱在线成人免费视频播放| 桃红色精品国产亚洲av| 亚洲精华国产精华精| netflix在线观看网站| 国内精品久久久久精免费| 最新美女视频免费是黄的| 亚洲天堂国产精品一区在线| 中文字幕人妻熟女乱码| 久久久久久久久久久久大奶| 免费在线观看亚洲国产| 日本三级黄在线观看| 中文字幕精品免费在线观看视频| 夜夜看夜夜爽夜夜摸| 无遮挡黄片免费观看| 欧美激情高清一区二区三区| 真人一进一出gif抽搐免费| 亚洲欧洲精品一区二区精品久久久| 看黄色毛片网站| 12—13女人毛片做爰片一| 国产一区二区三区综合在线观看| 国产精品秋霞免费鲁丝片| 中文亚洲av片在线观看爽| 最近最新中文字幕大全免费视频| 国产伦一二天堂av在线观看| 国产单亲对白刺激| 韩国精品一区二区三区| 中文字幕人妻熟女乱码| 日本黄色视频三级网站网址| 麻豆av在线久日| 欧美绝顶高潮抽搐喷水| 黑丝袜美女国产一区| 久久人人爽av亚洲精品天堂| 精品国产美女av久久久久小说| 狂野欧美激情性xxxx| 一进一出抽搐gif免费好疼| 18禁裸乳无遮挡免费网站照片 | 免费一级毛片在线播放高清视频 | 久久国产亚洲av麻豆专区| 宅男免费午夜| 欧美国产日韩亚洲一区| 精品一品国产午夜福利视频| 亚洲欧美精品综合久久99| 久久欧美精品欧美久久欧美| 成年人黄色毛片网站| 欧美不卡视频在线免费观看 | 午夜免费观看网址| 日韩欧美在线二视频| 久久人人精品亚洲av| 亚洲一卡2卡3卡4卡5卡精品中文| 一区二区三区激情视频| 一级片免费观看大全| 香蕉久久夜色| 亚洲av电影在线进入| 91老司机精品| 免费在线观看完整版高清| 亚洲av电影不卡..在线观看| 91成人精品电影| 午夜老司机福利片| 欧美+亚洲+日韩+国产| 91国产中文字幕| 老司机在亚洲福利影院| 色播亚洲综合网| 国产成+人综合+亚洲专区| av有码第一页| 黄色毛片三级朝国网站| 亚洲av美国av| 制服人妻中文乱码| 久久中文看片网| 亚洲第一电影网av| 美女午夜性视频免费| 无人区码免费观看不卡| 一进一出抽搐gif免费好疼| 两性夫妻黄色片| 黄色成人免费大全| a级毛片在线看网站| 亚洲最大成人中文| 亚洲av熟女| 国产伦一二天堂av在线观看| 中文字幕另类日韩欧美亚洲嫩草| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久午夜电影| 精品久久久久久久人妻蜜臀av | 久久这里只有精品19| 老汉色∧v一级毛片| 国产色视频综合| 一边摸一边做爽爽视频免费| a在线观看视频网站| 少妇 在线观看| 亚洲国产精品成人综合色| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美一区二区三区黑人| 亚洲自偷自拍图片 自拍| 免费高清在线观看日韩| 美国免费a级毛片| 嫩草影院精品99| 亚洲中文日韩欧美视频| 国产一区二区激情短视频| 国产99久久九九免费精品| 国产精品免费一区二区三区在线| 九色国产91popny在线| 国产av一区二区精品久久| 国产激情久久老熟女| 国产精品香港三级国产av潘金莲| 日韩欧美在线二视频| 日韩成人在线观看一区二区三区| 日韩高清综合在线| 亚洲av成人一区二区三| 国产精品野战在线观看| 亚洲中文字幕日韩| 国产av精品麻豆| 久久久久久久久免费视频了| avwww免费| 别揉我奶头~嗯~啊~动态视频| 精品福利观看| 精品国产乱码久久久久久男人| 夜夜躁狠狠躁天天躁| 久久久国产欧美日韩av| 麻豆国产av国片精品| 757午夜福利合集在线观看| 十八禁网站免费在线| 校园春色视频在线观看| 午夜福利高清视频| 午夜福利视频1000在线观看 | 男女午夜视频在线观看| 亚洲av美国av| 日韩欧美一区视频在线观看| 国产精品 国内视频| 啦啦啦韩国在线观看视频| 女人高潮潮喷娇喘18禁视频| 久久中文字幕一级| 亚洲性夜色夜夜综合| 精品午夜福利视频在线观看一区| 欧美不卡视频在线免费观看 | 欧美+亚洲+日韩+国产| 男女下面进入的视频免费午夜 | 成人免费观看视频高清| 午夜福利成人在线免费观看| 老司机靠b影院| 啪啪无遮挡十八禁网站| 精品一区二区三区av网在线观看| 色在线成人网| 青草久久国产| 亚洲专区字幕在线| 亚洲,欧美精品.| 非洲黑人性xxxx精品又粗又长| www日本在线高清视频| 亚洲精品在线观看二区| 国产又爽黄色视频| 如日韩欧美国产精品一区二区三区| 午夜福利影视在线免费观看| 亚洲精品美女久久av网站| 亚洲第一电影网av| 国产一区在线观看成人免费| 在线观看免费日韩欧美大片| 人妻丰满熟妇av一区二区三区| 久久精品国产清高在天天线| 亚洲成国产人片在线观看| 亚洲av成人不卡在线观看播放网| 午夜福利在线观看吧| 香蕉久久夜色| 欧美激情 高清一区二区三区| 99久久国产精品久久久| 午夜激情av网站| 亚洲av电影在线进入| 久久久久久久久免费视频了| 久久人人精品亚洲av| 欧美午夜高清在线| 亚洲av成人av| 麻豆国产av国片精品| 亚洲自偷自拍图片 自拍| av在线播放免费不卡| 1024视频免费在线观看| 国产精品久久电影中文字幕| 午夜成年电影在线免费观看| 国产亚洲av高清不卡| 一夜夜www| 我的亚洲天堂| 日韩免费av在线播放| 亚洲一码二码三码区别大吗| 91老司机精品| 亚洲国产高清在线一区二区三 | 亚洲片人在线观看| 最近最新中文字幕大全电影3 | 俄罗斯特黄特色一大片| 国产精品国产高清国产av| cao死你这个sao货| 国产精品亚洲一级av第二区| 午夜精品久久久久久毛片777| 97碰自拍视频| 亚洲欧美激情综合另类| avwww免费| 成人18禁高潮啪啪吃奶动态图| 中文字幕人妻丝袜一区二区| 国产亚洲精品综合一区在线观看 | 久久久水蜜桃国产精品网| 看片在线看免费视频| 丝袜美腿诱惑在线| 手机成人av网站| 丝袜人妻中文字幕| 日本免费一区二区三区高清不卡 | 国产蜜桃级精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 久久久水蜜桃国产精品网| 国产熟女xx| 欧美成人性av电影在线观看| 亚洲成人免费电影在线观看| 国产成人av教育| 午夜福利18| 亚洲国产中文字幕在线视频| 亚洲,欧美精品.| 黄色毛片三级朝国网站| 欧美成人性av电影在线观看| 亚洲中文字幕日韩| 久久久精品国产亚洲av高清涩受| av在线播放免费不卡| 国产日韩一区二区三区精品不卡| 琪琪午夜伦伦电影理论片6080| 熟妇人妻久久中文字幕3abv| 视频在线观看一区二区三区| 悠悠久久av| 国产区一区二久久| 丰满的人妻完整版| 99精品在免费线老司机午夜| 男女之事视频高清在线观看| 欧美成人午夜精品| 中文字幕人妻熟女乱码| 成人三级黄色视频| 最近最新中文字幕大全电影3 | 搡老妇女老女人老熟妇| 少妇 在线观看| 欧美人与性动交α欧美精品济南到| 亚洲成人久久性| 性欧美人与动物交配| 久久精品国产清高在天天线| 色综合婷婷激情| 国产精品久久久久久人妻精品电影| 一二三四在线观看免费中文在| 99香蕉大伊视频| 午夜免费成人在线视频| 色尼玛亚洲综合影院| 一二三四社区在线视频社区8| 韩国av一区二区三区四区| 色播亚洲综合网| 老司机靠b影院| 亚洲欧美激情综合另类| 国产av一区二区精品久久| 亚洲精品av麻豆狂野| 免费看美女性在线毛片视频| 国产亚洲av高清不卡| 久久久精品国产亚洲av高清涩受| 在线观看免费视频网站a站| 99久久综合精品五月天人人| 涩涩av久久男人的天堂| 精品午夜福利视频在线观看一区| 国产av一区二区精品久久| 一级毛片精品| 99国产精品一区二区蜜桃av| 国产成人欧美| 国产精品秋霞免费鲁丝片| 成人国产一区最新在线观看| 日本一区二区免费在线视频| 日本欧美视频一区| 国产精品 欧美亚洲| 精品国产美女av久久久久小说| 伦理电影免费视频| videosex国产| 亚洲午夜精品一区,二区,三区| 精品欧美国产一区二区三| 亚洲少妇的诱惑av| 国产人伦9x9x在线观看| 国产三级黄色录像| 美女国产高潮福利片在线看| 午夜精品久久久久久毛片777| 丁香欧美五月| 精品第一国产精品| 成在线人永久免费视频| 午夜两性在线视频| 女人精品久久久久毛片| 婷婷六月久久综合丁香| 色综合婷婷激情| 国产精品久久视频播放| 久久久久亚洲av毛片大全| 国内久久婷婷六月综合欲色啪| 亚洲精品粉嫩美女一区| 成人欧美大片| 成人免费观看视频高清| 成人国语在线视频| 91九色精品人成在线观看| 亚洲精品粉嫩美女一区| 我的亚洲天堂| 99精品欧美一区二区三区四区| 我的亚洲天堂| 亚洲av成人av| 50天的宝宝边吃奶边哭怎么回事| 亚洲熟妇中文字幕五十中出| 国内精品久久久久精免费| 长腿黑丝高跟| 一级毛片精品| 免费高清在线观看日韩| 香蕉久久夜色| 午夜亚洲福利在线播放| 又黄又爽又免费观看的视频| 老鸭窝网址在线观看| 亚洲精品一卡2卡三卡4卡5卡| 一个人观看的视频www高清免费观看 | 一个人免费在线观看的高清视频| netflix在线观看网站| 国产精品精品国产色婷婷| 757午夜福利合集在线观看| 欧美色视频一区免费| 男人舔女人下体高潮全视频| 国产精品野战在线观看| 中亚洲国语对白在线视频| 国产精品影院久久| 最新美女视频免费是黄的| 少妇裸体淫交视频免费看高清 | svipshipincom国产片| 精品久久久久久久毛片微露脸| 国产亚洲av高清不卡| 性欧美人与动物交配| 免费高清视频大片| 老熟妇乱子伦视频在线观看| 亚洲人成伊人成综合网2020| 91麻豆精品激情在线观看国产| 91在线观看av| 91国产中文字幕| 国产真人三级小视频在线观看| 91精品国产国语对白视频| 久久亚洲真实| 男女午夜视频在线观看| 一级毛片女人18水好多| 亚洲五月天丁香| 男女午夜视频在线观看| 怎么达到女性高潮| 亚洲五月天丁香| 久久久国产欧美日韩av| 99精品在免费线老司机午夜| 久久天躁狠狠躁夜夜2o2o| av中文乱码字幕在线| 精品第一国产精品| 男女下面插进去视频免费观看| 男人舔女人的私密视频| 欧美大码av| 国产成人精品无人区|