Alaa Emara Rabee1*
Ibrahim Abou-souliman2Ahmed I. Yousif3
Mebarek Lamara4*Mohamed A. El-Sherbieny3Eman Ali Elwakeel5
Ahmed M. Sallam2- 1Animal and Poultry Nutrition Department, Desert Research Center, Cairo, Egypt
- 2Animal and Poultry Breeding Department, Desert Research Center, Cairo, Egypt
- 3Animal Production Research Institute, Agricultural Research Center, Cairo, Egypt
- 4Forest Research Institute, University of Quebec in Abitibi-Temiscamingue, Rouyn-Noranda, QC, Canada
- 5Department of Animal and Fish Production, Faculty of Agriculture, University of Alexandria, Alexandria, Egypt
Introduction: Improving feed efficiency (FE) is a significant goal in animal breeding programs. Variations in FE and its relationship with rumen microbiota remain poorly understood across different goat breeds.
Methods: This study assessed the influence of the rumen microbiome and host genome on FE in 10 Shami (SH) goats and 10 Zaraibi (ZA) goats, all of which were fed the same diet. The animals were genotyped using the Illumina 65KSNP chip arrays v2, and their rumen bacteria and archaea were investigated using PCR-amplicon sequencing of the 16S rRNA gene.
Results: The results showed that the ZA goats exhibited higher FE than the SH goats (p < 0.05) and a greater relative abundance (p < 0.05) of rumen bacterial groups that are involved in the degradation of cellulose and hemicelluloses, such as Bacteroidota and Fibrobacterota, along with genera such as Prevotella, Lachnospiraceae AC2044 group, Lachnospiraceae NK3A20 group, and Succiniclasticum, which are linked to carbohydrate and nitrogen metabolism pathways. In addition, breed differences were found to affect the proportions of milk fatty acids. The association analysis identified 26 genome windows containing several putative candidate genes, such as TMEM241, AP4S1, FTO, HYAL2, BBS2, CD52, CRYBG2, PIGV, WDTC1, EEF1A2, GBA2, FNIP1, ACSL6, STARD10, VPS26B, ACAD8, GLB1L3, NRN1L, LCAT, and SLC7A6. These genes contributed to FE traits in Egyptian goats, as they are involved in obesity, metabolism, and the transport of energy, vitamins, fatty acids, proteins, and lipids through diverse biological pathways.
Discussion: This study suggests that specific genetic markers and rumen microbial traits could be used to identify high-efficient individuals in Egyptian goat breeds, and improving breeding strategies for FE.
Introduction
Increased milk and meat production in ruminant animals, such as goats, is associated with the depletion of feed resources and methane (CH4) emissions, highlighting the importance of enhancing productivity while mitigating environmental impacts (Xue et al., 2022).
Feed efficiency (FE) is closely associated with microbial fermentation, which depends on high-efficient microbial groups, including bacteria, archaea, protozoa, and fungi (Li and Guan, 2017; Løvendahl et al., 2018).
These microbial groups work together to break down and transform the ingested animal diet into microbial protein (MP) and volatile fatty acids (VFAs), which fulfill the host animal’s energy and protein needs (Firkins and Yu, 2015). Also, microbial fermentation generates hydrogen (H2) and carbon dioxide (CO2) that are used, besides other substrates, such as acetic and formic acid and methanol, as substrates to produce methane by rumen archaea (Difford et al., 2018). Methane represents a loss in the animal’s gross energy intake and contribute to greenhouse gas emissions (Rabee et al., 2022).
Understanding the relationship between the rumen microbiome and FE can improve current breeding programs by facilitating the selection of highly efficient animals (Xue et al., 2022).
Additionally, recent studies demonstrated the association between the host genome and the rumen microbial community in cattle (Li et al., 2019). Therefore, effective strategies to improve FE and reduce methane emissions in livestock should focus on both animal genomes and rumen microbiomes (Difford et al., 2018).
Animals with higher FE exhibited lower methane emissions, a lower relative abundance of Methanobrevibacter, and higher levels of Candidatus Methanomethylophilus archeae and propionic acid-producing bacteria (Tapio et al., 2017; Bharanidharan et al., 2018; McLoughlin et al., 2020). Previous studies on rumen microbiota indicated that the bacterial community in goats was affiliated mainly with the phyla Bacteroidetes and Firmicutes, and the dominant genera were Prevotella, Butyrivibrio, and Ruminococcus (Giger-Reverdin et al., 2020; Luo et al., 2022), whereas the gut archaeal community was dominated by the genera Methanobrevibacter and Methanosphaera (Luo et al., 2022).
Feed constitutes the largest expense in livestock production, making selective breeding a key strategy to reduce costs (Singh et al., 2022; Khanal et al., 2023). Thus, identifying genetically superior animals with high FE has become integral to breeding programs. The recent advances in high-throughput sequencing technologies have helped breeders collect adjunct information about the genetic potentiality of their animals to be included in genomic selection (Weigel, 2017). The best animals (i.e., animals that eat less feed while maintaining similar production to their herd mates) are then selected for breeding to transmit their genes of high FE to the next generations (Khanal et al., 2023). Genome-wide association studies (GWAS) use hundreds of thousands of genetic markers, typically single-nucleotide polymorphisms (SNPs), that spread out across the entire genome to identify genetic markers for the trait of interest, which supports genomic selection in livestock (Desire et al., 2018; Moreira et al., 2019; Taussat et al., 2020; Barría et al., 2021; Uffelmann et al., 2021; Khanal et al., 2023).
Zaraibi (also known as Egyptian Nubian) and Shami (also known as Damascus) goats are the most important goat breeds in Egypt due to their high potential for fertility and milk production under harsh conditions (Eid et al., 2020; Almasri et al., 2023). Therefore, animals of the two breeds were involved in genetic improvement programs for exotic breeds in several countries (Güney et al., 2006; Tatar et al., 2021). The rumen microbiome and host genome of goat breeds in arid regions received less attention. Furthermore, few studies explored the association between the rumen microbiome, host genome, and FE in ruminant animals. Therefore, this study aimed to investigate the relationship between FE, rumen microbiota, and host genome in Shami and Zaraibi goat breeds.
Materials and methods
Ethics statement
All animal procedures included in the current study were approved by the Animal Breeding Ethics Committee at the Desert Research Center (DRC) in Egypt (reference number: AB/NO2022) and the Research Ethics Review Committee, Faculty of Agriculture, University of Alexandria, Egypt (Reference: Alex. Agri. 082305307). All methods and protocols in this study comply with the ARRIVE 2.00 guidelines. The experiment did not include animal euthanasia, and all animals were released to the goat herd after the end of the experiment.
Animals and sampling
The experiment was conducted at the Agriculture Research Station, Animal Production Research Institute, Agricultural Research Center (ARC), Sakha, Kafr El Sheikh, Egypt. A total of 20 lactating goats were selected to represent both breeds, ensuring sufficient statistical power for comparative analysis. The animals comprised Shami goats (SH) with an average body weight of 40.6 ± 1.39 kg (n = 10) and Zaraibi goats (ZA) with an average body weight of 25.3 ± 1.57 kg (n = 10). The experiment lasted 45 days, with a 15-day adaptation period followed by 30 days of data collection. The adaptation period aimed to adapt the animals to individual housing conditions and feeding systems. All the animals used in this study were the offspring from the Sakha Agriculture Research Station goat herd in Kafr El Sheikh, Egypt. The body weight of the animals was recorded at the beginning of the experiment, and the animals were housed individually and had free access to drinking water. All animals in both groups received the same diet throughout the experiment: fresh Egyptian clover (Trifolium alexandrinum) ad libitum and commercial concentrates (2.5% of live weight per head/day). The quantities of offered and refused clover and concentrates feed mixture were estimated daily for each goat throughout the experiment. Subsequently, dry matter intake (DMI) was calculated as the difference between the offered and refused feed. After 20 days, approximately 200 mL of rumen fluid sample was collected from each doe via a stomach tube before morning feeding. The rumen fluid samples were filtered through a two-layer cheesecloth. The pH of rumen samples was measured using a digital pH meter (WPA CD70, ADWA, Szeged, Hungary). Then, the samples were separated into two portions, which were frozen at −20°C for DNA extraction and to analyze rumen ammonia (NH3-N) and volatile fatty acids (VFA). Milk yield was measured individually on days 20, 35, and 45 using twice-daily hand milking, morning and evening. Approximately 100 mL of representative milk samples were collected to conduct milk chemical composition. At the end of the experiment, all animals were released back to the goat herd without euthanasia.
FE calculation
The gross FE was calculated using the equation FE = milk production (kg/d) divided by DMI (kg/d). Fat-corrected milk for 3.5% fat content was calculated using the formula: milk yield (MY) 3.5% = (0.432 + 0.1625 × % milk fat) × milk yield, kg/d (Sklan et al., 1992). Subsequently, adjusted FE was calculated using the formula: FE = 3.5% fat-corrected milk yield (kg)/dry matter intake (kg). Milk net energy (Milk NE) was calculated using the equation described in National Research Council (NRC) (2001): Milk NE (Mcal of NE_L/d) = Milk Production x (0.0929 * Fat % + 0.0563 * Protein +0.0395 * Lactose %). Finally, FE for lactation was calculated using the formula: Mcal/kg = Milk NE (Mcal /d) / DMI /d.
Chemical analyses
Diets: Fresh Egyptian clover, concentrate feed mixture, and refused Egyptian clover were analyzed (replicate for every sample”) according to the method of AOAC (1997) to measure dry matter (DM, method 930.15), crude protein (CP, method 954.01), crude fiber (CF, method 962.09), ether extract (EE, method 920.39), and ash contents (method 942.05; Table 1). Nitrogen-free extract (NFE) was estimated by the difference from the sum of the protein, fat, ash, and crude fiber content.
Table 1. Chemical composition of the experimental diets.
Chemical composition of milk: The percentage of milk protein, fat, lactose, and total solid, in addition to somatic cell count (SCC, cells/ml), were analyzed using a MilkoScan (130 A/SN. Foss Electric, Hilleroed, Denmark) using three aliquots (6 mL/aliquot) for each sample. Before the testing, milk samples were heated to 40°C and homogenized by vortexing for 20 s.
Fatty acids in milk: The total lipid content of milk was extracted using the Folch et al. (1957) method. Briefly, 2 mL of milk sample were transferred into 15 mL screw cap tubes. Then, 6 mL of a mixture containing chloroform and methanol in a ratio of 2:1 were added to the tubes, followed by vortexing for 3 min. Next, 2 mL of deionized water was added to the tubes and vortexed for 3 min. The tubes were then centrifuged for 30 min at 5000xg. The lower phase containing the extracted lipids was transferred to clean tubes and allowed to dry at room temperature. To methylate the fatty acids, sodium methoxide (2 M) was used following the method suggested by Kramer et al. (1997). The extracted lipids dissolved in 1 mL of hexane were mixed with 200 μL of sodium methoxide and vortexed. The mixture was kept for 10 min at room temperature, and then the clear top layer was transferred to gas chromatography (GC) vials. The fatty acid methyl esters (FAME) were analyzed using GC with a mass spectrometer detector (Trac 1,300, Thermo Fisher Scientific, Waltham, United States) and a TG-5MS Zebron capillary column. The FAME compounds were identified using AMDIS software1 based on their retention times matching the NIST library database.
Rumen fermentation: For volatile fatty acid (VFA) and ammonia analysis, 1 mL of rumen fluid samples were transferred into 1.5 mL Eppendorf tubes. The samples were then acidified by 200 μL of meta-phosphoric acid 25% (w/v) and stored at −20°C for later analysis. Upon thawing, samples were centrifuged at 30,000 × g (15,000 rpm, JA-17 rotor) for 20 min, and then the supernatant was used for VFA and ammonia determination. A 750 μL of supernatant was transferred to GC vials for VFA analysis by injecting 1 μL into gas chromatography (TRACE 1300, Thermo Fisher Scientific, Waltham, United States) using a capillary column (TR-FFAP 30 m × 0.53 mmI D × 0.5 μm). A standard with known concentrations of VFA was used for calibration. The other 250 μL of supernatant was used for ammonia determination colorimetrically using an ammonia assay kit (Biodiagnostic Company, Dokki, Giza, Egypt). The predicted methane was determined using the following equation: Methane yield = 316/propionate +4.4, according to Williams et al. (2019).
Rumen microbial community
DNA extraction and PCR amplification: Total microbial DNA was extracted from 500 μL of rumen sample. Briefly, the sample was centrifuged at 13,000 rpm. According to the manufacturer’s instructions, DNA was extracted from the precipitated solid material using a QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany). The quantity and quality of DNA were assessed using agarose gel electrophoresis and a Nanodrop spectrophotometer 2000 (Thermo Scientific, Massachusetts, United States). The V4 region of the 16S rRNA gene was amplified with the primers 515F and 926R using the following PCR conditions: 94°C for 3 min; 35 cycles of 94°C for 45 s, 50°C for 60 s, and 72°C for 90 s; and a final extension at 72°C for 10 min. The rumen archaeal community was studied using primers Ar915aF and Ar1386R, and the PCR amplification was conducted under the following conditions: 95°C for 5 min; 30 cycles of 95°C for 20 s, 55°C for 15 s, 72°C for 5 min, and a final extension at 72°C for 10 min. PCR amplicons were purified and sequenced using the Illumina MiSeq system.
Bioinformatics and statistical analyses
The generated paired-end sequence reads were analyzed using the DADA2 pipeline in the R platform (Callahan et al., 2016). The fastq files of sequence reads were demultiplexed, and their quality was evaluated. Then, the sequences were filtered, trimmed, and dereplicated, followed by merging read 1 and read 2 together to get denoised sequences. Generate denoised amplicon sequence variants (ASVs), and chimeric ASVs were then removed. Taxonomic assignment of ASVs was performed using the “assign taxonomy” and “addSpecies” functions, as well as the SILVA reference database (version 138). Alpha diversity metrics were measured, including observed ASVs, Chao1, Shannon diversity, and inverse Simpson diversity indices.
Additionally, the beta diversity of the bacterial and archaeal communities was determined as principal coordinate analysis (PCoA) using Bray–Curtis dissimilarity, and figures were created using the phyloseq and ggplot R packages. The differences in alpha diversity indices, relative abundances of bacterial and archaeal phyla and dominant bacterial and archaeal genera, feed intake, rumen fermentation parameters, milk yield and composition, and FE were examined using an unpaired t-test at p < 0.05. Function prediction of microbial communities associated with Shami (SH) and Zaraibi (ZA) goat breeds was conducted based on the 16S rRNA data using PICRUSt2 (Douglas et al., 2020) based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The false discovery rate (FDR) method was used for multiple-comparison correction, and p-values below 0.05 were considered significant.
Blood sampling and genotyping
For animal genotyping, blood samples were collected using vacutainer tubes containing EDTA from each animal’s jugular vein. Samples were directly transferred in an icebox to the Molecular Genetics laboratory at DRC for DNA extraction. Genomic DNA was extracted using a Puregen Core Genomic DNA extraction from blood (Qiagen®, Germany) according to the manufacturer’s protocol. The quantity and quality of extracted DNA were assessed using a Nanodrop spectrophotometer 2000 (Thermo Scientific, Massachusetts, United States). High-quality DNA samples (≥ 50 ng/μL) were used for genotyping at the Research Institute for Farm Animal Biology, Dummerstorf, Germany, using the Illumina®Inc. Goat_IGGC_65K_v2 Infinium HD array (Illumina, San Diego, CA, United States). SNP locations reported in this paper are based on the latest goat genome version of Capra hircus available from the National Center for Biotechnology Information (NCBI) database (ARS1, NCBI). The genotyping BeadChip contained 59,727 SNPs, evenly distributed throughout the entire genome. Genotype calling was performed using GenomeStudio software (Illumina Inc., San Diego, CA, United States) according to the manufacturer’s protocols. The quality control of genotyped SNPs was performed using PLINK v1.9 software (Chang et al., 2015), with the following filtering criteria: (i) SNPs showing significant deviation from Hardy–Weinberg Equilibrium (HWE) were excluded (p < 10−6); (ii) SNPs with a minor allele frequency (MAF) ≤ 0.01 were removed; (iii) markers with a genotype call rate < 99% and individuals with a call rate < 90% were filtered out. Additionally, SNPs mapped to unknown chromosomal positions or duplicate positions on the same chromosome were excluded to ensure data integrity for subsequent analyses.
Genome-wide association analysis
Genetic variance explained by markers. The single-step GBLUP implemented in the BLUPF90 family (Legarra et al., 2014) was used to estimate the SNP effects from genomic estimated breeding values (GEBVs) of genotyped animals using the postGSf90 software of the BLUPF90 package (Aguilar et al., 2010). SNP effects were calculated as: û = DZ’ [ZDZ’]−1 âg, where û is the vector of the SNP effect; D is the diagonal matrix for weighting factors of the SNP effect; Z is the matrix of genotypes, and âg is the vector of breeding values predicted for genotyped animals. The variance explained by each SNP was calculated as σ2 = û2 2p(1 − p), where û is the SNP effect described above, and p is the allele frequency of the SNP (Zhang et al., 2010). The percentage of genetic variance explained by a window segment of 5 adjacent SNPs was calculated as: (Var (ai)) / (σ2a) x100%, where ai is the genetic value of the i-th region that consists of 5 adjacent SNPs, and σ2a is the total genetic variance (Wang et al., 2014). Manhattan plots of minus log10 of SNP p-values versus chromosomal location were drawn using the qqman package in R (Turner, 2018).
Functional annotation, candidate genes, and gene enrichment analysis
For genome windows above the threshold (i.e., explained >0.1% of the total genetic variance), functional annotation of genes was obtained from BioMart at the Ensembl Genome Browser2 (Kinsella et al., 2011). Gene functions and protein domains were identified by the UniProt OMIA (Online Mendelian Inheritance in Animals) and the GeneCards databases (Safran et al., 2021). The genes that overlapped with the identified genomic interval of the candidate genome were considered candidate genes for FE in goats and were enriched using ShinyGO v. 0.77 (Ge et al., 2020) software. The analysis was based on gene ontology (GO; Ashburner et al., 2000) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Kanehisa et al., 2023) against the goat gene set ontologies. The program’s default parameters were applied, and the results were adjusted to a false discovery rate (FDR < 0.05).
Results
Feed intake. Animals in this study were fed the same diet, which consisted of a concentrate feed mixture and fresh Egyptian clover. ZA goats exhibited significantly lower body weight compared to the SH goats (p < 0.05; Table 2). Goat breed significantly affected feed intake (g/kg^0.75, metabolic body weight), with ZA goats consuming more DM, CP, CF, EE, and NFE than the SH goats (p < 0.05).
Table 2. Feed intake (g/kg 0.75), milk yield and composition, and FE parameters in SH and ZA goat breeds.
Daily milk yield (DMY), chemical composition, and FE. The estimated average DMY was higher in the ZA goats than in the SH goats. However, the difference was not significant (p > 0.05), and 3.5% of fat-corrected milk followed the same trend (Table 2). Furthermore, the percentages of milk fat, lactose, total solids (TS), solids not fat (SNF), and somatic cell count (SCC) did not exhibit significant differences between the studied goat breeds, but milk protein was higher in the ZA goats (2.55%) than in the SH goats (2.29%; p < 0.05). Based on milk yield and DMI, ZA goats exhibited greater FE expressed as gross FE and adjusted FE (Table 2; p < 0.05). Moreover, milk net energy was higher in the ZA goats without significant difference, and FE for lactation was higher in the ZA goats compared to the SH goats (p < 0.05; Table 2).
Milk fatty acids: The most abundant fatty acids identified were palmitic (C16:0), oleic (C18:1-cis), stearic (C18:0), myristic (C14:0), capric (C10:0), lauric (C12:0), and linoleic (C18:2-cis-9.12; Table 3). The proportions of certain fatty acids in milk were influenced by the breed of goat. The SH goats exhibited higher percentages of caproic (C6:0), myristic (C14:0), palmitic (C16:0), and total saturated fatty acids (∑ SFA; p < 0.05). In contrast, the ZA goats exhibited higher percentages of tridecylic (C13:0), heptadecanoic (C17:0), stearic (C18:0), behenic (C22:0), pentadecenoic (C15:1, cis-10), oleic (C18:1, cis-9), elaidic (C18:1, trans-9), rumelenic (C18:2, cis 9, trans 11), mono-unsaturated fatty acids (∑ MUFA), and polyunsaturated fatty acids (∑ PUFA; p < 0.05; Table 3).
Table 3. Milk fatty acids (%) in Shami (SH) and Zaraibi (ZA) goat breeds.
Rumen fermentation: The concentrations of acetic, propionic, and total VFA were higher in the ZA goats (p < 0.05; Table 4). Additionally, rumen ammonia, isovaleric, and predicted methane were higher in the SH goats (p < 0.05; Table 4).
Table 4. Rumen fermentation parameters in Shami (SH) and Zaraibi (ZA) goat breeds.
Microbial community
Diversity of bacterial community: The total high-quality non-chimeric reads in the SH and ZA goats were 570,677 and 666,870, respectively, and the difference in the average reads number between the goat groups was not significant (p > 0.05; Table 5). Alpha diversity was calculated as observed ASVs, Chao1, Shannon, and inverse Simpson. (Table 5) showed that the ZA goat group had numerically lower alpha diversity metrics without significant differences (p > 0.05). Principal coordinate analysis (PCoA) based on the Bray-Curtis distance was conducted to show the similarity between microbial communities (Figure 1). The results revealed that bacterial communities were clustered separately based on the goat breeds.
Table 5. Alpha diversity indexes of rumen bacteria and relative abundances (%) of bacterial phyla in the rumen of SH and ZA goat breeds.
Figure 1. Principal coordinates analysis (PCoA) of bacterial community. PCoA of rumen bacteria in lactating goats based on Bray-Curtis dissimilarity. The analysis was conducted between two goats’ breeds: red circles for Shami breed (SH), blue circles for Zaraibi breed (ZA).
Structure of bacterial community: At the phylum level, the bacterial community in the rumen of goats was affiliated with eight bacterial phyla, including Bacteroidota, Cyanobacteria, Desulfobacterota, Fibrobacterota, Firmicutes, Planctomycetota, Spirochaetota, and Verrucomicrobiota, which were shared between all rumen samples (Table 5). In addition, one bacterial phylum, Synergistota, was observed only in the SH goats. The goat breed affected the relative abundance of bacterial groups (Table 5). Bacteroidota and Firmicutes were the most abundant bacterial phyla, representing more than 98% of the bacterial community. The Bacteroidota phylum dominated the bacterial community and accounted for 85.06 and 87.54% in the SH and ZA goats, respectively. This phylum was classified mainly into the families Prevotellaceae, Rikenellaceae, Bacteroidales RF16 group, Bacteroidales p-251-o5, Bacteroidales F082, Muribaculaceae, and Bacteroidales BS11 gut group (Table 6). Except for the family Prevotellaceae, other families showed higher relative abundance in the SH goats. Family Prevotellaceae represented the most abundant family and accounted for 60.38% in the SH goats and 66.57% in the ZA goats with a significant difference (p < 0.05); this family was classified mainly into two genera, Prevotella and Prevotellaceae NK3B31 group, that were higher in the ZA goats. Family Rikenellaceae was dominated by the Rikenellaceae RC9 gut group, which was higher in the SH goats (Table 6).
Table 6. Relative abundances (%) of dominant bacterial families (F) and genera (G) in the rumen of Shami (SH) and Zaraibi (ZA) goat breeds.
Phylum Firmicutes accounted for 14.77% of the SH goats and 11.03% of the ZA goats, with a significant difference (p < 0.05; Table 5). This phylum was classified mainly into four families, Ruminococcaceae, Oscillospiraceae, Christensenellaceae, and Lachnospiraceae, which were higher in the SH goats (p < 0.05; Table 6). Family Oscillospiraceae was affiliated mainly with the genus NK4A214 group, family Christensenellaceae was affiliated mainly with the genus Christensenellaceae R-7 group, and family Lachnospiraceae was affiliated mainly with the Lachnospiraceae NK3A20 group that was higher in the ZA goats (p < 0.05) in addition to the genus Butyrivibrio that was higher in the SH goats (Table 6). Phylum Fibrobacterota and Spirochaetota were more abundant in the ZA goats compared to the SH goats (p < 0.05; Table 5).
Diversity of archaeal community: After data processing and the chimeras were removed, a total of 40,313 remained in group SH, and 61,061 reads remained in the ZA goats without significant difference (p > 0.05; Table 7). Alpha diversity, including observed ASVs, Chao1, Shannon, and Invisimpsone, was similar between the goat breeds, and the ZA goats showed numerically higher observed ASVs, Chao1, and Shannon, and the SH goats exhibited higher Invisimpsone (Table 7; p > 0.05). Principal coordinate analysis (PCoA) based on the Bray-Curtis of archaeal community across the goat breeds demonstrated that the samples were clustered based on the animal breeds (Figure 2).
Table 7. Alpha diversity indexes of rumen archaea and relative abundances (%) of archaeal genera in the rumen of Shami (SH) and Zaraibi (ZA) goat breeds.
Figure 2. Principal coordinates analysis (PCoA) of archaeal community. The analysis was conducted between two goat breeds: red circles for the Shami breed (ASH) and blue circles for the Zaraibi breed (ZA).
Structure of archaeal community: The taxonomic analysis of the archaeal community revealed that the ASVs were assigned to two phyla, Thermoplasmatota, which was further classified to the family Methanomethylophilaceae, and phylum Euryarchaeota was affiliated with the family Methanobacteriaceae (Table 7). Goat breed affected the archaeal community, and the family Methanomethylophilaceae had two genera, Candidatus Methanomethylophilus, that accounted for 18.45% in group Shami goats (ASH) and 64.60% in Zaraibi goats (AZA) group with a significant difference (p < 0.05), and unclassified Methanomethylophilaceae that was higher in the SH goats compared to AZA. The family Methanobacteriaceae was further categorized into two genera. Methanobrevibacter constituted 61.44% of the archaeal community in the SH goats and 23.04% in the AZA group, showing a significant difference (p < 0.05). Additionally, the genus Methanosphaera was more prevalent in the SH goats compared to AZA (Table 7).
Function prediction of microbial community: The metabolic pathways of rumen microbial communities in Shami (SH) and Zaraibi (ZA) goats were performed using PICRUSt2. Principal components analysis (PCA) of PICRUSt2 function prediction was generated using the relative abundance of metabolic pathways (Figure 3), which revealed that samples were clustered based on animal breeds. The samples of the SH goats revealed significantly higher (q < 0.05) relative abundances of metabolic pathways related to heme biosynthesis-II, L-arginine degradation (ARGORNPROST-PWY), P562-PWY, polyamine biosynthesis I (POLYAMSYN-PWY), ARG + POLYAMINE-SYN, L-histidine degradation I (HISDEG-PWY), biotin biosynthesis-PWY, the glyoxylate cycle (GLYOXYLATE-BYPASS), pyruvate fermentation to acetone (PWY-6588), the glyoxylate cycle and TCA cycle (GLYCOLYSIS-TCA-GLYOX-BYPASS), anaerobic gondoate biosynthesis (PWY-7663), coenzyme B biosynthesis (TCA-GLYOX-BYPASS), and the tricarboxylic acid (TCA) cycle IV (2-oxoglutarate decarboxylase; P105-PWY; Figure 4). Additionally, the samples of group ZA revealed significantly higher (q < 0.05) relative abundances of metabolic pathways related to ascorbate and aldarate metabolism, carbohydrate metabolism, nitrogen metabolism, sulfur metabolism, lactose and galactose degradation (LACTOSECAT-PWY), amino acid metabolism (BRANCHED-CHAIN-AA-SYN-PWY), L-ornithine biosynthesis I (GLUTORN-PWY), glycolysis III (ANAGLYCOLYSIS-PWY), the sucrose degradation pathway (PWY-5384), nucleoside and nucleotide degradation (PWY-5532), GDP-d-glycero-alpha-d-manno-heptose biosynthesis (PWY-6478), and the tricarboxylic acid cycle citric acid cycle (PWY-7254; Figure 4).
Figure 3. Principal components analysis of PICRUSt2 functional prediction of microbial communities in the rumen of Shami and Zaraibi goat breeds. Orange squares refer to ZA samples, and blue circles refer to SH’s samples.
Figure 4. Effect of animal species on the relative abundances of metabolic pathways of rumen microbial communities in the rumen of SH and ZA goat breeds. Orange pars refer to ZA’s samples, and blue pars refer to SH’s samples.
Genetic variance explained by markers: A significant genome window was considered to contribute to the genetic variance of FE traits if it accounted for more than 0.1%. A total of 26 genome windows, each explaining >0.1% of the total genetic variance explained by all windows, with a total of 2.32% of genetic variance (Figure 5). These windows were located on chromosomes 1 (6, 110, and 132 Mb), 2 (9 and 116 Mb), 3 (24 Mb), 7 (88 and 94 Mb), 8 (60 Mb), 10 (76 Mb), 11 (94 and 95 Mb), 12 (47 Mb), 13 (53 Mb), 15 (31 Mb), 18 (21, 23, 25, and 36 Mb), 21 (25 and 41 Mb), 22 (49 Mb), 24 (33 Mb), 25 (1 and 22 Mb), and 26 (3 Mb; Table 8; Supplementary file 1).
Figure 5. Manhattan plots of genome-wide association results for gross FE (A), adjusted FE (B), milk net efficiency (C), and FE for lactation (D) in Egyptian goats. Each dot represents a genome window. The scale of the y-axis represents the significance of the percentage of genetic variance explained by the genome window, and chromosomes are displayed on the x-axis.
Table 8. Percentages of variance explained by genome windows and the annotated genes for FE traits in goats.
Gene enrichment analysis: The gene enrichment analysis identified 26 significant (q < 0.05) biological pathways using a list of the identified candidate genes in the current study (Supplementary Excel file S1). The identified candidate genes were significantly enriched in biological pathways related to the regulation of the transport process of many cellular components such as ions, lipids, organic substances, amino acids, and the general transporter activities mechanism (Figure 6).
Figure 6. Gene ontology (GO) enrichment analysis of biological pathways for the list of candidate genes resulted from genome-wide association analysis for FE in Egyptian goats. False discovery rate (FDR) < 0.05.
Discussion
FE represents an animal’s capacity to convert feed components into products such as meat and milk through rumen microbial fermentation. This efficiency directly impacts the profitability of livestock production enterprises (Fischer et al., 2022). Therefore, investigating the relationship between the rumen microbiome, animal genotype, and FE is crucial for genetic improvement programs in livestock species comprising goats. In the current study, goats’ diets consisted of a mixture of feed concentrates and fresh clover, which is a high-quality roughage. Estimates of milk yield and feed intake values were in agreement with the values reported by Ghoneem and El-Tanany (2023). On the other hand, milk yield in the current study (Table 2) was higher than those obtained by Alsheikh (2013) and Kholif et al. (2020) for Shami and Zaribi goat breeds. To the best of our knowledge, no studies compared the FE, feed intake, and milk yield between the SH and ZA goats under the same husbandry conditions. However, the studies conducted by Ghoneem and El-Tanany (2023) and Kholif et al. (2020) reported estimates of FE in the SH and ZA goats similar to our study.
ZA goats achieved higher FE (0.96 kg milk/kg DM) than the SH goats (0.69; p < 0.05), likely due to their superior feed utilization and adaptation to the arid conditions of southern Egypt, their region of origin (Alsheikh, 2013). Similarly, Knights and Garcia (1997) reviewed Zaraibi goats as producing higher milk than Shami goats.
Variations in feed intake, milk yield, and composition (Tables 2, 3), and FE between breeds were reported in different ruminant species (Knights and Garcia, 1997; Paz et al., 2016; Currò et al., 2019; Gustavsson et al., 2014; Niero et al., 2023). Bharanidharan et al. (2018) reported that high-efficiency dairy cows showed higher feed intake and higher milk yield. Our results showed that ZA goats exhibited higher levels of milk protein and unsaturated fatty acids (∑ PUFA), which is consistent with the higher abundance of Prevotella observed. The breed-specific differences in milk fatty acid profiles observed in this study align with previous findings in goats and sheep (Idamokoro et al., 2019). This could be attributed to variations between animal species and breeds in metabolic pathways of lipids in mammary glands and rumen microorganisms (Conte et al., 2022). Therefore, it is necessary to understand the differences in milk fat profiles in different livestock breeds.
Zhang et al. (2023) and Sasson et al. (2017) explained that host genome and rumen microbiome impact milk protein, which could be altered due to a higher abundance of rumen bacteria: Prevotella, Bacteroidales, Clostridiales, Flavefaciens, and Ruminococcus. Additionally, the proportion of polyunsaturated fatty acids (PUFA) is affected mainly by rumen microorganisms through rumen biohydrogenation (Shingfield et al., 2010; Huws et al., 2011). In this regard, Hu et al. (2020) identified a positive correlation between the higher Prevotella and higher unsaturated fatty acids (UFA) in yaks, as Prevotella play a vital role in energy harvesting in the rumen, and it might provide precursors for UFA synthesis, which demonstrates the higher UFA in ZA compared to the SH goats in this study. A similar conclusion was obtained by Huws et al. (2011), who studied the bacteria that play a predominant role in ruminal biohydrogenation.
Rumen ecosystem and FE
The rumen bacterial community in goats was dominated by the phyla Bacteroidota, Firmicutes, Spirochaetota, and Fibrobacterota (Table 5), aligning with findings from previous studies on goats (Giger-Reverdin et al., 2020; Luo et al., 2022). The Firmicutes phylum was higher in the SH goats with lower FE, which was also reported in low-efficiency beef cattle by Brooke et al. (2019).
Host genetics significantly influence the rumen microbiome, with approximately 35% of microbial taxa being heritable (Li et al., 2019). This implies that the difference in the microbial communities between the two goat breeds in the current study is expected. In this context, our PCoA analysis (Figure 1) showed that the bacterial communities in SH and ZA breeds were distinct, which agrees with the findings of Paz et al. (2016) in Holstein and Jersey breeds. Similar results were obtained by Mani et al. (2021) in the rumen of Damara and Meatmaster sheep breeds in South Africa. In the current study, ZA showed numerically low alpha bacterial diversity compared to the SH goats (p > 0.05; Table 5), which is similar to other findings in dairy cows (Shabat et al., 2016).
The variation in the composition of rumen microbiota results in changes in rumen fermentation efficiency, wherever the ZA goats had lower ammonia than the SH goats (24.51 vs. 37.10 mmol/L, p < 0.05). Furthermore, the ZA goats showed higher total VFA than the SH goats (82.76 vs. 72.83 mM, p < 0.05) as well as higher acetic and propionic (Table 4; p < 0.05). VFAs represent the main energy source of the host animal (Paz et al., 2016). Nitrogen and carbohydrate metabolism are the main microbial functions in the rumen (Xue et al., 2022). Higher VFA production in ZA was associated with higher fiber-degrading bacteria and higher carbohydrate metabolism pathways (Table 4; Figure 4), which agrees with Tapio et al. (2017) and Xue et al. (2022), indicating that Zaraibi goats have higher efficiency in carbohydrate metabolism. This finding is supported by higher relative abundances of Fibrobacters, Prevotellaceae, and Lachnospiraceae in the rumen of the ZA goats. These bacterial groups have essential roles in complex carbohydrates and protein metabolism (Thoetkiattikul et al., 2013; Paz et al., 2016; Neumann et al., 2018). Xue et al. (2022) reported that pathways related to carbohydrate metabolism could be used as prediction markers to differentiate efficient and inefficient rumen microbiomes or animals, which helps the future selection of high-efficiency animals. These findings were supported by Wang et al. (2023), who reported higher carbohydrate metabolism pathways in goats with higher growth performance. On the other hand, the ZA goats showed higher relative abundances of branched-chain amino acids (BCAA) pathways. The increase in the relative abundance of branched-chain amino acids (BCAA) pathways could refer to a higher concentration of BCAA in the rumen. Zhang H. L. et al. (2013) reported that BCAA stimulates fiber-degrading microorganisms and VFA production.
Furthermore, the increment in acetic acid production was accompanied by an increment in the PWY-7254 pathway (tricarboxylic acid cycle citric acid cycle), which is involved in the synthesis of precursors of acetic acid (De Sales-Millán et al., 2024). Group ZA showed lower rumen ammonia and higher relative abundances of pathways related to nitrogen metabolism, which could indicate the efficient use of nitrogen in microbial protein synthesis and the decline in nitrogen loss (Firkins et al., 2007). Higher synthesis of microbial protein increases the protein supply to the host animals Zhang H. L. et al. (2013). Furthermore, the ZA goats showed higher sulfur metabolism, which is linked to lower methane production (Wu et al., 2021), as sulfate-reducing bacteria are able to compete with methanogens for H2 in the rumen, which inhibits the methanogenesis (Zhao and Zhao, 2022).
Genus Prevotella is a major player in ruminal metabolism and predominated the rumen microbiome in several ruminant species (Betancur-Murillo et al., 2022). A higher relative abundance of genus Prevotella was associated with the highly efficient breed (i.e., ZA; Tables 2, 6). This finding is supported by Brooke et al. (2019), who reported that the Prevotella is a microbial marker identifying cattle with higher FE. Genus Prevotella includes members that utilize a wide range of substrates, such as hemicellulose and protein (Matsui et al., 2000). Myer et al. (2015) found a higher relative abundance of Prevotella in the gut of steers with high daily gain and intake. De Vadder et al. (2016) indicated that the genus Prevotella produces succinate, which is the precursor to propionate. The propionate synthesis uses hydrogen molecules, reducing its availability for methane production by methanogenic archaea (Betancur-Murillo et al., 2022). Thus, Prevotella has the potential to be used as an anti-methanogenic agent.
The members of the family Lachnospiraceae have fibrolytic and cellulolytic activities and were associated with higher milk production in Holstein cows (Thoetkiattikul et al., 2013; Paz et al., 2016). In this study, the candidate genus Lachnospiraceae NK3A20 was higher in the ZA goats; this genus was reported to produce butyric acid that promotes the development of the rumen (Huang et al., 2021). Boeckaert et al. (2009) indicated that the Lachnospiraceae NK3A20 group plays an important role in the efficiency of energy utilization. Genus Succiniclasticum is one of the dominant bacteria in the rumen of heifers (Liu et al., 2017), and it converts the succinate to propionate (van Gylswyk, 1995). Shabat et al. (2016) explained that cattle with higher propionate-producing rumen bacteria, such as Succiniclasticum and Prevotella, could achieve efficient feed utilization.
The rumen of the SH goat breed had a significantly higher proportion of the genus Christensenellaceae R-7 group, Mogibacterium. Christensenellaceae R-7 group was associated with FE, rumen digestion and absorption of nutrients, fiber digestion, and protein metabolism, and this genus produces acetic and butyric acids (Bach et al., 2019; Huang et al., 2021). Genus Mogibacterium was higher in the low FE group (i.e., SH). A previous study reported that Mogibacterium was associated with higher methane-emitting cattle (Wallace et al., 2015). Furthermore, Mogibacterium cannot degrade carbohydrates for energy, and it had a higher proportion in lower weight gain in steers (McLoughlin et al., 2020).
Furthermore, the rumen of the SH goats showed a numerically higher proportion than that of the family Ruminococcaceae. This family was linked to lower FE in beef cattle (Brooke et al., 2019), which supports the current findings. The function prediction of the rumen microbiome of the SH goats is enriched with hem and biotin synthesis pathways and L-arginine degradation. Hem is an important component of hemoglobin and a cofactor for many enzymes; therefore, it has an important role in oxygen transportation and various physiological processes (Yang et al., 2023). Arginine plays a key role in urea cycle regulation and protein synthesis, as ammonia is the main product of arginine degradation (Chacher et al., 2012). Biotin vitamin acts as a cofactor responsible for carbon dioxide transfer in carboxylases and was enriched in the rumen of cows with high milk yield and protein (Xue et al., 2022).
Rumen archaea: Methane emission from farm animals is one of the main determiners of FE. The composition of the rumen methanogens was associated with the animal breed, methane emission, and FE (Bharanidharan et al., 2018; McLoughlin et al., 2023). Archaea represent 1–2% of the microbial community and are considered the sole producer of methane in the rumen (Rabee et al., 2022; McLoughlin et al., 2023). Methane represents a 2–12% loss in FE to the host (Johnson and Johnson, 1995). PCoA analysis demonstrated (Figure 2) that the diversity of the archaeal community was influenced by the breed. The previous studies indicated that the abundance of methanogenic archaea in the rumen was associated with methane emissions and FE (Wallace et al., 2015; Zhou et al., 2010).
Methanobrevibacter accounted for 61.44% of the archaeal community in the SH goats, significantly more than in the ZA goats (Table 7; p < 0.05), which correlates with the higher predicted methane emissions in the SH goats. Meanwhile, ZA goats demonstrated a higher relative abundance of Candidatus Methanomethylophilus (64.60%) than the SH goats (18.45%). Candidatus Methanomethylophilus is an H2-dependent methylotrophic methanogen that derives its energy from the metabolism of methanol and methylamine (Borrel et al., 2012). This genus was associated with improved FE, higher daily gain, and lower methane emissions in different sheep breeds and Charolais steers (McLoughlin et al., 2023), which supports the current findings. Methanobrevibacter is the main methane producer in the rumen (Tapio et al., 2017). It uses hydrogen molecules, besides other substrates, such as acetate and formate, to produce methane (Jeyanathan et al., 2011). Bharanidharan et al. (2018) reported that Hanwoo steers (low efficiency) showed higher methane emissions and a higher abundance of Methanobrevibacter compared to Holstein steers (high efficiency), which supports our findings. The same conclusion was reported by Zhou et al. (2009) on beef steers. Previous studies reported that Methanobrevibacter uses acetate in methane production, leading to higher CH4 energy loss in low-FE animals (Nkrumah et al., 2006; Hegarty et al., 2007). Propionic acid production suppresses methane production by consuming the hydrogen from the rumen ecosystem. Furthermore, sulfur metabolism could reduce the availability of H2 for methane production (Zhao and Zhao, 2022). These findings explain the relationship between higher efficiency in the ZA goats and lower Methanobrevibacter and higher Prevotella as the main propionic producer in the rumen (Tapio et al., 2017; Bharanidharan et al., 2018). Thus, microbial communities in goats varied according to breed. Higher FE in the ZA goats could be attributed to higher relative abundances of lignocellulolytic bacteria that affect feed utilization positively. Moreover, a lower relative abundance of major methane-producing archaea might decline the methane energy emission (Difford et al., 2018; Gonzalez-Recio et al., 2018).
Variance explained by markers. FE is a quantitative trait that may be affected by a few genes with large or modest effects. Otherwise, it is affected by multiple genes with small effects (Zhang et al., 2020). In our study, we identified 26 genome windows. Each explained >0.1% of the total genetic variance explained by all windows, with a total of 1.89, 1.98, 2.39, and 1.78% of genetic variance for gross FE (G), adjusted FE (AFE), milk net energy (MNE), and FE for lactation (FEL), respectively. This suggests that the trait is likely controlled by multiple SNPs with small effects (Zhang et al., 2020). Additionally, most of the identified genes seem to have pleiotropic effects and should be considered in selective breeding programs (Sallam, 2023). Alternatively, there is a high genetic correlation between the studied traits that were undertaken to represent FE (e.g., G, AFE, MNE, and FEL). Thus, they shared common genome windows with the same positional candidate genes.
Our GWAS analysis identified multiple genomic regions that may be significantly associated with FE in Egyptian goats. The 132 megabases (Mb) on Capra hircus autosome (CHI) 1 was the most interesting genomic region that contributed to FE with approximately 0.26, 0.28, 0.27, and 0.3% of the total genetic variance explained for GFE, AFE, MNE, and FEL, respectively. This genomic region harbored the propionyl-CoA carboxylase subunit beta (PCCB) gene, which was involved in pathways related to the metabolism of water-soluble vitamins. Additionally, it is one of the two subunits of the biotin-dependent propionyl-CoA carboxylase (PCC), a mitochondrial enzyme involved in the catabolism of odd-chain fatty acids, branched-chain amino acids isoleucine, threonine, methionine, valine, and other metabolites (Wongkittichote et al., 2017).
An interesting genomic region was also located at 95 Mb on CHI11, which contributed to the genetic variance of G, AFE, MNE, and FEL with approximately 0.2, 0.16, 0.17, and 0.16%, respectively. This genomic region harbored the Golgin A1 (GOLGA1) gene, which is involved in vesicular trafficking at the Golgi apparatus level and endosome-to-Golgi trafficking (Wong and Munro, 2014), and the GTPase activating protein and VPS9 domains 1 (GAPVD1) gene, which enables GTPase activating protein binding activity and is involved in the regulation of protein transport. Additionally, the RABGAP1 (RAB GTPase activating protein 1) and DENND1A (DENN domain-containing 1A) genes are involved in pathways related to vesicle-mediated transport and trafficking regulation (Yoshimura et al., 2010; Samarelli et al., 2020). Thus, they were putative candidate genes for FE.
Two potential candidate genes were identified on CHI2 (116 Mb): the Metaxin 2 (MTX2) gene, which is involved in the transport of proteins into the mitochondrion (Del Castillo et al., 2021), and the Alkylglycerone Phosphate Synthase (AGPS) gene, which encodes a protein that catalyzes the second step of ether lipid biosynthesis (Chornyi et al., 2023). Similarly, the TMEM241 (transmembrane protein 241) was identified on CHI24 (33 Mb) and was involved in carbohydrate transport and transmembrane transport (Rodríguez et al., 2016). On CHI21 (40 Mb), the SCFD1 (Sec1 Family Domain Containing 1) gene was observed. It is involved in the regulation of protein transport (Hou et al., 2017), in addition to the Adaptor Related Protein Complex 4 Subunit Sigma 1 (AP4S1) gene, which was included in pathways related to transporter activity and obsolete protein transporter activity (Dell'Angelica et al., 1999). On CHI12 (35 Mb), the COMM Domain Containing 6 (COMMD6) gene, which is included in pathways related to protein metabolism, and TBC1D4 (TBC1 Domain Family Member 4), which is considered an insulin transport.
Similarly, the Enoyl-CoA Delta Isomerase 1 (ECI1) gene, involved in fatty acid metabolism pathways, was identified on CHI25, which spans 1 Mb. Important genes associated with fat and obesity metabolism were located on CHI18 (23 Mb), including the Iroquois Homeobox 3 (IRX3) gene, which acts as a regulator of energy metabolism and is linked to pathways involving the FTO obesity variant mechanism. Additionally, the Alpha-Ketoglutarate Dependent Dioxygenase (FTO) gene, which is known as the fat mass and obesity-associated gene, was also identified (Frayling et al., 2007).
On CHI22, which encompasses 49 Mb, the CACNA2D2 (Calcium Voltage-Gated Channel Auxiliary Subunit Alpha2delta 2) gene was found. This gene is involved in the regulation of calcium channels that control the entry of calcium ions into the cell (Bartuzi et al., 2013). Also identified on this chromosome is the TMEM115 gene, which a role in retrograde transport of proteins from the Golgi to the endoplasmic reticulum (Ong et al., 2014).
CHI8 spans 60 Mb and includes the Glucosylceramidase beta 2 (GBA2) gene, which plays a role in carbohydrate transport and metabolism (Woeste and Wachten, 2017). On CHI15, which spans 31 Mb, the STARD10 (StAR-related lipid transfer domain containing 10) gene is mapped and is predicted to be involved in lipid transport (Floris et al., 2019). The Acyl-CoA dehydrogenase family member 8 (ACAD8) gene, encoding a member of the acyl-CoA dehydrogenase family, catalyzes the dehydrogenation of acyl-CoA derivatives in fatty acids and amino acid metabolism (Zhuang et al., 2022). Additionally, the Galactosidase Beta 1 Like 3 (GLB1L3) gene, also located within this region, is predicted to be involved in carbohydrate metabolic processes (Nicoli et al., 2021).
The 36 megabases of CHI18 harbored several putative candidate genes: (1) the NRN1L (Neuritin 1-like) gene, which is involved in protein metabolism pathways (Zhang W. J. et al., 2013), and (2) the Lecithin-Cholesterol acyltransferase (LCAT) gene, which encodes an extracellular cholesterol esterifying enzyme essential for cholesterol transport (Kris-Etherton et al., 2023).
Lastly, multiple members of the SLC (solute carrier family) gene family transporters (e.g., SLC12A4 12, SLC7A6, SLC17A9, SLCO4A1, SLC38A3, and SLC35G2) were identified in most of the significant genome windows. SLC transporters participate in many important physiological functions, such as nutrient supply, and are highly expressed in different organs, including the kidney, brain, liver, gut, and heart (Zhang et al., 2019). Moreover, our gene enrichment analysis showed that the candidate genes identified in the current study were enriched in biological pathways related to the regulation of the transport process of many cellular components such as ions, lipids, organic substances, amino acids, and the general transporter activities mechanism. These findings support our hypothesis of the significant contribution of the animal genotype in FE in Egyptian goats. Thus, the findings of this study indicated that highly efficient animals have a higher representation of animal gene markers and microbial markers related to the metabolism of carbohydrates, protein, and lipids. Therefore, this information could be a useful aid in selective breeding programs that aim at improving FE in goats.
Conclusion
The ZA goats showed higher FE compared to the SH goats. These differences in FE may result from variations in rumen microbiota, with ZA goats demonstrating enhanced fiber utilization and reduced methane emissions. Furthermore, association analysis revealed 26 genome windows containing several putative candidate genes that significantly contribute to FE traits in Egyptian goats. These genes are involved in pathways related to cellular component metabolism and transport, making them promising candidates for enhancing FE in Egyptian goats. In addition, variations were observed between the goat breeds in terms of milk protein content and certain milk fatty acids.
Future studies should leverage genetic markers and microbial insights to improve goat breeding, enhance FE, and optimize milk quality in arid regions while also aiming to also reduce methane emissions.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found at: https://www.ncbi.nlm.nih.gov/sra/PRJNA1047779. The host genome datasets generated and/or analyzed during the current study are available for the corresponding breeds in the OSF repository using the following link: https://osf.io/7q6j8/ and can be downloaded upon request from the corresponding author.
Ethics statement
The animal study was approved by the Animal Breeding Ethics Committee at the Desert Research Center (DRC) in Egypt (Project ID: 43213) with approval reference number AB/NO2022. Research Ethics Review Committee, Faculty of Agriculture, University of Alexandria, Egypt, approved the experiment (Reference: Alex. Agri. 082305307). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
AR: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. IA-s: Writing – review & editing, Data curation, Formal analysis, Methodology, Validation, Writing – original draft. AY: Conceptualization, Formal analysis, Methodology, Project administration, Writing – original draft, Writing – review & editing, Data curation. ML: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. ME-S: Conceptualization, Data curation, Formal analysis, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. EE: Data curation, Formal analysis, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing, Investigation. AS: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Validation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by grants from the Science, Technology & Innovation Funding Authority (STDF), and the Egyptian Ministry for Scientific Research within the activities of the ongoing research project with reference number “43213.”
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1492742/full#supplementary-material
Footnotes
References
Aguilar, I., Misztal, I., Johnson, D. L., Legarra, A., Tsuruta, S., and Lawlor, T. J. (2010). Hot topic: a unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J. Dairy Sci. 93, 743–752. doi: 10.3168/jds.2009-2730
Almasri, O., Abou-Bakr, S., and Ibrahim, M. (2023). Lactation curve and milk production traits of syrian damascus goats. Egyptian J. Animal Production 60, 7–16. doi: 10.21608/ejap.2023.177189.1051
Alsheikh, S. M. (2013). Influence of age and live body weight on daily milk yield of Zaraibi and Shami goats in Sinai. Egypt. Annals of Agricul. Sci. 58, 1–3. doi: 10.1016/j.aoas.2013.01.001
AOAC (1997). Association of Official Analytical Chemists. International Official Methods of Analysis. 16th Edition, Arlington: AOAC.
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature genetics 25, 25–29. doi: 10.1038/75556
Bach, A., López-García, A., González-Recio, O., Elcoso, G., Fàbregas, F., Chaucheyras-Durand, F., et al. (2019). Changes in the rumen and colon microbiota and effects of live yeast dietary supplementation during the transition from the dry period to lactation of dairy cows. J. Dairy Sci. 102, 6180–6198. doi: 10.3168/jds.2018-16105
Barría, A., Benzie, J. A. H., Houston, R. D., De Koning, D. J., and de Verdal, H. (2021). Genomic selection and genome-wide association study for feed-efficiency traits in a farmed Nile Tilapia (Oreochromis niloticus) population. Front. Genet. 12:737906. doi: 10.3389/fgene.2021.737906
Bartuzi, P., Hofker, M. H., and van de Sluis, B. (2013). Tuning NF-κB activity: a touch of COMMD proteins. Biochim. Biophys. Acta 1832, 2315–2321. doi: 10.1016/j.bbadis.2013.09.014
Betancur-Murillo, C. L., Aguilar-Marín, S. B., and Jovel, J. (2022). Prevotella: a key player in ruminal metabolism. Microorganisms 11:1. doi: 10.3390/microorganisms11010001
Bharanidharan, R., Arokiyaraj, S., Kim, E. B., Lee, C. H., Woo, Y. W., Na, Y., et al. (2018). Ruminal methane emissions, metabolic, and microbial profile of Holstein steers fed forage and concentrate, separately or as a total mixed ration. PLoS One 13:e0202446. doi: 10.1371/journal.pone.0202446
Boeckaert, C., Morgavi, D. P., Jouany, J. P., Maignien, L., Boon, N., and Fievez, V. (2009). Role of the protozoan Isotricha prostoma, liquid-, and solid-associated bacteria in rumen biohydrogenation of linoleic acid. Animal: Int. J. Animal Biosci. 3, 961–971. doi: 10.1017/S1751731109004285
Borrel, G., Harris, H. M., Tottey, W., Mihajlovski, A., Parisot, N., Peyretaillade, E., et al. (2012). Genome sequence of "Candidatus Methanomethylophilus alvus" Mx1201, a methanogenic archaeon from the human gut belonging to a seventh order of methanogens. J. Bacteriol. 194, 6944–6945. doi: 10.1128/JB.01867-12
Brooke, C. G., Najafi, N., Dykier, K. C., and Hess, M. (2019). Prevotella copri, a potential indicator for high feed efficiency in western steers. Animal Sci. J. = Nihon chikusan Gakkaiho 90, 696–701. doi: 10.1111/asj.13197
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Chacher, B., Wang, D. M., Liu, H. Y., and Liu, J. X. (2012). Degradation of L-arginine and N-carbamoyl glutamate and their effect on rumen fermentation in vitro. Ital. J. Anim. Sci. 23, 989–1001. doi: 10.4081/ijas.2012.e68
Chang, C. C., Chow, C. C., Tellier, L. C., Vattikuti, S., Purcell, S. M., and Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience 4:7. doi: 10.1186/s13742-015-0047-8
Chornyi, S., Ofman, R., Koster, J., and Waterham, H. R. (2023). The origin of long-chain fatty acids required for de novo ether lipid/plasmalogen synthesis. J. Lipid Res. 64:100364. doi: 10.1016/j.jlr.2023.100364
Conte, G., Palombo, V., Serra, A., Correddu, F., D'Andrea, M., Macciotta, N. P. P., et al. (2022). Study of the fatty acid profile of Milk in different sheep breeds: evaluation by multivariate factorial analysis. Animals: Open Access J. MDPI 12:722. doi: 10.3390/ani12060722
Currò, S., Manuelian, C. L., De Marchi, M., Claps, S., Rufrano, D., and Neglia, G. (2019). Effects of breed and stage of lactation on Milk fatty acid composition of Italian goat breeds. Animals: Open Access J. MDPI 9:764. doi: 10.3390/ani9100764
De Sales-Millán, A., Reyes-Ferreira, P., Aguirre-Garrido, J. F., Corral-Guillé, I., Barrientos-Ríos, R., and Velázquez-Aragón, J. A. (2024). Comprehensive analysis of gut microbiota composition and functional metabolism in children with autism Spectrum disorder and Neurotypical children: implications for sex-based differences and metabolic dysregulation. Int. J. Mol. Sci. 25:6701. doi: 10.3390/ijms25126701
De Vadder, F., Kovatcheva-Datchary, P., Zitoun, C., Duchampt, A., Bäckhed, F., and Mithieux, G. (2016). Microbiota-Produced Succinate Improves Glucose Homeostasis via Intestinal Gluconeogenesis. Cell Metab. 24, 151–157. doi: 10.1016/j.cmet.2016.06.013
Del Castillo, U., Norkett, R., Lu, W., Serpinskaya, A., and Gelfand, V. I. (2021). Ataxin-2 is essential for cytoskeletal dynamics and neurodevelopment in Drosophila. iScience 25:103536. doi: 10.1016/j.isci.2021.103536
Dell'Angelica, E. C., Mullins, C., and Bonifacino, J. S. (1999). AP-4, a novel protein complex related to clathrin adaptors. J. Biol. Chem. 274, 7278–7285. doi: 10.1074/jbc.274.11.7278
Desire, S., Mucha, S., Coffey, M., Mrode, R., Broadbent, J., and Conington, J. (2018). “Deriving genomic breeding values for feed intake and body weight in dairy goats.” In Proceedings of the world congress on genetics applied to livestock production, volume species - caprine: 818.
Difford, G. F., Plichta, D. R., Løvendahl, P., Lassen, J., Noel, S. J., Højberg, O., et al. (2018). Host genetics and the rumen microbiome jointly associate with methane emissions in dairy cows. PLoS Genet. 14:e1007580. doi: 10.1371/journal.pgen.1007580
Douglas, G. M., Maffei, V. J., Zaneveld, J. R., Yurgel, S. N., Brown, J. R., Taylor, C. M., et al. (2020). PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 38, 685–688. doi: 10.1038/s41587-020-0548-6
Eid, J. I., Teleb, D. F., Mohamed, S. A., and El-Ghor, A. A. (2020). DGAT1 polymorphism in Egyptian Zaraibi goat breed and their association with milk yield and composition. J. Basic and Applied Zoology 81, 1–7. doi: 10.1186/s41936-020-00176-w
Firkins, J. L., and Yu, Z. (2015). RUMINANT NUTRITION SYMPOSIUM: how to use data on the rumen microbiome to improve our understanding of ruminant nutrition. J. Anim. Sci. 93, 1450–1470. doi: 10.2527/jas.2014-8754
Firkins, J. L., Yu, Z., and Morrison, M. (2007). Ruminal nitrogen metabolism: perspectives for integration of microbiology and nutrition for dairy. J. Dairy Sci. 90, E1–E16. doi: 10.3168/jds.2006-518
Fischer, A., Dai, X., and Kalscheur, K. F. (2022). Feed efficiency of lactating Holstein cows is repeatable within diet but less reproducible when changing dietary starch and forage concentrations. Animal 16:100599. doi: 10.1016/j.animal.2022.100599
Floris, A., Luo, J., Frank, J., Zhou, J., Orrù, S., Biancolella, M., et al. (2019). Star-related lipid transfer protein 10 (STARD10): a novel key player in alcohol-induced breast cancer progression. J. Experimental Clin. Cancer Res.: CR 38:4. doi: 10.1186/s13046-018-1013-y
Folch, J., Lees, M., and Sloane Stanley, G. H. (1957). A simple method for the isolation and purification of total lipides from animal tissues. J. Biol. Chem. 226, 497–509. doi: 10.1016/S0021-9258(18)64849-5
Frayling, T. M., Timpson, N. J., Weedon, M. N., Zeggini, E., and Freathy, R. M. (2007). A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science (New York, N.Y.) 316, 889–894. doi: 10.1126/science.1141634
Ge, S. X., Jung, D., and Yao, R. (2020). ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics (Oxford, England) 36, 2628–2629. doi: 10.1093/bioinformatics/btz931
Ghoneem, W. M. A., and El-Tanany, R. R. A. (2023). Impact of natural betaine supplementation on rumen fermentation and productive performance of lactating Damascus goats. Trop. Anim. Health Prod. 55:123. doi: 10.1007/s11250-023-03524-4
Giger-Reverdin, S., Domange, C., Broudiscou, L. P., Sauvant, D., and Berthelot, V. (2020). Rumen function in goats, an example of adaptive capacity. J. Dairy Res. 87, 45–51. doi: 10.1017/S0022029920000060
Gonzalez-Recio, O., Zubiria, I., García-Rodríguez, A., Hurtado, A., and Atxaerandio, R. (2018). Short communication: signs of host genetic regulation in the microbiome composition in 2 dairy breeds: Holstein and Brown Swiss. J. Dairy Sci. 101, 2285–2292. doi: 10.3168/jds.2017-13179
Güney, O., Torun, O., Özuyanık, O., and Darcan, N. (2006). Milk production, reproductive and growth performances of Damascus goats under northern Cyprus conditions. Small Rumin. Res. 65, 176–179. doi: 10.1016/j.smallrumres.2005.07.026
Gustavsson, F., Buitenhuis, A. J., Johansson, M., Bertelsen, H. P., Glantz, M., Poulsen, N. A., et al. (2014). Effects of breed and casein genetic variants on protein profile in milk from Swedish red, Danish Holstein, and Danish Jersey cows. J. Dairy Sci. 97, 3866–3877. doi: 10.3168/jds.2013-7312
Hegarty, R. S., Goopy, J. P., Herd, R. M., and McCorkell, B. (2007). Cattle selected for lower residual feed intake have reduced daily methane production. J. Anim. Sci. 85, 1479–1486. doi: 10.2527/jas.2006-236
Hou, N., Yang, Y., Scott, I. C., and Lou, X. (2017). The sec domain protein Scfd1 facilitates trafficking of ECM components during chondrogenesis. Dev. Biol. 421, 8–15. doi: 10.1016/j.ydbio.2016.11.010
Hu, R., Zou, H., Wang, H., Wang, Z., Wang, X., Ma, J., et al. (2020). Dietary energy levels affect rumen bacterial populations that influence the intramuscular fat fatty acids of fattening yaks (Bos grunniens). Animals: Open Access J. MDPI 10:1474. doi: 10.3390/ani10091474
Huang, C., Ge, F., Yao, X., Guo, X., Bao, P., Ma, X., et al. (2021). Microbiome and metabolomics reveal the effects of different feeding systems on the growth and ruminal development of yaks. Front. Microbiol. 12:682989. doi: 10.3389/fmicb.2021.682989
Huws, S. A., Kim, E. J., Lee, M. R., Scott, M. B., Tweed, J. K., Pinloche, E., et al. (2011). As yet uncultured bacteria phylogenetically classified as Prevotella, Lachnospiraceae incertae sedis and unclassified Bacteroidales, Clostridiales and Ruminococcaceae may play a predominant role in ruminal biohydrogenation. Environ. Microbiol. 13, 1500–1512. doi: 10.1111/j.1462-2920.2011.02452.x
Idamokoro, E. M., Muchenje, V., Afolayan, A. J., and Hugo, A. (2019). Comparative fatty-acid profile and atherogenicity index of milk from free grazing Nguni, Boer and non-descript goats in South Africa. Pastoralism 9, 1–8. doi: 10.1186/s13570-019-0138-8
Jeyanathan, J., Kirs, M., Ronimus, R. S., Hoskin, S. O., and Janssen, P. H. (2011). Methanogen community structure in the rumens of farmed sheep, cattle and red deer fed different diets. FEMS Microbiol. Ecol. 76, 311–326. doi: 10.1111/j.1574-6941.2011.01056.x
Johnson, K. A., and Johnson, D. E. (1995). Methane emissions from cattle. J. Anim. Sci. 73, 2483–2492. doi: 10.2527/1995.7382483x
Kanehisa, M., Furumichi, M., Sato, Y., Kawashima, M., and Ishiguro-Watanabe, M. (2023). KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 51, D587–D592. doi: 10.1093/nar/gkac963
Khanal, P., Johnson, J., Gouveia, G., Ross, P., and Deeb, N. (2023). Genomic evaluation of feed efficiency in US Holstein heifers. J. Dairy Sci. 106, 6986–6994. doi: 10.3168/jds.2023-23258
Kholif, A. E., Gouda, G. A., and Hamdon, H. A. (2020). Performance and Milk composition of Nubian goats as affected by increasing level of Nannochloropsis oculata microalgae. Animals: Open Access J. MDPI 10:2453. doi: 10.3390/ani10122453
Kinsella, R. J., Kähäri, A., Haider, S., Zamora, J., Proctor, G., Spudich, G., et al. (2011). Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database: J. Biological Databases and Curation 2011:bar030. doi: 10.1093/database/bar030
Knights, M., and Garcia, G. W. (1997). The status and characteristics of the goat (Capra hircus) and its potential role as a significant milk producer in the tropics: a review. Small Rumin. Res. 26, 203–215. doi: 10.1016/S0921-4488(96)00977-7
Kramer, J. K., Fellner, V., Dugan, M. E., Sauer, F. D., Mossoba, M. M., and Yurawecz, M. P. (1997). Evaluating acid and base catalysts in the methylation of milk and rumen fatty acids with special emphasis on conjugated dienes and total trans fatty acids. Lipids 32, 1219–1228. doi: 10.1007/s11745-997-0156-3
Kris-Etherton, P. M., Sanders, L., Lawler, O., Riley, T., and Maki, K. (2023). Hyperlipidemia. In encyclopedia of human nutrition. Academic Press. 4, 361–379.
Legarra, A., Christensen, O. F., Aguilar, I., and Misztal, I. (2014). Single step, a general approach for genomic selection. Livest. Sci. 166, 54–65. doi: 10.1016/j.livsci.2014.04.029
Li, F., and Guan, L. L. (2017). Metatranscriptomic profiling reveals linkages between the active rumen microbiome and feed efficiency in beef cattle. Appl. Environ. Microbiol. 83, e00061–e00017. doi: 10.1128/AEM.00061-17
Li, F., Li, C., Chen, Y., Liu, J., Zhang, C., Irving, B., et al. (2019). Host genetics influence the rumen microbiota and heritable rumen microbial features associate with feed efficiency in cattle. Microbiome 7:92. doi: 10.1186/s40168-019-0699-1
Liu, C., Li, X. H., Chen, Y. X., Cheng, Z. H., Duan, Q. H., Meng, Q. H., et al. (2017). Age-related response of rumen microbiota to mineral salt and effects of their interactions on enteric methane emissions in cattle. Microb. Ecol. 73, 590–601. doi: 10.1007/s00248-016-0888-4
Løvendahl, P., Difford, G. F., Li, B., Chagunda, M. G. G., Huhtanen, P., Lidauer, M. H., et al. (2018). Review: selecting for improved feed efficiency and reduced methane emissions in dairy cattle. Animal: Int. J. Animal Biosci. 12, s336–s349. doi: 10.1017/S1751731118002276
Luo, T., Li, Y., Zhang, W., Liu, J., and Shi, H. (2022). Rumen and fecal microbiota profiles associated with immunity of young and adult goats. Front. Immunol. 13:978402. doi: 10.3389/fimmu.2022.978402
Mani, S., Aiyegoro, O. A., and Adeleke, M. A. (2021). Characterization of rumen microbiota of two sheep breeds supplemented with direct-fed lactic acid Bacteria. Front. Vet. Sci. 7:570074. doi: 10.3389/fvets.2020.570074
Matsui, H., Ogata, K., Tajima, K., Nakamura, M., Nagamine, T., Aminov, R. I., et al. (2000). Phenotypic characterization of polysaccharidases produced by four Prevotella type strains. Curr. Microbiol. 41, 45–49. doi: 10.1007/s002840010089
McLoughlin, S., Spillane, C., Campion, F. P., Claffey, N., Sosa, C. C., McNicholas, Y., et al. (2023). Breed and ruminal fraction effects on bacterial and archaeal community composition in sheep. Sci. Rep. 13:3336. doi: 10.1038/s41598-023-28909-1
McLoughlin, S., Spillane, C., Claffey, N., Smith, P. E., O'Rourke, T., Diskin, M. G., et al. (2020). Rumen microbiome composition is altered in sheep divergent in feed efficiency. Front. Microbiol. 11:1981. doi: 10.3389/fmicb.2020.01981
Moreira, G. C. M., Poleti, M. D., Pértille, F., Boschiero, C., Cesar, A. S. M., Godoy, T. F., et al. (2019). Unraveling genomic associations with feed efficiency and body weight traits in chickens through an integrative approach. BMC Genet. 20:83. doi: 10.1186/s12863-019-0783-3
Myer, P. R., Smith, T. P., Wells, J. E., Kuehn, L. A., and Freetly, H. C. (2015). Rumen microbiome from steers differing in feed efficiency. PLoS One 10:e0129174. doi: 10.1371/journal.pone.0129174
National Research Council (NRC) (2001). Nutrient Requirements of Dairy Cattle. 7th Edition, Washington DC, USA: National Academic Press.
Neumann, A. P., Weimer, P. J., and Suen, G. (2018). A global analysis of gene expression in Fibrobacter succinogenes S85 grown on cellulose and soluble sugars at different growth rates. Biotechnol. Biofuels 11:295. doi: 10.1186/s13068-018-1290-x
Nicoli, E. R., Annunziata, I., d'Azzo, A., Platt, F. M., Tifft, C. J., and Stepien, K. M. (2021). GM1 Gangliosidosis-A Mini-Review. Front. Genet. 12:734878. doi: 10.3389/fgene.2021.734878
Niero, G., Thomas, S. A., Mouratidou, K., Visentin, G., De Marchi, M., Penasa, M., et al. (2023). Lactoferrin concentration in bovine milk: validation of radial immunodiffusion technique, sources of variation, and association to udder health status. Ital. J. Anim. Sci. 22, 230–238. doi: 10.1080/1828051X.2023.2180440
Nkrumah, J. D., Okine, E. K., Mathison, G. W., Schmid, K., Li, C., Basarab, J. A., et al. (2006). Relationships of feedlot feed efficiency, performance, and feeding behavior with metabolic rate, methane production, and energy partitioning in beef cattle. J. Anim. Sci. 84, 145–153. doi: 10.2527/2006.841145x
Ong, Y. S., Tran, T. H., Gounko, N. V., and Hong, W. (2014). TMEM115 is an integral membrane protein of the Golgi complex involved in retrograde transport. J. Cell Sci. 127, 2825–2839. doi: 10.1242/jcs.136754
Paz, H. A., Anderson, C. L., Muller, M. J., Kononoff, P. J., and Fernando, S. C. (2016). Rumen bacterial community composition in Holstein and Jersey cows is different under same dietary condition and is not affected by sampling method. Front. Microbiol. 7:1206. doi: 10.3389/fmicb.2016.01206
Rabee, A. E., Kewan, K. Z., El Shaer, H. M., Lamara, M., and Sabra, E. A. (2022). Effect of olive and date palm by-products on rumen methanogenic community in Barki sheep. AIMS Microbiol. 8, 26–41. doi: 10.3934/microbiol.2022003
Rodríguez, A., Gonzalez, L., Ko, A., Alvarez, M., Miao, Z., Bhagat, Y., et al. (2016). Molecular characterization of the lipid genome-wide association study signal on chromosome 18q11.2 implicates HNF4A-mediated regulation of the TMEM241 gene. Arterioscler. Thromb. Vasc. Biol. 36, 1350–1355. doi: 10.1161/ATVBAHA.116.307182
Safran, M., Rosen, N., Twik, M., BarShir, R., Stein, T. I., Dahary, D., et al. (2021). “The GeneCards suite,” in Practical Guide to Life Science Databases. eds. I. Abugessaisa and T. Kasukawa (21-01 Gateway East, Singapore: Springer), 27–56. doi: 10.1007/978-981-16-5812-9_2
Sallam, A. M. (2023). Effect of genetic polymorphisms in LALBA and prolactin genes on milk traits in Barki ewes. Small Rumin. Res. 226:107041. doi: 10.1016/j.smallrumres.2023.107041
Samarelli, A. V., Ziegler, T., Meves, A., Fässler, R., and Böttcher, R. T. (2020). Rabgap1 promotes recycling of active β1 integrins to support effective cell migration. J. Cell Sci. 133:jcs243683. doi: 10.1242/jcs.243683
Sasson, G., Kruger Ben-Shabat, S., Seroussi, E., Doron-Faigenboim, A., Shterzer, N., Yaacoby, S., et al. (2017). Heritable bovine rumen Bacteria are phylogenetically related and correlated with the Cow's capacity to harvest energy from its feed. MBio 8, e00703–e00717. doi: 10.1128/mBio.00703-17
Shabat, S. K., Sasson, G., Doron-Faigenboim, A., Durman, T., Yaacoby, S., Berg Miller, M. E., et al. (2016). Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 10, 2958–2972. doi: 10.1038/ismej.2016.62
Shingfield, K. J., Bernard, L., Leroux, C., and Chilliard, Y. (2010). Role of trans fatty acids in the nutritional regulation of mammary lipogenesis in ruminants. Animal: Int. J. Animal Biosci. 4, 1140–1166. doi: 10.1017/S1751731110000510
Singh, M. K., Dige, M. S., Pourouchottamane, R., Kumar, A., and Gowane, G. R. (2022). Influences of maternal factors on the estimate of genetic parameters for goat feed efficiency traits. Trop. Anim. Health Prod. 54:376. doi: 10.1007/s11250-022-03355-9
Sklan, D., Ashkenazi, R., Braun, A., Devorin, A., and Tabori, K. (1992). Fatty acids, calcium soaps of fatty acids, and cottonseeds fed to high yielding cows. J. Dairy Sci. 75, 2463–2472. doi: 10.3168/jds.S0022-0302(92)78008-4
Tapio, I., Snelling, T. J., Strozzi, F., and Wallace, R. J. (2017). The ruminal microbiome associated with methane emissions from ruminant livestock. J. Animal Sci. Biotechnol. 8, 7–11. doi: 10.1186/s40104-017-0141-0
Tatar, A. M., Tuncer, S. S., and Sireli, H. D. (2021). Comparison of yield characteristics of Damascus and Kilis goats in dry climatic conditions. Austral J. Vet. Sci. 51, 61–66. doi: 10.4067/S0719-81322019000200061
Taussat, S., Boussaha, M., Ramayo-Caldas, Y., Martin, P., Venot, E., Cantalapiedra-Hijar, G., et al. (2020). Gene networks for three feed efficiency criteria reveal shared and specific biological processes. Genetics, selection, evolution: GSE 52:67. doi: 10.1186/s12711-020-00585-z
Thoetkiattikul, H., Mhuantong, W., Laothanachareon, T., Tangphatsornruang, S., Pattarajinda, V., Eurwilaichitr, L., et al. (2013). Comparative analysis of microbial profiles in cow rumen fed with different dietary fiber by tagged 16S rRNA gene pyrosequencing. Curr. Microbiol. 67, 130–137. doi: 10.1007/s00284-013-0336-3
Turner, S. D. (2018). Qqman: an R package for visualizing GWAS results using QQ and Manhattan plots. J. Open Source Software 3:731. doi: 10.21105/joss.00731
Uffelmann, E., Huang, Q. Q., Munung, N. S., de Vries, J., Okada, Y., Martin, A. R., et al. (2021). Genome-wide association studies. Nature Rev. Methods Primers 1:59. doi: 10.1038/s43586-021-00056-9
van Gylswyk, N. O. (1995). Succiniclasticum ruminis gen. Nov., sp. nov., a ruminal bacterium converting succinate to propionate as the sole energy-yielding mechanism. Int. J. Syst. Bacteriol. 45, 297–300. doi: 10.1099/00207713-45-2-297
Wallace, R. J., Rooke, J. A., McKain, N., Duthie, C. A., Hyslop, J. J., Ross, D. W., et al. (2015). The rumen microbial metagenome associated with high methane production in cattle. BMC Genomics 16:839. doi: 10.1186/s12864-015-2032-0
Wang, D., Chen, L., Tang, G., Yu, J., Chen, J., Li, Z., et al. (2023). Multi-omics revealed the long-term effect of ruminal keystone bacteria and the microbial metabolome on lactation performance in adult dairy goats. Microbiome 11:215. doi: 10.1186/s40168-023-01652-5
Wang, H., Huang, J., Xu, X., and Xiao, Y. (2014). Damage attack on complex networks. Physica A: Statistical Mechanics and its Applications 408, 134–148. doi: 10.1016/j.physa.2014.04.001
Weigel, K. A. (2017). Genomic selection of dairy cattle: a review of methods, strategies, and impact. J. Animal Breeding and Genomics 1, 1–15. doi: 10.12972/jabng.20170001
Williams, S. R. O., Hannah, M. C., Jacobs, J. L., Wales, W. J., and Moate, P. J. (2019). Volatile fatty acids in ruminal fluid can be used to predict methane yield of dairy cows. Animals: Open Access J. MDPI 9:1006. doi: 10.3390/ani9121006
Woeste, M. A., and Wachten, D. (2017). The enigmatic role of GBA2 in controlling locomotor function. Front. Mol. Neurosci. 10:386. doi: 10.3389/fnmol.2017.00386
Wong, M., and Munro, S. (2014). Membrane trafficking. The specificity of vesicle traffic to the Golgi is encoded in the golgin coiled-coil proteins. Science (New York, N.Y.) 346:1256898. doi: 10.1126/science.1256898
Wongkittichote, P., Ah Mew, N., and Chapman, K. A. (2017). Propionyl-CoA carboxylase - a review. Mol. Genet. Metab. 122, 145–152. doi: 10.1016/j.ymgme.2017.10.002
Wu, H., Li, Y., Meng, Q., and Zhou, Z. (2021). Effect of high sulfur diet on rumen fermentation, microflora, and epithelial barrier function in steers. Animals: Open Access J. MDPI 11:2545. doi: 10.3390/ani11092545
Xue, M. Y., Xie, Y. Y., Zhong, Y., Ma, X. J., Sun, H. Z., and Liu, J. X. (2022). Integrated meta-omics reveals new ruminal microbial features associated with feed efficiency in dairy cattle. Microbiome 10:32. doi: 10.1186/s40168-022-01228-9
Yang, Q., Zhao, J., Zheng, Y., Chen, T., and Wang, Z. (2023). Microbial synthesis of Heme b: biosynthetic pathways, current strategies, detection, and future prospects. Molecules (Basel, Switzerland) 28:3633. doi: 10.3390/molecules28083633
Yoshimura, S., Gerondopoulos, A., Linford, A., Rigden, D. J., and Barr, F. A. (2010). Family-wide characterization of the DENN domain Rab GDP-GTP exchange factors. J. Cell Biol. 191, 367–381. doi: 10.1083/jcb.201008051
Zhang, H. L., Chen, Y., Xu, X. L., and Yang, Y. X. (2013). Effects of branched-chain amino acids on in vitro ruminal fermentation of wheat straw. Asian Australas. J. Anim. Sci. 26, 523–528. doi: 10.5713/ajas.2012.12539
Zhang, Z., Liu, J., Ding, X., Bijma, P., de Koning, D. J., and Zhang, Q. (2010). Best linear unbiased prediction of genomic breeding values using a trait-specific marker-derived relationship matrix. PLoS One 5:e12648. doi: 10.1371/journal.pone.0012648
Zhang, C., Wang, M., Liu, H., Jiang, X., Chen, X., Liu, T., et al. (2023). Multi-omics reveals that the host-microbiome metabolism crosstalk of differential rumen bacterial enterotypes can regulate the milk protein synthesis of dairy cows. J. Animal Sci. Biotechnol. 14, 63–18. doi: 10.1186/s40104-023-00862-z
Zhang, F., Wang, Y., Mukiibi, R., Chen, L., Vinsky, M., Plastow, G., et al. (2020). Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: I: feed efficiency and component traits. BMC Genomics 21:36. doi: 10.1186/s12864-019-6362-1
Zhang, Y., Zhang, Y., Sun, K., Meng, Z., and Chen, L. (2019). The SLC transporter in nutrient and metabolic sensing, regulation, and drug development. J. Mol. Cell Biol. 11, 1–13. doi: 10.1093/jmcb/mjy052
Zhang, W. J., Zhang, W. G., Zhang, P. Y., Cao, X. M., He, A. L., Chen, Y. X., et al. (2013). The expression and functional characterization associated with cell apoptosis and proteomic analysis of the novel gene MLAA-34 in U937 cells. Oncol. Rep. 29, 491–506. doi: 10.3892/or.2012.2129
Zhao, Y., and Zhao, G. (2022). Decreasing ruminal methane production through enhancing the sulfate reduction pathway. Animal Nutrit. (Zhongguo xu mu shou yi xue hui) 9, 320–326. doi: 10.1016/j.aninu.2022.01.006
Zhou, M., Hernandez-Sanabria, E., and Guan, L. L. (2009). Assessment of the microbial ecology of ruminal methanogens in cattle with different feed efficiencies. Appl. Environ. Microbiol. 75, 6524–6533. doi: 10.1128/AEM.02815-08
Zhou, M., Hernandez-Sanabria, E., and Guan, L. L. (2010). Characterization of variation in rumen methanogenic communities under different dietary and host feed efficiency conditions, as determined by PCR-denaturing gradient gel electrophoresis analysis. Appl. Environ. Microbiol. 76, 3776–3786. doi: 10.1128/AEM.00010-10
Zhuang, D. Y., Ding, S. X., Wang, F., Yang, X. C., Pan, X. L., Bao, Y. W., et al. (2022). Identification of six novel variants of ACAD8 in Isobutyryl-CoA dehydrogenase deficiency with increased C4 carnitine using tandem mass spectrometry and NGS sequencing. Front. Genet. 12:791869. doi: 10.3389/fgene.2021.791869
Keywords: goat breeds, feed efficiency, caprine genome, rumen fermentation, bacteria and archaea
Citation: Rabee AE, Abou-souliman I, Yousif AI, Lamara M, El-Sherbieny MA, Elwakeel EA and Sallam AM (2025) Variations in rumen microbiota and host genome impacted feed efficiency in goat breeds. Front. Microbiol. 16:1492742. doi: 10.3389/fmicb.2025.1492742
Edited by:
Anusorn Cherdthong, Khon Kaen University, ThailandReviewed by:
Maharach Matra, Mahasarakham University, ThailandRittikeard Prachumchai, Rajamangala University of Technology Thanyaburi, Thailand
Copyright © 2025 Rabee, Abou-souliman, Yousif, Lamara, El-Sherbieny, Elwakeel and Sallam. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alaa Emara Rabee, UmFiZWVfYV9tQHlhaG9vLmNvbQ==; YWxhYS5yYWJlZUBkcmMuZ292LmVn; Mebarek Lamara, bWViYXJlay5sYW1hcmFAdXFhdC5jYQ==