Taxonomic and predicted functional signatures reveal linkages between the rumen microbiota and feed efficiency in dairy cattle raised in tropical areas

Ruminants digest plant biomass more efficiently than monogastric animals due to their symbiotic relationship with a complex microbiota residing in the rumen environment. What remains unclear is the relationship between the rumen microbial taxonomic and functional composition and feed efficiency (FE), especially in crossbred dairy cattle (Holstein x Gyr) raised under tropical conditions. In this study, we selected twenty-two F1 Holstein x Gyr heifers and grouped them according to their residual feed intake (RFI) ranking, high efficiency (HE) (n = 11) and low efficiency (LE) (n = 11), to investigate the effect of FE on the rumen microbial taxa and their functions. Rumen fluids were collected using a stomach tube apparatus and analyzed using amplicon sequencing targeting the 16S (bacteria and archaea) and 18S (protozoa) rRNA genes. Alpha-diversity and beta-diversity analysis revealed no significant difference in the rumen microbiota between the HE and LE animals. Multivariate analysis (sPLS-DA) showed a clear separation of two clusters in bacterial taxonomic profiles related to each FE group, but in archaeal and protozoal profiles, the clusters overlapped. The sPLS-DA also revealed a clear separation in functional profiles for bacteria, archaea, and protozoa between the HE and LE animals. Microbial taxa were differently related to HE (e.g., Howardella and Shuttleworthia) and LE animals (e.g., Eremoplastron and Methanobrevibacter), and predicted functions were significatively different for each FE group (e.g., K03395-signaling and cellular process was strongly related to HE animals, and K13643-genetic information processing was related to LE animals). This study demonstrates that differences in the rumen microbiome relative to FE ranking are not directly observed from diversity indices (Faith's Phylogenetic Diversity, Pielou's Evenness, Shannon's diversity, weighted UniFrac distance, Jaccard index, and Bray-Curtis dissimilarity), but from targeted identification of specific taxa and microbial functions characterizing each FE group. These results shed light on the role of rumen microbial taxonomic and functional profiles in crossbred Holstein × Gyr dairy cattle raised in tropical conditions, creating the possibility of using the microbial signature of the HE group as a biological tool for the development of biomarkers that improve FE in ruminants.

Genome-Wide Analyses Reveal the Genetic Architecture and Candidate Genes of Indicine, Taurine, Synthetic Crossbreds, and Locally Adapted Cattle in Brazil

Cattle population history, breeding systems, and geographic subdivision may be reflected in runs of homozygosity (ROH), effective population size (N_e), and linkage disequilibrium (LD) patterns. Thus, the assessment of this information has become essential to the implementation of genomic selection on purebred and crossbred cattle breeding programs. In this way, we assessed the genotype of 19 cattle breeds raised in Brazil belonging to taurine, indicine, synthetic crossbreds, and Iberian-derived locally adapted ancestries to evaluate the overall LD decay patterns, N_e, ROH, and breed composition. We were able to obtain a general overview of the genomic architecture of cattle breeds currently raised in Brazil and other tropical countries. We found that, among the evaluated breeds, different marker densities should be used to improve the genomic prediction accuracy and power of genome-wide association studies. Breeds showing low N_e values indicate a recent inbreeding, also reflected by the occurrence of longer ROH, which demand special attention in the matting schemes to avoid extensive inbreeding. Candidate genes (e.g., ABCA7, PENK, SPP1, IFNAR1, IFNAR2, SPEF2, PRLR, LRRTM1, and LRRTM4) located in the identified ROH islands were evaluated, highlighting biological processes involved with milk production, behavior, rusticity, and fertility. Furthermore, we were successful in obtaining the breed composition regarding the taurine and indicine composition using single-nucleotide polymorphism (SNP) data. Our results were able to observe in detail the genomic backgrounds that are present in each breed and allowed to better understand the various contributions of ancestor breeds to the modern breed composition to the Brazilian cattle.

Lipopolysaccharide triggers different transcriptional signatures in taurine and indicine cattle macrophages: Reactive oxygen species and potential outcomes to the development of immune response to infections

Macrophages are classified upon activation as classical activated M1 and M2 anti-inflammatory regulatory populations. This macrophage polarization is well characterized in humans and mice, but M1/M2 profile in cattle has been far less explored. Bos primigenius taurus (taurine) and Bos primigenius indicus (indicine) cattle display contrasting levels of resistance to infection and parasitic diseases such as C57BL/6J and Balb/c murine experimental models of parasite infection outcomes based on genetic background. Thus, we investigated the differential gene expression profile of unstimulated and LPS stimulated monocyte-derived macrophages (MDMs) from Holstein (taurine) and Gir (indicine) breeds using RNA sequencing methodology. For unstimulated MDMs, the contrast between Holstein and Gir breeds identified 163 Differentially Expressed Genes (DEGs) highlighting the higher expression of C-C chemokine receptor type five (CCR5) and BOLA-DQ genes in Gir animals. LPS-stimulated MDMs from Gir and Holstein animals displayed 1,257 DEGs enriched for cell adhesion and inflammatory responses. Gir MDMs cells displayed a higher expression of M1 related genes like Nitric Oxide Synthase 2 (NOS2), Toll like receptor 4 (TLR4), Nuclear factor NF-kappa-B 2 (NFKB2) in addition to higher levels of transcripts for proinflammatory cytokines, chemokines, complement factors and the acute phase protein Serum Amyloid A (SAA). We also showed that gene expression of inflammatory M1 population markers, complement and SAA genes was higher in Gir in buffy coat peripheral cells in addition to nitric oxide concentration in MDMs supernatant and animal serum. Co-expression analyses revealed that Holstein and Gir animals showed different transcriptional signatures in the MDMs response to LPS that impact on cell cycle regulation, leukocyte migration and extracellular matrix organization biological processes. Overall, the results suggest that Gir animals show a natural propensity to generate a more pronounced M1 inflammatory response than Holstein, which might account for a faster immune response favouring resistance to many infection diseases.