Rumen microbiotas are known to influence the fat deposition (FD) in sheep, but controversy over causality remains unresolved. Here, we performed microbiome-wide association studies (MWAS), microbiome genome-wide association analysis (mbGWAS) and bidirectional mendelian randomization (MR) analyses on 1,150 sheep with genotype data from whole-genome resequencing, 16S rRNA sequencing and multilevel FD-traits data. We quantified the proportion of individual variation in FD-traits explained by host genetics, rumen microbiota, and their interaction effects. We identified 32 rumen microbiota biomarkers including Bifidobacterium that were associated with FD-traits (P_adj <0.05). Further, utilizing five MR methods, we identified eight causal associations between marker genera and FD-traits (P_adj <0.05), including Butyrivibrio, Olsenella, p-2534-18B5 gut group, Prevotellaceae UCG-003, and Pseudobutyrivibrio causing forward causal effects on FD, and changes in Flexilinea and Suttonella induced by FD. To our knowledge, this is the inaugural attempt to employ MR in sheep to investigate the causal relationships between gastrointestinal microbiota and complex phenotypes, underscoring the potential for developing interventions related to adipose deposition in sheep from the perspective of the rumen microbiome.

Mycoplasmal pneumonia of sheep (MPS), caused by Mesomycoplasma (Mycoplasma) ovipneumoniae, profoundly impacts ovine productivity and survival. Although gut–lung microbiota interactions are increasingly recognized in respiratory diseases, whether similar crosstalk occurs between the lung and rumen microbiota in MPS-affected sheep remains unknown. To investigate alterations in the lung and rumen microbiota of sheep with MPS, the crosstalk between these microbial communities, and their impacts on growth phenotypes. From a cohort of 414 naturally infected six-month-old male Hu sheep, we selected 10 individuals with severe pulmonary pathology and 10 healthy controls for detailed phenotypic and microbiome analyses. Assessment of 359 phenotypic traits revealed that MPS significantly impairs feed efficiency and growth rate (p < 0.05). Through 16S rRNA gene sequencing, we found that MPS significantly altered the pulmonary microbiota community structure (p < 0.01), with a noticeable impact on the rumen microbiota composition (p = 0.059). Succinivibrionaceae_UCG-001 was significantly depleted in both the rumen and lungs of diseased sheep (p < 0.05) and strongly associated with reduced average daily feed intake (p < 0.05). In addition, pulmonary Pasteurella and ruminal Succinivibrionaceae_UCG-002 were significantly enriched in MPS-affected sheep, showed a strong positive correlation (p < 0.05), and were both negatively associated with feed efficiency (p < 0.05). Notably, Pasteurella multocida subsp. gallicida may act as a keystone species influencing feed efficiency. These findings point to a previously unrecognized rumen-lung microbial axis that may modulate host productivity in sheep affected by MPS. This work provides new insights into the pathogenesis of MPS and offers potential targets for therapeutic intervention and management.

Population stratification based on gut microbiota composition has revealed several enterotypes in humans and animals, providing valuable tools for studying the gut microbiota landscape, which is crucial for animal health and production. However, knowledge about rumen enterotype identification in sheep, its influencing factors, and its association with growth performance and host genetics remains limited. Here, we investigated host genetic effects and phenotypic landscapes of rumen bacterial enterotypes in a large sheep population. Ruminal contents from 1150 healthy sheep were analyzed using 16S rRNA gene sequencing and genus-level clustering, complemented by extensive phenotypic data covering 47 traits spanning growth, feed efficiency, meat yield, and ruminal fermentation, along with whole-genome resequencing data. We identified two distinct enterotypes: Enterotype 1 (E1), a mixture of multiple genera, and Enterotype 2 (E2), dominated by Prevotella. E2 sheep exhibit superior growth and meat production performance, but lower feed efficiency and increased fat deposition. Two-part beta-regression models and co-occurrence network analyses revealed the extensive impact of enterotypes on microbial community structure, with E1 displaying a higher frequency of unique bacterial interactions. The estimated heritability of the enterotype was 0.47, and a GWAS identified five key genetic markers associated with rumen enterotype, localized to two candidate genes: CHODL and ENPP6. These markers significantly influence 58 ruminal bacterial genera, including key taxa and driving genus. Overall, our data provide new insights into sheep rumen-enterotype characteristics, contributing to a better understanding of microbial interactions that are crucial for improving ruminant growth performance.