Metatranscriptomic Profiling Reveals the Effect of Breed on Active Rumen Eukaryotic Composition in Beef Cattle With Varied Feed Efficiency

Abstract

Exploring the compositional characteristics of rumen eukaryotic community can expand our understanding of their role in rumen function and feed efficiency. In this study, we applied metatranscriptomics to characterize the active rumen eukaryotic community (protozoa and fungi) in beef cattle (n = 48) of three breeds [Angus (AN), Charolais (CH), and Kinsella Composite (KC)] and with divergent residual feed intake (RFI). The composition of active rumen eukaryotic microbiota was evaluated based on enriched 18S rRNAs from the metatranscriptomic datasets. At the phylum level, a total of four protozoal taxa (Ciliophora, Parabasalia, unclassified SAR, and unclassified Alveolata), six fungal taxa (Neocallimastigomycota, Basidiomycota, unclassified Fungi, Mucoromycota, Ascomycota, and Chytridiomycota), and one sister group of kingdom Fungi (unclassified Opisthokonta) were detected with relative abundances higher than 0.01% and in at least 50% of animals within each breed. Among these, Ciliophora, Parabasalia, unclassified Opisthokonta, and Neocallimastigomycota were the top four active eukaryotic phyla. At the genus level, a total of 8 ciliated protozoa, 5 flagellated protozoa, 5 anaerobic fungi, and 10 aerobic fungi taxa were detected, with unclassified Trichostomatia, Tetratrichomonas, unclassified Neocallimastigaceae, and Pleurotus being the most predominant taxa of ciliated protozoa, flagellated protozoa, anaerobic fungi, and aerobic fungi, respectively. Differential abundance analysis revealed that breed had a significant effect on the phylogenetic lineages of rumen eukaryotes, and seven fungal taxa were more abundant (linear discriminant analysis score > 2 with P < 0.05) in the rumen of KC steers than in the rumen of AN and CH steers. Although principal coordinate analysis (PCoA) revealed that the ruminal active eukaryotic profiles were not distinguishable between high- and low-RFI groups, the diversity indices, including Faith’s phylogenetic diversity (PD), observed operational taxonomic units (OTUs), and Shannon index of rumen eukaryotes were higher in low-RFI steers than those in high-RFI steers. Meanwhile, the abundance of genus Entodinium and the kingdom Fungi was higher in low-RFI steers than that in high-RFI steers. This information on active rumen eukaryotic microbiota and identified differential abundance of taxa between high- and low-RFI animals suggests the possibility of improving feed efficiency through altering rumen eukaryotic microbiota.