Comparison of MicroRNA Transcriptomes Reveals the Association between MiR-148a-3p Expression and Rumen Development in Goats

Identification of candidate genes associated with milk production and mastitis based on transcriptome-wide association study

Genetic research for the assessment of mastitis and milk production traits simultaneously has a long history. The main issue that arises in this context is the known existence of a positive correlation between the risk of mastitis and lactation performance due to selection. The transcriptome-wide association study (TWAS) approach endeavors to combine the expression quantitative trait loci and genome-wide association study summary statistics to decode complex traits or diseases. Accordingly, we used the farmgtex project results as a complete bovine database for mastitis and milk production. The results of colocalization and TWAS approaches were used for the detection of functional associated candidate genes with milk production and mastitis traits on multiple tissue-based transcriptome records. Also, we used the david database for gene ontology to identify significant terms and associated genes. For the identification of interaction networks, the genemania and string databases were used. Also, the available z-scores in TWAS results were used for the calculation of the correlation between tissues. Therefore, the present results confirm that LYNX1, DGAT1, C14H8orf33, and LY6E were identified as significant genes associated with milk production in eight, six, five, and five tissues, respectively. Also, FBXL6 was detected as a significant gene associated with mastitis trait. CLN3 and ZNF34 genes emerged via both the colocalization and TWAS approaches as significant genes for milk production trait. It is expected that TWAS and colocalization can improve our perception of the potential health status control mechanism in high-yielding dairy cows.

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.

MicroRNAs in Ruminants and Their Potential Role in Nutrition and Physiology

The knowledge of how diet choices, dietary supplements, and feed intake influence molecular mechanisms in ruminant nutrition and physiology to maintain ruminant health, is essential to attain. In the present review, we focus on the role of microRNAs in ruminant health and disease; additionally, we discuss the potential of circulating microRNAs as biomarkers of disease in ruminants and the state of technology for their detection, also considering the major difficulties in the transition of biomarker development from bench to clinical practice. MicroRNAs are an inexhaustible class of endogenous non-protein coding small RNAs of 18 to 25 nucleotides that target either the 3' untranslated (UTR) or coding region of genes, ensuring a tight post-transcriptionally controlled regulation of gene expression. The development of new "omics" technologies facilitated a fresh perspective on the nutrition-to-gene relationship, incorporating more extensive data from molecular genetics, animal nutrition, and veterinary sciences. MicroRNAs might serve as important regulators of metabolic processes and may present the inter-phase between nutrition and gene regulation, controlled by the diet. The development of biomarkers holds the potential to revolutionize veterinary practice through faster disease detection, more accurate ruminant health monitoring, enhanced welfare, and increased productivity. Finally, we summarize the latest findings on how microRNAs function as biomarkers, how technological paradigms are reshaping this field of research, and how platforms are being used to identify novel biomarkers. Numerous studies have demonstrated a connection between circulating microRNAs and ruminant diseases such as mastitis, tuberculosis, foot-and-mouth disease, fasciolosis, and metabolic disorders. Therefore, the identification and analysis of a small number of microRNAs can provide crucial information about the stage of a disease, etiology, and prognosis.

LncRNA-mRNA modules involved in goat rumen development: Insights from genome-wide transcriptome profiling

The rumen is an essential digestive and absorption organ of ruminants. During fetal life, lactation, and post-weaning period, goat rumen undergoes drastic morphological and metabolic-functional changes triggered by potential regulated genes and non-coding RNA molecules. As the essential regulatory factors, long non-coding RNAs (lncRNAs) have vital functions in various biological activities. However, their roles during rumen development are still poorly explored in ruminants. To explore the genome-wide expression profiles of lncRNAs and mRNAs in the goat rumens, we generated 5,007 lncRNAs and 19,738 mRNAs identified during the fetal and prepubertal stages by the high-throughput RNA sequencing. Notably, 365 lncRNAs and 2,877 mRNAs were considered to be differentially expressed. The weighted gene co-expression network analysis and functional analysis were performed to explore the regulatory roles of those differentially expressed molecules. The cis-and trans-target genes of differently expressed lncRNAs were enriched for pathways related to focal adhesion, cGMP-PKG signaling pathway, alpha-linolenic acid metabolism, arachidonic acid metabolism, and fat digestion and absorption. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes analyses showed that the differently expressed genes mainly participated in mitotic cytokinesis, desmosome, fatty acid degradation, cell adhesion molecules, and fatty acid metabolism. The prediction of lncRNA-mRNA interaction networks further revealed transcripts potentially involved in rumen development. The present study profiles a global overview of lncRNAs and mRNAs during rumen development. Our findings provide valuable resources for genetic regulation and molecular mechanisms of rumen development in ruminants.

Characterization of presence and activity of microRNAs in the rumen of cattle hints at possible host-microbiota cross-talk mechanism

MicroRNAs (miRNAs), as important post-transcriptional regulators, are ubiquitous in various tissues. The aim of this exploratory study was to determine the presence of miRNAs in rumen fluid, and to investigate the possibility of miRNA-mediated cross-talk within the ruminal ecosystem. Rumen fluid samples from four cannulated Holstein cows were collected during two feeding regimes (forage and high-grain diet) and DNA and RNA were extracted for amplicon and small RNA sequencing. Epithelial biopsies were simultaneously collected to investigate the co-expression of miRNAs in papillae and rumen fluid. We identified 377 miRNAs in rumen fluid and 638 in rumen papillae, of which 373 were shared. Analysis of microbiota revealed 20 genera to be differentially abundant between the two feeding regimes, whereas no difference in miRNAs expression was detected. Correlations with at least one genus were found for 170 miRNAs, of which, 39 were highly significant (r > |0.7| and P < 0.01). Both hierarchical clustering of the correlation matrix and WGCNA analysis identified two main miRNA groups. Putative target and functional prediction analysis for the two groups revealed shared pathways with the predicted metabolic activities of the microbiota. Hence, our study supports the hypothesis of a cross-talk within the rumen at least partly mediated by miRNAs.

Developmental Alterations of Colonic microRNA Profiles Imply Potential Biological Functions in Kid Goats

The colon is a crucial digestive organ of the hind gut in ruminants. The bacterial diversity and mucosal immune maturation in this region are related to age. However, whether the microRNA expression in the colon of goats is affected by age is still unclear. In the current study, we analyzed the transcriptomes of colon microRNAs during preweaning (Day 10 and Day 25) and postweaning (Day 31). A total of 1572 microRNAs were identified in the colon tissues. Of these, 39 differentially expressed microRNAs (DEmiRNAs) and 88 highly expressed microRNAs (HEmiRNAs) were screened. The target genes regulated by the DEmiRNAs and HEmiRNAs were commonly enriched in the MAPK signaling pathway, Wnt signaling pathway, Hippo signaling pathway, cell adhesion molecules, focal adhesion, and adherens junction. Remarkably, the targeted genes of the DEmiRNAs were highly enriched for the prevention of microbial invasion via the Erbb−MAPK network while the targeted genes of HEmiRNAs contributed to the permeable barrier maintenance and cell damage surveillance. Additionally, there were eight different expression profiles of 87 dynamic miRNAs, in which approximately half of them were affected by age. Taken together, our study reveals the different roles of DEmiRNAs, HEmiRNAs, and dynamic microRNAs in the development of the colon and gives new insights into the regulatory mechanism of colon development in goats.

Bovine rumen epithelial miRNA-mRNA dynamics reveals post-transcriptional regulation of gene expression upon transition to high-grain feeding and phytogenic supplementation

The rumen epithelium has a pivotal role in nutrient uptake and host health. This study aimed to explore the role of microRNAs (miRNAs) in the epithelial transcriptome during diet transition from forage to high-grain feeding and the modulation through supplementation with a phytogenic feed additive. Rumen biopsies were collected from 9 ruminally-cannulated non-lactating Holstein cows fed a baseline forage diet (FD) and then transitioned to high-grain feeding (HG; 65% concentrate on a dry matter basis). Cows were randomly allocated into a control group (CON, n = 5) and a group supplemented with a phytogenic feed additive (PHY, n = 4). MiRNA and mRNA sequencing was performed in parallel and transcripts were analyzed for differential expression, pathway enrichment analysis, and miRNA-mRNA interaction networks. We identified 527 miRNAs shared by all samples of the rumen epithelium, from which, bta-miR-21-5p, bta-miR-143 and bta-miR-24-3p were the most expressed. Six miRNAs were differentially expressed between CON and PHY and 8 miRNAs between FD and HG feeding, which were mainly associated with fat metabolism. Transcriptome analysis identified 9481 differentially expressed genes (DEGs) between FD and HG, whereas PHY supplementation resulted in 5 DEGs. DEGs were mainly involved in epithelium development and morphogenesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with tricarboxylic acid and short chain fatty acid (SCFA) metabolism were enriched in DEGs between diets. MiRNA target prediction and anti-correlation analysis was used to construct networks and identify DEGs targeted by DE miRNAs responsive to diet or PHY. This study allowed the identification of potential miRNA regulation mechanisms of gene expression during transition from FD to HG feeding and phytogenic supplementation, evidencing a direct role of miRNAs in host responses to nutrition.

Transcriptomic Profiling of Circular RNAs in the Goat Rumen During Fetal and Prepubertal Period

Circular RNAs (circRNAs) are key regulatory factors with vital functions in various biological activities. However, little has been reported concerning the genetic regulation of circRNAs during rumen development in goats. The aim of this study was to identify the genome-wide expression profiles of circRNAs in the rumen of goats during fetal development and before and after weaning. Histological morphology showed that from the fetal period (days 60 and 135 of gestation) to the prepuberal period (days 60 and 150 of age) the rumen papilla developed gradually, and the thickness of the rumen muscular layer increased. A total of 11,149 circRNAs were identified in the four development stages by RNA-sequencing. From this, 1,518 were differentially expressed circRNAs (DECs). Fifty-eight DECs were up-regulated from 60 to 135 days of gestation, and 93 from day 135 of pregnancy to 30 days after birth. A large proportion (598) of DECs were down-regulated from day 135 of gestation to 30 days after birth. The expression levels of six randomly selected circRNAs were validated by qPCR, and their back-splicing junction (BSJ) sites were also confirmed. Ontology and pathway analyses revealed that the parental genes of DECs were mainly involved in the signaling pathways related to cell proliferation and apoptosis. The interaction network of circRNAs with their target miRNAs showed its involvement in cell proliferation and apoptosis signaling pathways. In conclusion, we identified the genome-wide expression profiles of circRNAs in the rumen of goats during fetal development and before and after weaning. These results provide a basis for further study on the regulatory effect of circRNAs on the development of rumen tissues.

Genome-Wide Identification of Reference Genes for Reverse-Transcription Quantitative PCR in Goat Rumen

Simple Summary

The rumen plays an essential role as a digestive organ and serves as the primary site of energy substrate absorption for the productive ruminants. Understanding gene expression profiles is necessary to explore the intrinsic regulatory mechanisms of rumen development in goats. The selection of suitable reference genes (RGs) was the primary assay before the real-time quantitative PCR (RT-qPCR). We identified sixteen genome-wide candidate RGs for normalization of gene expression assessments in goat rumen tissues. We demonstrate that the RGs selected (RPS4X and RPS6) were more stably expressed than the commonly used HKGs (ACTB and GAPDH) in goat rumen tissues, suggesting that the ribosomal protein gene family may be another source for the RG pool.

Abstract

As the largest chamber of the ruminant stomach, the rumen not only serves as the principal absorptive surface and nutrient transport pathway from the lumen into the animal, but also plays an important short-chain fatty acid (SCFA) metabolic role in addition to protective functions. Accurate characterization of the gene expression profiles of genes of interest is essential to the exploration of the intrinsic regulatory mechanisms of rumen development in goats. Thus, the selection of suitable reference genes (RGs) is an important prerequisite for real-time quantitative PCR (RT-qPCR). In the present study, 16 candidate RGs were identified from our previous transcriptome sequencing of caprine rumen tissues. The quantitative expressions of the candidate RGs were measured using the RT-qPCR method, and the expression stability of the RGs was assessed using the geNorm, NormFinder, and BestKeeper programs. GeNorm analysis showed that the M values were less than 0.5 for all the RGs except GAPT4, indicating that they were stably expressed in the rumen tissues throughout development. RPS4X and RPS6 were the two most stable RGs. Furthermore, the expressions of two randomly selected target genes (IGF1 and TOP2A), normalized by the selected most stable RGs (RPS4X and RPS6), were consistent with the results of RNA sequencing, while the use of GAPDH and ACTB as RGs resulted in altered profiles. Overall, RPS4X and RPS6 showed the highest expression stability and the lowest coefficients of variation, and could be used as the optimal reference combination for quantifying gene expression in rumen tissues via RT-qPCR analysis.

Identification of microRNAs in Silver Carp (Hypophthalmichthys molitrix) Response to Hypoxia Stress

Simple Summary

Hypoxia stress is one of the main problems in silver carp (Hypophthalmichthys molitrix) culture. Severe hypoxia stress can lead to damage and even death of silver carp. Therefore, it is very important to explore how silver carp adapt to and respond to hypoxia stress. MicroRNAs play an important role in a series of important life activities in organisms. In this study, the differentiallyexpressed miRNAs were screened from a mixed pool of liver, brain, heart and gill of silver carp under different levels of hypoxia stress by high-throughput sequencing. Our findings provided new insights to further study the miRNA regulatory mechanism and molecular characteristics of anoxic response in silver carp.

Abstract

Hypoxia is one of the serious stresses in fish culture, which can lead to physical and morphological changes, and cause injury and even death to fish. Silver carp (Hypophthalmichthys molitrix) is an important economic fish and widely distributed in China. MicroRNA is a kind of endogenous non-coding single-stranded small RNA, which is involved in cell development, and immune response and gene expression regulation. In this study, silver carp were kept in the closed containers for hypoxia treatment by spontaneous oxygen consumption. The samples of heart, brain, liver and gill were collected, and the total RNAs extracted separately from the four tissues were mixed in equal amounts according to the concentration. Afterwards, the RNA pool was constructed for high-throughput sequencing, and based on the small RNA sequencing, the differentially expressed microRNAs were identified. Furthermore, their target gene prediction and enrichment analyses were carried out. The results showed that a total of 229 known miRNAs and 391 putative novel miRNAs were identified, which provided valuable resources for further study on the regulatory mechanism of miRNAs in silver carp under hypoxia stress. The authors verified 16 differentially expressed miRNAs by qRT-PCR, and the results were consistent with small RNA sequencing (sRNA-seq). The predicted target genes number of differentially expressed miRNAs was 25,146. GO and KEGG functional enrichment analysis showed that these target genes were mainly involved in the adaption of hypoxia stress in silver carp through biological regulation, catalytic activity and apoptosis. This study provides references for further study of interaction between miRNAs and target genes, and the basic data for the response mechanism under hypoxia stress in silver carp.

Differentially Expressed MiRNAs of Goat Submandibular Glands Among Three Developmental Stages Are Involved in Immune Functions

Submandibular glands (SMGs) are one of the primary components of salivary glands in goats. The proteins and biologically active substances secreted by the SMGs change with growth and development. Our previous studies showed that most of the differentially expressed genes in the SMGs of goats at different developmental stages are involved in immune-related signaling pathways, but the miRNA expression patterns in the same tissues are unknown. The aim of this study was to reveal the expression profile of miRNAs at three different developmental stages, detect differentially expressed miRNAs (DE miRNAs) and predict disease-related DE miRNAs. SMG tissue samples were collected from groups of 1-month-old kids, 12-month-old maiden goats and 24-month-old adult goats (three samples from each group), and high-throughout transcriptome sequencing was conducted. A total of 178, 241 and 7 DE miRNAs were discovered between 1-month-old kids and 12-month-old maiden goats, between 1-month-old kids and 24-month-old adult goats, and between 12-month-old maiden goats and 24-month-old adult goats, respectively. Among these DE miRNAs, 88 DE miRNAs with medium or high expression levels (TPM ≥50) were classified into five expression pattern clusters. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that some of the predicted target genes of the DE miRNAs in the five clusters were enriched in disease-related GO terms and pathways. MiRNA target genes in significant pathways were significantly enriched in Hepatitis B (FDR = 9.03E-10) and Pathways in cancer (FDR = 4.2E-10). Further analysis was performed with a PPI network, and 10 miRNAs were predicted to play an important role in the occurrence and prevention of diseases during the growth and development of goats.