INTRODUCTION
Ruminants rely on a diverse and robust rumen microbiota for the efficient digestion of plant fibers, which are converted into volatile fatty acids (VFAs), microbial proteins, ammonia and other metabolites essential for growth [
1]. The rumen contains a diverse array of microorganisms, such as bacteria, protozoa, and fungi, which co-evolve with the host to influence the phenotype and carry out vital functions in metabolism, digestion, and immunity [
2]. Diet, environment, and host genetics were the primary drivers of the observed variations in the rumen microbiota community structure. Animal performance traits, such as average daily gain and feed efficiency have also been associated with specific microbial states in the rumen, indicating that microbiota can influence livestock productivity [
3].
Gender and age, fundamental characteristics of individual organisms, significantly influence the composition of the gut microbiota community [
4]. Gender differences may lead to variations in hormone levels, dietary preferences, and metabolic functions, thereby impacting the structure of the gut microbiota [
5]. Studies have indicated that females tend to have higher intestinal microbial diversity and a greater abundance of antibiotic-resistance genes than males [
6]. Additionally, host phenotypes, such as body mass index (BMI) and diastolic blood pressure, inversely correlate with the core microbiota, which is influenced by gender [
7]. Notably, research conducted on non-obese diabetic (NOD) mouse models revealed that transplanting male colonic contents profoundly reshaped the gut microbiota of female NOD mice. This intervention alleviated pancreatic inflammation and reduced insulin antibody production [
8], suggesting a gender-specific role of the gut microbiota in modulating the occurrence of host diseases. Furthermore, age significantly influences the dynamic assembly of the gut microbiota. Changes in rumen structure and physiological characteristics correlate with the development of the rumen microbiota [
9]. Studies have demonstrated significant changes in the composition of rumen microbiota between 1-day-old and 3-day-old calves, despite their identical diets and nearly identical symptoms [
10]. Following calf birth, aerobic bacteria in the rumen rapidly deplete the available oxygen, prompting a shift from aerobic to anaerobic conditions in the microbial ecological environment. It has been found that age-related microbial changes are closely related to host inflammation [
11]. Consequently, changes in ecological niches lead to a restructuring of the microbial community by the rumen microbiota. The initial bacterial communities colonizing the rumen differ significantly from those found in mature animals. Increasing age contributes to enhanced complexity and stability within these communities [
10]. Additionally, females exhibited higher carcass fat and back fat content compared to males, and intramuscular fat content significantly increased with age [
12]. Studies have indicated that tail fat deposition was negatively correlated with the abundance of Verrucomicrobiota and
Bacteroides in the colon of Small-Tailed Han sheep [
13], while the abundance of
Moryella in the rumen was positively correlated with intramuscular fat content in Angus× Simmental steers [
14]. These studies have demonstrated that the gut microbiota is influenced by the gender and age of the host, and is closely associated with host phenotypes such as fat deposition. However, our understanding of the gender and age disparities in the composition and functionality of rumen microbiota among ruminants is still limited due to various constraints, such as animal breeds, feeding management practices, geographical locations, and environmental conditions.
Therefore, we conducted a comparative analysis of Qinchuan cattle with uniform genetic backgrounds and consistent feeding management practices, varying only in gender and age, to explore their respective impacts on ruminal fermentation function and microbial community structure. Our study lays the groundwork for improving ruminant production efficiency by strategically regulating rumen microbiota. It also offers valuable insights into conserving and strategically utilizing genetic resources in Chinese local yellow cattle.
MATERIALS AND METHODS
Animal care
The management and handling of experimental animals are carried out in accordance with the Animal Welfare and Ethical Rules of the Laboratory Animal Management Committee of Northwest A&F University of Science and Technology in China and have been approved (protocol number: NWAFUCAST2018-168).
Experiment design and sample collection
The study was conducted at the Qinchuan Cattle Breeding Farm of the National Beef Cattle Improvement Center. A total of 38 QinChuan beef cattle, aged 3 months (6 males and 6 females), 9 months (7 males and 7 females), and 15 months (6 males and 6 females), with good health and consistent feeding practices, were selected. The experimental animals, grouped by gender and age, were all maintained on a same diet for three months. The diet consisted of corn silage, wheat straw, and concentrate, fed twice daily at 6:00 AM and 4:00 PM. The feed was manually mixed according to the formulation and provided to the experimental animals. The feeding amount was set to 110% of the previous intake to ensure feed remained each day. The concentrate was formulated according to the Beef Cattle Feeding Standard (NY/T 81-2004). The composition and nutrient levels of the diet were detailed in
Table 1. Clean and ample drinking water was provided, and standard immunization procedures were followed.
Rumen fluid and blood samples were collected when the cattle were 6 months old (QC6), 12 months old (QC12), and 18 months old (QC18). Rumen fluid was collected two hours after the morning feeding using an oral rumen tube, then filtered through gauze, and divided into 5 mL centrifuge tubes. Samples were immediately frozen in liquid nitrogen and stored at −80°C for subsequent determination of rumen fermentation performance and microbial 16S rRNA sequencing. Blood samples were obtained from the jugular vein of cattle on an empty stomach using conventional vacutainer tubes. The samples were then incubated at 25°C for 4–6 h, and the serum was subsequently separated by centrifugation at 3,000×g for 15 min. Serum samples, each 500 μL, were divided into 1.5 mL Eppendorf tubes and stored at −80°C for the analysis of serum metabolites.
Chemical analysis
The chemical composition of the feed sample was determined after drying at 65°C for 72 h and crushing it through a 1 mm screen. Nitrogen content was determined using the Kjeldahl method with copper sulphate and potassium sulphate (1:10, w/w) as a catalyst (UDK159, VELP, Usmate, MB, Italy). Crude protein (CP) content was calculated as N×6.25. The content of neutral detergent fiber (NDF) and acid detergent fiber (ADF) was analyzed using a fibre analyser (F800, Hanon, Jinan, China), following the method described by Van Soest et al [
15]. Ash content was measured using a muffle furnace at 550°C for 5 h, with preliminary ashing in an electric heating panel (F47910-33, Thermo Scientific, Waltham, MA, USA). Organic matter (OM) content was calculated by subtracting the ash content from 100%. The ether extract (EE) content was determined using an automatic Soxhlet extractor (SOX606, Hanon).
The concentration of VFAs in the rumen fluid was determined using gas chromatography-flame ionization detection (GC-FID) (internal standard method, GC-6850, Agilent, Santa Clara, CA, USA), with crotonic acid as serving as the internal standard. Ammonia nitrogen (NH3-N) concentration was measured by using an Ultraviolet-visible spectrophotometer (colorimetric method; Cary 60, Agilent). The concentration of serum metabolites was assessed using an automatic chemistry analyzer (BK-400, Biobase, Jinan, China).
DNA extraction and high-throughput sequencing
The TIANamp Stool DNA Kit (Tiangen Biochemical Technology Co., Ltd., Beijing, China) was used to extract the total DNA of bacterial in rumen fluid. The V3–V4 region of the 16S rRNA genes in the total DNA was amplified using specific primers (F: ACTCCTACGGGAGGCAGCA; R: GGACTACHVGGGTWTCTAAT), with a sequencing adapter added to the end of the primer before PCR amplification. PCR products were then purified, quantified, and normalized to construct a sequencing library. Following quality assessment, the amplified products were sequenced on the Illumina NovaSeq 6000 platform (Illumina, San Diego, CA, USA), with the sequencing procedures and raw data quality control entrusted to Biomic Biotechnology Co., Ltd. (Beijing, China).
Bioinformatics analysis
After filtering, removing primer sequences and chimeras, the DADA2 algorithm [
16] within QIIME2 2020.6 [
17] was used to denoise the data, applying a threshold of 0.005% to filter amplicon sequence variants (ASVs). The SILVA database (Release138,
http://www.arb-silva.de) served as the reference database, and the Naïve Bayes classifier annotated feature sequences, classified species, and filtered features with less than 2% sequence abundance. Alpha diversity (Chao1 and Shannon index) and beta diversity were assessed using QIIME2 2020.6 [
17]. Principal Coordinates Analysis (PCoA) based on Bray-Curtis dissimilarity visualized differences in the bacterial community among groups, followed by permutational multivariate analysis of variance (PERMANOVA) to test the significance of these differences. Gender and age differences in rumen microorganisms at the phylum and genus levels were analyzed using STAMP (v2.1.3) [
18] via Welch’s t-test, with unclassified taxa and relative abundances less than 0.1% being filtered out. Furthermore, functions from the Kyoto Encyclopedia of Genes and Genomes (KEGG) were predicted across different genders and ages using PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) [
19]. Redundancy analysis (RDA) in CANOCO5 was utilized to examine the correlations among serum metabolites, rumen ammonia and VFAs production and dominant bacteria in the rumen.
Statistics and Analysis
To investigate the changes in serum metabolites, rumen ammonia, and VFAs production variations in gender and age, a general linear model (GLM) procedure in SPSS 26.0 was employed. The analysis considered the fixed effects of gender and age, as well as the interaction between gender and age. The Pearson correlation coefficients (r) were also calculated using SPSS 26.0. Data were expressed as mean±standard error. Differences between groups were considered significant at p<0.05, while a p-value between 0.05 and 0.1 indicated a significant trend.
DISCUSSION
The production of ammonia nitrogen through the microbial deamination of feed proteins is essential for the digestion and utilization of nutrients in the rumen [
20]. In this study, rumen ammonia nitrogen concentrations in Qinchuan cattle showed no significant differences across genders or ages. However, studies on Holstein cattle indicated a significantly higher rumen ammonia nitrogen concentration in males compared to females [
21]. This is potentially due to differences in breed and diet composition.
VFAs are the primary energy source for ruminants, produced through the sequential degradation of feed nutrients by ruminal microbes and providing over 70% of their energy needs [
1]. Acetate, propionate, and butyrate account for approximately 95% of the total VFAs. The proportions of acetate and propionate were significantly influenced by dietary composition, while the proportions of butyrate and valerate were less affected by diet [
22]. Increasing the proportion of non-structural carbohydrates (NFC) in the diet typically decrease the acetate ratio and increase the propionate ratio [
23]. In this study, all experimental cattle were fed the same diet, and no significant gender differences were observed in the concentrations of acetate and propionate. However, the proportion of acetate showed a trend toward increasing with age, potentially due to variations in dry matter intake (DMI). The acetate/propionate ratio did not show significant differences across genders or ages in Qinchuan cattle, suggesting similar ruminal fermentation patterns across these groups. Moreover, butyrate concentrations were significantly higher in females than in males, indicating clear gender differences. The concentrations of isobutyrate, butyrate, and isovalerate also increased with age (
Table 1). Butyrate is produced by ruminal microbes through the fermentation of carbohydrates in the feed, such as insoluble cellulose and resistant starch, and is absorbed by the epithelial cells within the rumen. Most of the butyrate was subsequently converted into ketone bodies (β-HB), which serve as substrates for the synthesis of body fat [
24]. An increase in butyrate concentration stimulated the expansion of the surface of rumen epithelial cells, and promoted the oxidation of short-chain fatty acids in the ketogenesis pathway [
25]. Variations in butyrate concentration by gender and age may be attributed to differences in fat deposition capacities between males and females and to varying rates of fat deposition among different age groups of Qinchuan cattle. Isobutyrate and isovalerate, branched-chain VFAs produced from the fermentation of branched-chain amino acids, are crucial for the degradation of structural carbohydrates and for microbial protein synthesis in the rumen [
26]. Supplementing branched-chain VFAs has been shown to enhance the production of ruminal microbes and increase cellulase activity in sheep [
27]. In this study, the increase in branched-chain fatty acid concentrations in QC
18 was likely related to the rise in total VFA concentrations.
In general, greater microbial community diversity correlates with enhanced stability of the community structure [
28]. Current research suggested that sex hormones were a primary factor causing differences in gut microbiota between genders [
29]. Research involving mice has shown that male mice had lower gut microbial diversity compared to female mice, influenced by hormonal variations. Castration eliminated these gender differences, while hormone replacement reinstated them, suggesting a direct hormonal effect on microbial composition [
30]. However, research on older humans, such as postmenopausal women and elderly men, showed no significant gender differences in gut microbial diversity [
31], suggesting that the observed gender differences in younger populations were influenced by fluctuating sex hormone levels. In our study, we observed no significant differences in the rumen microbial diversity between male and female Qinchuan cattle. This might be due to rumen microbial ecosystem in cattle might have responded differently to hormonal changes compared to the human or murine gut. The unique physiology of ruminants could lead to distinct interactions between hormones and microbiota. However, since the levels of sex hormones were not measured in this study, further research is needed to provide supporting evidence to elucidate the impact of sex hormones on rumen microbial diversity. Furthermore, rumen microbial diversity varied across different developmental stages. Studies on dairy cows have shown that both rumen microbial diversity and community similarity increased with age [
10], which was consistent with our findings, indicating that age was one of the most important factors in the stability of rumen microbial community structure.
Firmicutes and Bacteroidetes were the two dominant phyla with the highest abundance in the rumen [
32]. In this study, significant differences were observed in the abundance of these dominant phyla between males and females. The abundance of Bacteroidetes was significantly lower, whereas Firmicutes was significantly higher in females than in males (
Figure 3). The ratio of Firmicutes to Bacteroidetes is crucial for maintaining normal gut homeostasis, and an increased ratio has been closely associated with the development of human obesity [
33]. In addition, increased abundance of Firmicutes have been shown to be associated with higher feed efficiency and greater body fat accumulation in cattle [
34,
35]. This suggested that Firmicutes play a significant role in the regulation of body fat and the development of obesity. Generally, females have a higher percentage of body fat than males [
12], which may be related to the increase in Firmicutes in females. Furthermore, Desulfobacteria is one of the main sulfate-reducing bacteria in the rumen, and the accumulation of its metabolites in the rumen epithelium may induce intestinal inflammation [
36]. Akkermansia muciniphila, a member of Verrucomicrobia found in mammalian intestines, colonizes the intestinal mucosa and exhibits anti-inflammatory properties that protect the host from intestinal pathogens [
37]. In our study, the decrease in Desulfobacteria in males may be related to the increase in Verrucomicrobia. Notably, there were also significant differences in Desulfobacteria and Verrucomicrobia among different age groups. Compared to QC
6, Desulfobacteria significantly increased in QC
12, while Verrucomicrobia exhibited a significant increase in QC
18. These findings indicate dynamic changes in gut microbiota across different developmental stages of the host and highlight the beneficial role of Verrucomicrobia in promoting intestinal health.
Significant gender differences were observed in the abundance of
Prevotella,
Ruminococcus,
Saccharofermentans,
Succiniclasticum and
NK4A214_group.
Prevotella, which is widely present in the mammalian gastrointestinal tract, plays a crucial role in degrading carbohydrates through the fermentation of various plant polysaccharides [
38].
Prevotella primarily produces propionate as a fermentation end-product. Once absorbed through the rumen wall, propionate serves as an important precursor for hepatic gluconeogenesis, playing a critical role in maintaining GLU homeostasis [
39]. However, in our study, the propionate concentration in the rumen of male Qinchuan cattle did not show a significant change with the increase in
Prevotella abundance. This might be related to the higher GLU utilization efficiency in males, meaning that the propionate produced by
Prevotella was quickly converted to GLU and were utilized [
40]. The significant enrichment of microbial functions related to glycan biosynthesis and metabolism in males also supported the finding. Butyrate was closely related to feed efficiency in ruminants, and researchers have observed a significant increase in butyrate concentration in Holstein Friesian dairy cows with high feed efficiency [
41]. The higher TG levels were consistent with the higher body fat percentage and fat storage capacity typically observed in females [
42]. Furthermore, females oxidize a greater proportion of lipids (EFA) for energy compared to carbohydrates and proteins than males. Additionally, females store more circulating free fatty acid (FFA) into adipose tissue during rest periods than males [
43]. This aligns with our findings that the rumen microbiota in male Qinchuan cattle were significantly enriched in carbohydrate metabolism, while females were enriched in lipid metabolism. Moreover,
Ruminococcus exhibited a significant positive correlation with the concentration of butyrate, and the levels of serum TG and NEFA (
Figure 6). The higher levels of butyrate, TG, NEFA and increased abundance of
Ruminococcus in females indicated that gender differences in feed efficiency and fat metabolism may be partly attributed to variations in the rumen microbiome composition, where
Ruminococcus play key roles.
Ruminococcus,
UCG-005,
NK4A214_group,
Saccharofermentans, members of the Ruminococcaceae family, were famous for their robust fiber decomposition capabilities and were also primary butyrate producers [
44]. The significant increase in these microbes in females suggested enhanced fiber degradation capabilities, enabling increased fermentation of plant fibers, which aligned with findings in Tibetan sheep [
45]. In addition, we identified age-related differences in microbial abundance. In QC
6, there was a notable increase in the abundance of
Pseudobutyrivibrio, which is capable of fermenting various carbohydrates and producing butyrate as a key end product [
46]. However, the concentration of butyrate in QC
6 did not significantly improve. This aligns with previous findings, showing that variations in
Pseudobutyrivibrio abundance did not correspond with changes in butyrate concentration [
47]. The rumen microbial community exhibits significant functional redundancy, with multiple microorganisms capable of degrading the same substrates [
48]. Alterations in community composition often do not lead to substantial changes in VFA concentrations [
48,
49]. The fluctuations in abundance of
Pseudobutyrivibrio may not be sufficient to influence butyrate levels in the rumen. Notably, we did not observe any significant changes in the abundance of other butyrate-producing bacteria, suggesting that further investigation is necessary to fully understand these findings.