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

    Preweaning period is a critical window for rumen microbial regulation of average daily gain in Holstein heifer calves

    2024-03-14 13:20:08ShengyangXuChongJiaoQiyuDiaoandYanTu

    Shengyang Xu ,Chong Jiao ,Qiyu Diao and Yan Tu*

    Abstract Background Rumen bacterial groups can affect growth performance,such as average daily gain (ADG),feed intake,and efficiency.The study aimed to investigate the inter-relationship of rumen bacterial composition,rumen fermentation indicators,serum indicators,and growth performance of Holstein heifer calves with different ADG.Twelve calves were chosen from a trail with 60 calves and divided into higher ADG (HADG,high pre-and post-weaning ADG,n=6)and lower ADG (LADG,low pre-and post-weaning ADG,n=6) groups to investigate differences in bacterial composition and functions and host phenotype.Results During the preweaning period,the relative abundances of propionate producers,including g_norank_f_Butyricicoccaceae,g_Pyramidobacter,and g_norank_f_norank_o_Clostridia_vadinBB60_group,were higher in HADG calves (LDA > 2,P<0.05).Enrichment of these bacteria resulted in increased levels of propionate,a gluconeogenic precursor,in preweaning HADG calves (adjusted P<0.05),which consequently raised serum glucose concentrations(adjusted P<0.05).In contrast,the relative abundances of rumen bacteria in post-weaning HADG calves did not exert this effect.Moreover,no significant differences were observed in rumen fermentation parameters and serum indices between the two groups.Conclusions The findings of this study revealed that the preweaning period is the window of opportunity for rumen bacteria to regulate the ADG of calves.

    Keywords Average daily gain,Calves,Preweaning period,Rumen microbiota

    Background

    Dairy cows provide high-quality milk to meet the growing demand for nutrients by the population.Therefore,rearing healthy calves for the dairy industry is critical as they form the reserve force for dairy cows.The rearing conditions and growth periods of calves are crucial phases during which any environmental changes can have long-term effects on the growth,health,and milk yield in adult life [1].A previous study has shown that preweaning liquid feed and solid feed dry matter intake synergistically improve milk,fat,and protein production.In addition,the study also revealed a positive association between first-lactation performance and preweaning ADG [2].A meta-analysis showed that for every additional 100 g/d increase in preweaning ADG,first-lactation production was expected to increase by 155 kg,and this relationship did not appear to have a plateau[3].The analysis using a test-day model revealed that the preweaning ADG accounted for 22% of the variation in first lactation milk yield [4].Together,these studies indicate that an increased preweaning ADG in calves can improve adult lactation performance.

    Rumen bacterial composition can affect feed efficiency traits,such as ADG,average daily feed intake,and gain:feed ratio,and has been shown to explain 25.3% of the variation for heifer ADG [5].Xie et al.[6] suggested that some rumen Lachnospiraceae species involved in carbohydrate metabolism contribute to the variation in residual feed intake.The findings highlighted the potential of the modulation of ruminal microbes or microbial functions to improve the feed efficiency of dairy cows [6].Xue et al.[7] demonstrated that the rumen microbiota and its metabolites contributed to different milk protein yields.These studies indicated that rumen microbes can influence the host phenotype.

    The microbial colonization of the gastrointestinal tract occurs before or during delivery,and methanogens and fibrolytic bacteria have been detected within 20 min after birth [8].Before weaning or exposure to solid feed,cellulolytic bacterium and other important bacterial species for the proper functioning of the adult rumen were detected on the first day of age [9].Afterward,the rumen ecosystem undergoes a rapid change by d 3 of age,reflected by an increased abundance of anaerobic genera,a decline in aerobic and facultative anaerobic taxa and an establishment of a new anaerobic environmental niche[9].The rumen microbes play important roles in the fermentation of the diet and energy supply of the host,which produces volatile fatty acids (VFA) to meet 70%of the daily energy need of a ruminant [10].Few studies have shown the relationships between ruminal bacteria,VFA and ADG in pre-and post-weaning calves.Therefore,understanding the relationship between ruminal bacteria and ADG in calves could help modulate ruminal bacteria to improve the preweaning ADG of calves.

    In this study,we hypothesized that calves with the same dry matter intake (DMI) and diet structure exhibit different rumen fermentation driven by rumen bacteria,leading to variations in energy levels and ultimately affecting the preweaning ADG of the calves.To test this hypothesis,the study aimed to investigate the rumen bacterial composition,rumen fermentation indicators,serum indicators,and growth performance of calves with different ADG.The goal of this study is to provide new insights into calf management.

    Materials and methods

    Animal management,diet and experimental design

    This study strictly complied with the requirements of the Animal Ethics Committee of the Chinese Academy of Agricultural Sciences (Beijing,China).Sixty Holstein heifer calves (age,14 ± 3 d;body weight,41.25 ± 2.36 kg) were kept in individual hutches till they were 90-day-old.Pasteurized milk (6 L/d) was provided to 14-60-day-old calves (preweaning) twice daily at 0630 and 1700 h.All calves were gradually weaned by reducing the dose of pasteurized milk provided each day by 1 L/d (6,5,4,3,2,1,and 1 L/d,respectively) at the age of 61-67 d (weaning period).The claves of 68 to 90-day-old (postweaning) relied solely on solid feed to meet their nutritional needs.All calves had free access to starter feed and clean water throughout the experimental period (Additional file 1: Table S1).Oat hay was fed to all calves at 80 g/d (30-60-day-old) and 200 g/d(61-90-day-old).

    The body weight (BW) of 14,60,67,and 90-dayold calves was recorded to determine the preweaning,weaning,and postweaning ADG.The initial and left-over amounts of milk,starter feed,and oat hay were recorded to calculate the DMI.The calves were grouped based on their pre-and post-weaning ADG:those with both higher pre-and post-weaning ADG were assigned to the higher ADG group (HADG;n=6),while those with both lower pre-and post-weaning ADG were assigned to the lower ADG group (LADG,n=6).We ensured that the birth weight,DMI and BW of the 14-day-old calves did not differ significantly between the two groups during the grouping process to eliminate the effect of feed and initial weight on ADG.Both groups of calves were descendants of the same bull,aiming to minimize the potential interference of genetic factors.

    Sample collection and measurement

    An oral stomach tube was used to collect rumen fluid before morning feeding in 45 and 75-day-old calves [11].The rumen fluid was collected in two tubes and stored at -80 °C for microbial DNA extraction and VFA concentration determination.The VFA concentration was determined using gas chromatography [12].Blood samples(10 mL) of each calve (30,60,and 90-day-old) were collected from the jugular vein before the morning feeding.The blood samples were then centrifuged at 2,000 ×gfor 10 min to collect the serum,which was placed in tubes and frozen at -80 °C for further analyses.The serum parameters,including total protein (TP),albumin (ALB),glucose (GLU),blood urea nitrogen (BUN),glucose(GLU) and immunoglobulins (IgA,IgG,and IgM),were analyzed using a KHB-1280 automatic biochemical analyzer (Kehua Biological Engineering Co.,Ltd.,Shanghai,China).The levels of globulin (GLB) were also calculated.The concentration of growth hormone (GH),insulin-like growth factor 1 (IGF-1) and insulin (INS) in serum were determined using bovine ELISA kits (Jinhaikeyu Biological Technology Development Co.,Ltd.,Beijing,China).

    DNA extraction,16S rRNA gene amplification and sequencing

    The microbial DNA was extracted from rumen liquid samples using the E.Z.N.A.?soil DNA Kit (Omega Biotek,Norcross,GA,USA) following the manufacturer’s instructions.The 1.0% agarose gel electrophoresis and a NanoDrop?ND-2000 spectrophotometer (Thermo Scientific Inc.,Waltham,MA,USA) were used to determine the quality and concentration of DNA,respectively.An ABI GeneAmp?9700 PCR thermocycler (Applied Biosystems,Carlsbad,CA,USA) was used to amplify the hypervariable region V3-V4 of the bacterial 16S rRNA gene with primer pairs 338F (5’-ACT CCT ACG GGA GGC AGC AG-3’) and 806R (5’-GGA CTA CHVGGG TWT CTAAT-3’) [13].The PCR reaction mixture comprised 4 μL 5× Fast Pfu buffer,2 μL 2.5 mmol/L dNTPs,0.8 μL each primer (5 μmol/L),0.4 μL Fast Pfu polymerase,10 ng template DNA,and ddH2O to a final volume of 20 μL.The PCR amplification cycling parameters were 95 °C for 3 min,27 cycles of 95 °C for 30 s,55 °C for 30 s and 72 °C for 45 s with a final extension at 72 °C for 10 min and end at 4 °C.All samples were amplified in triplicate.The 2% agarose gel was used to extract the PCR product.The PCR product was purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences,Union City,CA,USA) and quantified using Quantus?Fluorometer (Promega,Madison,WI,USA).Purified amplicons were pooled and paired-end sequenced on an Illumina MiSeq PE300 platform (Illumina,San Diego,CA,USA) at Majorbio Bio-Pharm Technology Co.,Ltd.(Shanghai,China).

    The raw FASTQ files were demultiplexed using an in-house Perl script.The raw data was quality-filtered using fastp version 0.19.6 [14] and merged using FLASH version 1.2.7 [15] according to the following criteria: (i)the 300 bp reads were truncated at any site receiving an average quality score of < 20 over a 50 bp sliding window,and the truncated reads shorter than 50 bp and reads containing ambiguous characters were discarded;(ii) overlapping sequences longer than 10 bp could be assembled,the maximum mismatch ratio of the overlap region was set at 0.2,and unassembled reads were discarded;(iii) samples were distinguished based on the barcode and primers,and the sequence direction was adjusted,exact barcode matching,2 nucleotide mismatch in primer matching.Then the optimized sequences were used to cluster operational taxonomic units (OTUs) using UPARSE 7.1 with a 97% sequence similarity level [16].The RDP Classifier version 2.2 was used to analyze the taxonomy of each OTU representative sequence against the 16S rRNA gene database using a confidence threshold of 0.7 [17].Metagenomic functions were predicted through the utilization of Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2),in conjunction with the KEGG (Kyoto Encyclopedia of Genes and Genomes) database [18].

    Statistical analysis

    The ADG,DMI,serum parameters,rumen fermentation parameters,and bacterial taxonomic and functional data were analyzed using SAS software (version 9.4,SAS Institute Inc.,Cary,NC,USA).A mixed procedure was used as follows:

    whereYjikis the dependent variable,μis the overall mean,Giis the effect of theithgroup,Pj(Dj) is the effect of thejthperiod (day),GPij(GDij) is the interaction between group and period (group and day),Ckis the random effect ofkthcalf,andeijkis the residual error.

    The independent samplet-test was performed to analyze birth weight and BW of 14-day-age using the SAS software.The taxonomic and functional data were analyzed on the online platform of the Majorbio Cloud Platform [19].Rumen microbial data’s alpha diversity between groups was assessed using the Wilcoxon rank-sum test,with a false discovery rate (FDR)adjustedP-value <0.05 indicating significance.Principal co-ordinates analysis (PCoA) was applied to assess beta diversity using the Bray-Curtis distance algorithm,while intergroup differences were tested using the adonis method.The taxonomic and functional data were also compared using Linear Discriminant Analysis Effect Size (LEfSe);significance was determined with aP-value < 0.05 and a LDA score > 2.The LEfSe was also used to compare the KEGG enzymes with |log2Fold Change| > 1,LDA > 2 andP-value <0.05 considered significant.Spearman’s rank correlations were employed to unveil the associations among rumen bacteria at the genus level,with a FDR adjustedP-value <0.05 and a Spearman’s |r|> 0.50 deemed as statistically significant correlations.The Gephi software (version 0.10) was used for visualization.Spearman’s rank correlations were also constructed between ADG,GLU30 (serum glucose in 30-day-old calves),GLU60 (serum glucose in 60-day-old calves),GH30 (serum growth hormone in 30-day-old calves),GFI-1 30 (serum insulin-like growth factor 1 in 30-day-old calves),INS30 (serum insulin in 30-day-old calves),propionate,total VFA,microorganism and KEGG enzymes during the preweaning period(Spearman’s |r|> 0.50 and FDR adjustedP-value < 0.05).Multiplex networks were visualized using Gephi software (version 0.10).

    Results

    Intake,growth performance and serum parameters

    The HADG group had significantly higher ADG than the LADG group during the pre-and post-weaning periods(P<0.05) (Fig.1A).However,the total DMI,stater DMI,birth weight and BW of 14-day-old calves did not differ between the two groups (P>0.05;Fig.1A).Furthermore,the milk and forage given in a quantified amount and were fully consumed by the calves,which resulted in similar forage DMI,and milk DMI between the two groups during the tested periods (P>0.05;Fig.1A).These findings suggest that despite consistent birth weight,BW at 14 days of age,total DMI,and diet composition,variations in microbial-mediated ruminal fermentation might cause the differences in ADG.

    Fig.1 Phenotypes of calves in the HADG and LADG groups at different periods.A Differences in feed intake,body weight and ADG between HADG and LADG calves in different periods.B Differences in serum indices of HADG and LADG calves in different periods.HADG,higher average daily gain group;LADG,lower average daily gain group;* represents significantly different between the two groups (P<0.05).# represents a tendency to differ between the two groups (0.05<P<0.10)

    As shown in Fig.1B,no changes in other serum parameters were observed between the two groups (P>0.05).However,the HADG calves exhibited elevated concentrations of serum GLU30,GLU60,GH30,and IGF-1 30 during the preweaning period compared to those in LADG calves (P<0.05).Additionally,INS30 also showed a tendency to be higher in the serum of HADG calves during the preweaning period (P<0.10).These findings collectively suggest an increase in hormones related to growth and an augmented energy supply in HADG calves(Fig.1B).

    Ruminal fermentation parameters concerning VFA

    Remarkably,higher concentrations of total VFA and propionate were detected in the rumen of HADG calves than those in LADG calves during the preweaning period (P<0.05) (Fig.2).However,no significant differences in fermentation parameters were observed between the two groups during the postweaning period(P>0.05) (Fig.2).

    Fig.2 Rumen fermentation in HADG and LADG groups.HADG,higher average daily gain group;LADG,lower average daily gain group;* represents significantly different between the two groups (P<0.05)

    Bacterial diversities and composition

    Alpha diversity calculations revealed no marked divergence of the Shannon and Chao indices (adjustedP>0.05),indicating unchanged bacterial richness and evenness between the two groups at any stage (Fig.3AD).The principal coordinate analysis (PCoA) showed no clear separation between the HADG and LADG groups at any time point (Fig.3E and F).

    Fig.3 Alpha and beta diversity between HADG and LADG groups in different periods.A and B Alpha diversity (Chao and Shannon indices)of rumen bacteria between HADG and LADG groups during the preweaning period.C and D Alpha diversity (Chao and Shannon indices) of rumen bacteria between HADG and LADG groups during the postweaning period.E Beta diversity of a principal coordinate analysis (PCoA) based on Bray-Curtis dissimilarity matrices between HADG and LADG groups during the preweaning period.F Beta diversity of a principal coordinate analysis(PCoA) based on Bray-Curtis dissimilarity matrices between HADG and LADG groups during the postweaning period.HADG,higher average daily gain group;LADG,lower average daily gain group

    At the phylum level,the predominant phyla were Firmicutes (61.70%),Bacteroidota (21.42%),Actinobacteriota (15.30%) and Patescibacteria (0.72%) across the 12 samples from the preweaning period (Fig.4A;Additional file 1: Table S2).Meanwhile,none of the phyla showed significant differences between the HADG and LADG groups.At the genus level,the dominant genera across the 12 samples before weaning wereLachnospiraceae_NK3A20_group(15.09%),Olsenella(14.82%),Ruminococcus(14.18%),Prevotella(13.61%) andShuttleworthia(3.38%) (Fig.4B;Additional file 1: Table S3).Several genera that differed significantly between the two groups were identified;the abundance of most of these genera was higher in the HADG group than that in the LADG group.Specifically,the relative abundances ofnorank_f_Butyricicoccaceae,Pyramidobacter,norank_f_norank_o_Clostridia_vadinBB60_group,Blautia,andUCG-004were enriched in the HADG group,whereas the relative abundances ofSolibacillusandAtopobiumwere enriched in the LADG group (LDA > 2,P<0.05)(Fig.4C;Additional file 1: Table S4).The support provided by the mixed model further enhances the quality of our results.The results showed that the relative abundances ofnorank_f_Butyricicoccaceaeandnorank_f_norank_o_Clostridia_vadinBB60_groupwere found to be enriched in the HADG group,while the relative abundances ofSolibacillusandAtopobiumwere enriched in the LADG group (P<0.05).Pyramidobacter,Blautia,andUCG-004also exhibited a tendency to be higher in the HADG group (P<0.10) (Additional file 2: Fig.S1).

    Fig.4 Bacterial composition and differential taxa between HADG and LADG groups in different periods.A and B The relative abundances of the rumen bacterial community at the phylum and genus levels during the preweaning period.D and E The relative abundances of the rumen bacterial community at the phylum and genus levels during the postweaning period.C and F The significantly different bacterial taxa between HADG and LADG groups during the preweaning (C) and postweaning period (F).HADG,higher average daily gain group;LADG,lower average daily gain group

    Next,we focused on the samples (n=12) from the postweaning period.The predominant phyla were Firmicutes(77.01%),Actinobacteriota (18.03%),and Bacteroidota(0.72%) (Fig.4D;Additional file 1: Table S5).No significantly different taxa were found at the phylum level.At the genus level,the dominant genera across the 12 samples from postweaning includedErysipelotrichaceae_UCG-002(20.73%),Acetitomaculum(17.99%),Olsenella(15.75%) andnorank_f__norank_o__Clostridia_UCG-014(6.77%) (Fig.4E;Additional file 1: Table S6).However,the abundance of none of these genera differed significantly between the HADG and LADG groups.On the contrary,the abundance ofLachnospiraandLachnoclostridiumwas significantly higher in the LADG group than that in the HADG group during the postweaning period(LDA > 2,P<0.05) (Fig.4F;Additional file 1: Table S7).The mixed model also confirms these results (Additional file 2: Fig.S1).

    Microbial interactions

    Genera enriched in the rumen of HADG calves had more connections in the network during the preweaning period,suggesting a richer metabolic association and a more critical role in the network.Conversely,the higher abundance ofSolibacillusin the LADG group during the preweaning period was found to be on the network’s periphery with fewer connections,suggesting a minor role in the overall network (Fig.5).

    Fig.5 Network analysis to reveal microbial interactions during the preweaning period.The network analysis showed the degree of correlation between the bacteria at the genus level (Spearman’s |r|> 0.50 and adjusted P<0.05).The node colors represent the phylum classification of the genera.Lines between two nodes represent the correlation,with a red line indicating a positive correlation and a blue line indicating a negative correlation.The green text represents the significantly higher genera in the HADG group;the blue text represents the significantly higher genera in the LADG group (LDA > 2,P<0.05).HADG,higher average daily gain group;LADG,lower average daily gain group

    During the postweaning period,onlyLachnospiraappeared to be located at the edge of the network.Negative correlations were observed betweenLachnospiraand other taxa (Additional file 2: Fig.S2).

    Functions of microbiota

    The software PICRUSt2 was used to predict microbial function.At KEGG pathway level 1,seven categories,including Metabolism (47.07%),Genetic Information Processing (20.65%),Environmental Information Processing (14.72%),Unclassified (13.03%),Cellular Processes (3.00%),Organismal Systems (0.79%),and Human Diseases (0.74%),were identified in the 12 samples during the preweaning period (Fig.6A).Among these,carbohydrate (11.34%),amino acid (9.57%),and energy(5.47%) metabolism were the most active in the metabolism pathway (Fig.6B).The top 10 pathways for carbohydrate and amino acid metabolism were further observed at KEGG pathway level 3 (Fig.6C).However,only eight pathways at KEGG level 3 were detected in the energy metabolism pathway (Fig.6C).Additionally,several KEGG enzymes related to carbohydrate metabolism and propanoate metabolism were also significantly enriched in preweaning HADG calves (log2FC > 1,LDA > 2 andP<0.05) (Fig.6D;Additional file 1: Table S8).The mixed model can yield identical results as well (P<0.05) (Additional file 2: Fig.S3).

    Fig.6 Functional predictions of rumen bacteria between HADG and LADG groups during preweaning period.A-C Relative abundance of function KEGG pathway level 1,level 2 and level 3 (top 10 of carbohydrate metabolism,top 10 of amino acid metabolism,and top 8 of energy metabolism).D Significantly different KEGG enzymes related to carbohydrate and propionate metabolisms between HADG and LADG groups (LDA > 2,P<0.05).HADG,higher average daily gain group;LADG,lower average daily gain group

    The KEGG enzymes related to carbohydrate and propionate metabolisms in the rumen of the preweaning calves were summarized in Fig.7.Among the KEGG enzymes involved in starch and sucrose metabolism pathway,one enzyme,EC 3.2.1.93,was more abundant in the HADG calves than those in LADG group (log2FC > 1,LDA > 2 andP<0.05;Fig.6D).EC 1.2.1.59 in the glycolysis pathway,EC 1.1.1.140 in the fructose and mannose metabolism pathway,EC 1.2.7.5 in pentose phosphate pathway,and EC 2.8.3.1 in the propanoate metabolism pathway were more abundant in the HADG calves (log2FC > 1,LDA > 2 andP<0.05;Fig.6D).

    Next,we focused on the postweaning calves.At the KEGG pathway level 1,6 categories were observed in the 12 samples during the postweaning period,including Metabolism(76.95%),Genetic Information Processing (8.57%),Environmental Information Processing (6.08%),Cellular Processes(3.79%),Human Diseases (3.00%),and Organismal Systems(1.61%) (Additional file 1: Table S9).In the Metabolism pathway,global and overview maps (41.13%),carbohydrate metabolism (10.81%) and amino acid metabolism (7.09%)were the most active (Additional file 1: Table S10).The top 10 pathways at the KEGG pathway level 3 were involved with carbohydrate and amino acid metabolisms,wherein 8 pathways were identified under the global and overview maps pathway (Additional file 1: Table S11).None of the KEGG pathways differed significantly.Furthermore,no significant differences were observed in the KEGG enzymes related to carbohydrate and propionate metabolisms.These results suggest that HADG calves lose the advantages imparted by carbohydrate and propionate metabolisms during the postweaning period.

    Relationships between phenotypes,microorganisms,and microbial functions

    We investigated the relationships between various indicators by constructing spearman’s rank correlations between ADG,GLU30,GLU60,GH30,INS30,IGF-1 30,propionate,total VFA,microorganisms,and KEGG enzymes during the preweaning period.The correlations were visualized by multiplex networks (Fig.8;Additional file 1: Table S12).Norank_f_Butyricicoccaceaepositively correlated with thePyramidobacterandnorank_f_norank_o_Clostridia_vadinBB60_group(adjustedP<0.05).In addition,we identified a few relationships between bacteria and KEGG enzymes.Specifically,the relative abundancesnorank_f_Butyricicoccaceaehad positive correlations with the relative abundances of KEGG enzymes related to carbohydrate and propionate metabolisms (EC 1.1.1.140,EC 1.2.1.59,EC 2.8.3.1,and EC 3.2.1.93) (adjustedP<0.05).The relative abundance ofnorank_f_norank_o_Clostridia_vadinBB60_groupshowed positive correlations with the relative abundances of EC 1.1.1.140,EC 1.2.1.59,and EC 3.2.1.93 (adjustedP<0.05).ThePyramidobactergenus showed positive relationships with EC 1.1.1.140,EC 2.8.3.1,and EC 3.2.1.93 (adjustedP<0.05).The concentrations of total VFA and propionate were positively correlated with EC 1.2.1.59,EC 3.2.1.93,and EC 1.2.7.5 (adjustedP<0.05).Propionate is the main precursor for ruminant liver gluconeogenesis,and glucose is the main energy supplier for living organisms.In our study,total VFA and propionate were positively correlated with INS30,GLU30 and GLU60(adjustedP<0.05),wherein serum GLU30 and INS30 showed a positive correlation with ADG (adjustedP<0.05).

    Discussion

    In this study,we systematically compare ruminal bacteria and host phenotypes of calves with different ADG in pre-and post-weaning periods.Our study sheds light on the role of ruminal bacteria in driving variation in preweaning ADG.Our data indicate that rumen bacteria play a pivotal role in propionate and total VFA production,leading to elevated serum glucose and insulin levels,ultimately contributing to a higher preweaning ADG.

    The differences in total DMI were likely responsible for most of the improvement in preweaning ADG,as the higher DMI can provide the necessary energy and protein needed for growth [20].Additionally,diet structure can also affect the growth performance of calves,with liquid feeds being shown to meet the energy needs of calves better than solid feeds [21].A previous study has also shown that increased milk intake positively affects calf growth during preweaning [22].However,in our study,we did not observe any differences in total DMI or diet structure between the two groups,suggesting other factors are responsible for the difference in preweaning ADG.

    Further analysis of serum indicators revealed that serum GLU,INS,GH,and IGF-1 were higher in HADG calves than that in LADG calves during the preweaning period,and the these differences between the two groups disappeared during the postweaning period.Glucose is an essential energy-producing nutrient for all mammals and is important for the survival and growth of cattle[23].Blood glucose concentration in ruminants is closely related to the glucose turnover rate,meaning that a high blood glucose level can lead to increased glucose utilization [24].Insulin is an important metabolic hormone that regulates glucose concentration and stimulates the uptake and utilisation of glucose by tissues [25].Growth hormone and its closely associated hormone,insulin-like growth factor 1,play a pivotal role in stimulating bodily growth and advancement by exerting their effects on vital metabolic organs such as the liver,skeletal muscle,and bones [26].Hence,elevated levels of glucose,growth hormone,insulin,and insulin-like growth factor 1 in serum potentially supply greater energy to calves,enhancing glucose utilization and potentially fostering bodily growth and development.This may ultimately contribute to variations in preweaning ADG.

    The reasons for the differences in serum parameters still need to be further explored.In terms of rumen fermentation indicators,we found important clues.During the preweaning period,the HADG calves had higher total VFA and propionate concentrations.Approximately 70% of the host’s energy needs can be met by VFA energy produced by rumen microbes [27].Cattle primarily obtain glucose from ruminal VFAs through gluconeogenesis [28].Undoubtedly,higher total VFA concentrations are beneficial for increasing host energy levels and promoting host growth and development.Propionate,valerate and isobutyrate can be used as precursors of hepatic glucose production for the net synthesis of glucose and as the main source of glucose and energy for ruminants[27,29,30].Propionate is the main precursor substance for hepatic gluconeogenesis,with a contribution rate of 60%-74% [28].In this study,the accumulation of propionate in the rumen of HADG calves during the preweaning period resulted in elevated serum glucose.

    Ruminal bacteria are responsible for driving ruminal fermentation,which produces VFAs [31].Hence,analyzing rumen bacteria is necessary to investigate the reasons for changes in VFAs.Similar to previous studies,Firmicutes and Bacteroidota were the two major bacterial phyla detected in preweaning calves [32,33].Firmicutes have been shown to be dominant in high-grain diets,while the Bacteroidota are dominant in hay diets [34].In this study,preweaning calves were primarily fed starter and less hay,which may have led to the dominance of the Firmicutes phylum.Almost all differential genera(norank_f_Butyricicoccaceae,Pyramidobacter,norank_f_norank_o_Clostridia_vadinBB60_group,Blautia,andUCG-004) had higher abundance in the rumen of HADG calves during the preweaning period.

    The Butyricicoccaceae family is an important butyrate producer [35].Our study demonstrated thatnorank_f_Butyricicoccaceaelikely possess an acid-producing function,leading to differences in rumen fermentation.Pyramidobactercontains many acid-producing strains that produce acetate,isovalerate,propionate,isobutyrate,succinate,and phenylacetic acid as end products [36,37].A Previous study revealed thatPyramidobactermay synthesize thiamine (vitamin B1) in the rumen to meet the needs of ruminants [38].Additionally,thePyramidobactergenus was found to be associated with rumen development [39].In short,thePyramidobactergenus not only can produce acid but also may produce other substances beneficial to the host,which may promote the growth and development of the host.

    The Clostridiales vadinBB60 group is known to consist of conditional pathogens [40].Under normal conditions,a balance is maintained between conditional pathogens and their hosts and other flora through factors such as nutritional competition and mutual restraint of metabolites.However,this equilibrium relationship is disrupted under certain conditions,and the conditional pathogens may show a negative impact.This suggests that conditional pathogens may also perform beneficial functions.Zhang et al.[41] suggested that the unclassified Clostridiales vadinBB60 group was positively correlated with propionate.In this study,the enrichment ofnorank_f_norank_o_Clostridia_vadinBB60_groupmay also be a reason for the increased propionate concentration.

    Blautiais a genus of probiotic bacteria widely distributed in mammalian feces and intestines and expresses enzymes related to propionate production [42,43].Blautiais also positively related to plasma concentrations of amino acids,biogenic amines,and sphingomyelins,as well as ADG in calves [44].In terms of microbial networks,norank_f_Butyricicoccaceae,Pyramidobacter,norank_f_norank_o_Clostridia_vadinBB60_group,andBlautiahave more connections with other bacterial taxa,especiallyPyramidobacter,occupies a relatively important position in the network.Enrichment of these four genera in the rumen of HADG calves may have increased the metabolic efficiency of the rumen flora [45],resulting in a greater impact on the rumen flora.

    The relative abundance ofSolibacilluswas increased in the LADG calves during the preweaning period.A higher predisposition to diarrhea may be associated with theSolibacillus[46].This indicates that the enrichment of theSolibacillusis detrimental to the growth and development of the host.Nevertheless,the network analysis identifiedSolibacillus at the edge of the microbial network indicating its less connection to other taxa and less influence on the network.

    During the postweaning period,Firmicutes and Actinobacteria emerged as the dominant bacterial phyla.On the contrary,previous studies have reported Bacteroidetes and Firmicutes as the main bacterial phylum in postweaning calves [47].This discrepancy could be attributed to the different dietary structures in the studies,which may have contributed to the enrichment of Actinobacteria.All the differential genera (LachnospiraandLachnoclostridium) were enriched in postweaning LADG calves,suggesting that the HADG calves have lost the advantage of acidogenic bacteria,as evidenced by ruminal fermentation parameters during the postweaning period.Lachnospiracomprises many pectin degradation bacteria and has been shown to be enriched in the rumen of grass-fed steers [48,49].Lachnoclostridiumconsists of many beneficial bacteria [50].However,an increase in the relative abundances ofLachnospiraandLachnoclostridiumwas insufficient to alter ruminal fermentation in LADG calves.

    The KEGG enzymes predicted by PICRUSt2 provide a way to evaluate the functions of rumen microorganisms.During the preweaning period,five enzymes (EC 2.8.3.1,EC 1.2.7.5,EC 1.1.1.140,EC 1.2.1.59,and EC 3.2.1.93)related to carbohydrate and propionate metabolisms were enriched in the HADG calves.Propionyl-CoA carboxylase (EC 2.8.3.1),a key enzyme,played a vital role in propionate formation via propanoate metabolism[51].Aldehyde ferredoxin oxidoreductase (EC 1.2.7.5)is involved in the pentose phosphate pathway,and its catalyzed downstream metabolites are involved in the Glycolysis.Sorbitol-6-phosphate 2-dehydrogenases (EC 1.1.1.140) catalyze the interconversion ofD-sorbitol-6-phosphate toD-fructose 6-phosphate,which can be utilized by bacteria to utilize sorbitol [52].Glyceraldehyde-3-phosphate dehydrogenase (EC 1.2.1.59) promote the degradation of carbohydrates in the glycolysis pathway [53].Phosphotrehalase (EC 3.2.1.93) promotes sucrose degradation in the starch and sucrose metabolism.Overall,the enrichment of the five KEGG enzymes may promote the degradation of carbohydrates in the fructose and mannose metabolism,starch and sucrose metabolism,pentose phosphate pathway,and glycolysis pathway.They also support propanoate metabolism by synthesizing lactate from pyruvate,resulting in an increased propionate concentration.However,during the postweaning period,the two groups of calves did not differ significantly in KEGG enzymes involved in carbohydrate degradation and propionate metabolism.

    Multiplex networks revealed that rumen microbes regulate ADG during the preweaning period.The abundance of three genera (norank_f_Butyricicoccaceae,Pyramidobacter,andnorank_f_norank_o_Clostridia_vadinBB60_group) led to the enrichment of three KEGG enzymes(EC 1.2.1.59,EC 1.2.7.5,and EC 3.2.1.93) related to carbohydrate degradation,which in turn increased the concentrations of propionate and total VFA.Propionate can be used as a substrate for hepatic gluconeogenesis to produce glucose and promote serum glucose concentrations.Higher serum glucose concentrations not only provide calves with more energy but also increase the efficiency of serum glucose utilization in calves,ultimately resulting in higher ADG.Concurrently,heightened serum glucose levels prompt the body to generate insulin,thereby further amplifying glucose utilization and providing the body with increased energy resources.

    It is important to emphasize that the HADG calves maintained higher ADG even during the postweaning period,but this was not due to differences in rumen fermentation driven by microorganisms.In other words,the window of microbial regulation of ADG may disappear during the postweaning period,but the advantage of higher ADG caused by microorganisms during the preweaning period will persist for a long time.

    Conclusion

    In conclusion,this study highlights the preweaning period as a critical time window for rumen bacteria to regulate ADG (Fig.9).Specifically,the abundance of three key genera (norank_f_Butyricicoccaceae,Pyramidobacter,andnorank_f_norank_o_Clostridia_vadinBB60_group) with propionate-producing and carbohydrate-degrading abilities in preweaning calves led to higher concentrations of propionate and total VFA.The bacteria-driven production of propionate and total VFA were strongly correlated with preweaning serum glucose and insulin,which in turn were positively associated with preweaning ADG.These correlation networks are crucial in regulating multiple host phenotypes during preweaning,particularly preweaning ADG.It is important to note that although the window period for rumen bacteria to regulate ADG disappears after weaning,the beneficial effects of bacteria on ADG can persist in the postweaning period.Further studies are needed to confirm these findings and to explore the potential long-term effects of these differences on the health and productivity of the animals.

    Fig.9 The preweaning period is the window of opportunity for rumen bacteria to regulate the average daily gain (ADG).VFA,volatile fatty acids

    Abbreviations

    ADG Average daily gain

    ALB Albumin

    BUN Blood urea nitrogen

    DMI Dry matter intake

    GH30 Serum growth hormone in 30-day-old calves

    GLB Globulin

    GLU30 Serum glucose in 30-day-old calves

    GLU60 Serum glucose in 60-day-old calves

    HADG Higher average daily gain

    IGF-1 30 Serum insulin-like growth factor 1 in 30-day-old calves

    INS30 Serum insulin in 30-day-old calves

    KEGG Kyoto Encyclopedia of Genes and Genomes

    LADG Lower average daily gain

    TP Total protein

    VFA Volatile fatty acids

    Supplementary Information

    The online version contains supplementary material available at https://doi.org/10.1186/s40104-023-00934-0.

    Additional file 1: Table S1.The composition and nutrient levels of starter diets (DM basis),%.Table S2.The relative abundances of the top 10 phyla between HADG and LADG groups during the preweaning period.Table S3.The relative abundances of the top 10 genera between HADG and LADG groups during the preweaning period.Table S4.The relative abundances of significantly different genera between HADG and LADG groups during the preweaning period.Table S5.The relative abundances of the top 9 phyla between HADG and LADG groups during the postweaning period.Table S6.The relative abundances of the top 10 genera between HADG and LADG groups during the postweaning period.Table S7.The relative abundances of significantly different genera between HADG and LADG groups during the postweaning period.Table S8.The relative abundances of significantly different KEGG enzymes between HADG and LADG groups during the preweaning period.Table S9.The relative abundances of KEGG pathways level 1 between HADG and LADG groups during the postweaning period.Table S10.The relative abundances of KEGG pathways level 2 between HADG and LADG groups during the postweaning period.Table S11.The relative abundances of KEGG pathways level 3 between HADG and LADG groups during the postweaning period.Table S12.Correlation analysis to reveal interactions among the significantly different bacterial genera,enzymes,fermentation parameters,serum indicators and ADG during the preweaning period.

    Additional file 2: Fig.S1.The mixed model reveals significantly different bacterial genera.Fig.S2.Network analysis to reveal microbial interactions during the postweaning period.Fig.S3.The mixed model reveals significantly different KEGG enzymes.

    Acknowledgements

    The authors would like to thank reviewers for their kind suggestions and comments on this work.

    Authors’ contributions

    SX contributed to study design,data analyses,and writing the manuscript.CJ contributed to sample collection and sample processing.QD and YT contributed to project administration,study design,and supervision.

    Funding

    This study was funded by National Key R&D Program of China(2022YFA1304204),Agricultural Science and Technology Innovation Program(CAAS-ASTIP-2017-FRI-04) and Beijing Innovation Consortium of livestock Research System (BAIC05-2023).

    Availability of data and materials

    All data involved in this study are included in this article and its Supplementary files.

    Declarations

    Ethics approval and consent to participate

    The Animals involved in this experiment were cared for according to the guidelines of the Animal Ethics Committee of the Chinese Academy of Agricultural Sciences.The committee reviewed and approved the experiments and procedures involving the use of Animals.

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Received:15 May 2023 Accepted:1 September 2023

    久久久欧美国产精品| 国产精品精品国产色婷婷| 国产成人福利小说| 国产精品一区二区性色av| 特级一级黄色大片| 国产成人一区二区在线| 秋霞在线观看毛片| 国产精品蜜桃在线观看| 午夜精品国产一区二区电影 | 99九九线精品视频在线观看视频| 一区二区三区四区激情视频| 亚洲国产高清在线一区二区三| 人人妻人人爽人人添夜夜欢视频 | 精品一区二区免费观看| 自拍欧美九色日韩亚洲蝌蚪91 | 国产大屁股一区二区在线视频| 国产精品一区www在线观看| 一区二区三区免费毛片| av国产精品久久久久影院| 一区二区三区免费毛片| 在线观看国产h片| 激情五月婷婷亚洲| 午夜激情久久久久久久| 爱豆传媒免费全集在线观看| 亚洲自拍偷在线| 国产真实伦视频高清在线观看| 少妇人妻一区二区三区视频| 欧美日韩综合久久久久久| 身体一侧抽搐| 婷婷色麻豆天堂久久| 亚洲熟女精品中文字幕| 国产女主播在线喷水免费视频网站| 最近中文字幕2019免费版| 波野结衣二区三区在线| 亚洲人成网站高清观看| 国产色爽女视频免费观看| 狂野欧美激情性xxxx在线观看| 免费观看性生交大片5| 亚洲精品国产色婷婷电影| av国产精品久久久久影院| 免费大片18禁| 中国国产av一级| 一级a做视频免费观看| 免费高清在线观看视频在线观看| 久久热精品热| 国产精品久久久久久久电影| 人人妻人人澡人人爽人人夜夜| tube8黄色片| 亚洲欧美成人精品一区二区| 久久精品国产鲁丝片午夜精品| 精品久久久久久久人妻蜜臀av| 欧美高清成人免费视频www| 国产毛片在线视频| 国产毛片a区久久久久| videos熟女内射| 国产精品人妻久久久久久| 夜夜爽夜夜爽视频| 亚洲美女搞黄在线观看| 亚洲精品日本国产第一区| 国产精品一二三区在线看| 成年人午夜在线观看视频| 精品久久久久久久末码| 亚洲欧洲日产国产| 亚洲婷婷狠狠爱综合网| 国产亚洲午夜精品一区二区久久 | 51国产日韩欧美| 婷婷色av中文字幕| 综合色丁香网| 麻豆成人午夜福利视频| 精品久久国产蜜桃| 搡老乐熟女国产| 最近最新中文字幕免费大全7| 少妇的逼好多水| 三级国产精品欧美在线观看| 国产成人精品久久久久久| 日产精品乱码卡一卡2卡三| 免费大片黄手机在线观看| 日本免费在线观看一区| av网站免费在线观看视频| 亚洲精品一二三| 日韩视频在线欧美| 十八禁网站网址无遮挡 | 舔av片在线| 国产av码专区亚洲av| 亚洲成人精品中文字幕电影| 久久精品国产自在天天线| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品乱码久久久久久按摩| 99热国产这里只有精品6| 99久久九九国产精品国产免费| 午夜福利高清视频| 国产永久视频网站| 欧美日韩在线观看h| www.色视频.com| 国产视频首页在线观看| 嫩草影院新地址| 日日摸夜夜添夜夜爱| 美女cb高潮喷水在线观看| 欧美丝袜亚洲另类| 欧美xxxx性猛交bbbb| 99精国产麻豆久久婷婷| av网站免费在线观看视频| 久久久久久久久久久免费av| 只有这里有精品99| av国产免费在线观看| 国产伦精品一区二区三区视频9| 91久久精品国产一区二区三区| 日韩成人av中文字幕在线观看| 国语对白做爰xxxⅹ性视频网站| 69人妻影院| 亚洲欧美日韩另类电影网站 | 少妇裸体淫交视频免费看高清| 丰满人妻一区二区三区视频av| 日本av手机在线免费观看| 我要看日韩黄色一级片| 国产女主播在线喷水免费视频网站| 久久久亚洲精品成人影院| 日日撸夜夜添| av在线天堂中文字幕| 久久久国产一区二区| 久久99热这里只频精品6学生| 97在线视频观看| 欧美成人a在线观看| 简卡轻食公司| 亚洲精品影视一区二区三区av| 麻豆精品久久久久久蜜桃| 美女cb高潮喷水在线观看| 国产精品久久久久久精品古装| 国产91av在线免费观看| 成人一区二区视频在线观看| 国产日韩欧美亚洲二区| 中文资源天堂在线| 免费观看av网站的网址| 日韩在线高清观看一区二区三区| 国产男女超爽视频在线观看| 国产色婷婷99| 三级男女做爰猛烈吃奶摸视频| 亚洲色图av天堂| 亚洲成人精品中文字幕电影| 熟女av电影| 黄色配什么色好看| 久久精品夜色国产| 大又大粗又爽又黄少妇毛片口| 国产亚洲一区二区精品| 日韩视频在线欧美| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩东京热| 天堂中文最新版在线下载 | 国产精品国产三级国产av玫瑰| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品一区二区免费观看| 男女边吃奶边做爰视频| 天天躁夜夜躁狠狠久久av| 街头女战士在线观看网站| 久久99蜜桃精品久久| 九九久久精品国产亚洲av麻豆| 日韩制服骚丝袜av| 午夜福利视频精品| 亚洲va在线va天堂va国产| 国产久久久一区二区三区| 亚洲av中文字字幕乱码综合| 国产精品久久久久久久电影| 亚洲精品自拍成人| 国产精品一及| 国产亚洲91精品色在线| 免费观看a级毛片全部| 久久影院123| 中文天堂在线官网| 视频区图区小说| 国产av国产精品国产| 亚洲四区av| 国产精品一区www在线观看| 乱码一卡2卡4卡精品| 亚洲人成网站在线观看播放| 小蜜桃在线观看免费完整版高清| 99热6这里只有精品| 中文字幕久久专区| 少妇高潮的动态图| 看黄色毛片网站| 高清毛片免费看| 少妇的逼水好多| 久久人人爽人人片av| 久热久热在线精品观看| 国产淫片久久久久久久久| 亚洲精品第二区| 午夜爱爱视频在线播放| 久久精品国产鲁丝片午夜精品| 亚洲精品456在线播放app| 成年免费大片在线观看| av在线播放精品| 18禁在线无遮挡免费观看视频| 在线免费观看不下载黄p国产| 夜夜爽夜夜爽视频| 国产成人a区在线观看| 男人狂女人下面高潮的视频| 国产精品99久久久久久久久| 制服丝袜香蕉在线| 日本wwww免费看| 精华霜和精华液先用哪个| 欧美日韩视频精品一区| 视频中文字幕在线观看| 亚洲国产最新在线播放| 97在线人人人人妻| 亚洲av不卡在线观看| 97超碰精品成人国产| 成人亚洲欧美一区二区av| 国产成人午夜福利电影在线观看| 国产av码专区亚洲av| 国产午夜福利久久久久久| 天天躁夜夜躁狠狠久久av| 国产熟女欧美一区二区| 国产一区二区三区av在线| 亚洲精品国产成人久久av| 精品一区在线观看国产| 国产 精品1| 亚洲欧美清纯卡通| 中文在线观看免费www的网站| 黄色欧美视频在线观看| 精品国产一区二区三区久久久樱花 | 欧美最新免费一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人中文字幕在线播放| 国产成年人精品一区二区| 一个人观看的视频www高清免费观看| 一级黄片播放器| 国产日韩欧美亚洲二区| 日本一本二区三区精品| 交换朋友夫妻互换小说| 国内揄拍国产精品人妻在线| 卡戴珊不雅视频在线播放| 成人毛片a级毛片在线播放| 搡女人真爽免费视频火全软件| 亚洲精品亚洲一区二区| 久久影院123| 听说在线观看完整版免费高清| 香蕉精品网在线| 大香蕉久久网| 少妇人妻一区二区三区视频| 亚洲aⅴ乱码一区二区在线播放| 五月天丁香电影| 91久久精品国产一区二区成人| 国产精品一二三区在线看| 亚洲欧美日韩东京热| 18+在线观看网站| 在线天堂最新版资源| 亚洲欧美精品自产自拍| 国国产精品蜜臀av免费| 大香蕉97超碰在线| 亚洲av中文字字幕乱码综合| 国产永久视频网站| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜爱| 最近最新中文字幕大全电影3| 伦理电影大哥的女人| 最近最新中文字幕免费大全7| 亚州av有码| 午夜激情福利司机影院| 国产精品人妻久久久久久| 精品人妻视频免费看| 一二三四中文在线观看免费高清| 久久精品国产亚洲av天美| 成人免费观看视频高清| 亚洲国产精品专区欧美| 一区二区三区乱码不卡18| 国产色婷婷99| 纵有疾风起免费观看全集完整版| 嫩草影院精品99| 又爽又黄无遮挡网站| 日本熟妇午夜| 国内揄拍国产精品人妻在线| 婷婷色综合www| 国产免费又黄又爽又色| 在线观看av片永久免费下载| 狂野欧美激情性xxxx在线观看| 国产精品蜜桃在线观看| 免费av不卡在线播放| 国产精品国产三级国产av玫瑰| 最近2019中文字幕mv第一页| 欧美变态另类bdsm刘玥| 国产精品偷伦视频观看了| 久久久久久久精品精品| 久久精品国产亚洲av天美| 国产午夜福利久久久久久| 国国产精品蜜臀av免费| 极品教师在线视频| 免费观看av网站的网址| 高清av免费在线| 又爽又黄无遮挡网站| 亚洲av二区三区四区| 国产男人的电影天堂91| 亚洲国产色片| 国产 一区 欧美 日韩| 超碰97精品在线观看| 久久女婷五月综合色啪小说 | 六月丁香七月| 99热这里只有是精品50| av网站免费在线观看视频| 青春草亚洲视频在线观看| 性色av一级| 天天一区二区日本电影三级| 97人妻精品一区二区三区麻豆| 成人无遮挡网站| 国产欧美日韩精品一区二区| 国产一区二区三区av在线| 欧美成人a在线观看| 成人欧美大片| 亚洲精品视频女| 欧美区成人在线视频| 久久久国产一区二区| 女人被狂操c到高潮| 99精国产麻豆久久婷婷| 亚洲国产欧美在线一区| 久久精品久久久久久久性| 五月开心婷婷网| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 秋霞伦理黄片| 日本与韩国留学比较| 97热精品久久久久久| 日韩一区二区视频免费看| 亚洲av国产av综合av卡| 久久精品久久精品一区二区三区| 国产一区有黄有色的免费视频| 神马国产精品三级电影在线观看| 舔av片在线| 国产在线一区二区三区精| 中文字幕av成人在线电影| 天美传媒精品一区二区| 又黄又爽又刺激的免费视频.| 国产极品天堂在线| 嫩草影院入口| 日韩视频在线欧美| 亚洲天堂av无毛| 久久韩国三级中文字幕| 国产高清三级在线| 久久99蜜桃精品久久| 国产在视频线精品| 另类亚洲欧美激情| 少妇人妻 视频| 国产亚洲精品久久久com| 涩涩av久久男人的天堂| 免费高清在线观看视频在线观看| 最后的刺客免费高清国语| 免费人成在线观看视频色| 在线观看一区二区三区| 水蜜桃什么品种好| 亚洲丝袜综合中文字幕| 男的添女的下面高潮视频| 激情 狠狠 欧美| 五月天丁香电影| 欧美激情国产日韩精品一区| 国产成年人精品一区二区| 啦啦啦啦在线视频资源| 国产毛片a区久久久久| 免费不卡的大黄色大毛片视频在线观看| 干丝袜人妻中文字幕| 亚洲第一区二区三区不卡| 激情五月婷婷亚洲| 中文字幕免费在线视频6| 欧美精品人与动牲交sv欧美| 嫩草影院入口| 午夜视频国产福利| 99久久中文字幕三级久久日本| 干丝袜人妻中文字幕| 成年女人看的毛片在线观看| 少妇的逼好多水| 美女高潮的动态| 夜夜看夜夜爽夜夜摸| 国产黄片视频在线免费观看| 青春草视频在线免费观看| 3wmmmm亚洲av在线观看| 乱系列少妇在线播放| 国产精品一区二区在线观看99| 亚洲av免费在线观看| 欧美日韩亚洲高清精品| 日韩人妻高清精品专区| 午夜福利视频精品| 欧美日本视频| 18禁裸乳无遮挡动漫免费视频 | 免费播放大片免费观看视频在线观看| 水蜜桃什么品种好| 国产有黄有色有爽视频| 综合色丁香网| 欧美潮喷喷水| 久久久成人免费电影| 又爽又黄无遮挡网站| 日韩欧美 国产精品| av国产久精品久网站免费入址| 欧美老熟妇乱子伦牲交| 国产精品人妻久久久影院| 久久精品久久久久久久性| 午夜福利网站1000一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 国产在线男女| 丰满人妻一区二区三区视频av| 欧美成人精品欧美一级黄| 男女边摸边吃奶| 久久国内精品自在自线图片| 成人免费观看视频高清| 一级二级三级毛片免费看| 联通29元200g的流量卡| 国产成人91sexporn| 少妇 在线观看| 成人毛片60女人毛片免费| 国产高清有码在线观看视频| 人人妻人人看人人澡| 亚洲国产欧美在线一区| 中文天堂在线官网| 亚洲精品乱码久久久久久按摩| 亚洲第一区二区三区不卡| 菩萨蛮人人尽说江南好唐韦庄| av国产久精品久网站免费入址| 久久久久久久大尺度免费视频| 久久女婷五月综合色啪小说 | av黄色大香蕉| 日本wwww免费看| 日本三级黄在线观看| av在线蜜桃| 亚洲自拍偷在线| 能在线免费看毛片的网站| 搞女人的毛片| 久久久久久久午夜电影| 99热这里只有是精品50| 国产精品秋霞免费鲁丝片| 99热全是精品| 亚洲第一区二区三区不卡| 好男人视频免费观看在线| 在线精品无人区一区二区三 | 精品人妻熟女av久视频| 一级毛片黄色毛片免费观看视频| 国产高清有码在线观看视频| 欧美高清成人免费视频www| 日韩av免费高清视频| 日本一本二区三区精品| 久久久久久伊人网av| 18禁动态无遮挡网站| 精品人妻视频免费看| 秋霞在线观看毛片| 国内少妇人妻偷人精品xxx网站| 性色avwww在线观看| 成人免费观看视频高清| 欧美精品国产亚洲| 欧美3d第一页| 免费观看av网站的网址| 精品久久久精品久久久| 久久久久久久久久人人人人人人| 嫩草影院新地址| 在线精品无人区一区二区三 | 激情五月婷婷亚洲| 国产伦精品一区二区三区四那| 国产精品人妻久久久影院| 97在线人人人人妻| 国产伦精品一区二区三区视频9| 成人亚洲精品av一区二区| 欧美精品人与动牲交sv欧美| 91精品伊人久久大香线蕉| 成人高潮视频无遮挡免费网站| 国产精品精品国产色婷婷| 日韩欧美一区视频在线观看 | 老司机影院毛片| 国产视频首页在线观看| 国产乱人偷精品视频| 婷婷色综合大香蕉| 中文乱码字字幕精品一区二区三区| 国产免费福利视频在线观看| 2021少妇久久久久久久久久久| 成人亚洲欧美一区二区av| 一级毛片久久久久久久久女| 波多野结衣巨乳人妻| 免费不卡的大黄色大毛片视频在线观看| 国产爽快片一区二区三区| 亚洲av国产av综合av卡| 国产成人aa在线观看| 天天躁夜夜躁狠狠久久av| 国产淫片久久久久久久久| 久久久久国产网址| 丰满人妻一区二区三区视频av| 又爽又黄a免费视频| 日韩av不卡免费在线播放| 亚洲精品国产色婷婷电影| 精品久久久噜噜| 黑人高潮一二区| 丝瓜视频免费看黄片| 日韩中字成人| 久久久久久久大尺度免费视频| 日韩欧美精品v在线| 色视频www国产| 亚洲真实伦在线观看| 午夜激情福利司机影院| 日韩av不卡免费在线播放| 在线 av 中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 国产亚洲最大av| 大又大粗又爽又黄少妇毛片口| 国产一区亚洲一区在线观看| 精品亚洲乱码少妇综合久久| 日韩精品有码人妻一区| 99热这里只有精品一区| 亚洲天堂av无毛| 国产精品伦人一区二区| 街头女战士在线观看网站| 女人被狂操c到高潮| 日本午夜av视频| 亚洲人成网站在线播| 成人国产麻豆网| 老师上课跳d突然被开到最大视频| 欧美少妇被猛烈插入视频| 天天一区二区日本电影三级| 一本一本综合久久| 久久国产乱子免费精品| 免费黄频网站在线观看国产| 99九九线精品视频在线观看视频| 人体艺术视频欧美日本| 人妻系列 视频| 国产精品av视频在线免费观看| 成年女人在线观看亚洲视频 | 免费观看av网站的网址| 啦啦啦啦在线视频资源| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品一二三| av国产久精品久网站免费入址| 久久精品人妻少妇| 国产欧美日韩一区二区三区在线 | 成人午夜精彩视频在线观看| 久久韩国三级中文字幕| 在线观看av片永久免费下载| 免费黄网站久久成人精品| 大码成人一级视频| 国产成人freesex在线| 亚洲av免费高清在线观看| 欧美日韩在线观看h| 中文在线观看免费www的网站| 久久久精品94久久精品| 色综合色国产| 一本一本综合久久| 欧美日韩精品成人综合77777| 久久久久久久亚洲中文字幕| 深爱激情五月婷婷| 大陆偷拍与自拍| 亚洲av免费高清在线观看| 国内精品美女久久久久久| 人妻少妇偷人精品九色| 日韩av在线免费看完整版不卡| 久久久久久久大尺度免费视频| 国产一区亚洲一区在线观看| 国产成人freesex在线| 久久久久网色| 精品久久久久久久久亚洲| 联通29元200g的流量卡| 亚洲色图综合在线观看| 国产免费一级a男人的天堂| 韩国高清视频一区二区三区| 日韩三级伦理在线观看| 成人高潮视频无遮挡免费网站| av女优亚洲男人天堂| 欧美日韩精品成人综合77777| 免费人成在线观看视频色| 亚洲aⅴ乱码一区二区在线播放| 国内少妇人妻偷人精品xxx网站| 欧美成人午夜免费资源| 搡老乐熟女国产| 在线a可以看的网站| 午夜激情福利司机影院| 99久久中文字幕三级久久日本| 伦精品一区二区三区| 色婷婷久久久亚洲欧美| 久久久久精品久久久久真实原创| 小蜜桃在线观看免费完整版高清| 午夜老司机福利剧场| 2022亚洲国产成人精品| 真实男女啪啪啪动态图| 啦啦啦中文免费视频观看日本| 国产成人精品久久久久久| 特大巨黑吊av在线直播| 97在线视频观看| 国产老妇女一区| 女人被狂操c到高潮| 国产av不卡久久| 欧美激情久久久久久爽电影| 99久国产av精品国产电影| 日韩伦理黄色片| 免费看不卡的av| 亚洲国产最新在线播放| 少妇人妻一区二区三区视频| 成年女人在线观看亚洲视频 | 欧美精品人与动牲交sv欧美| 久久久久久久久久久丰满| 国产精品精品国产色婷婷| 大话2 男鬼变身卡| 亚洲精品乱码久久久v下载方式| 嫩草影院精品99| 亚洲欧美精品自产自拍| 99热这里只有精品一区| 嫩草影院精品99| 身体一侧抽搐| 久久精品国产亚洲av天美| 深爱激情五月婷婷| 国产高潮美女av| 美女视频免费永久观看网站| 97在线视频观看| 夜夜看夜夜爽夜夜摸| 中文字幕制服av| 舔av片在线| 国产精品一区二区性色av| 亚洲图色成人| 欧美日韩视频高清一区二区三区二| 亚洲第一区二区三区不卡| 国产爽快片一区二区三区| 一级二级三级毛片免费看| 欧美精品一区二区大全| 国产日韩欧美亚洲二区| 91在线精品国自产拍蜜月| 亚洲激情五月婷婷啪啪| 97在线人人人人妻| 啦啦啦中文免费视频观看日本| 亚洲av国产av综合av卡|