Establishment and transcriptomic analyses of a cattle rumen epithelial primary cells (REPC) culture by bulk and single-cell RNA sequencing to elucidate interactions of butyrate and rumen development

An integrated microbiome- and metabolome-genome-wide association study reveals the role of heritable ruminal microbial carbohydrate metabolism in lactation performance in Holstein dairy cows

BACKGROUND

Despite the growing number of studies investigating the connection between host genetics and the rumen microbiota, there remains a dearth of systematic research exploring the composition, function, and metabolic traits of highly heritable rumen microbiota influenced by host genetics. Furthermore, the impact of these highly heritable subsets on lactation performance in cows remains unknown. To address this gap, we collected and analyzed whole-genome resequencing data, rumen metagenomes, rumen metabolomes and short-chain fatty acids (SCFAs) content, and lactation performance phenotypes from a cohort of 304 dairy cows.

RESULTS

The results indicated that the proportions of highly heritable subsets (h2 ≥ 0.2) of the rumen microbial composition (55%), function (39% KEGG and 28% CAZy), and metabolites (18%) decreased sequentially. Moreover, the highly heritable microbes can increase energy-corrected milk (ECM) production by reducing the rumen acetate/propionate ratio, according to the structural equation model (SEM) analysis (CFI = 0.898). Furthermore, the highly heritable enzymes involved in the SCFA synthesis metabolic pathway can promote the synthesis of propionate and inhibit the acetate synthesis. Next, the same significant SNP variants were used to integrate information from genome-wide association studies (GWASs), microbiome-GWASs, metabolome-GWASs, and microbiome-wide association studies (mWASs). The identified single nucleotide polymorphisms (SNPs) of rs43470227 and rs43472732 on SLC30A9 (Zn2+ transport) (P < 0.05/nSNPs) can affect the abundance of rumen microbes such as Prevotella_sp., Prevotella_sp._E15-22, Prevotella_sp._E13-27, which have the oligosaccharide-degradation enzymes genes, including the GH10, GH13, GH43, GH95, and GH115 families. The identified SNPs of chr25:11,177 on 5s_rRNA (small ribosomal RNA) (P < 0.05/nSNPs) were linked to ECM, the abundance alteration of Pseudobutyrivibrio_sp. (a genus that was also showed to be linked to the ECM production via the mWASs analysis), GH24 (lysozyme), and 9,10,13-TriHOME (linoleic acid metabolism). Moreover, ECM, and the abundances of Pseudobutyrivibrio sp., GH24, and 9,10,13-TRIHOME were significantly greater in the GG genotype than in the AG genotype at chr25:11,177 (P < 0.05). By further the SEM analysis, GH24 was positively correlated with Pseudobutyrivibrio sp., which was positively correlated with 9,10,13-triHOME and subsequently positively correlated with ECM (CFI = 0.942).

CONCLUSION

Our comprehensive study revealed the distinct heritability patterns of rumen microbial composition, function, and metabolism. Additionally, we shed light on the influence of host SNP variants on the rumen microbes with carbohydrate metabolism and their subsequent effects on lactation performance. Collectively, these findings offer compelling evidence for the host-microbe interactions, wherein cows actively modulate their rumen microbiota through SNP variants to regulate their own lactation performance. Video Abstract.

Genome-Wide Association Study of Conformation Traits in Brazilian Holstein Cattle

The linear conformation of animals exerts an influence on health, reproduction, production, and welfare, in addition to longevity, which directly affects the profitability of milk-producing farms. The objectives of this study were (1) to perform genome-wide association studies (GWASs) of conformation traits, namely the Rump, Feet and Legs, Mammary System, Dairy Strength, and Final Classification traits, and (2) to identify genes and related pathways involved in physiological processes associated with conformation traits in Brazilian Holstein cattle. Phenotypic and genotypic data from 2339 Holstein animals distributed across the states of Rio Grande do Sul, Paraná, São Paulo, and Minas Gerais were used. The genotypic data were obtained with a 100 K SNP marker panel. The single-step genome-wide association study (ssGWAS) method was employed in the analyses. Genes close to a significant SNP were identified in an interval of 100 kb up- and downstream using the Ensembl database available in the BioMart tool. The DAVID database was used to identify the main metabolic pathways and the STRING program was employed to create the gene regulatory network. In total, 36 significant SNPs were found on 15 chromosomes; 27 of these SNPs were linked to genes that may influence the traits studied. Fourteen genes most closely related to the studied traits were identified, as well as four genes that showed interactions in important metabolic pathways such as myogenesis, adipogenesis, and angiogenesis. Among the total genes, four were associated with myogenesis (TMOD2, TMOD3, CCND2, and CTBP2), three with angiogenesis (FGF23, FGF1, and SCG3), and four with adipogenesis and body size and development (C5H12orf4, CCND2, EMILIN1, and FGF6). These results contribute to a better understanding of the biological mechanisms underlying phenotypic variability in conformation traits in Brazilian Holstein cattle.

Advances in single-cell transcriptomics in animal research

Understanding biological mechanisms is fundamental for improving animal production and health to meet the growing demand for high-quality protein. As an emerging biotechnology, single-cell transcriptomics has been gradually applied in diverse aspects of animal research, offering an effective method to study the gene expression of high-throughput single cells of different tissues/organs in animals. In an unprecedented manner, researchers have identified cell types/subtypes and their marker genes, inferred cellular fate trajectories, and revealed cell‒cell interactions in animals using single-cell transcriptomics. In this paper, we introduce the development of single-cell technology and review the processes, advancements, and applications of single-cell transcriptomics in animal research. We summarize recent efforts using single-cell transcriptomics to obtain a more profound understanding of animal nutrition and health, reproductive performance, genetics, and disease models in different livestock species. Moreover, the practical experience accumulated based on a large number of cases is highlighted to provide a reference for determining key factors (e.g., sample size, cell clustering, and cell type annotation) in single-cell transcriptomics analysis. We also discuss the limitations and outlook of single-cell transcriptomics in the current stage. This paper describes the comprehensive progress of single-cell transcriptomics in animal research, offering novel insights and sustainable advancements in agricultural productivity and animal health.

Exploring the effects of Hippo signaling pathway on rumen epithelial proliferation

BACKGROUND

The current understanding to the mechanism of rumen development is limited. We hypothesized that the Hippo signaling pathway controlled the proliferation of rumen epithelium (RE) during postnatal development. In the present study, we firstly tested the changes of the Hippo signaling pathway in the RE during an early growing period from d5 to d25, and then we expanded the time range to the whole preweaning period (d10-38) and one week post weaning (d45). An in vitro experiment was also carried out to verify the function of Hippo signaling pathway during RE cell proliferation.

RESULTS

In the RE of lambs from d5 to d25, the expression of baculoviral IAP repeat containing (BIRC3/5) was increased, while the expressions of large tumor suppressor kinase 2 (LATS2), TEA domain transcription factor 3 (TEAD3), axin 1 (AXIN1), and MYC proto-oncogene (MYC) were decreased with rumen growth. From d10 to d38, the RE expressions of BIRC3/5 were increased, while the expressions of LATS2 and MYC were decreased, which were similar with the changes in RE from d5 to d25. From d38 to d45, different changes were observed, with the expressions of LATS1/2, MOB kinase activator 1B (MOB1B), and TEAD1 increased, while the expressions of MST1 and BIRC5 decreased. Correlation analysis showed that during the preweaning period, the RE expressions of BIRC3/5 were positively correlated with rumen development variables, while LAST2 was negatively correlated with rumen development variables. The in vitro experiment validated the changes of LATS2 and BIRC3/5 in the proliferating RE cells, which supported their roles in RE proliferation during preweaning period.

CONCLUSIONS

Our results suggest that the LATS2-YAP1-BIRC3/5 axis participates in the RE cell proliferation and promotes rumen growth during the preweaning period.

Genome-Wide Acetylation Modification of H3K27ac in Bovine Rumen Cell Following Butyrate Exposure

Butyrate contributes epigenetically to the changes in cellular function and tissue development of the rumen in ruminant animals, which might be achieved by its genetic or epigenetic regulation of gene expression. To explore the role of butyrate on bovine rumen epithelial function and development, this study characterized genome-wide H3K27ac modification changes and super-enhancer profiles in rumen epithelial primary cells (REPC) induced with butyrate by ChIP-seq, and analyzed its effects on gene expression and functional pathways by integrating RNA-seq data. The results showed that genome-wide acetylation modification was observed in the REPC with 94,675 and 48,688 peaks in the butyrate treatment and control group, respectively. A total of 9750 and 5020 genes with increased modification (H3K27ac-gain) and decreased modification (H3K27ac-loss) were detected in the treatment group. The super-enhancer associated genes in the butyrate-induction group were involved in the AMPK signaling pathway, MAPK signaling pathway, and ECM-receptor interaction. Finally, the up-regulated genes (PLCG1, CLEC3B, IGSF23, OTOP3, ADTRP) with H3K27ac gain modification by butyrate were involved in cholesterol metabolism, lysosome, cell adhesion molecules, and the PI3K-Akt signaling pathway. Butyrate treatment has the role of genome-wide H3K27ac acetylation on bovine REPC, and affects the changes in gene expression. The effect of butyrate on gene expression correlates with the acetylation of the H3K27ac level. Identifying genome-wide acetylation modifications and expressed genes of butyrate in bovine REPC cells will expand the understanding of the biological role of butyrate and its acetylation.

Investigation of rumen long noncoding RNA before and after weaning in cattle

Background

This study aimed to identify long non-coding RNA (lncRNA) from the rumen tissue in dairy cattle, explore their features including expression and conservation levels, and reveal potential links between lncRNA and complex traits that may indicate important functional impacts of rumen lncRNA during the transition to the weaning period.

Results

A total of six cattle rumen samples were taken with three replicates from before and after weaning periods, respectively. Total RNAs were extracted and sequenced with lncRNA discovered based on size, coding potential, sequence homology, and known protein domains. As a result, 404 and 234 rumen lncRNAs were identified before and after weaning, respectively. However, only nine of them were shared under two conditions, with 395 lncRNAs found only in pre-weaning tissues and 225 only in post-weaning samples. Interestingly, none of the nine common lncRNAs were differentially expressed between the two weaning conditions. LncRNA averaged shorter length, lower expression, and lower conservation scores than the genome overall, which is consistent with general lncRNA characteristics. By integrating rumen lncRNA before and after weaning with large-scale GWAS results in cattle, we reported significant enrichment of both pre- and after-weaning lncRNA with traits of economic importance including production, reproduction, health, and body conformation phenotypes.

Conclusions

The majority of rumen lncRNAs are uniquely expressed in one of the two weaning conditions, indicating a functional role of lncRNA in rumen development and transition of weaning. Notably, both pre- and post-weaning lncRNA showed significant enrichment with a variety of complex traits in dairy cattle, suggesting the importance of rumen lncRNA for cattle performance in the adult stage. These relationships should be further investigated to better understand the specific roles lncRNAs are playing in rumen development and cow performance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08758-4.

Transcriptome Analysis of Bovine Rumen Tissue in Three Developmental Stages

Rumen development is a crucial physiological challenge for ruminants. However, the molecular mechanism regulating rumen development has not been clearly elucidated. In this study, we investigated genes involved in rumen development in 13 rumen tissues from three developmental stages (birth, youth, and adult) using RNA sequencing. We identified that 6,048 genes were differentially expressed among three developmental stages. Using weighted correlation network analysis, we found that 12 modules were significantly associated with developmental stages. Functional annotation and protein–protein interaction (PPI) network analysis revealed that CCNB1, CCNB2, IGF1, IGF2, HMGCL, BDH1, ACAT1, HMGCS2, and CREBBP involved in rumen development. Integrated transcriptome with GWAS information of carcass weight (CW), stomach weight (SW), marbling score (MS), backfat thickness (BFT), ribeye area (REA), and lean meat weight (LMW), we found that upregulated DEGs (fold change 0∼1) in birth–youth comparison were significantly enriched with GWAS signals of MS, downregulated DEGs (fold change >3) were significantly enriched with GWAS signals of SW, and fold change 0∼1 up/downregulated DEGs in birth–adult comparison were significantly enriched with GWAS signals of CW, LMW, REA, and BFT. Furthermore, we found that GWAS signals for CW, LMW, and REA were enriched in turquoise module, and GWAS signals for CW was enriched in lightgreen module. Our study provides novel insights into the molecular mechanism underlying rumen development in cattle and highlights an integrative analysis for illustrating the genetic architecture of beef complex traits.

Integration of Long Non-Coding RNA and mRNA Profiling Reveals the Mechanisms of Different Dietary NFC/NDF Ratios Induced Rumen Development in Calves

The aim of the present study was to explore the effects of dietary non-fibrous carbohydrate to neutral detergent fiber (NFC/NDF) ratios on rumen development of calves, and to investigate the mechanisms by integrating of lncRNA and mRNA profiling. Forty-five weaned Charolais hybrid calves [body weight = 94.38 ± 2.50 kg; age = 70 ± 2.69 d] were randomly assigned to 1 of 3 treatment groups with different dietary NFC/NDF ratios: 1.10 (H group), 0.94 (M group) and 0.60 (L group), respectively. The ventral sac of the rumen was sampled for morphological observation and transcriptional sequencing. The average daily gain of calves in the high NFC/NDF ratio group was significantly higher than that in other groups (p < 0.05). Papillae width was largest in high NFC/NDF ratio group calves (p < 0.05). Identified differentially expressed genes that were significantly enriched in pathways closely related to rumen epithelial development included focal adhesion, Wingless-int signaling pathway, thyroid hormone signaling pathway, regulation of actin cytoskeleton and cGMP-PKG signaling pathway. The lncRNA-mRNA network included XLOC_068691 and MOAB, XLOC_023657 and DKK2, XLOC_064331 and PPP1R12A which we interpret to mean they have important regulatory roles in calve rumen development. These findings will serve as a theoretical basis for further analysis of the molecular genetic mechanism of dietary factors affecting rumen development in calves.

Functional annotation of regulatory elements in cattle genome reveals the roles of extracellular interaction and dynamic change of chromatin states in rumen development during weaning

We profiled landscapes of bovine regulatory elements and explored dynamic changes of chromatin states in rumen development during weaning. The regulatory elements (15 chromatin states) and their coordinated activities in cattle were defined through genome-wide profiling of four histone modifications, CTCF-binding, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial tissues. Each chromatin state presented specific enrichment for sequence ontology, methylation, trait-associated variants, transcription, gene expression-associated variants, selection signatures, and evolutionarily conserved elements. During weaning, weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states and occurred in tandem with significant variations in gene expression and DNA methylation, significantly associated with stature, production, and reproduction economic traits. By comparing with in vitro cultured epithelial cells and in vivo rumen tissues, we showed the commonness and uniqueness of these results, especially the roles of cell interactions and mitochondrial activities in tissue development.

Transcriptional Reprogramming in Rumen Epithelium during the Developmental Transition of Pre-Ruminant to the Ruminant in Cattle

Simple Summary

The rumen is the critical organ mediating nutrient uptake and use in cattle. Health development is essential to ensure animal feed efficiency. In this report, we present an analysis of gene expression dynamic in rumen epithelium during the transition from pre-ruminant to ruminant in cattle fed with hay or concentrated diets at weaning. The global shifts in gene expression and correlated transcription factors activities indicate transcriptional reprogramming during weaning. Transcriptional reprogramming in rumen epithelial tissue reflects critical nutrient-gene interactions occurring during the developmental progression. The results unveiled that nutrient-gene interactions compel transcriptional reprogramming. Our findings also suggest that this transcriptional reprogramming is the molecular basis of the transitional development of pre-ruminant to the ruminant in cattle.

Abstract

We present an analysis of transcriptomic dynamics in rumen epithelium of 18 Holstein calves during the transition from pre-rumination to rumination in cattle-fed hay or concentrated diets at weaning. Three calves each were euthanized at 14 and 42 d of age to exemplify preweaning, and six calves each were provided diets of either milk replacer and grass hay or calf starter to introduce weaning. The two distinct phases of rumen development and function in cattle are tightly regulated by a series of signaling events and clusters of effectors on critical pathways. The dietary shift from liquid to solid feeds prompted the shifting of gene activity. The number of differentially expressed genes increased significantly after weaning. Bioinformatic analysis revealed gene activity shifts underline the functional transitions in the ruminal epithelium and signify the transcriptomic reprogramming. Gene ontogeny (GO) term enrichment shows extensively activated biological functions of differentially expressed genes in the ruminal epithelium after weaning were predominant metabolic functions. The transcriptomic reprogramming signifies a correlation between gene activity and changes in metabolism and energy production in the rumen epithelium, which occur at weaning when transitioning from glucose use to VFA use by epithelium during the weaning.

Single-cell transcriptomic analyses of dairy cattle ruminal epithelial cells during weaning

Using the 10× Genomics Chromium Controller, we obtained scRNA-seq data of 5064 and 1372 individual cells from two Holstein calf ruminal epithelial tissues before and after weaning, respectively. We detected six distinct cell clusters, designated their cell types, and reported their marker genes. We then examined these clusters' underlining cell types and relationships by performing cell cycle, pseudotime trajectory, regulatory network, weighted gene co-expression network and gene ontology analyses. By integrating these cell marker genes with Holstein GWAS signals, we found they were enriched for animal production and body conformation traits. Finally, we confirmed their cell identities by comparing them with human and mouse stomach epithelial cells. This study presents an initial effort to implement single-cell transcriptomic analysis in cattle, and demonstrates ruminal tissue epithelial cell types and their developments during weaning, opening the door for new discoveries about tissue/cell type roles in complex traits at single-cell resolution.