Transcriptome Analysis of Three Sheep Intestinal Regions reveals Key Pathways and Hub Regulatory Genes of Large Intestinal Lipid Metabolism

Integrated multi-omics reveals novel microbe-host lipid metabolism and immune interactions in the donkey hindgut

Evidence has shown that gut microbiota play a key role in host metabolism and health; however, little is known about the microbial community in the donkey hindgut as well as the interactions that occur between gut microbes and the host. This study aimed to explore the gut microbiome differences by analyzing the microbial community and differentially expressed genes (DEGs) related to lipid metabolism and the immune system along the donkey hindgut. The hindgut tissues (cecum, ventral colon, and dorsal colon) were separated, and the contents of each section were collected from six male donkeys for multi-omics analysis. There were significant differences in terms of dominant bacteria among the three sections, especially between the cecum and dorsal colon sites. For instance, species belonging to Prevotella and Treponema were most abundant in the cecum, while the Clostridiales_bacterium, Streptococcus_equinus, Ruminococcaceae_bacterium, etc., were more abundant in the dorsal colon. Apart from propionate, the concentrations of acetate, isobutyrate, valerate and isovalerate were all lower in the cecum than in the dorsal colon (p < 0.05). Furthermore, we identified some interesting DEGs related to lipid metabolism (e.g., ME1, MBOAT1, ACOX1, ACOX2 and LIPH) and the immune system (e.g., MUC3B, mucin-2-like, IL17RC, IL1R2, IL33, C1QA, and MMP9) between the cecum and dorsal colon and found that the PPAR pathway was mainly enriched in the cecum. Finally, we found a complex relationship between the gut microbiome and gene expression, especially with respect to the immune system, and combined with protein-protein interaction (PPI) data, suggesting that the PPAR pathway might be responsible, at least in part, for the role of the hindgut microbiota in the donkeys' gut homeostasis. Our data provide an in-depth understanding of the interaction between the microbiota and function in the healthy equine hindgut and may also provide guidance for improving animal performance metrics (such as product quality) and equine welfare.

Expression characteristics of piRNAs in ovine luteal phase and follicular phase ovaries

PIWI-interacting RNAs (piRNAs), as a novel class of small non-coding RNAs that have been shown to be indispensable in germline integrity and stem cell development. However, the expressed characteristics and regulatory roles of piRNAs during different reproductive phases of animals remain unknown. In this study, we investigated the piRNAs expression profiles in ovaries of sheep during the luteal phase (LP) and follicular phase (FP) using the Solexa sequencing technique. A total of 85,219 and 1,27,156 piRNAs tags were identified in ovine ovaries across the two phases. Most expressed piRNAs start with uracil. piRNAs with a length of 24 nt or 27–29 nts accounted for the largest proportion. The obvious ping-pong signature appeared in the FP ovary. The piRNA clusters in the sheep ovary were unevenly distributed on the chromosomes, with high density on Chr 3 and 1. For genome distribution, piRNAs in sheep ovary were mainly derived from intron, CDS, and repeat sequence regions. Compared to the LP ovary, a greater number of expressed piRNA clusters were detected in the FP ovary. Simultaneously, we identified 271 differentially expressed (DE) piRNAs between LP and FP ovaries, with 96 piRNAs upregulated and 175 piRNAs downregulated, respectively. Functional enrichment analysis (GO and KEGG) indicated that their target genes were enriched in reproduction-related pathways including oocyte meiosis, PI3K-Akt, Wnt, and TGF-β signaling pathways. Together, our results highlighted the sequence and expression characteristics of the piRNAs in the sheep ovary, which will help us understand the roles of piRNAs in the ovine estrus cycle.

Genomic sequencing analysis reveals copy number variations and their associations with economically important traits in beef cattle

Copy number variation (CNV) represents a major source of genetic variation, which may have potentially large effects, including alternating gene regulation and dosage, as well as contributing to gene expression and risk for normal phenotypic variability. We carried out a comprehensive analysis of CNV based on whole genome sequencing in Chinese Simmental beef cattle. Totally, we found 9313 deletion and 234 duplication events, covering 147.5 Mb autosomal regions. Within them, 257 deletion events of high frequency overlapped with 193 known RefGenes. Among these genes, we observed several genes were related to economically important traits, like residual feed intake, immune responding, pregnancy rate and muscle differentiation. Using a locus-based analysis, we identified 11 deletions and 1 duplication, which were significantly associated with three traits including carcass weight, tenderloin and longissimus muscle area. Our sequencing-based study provided important insights into investigating the association of CNVs with important traits in beef cattle.

Profiling of gene expression in methicillin-resistant Staphylococcus aureus in response to cyclo-(L-Val-L-Pro) and chloramphenicol isolated from Streptomyces sp., SUK 25 reveals gene downregulation in multiple biological targets

Chloramphenicol (CAP) and cyclo-(L-Val-L-Pro) were previously isolated from Streptomyces sp., SUK 25 which exhibited a high potency against methicillin-resistant Staphylococcus aureus (MRSA). This study aimed to profile gene expression of MRSA treated with CAP and cyclo-(L-Val-L-Pro) compounds using DNA microarray. Treatment of MRSA with CAP resulted in upregulation of genes involved in protein synthesis, suggesting the coping mechanism of MRSA due to the inhibition of protein synthesis effect from CAP. Most upregulated genes in cyclo-(L-Val-L-Pro) were putative genes with unknown functions. Interestingly, genes encoding ribosomal proteins, cell membrane synthesis, DNA metabolism, citric acid cycle and virulence were downregulated in MRSA treated with cyclo-(L-Val-L-Pro) compound, suggesting the efficacy of this compound in targeting multiple biological pathways. Contrary to CAP, with only a single target, cyclo-(L-Val-L-Pro) isolated from this study had multiple antimicrobial targets that can delay antibiotic resistance and hence is a potential antimicrobial agent of MRSA.

Integrated Hypothalamic Transcriptome Profiling Reveals the Reproductive Roles of mRNAs and miRNAs in Sheep

Early studies have provided a wealth of information on the functions of microRNAs (miRNAs). However, less is known regarding their functions in the hypothalamus involved in sheep reproduction. To explore the potential roles of hypothalamic messenger RNAs (mRNAs) and miRNAs in sheep without FecB mutation, in total, 172 and 235 differentially expressed genes (DEGs) and 42 and 79 differentially expressed miRNAs (DE miRNAs) were identified in polytocous sheep in the follicular phase versus monotocous sheep in the follicular phase (PF vs. MF) and polytocous sheep in the luteal phase versus monotocous sheep in the luteal phase (PL vs. ML), respectively, using RNA sequencing. We also identified several key mRNAs (e.g., POMC, GNRH1, PRL, GH, TRH, and TTR) and mRNA–miRNAs pairs (e.g., TRH co-regulated by oar-miR-379-5p, oar-miR-30b, oar-miR-152, oar-miR-495-3p, oar-miR-143, oar-miR-106b, oar-miR-218a, oar-miR-148a, and PRL regulated by oar-miR-432) through functional enrichment analysis, and the identified mRNAs and miRNAs may function, conceivably, by influencing gonadotropin-releasing hormone (GnRH) activities and nerve cell survival associated with reproductive hormone release via direct and indirect ways. This study represents an integral analysis between mRNAs and miRNAs in sheep hypothalamus and provides a valuable resource for elucidating sheep prolificacy.

Effect of Divergent Selection for Intramuscular Fat Content on Muscle Lipid Metabolism in Chickens

Simple Summary

Intramuscular fat is an important factor affecting meat quality and consumer acceptance. Appropriate increases in the intramuscular fat content contribute to the improvement of meat quality, and genetic selection is an effective method to increase the intramuscular fat content in chickens. In this study, chicken lines divergently selected for their intramuscular fat content were used to investigate the mechanisms behind differential intramuscular fat deposition. These results found in this study may contribute to the improvement of meat quality in chickens.

Abstract

Intramuscular fat (IMF)—an important factor affecting meat quality—can be appropriately increased by genetic selection. Chicken lines divergently selected for IMF content were used in this study to investigate the mechanisms behind differential IMF deposition. Sixty 15th generation chickens were genotyped using the IASCHICK 55K single nucleotide polymorphism (SNP) chip. After quality control, 59 chickens and 36,893 SNPs were available for subsequent analysis. Population structure assessment indicated that the lines were genetically differentiated. Based on the top 1% paired fixation index values, three pathways were significantly (p < 0.05) enriched, and nine genes were considered candidate genes for differential IMF deposition. Differences between the lines in the expressions of representative genes involved in the above pathways were detected in 16th generation chickens. This study suggests that genetic selection for increased IMF in the pectoralis major muscle may enhance fatty acid synthesis, transport, and esterification, and reduce triglyceride hydrolysis. The peroxisome proliferator-activated receptor (PPAR) signaling pathway, glycerolipid metabolism, and fatty acid degradation pathway may have contributed to the differences in IMF deposition between the lines. These results contribute to the understanding of the genetic mechanisms behind IMF deposition, and the improvement of chicken meat quality.

Transcriptome analysis reveals mechanism underlying the differential intestinal functionality of laying hens in the late phase and peak phase of production

Background

The compromised performance of laying hens in the late phase of production relative to the peak production was thought to be associated with the impairment of intestinal functionality, which plays essential roles in contributing to their overall health and production performance. In the present study, RNA sequencing was used to investigate differences in the expression profile of intestinal functionality-related genes and associated pathways between laying hens in the late phase and peak phase of production.

Results

A total of 104 upregulated genes with 190 downregulated genes were identified in the ileum (the distal small intestine) of laying hens in the late phase of production compared to those at peak production. These upregulated genes were found to be enriched in little KEGG pathway, however, the downregulated genes were enriched in the pathways of PPAR signaling pathway, oxidative phosphorylation and glutathione metabolism. Besides, these downregulated genes were mapped to several GO clusters in relation to lipid metabolism, electron transport of respiratory chain, and oxidation resistance. Similarly, there were lower activities of total superoxide dismutase, glutathione S-transferase and Na⁺/K⁺-ATPase, and reductions of total antioxidant capacity and ATP level, along with an elevation in malondialdehyde content in the ileum of laying hens in the late phase of production as compared with those at peak production.

Conclusions

The intestine of laying hens in the late phase of production were predominantly characterized by a disorder of lipid metabolism, concurrent with impairments of energy production and antioxidant property. This study uncovers the mechanism underlying differences between the intestinal functionality of laying hens in the late phase and peak phase of production, thereby providing potential targets for the genetic control or dietary modulation of intestinal hypofunction of laying hens in the late phase of production.

A Transcriptomic Study of the Tail Fat Deposition in Two Types of Hulun Buir Sheep According to Tail Size and Sex

Simple Summary

Based on tail types, Hulun Buir sheep were divided into two lines including small and big fat-tailed, but these two lines have similar genetic background. In this study, we investigated the morphology and transcription level differences of tail fat between these two lines. The RNA-seq analyses indicated several differentially expressed genes when compared between sexes or two tail sizes. Interestingly, we also found an obvious sex difference in the fat metabolism in Hulun Buir sheep. Two different co-expression networks were only shown either in male or in female sheep. Our findings will provide theoretical background in understanding the genetic mechanism of fat deposition in sheep.

Abstract

Hulun Buir sheep of similar genetic background were divided into two lines based on tail types: Small- and big fat-tailed. To explore the molecular mechanism of fat deposition in sheep tails, we firstly evaluated the morphology and transcription level differences of tail fat between these two lines. RNA-Seq technology was used to identify differentially expressed genes (DEGs) in phenotypic extremes of tail sizes. Five comparisons were performed taking into account two factors, sex and tail type. We screened out 373 DEGs between big-tailed and small-tailed Hulun Buir sheep, and 775 and 578 DEGs between two types of tails in male and female sheep, respectively. The results showed an obvious sex difference in the fat metabolism in sheep based on gene ontology (GO), pathway, and network analyses. Intriguingly, there were two different co-expression networks only respectively shown in male and female sheep, which were insulin-related network acting on upstream pathways and PPARG-related network effect in downstream pathways. Furthermore, these two networks were linked by a classic pathway of regulating adipogenesis. This is the first study to investigate the sex differences of fat metabolism in domestic animals, and it demonstrates a new experimental way to study fat metabolism. Our findings will provide theoretical background in understanding the tail-size phenotype in sheep and can be exploited in breeding small-tailed sheep.

Deep transcriptome analysis using RNA-Seq suggests novel insights into molecular aspects of fat-tail metabolism in sheep

Fat-tail content of sheep breeds is varied and the molecular mechanisms regulating fat-tail development have not been well characterized. Aiming at better identifying the important candidate genes and their functional pathways contributing to fat deposition in the tail, a comparative transcriptome analysis was performed between fat- (Lori-Bakhtiari) and thin-tailed (Zel) Iranian sheep breeds using RNA-seq. The experiment was conducted on six male lambs (three lambs per each breed) at seven months of age. Four different combinations of aligners and statistical methods including Hisat2 + edgeR, Hisat2 + DESeq2, STAR + edgeR and STAR + DESeq2 were used to identify the differentially expressed genes (DEGs). The DEGs were selected for functional enrichment analysis and protein-protein interaction (PPI) network construction. Module analysis was also conducted to mine the functional sub-networks from the PPI network. In total, 264 genes including 80 up- and 184 down-regulated genes were identified as DEGs. The RNA-Seq results were validated by Q-RT-PCR. Functional analysis of DEGs and the module analysis of PPI network demonstrated that in addition to pathways affecting lipid metabolism, a series of enriched functional terms related to "response to interleukin", "MAPK signaling pathways", "Wnt signaling pathway", "ECM-receptor interaction", "regulation of actin cytoskeleton", and "response to cAMP" might contribute to the deposition of fat in tails of sheep. Overall results using RNA-Seq analysis characterized important candidate genes involved in the fatty acid metabolism and regulation of fat deposition, suggesting novel insights into molecular aspects of fat-tail metabolism in sheep. Selected DEGs should be further investigated as potential markers associated with the fat-tail development in sheep breeds.

Genome-wide association study between CNVs and milk production traits in Valle del Belice sheep

Copy number variation (CNV) is a major source of genomic structural variation. The aim of this study was to detect genomic CNV regions (CNVR) in Valle del Belice dairy sheep population and to identify those affecting milk production traits. The GO analysis identified possible candidate genes and pathways related to the selected traits. We identified CNVs in 416 individuals genotyped using the Illumina OvineSNP50 BeadChip array. The CNV association using a correlation-trend test model was examined with the Golden Helix SVS 8.7.0 tool. Significant CNVs were detected when their adjusted p-value was <0.01 after false discovery rate (FDR) correction. We identified 7,208 CNVs, which gave 365 CNVRs after aggregating overlapping CNVs. Thirty-one CNVRs were significantly associated with one or more traits included in the analysis. All CNVRs, except those on OAR19, overlapped with quantitative trait loci (QTL), even if they were not directly related to the traits of interest. A total of 222 genes were annotated within the significantly associated CNVRs, most of which played important roles in biological processes related to milk production and health-related traits. Identification of the genes in the CNVRs associated with the studied traits will provide the basis for further investigation of their role in the metabolic pathways related to milk production and health traits.

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

One of the most significant physiological challenges to neonatal and juvenile ruminants is the development and establishment of the rumen. Using a subset of RNA-Seq data from our high-resolution atlas of gene expression in sheep (Ovis aries) we have provided the first comprehensive characterization of transcription of the entire gastrointestinal (GI) tract during the transition from pre-ruminant to ruminant. The dataset comprises 164 tissue samples from sheep at four different time points (birth, one week, 8 weeks and adult). Using network cluster analysis we illustrate how the complexity of the GI tract is reflected in tissue- and developmental stage-specific differences in gene expression. The most significant transcriptional differences between neonatal and adult sheep were observed in the rumen complex. Comparative analysis of gene expression in three GI tract tissues from age-matched sheep and goats revealed species-specific differences in genes involved in immunity and metabolism. This study improves our understanding of the transcriptomic mechanisms involved in the transition from pre-ruminant to ruminant by identifying key genes involved in immunity, microbe recognition and metabolism. The results form a basis for future studies linking gene expression with microbial colonization of the developing GI tract and provide a foundation to improve ruminant efficiency and productivity through identifying potential targets for novel therapeutics and gene editing.

Transcriptome analysis reveals differential gene expression in intramuscular adipose tissues of Jinhua and Landrace pigs

Meat is a rich source of protein, fatty acids and carbohydrates for human needs. In addition to necessary nutrients, high fat contents in pork increase the tenderness and juiciness of the meat, featuring diverse application in various dishes. This study investigated the transcriptomic profiles of intramuscular adipose tissues in Jinhua and Landrace pigs by employing advanced RNA sequencing. Results showed significant interesting to note that there were significant differences in the expression of genes. 1,632 genes showed significant differential expression, 837 genes were up-regulated and 195 genes were down-regulated. Variations in genes responsible for cell aggregation, extracellular matrix formation, cellular lipid catabolic process, and fatty acid binding strongly supported that both pig breeds feature variable fat and muscle metabolism. Certain differentially expressed genes are included in the pathway of mitogen-activated protein kinase signaling pathway, Ras signaling pathway and insulin pathway. Results from real-time quantitative polymerase chain reaction also validated the differential expression of 17 mRNAs between meats of the two pig breeds. Overall, these findings reveal significant differences in fat and protein metabolism of intramuscular adipose tissues of two pig breeds at the transcriptomic level and suggest diversification at the genetic level between breeds of the same species.

Integration of lncRNA and mRNA Transcriptome Analyses Reveals Genes and Pathways Potentially Involved in Calf Intestinal Growth and Development during the Early Weeks of Life

A better understanding of the factors that regulate growth and immune response of the gastrointestinal tract (GIT) of calves will promote informed management practices in calf rearing. This study aimed to explore genomics (messenger RNA (mRNA)) and epigenomics (long non-coding RNA (lncRNA)) mechanisms regulating the development of the rumen and ileum in calves. Thirty-two calves (≈5-days-old) were reared for 96 days following standard procedures. Sixteen calves were humanely euthanized on experiment day 33 (D33) (pre-weaning) and another 16 on D96 (post-weaning) for collection of ileum and rumen tissues. RNA from tissues was subjected to next generation sequencing and 3310 and 4217 mRNAs were differentially expressed (DE) between D33 and D96 in ileum and rumen tissues, respectively. Gene ontology and pathways enrichment of DE genes confirmed their roles in developmental processes, immunity and lipid metabolism. A total of 1568 (63 known and 1505 novel) and 4243 (88 known and 4155 novel) lncRNAs were detected in ileum and rumen tissues, respectively. Cis target gene analysis identified BMPR1A, an important gene for a GIT disease (juvenile polyposis syndrome) in humans, as a candidate cis target gene for lncRNAs in both tissues. LncRNA cis target gene enrichment suggested that lncRNAs might regulate growth and development in both tissues as well as posttranscriptional gene silencing by RNA or microRNA processing in rumen, or disease resistance mechanisms in ileum. This study provides a catalog of bovine lncRNAs and set a baseline for exploring their functions in calf GIT development.

Identification and characterization of microRNAs in the intestinal tissues of sheep (Ovis aries)

Sheep are small ruminants, and their long intestines exhibit high digestive and absorptive capacity in many different rearing conditions; however, the genetic bases of this characteristic remains unclear. MicroRNAs (miRNAs) play a major role in maintaining both intestinal morphological structure as well as in regulating the physiological functions of this organ. However, no study has reported on the miRNA expression profile in the intestinal tissues of sheep. Here, we analyzed and identified the miRNA expression profile of three different intestinal tissues (i.e., duodenum, cecum, and colon) of sheep (Ovis aries) using high-throughput sequencing and bioinformatic methods. In total, 106 known miRNAs were identified, 458 conserved miRNAs were detected, 192 unannotated novel miRNAs were predicted, and 195 differentially expressed miRNAs were found between the different tissues. Additionally, 3,437 candidate target genes were predicted, and 17 non-redundant significantly enriched GO terms were identified using enrichment analysis. A total of 99 candidate target genes were found to significantly enriched in 4 KEGG biological pathways. A combined regulatory network was constructed based on 92 metabolism-related candidate target genes and 65 differentially expressed miRNAs, among which 7 miRNAs were identified as hub miRNAs. Via these mechanisms, miRNAs may play a role in maintaining intestinal homeostasis and metabolism. This study helps to further explain the mechanisms that underlie differences in tissue morphology and function in three intestinal segments of sheep.