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

    Gut microbiota contributes to the development of endometrial glands in gilts during the ovary-dependent period

    2021-12-17 11:51:20BaoyangXuWenxiaQinYiqinYanYimeiTangShuyiZhouJunchengHuangChunlinXieLibaoMaandXianghuaYan

    Baoyang Xu,Wenxia Qin,Yiqin Yan,Yimei Tang,Shuyi Zhou,Juncheng Huang,Chunlin Xie,Libao Ma*and Xianghua Yan*

    Abstract

    Background:The hyper-prolificacy Meishan gilts achieved a superior endometrial gland development(EGD)than white crossbred gilts during the ovary-independent period(before 60 d of age).Then,the EGD continues under the management of ovary-derived steroid hormones that regulated by gut microbiota(after 60d of age).However,whether Meishan gilts’superiority in EGD lasting to the ovary-dependent period(after 60 d of age)and the role of gut microbiota in this period both remain unclear.

    Methods:Meishan gilts and Landrace x Yorkshire(LxY)gilts were raised under the same housing and feeding conditions until sexual maturity and then we compared their EGD and gut microbiota.Meanwhile,we transplanted fecal microbiota from Meishan gilts to L×Y gilts to explore the role of gut microbiota in EGD.We sampled plasma every 3 weeks and collected the uterus,ovary,liver,and rectal feces after the sacrifice.We then determined the hormone concentrations and expressions of the EGD-related genes.We also profiled the gut microbiota using 16S rDNA sequencing and metabolites of plasma and liver tissue using untargeted metabolomics.Finally,the correlation analysis and significant test was conducted between FMT-shifted gut microbes and EGD-related indices.

    Results:Meishan gilts have larger endometrial gland area(P<0.001),longer uterine horn length(P<0.01)but lighter uterine horn weight(P<0.05),a distinctive gut microbiota compared with L×Y gilts.Fecal microbiota transplantation(FMT)increased endometrial gland area(P<0.01).FMT markedly shifted the metabolite profiles of both liver and plasma,and these differential metabolites enriched in steroid hormone biosynthesis pathway.FMT increased estradiol and insulin-like growth factor 1 but decreased progesterone dynamically.FMT also increased the expression of the EGD-related genes estrogen receptor 1 gene,epithelial cadherin,and forkhead box protein A2.There is a significant correlation between FMT-shifted gut microbes and EGD-related indices.

    Conclusion:Sexually matured Meishan gilts achieved a superior EGD than LxY gilts.Meanwhile,gut microbiota contribute to the EGD potentially via regulating of steroid hormones during the ovary-dependent period.

    Keywords:Endometrial gland,Fecal microbiota transplantation,Gut microbiota,Meishan gilts,Steroid hormones,Untargeted metabolomics,Uterus

    Background

    Mammalian uteri contain glands in the endometrium that synthesize or transport and secrete substances essential for survival and development of the conceptus(embryo/fetus and associated extraembryonic membranes),thus contributing to female prolificacy[1,2].Endometrial adenogenesis is primarily a postnatal event in sheep,pigs,and rodents[3].At birth,no endometrial glands are present in uterus and then within 7d,endometrial glands begin to differentiate,and they continue to differentiate during the next several months[3].Neonatal porcine uterine development occurs in an ovaryindependent manner before,and in an ovary-dependent manner after,day 60 of neonatal life[4,5].In ovarydependent development period,ovarian hormones estradiol and progesterone play crucial roles in endometrial glands development[5].Estradiol is a primary mitogen for uterine epithelium and its increase around estrous induces proliferation of both uterine luminal epithelium and glandular epithelium[5].Neonatal estradiol exposure alters uterine morphology and endometrial transcriptional activity in prepubertal gilts[6].Progesterone can inhibit estradiol actions in stimulation of uterine epithelial proliferation through epithelial progesterone receptors[5,7].This relationship between estradiol and progesterone may regulate the endometrial gland development and uterine function.

    Meishan pigs as a Chinese indigenous breed is characterized by its prolificacy,producing an average of three to four more piglets per litter than European breeds[8].However,the underlying mechanism by which Meishan pigs achieve superior litter size remains not fully understood.Christenson et al.found that Meishan gilts’uterus secreted more abundance of endometrial proteins,which increased in association with endometrial gland development before day 60 of neonatal life[9].Further histological analysis validated that Meishan gilts achieved a superior endometrial gland development than white crossbred gilts at 60d of age[9].However,the further comparative morphological study of their endometrial gland development after 60 d of age remains undetermined.A recent study showed that the Meishan sows of high-prolificacy both have increased fecal microbiota diversity and fecal steroid hormones estradiol and progesterone levels than low-prolificacy ones,which may contribute to the improvement of sows’reproductive performance[10].Interestingly,microbiota that harbor inside the host intestine function as a vital endocrine organ get involved in regulating steroid hormones estradiol and progesterone[11,12].

    Given that the vital role of steroid hormones estradiol and progesterone in endometrial glands development during the ovary-dependent period[7,13,14],we hypothesis that gut microbiota get involved in the superior endometrial gland’s development of prolific Meishan gilts than white crossbred gilts during the ovarydependent period.Fecal microbiota transplantation(FMT)involves the transfer of donor fecal microbiome to a recipient in order to establish a more desirable or donor-like gut microbiome[15,16].FMT provides an effective approach to link gut microbiota and host phenotype[15,17].To test our hypothesis,we compared the uterine development and gut microbiota between sexually matured Meishan gilts and crossbred Landrace× Yorkshire(L×Y)gilts under the same housing and feeding conditions.Meanwhile,for exploring the role of gut microbiota in endometrial glands development,we transplanted fecal microbiota of Meishan gilts to recipient L×Y gilts.

    Methods

    Experimental design and animals

    All experiments involving swine were carried out under the recommendations of the Guide for the Care and Use of Laboratory Animals Monitoring Committee of Hubei Province,China,and the protocol was approved by the Scientific Ethics Committee of Huazhong Agricultural University(approval numbers HZAUSW-2018-015).All efforts were made to minimize animal suffering.Eight Meishan gilts of age 95±7d from the National Conservation Farm of Meishan Pigs(Jiading,Shanghai)and 14 Landrace× Yorkshire(L×Y)gilts of age 90±4d from a commercial farm were housed in National Engineering and Technology Research Center for Livestock and fed with diets according to NRC(2012).Experimental design:L×Y gilts receive 20mL saline(Ctrl)or fecal suspension from Meishan gilts(FMT)and these Meishan gilts receive the same volume of saline(MS)every other day until the puberty onset.For reducing the anaerobic bacteria die,we prepared fecal slurry with O2?free saline immediately after collecting fresh feces.We then regulated the concentration of fecal slurry after counting the live microbes using optical microscopy combined with methylene blue staining[18].All gilts were sacrificed at 24h after detection of the third estrus(excluded a Meishan gilts because of its health concerns).Blood samples(10mL)were collected and then centrifugate for 10min at 3,500×g at 4 °C to get plasma samples.Immediately after slaughter,intestinal contents from rectum were collected and placed in liquid nitrogen.We determined the weight and length of uterine horn.Uterus tissue was removed from the middle portion of the uterus and rinsed with iced cold PBS,then fixed in 4% paraformaldehyde solution for measuring uterus morphology or placed in liquid nitrogen.

    Fecal microbiota transplantation(FMT)experiment

    The fecal suspension was prepared using the protocol previously described with minor optimization [19].Briefly,fresh feces samples were obtained from Meishan gilts and immediately homogenized in sterile and O2-free saline solution.Then fecal slurry was passed through the sterilized gauze and then a 0.224-mm stainless cell strainer to remove the particles.We used optical microscopy combined with methylene blue staining to count the live microbes in the fecal slurry.The Meishan gilts’neighbor L×Y gilts received sterile saline(20mL)or fecal microbiota suspensions(107CFU/mL,20mL)every other day from the age of 90 d to the age of puberty onset(about 200 d of age).

    Evaluation of endometrial gland development

    We evaluated endometrial gland development by examining hematoxylin-and eosin-stained(H&E)uterine wall sections(10μm)using computer-assisted morphometry.For each section,six areas of the endometrium from lumen to myometrium were randomly selected and measured.Then,the total area of glands within the endometrial area was measured by encircling each gland within the previously outlined endometrial area.Gland development for that section was then expressed as the ratio of gland area to endometrial area[20].

    Gut microbiota profiling

    The total genomic DNA of fecal bacteria(stool from the rectum)was extracted using the protocol of the Repeated Bead Beating Plus Column Method[21].The integrity of DNA was assessed by agarose gel electrophoresis.The genomic DNA was used as a template for PCR amplification.Universal primers 338F and 806R were used for PCR amplification of the V3–V4 hypervariable regions of 16S rRNA genes(338F,5′-ACTC CTACGGGAGGCAGCA-3′,806R,5′-GGACTACHVG GGTWTCTAAT-3′)[22].The generated DNA pool was then sequenced on the Illumina HiSeq system with the sequencing strategy PE300 and then the sequencing data were analyzed using the Quantitative Insights Into Microbial Ecology software package[23].

    To obtain more accurate and reliable results in subsequent bioinformatics analysis,raw data were cleansed by the in-house procedure[22].Then paired end reads with overlap were merged to tags.The high-quality pairedend reads were combined to tags based on overlaps with FLASH[24].The tags were clustered to OTU(Operational Taxonomic Unit) by scripts of software USEARCH(v7.0.1090,[25]).After that,they were clustered into OTU with a 97% threshold using UPARSE,and the OTU unique representative sequences were obtained;Chimeras were filtered out by UCHIME(v4.2.40);OTU representative sequences were taxonomically classified using Ribosomal Database Project(RDP)Classifier v.2.2 trained on the database Greengene_2013_5_99[26]with 0.6 confidence values as cutoff.The alpha diversity indices including observed species value,Chao1 value,ACE value,Shannon,and Simpson value are calculated by Mothur(v1.31.2,[27])with the corresponding rarefaction curve was drawn by software R(v3.6.0).Phylogenetic beta diversity measures such as unweighted UniFrac distance metrics analysis and principalcomponent analysis(PCoA)were done using the QIIME 2[23].Functional analysis of gut microbiota was predicted by Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2(PICRUSt2,[28]).

    Gut microbial metabolites

    We determined the microbial Metabolites short-chain fatty acids(SCFAs)(acetate,propionate,butyrate)in the diegesta of colon using gas chromatography with a modification of the previous method[29].In brief,1g of the digesta samples was weighed into a 2-mL centrifuge tube with 1 mL of methanol added.After being vortexed for 30 s,the sample was centrifuged(12,000× g)at 4°C for 10min.The supernatant was transferred into centrifuge tubes(2 mL)and mixed with 0.2mL 25% metaphosphoric acid.After 30min at 4°C,the tubes were centrifuged(12,000× g)again at 4°C for 10 min.To quantify SCFAs,a calibration curve for the concentration range of 0.015–1 mg/mL was constructed.SCFAs measurements were performed following a recently published protocol[30].

    Quantitative enzyme-linked immunosorbent assay(ELISA)

    Steroid hormones estradiol,progesterone,and insulinlike growth factor 1(IGF-1)in plasma and tissue were quantified by ELISA(Cusabio,Wuhan,China).We collected plasma samples using EDTA as an anticoagulant and then the samples were centrifuged for 15min at 1,500× g and 4°C.For tissue homogenates,100 mg tissue was rinsed with 1× PBS,homogenized in 1mL of 1×PBS and stored overnight at ?20°C.After two freezethaw cycles were performed to break the cell membranes,the homogenates were centrifuged for 5 min at 5,000× g,and 4°C.The supernatant was removed and assayed immediately.Alternatively,aliquot and store samples at ?20°C or?80 °C.Centrifuge the sample again after thawing before the assay.Avoid repeated freeze-thaw cycles.All samples were measured in 3 replicates following the recommended procedures in the instruction.

    Quantitative real-time PCR

    qPCR was performed with the PowerUp? SYBR? Green Master Mix(Applied Biosystems)on a CFX384 Real-Time PCR system(Bio-Rad)and changes in gene expression were calculated relative to GAPDH.RT-qPCR primer pairs are:estrogen receptor 1 gene (ESR1)5′-AGCACCCTGAAGTCTCTGGA-3′, 5′-TGTGCC TGAAGTGAGACAGG-3′;epithelial cadherin(CDH1)5′-ATGTGCACGTATGCGACTGT-3′,5′-GGAACT TGCAATCCTGCTTC-3;forkhead box protein A2(FOXA2)5′-ATGCTGGGAGCGGTGAAGAT-3′,5′-AGCGAGTGGCGGATGGAGTT-3′;glyceraldehyde-3-phosphate dehydrogenase(GAPDH)5′-T CGGAGTGA ACGGATTTGGC-3′,5′-TGCCGTGGGTGGAATCAT AC-3′.

    Hematoxylin and eosin(H&E)staining

    Briefly,ovaries were preserved in 4% paraformaldehyde overnight at 4°C(Sigma-Aldrich).The tissues were subsequently embedded in paraffin wax(Fisher Scientific)following immersion in a graded series of alcohols(70–100%).Embedded tissue was sectioned(5 mm)using a rotary microtome.We adhered these sections to microscope slides and then dried them at 56°C for 24 h.Next,we passed slides through a series of the clearing agent xylene and rehydrated in a graded series of ethanol(100%,100%,100%,95%,80%).After a brief wash in distilled water,we incubated the slides with hematoxylin solution.The sections were then washed with running tap water to remove excess hematoxylin.Then,we differentiated the sections in 1% acid alcohol for 30 s and then washed them with running tap water for 1 min.This step was followed by an incubation in the eosin counterstain,subsequent dehydration in a graded series of ethanol(80%,95%,95%,95%,100%,100%,100%),and immersion in xylene.

    Untargeted metabolomics profiling

    Sample preparation.Plasma metabolite extraction was performed by mixing 40 μL thawed plasma and 120 μL ice-cold methanol,vortexing for 1 min,10min standing at room temperature,and keeping at ?20°C for 2 h.After centrifuged at 4,000 × g and 4 °C for 30min,25 μL of the supernatant was transferred into a glass vial containing 225 μL of 50% methanol.Take 50 μL of each sample and mix it into a quality control(QC)sample and divide all samples into 96-well plates in the order of loading and put 60μL in each well.Liver metabolite extraction was performed by mixing 25 mg thawed liver tissue and 800μL of chilled methanol/water(1:1)solution.Add two small steel balls to each EP tube,place the sample in TissueLyser and grind.The parameter is set to 50 Hz for 3 min.After grinding,remove the steel ball,place the EP tube in the ?20°C refrigerator for 2 h.After centrifuged at 30,000 × g and 4°C for 15min,carefully take 650 mL of each sample into a new EP tube.Then centrifuged at 25,000× g and 4°C for 20min,take 550 mL of each sample in a new EP tube.Take 1 mL of acetonitrile through the column with a positive pressure extractor,then take 1mL of 100% methanol through the column at a constant rate,and finally take 1 mL of Milli-Q pure water through the column at a constant rate.Draw 400 mL supernatant after centrifugation for each sample,transfer to the activated and equilibrated SPE column,after repeating once,discard the supernatant,add 400mL 100% acetonitrile to elute,keep eluted solution in a new EP tube.Take 40μL of all samples and mix them into QC.Dispense all the samples into 96-well plates with 60 μl of each sample in each well.

    Metabolomics analysis.For plasma and liver tissue sample,all samples were gained by the Liquid Chromatography Mass Spectrometry(LC-MS)system followed machine orders.First,all chromatographic separations were performed using an ultra-performance liquid chromatography(UPLC)system(Waters,UK).For the reversed phase separation,plasma sample used an ACQUITY UPLC BEH C18 column(100 mm×2.1mm,1.7μm,Waters,UK)and liver tissue sample used an ACQUITY UPLC HSS T3 column(100 mm×2.1mm,1.8μm,Waters,UK).The column oven was maintained at 50°C.The flow rate was 0.4 mL/min and the mobile phase comprised solvent A(water+0.1% formic acid)and solvent B(acetonitrile+0.1% formic acid).Gradient elution conditions were set:0–2 min,100% phase A;2–11min,0 to 100%B;11–13min,100%B;13–15 min,0 to 100%A.The injection volume for each sample was 10μL.A high-resolution tandem mass spectrometer Xevo G2 XS QTOF(Waters,UK)was used to detect metabolites eluted form the column.The Q-TOF was operated in both positive and negative ion modes.For positive ion mode,the capillary and sampling cone voltages were set at 3.0kV and 40.0 V,respectively.For negative ion mode,the capillary and sampling cone voltages were set at 2.0kV and 40.0 V,respectively.The mass spectrometry data were gained in Centroid MSE mode.The TOF mass range was from 50 to 1200Da and the scan time was 0.2s.For the MS/MS detection,all precursors were fragmented using 20–40 eV,and the scan time was 0.2s.During the acquisition,the LE signal was gained every 3 s to calibrate the mass accuracy.In order to test the stability of the LC-MS during the whole acquisition,a quality control sample(Pool of all samples)was gained after every 10 samples.Peak extraction is mainly achieved through the commercial software Progenesis QI(version 2.2),including peak alignment,peak extraction,normalization,deconvolution,and compound identification.Based on QC sample information,local polynomial regression fitting signal correction(Quality control–based robust LOESS signal correction,QCRSC)is performed on the real sample signal[31].

    Statistical analysis

    Results were presented as mean±SEM.Experimental data were analyzed by one-way analysis of variance tests,followed by Fisher’s least significant difference and the Duncan multiple comparison test with GraphPad 8.0 software.Significance was presented as*P<0.05,and**P<0.01,whereas P values between 0.05 and 0.10 were considered as indicative of a trend.

    Results

    Meishan gilts have distinct uterine characters and gut microbiota compared with LxY gilts

    Under the same housing and feeding conditions,Meishan gilts have larger(P<0.001)endometrial gland area normaled to the endometrial area(Fig.1a-b)and longer(P<0.01 for left uterine horn and P<0.001 for right uterine horn)uterine horn length(Fig.1c)but lighter(P<0.05)uterine horn weight(Fig.1d)compared with L×Y gilts at the age of third estrus.The beta diversity of gut microbiota showed by PCoA showed obvious difference between Meishan gilts and L×Y gilts(Fig.1e).Meishan gilts had superior alpha diversity indices including observed species(P<0.01),Chao1(P<0.01),ACE(P<0.01),Shannon(P<0.01)and Simpson(P<0.001)compared with L×Y gilts(Fig.1f).The alpha diversity index Good coverage suggested that all the sequencing data were of high quality and could cover the most species.Further taxon analysis showed that Meishan gilts’gut harbored more bacteria of phyla Actinobacteria(P<0.01),Firmicutes(P<0.01),Fibrobacteres(P<0.001),Lentisphaerae(P<0.05),and Synergistetes(P<0.01)but fewer bacteria of phyla Bacteroidetes(P<0.001),and WPS-2(P<0.01)compared with L×Y gilts(Supplemental Figure S1a).At genus level,Meishan gilts’gut harbored more Bacteroides(P<0.01),Bifidobacterium(P<0.001),Fibrobacter(P<0.001),Oscillospira(P<0.05),SMB53(P<0.05),and Turicibacter(P<0.05)but less Anaerovibrio(P<0.05),Bulleidia(P<0.01),Butyricicoccus(P<0.01),Coprococcus(P<0.05),Dorea(P<0.05),Gemmiger(P<0.05),Lachnospira(P<0.01),Pseudobutyrivibrio(P<0.01),and Streptococcus(P<0.05)compared with L×Y gilts(Supplemental Figure S1b).

    Fig.1 Meishan gilts have distinct uterine characters and gut microbiota compared with LxY gilts.a Representative H&E-stained uterine sections.Scale bar,1,000μm and 200μm.b Endometrial gland area normalized to endometrial area.Length(c)and weigh(d)of Uterine horn.e PCoA based on weighted UniFrac distance show a distinct fecal microbial composition inter-groups.Each data point represents a sample from a distinct sow projected onto the first two principal coordinates(percent variation explained by each PCo is shown in parentheses).f Species diversity analysis of gut microbiota with alpha diversity indices including Observed species,Chao1,ACE,Shannon,Simpson and Good coverage.*P<0.05,**P<0.01,***P<0.001,data were shown as means±SEM;n=7

    FMT shifted gut microbiota and increased the endometrial gland area in recipient LxY gilts

    To explore the association between the gut microbiota and uterine development,we transplanted the fecal microbiota of Meishan gilts to L×Y gilts.PCoA showed that FMT distinctly shifted the gut microbiota composition of the recipient L×Y gilts and made it closed to that of Meishan gilts(Fig.2a).The top two principalcomponent(PC1 and PC2)suggested that FMT explained 60.49% of the variation(Fig.2a).FMT recipient gilts had superior alpha diversity indices including observed species(P<0.01),Chao1(P<0.01),ACE(P<0.05),Shannon(P<0.05)and Simpson(P<0.01)compared with control L×Y gilts(Fig.2b).The alpha diversity index Good coverage suggested that all the sequencing data were of high quality and could cover the most species.Further taxon analysis showed that FMT recipient gilts’gut harbored more bacteria of phyla Firmicutes(P<0.001),Fibrobacteres(P<0.001),but fewer bacteria of phyla Bacteroidetes(P<0.001)compared with control L×Y gilts(Supplemental Figure S2a).At genus level,FMT recipient gilts’gut harbored more Bifidobacterium(P<0.01)and Fibrobacter(P<0.01)compared with control L×Y gilts(Supplemental Figure S2b).Morphological analysis(H&E)of uterine showed that FMT recipient L×Y gilts had increased(P<0.01)endometrial gland area(Fig.3c-d).FMT had no effect on the uterine horn length(Fig.3e)but had a trend(P=0.069 for left uterine and P=0.067 for right uterine)to reduce the uterine horn weight(Fig.3f).

    Fig.2 FMT shifted gut microbiota and increased the endometrial gland area in recipient L×Y gilts.a PCoA based on weighted UniFrac distance show a distinct fecal microbial composition intergroups.Each data point represents a sample from a distinct sow projected onto the first two principal coordinates(percent variation explained by each PCo is shown in parentheses).b Species diversity analysis of gut microbiota with alpha diversity indices including Observed species,Chao1,ACE,Shannon,Simpson and Good coverage.c Representative H&E-stained uterine sections.Scale bar,1,000μm and 200μm.d Endometrial gland area normalized to endometrial area.Length(e)and weigh(f)of Uterine horn.*P<0.05,**P<0.01,***P<0.001,data were shown as means±SEM;n=7

    Untargeted metabolomics revealed metabolism shift induced by FMT in liver and plasma

    To link the shifted gut microbiota and increased endometrial gland area,we determined the metabolites profiles of liver and plasma using untargeted metabolomics.Partial least squares discriminant analysis(PLS-DA)showed that the detected metabolites of both pos and neg mode were obviously shifted by FMT in liver(Fig.3a)and plasma(Fig.3b).The volcano plots showed the significantly differential metabolites(red dot)under both pos and neg mode in liver(Fig.3c)and plasma(Fig.3d).We further conducted KEGG pathway enrichment analysis for the different metabolites in both liver and plasma.For liver,the FMT-induced differential metabolites enriched in pathways of steroid hormone biosynthesis,steroid biosynthesis,ovarian steroidogenesis,linoleic acid metabolism,folate biosynthesis,and biosynthesis of unsaturated fatty acids(Fig.3e).And for plasma, the FMT-induced differential metabolites enriched in pathways of steroid hormone biosynthesis,primary bile acid biosynthesis,fatty acid metabolism,fatty acid degradation,estrogen signaling pathway,butanoate metabolism,and alpha-Linolenic acid metabolism(Fig.3f).Noteworthily,steroid hormone biosynthesis pathway was enriched by FMT-induced differential metabolites from both liver and plasma,which indicated that the circulating steroid hormone profile has been shifted by FMT.

    Fig.3 FMT shifted the metabolite profiles of liver and plasma in recipient L×Y gilts.PLS-DA of the detected metabolites of both pos and neg mode in liver(a)and plasma(b).The volcano plots of the significantly differential metabolites(red dot)under both pos and neg mode in liver(c)and plasma(d).KEGG pathway enrichment analysis for the different metabolites in both liver(e)and plasma(f).n=7

    FMT shifted endometrial gland development-related hormones in recipient LxY gilts

    Given that FMT shifted the circulating steroid hormone profile,we then determined the concentration of steroid hormones estradiol and progesterone,and IGF-1 in plasma.FMT increased plasma estradiol concentration on 132(P<0.05),153(P <0.05),174(P <0.05),and 195(P<0.05)d of age(Fig.4a),but decreased plasma progesterone concentration on 111(P<0.05),132(P=0.057),and 153(P<0.05)d of age(Fig.4b).Estradiol and progesterone can regulate of IGF-1 gene expression in pig uterus[32].We therefore determined the IGF-1 concentration in plasma and uterine tissue.The results showed that FMT increased IGF-1 on 111(P<0.05),132(P<0.05),153(P=0.07),174(P<0.05),and 195(P<0.01)d of age(Fig.4c),and in uterine tissue(P<0.05)(Fig.4f).The enhanced secretion of estradiol(P<0.01)and progesterone(P<0.001)were determined in ovary tissue(Fig.4d-e).Correspondingly,FMT elevated(P<0.01)the mRNA expression level of estrogen receptor 1 gene(Esr1)(Fig.4g).Besides,FMT elevated the mRNA expression level of uterine development-related epithelial cadherin(CDH1)(P<0.01)and forkhead box protein A2(FOXA2)(P<0.05)(Fig.4g).

    Fig.4 FMT shifted uterine gland development-related hormones and gene expression in recipient LxY gilts.Dynamic shifts of uterine gland development-related hormones estradiol(a),progesterone(b),and IGF-1(c).The concentrations of estradiol(d)and progesterone(e)in ovary tissue;and IGF-1 concentration(f)in uterine tissue.The relative mRNA expression of endometrial gland development-related genes(g).*P<0.05,**P<0.01,***P<0.001,data were shown as means±SEM;n=7

    FMT promoted carbohydrate metabolism and digestive system functions of gut microbiota in recipient LxY gilts

    To explore the underlying mechanism by which the shifted gut microbiota contributes to the endometrial gland development,we determined the function shift of gut microbiota using PICRUSt2.The results suggested that FMT enriched KEGG functions of carbohydrate metabolism(P<0.01)and digestive system(P<0.05)in gut microbiota of recipient L×Y gilts(Fig.5a).This enrichment effect of FMT was confirmed by the increased concentrations of SCFAs propionate(P<0.01)and butyrate(P<0.05)but not acetate(Fig.5b-d).

    Fig.5 FMT promoted carbohydrate metabolism and digestive system functions of gut microbiota in recipient L×Y gilts.KEGG functions of gut microbiota(a)and the corresponding concentrations of SCFAs(b).The differential KEGG function was indicated in red.*P<0.05,**P<0.01,data were shown as means±SEM;n=7

    FMT-shifted gut microbes were correlated to the indices of uterine development

    We finally determined the correlation between FMT-shifted gut microbes and uterine development-related indices using spearman’s rank correlation coefficient and significance test.The results in Fig.6 showed that Lentisphaerae(P<0.05),Bifidobacterium(P<0.05),and Fibrobacter(P <0.05) were positively correlated to endometrial gland area.Firmicutes(P <0.01) and Fibrobacter(P<0.05),were positively correlated to estradiol concentration;Bacteroidetes was significantly negatively correlated to estradiol(P <0.01).For IGF-1,Bacteroidetes was negatively correlated to IGF-1 concentrations(P<0.01);Fibrobacteres(P<0.05),Firmicutes(P<0.001),Bifidobacterium(P<0.05),and Fibrobacter(P<0.05)were positively correlated to IGF-1 concentration.For progesterone,Bacteroidetes was positively correlated to progesterone concentration(P<0.01);Fibrobacteres(P<0.001),Firmicutes(P<0.01),Bifidobacterium(P<0.01),and Fibrobacter(P<0.01)were negatively correlated to progesterone concentration.

    Fig.6 FMT-shifted gut microbes were correlated to the indices of uterine development.The correlation between FMT-shifted gut microbes and uterine development-related indices using spearman’s rank correlation coefficient and significance test(n=7).*P<0.05,**P<0.01,***P<0.001;n=7

    Discussion

    Porcine uterine development plays a crucial role in determining sows’prolificacy[4,5].Earlier evidences suggested that porcine uterine development occurs in an ovary-independent manner before,and in an ovarydependent manner after,60 d of age[4,5].Before 60 d of age,Christenson et al.found that Meishan gilts achieved a superior endometrial gland development than white crossbred gilts[9].Here,we further compared the uterine development between Meishan gilts and L×Y gilts during ovary-dependent development period(after 60d of age),and found that Meishan gilts have increased endometrial gland area and uterine horn length than L×Y gilts.Thus,the superior uterine development of Meishan gilts than white crossbred L×Y gilts lasted from neonate to mature.In ovary-dependent development period,ovary-derived hormones estradiol and progesterone play crucial roles in endometrial glands development[5].Recent evidences suggest that gut microbiota plays a role in the regulation of ovarian hormones estradiol and progesterone[12,33,34].We found that Meishan gilts had a distinctive structure of gut microbiota and superior alpha diversity including observed species,Chao1,ACE,Shannon and Simpson compared with L×Y gilts under the same housing and feeding conditions.These conservative differences of gut microbiota between Meishan gilts and L×Y gilts indicated that host genetics can shape gut microbiota,which had been also reported in human[35].Compared well with our results,Yuan et al.reported patients with endometriosis,a chronic and estrogen-dependent disease characterized by endometrial glands and stroma outside the uterus,had a disordered gut microbiota of reduced alpha diversity[36].These increased alpha diversity indices had been also determined in Meishan sows whose reproductive performance were promoted by dietary crude fiber supplementation[10].Further taxon analysis showed that Meishan gilts’gut harbored more bacteria of phyla Actinobacteria Firmicutes,Fibrobacteres,Lentisphaerae,and Synergistetes but fewer bacteria of phyla Bacteroidetes and WPS-2 compared with L×Y gilts.Yang et al.also reported a similar difference of gut microbiota between Meishan sows and Duroc sows except Bacteroidetes[37].These inconsistent results may because of the detect method(real-time PCR analysis vs.16S rDNA high-throughput sequencing)or animal conditions(age etc.).Consistent with our results,Jiang et al.also determined an increased Bacteroidetes in fertility-improved Meishan sows by dietary crude fiber supplementation[10].Among the changed bacteria at genus levels,the increased Fibrobacter,but decreased Streptococcus in Meishan gilts were comparable to the earlier reports[10,38].Taken together,our data firstly suggested that sexually mature Meishan gilts achieved superior endometrial gland area and uterine horn length than L×Y gilts,which may be regulated at least partly by their distinctive gut microbiota.

    We then explored the association between gut microbiota and uterine development by transplanting the fecal microbiota of Meishan gilts to L×Y gilts.As reported in other FMT experiment[39,40],FMT had a lasted shifteffect on the gut microbiota of recipient L×Y gilts and made them similar to that of the donor Meishan gilts.Interestingly,FMT from Meishan gilts increased endometrial gland area in recipient L×Y gilts.To further understand the association between the shifted gut microbiota and increased endometrial gland area,we profiled the metabolites of liver and plasma and found that FMT markedly shifted the metabolite profiles of both liver and plasma,and these FMT-induced differential metabolites from both liver and plasma enriched in steroid hormone biosynthesis pathway.Compared well with the results of metabolomics,we found that FMT increased estradiol concentrations,but decreased progesterone concentration in plasma during the trial dynamically.Ovary-derived steroid hormones of estradiol and progesterone play crucial roles in endometrial glands development[7].Estradiol is a primary mitogen for uterine epithelium and can induce proliferation of both uterine luminal epithelium and glandular epithelium[7].Progesterone reversely inhibit estradiol actions in stimulation of uterine epithelial proliferation through epithelial progesterone receptors[7,13].Therefore,FMT induced an increase of plasma estradiol concentrations,but a decrease of plasma progesterone concentration during our trial dynamically may account for its effect on promoting endometrial gland area.The elevated estradiol concentration but decreased progestogen in ovary tissue showed that FMT could regulate steroid hormones synthesis of ovary.Besides,estradiol and progestogen affect endometrial IGF-I expression and secretion [32].In our study, FMT elevated IGF-1 concentration,and increased expression of ESR1 gene in uterine tissue,suggesting that gut microbiota affect endometrial IGF-I expression and secretion by regulating steroid hormones estradiol and progestogen.CDH1 and FOXA2 are essential for endometrial gland differentiation and development[41,42].The elevated expression of both CDH1 and FOXA2 genes further confirmed that FMT promoted endometrial gland differentiation and development in recipient L×Y gilts.Steroid hormones regulate gene expression of CDH1 during the peri-implantation period of pregnancy in pigs[43].In ovariectomized mice,estradiol suppressed Foxa2 mRNA expression at 48h but not 4 h and 24h after treatment[44].This suppression effect of estradiol was sped up by co-treatment with progestogen[44].However,the role of steroid hormones on the expression of CDH1 and FOXA2 genes during pre-puberty gilts remain unknown.Overall,FMT shifted the steroid hormones profiles,which may contribute to its effect on promoting endometrial gland development.

    Given that gut microbiota gets involved in regulation of ovarian hormones estradiol and progesterone[12,33,34],we then identified the FMT-shifted bacteria.FMT increased Firmicutes,Fibrobacteres,Bifidobacterium, and Fibrobacter, but deceased Bacteroidetes.Nuriel-Ohayon et al.found that Bifidobacterium abundance increased in the gut during pregnancy in women and mice,and progesterone supplementation increases Bifidobacterium abundance in mice and in vitro[12].Consistently,these FMT-shifted bacteria along with increased estradiol and progesterone had been also determined in Meishan sows whose reproductive performance had been promoted by dietary crude fiber supplementation[10],which suggested the potential interaction between gut microbiota and progesterone.More directly,gut microbiota gets involved in enterohepatic circulation of estradiol by encoding β-glucuronidase,which plays an important role in cleaving conjugated estradiol[38,45].The gut microbiome includes an ‘estrobolome’that the aggregate of enteric bacterial genes whose products can cleave conjugated estradiol[45].Notably,these FMT-shifted bacteria belong to the estrobolome.Taken together,FMT-shifted gut microbiota may promoted endometrial gland area by regulating steroid hormones.Furthermore,We found that FMT enriched KEGG functions of carbohydrate metabolism and digestive system.These function shifts of gut microbiota compared well with the increased SCFAsproducing bacteria,including Fibrobacter.In addition,these enrichment effects of FMT was confirmed by the increased concentrations of fecal SCFAs,including propionate and butyrate but not acetate.Although there is no direct evidence of SCFAs regulating endometrial gland development,the elevated fecal SCFAs get well along with the hyper-prolificacy of sows had been reported in Meishan sows and crossed white sows[10,38,46].Recently studies indicated that these elevated SCFAs may promote sows’prolificacy indirectly.Administration of SCFAs to mice alleviate stress-induced brain–gut axis alterations[47].Kimura-Todani et al.also reported that acetate exerted an anxiolytic effect on the host[48].However,we haven’t observed a significant increase of acetate concentration in distal colon,which may because of the acetate could be absorbed in colon[49].In modern pig production system,sows suffer from postpartum stress even anxiety,which results in decreased maternal performance of lactating sows and increased preweaning piglet mortality[50,51].Given that the number of pigs weaned per sow per year is used to determine sows’prolificacy[52],gut microbial metabolites SCFAs may contribute to the hyper-prolificacy of sows by ameliorating sows’postpartum stress.Besides,the significant correlation between FMT-shifted gut microbes and uterine development-related indices further supports our hypotheses that gut microbiota contribute to endometrial gland development during the ovary-dependent period.

    Conclusion

    In conclusion,this study presented that sexually mature Meishan gilts achieved a superior uterine development than L×Y gilts under the same housing and feeding conditions.Meanwhile,Meishan gilts harbored a distinctive gut microbiota that superior in alpha diversity when compared with L×Y gilts.Transplanting fecal microbiota from Meishan gilts made the gut microbiota of recipient L×Y gilts similar to that of the donors lastingly and promote endometrial gland development,simultaneously.The FMT-induced profiler shifts of steroid hormones which potentially driven by gut microbiota contribute to the promotion of endometrial gland development.Therefore,our data suggests that gut microbiota contribute to the endometrial gland development during the ovary-dependent period.

    Abbreviations

    SCFAs:Short chain fatty acids;ELISA:Enzyme-linked immunosorbent assay;IGF-1:Insulin-like growth factor 1;EDTA:Ethylenediaminetetraacetic acid;ESR1:Estrogen receptor 1 gene;CDH1:Cadherin-1;FOXA2:Forkhead box protein A2;H&E:Hematoxylin and eosin;LC-MS:Liquid chromatography mass spectrometry;PICRUSt2:Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2;OTU:Operational taxonomic units

    Supplementary Information

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

    Additional file 1.Supplemental Table S1.Ingredient composition and chemical analysis of the basal diet for tilapia.

    Additional file 2:Supplemental Figure S1.Relative abundance of gut microbes in Meishan and LxY gilts.

    Additional file 3:Supplemental Figure S2.FMT shifted the composition of gut microbiota in recipient LxY gilts.

    Acknowledgements

    We thank Yan lab’s members for helpful discussions and critical reading of the manuscript.

    Authors’contributions

    B.X.,L.M.,and X.Y.designed the research.B.X.,Y.Y.,W.Q.,Y.T.,S.Z.,and J.H.performed the research.B.X.,Y.Y.,W.Q.,Y.T.,S.Z.,and X.Y.analyzed the data.B.X.,and X.Y.wrote the paper with the help of all authors.All authors read and approved the final version of the manuscript.

    Funding

    This work was supported by the Natural Science Foundation of China(31730090 and 31925037),Hubei Provincial Natural Science Foundation of China(2018CFA020).

    Availability of data and materials

    All data generated or analyzed during this study can be made available by the corresponding author upon reasonable request.

    Declarations

    Ethics approval

    The animal model and experimental procedures used in this experiment were approved by the Scientific Ethics Committee of Huazhong Agricultural University(approval numbers HZAUSW-2018-015).

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Author details

    1State Key Laboratory of Agricultural Microbiology,College of Animal Sciences and Technology,Huazhong Agricultural University,Wuhan 430070,Hubei,China.2The Cooperative Innovation Center for Sustainable Pig Production,Wuhan 430070,Hubei,China.3Hubei Provincial Engineering Laboratory for Pig Precision Feeding and Feed Safety Technology,Wuhan 430070,Hubei,China.

    Received:29 October 2020 Accepted:1 March 2021

    久久亚洲精品不卡| 国产毛片a区久久久久| 毛片一级片免费看久久久久| 最好的美女福利视频网| 青春草国产在线视频 | 久久久久久伊人网av| 国产成人a区在线观看| 一级二级三级毛片免费看| 国产亚洲5aaaaa淫片| 国产亚洲精品av在线| 精品无人区乱码1区二区| 超碰av人人做人人爽久久| 亚洲精品成人久久久久久| 伦精品一区二区三区| 在线免费观看不下载黄p国产| 久久精品国产亚洲av香蕉五月| 一级黄色大片毛片| 久久精品国产鲁丝片午夜精品| 简卡轻食公司| 婷婷六月久久综合丁香| 亚洲色图av天堂| 国内精品一区二区在线观看| 美女黄网站色视频| 国产午夜精品久久久久久一区二区三区| eeuss影院久久| 九色成人免费人妻av| 日日干狠狠操夜夜爽| 久99久视频精品免费| 成年免费大片在线观看| 欧美一级a爱片免费观看看| 欧美不卡视频在线免费观看| 国产精品1区2区在线观看.| 久久久久久久久久成人| 尾随美女入室| 亚洲国产精品久久男人天堂| 国产成人a∨麻豆精品| 国产私拍福利视频在线观看| 婷婷六月久久综合丁香| 看非洲黑人一级黄片| 国产亚洲91精品色在线| 一夜夜www| 午夜激情欧美在线| 老女人水多毛片| 欧美高清成人免费视频www| 丝袜喷水一区| 看片在线看免费视频| 国产亚洲精品久久久久久毛片| 亚洲精品自拍成人| 亚洲电影在线观看av| 亚洲成人av在线免费| 可以在线观看毛片的网站| 男人的好看免费观看在线视频| 国产精品一区二区三区四区免费观看| 国产真实乱freesex| 国产一区二区激情短视频| 人人妻人人澡人人爽人人夜夜 | 国产精品一二三区在线看| 天天躁日日操中文字幕| 91久久精品国产一区二区三区| 一区福利在线观看| 亚洲精品乱码久久久久久按摩| 国产成人精品久久久久久| 亚洲无线观看免费| 内射极品少妇av片p| 日韩欧美在线乱码| 插逼视频在线观看| 人妻制服诱惑在线中文字幕| 国产精品1区2区在线观看.| 精品国内亚洲2022精品成人| 亚洲成人久久爱视频| 中文字幕熟女人妻在线| 人人妻人人澡欧美一区二区| 亚洲国产高清在线一区二区三| av.在线天堂| 热99re8久久精品国产| 麻豆av噜噜一区二区三区| 久久久国产成人精品二区| 啦啦啦观看免费观看视频高清| 亚洲性久久影院| 国产高潮美女av| 亚洲美女视频黄频| 18禁黄网站禁片免费观看直播| 国产一区亚洲一区在线观看| 日本爱情动作片www.在线观看| 欧美色视频一区免费| 午夜福利成人在线免费观看| 日韩在线高清观看一区二区三区| 人妻夜夜爽99麻豆av| 亚洲激情五月婷婷啪啪| 菩萨蛮人人尽说江南好唐韦庄 | 日本黄大片高清| 亚洲国产精品sss在线观看| 国产一区二区亚洲精品在线观看| 免费搜索国产男女视频| 免费搜索国产男女视频| 亚洲七黄色美女视频| 久久这里有精品视频免费| av黄色大香蕉| 国产欧美日韩精品一区二区| 国产亚洲精品av在线| 婷婷亚洲欧美| 免费av不卡在线播放| 一区二区三区高清视频在线| 国产午夜精品久久久久久一区二区三区| 久久精品国产清高在天天线| 舔av片在线| 欧美zozozo另类| 亚洲国产精品成人久久小说 | 女同久久另类99精品国产91| 日日啪夜夜撸| 亚洲av男天堂| 99九九线精品视频在线观看视频| 春色校园在线视频观看| 干丝袜人妻中文字幕| 啦啦啦韩国在线观看视频| 亚洲人成网站在线播放欧美日韩| 精品国产三级普通话版| 青春草亚洲视频在线观看| 又爽又黄无遮挡网站| 中文亚洲av片在线观看爽| 亚洲国产欧洲综合997久久,| 久久久久网色| 能在线免费观看的黄片| 97人妻精品一区二区三区麻豆| 22中文网久久字幕| 午夜精品一区二区三区免费看| 欧美三级亚洲精品| 久久久久久伊人网av| 亚洲自偷自拍三级| 毛片一级片免费看久久久久| 久久99热6这里只有精品| 国产精品电影一区二区三区| 国产精品久久久久久久电影| 一区福利在线观看| 校园人妻丝袜中文字幕| 精品无人区乱码1区二区| 中文字幕久久专区| 欧美区成人在线视频| 黄片wwwwww| 深爱激情五月婷婷| 国产精品一区www在线观看| 乱码一卡2卡4卡精品| 老熟妇乱子伦视频在线观看| 久久人人精品亚洲av| 国国产精品蜜臀av免费| 亚洲av免费高清在线观看| 精品人妻一区二区三区麻豆| av福利片在线观看| 久久精品人妻少妇| 国产精品免费一区二区三区在线| 又粗又硬又长又爽又黄的视频 | 超碰av人人做人人爽久久| 性色avwww在线观看| 国产日本99.免费观看| 一本精品99久久精品77| 搞女人的毛片| 三级国产精品欧美在线观看| 一区福利在线观看| 一级二级三级毛片免费看| 一本一本综合久久| 高清在线视频一区二区三区 | 一级毛片aaaaaa免费看小| 亚洲经典国产精华液单| 亚洲乱码一区二区免费版| 综合色av麻豆| 亚洲精品久久国产高清桃花| 日本一本二区三区精品| 美女大奶头视频| 日韩欧美精品免费久久| 成人永久免费在线观看视频| 久久久成人免费电影| 久久久国产成人免费| 国产成人一区二区在线| 长腿黑丝高跟| 亚洲国产精品久久男人天堂| 国产亚洲5aaaaa淫片| 夜夜爽天天搞| 久久亚洲国产成人精品v| 99久久精品热视频| 麻豆乱淫一区二区| 美女脱内裤让男人舔精品视频 | 99热这里只有是精品50| 黄片wwwwww| 精品欧美国产一区二区三| 免费av观看视频| 最新中文字幕久久久久| 韩国av在线不卡| av天堂中文字幕网| 少妇裸体淫交视频免费看高清| 国产在视频线在精品| 成人毛片a级毛片在线播放| 免费大片18禁| 欧美潮喷喷水| 日韩制服骚丝袜av| 男人狂女人下面高潮的视频| 国产精品综合久久久久久久免费| videossex国产| 赤兔流量卡办理| 久久这里只有精品中国| 免费看av在线观看网站| 国内少妇人妻偷人精品xxx网站| 床上黄色一级片| 久久韩国三级中文字幕| 国产黄片视频在线免费观看| 国产亚洲av嫩草精品影院| 久久久a久久爽久久v久久| 亚洲人成网站在线播放欧美日韩| 欧美日本亚洲视频在线播放| av免费观看日本| 久久精品国产自在天天线| 18禁在线播放成人免费| 亚洲欧美中文字幕日韩二区| 国产综合懂色| 特级一级黄色大片| 亚洲精品乱码久久久久久按摩| 哪里可以看免费的av片| 成人无遮挡网站| 精品久久久久久成人av| 精品少妇黑人巨大在线播放 | 免费观看a级毛片全部| 久久久精品欧美日韩精品| 在线免费十八禁| 尤物成人国产欧美一区二区三区| av天堂中文字幕网| av免费在线看不卡| 国产精华一区二区三区| 亚洲久久久久久中文字幕| 淫秽高清视频在线观看| 精品一区二区三区视频在线| 中文在线观看免费www的网站| 日韩欧美在线乱码| 亚洲成人久久性| 高清日韩中文字幕在线| 嫩草影院精品99| 最近手机中文字幕大全| 国产 一区 欧美 日韩| 免费av不卡在线播放| 日韩精品有码人妻一区| 亚洲av熟女| 欧洲精品卡2卡3卡4卡5卡区| 亚洲色图av天堂| 精品一区二区免费观看| 天天躁日日操中文字幕| 日韩欧美三级三区| 午夜老司机福利剧场| 久久久a久久爽久久v久久| a级一级毛片免费在线观看| 亚洲中文字幕日韩| or卡值多少钱| 国产 一区 欧美 日韩| 国产 一区精品| 人妻少妇偷人精品九色| 寂寞人妻少妇视频99o| 少妇猛男粗大的猛烈进出视频 | 国产精品国产三级国产av玫瑰| 一级av片app| 日韩欧美 国产精品| 91久久精品国产一区二区成人| 国产成人午夜福利电影在线观看| 有码 亚洲区| 国产av一区在线观看免费| 欧美日本亚洲视频在线播放| 亚洲人与动物交配视频| 色哟哟哟哟哟哟| 毛片一级片免费看久久久久| 九草在线视频观看| 超碰av人人做人人爽久久| 啦啦啦观看免费观看视频高清| 一区福利在线观看| 免费看av在线观看网站| 蜜臀久久99精品久久宅男| 久久九九热精品免费| 九草在线视频观看| 天美传媒精品一区二区| 久久久久久久久中文| 免费电影在线观看免费观看| 在线播放国产精品三级| 26uuu在线亚洲综合色| 午夜激情福利司机影院| 在线观看av片永久免费下载| 久久午夜亚洲精品久久| 99视频精品全部免费 在线| 校园人妻丝袜中文字幕| 国产伦一二天堂av在线观看| 我要看日韩黄色一级片| 亚洲av中文字字幕乱码综合| 精品久久国产蜜桃| 又黄又爽又刺激的免费视频.| 五月伊人婷婷丁香| 嘟嘟电影网在线观看| 日韩精品青青久久久久久| 亚洲av免费高清在线观看| 变态另类成人亚洲欧美熟女| 亚洲高清免费不卡视频| 亚洲色图av天堂| 亚洲av一区综合| 麻豆一二三区av精品| 久久久久久久久久黄片| 国产亚洲av片在线观看秒播厂 | 色综合色国产| 两个人视频免费观看高清| 日韩在线高清观看一区二区三区| 国产视频首页在线观看| 18+在线观看网站| 哪里可以看免费的av片| 99热精品在线国产| 中文精品一卡2卡3卡4更新| 99久久人妻综合| 亚洲精品日韩av片在线观看| 国产单亲对白刺激| 日本黄大片高清| 18禁在线播放成人免费| 欧洲精品卡2卡3卡4卡5卡区| 免费看光身美女| 国产爱豆传媒在线观看| 亚洲av.av天堂| 性色avwww在线观看| 亚洲国产精品国产精品| 国产午夜福利久久久久久| 国产一区二区三区在线臀色熟女| 一级av片app| av福利片在线观看| 国产男人的电影天堂91| 麻豆乱淫一区二区| 黄色日韩在线| 中文精品一卡2卡3卡4更新| 99热全是精品| 国产精品一区二区三区四区免费观看| 我的女老师完整版在线观看| 51国产日韩欧美| 成人高潮视频无遮挡免费网站| 日韩精品有码人妻一区| 亚洲av免费高清在线观看| 亚洲精华国产精华液的使用体验 | 蜜桃亚洲精品一区二区三区| 变态另类成人亚洲欧美熟女| 国产精品三级大全| 亚洲久久久久久中文字幕| 亚洲激情五月婷婷啪啪| 99热6这里只有精品| 日韩人妻高清精品专区| 热99re8久久精品国产| 欧美一区二区精品小视频在线| 久久国内精品自在自线图片| av在线蜜桃| 禁无遮挡网站| 麻豆久久精品国产亚洲av| 国内精品宾馆在线| 精品一区二区三区人妻视频| 在线免费十八禁| 我的老师免费观看完整版| 免费观看的影片在线观看| 99视频精品全部免费 在线| 午夜福利在线观看吧| 一个人看的www免费观看视频| 91久久精品国产一区二区成人| 黄色配什么色好看| 欧美成人a在线观看| 岛国毛片在线播放| 亚洲婷婷狠狠爱综合网| 熟女电影av网| 国产亚洲精品久久久com| av.在线天堂| 禁无遮挡网站| 一区二区三区高清视频在线| 午夜激情福利司机影院| 日本与韩国留学比较| 啦啦啦啦在线视频资源| 哪个播放器可以免费观看大片| 男人狂女人下面高潮的视频| 在线观看一区二区三区| 精品久久久久久久久久久久久| 日韩欧美一区二区三区在线观看| 精品人妻一区二区三区麻豆| 亚洲18禁久久av| 精品久久久噜噜| 五月玫瑰六月丁香| 天堂av国产一区二区熟女人妻| 久久精品综合一区二区三区| 国产亚洲精品av在线| 国产精品精品国产色婷婷| 欧美日本视频| 国产黄片视频在线免费观看| 网址你懂的国产日韩在线| 卡戴珊不雅视频在线播放| 国产毛片a区久久久久| 亚洲国产精品成人综合色| 三级经典国产精品| 97热精品久久久久久| 熟女电影av网| av女优亚洲男人天堂| 久久热精品热| 桃色一区二区三区在线观看| eeuss影院久久| 亚洲人成网站在线播| 久久精品国产亚洲av香蕉五月| 99久久成人亚洲精品观看| 国产精品免费一区二区三区在线| 免费观看精品视频网站| 欧美一区二区国产精品久久精品| 国产亚洲精品av在线| 欧美日韩乱码在线| 国产黄片美女视频| 午夜福利视频1000在线观看| 22中文网久久字幕| 人妻系列 视频| 2021天堂中文幕一二区在线观| 欧美丝袜亚洲另类| 村上凉子中文字幕在线| 国产91av在线免费观看| av免费在线看不卡| 亚洲丝袜综合中文字幕| 99九九线精品视频在线观看视频| 亚洲国产精品合色在线| 毛片一级片免费看久久久久| 成年女人看的毛片在线观看| 可以在线观看的亚洲视频| 精品熟女少妇av免费看| 夜夜爽天天搞| 成年女人看的毛片在线观看| 国产真实乱freesex| 久久精品人妻少妇| 久久草成人影院| 热99在线观看视频| 亚洲在线自拍视频| 综合色av麻豆| 亚洲国产精品合色在线| 亚洲精品国产av成人精品| 身体一侧抽搐| 少妇熟女欧美另类| 精品久久久久久久久亚洲| 国内久久婷婷六月综合欲色啪| 国产私拍福利视频在线观看| 久久这里有精品视频免费| 黄色视频,在线免费观看| 国产成人freesex在线| 夜夜爽天天搞| 国产精品一区二区性色av| 日韩 亚洲 欧美在线| 人人妻人人澡欧美一区二区| 精品久久久久久成人av| 国产精品人妻久久久影院| 波多野结衣巨乳人妻| 国产精品女同一区二区软件| 白带黄色成豆腐渣| 在现免费观看毛片| 久久欧美精品欧美久久欧美| 欧美另类亚洲清纯唯美| 国产一区亚洲一区在线观看| 午夜免费男女啪啪视频观看| 观看美女的网站| 国产精品一区二区在线观看99 | 精品午夜福利在线看| 哪里可以看免费的av片| 国产精品电影一区二区三区| 99九九线精品视频在线观看视频| 伦理电影大哥的女人| 国产亚洲av嫩草精品影院| 久久精品人妻少妇| 欧美成人精品欧美一级黄| 免费av毛片视频| 大又大粗又爽又黄少妇毛片口| or卡值多少钱| 国产精品国产高清国产av| 国产成人影院久久av| 成人毛片a级毛片在线播放| 桃色一区二区三区在线观看| 久久这里只有精品中国| 成人综合一区亚洲| 我要看日韩黄色一级片| 成年版毛片免费区| 亚洲最大成人手机在线| 蜜臀久久99精品久久宅男| 国产91av在线免费观看| 波多野结衣高清作品| 97超碰精品成人国产| 69av精品久久久久久| 国产精品久久久久久久电影| 非洲黑人性xxxx精品又粗又长| 深夜a级毛片| 天堂影院成人在线观看| 欧美一区二区国产精品久久精品| 天天一区二区日本电影三级| 日本一二三区视频观看| 亚洲精品日韩在线中文字幕 | 亚洲国产欧美在线一区| 在线播放无遮挡| 毛片一级片免费看久久久久| 悠悠久久av| 啦啦啦啦在线视频资源| 亚洲av.av天堂| 久久久久久久亚洲中文字幕| 国产精品国产三级国产av玫瑰| 又爽又黄a免费视频| 亚洲国产精品sss在线观看| 中文字幕精品亚洲无线码一区| 午夜福利在线观看免费完整高清在 | 国产一区二区三区在线臀色熟女| 精品久久久噜噜| 久久中文看片网| 女的被弄到高潮叫床怎么办| 久久久午夜欧美精品| 一级黄片播放器| 成人特级黄色片久久久久久久| 欧美在线一区亚洲| 日韩,欧美,国产一区二区三区 | 毛片一级片免费看久久久久| 国产毛片a区久久久久| 久久6这里有精品| 欧美不卡视频在线免费观看| 99国产精品一区二区蜜桃av| 成人高潮视频无遮挡免费网站| 精品免费久久久久久久清纯| 国产精品乱码一区二三区的特点| 亚洲av第一区精品v没综合| 色5月婷婷丁香| 超碰av人人做人人爽久久| 搞女人的毛片| 欧美激情在线99| 网址你懂的国产日韩在线| videossex国产| 中文亚洲av片在线观看爽| 毛片女人毛片| 12—13女人毛片做爰片一| 亚洲成人久久爱视频| 国产av在哪里看| 日本黄色片子视频| 精华霜和精华液先用哪个| 久久久久久久久久久丰满| 精品午夜福利在线看| 国产精品99久久久久久久久| 日本欧美国产在线视频| 亚洲国产精品sss在线观看| 99热这里只有是精品50| 国产黄色小视频在线观看| 国产极品精品免费视频能看的| 色视频www国产| 国产精品,欧美在线| 尾随美女入室| 久久久久九九精品影院| 久久人人精品亚洲av| 男女那种视频在线观看| 日本av手机在线免费观看| 色哟哟哟哟哟哟| 女的被弄到高潮叫床怎么办| 亚洲精品粉嫩美女一区| 亚洲经典国产精华液单| 99热6这里只有精品| 亚洲综合色惰| 日日摸夜夜添夜夜添av毛片| 禁无遮挡网站| av免费观看日本| 国产精品电影一区二区三区| av卡一久久| 中文字幕av成人在线电影| 亚洲不卡免费看| a级毛片免费高清观看在线播放| av福利片在线观看| 99热这里只有是精品在线观看| 美女黄网站色视频| 精品少妇黑人巨大在线播放 | 欧美最黄视频在线播放免费| 欧美zozozo另类| 久久热精品热| 精品久久久久久久久久久久久| 欧美性感艳星| 日本-黄色视频高清免费观看| 99久国产av精品| 国产成人福利小说| 国产精品1区2区在线观看.| 久久精品国产亚洲网站| 高清日韩中文字幕在线| 欧美激情国产日韩精品一区| 国产淫片久久久久久久久| 成人国产麻豆网| 亚洲成av人片在线播放无| 亚洲va在线va天堂va国产| 欧美日韩在线观看h| 精品国内亚洲2022精品成人| 六月丁香七月| 久久99蜜桃精品久久| 国产色婷婷99| 男女视频在线观看网站免费| 亚洲av中文av极速乱| 真实男女啪啪啪动态图| 亚洲av熟女| 日韩av不卡免费在线播放| 乱码一卡2卡4卡精品| 校园人妻丝袜中文字幕| 久久热精品热| 老师上课跳d突然被开到最大视频| 亚洲av一区综合| 一进一出抽搐gif免费好疼| 欧美高清性xxxxhd video| 亚洲最大成人av| 亚洲第一区二区三区不卡| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品国产成人久久av| av.在线天堂| 校园春色视频在线观看| 狂野欧美激情性xxxx在线观看| 你懂的网址亚洲精品在线观看 | 白带黄色成豆腐渣| 欧美色欧美亚洲另类二区| 欧美精品国产亚洲| 波野结衣二区三区在线| 欧美丝袜亚洲另类| 国产成人精品婷婷| av天堂中文字幕网| 人妻夜夜爽99麻豆av| 日本在线视频免费播放| 亚洲av一区综合| 天天躁夜夜躁狠狠久久av| 国产大屁股一区二区在线视频| 两个人视频免费观看高清|