Comprehensive Transcriptional Profiling of the Gastrointestinal Tract of Ruminants from Birth to Adulthood Reveals Strong Developmental Stage Specific Gene Expression

The early life immune dynamics and cellular drivers at single-cell resolution in lamb forestomachs and abomasum

BACKGROUND

Four-chambered stomach including the forestomachs (rumen, reticulum, and omasum) and abomasum allows ruminants convert plant fiber into high-quality animal products. The early development of this four-chambered stomach is crucial for the health and well-being of young ruminants, especially the immune development. However, the dynamics of immune development are poorly understood.

RESULTS

We investigated the early gene expression patterns across the four-chambered stomach in Hu sheep, at 5, 10, 15, and 25 days of age. We found that forestomachs share similar gene expression patterns, all four stomachs underwent widespread activation of both innate and adaptive immune responses from d 5 to 25, whereas the metabolic function were significantly downregulated with age. We constructed a cell landscape of the four-chambered stomach using single-cell sequencing. Integrating transcriptomic and single-cell transcriptomic analyses revealed that the immune-associated module hub genes were highly expressed in T cells, monocytes and macrophages, as well as the defense-associated module hub genes were highly expressed in endothelial cells in the four-stomach tissues. Moreover, the non-immune cells such as epithelial cells play key roles in immune maturation. Cell communication analysis predicted that in addition to immune cells, non-immune cells recruit immune cells through macrophage migration inhibitory factor signaling in the forestomachs.

CONCLUSIONS

Our results demonstrate that the immune and defense responses of four stomachs are quickly developing with age in lamb's early life. We also identified the gene expression patterns and functional cells associated with immune development. Additionally, we identified some key receptors and signaling involved in immune regulation. These results help to understand the early life immune development at single-cell resolution, which has implications to develop nutritional manipulation and health management strategies based on specific targets including key receptors and signaling pathways.

Single-cell transcriptome and metagenome profiling reveals the genetic basis of rumen functions and convergent developmental patterns in ruminants

The rumen undergoes developmental changes during maturation. To characterize this understudied dynamic process, we profiled single-cell transcriptomes of about 308,000 cells from the rumen tissues of sheep and goats at 17 time points. We built comprehensive transcriptome and metagenome atlases from early embryonic to rumination stages, and recapitulated histomorphometric and transcriptional features of the rumen, revealing key transitional signatures associated with the development of ruminal cells, microbiota, and core transcriptional regulatory networks. In addition, we identified and validated potential cross-talk between host cells and microbiomes and revealed their roles in modulating the spatiotemporal expression of key genes in ruminal cells. Cross-species analyses revealed convergent developmental patterns of cellular heterogeneity, gene expression, and cell-cell and microbiome-cell interactions. Finally, we uncovered how the interactions can act upon the symbiotic rumen system to modify the processes of fermentation, fiber digestion, and immune defense. These results significantly enhance understanding of the genetic basis of the unique roles of rumen.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Establishment of SV40 Large T-Antigen-Immortalized Yak Rumen Fibroblast Cell Line and the Fibroblast Responses to Lipopolysaccharide

The yak lives in harsh alpine environments and the rumen plays a crucial role in the digestive system. Rumen-associated cells have unique adaptations and functions. The yak rumen fibroblast cell line (SV40T-YFB) was immortalized by introducing simian virus 40 large T antigen (SV40T) by lentivirus-mediated transfection. Further, we have reported the effects of lipopolysaccharide (LPS) of different concentrations on cell proliferation, extracellular matrix (ECM), and proinflammatory mediators in SV40T-YFB. The results showed that the immortalized yak rumen fibroblast cell lines were identified as fibroblasts that presented oval nuclei, a fusiform shape, and positive vimentin and SV40T staining after stable passage. Chromosome karyotype analysis showed diploid characteristics of yak (n = 60). LPS at different concentrations inhibited cell viability in a dose-dependent manner. SV40T-YFB treated with LPS increased mRNA expression levels of matrix metalloproteinases (MMP-2 and MMP-9), inflammatory cytokines (TNF-α, IL-1β, IL-6), and urokinase-type plasminogen activator system components (uPA, uPAR). LPS inhibits the expression of tissue inhibitors of metalloproteinases (TIMP-1 and TIMP-2), plasminogen activator inhibitor-2 (PAI-2), fibronectin (FN), anti-inflammatory factor IL-10, and collagen I (COL I) in SV40T-YFB. Overall, these results suggest that LPS inhibits cell proliferation and induces ECM degradation and inflammatory response in SV40T-YFB.

Microbiota-host crosstalk in the newborn and adult rumen at single-cell resolution

BACKGROUND

The rumen is the hallmark organ of ruminants, playing a vital role in their nutrition and providing products for humans. In newborn suckling ruminants milk bypasses the rumen, while in adults this first chamber of the forestomach has developed to become the principal site of microbial fermentation of plant fibers. With the advent of single-cell transcriptomics, it is now possible to study the underlying cell composition of rumen tissues and investigate how this relates the development of mutualistic symbiosis between the rumen and its epithelium-attached microbes.

RESULTS

We constructed a comprehensive cell landscape of the rumen epithelium, based on single-cell RNA sequencing of 49,689 high-quality single cells from newborn and adult rumen tissues. Our single-cell analysis identified six immune cell subtypes and seventeen non-immune cell subtypes of the rumen. On performing cross-species analysis of orthologous genes expressed in epithelial cells of cattle rumen and the human stomach and skin, we observed that the species difference overrides any cross-species cell-type similarity. Comparing adult with newborn cattle samples, we found fewer epithelial cell subtypes and more abundant immune cells, dominated by T helper type 17 cells in the rumen tissue of adult cattle. In newborns, there were more fibroblasts and myofibroblasts, an IGFBP3+ epithelial cell subtype not seen in adults, while dendritic cells were the most prevalent immune cell subtype. Metabolism-related functions and the oxidation-reduction process were significantly upregulated in adult rumen epithelial cells. Using 16S rDNA sequencing, fluorescence in situ hybridization, and absolute quantitative real-time PCR, we found that epithelial Desulfovibrio was significantly enriched in the adult cattle. Integrating the microbiome and metabolome analysis of rumen tissues revealed a high co-occurrence probability of Desulfovibrio with pyridoxal in the adult cattle compared with newborn ones while the scRNA-seq data indicated a stronger ability of pyroxidal binding in the adult rumen epithelial cell subtypes. These findings indicate that Desulfovibrio and pyridoxal likely play important roles in maintaining redox balance in the adult rumen.

CONCLUSIONS

Our integrated multi-omics analysis provides novel insights into rumen development and function and may facilitate the future precision improvement of rumen function and milk/meat production in cattle.

Integrated multi-omics of the gastrointestinal microbiome and ruminant host reveals metabolic adaptation underlying early life development

BACKGROUND

The gastrointestinal tract (GIT) microbiome of ruminants and its metabolic repercussions vastly influence host metabolism and growth. However, a complete understanding of the bidirectional interactions that occur across the host-microbiome axis remains elusive, particularly during the critical development stages at early life. Here, we present an integrative multi-omics approach that simultaneously resolved the taxonomic and functional attributes of microbiota from five GIT regions as well as the metabolic features of the liver, muscle, urine, and serum in sika deer (Cervus nippon) across three key early life stages.

RESULTS

Within the host, analysis of metabolites over time in serum, urine, and muscle (longissimus lumborum) showed that changes in the fatty acid profile were concurrent with gains in body weight. Additional host transcriptomic and metabolomic analysis revealed that fatty acid β-oxidation and metabolism of tryptophan and branched chain amino acids play important roles in regulating hepatic metabolism. Across the varying regions of the GIT, we demonstrated that a complex and variable community of bacteria, viruses, and archaea colonized the GIT soon after birth, whereas microbial succession was driven by the cooperative networks of hub populations. Furthermore, GIT volatile fatty acid concentrations were marked by increased microbial metabolic pathway abundances linked to mannose (rumen) and amino acids (colon) metabolism. Significant functional shifts were also revealed across varying GIT tissues, which were dominated by host fatty acid metabolism associated with reactive oxygen species in the rumen epithelium, and the intensive immune response in both small and large intestine. Finally, we reveal a possible contributing role of necroptosis and apoptosis in enhancing ileum and colon epithelium development, respectively.

CONCLUSIONS

Our findings provide a comprehensive view for the involved mechanisms in the context of GIT microbiome and ruminant metabolic growth at early life. Video Abstract.

Dynamics of rumen gene expression, microbiome colonization, and their interplay in goats

Background

Preweaned rumen development is vital for animal health and efficient fermentation. In this study, we integrated ruminal transcriptomic and metagenomic data to explore the dynamics of rumen functions, microbial colonization, and their functional interactions during the first 8 weeks of life in goats.

Results

The dynamic rumen transcriptomic and microbial profiles both exhibited two distinct phases during early rumen development. The differentially expressed genes of the rumen transcriptome between the two phases showed that the immune-related response was enriched in the first phase and nutrient-related metabolism was enriched in the second phase, whereas the differentially expressed genes of the rumen microbiome were enriched in bacteriocin biosynthesis and glycolysis/gluconeogenesis activities. The developmental shift in the rumen transcriptome (at d 21) was earlier than the feed stimulus (at d 25) and the shift in the rumen microbiome (at d 42). Additionally, 15 temporal dynamic rumen gene modules and 20 microbial modules were revealed by coexpression network analysis. Functional correlations between the rumen and its microbiome were primarily involved in rumen pH homeostasis, nitrogen metabolism and the immune response. Rumen gene modules associated with the microbial alpha diversity index were also enriched in the immune response process.

Conclusions

The present study touched the critical developmental process of rumen functions, microbial colonization and their functional interactions during preweaned development. Taken together, these results demonstrated that rumen development at the first phase is more likely a programmed process rather than stimulation from feed and the microbiome, while the shift of rumen metagenomes was likely regulated by both the diet and host. The intensive functional correlations between rumen genes and the microbiome demonstrated that synergistic processes occurred between them during early rumen development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07595-1.

A Mini-Atlas of Gene Expression for the Domestic Goat (Capra hircus)

Goats (Capra hircus) are an economically important livestock species providing meat and milk across the globe. They are of particular importance in tropical agri-systems contributing to sustainable agriculture, alleviation of poverty, social cohesion, and utilisation of marginal grazing. There are excellent genetic and genomic resources available for goats, including a highly contiguous reference genome (ARS1). However, gene expression information is limited in comparison to other ruminants. To support functional annotation of the genome and comparative transcriptomics, we created a mini-atlas of gene expression for the domestic goat. RNA-Seq analysis of 17 transcriptionally rich tissues and 3 cell-types detected the majority (90%) of predicted protein-coding transcripts and assigned informative gene names to more than 1000 previously unannotated protein-coding genes in the current reference genome for goat (ARS1). Using network-based cluster analysis, we grouped genes according to their expression patterns and assigned those groups of coexpressed genes to specific cell populations or pathways. We describe clusters of genes expressed in the gastro-intestinal tract and provide the expression profiles across tissues of a subset of genes associated with functional traits. Comparative analysis of the goat atlas with the larger sheep gene expression atlas dataset revealed transcriptional similarities between macrophage associated signatures in the sheep and goats sampled in this study. The goat transcriptomic resource complements the large gene expression dataset we have generated for sheep and contributes to the available genomic resources for interpretation of the relationship between genotype and phenotype in small ruminants.

Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.

A Gene Expression Atlas of the Domestic Water Buffalo (Bubalus bubalis)

The domestic water buffalo (Bubalus bubalis) makes a major contribution to the global agricultural economy in the form of milk, meat, hides, and draught power. The global water buffalo population is predominantly found in Asia, and per head of population more people depend upon the buffalo than on any other livestock species. Despite its agricultural importance, there are comparatively fewer genomic and transcriptomic resources available for buffalo than for other livestock species. We have generated a large-scale gene expression atlas covering multiple tissue and cell types from all major organ systems collected from three breeds of riverine water buffalo (Mediterranean, Pandharpuri and Bhadawari) and used the network analysis tool Graphia Professional to identify clusters of genes with similar expression profiles. Alongside similar data, we and others have generated for ruminants as part of the Functional Annotation of Animal Genomes Consortium; this comprehensive transcriptome supports functional annotation and comparative analysis of the water buffalo genome.