Genome-wide analysis of microRNAs identifies the lipid metabolism pathway to be a defining factor in adipose tissue from different sheep

AMPK Signaling Axis-Mediated Regulation of Lipid Metabolism: Ameliorative Effects of Sodium Octanoate on Intestinal Dysfunction in Hu Sheep

At the present stage, heavy metal pollution, led by environmental exposure to cadmium (Cd), has caused incalculable losses in animal husbandry. The potential value of caprylic acid as a medium- and long-chain fatty acid with a unique role in regulating lipid metabolism has attracted much attention. Our previous study found that octanoic acid levels were significantly reduced under Cd-exposed conditions in Hu Sheep, on the basis of which we investigated the protective effect of sodium octanoate, a derivative of octanoic acid, against Cd exposure in Hu Sheep in the present study. In this study, an animal model of Cd exposure in Hu Sheep was established. Comprehensive assessment of Cd-induced intestinal injury using hematoxylin and eosin (H&E) staining, immunostaining and carried out in-depth analyses combined with lipid metabolomics and transcriptomics. The results showed that Cd exposure triggered intestinal inflammation, barrier function damage and oxidative stress imbalance. Lipid metabolomics analysis showed that Cd exposure severely disrupted lipid metabolic processes, especially the glycerophospholipid metabolic pathway, suggesting that lipid metabolic disorders are closely related to intestinal injury. Notably, sodium octanoate could partially reverse the lipid metabolism abnormality by regulating the Adenosine 5'-monophosphate (AMP)-activated protein kinase (AMPK) signaling pathway, effectively alleviating the Cd toxicity, which provides a brand-new prevention and control strategy for Cd-induced intestinal injury in the livestock industry pollution-mediated disease.

Profile of miRNAs induced during sheep fat tail development and roles of four key miRNAs in proliferation and differentiation of sheep preadipocytes

Background

The fat tail of sheep is an adaptive trait that facilitates their adaptation to harsh natural environments. MicroRNAs (miRNAs) have been demonstrated to play crucial roles in the regulation of tail fat deposition.

Methods

In this study, miRNA-Seq was employed to investigate the expression profiles of miRNAs during different developmental stages of sheep fat tails and elucidate the functions of differentially expressed miRNAs (DE miRNAs).

Results

A total of 350 DE miRNAs were identified, among which 191, 60, 26, and 21 were significantly upregulated in tail fat tissues of fetal, lamb, hogget Altay sheep, and adult Xinjiang fine wool (XFW) sheep but downregulated in other stages. Furthermore, we predicted a set of candidate target genes (4,476) for the top 20 DE miRNAs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that they involve in several adipogenesis-related pathways. Subsequent investigations indicated that four DE miRNAs, miR-433-3p, miR-485-3p, miR-409-3p, and miR-495-3p, could suppress the expression of peroxisome proliferator-activated receptor gamma (PPARγ) and phosphoinositide-3-kinase regulatory subunit 3 (PIK3R3) and regulate the preadipocyte development in sheep. Meanwhile, the lipid metabolism-related genes, fatty acid-binding protein (FABP3), perilipin 1 (PLIN1), adiponectin C1Q and collagen domain containing (ADIPOQ), and lipoprotein lipase (LPL), were significantly downregulated (p < 0.01).

Conclusion

The expression patterns of miRNAs exhibited significant fluctuations during different development periods of the fat tail, and some of them may participate in the regulation of tail fat deposition by modulating the proliferation and differentiation of preadipocytes.

Comparative analysis of alternative splicing events in skeletal muscle of different sheep

This paper aims to investigate the relationship between genes with alternative splicing (AS) events and breed-specific differences in muscle development in two breeds of sheep. RNA-seq was utilized to identify genes with AS between Small-tailed Han sheep and Dorset sheep. The gene lists of differentially spliced genes were identified, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were conducted on these genes. In this study, 299 genes with 356 AS indicated significant differences between two diffrerent breeds. There are differences in 31 genes with 35 AS. Cassette, alt5' and alt3' exhibited the highest levels of enrichment across various significant levels. GO and KEGG enrichment analysis demonstrated a significant correlation between Wnt, TGF-beta, Notch and MAPK signaling pathways and the development of muscle in sheep. These findings indicate that genes with AS are linked to variations in muscle development in sheep. These results offer significant scientific and practical implications for improving the quality of sheep meat.

Comparison of alternative splicing (AS) events in adipose tissue of polled dorset versus small tail han sheep

Background

During the alternative splicing (AS), the exons of primary transcripts are spliced in various arrangements, resulting in structurally and functionally distinct mRNAs and proteins. This study aimed to examine genes with AS events from Small Tail Han sheep and Dorset sheep to explore the mechanism of adipose developments.

Methods

This study identified the genes with AS events in adipose tissues of two different sheep with next-generation sequencing. In this paper, genes with significantly different AS events were performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses.

Results

364 genes with 411 A S events showed significant differences in adipose tissues between the two breeds; 108 genes with 120 A S events were extremely significant differences between the two breeds. We identified several novel genes that are related with adipose growth and development. The results of KEGG and GO analysis indicated that oocyte meiosis, mitogen-activated protein kinase (Wnt), mitogen-activated protein kinase (MAPK) signaling pathway, etc. Were closely related to the adipose tissue developments.

Conclusions

This paper revealed that the genes with AS events are important for adipose tissues in sheep, exploring the mechanisms of AS events associated with adipose tissue developments in sheep of different breeds.

MicroRNA and circular RNA profiling in the deposited fat tissue of Sunite sheep

As the most typical deposited fat, tail fat is an important energy reservoir for sheep adapted to harsh environments and plays an important role as a raw material in daily life. However, the regulatory mechanisms of microRNA (miRNA) and circular RNA (circRNA) in tail fat development remain unclear. In this study, we characterized the miRNA and circRNA expression profiles in the tail fat of sheep at the ages of 6, 18, and 30 months. We identified 219 differentially expressed (DE) miRNAs (including 12 novel miRNAs), which exhibited a major tendency to be downregulated, and 198 DE circRNAs, which exhibited a tendency to be upregulated. Target gene prediction analysis was performed for the DE miRNAs. Functional analysis revealed that their target genes were mainly involved in cellular interactions, while the host genes of DE circRNAs were implicated in lipid and fatty acid metabolism. Subsequently, we established a competing endogenous RNA (ceRNA) network based on the negative regulatory relationship between miRNAs and target genes. The network revealed that upregulated miRNAs play a leading role in the development of tail fat. Finally, the ceRNA relationship network with oar-miR-27a_R-1 and oar-miR-29a as the core was validated, suggesting possible involvement of these interactions in tail fat development. In summary, DE miRNAs were negatively correlated with DE circRNAs during sheep tail fat development. The multiple ceRNA regulatory network dominated by upregulated DE miRNAs may play a key role in this developmental process.

Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes

BACKGROUND

After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.

REVIEW

Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.

CONCLUSIONS

Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.

Identification of potential miRNA-mRNA regulatory network and the key miRNAs in intramuscular and subcutaneous adipose

Intramuscular fat (IMF) is an important indicator for evaluating meat quality. Breeds with high IMF content are often accompanied by high subcutaneous fat (SCF), severely affecting the meat rate of pigs. Studying the mechanisms of miRNAs related to lipogenesis and lipid metabolism has important implications for pig breeding. We constructed two small RNA libraries from intramuscular and subcutaneous fat to evaluate the patterns of lipogenesis in Laiwu pig, a Chinese breed. A total of 286 differentially expressed miRNAs (DEmiRNAs), including 193 known miRNA and 93 novel miRNAs, were identified from two types of adipose. GO and KEGG enrichment analysis for DEmiRNAs showed that their target genes involved in many adipogenesis and lipid metabolism biological processes and signaling pathways, such as Wnt signaling pathway,MAPK signaling pathway, Hippo signaling pathway, PI3K-Akt signaling pathway, Melanogenesis, Signaling pathways regulating pluripotency of stem cells and so on. Then, we constructed a miRNA-mRNA interaction network to find out which miRNAs were the key miRNAs of regulation in Wnt signaling pathway. In this pathway, miR-331-3p, miR-339-5p, miR-874 and novel346_mature target PPARD, WNT10B, RSPO3, WNT2B. This study provides a theoretical basis for further understanding the post-transcriptional regulation mechanism of meat quality formation and predicting and treating diseases caused by ectopic fat.

Genome-Wide Analysis of microRNAs Identifies the Lipid Metabolism Pathway to Be a Defining Factor in Adipose Tissue From Different Sheep

microRNAs are a class of important non-coding RNAs, which can participate in the regulation of biological processes. In recent years, miRNA has been widely studied not only in humans and mice, but also in animal husbandry. However, compared with other livestock and poultry breeds, the study of miRNA in subcutaneous adipose tissue of sheep is not comprehensive. Transcriptome analysis of miRNAs in subcutaneous adipose tissue of Duolang sheep, and Small Tail Han sheep was performed using RNA-Seq technology. Differentially expressed miRNAs were screened between different breeds. Target genes were predicted, and then the joint analysis of candidate genes were conducted based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Finally, the RNA-Seq data were verified by real-time quantitative polymerase chain reaction (qRT-PCR). Herein, we identified 38 differentially expressed miRNAs (9 novel miRNAs and 29 known miRNAs). In addition, a total of 854 target genes were predicted by miRanda software. GO and KEGG pathway analysis demonstrated that regulation of lipolysis in adipocytes plays a key role in the deposition of subcutaneous adipose tissue in Duolang sheep and Small Tail Han sheep. The miRNAs might regulate fat deposits by regulating genes involved in regulation of lipolysis in adipocytes. Specifically, NC_ 040278.1_ 37602, oar-mir-493-3p, NC_ 040278.1_ 37521 and NC_ 040255.1_ 11627 might target PTGS2, AKT2, AKT3, and PIK3CA, respectively, and then play critical regulatory role. In conclusion, all the results provide a good idea for further revealing the mechanism of subcutaneous adipose tissue deposition and improving the meat production performance of sheep, and lay a foundation for promoting the development of animal husbandry.

MiR-128-1-5p regulates differentiation of ovine stromal vascular fraction by targeting the KLF11 5'-UTR

Fat content is an important index to evaluate the individual performance of livestock animals such as sheep for meat production purposes. Reducing the subcutaneous and visceral fat while increasing the intramuscular fat is a valuable goal to achieve for the meat production industry. Here, we investigated the effect of miR-128-1-5p on adipogenesis of subcutaneous fat by targeting 5'-UTR in KLF11, a rare mechanism where most miRNAs bind the 3'-UTR of mRNAs. A dual fluorescence reporter assay was conducted to validate the binding sites of miR-128-1-5p on 5'-UTR of KLF11 mRNA. Roles of miR-128-1-5p in KLF11 expression were measured through co-transfecting miRNA mimics with KLF11-expressing vectors (CDSs together with or without the 5'-UTR) into ovine stromal vascular fractions (SVF). Additionally, functional roles of miR-128-1-5p, and KLF11 in adipogenesis of ovine subcutaneous fat were investigated. Results showed that miR-128-1-5p targeted KLF11 5'-UTR, reduced the fluorescence activity of the dual fluorescent reporter vector, as well as KLF11 mRNA, and protein expression levels. During the differentiation of SVF, disturbing the expression of miR-128-1-5p and KLF11 changed the adipogenic differentiation of SVF as observed in the lipid formation, and adipogenic marker genes. This study indicates that miR-128-1-5p promotes the expression of lipogenic marker genes and the formation of lipid droplets by targeting KLF11 5'-UTR. Furthermore, overexpression, and inhibition of KLF11 indicate that KLF11 inhibited SVF differentiation. In summary, the 5'-UTR binding mechanism discovered in this study extends the understanding of miRNA functions. Key roles of miR-128-1-5p and KLF11 in the adipogenesis of sheep subcutaneous fat have potential values for improving the meat and/or fat ratio of domestic animals.

Emerging Roles of Non-Coding RNAs in the Feed Efficiency of Livestock Species

A global population of already more than seven billion people has led to an increased demand for food and water, and especially the demand for meat. Moreover, the cost of feed used in animal production has also increased dramatically, which requires animal breeders to find alternatives to reduce feed consumption. Understanding the biology underlying feed efficiency (FE) allows for a better selection of feed-efficient animals. Non-coding RNAs (ncRNAs), especially micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs), play important roles in the regulation of bio-logical processes and disease development. The functions of ncRNAs in the biology of FE have emerged as they participate in the regulation of many genes and pathways related to the major FE indicators, such as residual feed intake and feed conversion ratio. This review provides the state of the art studies related to the ncRNAs associated with FE in livestock species. The contribution of ncRNAs to FE in the liver, muscle, and adipose tissues were summarized. The research gap of the function of ncRNAs in key processes for improved FE, such as the nutrition, heat stress, and gut–brain axis, was examined. Finally, the potential uses of ncRNAs for the improvement of FE were discussed.

Signature selection analysis reveals candidate genes associated with production traits in Iranian sheep breeds

BACKGROUND

Sheep were among the first animals to be domesticated. They are raised all over the world and produce a major scale of animal-based protein for human consumption and play an important role in agricultural economy. Iran is one of the important locations for sheep genetic resources in the world. Here, we compared the Illumina Ovine SNP50 BeadChip data of three Iranian local breeds (Moghani, Afshari and Gezel), as a population that does not undergone artificial breeding programs as yet, and five other sheep breeds namely East Friesian white, East Friesian brown, Lacaune, DorsetHorn and Texel to detect genetic mechanisms underlying economical traits and daptation to harsh environments in sheep.

RESULTS

To identify genomic regions that have been targeted by positive selection, we used fixation index (Fst) and nucleotide diversity (Pi) statistics. Further analysis indicated candidate genes involved in different important traits such as; wool production included crimp of wool (PTPN3, NBEA and KRTAP20-2 genes), fiber diameter (PIK3R4 gene), hair follicle development (LHX2 gene), the growth and development of fiber (COL17A1 gene)), adaptation to hot arid environments (CORIN gene), adaptive in deficit water status (CPQ gene), heat stress (PLCB4, FAM107B, NBEA, PIK3C2B and USP43 genes) in sheep.

CONCLUSIONS

We detected several candidate genes related to wool production traits and adaptation to hot arid environments in sheep that can be applicable for inbreeding goals. Our findings not only include the results of previous researches, but also identify a number of novel candidate genes related to studied traits. However, more works will be essential to acknowledge phenotype- genotype relationships of the identified genes in our study.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differentially expressed transcripts in sheep tail fat

Background

Nowadays, both customers and producers prefer thin-tailed fat sheep. To effectively breed for this phenotype, it is important to identify candidate genes and uncover the genetic mechanism related to tail fat deposition in sheep. Accumulating evidence suggesting that post-transcriptional modification events of precursor-messenger RNA (pre-mRNA), including alternative splicing (AS) and alternative polyadenylation (APA), may regulate tail fat deposition in sheep. Differentially expressed transcripts (DETs) analysis is a way to identify candidate genes related to tail fat deposition. However, due to the technological limitation, post-transcriptional modification events in the tail fat of sheep and DETs between thin-tailed and fat-tailed sheep remains unclear.

Methods

In the present study, we applied pooled PacBio isoform sequencing (Iso-Seq) to generate transcriptomic data of tail fat tissue from six sheep (three thin-tailed sheep and three fat-tailed sheep). By comparing with reference genome, potential gene loci and novel transcripts were identified. Post-transcriptional modification events, including AS and APA, and lncRNA in sheep tail fat were uncovered using pooled Iso-Seq data. Combining Iso-Seq data with six RNA-sequencing (RNA-Seq) data, DETs between thin- and fat-tailed sheep were identified. Protein protein interaction (PPI) network, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were implemented to investigate the potential functions of DETs.

Results

In the present study, we revealed the transcriptomic complexity of the tail fat of sheep, result in 9,001 potential novel gene loci, 17,834 AS events, 5,791 APA events, and 3,764 lncRNAs. Combining Iso-Seq data with RNA-Seq data, we identified hundreds of DETs between thin- and fat-tailed sheep. Among them, 21 differentially expressed lncRNAs, such as ENSOART00020036299, ENSOART00020033641, ENSOART00020024562, ENSOART00020003848 and 9.53.1 may regulate tail fat deposition. Many novel transcripts were identified as DETs, including 15.527.13 (DGAT2), 13.624.23 (ACSS2), 11.689.28 (ACLY), 11.689.18 (ACLY), 11.689.14 (ACLY), 11.660.12 (ACLY), 22.289.6 (SCD), 22.289.3 (SCD) and 22.289.14 (SCD). Most of the identified DETs have been enriched in GO and KEGG pathways related to extracellular matrix (ECM). Our result revealed the transcriptome complexity and identified many candidate transcripts in tail fat, which could enhance the understanding of molecular mechanisms behind tail fat deposition.

Genome-wide transcriptome profiling uncovers differential miRNAs and lncRNAs in ovaries of Hu sheep at different developmental stages

Ovary development is an important determinant of the procreative capacity of female animals. Here, we performed genome-wide sequencing of long non-coding RNAs (lncRNAs) and mRNAs on ovaries of 1, 3 and 8 months old Hu sheep to assess their expression profiles and roles in ovarian development. We identified 37,309 lncRNAs, 45,404 messenger RNAs (mRNAs) and 330 novel micro RNAs (miRNAs) from the transcriptomic analysis. Six thousand, seven hundred and sixteen (6716) mRNAs and 1972 lncRNAs were significantly and differentially expressed in ovaries of 1 month and 3 months old Hu sheep (H1 vs H3). These mRNAs and target genes of lncRNAs were primarily enriched in the TGF-β and PI3K-Akt signalling pathways which are closely associated with ovarian follicular development and steroid hormone biosynthesis regulation. We identified MSTRG.162061.1, MSTRG.222844.7, MSTRG.335777.1, MSTRG.334059.16, MSTRG.188947.6 and MSTRG.24344.3 as vital genes in ovary development by regulating CTNNB1, CCNA2, CDK2, CDC20, CDK1 and EGFR expressions. A total of 2903 mRNAs and 636 lncRNAs were differentially expressed in 3 and 8 months old ovaries of Hu sheep (H3 vs H8); and were predominantly enriched in PI3K-Akt, progesterone-mediated oocyte maturation, estrogen metabolism, ovulation from the ovarian follicle and oogenesis pathways. These lncRNAs were also found to regulate FGF7, PRLR, PTK2, AMH and INHBA expressions during follicular development. Our result indicates the identified genes participate in the development of the final stages of follicles and ovary development in Hu sheep.

Quantitative proteomic analysis identified differentially expressed proteins with tail/rump fat deposition in Chinese thin- and fat-tailed lambs

Tail adipose as one of the important functional tissues can enhance hazardous environments tolerance for sheep. The objective of this study was to gain insight into the underlying development mechanisms of this trait. A quantitative analysis of protein abundance in ovine tail/rump adipose tissue was performed between Chinese local fat- (Kazakh, Hu and Lanzhou) and thin-tailed (Alpine Merino, Tibetan) sheep in the present study by using lable-free approach. Results showed that 3400 proteins were identified in the five breeds, and 804 were differentially expressed proteins, including 638 up regulated proteins and 83 down regulated proteins in the tail adipose tissues between fat- and thin-tailed sheep, and 8 clusters were distinguished for all the DEPs’ expression patterns. The differentially expressed proteins are mainly associated with metabolism pathways and peroxisome proliferator activated receptor signaling pathway. Furthermore, the proteomics results were validated by quantitative real-time PCR and Western Blot. Our research has also suggested that the up-regulated proteins ACSL1, HSD17β4, FABP4 in the tail adipose tissue might contribute to tail fat deposition by facilitating the proliferation of adipocytes and fat accumulation in tail/rump of sheep. Particularly, FABP4 highly expressed in the fat-tail will play an important role for tail fat deposition. Our study might provide a novel view to understanding fat accumulation in special parts of the body in sheep and other animals.

Identification and analysis of miRNAs in the normal and fatty liver from the Holstein dairy cow

MicroRNAs (miRNAs) are a class of non-coding short RNAs with ∼22 nts in length, which play important roles in the regulation of numerous biological processes in animals. In this study, two small RNA libraries from fatty (S01) and normal livers (S02) from Holstein Dairy Cow (HDC) were sequenced through deep sequencing. A total of 12,964,411 and 15,426,289 clean reads were obtained, representing 370 known and 182 novel miRNAs, respectively. The characterization, expression pattern, potential functions and target genes of these miRNAs were investigated. Analysis identified 66 upregulated and seven downregulated differentially expressed miRNAs (DIE-miRNAs). To verify the sequencing results, 10 DIE-miRNAs were selected for qRT-PCR, and the results were confirmed to be consistent with the miRNA sequencing. In addition, a total of 5,578 targets of the 73 DIE-miRNAs were predicted. GO analysis revealed that DIE-miRNAs targets are associated with cellular process, cell part and molecular transducer activity. KEGG pathway analysis showed that Arrhythmogenic right ventricular cardiomyopathy, Axon guidance, Ether lipid metabolism and Cocaine addiction were closely associated with liver metabolism. These findings will provide valuable information for further functional verification of miRNAs between normal and fatty liver, as might exploit new attractive miRNAs biomarkers for diseases detection in HDC.

Comparative DNA methylome analysis of estrus ewes reveals the complex regulatory pathways of sheep fecundity

BACKGROUND/AIMS

Sheep are important livestock with variant ovulation rate and fertility. Dorset sheep is a typical breed with low prolificacy, whereas Small Tail Han sheep with FecB mutation (HanBB) have hyperprolificacy. Our previous studies have revealed the gene expression difference between the ovaries from Dorset and HanBB sheep contributes to the difference of fecundity, however, what leads to these gene expression difference remains unclear. DNA methylation, an important epigenetic process, plays a crucial role in gene expression regulation.

METHODS

In the present study, we constructed a methylated DNA immunoprecipitation combined with high throughput sequencing (MeDIP-seq) strategy to investigate the differentially methylated genes between the Dorset and HanBB ovaries.

RESULTS

Our findings suggest the genes involved in immune response, branched-chain amino acid metabolism, cell growth and cell junction were differentially methylated in or around the gene body regions.

CONCLUSIONS

These findings provide prospective insights on the epigenetic basis of sheep fecundity.

Comparison Between the Gut Microbiota in Different Gastrointestinal Segments of Large-Tailed Han and Small-Tailed Han Sheep Breeds with High-Throughput Sequencing

Commensal microorganisms are essential to the normal development and function of many aspects of animal biology. However, the dynamic shift patterns of the microbiota of different gut segments in sheep and the correlation between fat type large-tailed phenotype and microbiota remain poorly unknown. This study therefore sought to assess the composition and distribution of the intestinal microbiome, and compared the difference of gut microbiota from different gastrointestinal segments within breeds and same intestinal sections between breeds. For these analyses, 16S rRNA V4 regions from 4 gut sections prepared from each of six individuals (3 from each breed) were sequenced to detect the microbiome composition in these samples. These analyses revealed the presence of 51,173 operational taxonomic units distributed across 24 phyla and 420 genera in these samples, with Firmicutes and Bacteroidetes being the most prevalent phyla of microbes present in these samples. Moreover, the bacterial composition showed distinct microbial communities in different gastrointestinal segments within breed, but showed similar and relative fixed bacterial abundance in the same intestinal segments from individuals of different breeds. We also found that only a few bacterial species (Lachnospiraceae, Akkermansia) were needed to distinguish between Small-tailed Han sheep (STH) and Large-tailed Han sheep (LTH) and their metabolic process maybe influence the fat type large-tailed phenotype formation in sheep. The functional profile analysis revealed that the environment information processing, genetic information processing, and metabolic pathways were enriched in all samples. The main functional roles of the gut microbiota were amino acid metabolism, replication and repair, carbohydrate metabolism, and membrane transport. Finally, our findings suggested that distinguished gut species between STH and LTH have relative fixed and the potential correlation is existing between the intestinal microorganisms and the large-tailed phenotype trait formation of sheep, which may offer clues for further investigation to detect the roles of intestinal microbiota in the metabolism and fat deposition in the tail of sheep.

miR-27a is an important adipogenesis regulator associated with differential lipid accumulation between intramuscular and subcutaneous adipose tissues of sheep

Micro ribonucleic acids (miRNAs) are crucial regulators for various biological processes. Despite important function in the proliferation and differentiation of preadipocytes, miRNA studies are limited in regional differences in adipogenesis. Here, we show that miR-27a plays an important role in regulating differential lipid accumulation between intramuscular (IM) and subcutaneous (SC) adipose tissues in sheep. Invivo, we observed that miR-27a expression in IM adipose tissue is more abundant than in SC adipose tissue. However, the expression of Peroxisome Proliferator-Activated Receptor Gamma (PPARG) and retinoid X receptor alpha (RXR alpha) in IM adipose tissue was significantly lower than that in SC adipose tissue. In the ovine preadipocyte differentiation model, we found that the expression of miR-27a was significantly decreased in differentiated ovine adipocytes. Overexpression of miR-27a significantly downregulated the expression of PPARG and RXR alpha and suppressed the accumulation of triglyceride but promoted the proliferation of ovine preadipocytes. Whereas, inhibition of miR-27a suppressed preadipocyte proliferation but enhanced PPARG and RXR alpha expression and lipid droplet formation. In addition, dual-luciferase activity assays showed that RXR alpha was a direct target of miR-27a. Thus, miR-27a enhances ovine preadipocytes proliferation and inhibits ovine preadipocytes differentiation through regulating the expression of target RXR alpha. Collectively, our study demonstrates the functional importance of miR-27a in ovine adipogenesis and provides novel insights into exploring regional differences in adipogenesis.

Identification and Profiling of Pituitary microRNAs of Sheep during Anestrus and Estrus Stages

MicroRNAs (miRNAs) are a class of small non-coding RNAs, molecules of 21 to 25 nucleotides in length, that regulate gene expression by binding to their target mRNA and play a significant role in animal development. The expression and role of miRNAs in regulating sheep estrus, however, remain elusive. Transcriptome analysis is helpful to understand the biological roles of miRNAs in the pituitary gland of sheep. A sheep's pituitary gland has a significant difference between estrus and anestrus states. Here, we investigate the expression profiles of sheep anterior pituitary microRNAs (miRNAs) in two states, estrus and anestrus, using Illumina HiSeq-technology. This study identified a total of 199 miRNAs and 25 differentially expressed miRNAs in the estrus and anestrus pituitary gland in sheep. Reverse transcription quantitative-PCR (RT-qPCR) analysis shows six differentially (p < 0.05) expressed miRNAs, that are miR-143, miR-199a, miR-181a, miR-200a, miR-218, and miR-221 in both estrus and anestrus states. miRNAs containing estrus-related terms and pathways regulation are enriched using enrichment analysis from gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Moreover, we also envisioned a miRNA-mRNA interaction network to understand the function of miRNAs involved in the pituitary gland regulatory network. In conclusion, miRNA expression profiles in sheep pituitary gland in the anestrus and estrus deliver a theoretical basis for the study of pituitary gland biology in sheep.

MicroRNA Expression Profile in Peripheral Blood Lymphocytes of Sheep Vaccinated with Nigeria 75/1 Peste Des Petits Ruminants Virus

Peste des petits ruminants (PPR) is one of the highly contagious transboundary viral diseases of small ruminants. Host microRNA (miRNA) expression patterns may change in response to virus infection, and it mainly works as a post-transcriptional moderator in gene expression and affects viral pathogenesis and replication. In this study, the change of miRNA expression profile in peripheral blood lymphocyte (PBMC) from sheep inoculated with PPR vaccine virus in vivo as well as primary sheep testicular (ST) cells inoculated with PPR vaccine virus in vitro were determined via deep sequencing technology. In PBMC cells, 373 and 115 differentially expressed miRNAs (DEmiRNAs) were identified 3 days and 5 days post inoculated (dpi), respectively. While, 575 DEmiRNAs were identified when comparing miRNA profiles on 5 dpi with 3 dpi. Some of the DEmiRNAs were found to change significantly via time-course during PPR vaccine virus inoculated. Similarly, in ST cells, 136 DEmiRNAs were identified at 3 dpi in comparison with mock-inoculation. A total of 12 DEmiRNAs were validated by real-time quantitative PCR (RT-qPCR). The oar-miR-150, oar-miR-370-3p and oar-miR-411b-3p were found common differentially expressed in both PPR vaccine virus-inoculated PBMC cells and ST cells. Targets prediction and functional analysis of the DEmiRNAs uncovered mainly gathering in antigen processing and presentation pathways, protein processing in endoplasmic reticulum pathways and cell adhesion molecules pathways. Our study supplies information about the DEmiRNAs in PPR vaccine virus-inoculated PBMC cells and ST cells, and provides clues for further understanding the function of miRNAs in PPR vaccine virus replication.

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.

Identification and Expression Analysis of Long Noncoding RNAs in Fat-Tail of Sheep Breeds

Emerging evidence suggests that long non-coding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. However, most studies have been focused on a few established model organisms and little is known about lncRNAs in fat-tail development in sheep. Here, the first profile of lncRNA in sheep fat-tail along with their possible roles in fat deposition were investigated, based on a comparative transcriptome analysis between fat-tailed (Lori-Bakhtiari) and thin-tailed (Zel) Iranian sheep breeds. Among all identified lncRNAs candidates, 358 and 66 transcripts were considered novel intergenic (lincRNAs) and novel intronic (ilncRNAs) corresponding to 302 and 58 gene loci, respectively. Our results indicated that a low percentage of the novel lncRNAs were conserved. Also, synteny analysis identified 168 novel lincRNAs with the same syntenic region in human, bovine and chicken. Only seven lncRNAs were identified as differentially expressed genes between fat and thin tailed breeds. Q-RT-PCR results were consistent with the RNA-Seq data and validated the findings. Target prediction analysis revealed that the novel lncRNAs may act in cis or trans and regulate the expression of genes that are involved in the lipid metabolism. A gene regulatory network including lncRNA-mRNA interactions were constructed and three significant modules were found, with genes relevant to lipid metabolism, insulin and calcium signaling pathway. Moreover, integrated analysis with AnimalQTLdb database further suggested six lincRNAs and one ilncRNAs as candidates of sheep fat-tail development. Our results highlighted the putative contributions of lncRNAs in regulating expression of genes associated with fat-tail development in sheep.

Expression analysis of lung miRNAs responding to ovine VM virus infection by RNA-seq

BACKGROUND

MicroRNAs (miRNAs) are short endogenous, single-stranded, noncoding small RNA molecules of approximately 22 nucleotides in length. They regulate gene expression posttranscriptionally by silencing mRNA expression, thus orchestrating many physiological processes. The Small Ruminant Lentiviruses (SRLV) group includes the Visna Maedi Virus (VMV) and Caprine Arthritis Encephalitis (CAEV) viruses, which cause a disease in sheep and goats characterized by pneumonia, mastitis, arthritis and encephalitis. Their main target cells are from the monocyte/macrophage lineage. To date, there are no studies on the role of miRNAs in this viral disease.

RESULTS

Using RNA-seq technology and bioinformatics analysis, the expression levels of miRNAs during different clinical stages of infection were studied. A total of 212 miRNAs were identified, of which 46 were conserved sequences in other species but found for the first time in sheep, and 12 were completely novel. Differential expression analysis comparing the uninfected and seropositive groups showed changes in several miRNAs; however, no significant differences were detected between seropositive asymptomatic and diseased sheep. The robust increase in the expression level of oar-miR-21 is consistent with its increased expression in other viral diseases. Furthermore, the target prediction of the dysregulated miRNAs revealed that they control genes involved in proliferation-related signalling pathways, such as the PI3K-Akt, AMPK and ErbB pathways.

CONCLUSIONS

To the best of our knowledge, this is the first study reporting miRNA profiling in sheep in response to SRLV infection. The known functions of oar-miR-21 as a regulator of inflammation and proliferation appear to be a possible cause of the lesions caused in the sheep's lungs. This miRNA could be an indicator for the severity of the lung lesions, or a putative target for therapeutic intervention.

RumimiR: a detailed microRNA database focused on ruminant species

The ever-increasing use of next-generation sequencing technologies to explore the genome has generated large quantities of data in recent years. Numerous publications have described several thousand sequences of microRNAs, all species included. A new database (RumimiR) has been created from the literature to provide a detailed description of microRNAs for three ruminant species: cattle, goats and sheep. To date, 2887, 2733 and 5095 unique microRNAs from bovine, caprine and ovine species, respectively, are included. In addition to the most recent reference genomic position and sequence of each microRNA, this database contains details about the animals, tissue origins and experimental conditions mentioned in the publications. Identity to human or mouse microRNA is also indicated. The RumimiR database allows data filtering by selecting microRNAs on the basis of defined criteria such as animal status or tissue origin. For ruminant studies, RumimiR supplements the widely used miRBase database, by using complementary criteria to allow browsing and filtering, and integrates all newly described published sequences. The principal goal of this database is to provide easy access to all the ruminant microRNAs described in the literature.

Porcine reproductive and respiratory syndrome virus counteracts type I interferon-induced early antiviral state by interfering IRF7 activity

Porcine reproductive and respiratory syndrome (PRRS) is an economically important disease with a significant impact on the pig industry. It is caused by PRRS virus (PRRSV), which predominantly infects and replicates in porcine pulmonary alveolar macrophages (PAMs). We pretreated PAMs with porcine interferon (IFN)-α to induce an antiviral state within the cells and subsequently infected them with highly pathogenic PRRSV. Changes in global gene expression in IFN-α-pretreated PAMs in response to PRRSV infection were determined by RNA-sequence analysis and confirmed by real-time PCR. We found that IRF7 and other antiviral interferon stimulating genes (ISG)s were suppressed by PRRSV infection. Further studies demonstrated that PRRSV could down-regulate the expression of IRF7 by the non-structure protein 7 (nsp7). In conclusion, PRRSV infection had a strong immunosuppressive effect of IFN. PRRSV nsp7 inhibits the expression of IRF7, thereby down-regulating the expression of IFN and downstream ISGs and facilitated the virus to replicate.

Early adipose deposits in sheep: comparative analysis of the perirenal fat transcriptome of Assaf and Churra suckling lambs

Adipose deposits influence the quality of ruminant carcasses, and in suckling lambs, internal types of adipose deposits represent a notable proportion of total fat. The aim of this study was to perform a comparative analysis of the perirenal fat transcriptomes of suckling lambs from two breeds with different growth and carcass characteristics. The perirenal fat tissue from 14 suckling lambs (Assaf, n = 8; Churra, n = 6) was used for the RNA-seq analysis. The functional enrichment analysis of the 670 highly expressed genes (>150 fragments per kilobase of exon per million fragments mapped) in the perirenal fat transcriptome of both breeds revealed that the majority of these genes were involved in energy processes. The expression of the UCP1 gene, a classical biomarker of brown fat, and the presence of multilocular adipocytes in the two breeds supported the presence of brown fat at the transition stage towards white fat tissue. The differential expression analysis performed identified 373 differentially expressed genes (DEGs) between the two compared breeds. Brown/white fat gene biomarkers were not included in the list of DEGs. In Assaf lambs, DEGs were enriched in Gene Ontology (GO) biological processes related to fatty-acid oxidation, whereas in Churra lambs, the majority of the significantly enriched GO terms were related to cholesterol synthesis, which suggests that upregulated DEGs in Assaf lambs are implicated in fat burning, whereas the Churra upregulated DEGs are linked to fat accumulation. These results can help to increase knowledge of the genes controlling early fat deposition in ruminants and shed light on fundamental aspects of adipose tissue growth.

Ovarian transcriptomic analysis reveals the alternative splicing events associated with fecundity in different sheep breeds

Alternative splicing (AS) is one of the most common mechanisms that accounts for the greater macromolecular and cellular complexity of higher eukaryotic organisms. This study focused on the splicing events in the ovaries of different sheep breeds, namely the Han and Dorset breeds. Of the groups studied, cassette splicing events accounted for the maximum number of the AS events with significant differences, whereas the splicing events that were mutually exclusive with introns accounted for the smallest proportion of splicing events. Greater than 1000 AS events with significant differences were identified between the Han BB and Dorset sheep. The number of AS events with significant differences between Han ++ and Dorset sheep, however, was fewer than that of the comparison of Han BB and Dorset sheep. Seven randomly selected genes with AS events were detected in this study and were validated by reverse transcription-polymerase chain reaction (RT-PCR). Furthermore, there are many genes which were common to the two genotype groups (Han BB and Dorset sheep, as well as the Han ++ and Dorset sheep). In addition, genes detected in the present study were involved in different pathways, including the pathways related with fertility or fecundity. The present study could provide the detailed understanding on the mechanisms of alternative splicing events associated with fecundity in different sheep breeds.

Comparative Transcriptome Profiling of mRNA and lncRNA Related to Tail Adipose Tissues of Sheep

The Lanzhou Fat-Tail sheep (LFTS, long fat-tailed sheep) is an endangered sheep breed in China with a fat tail compared to the traditional local varieties, Small Tail Han sheep (STHS, thin-tailed sheep) with a small tail, and Tibetan sheep (TS, short thin-tailed sheep) with a little tail. However, little is known regarding how tail fat deposition is regulated by long noncoding RNA (lncRNA). To evaluate the lncRNA and mRNA associated with tail fat deposition and development among these breeds, high-throughput RNA sequencing of three individuals each of LFTS, STHS, and TS were performed and analyzed in this study. RNA sequencing data from these three groups revealed 10 differentially expressed genes (DEGs) and 37 differentially expressed lncRNAs between the LFTS and STHS groups, 390 DEGs and 59 differentially expressed lncRNAs between the LFTS and TS groups, and 80 DEGs and 16 differentially expressed lncRNAs between the STHS and TS groups (p-value < 0.05 and fold change ≥ 2), respectively. Gene Ontology and pathway analysis of DEGs and target genes of differentially expressed lncRNAs revealed enrichment in fatty acid metabolism and fatty acid elongation-related pathways that contribute to fat deposition. Subsequently, the expression of 14 DEGs and 6 differentially expressed lncRNAs was validated by quantitative real-time PCR. Finally, two co-expression networks of differentially expressed mRNA and lncRNAs were constructed. The results suggested that some differentially expressed lncRNAs (TCONS_00372767, TCONS_00171926, TCONS_00054953, and TCONS_00373007) may play crucial roles as core lncRNAs in tail fat deposition processes. In summary, the present study extends the sheep tail fat lncRNA database and these differentially expressed mRNA and lncRNAs may provide novel candidate regulators for future genetic and molecular studies on tail fat deposition of sheep.

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.

An Integrated Analysis of miRNAs and Methylated Genes Encoding mRNAs and lncRNAs in Sheep Breeds with Different Fecundity

In our previous study, we investigated the regulatory relationship between lncRNAs, miRNA, and mRNAs in an effort to shed light onto the regulatory mechanisms involved in sheep fecundity. As an extension of this study, here, we aimed to identify potential regulators of sheep fecundity using a genome-wide analysis of miRNAs and the methylated genes encoding mRNAs and lncRNAs in the ovaries of Dorset sheep (low fecundity) and Small Tail Han ewes (high fecundity) with the genotype BB (Han BB) and the genotype ++ (Han ++) by performing RNA-Seq and MeDIP-Seq analyses. Methylated coding-non-coding gene co-expression networks for Han and Dorset sheep were constructed using the methylated genes encoding the differentially expressed mRNAs and lncRNAs identified in this study. In the Han BB vs. Dorset comparison, the lncRNAs TTC26 and MYH15 had the largest degree. Similarly, the lncRNA NYAP1 had the largest degree in the Han ++ vs. Dorset comparison. None of the methylated genes encoding lncRNAs were co-expressed with the methylated genes encoding mRNAs in the Han BB vs. Han ++ comparison. The methylated genes encoding lncRNAs identified here may play a vital regulatory role in sheep breeding. Our results suggest that miRNAs might play a key role in sheep prolificacy by regulating target genes related to thyroid hormone synthesis, and methylated genes encoding lncRNAs associated with tight junctions might contribute to the high breeding rate in Han sheep. These findings may contribute to a deeper understanding of sheep prolificacy.

Differential regulation of mRNAs and lncRNAs related to lipid metabolism in two pig breeds

Long non-coding RNAs (lncRNAs) can regulate lipid metabolism and adipogenesis. However, there is little research on the role of lncRNAs in fat deposition in pig. In this study, RNA-seq technology was used to analyze the gene expression profiles of subcutaneous adipose tissue in Laiwu (LW) and Large White (LY) pigs. Then, key lncRNAs and genes associated with lipid metabolism and adipogenic differentiation were identified. Fifty four lncRNAs and 482 known mRNAs were differentially expressed in the two pig breeds. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses revealed that differentially expressed genes and the target genes of differentially expressed lncRNAs were significantly enriched in PPAR signaling pathway and biological processes including fat cell differentiation and fatty acid metabolism. Key lncRNAs might regulate adipogenic differentiation and fatty acid metabolism by regulating genes involved in above signaling pathway and biological processes. Specifically, XLOC_014379, XLOC_011279, XLOC_064871, XLOC_019518 and XLOC_013639 might target SCD, LPIN1, TRIB3, EGR2 and FABP3, respectively, and then play critical regulatory role. These results are useful for understanding fat deposition in pig, breeding livestock with high quality meat, and preventing and treating lipid metabolic disease.

Global transcriptome analysis identifies differentially expressed genes related to lipid metabolism in Wagyu and Holstein cattle

Fat deposition of beef cattle varies between breeds. However, the regulation mechanism is still not elucidated completely at molecular level. In the present study, we comparatively analyzed transcriptome of subcutaneous adipose tissue between Wagyu and Holstein cattle with a significant difference in fat deposition to identify key genes associated with fat metabolism and adipogenesis by high-throughput RNA-seq technology. A total of 59,149,852 and 69,947,982 high quality reads were generated, respectively. With further analysis, 662 differentially expressed genes were identified. Gene Ontology and KEGG pathway analysis revealed that many differentially expressed genes were enriched in several biological processes and pathways relevant to adipogenesis and lipid metabolism, in which PPAR and fatty acid metabolism signaling pathways with related genes such as PPARγ, PLIN2 and ELOVL6 et al. play a critical role. Protein-protein interaction network analysis showed EGR1, FOS, SERPINE1, AGT, MMP2 may have great impact on adipocyte differentiation and adipogenesis. Moreover, potential alternative splicing events and single nucleotide polymorphisms (SNPs) were also identified. In summary, we comprehensively analyzed and discussed the transcriptome of subcutaneous adipose tissue of Wagyu and Holstein cattle, which might provide a theoretical basis for better understanding molecular mechanism of fat metabolism and deposition in beef cattle.

Identification and characterization of differentially expressed miRNAs in subcutaneous adipose between Wagyu and Holstein cattle

MicroRNAs (miRNAs) are important post-transcriptional regulators involved in animal adipogenesis, however, their roles in bovine fat deposition remain poorly understood. In the present study, we conducted a comparative RNA sequencing to identify the key miRNAs involved in beef lipid accumulation by comparing the backfat small RNA samples between Wagyu (high intramuscular fat) and Holstein (moderate intramuscular fat) cattle. Fifteen miRNAs such as bta-miR-142-3p, bta-miR-379, bta-miR-196a, bta-miR-196b, bta-miR-30f and bta-miR-2887 were identified to have a higher expression level in Wagyu cattle compared with Holstein, whereas bta-miR-320a, bta-miR-874 and bta-miR-1247-3p had a lower expression level in Wagyu. Furthermore, a total of 1345 potential target genes of differentially expressed miRNAs were predicted using bioinformatics tools, in which PPARα and RXRα were known to play a critical role in adipocyte differentiation and lipid metabolism. In conclusion, the present study constructed a high-throughput RNA sequencing screen and successfully identified miRNAs such as bta-miR-874, bta-miR-320a and bta-miR-196b which may affect beef fat deposition. The present findings may provide a theoretical foundation for the utilization of beef cattle germplasm resources.

Co-expression analysis and identification of fecundity-related long non-coding RNAs in sheep ovaries

Small Tail Han sheep, including the FecB^BFecB^B (Han BB) and FecB⁺ FecB⁺ (Han++) genotypes, and Dorset sheep exhibit different fecundities. To identify novel long non-coding RNAs (lncRNAs) associated with sheep fecundity to better understand their molecular mechanisms, a genome-wide analysis of mRNAs and lncRNAs from Han BB, Han++ and Dorset sheep was performed. After the identification of differentially expressed mRNAs and lncRNAs, 16 significant modules were explored by using weighted gene coexpression network analysis (WGCNA) followed by functional enrichment analysis of the genes and lncRNAs in significant modules. Among these selected modules, the yellow and brown modules were significantly related to sheep fecundity. lncRNAs (e.g., NR0B1, XLOC_041882, and MYH15) in the yellow module were mainly involved in the TGF-β signalling pathway, and NYAP1 and BCORL1 were significantly associated with the oxytocin signalling pathway, which regulates several genes in the coexpression network of the brown module. Overall, we identified several gene modules associated with sheep fecundity, as well as networks consisting of hub genes and lncRNAs that may contribute to sheep prolificacy by regulating the target mRNAs related to the TGF-β and oxytocin signalling pathways. This study provides an alternative strategy for the identification of potential candidate regulatory lncRNAs.

Genome-wide analysis of miRNAs in the ovaries of Jining Grey and Laiwu Black goats to explore the regulation of fecundity

Goat fecundity is important for agriculture and varies depending on the genetic background of the goat. Two excellent domestic breeds in China, the Jining Grey and Laiwu Black goats, have different fecundity and prolificacies. To explore the potential miRNAs that regulate the expression of the genes involved in these prolific differences and to potentially discover new miRNAs, we performed a genome-wide analysis of the miRNAs in the ovaries from these two goats using RNA-Seq technology. Thirty miRNAs were differentially expressed between the Jining Grey and Laiwu Black goats. Gene Ontology and KEGG pathway analyses revealed that the target genes of the differentially expressed miRNAs were significantly enriched in several biological processes and pathways. A protein-protein interaction analysis indicated that the miRNAs and their target genes were related to the reproduction complex regulation network. The differential miRNA expression profiles found in the ovaries between the two distinctive breeds of goats studied here provide a unique resource for addressing fecundity differences in goats.

Ovarian transcriptomic study reveals the differential regulation of miRNAs and lncRNAs related to fecundity in different sheep

miRNAs and lncRNAs, which represent one of the most highly expressed classes of ncRNAs in development, are attracting increasing interest. A variety of regulators is considered to be implicated in sheep species with different fecundity. However, interactions between miRNAs and lncRNAs and changes in the expression of regulatory lncRNAs in sheep fecundity have not yet been reported. To characterize the important roles of miRNAs and lncRNAs and elucidate their regulating networks in sheep prolificacy, a genome-wide analysis of miRNAs and lncRNAs from Small Tail Han sheep of genotypes FecB^BFecB^B (Han BB) and FecB⁺ FecB⁺ (Han++) and from Dorset sheep (Dorset) was performed. An integrated analysis of miRNAs and lncRNAs was performed to study the regulatory function of miRNAs and lncRNAs in fecundity, revealing significantly correlated patterns of expression. Dramatic changes of miRNAs and lncRNAs suggest their critical roles in sheep fecundity. In conclusion, this is the first study performing thorough investigations of regulatory relationships among lncRNAs, miRNA and mRNAs, which will provide a novel view of the regulatory mechanisms involved in sheep fecundity. These results may provide further insight into sheep fecundity and help us to improve sheep prolificacy.

Ovarian proteomic study reveals the possible molecular mechanism for hyperprolificacy of Small Tail Han sheep

Small Tail Han sheep is a widely bred farm animal in China which has attracted lots of attention due to their high prolificacy and year-round estrus. However, the molecular mechanism of its fecundity remains unrevealed. The FecB gene polymorphism has been found to be associated with the ovulation rate and litter size of sheep. In the present study, we constructed an iTRAQ-based quantitative proteomics analysis to compare the ovarian proteomes of FecB+FecB+ genotype Small Tail Han sheep ewes (Han ++), FecB^BFecB^B Han ewes (Han BB) and Dorset ewes (Dorset). Hundreds of differentially expressed proteins between each two groups were identified; GO and KEGG pathway analysis indicated that the expressions of those proteins involved in ribosome assembly, protein translation and mTOR pathway between Dorset and both Han groups were highly different. Between Han ++ and Han BB groups, higher level of protein expressions were related to mitochondrial oxidation functions such as oxidoreductase activity, cytochrome-c oxidase activity and electron carrier activity. This was identified in Han BB group, which may contribute to the elevated ovulation rate of Han BB ewes. In conclusion, our work provided a prospective understanding of the molecular mechanism for high prolificacy of Small Tail Han sheep.

Genome-wide transcriptome analysis in the ovaries of two goats identifies differentially expressed genes related to fecundity

The goats are widely kept as livestock throughout the world. Two excellent domestic breeds in China, the Laiwu Black and Jining Grey goats, have different fecundities and prolificacies. Although the goat genome sequences have been resolved recently, little is known about the gene regulations at the transcriptional level in goat. To understand the molecular and genetic mechanisms related to the fecundities and prolificacies, we performed genome-wide sequencing of the mRNAs from two breeds of goat using the next-generation RNA-Seq technology and used functional annotation to identify pathways of interest. Digital gene expression analysis showed 338 genes were up-regulated in the Jining Grey goats and 404 were up-regulated in the Laiwu Black goats. Quantitative real-time PCR verified the reliability of the RNA-Seq data. This study suggests that multiple genes responsible for various biological functions and signaling pathways are differentially expressed in the two different goat breeds, and these genes might be involved in the regulation of goat fecundity and prolificacy. Taken together, our study provides insight into the transcriptional regulation in the ovaries of 2 species of goats that might serve as a key resource for understanding goat fecundity, prolificacy and genetic diversity between species.