Genome-wide transcriptome analysis of mRNAs and microRNAs in Dorset and Small Tail Han sheep to explore the regulation of fecundity

DIA-Based Proteomic Analysis Reveals MYOZ2 as a Key Protein Affecting Muscle Growth and Development in Hybrid Sheep

Hybridization of livestock can be used to improve varieties, and different hybrid combinations produce unique breeding effects. In this study, male Southdown and Suffolk sheep were selected to hybridize with female Hu sheep to explore the effects of male parentage on muscle growth and the development of offspring. Using data-independent acquisition technology, we identified 119, 187, and 26 differentially abundant proteins (DAPs) between Hu × Hu (HH) versus Southdown × Hu (NH), HH versus Suffolk × Hu (SH), and NH versus SH crosses. Two DAPs, MYOZ2 and MYOM3, were common to the three hybrid groups and were mainly enriched in muscle growth and development-related pathways. At the myoblast proliferation stage, MYOZ2 expression decreased cell viability and inhibited proliferation. At the myoblast differentiation stage, MYOZ2 expression promoted myoblast fusion and enhanced the level of cell fusion. These findings provide new insights into the key proteins and metabolic pathways involved in the effect of male parentage on muscle growth and the development of hybrid offspring in sheep.

Identification and analysis of differentially expressed microRNAs in endometrium to explore the regulation of sheep fecundity

BACKGROUND

MicroRNAs (miRNAs) play an important regulatory role in mammalian reproduction. Currently, most studies are primarily concentrated on ovarian miRNAs, ignoring the influence of endometrial miRNAs on the fecundity of female sheep. To uncover potential regulators of sheep fecundity, RNA-seq was used to comparatively analyze miRNA expression profiles of endometrium between high prolificacy sheep (HP, litter size = 3) and low prolificacy sheep (LP, litter size = 1) with FecB genotype.

RESULTS

Firstly, genomic features of miRNAs from endometrium were analyzed. Furthermore, 58 differentially expressed (DE) miRNAs were found in the endometrium of Hu sheep with different litter size. A co-expression network of DE miRNAs and target genes has been constructed, and hub genes related litter size are included, such as DE miRNA unconservative_NC_019472.2_1229533 and unconservative_NC_019481.2_1637827 target to estrogen receptor α (ESR1) and unconservative_NC_019481.2_1637827 targets to transcription factor 7 (TCF7). Moreover, functional annotation analysis showed that the target genes (NRCAM and NEGR1) of the DE miRNAs were significantly enriched in cell adhesion molecules (CAMs) signaling pathway, which was related to uterine receptivity.

CONCLUSION

Taken together, this study provides a new valuable resource for understanding the molecular mechanisms underlying Hu sheep prolificacy.

Genome wide identification and characterization of fertility associated novel CircRNAs as ceRNA reveal their regulatory roles in sheep fecundity

Reproductive traits play a vital role in determining the production efficiency of sheep. Maximizing the production is of paramount importance for breeders worldwide due to the growing population. Circular RNAs (circRNAs) act as miRNA sponges by absorbing miRNA activity through miRNA response elements (MREs) and participate in ceRNA regulatory networks (ceRNETs) to regulate mRNA expression. Despite of extensive research on role of circRNAs as miRNA sponges in various species, their specific regulatory roles and mechanism in sheep ovarian tissue are still not well understood. In this study, we performed whole genome sequencing of circRNAs, miRNA and mRNA employing bioinformatic techniques on ovine tissues of two contrasting sheep breeds "Small tail Han (X_LC) and Dolang sheep (D_LC)", which results into identification of 9,878 circRNAs with a total length of 23,522,667 nt and an average length of 2,381.32 nt. Among them, 44 differentially expressed circRNAs (DECs) were identified. Moreover, correlation between miRNA-mRNA and lncRNA-miRNA provided us with to prediction of miRNA binding sites on nine differentially expressed circRNAs and 165 differentially expressed mRNAs using miRanda. miRNA-mRNA and lncRNA-miRNA pairs with negative correlation were selected to determine the ceRNA score along with positively correlated pairs from lncRNA and mRNA network. Integration of ceRNA score and positively correlated pairs exhibit a significant ternary relationship among circRNAs-miRNA-mRNA demonestrated by ceRNA, comprising of 50 regulatory pairs sharring common nodes and predicted potential differentially expressed circRNAs-miRNAs-mRNAs regulatory axis. Based on functional enrichment analysis shortlisted key ceRNA regulatory pairs associated with reproduction including circRNA_3257-novel579_mature-EPHA3, circRNA_8396-novel130_mature-LOC101102473, circRNA_4140- novel34_mature > novel661_mature-KCNK9, and circRNA_8312-novel339_mature-LOC101110545. Furthermore, expression profiling, functional enrichments and qRT-PCR analysis of key target genes infer their implication in reproduction and metabolism. ceRNA target mRNAs evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations following genomic organizations will provide new insights underlying molecular mechanisms of reproduction, and establish a solid foundation for future research. Graphical abstract summarizing the scheme of study.

A meta-analysis of mRNA expression profiling studies in sheep with different FecB genotypes

The FecB mutation in the sheep BMPRIB is strongly correlated with high ovulation traits but its mechanism remains unclear. This study explored differentially expressed genes (DEGs) and their associated molecular mechanisms that may be involved in FecB mutation-induced high ovulation from the perspective of the hypothalamic-pituitary-gonadal (HPG) axis by conducting a systematic review and meta-analysis. The PubMed, EMBASE, CNKI, WanFang, and CBM databases were searched for eligible articles published before August 2022, focusing on mRNA sequencing of different tissues in the HPG axis in sheep with different FecB genotypes. A total of 6555 DEGs were identified from the analysis of six published articles and experimental results from our laboratory. The DEGs were screened by vote-counting rank and robust rank aggregation. Among these, in the follicular phase, FKBP5, CDCA7 and CRABP1 were upregulated in the hypothalamus. INSM2 was upregulated, while LDB3 was downregulated in the pituitary. CLU, SERPINA14, PENK, INHA and STAR were upregulated, while FERMT2 and NPY1R were downregulated in the ovary. On the HPG axis, TAC1 was upregulated and NPNT was downregulated. Many DEGs were found in sheep with different FecB genotypes. The genes FKBP5, CDCA7, CRABP1, INSM2, LDB3, CLU, SERPINA14, PENK, INHA, STAR, FERMT2, NPY1R, TAC1 and NPNT, may be associated with FecB mutation-induced high ovulation in different tissues. These candidate genes will further improve the mechanism of multiple fertility traits induced by the FecB mutation from the perspective of the HPG axis.

Impact of FecB Mutation on Ovarian DNA Methylome in Small-Tail Han Sheep

Booroola fecundity (FecB) gene, a mutant of bone morphogenetic protein 1B (BMPR-1B) that was discovered in Booroola Merino, was the first prolificacy gene identified in sheep related to increased ovulation rate and litter size. The mechanism of FecB impact on reproduction is unclear.

METHODS

In this study, adult Han ewes with homozygous FecB(B)/FecB(B) mutations (Han BB group) and ewes with FecB(+)/FecB(+) wildtype (Han ++ group) were selected. Methylated DNA immunoprecipitation and high-throughput sequencing (MeDIP-seq) was used to identify differences in methylated genes in ovary tissue.

RESULTS

We examined differences in DNA methylation patterns between HanBB and Han ++ sheep. In both sheep, methylated reads were mainly distributed at the gene body regions, CpG islands and introns. The differentially methylated genes were enriched in neurotrophy in signaling pathway, Gonadotropin Releasing Hormone (GnRH) signaling pathway, Wnt signaling pathway, oocyte meiosis, vascular endothelial growth factor (VEGF) signaling pathway, etc. Differentially-methylated genes were co-analyzed with differentially-expressed mRNAs. Several genes which could be associated with female reproduction were identified, such as FOXP3 (forkhead box P3), TMEFF2 (Transmembrane Protein with EGF Like and Two Follistatin Like Domains 2) and ADAT2 (Adenosine Deaminase TRNA Specific 2).

CONCLUSIONS

We constructed a MeDIP-seq based methylomic study to investigate the ovarian DNA methylation differences between Small-Tail Han sheep with homozygous FecB mutant and wildtype, and successfully identified FecB gene-associated differentially-methylated genes. This study has provided information with which to understand the mechanisms of FecB gene-induced hyperprolificacy in sheep.

Integrative analysis of renal microRNA and mRNA to identify hub genes and pivotal pathways associated with cyclosporine-induced acute kidney injury in mice

Cyclosporine (CsA) is an immunosuppressive agent that often causes acute kidney injury (AKI) in children. The specific mechanisms underlying CsA-induced AKI are currently unknown. This study used an integrated network analysis of microRNA (miRNA) and mRNA expression profiles, biochemical and pathological analyses to further investigate these potential mechanisms of CsA-induced AKI. Small RNA sequence analysis identified 25 differentially expressed miRNAs, RNA sequencing analysis identified 4,109 differentially expressed mRNAs. We obtained a total of 4,367 target genes from the 25 differentially expressed miRNAs based on three algorithms, including the Mirdb, Mirtarbase, and TargetScan. 971 target genes overlapped between the 4,367 target genes and 4,109 differentially expressed mRNAs were identified for further bioinformatics analysis. Finally, 30 hub genes and two main modules were recognized. Functional enrichment analysis of 30 hub genes indicated that inflammation and epithelial-mesenchymal transition (EMT) related genes were mainly concentrated together. Pathway analysis revealed that the PI3K-Akt signaling pathway plays an integral role in CsA-induced AKI. Network analysis identified 3 important miRNAs, mmu-miR-17b-5p, mmu-miR-19b-3p, and mmu-mir-423-5p that may further promote the development of inflammatory responses and EMT by mediating a complex network of factors. Our research provides a clearer understanding the molecular mechanism of this specific drug-induced AKI by CsA use, which is useful for discovering potential targets for gene therapies, and drug development in CsA-induced AKI.

Integrative Analysis of miRNAs Involved in Fat Deposition in Different Pig Breeds

BACKGROUND

miRNAs are a set of small, noncoding RNAs that bind to partially complementary sequences on target mRNAs. This leads to the post-transcriptional regulation of gene expression. Many studies have shown that microRNAs play critical roles in adipose cell differentiation and fat metabolism. The aim of this study was to explore the regulatory functions of miRNAs in fat deposition for the prevention and therapy of lipid metabolism-related diseases.

METHODS

The significant differences in the fat deposition of Laiwu (LW) pigs and Large White (LY) pigs were studied. To investigate the genetic relationships of miRNAs that regulate fat deposition, we performed a genome-wide analysis of miRNAs derived from subcutaneous adipose tissue of LW and LY pigs using RNA-seq.

RESULTS

There were 39 known miRNAs and 56 novel miRNAs significantly differential expressed between the two breeds of pigs. In the analysis of the Gene Ontology and KEGG pathways, predicted targets of these differentially expressed miRNAs were involved in several fat-associated pathways, such as the peroxisome proliferator-activated receptor (PPAR), mitogen-activated protein kinases (MAPK) and Wnt signaling pathways. In addition, ssc-miR-133a-3p, ssc-miR-486 and ssc-miR-1 each had a great impact on the development of porcine subcutaneous fat through the PPAR signaling pathway.

CONCLUSIONS

We explored the role of differentially expressed miRNAs and studied the mechanisms of adipogenesis and fat deposition between two different pig breeds. In addition, these results also contribute to research relevant to human obesity.

Differential regulation of mRNAs and lncRNAs related to lipid metabolism in Duolang and Small Tail Han sheep

The function of long non-coding RNA (lncRNA) can be achieved through the regulation of target genes, and the deposition of fat is regulated by lncRNA. Fat has an important effect on meat quality. However, there are relatively few studies on lncRNAs in the subcutaneous adipose tissue of Duolang sheep and Small Tail Han sheep. In this study, RNA-Seq technology and bioinformatics methods were used to identify and analyze the lncRNA and mRNA in the subcutaneous adipose tissue of the two breeds of sheep. The results showed that 107 lnRNAs and 1329 mRNAs were differentially expressed. The differentially expressed genes and lncRNA target genes were significantly enriched in the biosynthesis of unsaturated fatty acids signaling pathway, fatty acid metabolism, adipocyte differentiation and other processes related to fat deposition. Among them, LOC105616076, LOC114118103, LOC105607837, LOC101116622, and LOC105603235 target FADS1, SCD, ELOVL6, HSD17B12 and HACD2, respectively. They play a key regulatory role in the biosynthesis of unsaturated fatty acids. This study lays a foundation for the study of the molecular mechanism of lncRNA on fat development, and has reference value for studying the differences in fat deposition between Duolang sheep and Small Tail Han sheep.

Transcriptome Sequencing to Identify Important Genes and lncRNAs Regulating Abdominal Fat Deposition in Ducks

Simple Summary

Abdominal fat deposition affects the quality of duck meat and the feed conversion ratio. Here, we performed transcriptome sequencing of the abdominal adipose tissue of ducks with high and low abdominal fat rate by RNA sequencing, exploring the key regulatory genes and lncRNAs related to abdominal fat deposition. As a result, several candidate genes, lncRNAs, and pathways related to abdominal fat deposition in ducks were retrieved. This study lays the foundations for exploring molecular mechanisms underlying the regulation of abdominal fat deposition in ducks, providing a theoretical reference for breeding high-quality meat-producing ducks.

Abstract

Abdominal fat deposition is an important trait in meat-producing ducks. F2 generations of 304 Cherry Valley and Runzhou Crested White ducks were studied to identify genes and lncRNAs affecting abdominal fat deposition. RNA sequencing was used to study abdominal fat tissue of four ducks each with high or low abdominal fat rates. In all, 336 upregulated and 297 downregulated mRNAs, and 95 upregulated and 119 downregulated lncRNAs were identified. Target gene prediction of differentially expressed lncRNAs identified 602 genes that were further subjected to Gene Ontology and KEGG pathway analysis. The target genes were enriched in pathways associated with fat synthesis and metabolism and participated in biological processes, including Linoleic acid metabolism, lipid storage, and fat cell differentiation, indicating that these lncRNAs play an important role in abdominal fat deposition. PPAPA, FOXO3, FASN, PNPLA2, FKBP5, TCF7L2, BMP2, FGF2, LIFR, ZBTB16, SIRT, GYG2, NCOR1, and NR3C1 were involved in the regulation of abdominal fat deposition. PNPLA2, TCF7L2, FGF2, LIFR, BMP2, FKBP5, GYG2, and ZBTB16 were regulated by the lncRNAs TCONS_00038080, TCONS_0033547, TCONS_00066773, XR_001190174.3, XR_003492471.1, XR_003493494.1, XR_001192142.3, XR_002405656.2, XR_002401822.2, XR_003497063.1, and so on. This study lays foundations for exploring molecular mechanisms underlying the regulation of abdominal fat deposition in ducks and provides a theoretical basis for breeding high-quality meat-producing ducks.

Exploring the effects of different male parent crossings on sheep muscles and related regulatory genes using mRNA-Seq

Objective

With improvements in living standards and increase in population globally, the demand for meat products has been increasing; improved meat production from livestock could effectively meet this demand. In this study, we examined the differences in the muscle traits of different male crossbred sheep and attempted to identify key genes that regulate these traits.

Methods

Dubo sheep × Small-tailed Han sheep (DP×STH) and Suffolk × Small-tailed Han sheep (SFK×STH) were selected to determine meat quality and production performance by Masson staining. Transcriptome sequencing and bioinformatic analysis were performed to identify differentially expressed genes (DEGs) related to meat quality. The presence of DEGs was confirmed by real-time polymerase chain reaction (qPCR).

Results

The production performance of SFK×STH sheep was better than that of DP×STH sheep, but the meat quality of DP×STH sheep was better than that of SFK×STH sheep. The muscle fiber diameter of DP×STH sheep was smaller than that of SFK×STH sheep. Twenty-two DEGs were identified. Among them, four Gene Ontology terms were related to muscle traits, and three DEGs were related to muscle or muscle fibers. There were no significant differences in the number of single nucleotide mutations and mutation sites in the different male parent cross combinations.

Conclusion

This study provides genetic resources for future sheep muscle development and cross-breeding research.

Interference With ACSL1 Gene in Bovine Adipocytes: Transcriptome Profiling of mRNA and lncRNA Related to Unsaturated Fatty Acid Synthesis

The enzyme long-chain acyl-CoA synthetase 1 (ACSL1) is essential for lipid metabolism. The ACSL1 gene controls unsaturated fatty acid (UFA) synthesis as well as the formation of lipid droplets in bovine adipocytes. Here, we used RNA-Seq to determine lncRNA and mRNA that regulate UFA synthesis in bovine adipocytes using RNA interference and non-interference with ACSL1. The corresponding target genes of differentially expressed (DE) lncRNAs and the DE mRNAs were found to be enriched in lipid and FA metabolism-related pathways, according to GO and KEGG analyses. The differentially expressed lncRNA- differentially expressed mRNA (DEL-DEM) interaction network indicated that some DELs, such as TCONS_00069661, TCONS_00040771, TCONS_ 00035606, TCONS_00048301, TCONS_001309018, and TCONS_00122946, were critical for UFA synthesis. These findings assist our understanding of the regulation of UFA synthesis by lncRNAs and mRNAs in bovine adipocytes.

Differential Expression and Functional Analysis of CircRNA in the Ovaries of Low and High Fecundity Hanper Sheep

Simple Summary

Litter size is an important trait affecting reproductive capacity and breeding economics in meat sheep. Consequently, revealing its molecular mechanism helps us understand multiple lambs from the genetic perspective. In this study, we provide a genome-wide expression profile of circular RNAs (circRNAs) expression in Hanper sheep, which is a new breed of meat sheep raised by cross and self-group breeding for 15 years. In this study, ovarian circular RNAs and miRNAs associated with high and low fertility Hanper sheep are identified during the follicular and luteal phases of the estrous cycle, and their potential biological functions are predicted through Gene Ontology (GO), KEGG, GSEA, STEM, WGCNA analysis.

Abstract

Litter size is a considerable quality that determines the production efficiency of mutton sheep. Therefore, revealing the molecular regulation of high and low fertility may aid the breeding process to develop new varieties of mutton sheep. CircRNAs are the important factors regulating follicular development, but their mechanism role in the regulation of litter size in Hanper sheep is not clear. In the present study, ovarian tissues from the follicular (F) or luteal phase (L) of Hanper sheep that were either consecutive monotocous (M) or polytocous were collected. Then, we performed transcriptome sequencing to screen for differentially expressed circRNAs (DE-circRNAs) and elucidate their function. In total, 4256 circRNA derived from 2184 host genes were identified in which 183 (146 were upregulated, while 37 were downregulated) were differentially expressed in monotocous sheep in the follicular phase versus polytocous sheep in the follicular phase (MF vs. PF). Moreover, 34 circRNAs (14 were upregulated, while 20 were downregulated) were differentially expressed in monotocous sheep in the luteal phase versus polytocous sheep in the luteal sheep (ML vs. PL). This was achieved through DE-circRNAs function enrichment annotation analysis by GESA, GO, and KEGG, which function through the EGF-EGFR-RAS-JNK, TGF-β and thyroid hormone signaling pathway to affect the litter size of Hanper sheep in MF vs. PF and ML vs. PL. STEM results showed that MAPK signaling pathways play a key role in MF vs. PF and ML vs. PL. Through WGCNA analysis, AKT3 was a core gene in MF vs. PF and ML vs. PL. Moreover, competitive endogenous RNA (ceRNA) network analysis revealed the target binding sites for miRNA such as oar-miR-27a, oar-miR-16b, oar-miR-200a/b/c, oar-miR-181a, oar-miR-10a/b, and oar-miR-432 in the identified DE-cirRNAs.

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.

Differential expression and prediction of function of lncRNAs in the ovaries of low and high fecundity Hanper sheep

Litter size is an important trait that determines the production efficiency of sheep bred for meat. Its detailed investigation can reveal the molecular mechanisms that control the fecundity of sheep and possibly accelerate the breeding process of new varieties of sheep that have high prolificacy. Long non-coding RNAs (lncRNAs) have proven to be an important factor in the regulation of follicular development. However, the mechanisms by which lncRNAs regulate litter size in sheep remain unclear. In the present study, ovarian tissues from the follicular (F) or luteal phase (L) of Hanper sheep that were either monotocous (M) or polytocous (P; FM, FP, LM and LP groups) were collected and sequenced to identify differentially expressed lncRNAs and predict their function. The results indicate that the number of up- and down-regulated lncRNAs in the follicular phase (FM vs. FP) was 95 and 111 and 109 and 49, respectively, in the luteal phase (LM vs. LP). The functional enrichment of the different lncRNAs coexpressed with mRNA was analysed. The results demonstrated that the KISS1-GnRH-LH/FSH-E2 and EGF-EGFR-RAS-PI3K signalling pathways promoted the initiation of the primordial period, follicular development and ovulation in the follicular phase (FM vs. FP). During the luteal phase (LM vs. LP), the production and development of the corpus luteum in ewes was influenced by the KITLG-KIT/FGF-FGFR/HGF-MET-RAS-ERK signalling pathway. STEM clustering functional enrichment analysis of the differentially expressed lncRNAs indicated that profile11 was principally enriched in the Cytokine-Jak-STAT, PDGF-PDGFR-PI3K and KITLG-KIT-RAS-ERK signalling pathways. By analysis of the differential expression of the lncRNAs and their expression in each group, lncRNAs Xist (loc101112291) and Gtl2 (loc101123329) were found to be highly expressed, suggesting that regulation of follicular development was mediated through methylation processes.

Genome-wide association study of four composite reproductive traits in Iranian fat-tailed sheep

Composite reproductive traits are a combination of growth and reproductive traits. They have the advantage of being better attuned to the market drivers since producers are paid on a per kilogram basis and not on a per head basis. In this study, 124 Lori---Bakhtiari ewes were genotyped using the medium-density Illumina Ovine SNP50 array. A genome-wide association study was performed on estimated breeding values of four composite reproductive traits and genetic parameters were also estimated. The traits were litter mean weight at birth, litter mean weight at weaning, total litter weight at birth and total litter weight at weaning. Several suggestive and associated single nucleotide polymorphisms (SNPs) were identified. Neighbouring the top SNPs there were five genes, inhibin β E subunit (INHBE), inhibin β C subunit (INHBC), testis expressed 12 (TEX12), β-carotene oxygenase 2 (BCO2) and WD Repeat Domain 70 (WDR70) identified as possible candidate genes for composite reproductive traits of the Lori-Bakhtiari sheep. These genes are in pathways known to be relevant to fertility and growth characteristics. The results provide new information for the functional annotation of genes associated with fertility traits and add new evidence towards a consensus of quantitative trait loci associated with reproductive traits in sheep.

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.

Long non-coding RNA366.2 controls endometrial epithelial cell proliferation and migration by upregulating WNT6 as a ceRNA of miR-1576 in sheep uterus

Long non-coding RNAs (lncRNAs) play an important regulatory role in mammalian fecundity. Currently, most studies are primarily concentrated on ovarian lncRNAs, ignoring the influence of uterine lncRNAs on the fecundity of female sheep. In this study, we found a higher density of uterine glands and endometrial microvessel density (MVD) in high prolificacy group of Hu sheep compared to low prolificacy groups (p < 0.05) as well as an increased level of serum placental growth factor (PLGF). Hundreds of differentially expressed (DE) lncRNAs were identified in Hu sheep with different fecundity by RNA sequencing (RNA-seq), and their targets were enriched in some signaling pathways involved in endometrial functions, such as the estrogen signaling pathway, nuclear factor kappa B (NF-κB) signaling pathway, oxytocin signaling pathway, and Wnt signaling pathway. Furthermore, the underlying mechanisms of competitive endogenous RNA (ceRNA) of lncRNA366.2-miR-1576- WNT6 were determined by bioinformatics analysis. Functionally, our results indicated that lncRNA366.2 promoted endometrial epithelial cell (EEC) proliferation, migration, and growth factor expression by sponging miR-1576 to upregulate WNT6 expression and activate the Wnt/β-catenin pathway. Taken together, our research indicated the regulatory mechanism of the lncRNA366.2-miR-1576-WNT6 in EEC proliferation and migration. Furthermore, this study provides a new theoretical reference for the identification of candidate genes related to fecundity.

Eukaryotic expression, Co-IP and MS identify BMPR-1B protein-protein interaction network

BACKGROUND

BMPR-1B is part of the transforming growth factor β super family and plays a pivotal role in ewe litter size. Functional loss of exon-8 mutations in the BMPR-1B gene (namely the FecB gene) can increase both the ewe ovulation rate and litter size.

RESULTS

This study constructed a eukaryotic expression system, prepared a monoclonal antibody, and characterized BMPR-1B/FecB protein-protein interactions (PPIs). Using Co-immunoprecipitation coupled to mass spectrometry (Co-IP/MS), 23 proteins were identified that specifically interact with FecB in ovary extracts of ewes. Bioinformatics analysis of selected PPIs demonstrated that FecB associated with several other BMPs, primarily via signal transduction in the ovary. FecB and its associated interaction proteins enriched the reproduction process via BMP2 and BMP4 pathways. Signal transduction was identified via Smads proteins and TGF-beta signaling pathway by analyzing the biological processes and pathways. Moreover, other target proteins (GDF5, GDF9, RhoD, and HSP 10) that interact with FecB and that are related to ovulation and litter size in ewes were identified.

CONCLUSIONS

In summary, this research identified a novel pathway and insight to explore the PPi network of BMPR-1B.

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.

Comparative Transcriptome Analysis Identifying the Different Molecular Genetic Markers Related to Production Performance and Meat Quality in Longissimus Dorsi Tissues of MG × STH and STH Sheep

Crossbred sheep have many prominent traits, such as excellent production performance and high-quality meat, when compared to local sheep breeds. However, the genetic molecular markers related to these characteristics remain unclear. The crossbred MG × STH (small-tailed Han sheep (STH) × Mongolian sheep (MG)) breed and the STH breed were selected to measure production performance and meat quality. We used 14 indexes of production performance and meat quality, which in the MG × STH population showed significant differences compared to the STH breed. Subsequently, the longissimusdorsi from the two sheep were subjected to comparative transcriptomic analyses to identify differentially expressed genes (DEGs) related to production performance and meat quality. A total of 874 DEGs were identified between the two sheep groups. A total of 110 unique DEGs related to sheep production performance and meat quality were selected as the candidate DEGs. We found 6 production-performance-related and 30 meat-quality-related DEGs through a correlation analysis, including SPARC, ACVRL1, FNDC5 and FREM1. The expression levels of 11 DEGs were validated by real-time PCR, and the results were in accordance with the results of the comparative transcriptomic and correlation analyses. These results will assist in understanding sheep heterosis and molecular marker-assisted selection.

Transcriptome profile of skeletal muscle at different developmental stages in Large White and Mashen pigs

From the perspectives of promoting individual growth and development, increasing pork yield, and improving feed utilization, it is desirable to screen candidate genes underlying pig muscle growth and regulation. In this study, we investigated transcriptome differences at 1, 90, and 180 d of age in Large White and Mashen pigs, characterized differentially expressed genes (DEGs), and screened candidate genes affecting skeletal muscle growth and development. RNA-seq was applied to analyze the transcriptome of the longissimus dorsi (LD) in the two breeds. In LD samples from the two breeds at three growth stages, 7215, 6332, 237, 3935, 3404, and 846 DEGs were obtained for L01 vs. L90, L01 vs. L180, L90 vs. L180, MS01 vs. MS90, MS01 vs. MS180, and MS90 vs. MS180, respectively. Significant tendencies in DEG expression could be grouped into eight profiles. Based on the functional analysis of DEGs, 16 candidate genes related to skeletal muscle growth and development were identified, including PCK2, GNAS, ADCY2, PRKAB1, PRKAB2, PRKAG1, PRKAG2, PHKA1, PHKA2, PHKG1, PHKG2, ITPR3, IGF1R, FGFR4, FGF1, and FGF18. The results of this study thus provide a theoretical basis for the mechanisms and candidate genes underlying skeletal muscle development in pigs.

Correlation analysis of candidate gene SNP for high-yield in Liaoning cashmere goats with litter size and cashmere performance

This study was designed to identify the relationship of four genes (GDF9, BMPR-IB, FecB and ESR) polymorphisms in the 3'UTR region with litter size and cashmere performance of Liaoning cashmere goats (LCG, n = 1140). The ESR C463T and T575G loci of LCG were genotyped. The results of correlation analysis showed that five effective single nucleotide polymorphisms (SNPs) loci (C47T, C94T, C299T, C463T and T575G) were found in the four genes. The lambing number of CC and CT genotypic individuals at FecB C94T locus was significantly higher than that of TT genotypic individuals (45.7 and 46.8%, respectively); the lambing number of CC genotypic individuals at ESR C463T locus was significantly higher than that of CT, TT genotypic individuals (9 and 15%, respectively); There was a positive correlation between CC genotype at C463T locus and cashmere fineness. In this study, the relationship between FecB C94T and ESR C463T loci C alleles and lambing number in LCG was preliminarily revealed. These results further confirmed that FecB and ESR genes may be significantly correlated with high fecundity of LCG.

Identification and characterization of mRNAs and lncRNAs in the uterus of polytocous and monotocous Small Tail Han sheep (Ovis aries)

Background

Long non-coding RNAs (lncRNAs) regulate endometrial secretion and uterine volume. However, there is little research on the role of lncRNAs in the uterus of Small Tail Han sheep (FecB++). Herein, RNA-seq was used to comparatively analyze gene expression profiles of uterine tissue between polytocous and monotocous sheep (FecB++) in follicular and luteal phases.

Methods

To identify lncRNA and mRNA expressed in the uterus, the expression of lncRNA and mRNA in the uterus of Small Tail Han sheep (FecB++) from the polytocous group (n = 6) and the monotocous group (n = 6) using RNA-sequencing and real-time polymerase chain reaction (RT-PCR). Identification of differentially expressed lncRNAs and mRNAs were performed between the two groups and two phases . Gene ontology (GO) and pathway enrichment analyses were performed to analyze the biological functions and pathways for the differentially expressed mRNAs. LncRNA-mRNA co-expression network was constructed to further analyses the function of related genes.

Results

In the follicular phase, 473 lncRNAs and 166 mRNAs were differentially expressed in polytocous and monotocous sheep; in the luteal phase, 967 lncRNAs and 505 mRNAs were differentially expressed in polytocous and monotocous sheep. GO and KEGG enrichment analysis showed that the differentially expressed lncRNAs and their target genes are mainly involved in ovarian steroidogenesis, retinol metabolism, the oxytocin signaling pathway, steroid hormone biosynthesis, and the Foxo signaling pathway. Key lncRNAs may regulate reproduction by regulating genes involved in these signaling pathways and biological processes. Specifically, UGT1A1, LHB, TGFB1, TAB1, and RHOA, which are targeted by MSTRG.134747, MSTRG.82376, MSTRG.134749, MSTRG.134751, and MSTRG.134746, may play key regulatory roles. These results offer insight into molecular mechanisms underlying sheep prolificacy.

Single nucleotide polymorphisms in BMP2 and BMP7 and the association with litter size in Small Tail Han sheep

Although it has been investigated for many years, the physiological processes regulating prolificacy in sheep remains unclear because of regulation by many genes. To better understand the effects of three single nucleotide polymorphisms (SNPs) comprising g.48462350C>T in BMP2, g.58171856C>G and g.58171886A>C in BMP7, a population genetic analysis was conducted using data obtained from genotyping in 768 sheep from six breeds (three polytocous and three monotocous). The results indicate that all the sheep breeds were considered to conform to the Hardy-Weinberg equilibrium (P > 0.05). The associations of these three SNPs with litter size in 384 Small Tail Han sheep were analyzed, therefore, and found to be correlated with fecundity as assessed by mean litter size (P < 0.05). Bioinformatic analysis indicated there was a transmembrane domain change that occurred after a mutation in BMP2 at g.48462350C>T, and changes involving transcription factors such as USF1, USF2 and INMS1 in the BMP7 promoter region might be involved in greater sheep prolificacy.

Integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old Small Tail Han Sheep ovaries reveals that oar-miR-432 downregulates RPS6KA1 expression

Small Tail Han Sheep are an excellent local sheep breed in China, and their outstanding reproductive performance is one of their very important biological characteristics. Clarifying the ovary development process of these ewes should provide a theoretical basis for improving their reproductive efficiency. In this study, we identified the differentially expressed (DE) microRNAs (miRNAs) in 2-, 6-, and 12-month-old small-tail Han sheep ovaries by constructing and analyzing the miRNA expression profiles. These findings clarify the molecular mechanisms regulating the excellent reproductive performance of small-tail Han ewes. We used RNA-Seq technology and bioinformatic to analyze these profiles. Eleven, 13, and 19 DE miRNAs were identified in the 2- vs 6-, 6- vs 12-, and 2- vs 12-month-old ovaries, respectively. In total, 54, 37, and 198 predicted target genes of these DE miRNAs were identified in these three groups, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that in the 2- vs 6-month-old ovaries, the target genes of DE known sheep miRNAs were involved in 102 GO terms and seven signaling pathways; in the 6- vs 12-month-old ovaries, the target genes of DE known sheep miRNAs were involved in 52 GO terms and three signaling pathways; and in the 2- vs 12-month-old ovaries, the target genes of DE known sheep miRNAs were involved in 88 GO terms and six signaling pathways. Three miRNA-target regulatory networks were constructed based on these DE miRNA-target interactions. Nine miRNAs were selected to confirm to the accuracy of the miRNA sequencing data with qRT-PCR. The site at which oar-miR-432 binds RPS6KA1 was determined with a dual-luciferase system. This is the first integrated analysis the expression profiles of miRNAs and their targets during ovarian development in small-tail Han sheep. These data clarify the molecular regulatory mechanisms underlying sheep ovarian development and identify biomarkers that influence the reproductive performance of small-tail Han ewes.

Pituitary Transcriptomic Study Reveals the Differential Regulation of lncRNAs and mRNAs Related to Prolificacy in Different FecB Genotyping Sheep

Long non-coding RNA (LncRNA) have been identified as important regulators in the hypothalamic-pituitary-ovarian axis associated with sheep prolificacy. However, their expression pattern and potential roles in the pituitary are yet unclear. To explore the potential mRNAs and lncRNAs that regulate the expression of the genes involved in sheep prolificacy, we used stranded specific RNA-seq to profile the pituitary transcriptome (lncRNA and mRNA) in high prolificacy (genotype FecB BB, litter size = 3; H) and low prolificacy sheep (genotype FecB B+; litter size = 1; L). Our results showed that 57 differentially expressed (DE) lncRNAs and 298 DE mRNAs were found in the pituitary between the two groups. The qRT-PCR results correlated well with the RNA-seq results. Moreover, functional annotation analysis showed that the target genes of the DE lncRNAs were significantly enriched in pituitary function, hormone-related pathways as well as response to stimulus and some other terms related to reproduction. Furthermore, a co-expression network of lncRNAs and target genes was constructed and reproduction related genes such as SMAD2, NMB and EFNB3 were included. Lastly, the interaction of candidate lncRNA MSTRG.259847.2 and its target gene SMAD2 were validated in vitro of sheep pituitary cells. These differential mRNA and lncRNA expression profiles provide a valuable resource for understanding the molecular mechanisms underlying Hu sheep prolificacy.

Integrated Analysis of miRNA and mRNA Expression Profiles Reveals Functional miRNA-Targets in Development Testes of Small Tail Han Sheep

Small Tail Han Sheep is a highly valued local breed in China because of their precocity, perennial estrus, and high fecundity. The average annual lambing rate of ewes is as high as 180–270%, the semen of ram has characteristics of high yield, high density, and good motility. To reveal the key miRNAs and miRNA-targets underlying testis development and spermatogenesis in male Small Tail Han Sheep, integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old testes was performed by RNA-seq technology and bioinformatics methods. The results showed that total of 153 known sheep miRNAs and 2712 novel miRNAs were obtained in 2-,6 - and 12-month-old Small Tail Han Sheep testes; 5, 1, and 4 differentially expressed (DE) known sheep miRNAs, and 132, 105, and 24 DE novel miRNAs were identified in 2- vs. 6-, 6- vs. 12-, and 2- vs. 12-month-old testes, respectively. We combined miRNA results of this study and the mRNA results obtained in our previous study to predict the target mRNAs of DE known sheep miRNAs; 131, 10, and 15 target mRNAs of DE known sheep miRNAs and 76, 1, and 11 DE miRNA–targets were identified in the three groups, respectively. GO and KEGG analyses showed that: in 2- vs. 6-month-olds, the target genes of DE known sheep miRNAs were involved in 100 biological processes and 11 signaling pathways; in 6- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 4 biological processes; and in 2- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 17 biological processes and 4 signaling pathways. Three miR–target regulatory networks were constructed based on these DE miRNA–targets. The key miRNA-Targets involved in testis development and spermatogenesis were screened. 6 known sheep miRNAs and 6 novel miRNAs were selected to validate the accuracy of miRNA sequencing data by qRT-PCR. The binding sites of oar-miR-379-5p with WNT8A was validated by a dual luciferase reporter gene detection system.

Differential expression of mRNA-miRNAs related to intramuscular fat content in the longissimus dorsi in Xinjiang brown cattle

In this study, we examined the role of mRNAs and miRNAs in variations in intramuscular fat content in the longissimus dorsi muscle in Xinjiang brown cattle. Two groups of Xinjiang brown cattle with extremely different intramuscular fat content in the longissimus dorsi were selected for combined of miRNA and mRNA analysis using an RNA-Seq. In total, 296 mRNAs and 362 miRNAs were significantly differentially expressed, including 155 newly predicted miRNAs, 275 significantly upregulated genes, 252 significantly upregulated miRNAs, 21 significantly downregulated genes and 110 significantly downregulated miRNAs. The combined miRNA and mRNA analysis identified 96 differentially expressed miRNAs and 27 differentially expressed mRNAs. In all, 47 upregulated miRNAs had a regulatory effect on 14 differentially downregulated target genes, and 49 downregulated miRNAs had a regulatory effect on 13 upregulated target genes. To verify the sequencing results, 10 differentially expressed genes (DEGs) and 10 differentially expressed miRNAs were selected for qRT-PCR. The qRT-PCR results confirmed the sequencing results. The results of this study shed light on the molecular regulation of bovine adipose tissue, which might help with the development of new strategies for improving meat quality and animal productivity in beef cattle to provide healthier meat products for consumers.

Genome -wide transcriptome profiling in ovaries of small-tail Han sheep during the follicular and luteal phases of the oestrous cycle

Ovarian function, the control of which is predominantly integrated by the hypothalamic-pituitary-ovarian axis, is pivotal for maintaining reproductive efficiency in sheep. To understand the regulatory mechanism of the follicular-luteal phase transition in prolific sheep, the genome-wide expression patterns of microRNAs and genes in ovarian tissue of Small Tail Han sheep were examined during the follicular and luteal phases of the oestrous cycle. Differentially expressed genes (n = 450) were identified and a total of 139 known and 72 novel miRNAs were identified in the two libraries. It is suggested that differential abundance of miR-200a, 200b and 200c may have an important role in the follicular-luteal transition. A miRNA-regulated gene expression network was created for exploring the regulation of the follicular-luteal transition and quantitative real-time PCR verified the reliability of the RNA-seq data.

Comparative mRNA and miRNA expression in European mouflon (Ovis musimon) and sheep (Ovis aries) provides novel insights into the genetic mechanisms for female reproductive success

Prolific breeds of domestic sheep (Ovis aries) are important genetic resources due to their reproductive performance, which is characterized by multiple lambs per birth and out-of-season breeding. However, the lack of a comprehensive understanding of the genetic mechanisms underlying the important reproductive traits, particularly from the evolutionary genomics perspective, has impeded the efficient advancement of sheep breeding. Here, for the first time, by performing RNA-sequencing we built a de novo transcriptome assembly of ovarian and endometrial tissues in European mouflon (Ovis musimon) and performed an mRNA–miRNA integrated expression profiling analysis of the wild species and a highly prolific domestic sheep breed, the Finnsheep. We identified several novel genes with differentially expressed mRNAs (e.g., EREG, INHBA, SPP1, AMH, TDRD5, and ZP2) between the wild and domestic sheep, which are functionally involved in oocyte and follicle development and fertilization, and are significantly (adjusted P-value < 0.05) enriched in the Gene Ontology (GO) terms of various reproductive process, including the regulation of fertilization, oogenesis, ovarian follicle development, and sperm–egg recognition. Additionally, we characterized 58 differentially expressed miRNAs and 210 associated target genes that are essential for the regulation of female reproduction cycles through specific regulatory networks [e.g., (miR-136, miR-374a, miR-9-5p)-(EREG, INHBA)]. Furthermore, our integrated mRNA and miRNA expression profiling analysis elucidated novel direct and indirect miRNA/mRNA causal regulatory relationships related to the reproductive traits of the Ovis species. This study provides in-depth insights into the genomic evolution underlying the reproductive traits of the Ovis species and valuable resources for ovine genomics.

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.

Metabolic Effects of FecB Gene on Follicular Fluid and Ovarian Vein Serum in Sheep (Ovis aries)

The FecB gene has been discovered as an important gene in sheep for its high relationship with the ovulation rate, but its regulatory mechanism remains unknown. In the present study, liquid chromatography-mass spectrometry (LC-MS) and gas chromatography-mass spectrometry (GC-MS) techniques were adopted to detect the metabolic effects of FecB gene in follicular fluid (FF) and ovarian vein serum (OVS) in Small Tail Han (STH) sheep. ANOVA and random forest statistical methods were employed for the identification of important metabolic pathways and biomarkers. Changes in amino acid metabolism, redox environment, and energy metabolism were observed in FF from the three FecB genotype STH ewes. Principal component analysis (PCA) and hierarchical clustering analysis (HCA) showed that metabolic effects of FecB gene are more pronounced in FF than in OVS. Therefore, the difference of the metabolic profile in FF is also affected by the FecB genotypes. In Spearman correlation analysis, key metabolites (e.g., glucose 6-phosphate, glucose 1-phosphate, aspartate, asparagine, glutathione oxidized (GSSG), cysteine-glutathione disulfide, γ-glutamylglutamine, and 2-hydrosybutyrate) in ovine FF samples showed a significant correlation with the ovulation rate. Our findings will help to explain the metabolic mechanism of high prolificacy ewes and benefit fertility identification.

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.

RNA-Seq analysis on chicken taste sensory organs: An ideal system to study organogenesis

RNA-Seq is a powerful tool in transcriptomic profiling of cells and tissues. We recently identified many more taste buds than previously appreciated in chickens using molecular markers to stain oral epithelial sheets of the palate, base of oral cavity, and posterior tongue. In this study, RNA-Seq was performed to understand the transcriptomic architecture of chicken gustatory tissues. Interestingly, taste sensation related genes and many more differentially expressed genes (DEGs) were found between the epithelium and mesenchyme in the base of oral cavity as compared to the palate and posterior tongue. Further RNA-Seq using specifically defined tissues of the base of oral cavity demonstrated that DEGs between gustatory (GE) and non-gustatory epithelium (NGE), and between GE and the underlying mesenchyme (GM) were enriched in multiple GO terms and KEGG pathways, including many biological processes. Well-known genes for taste sensation were highly expressed in the GE. Moreover, genes of signaling components important in organogenesis (Wnt, TGFβ/ BMP, FGF, Notch, SHH, Erbb) were differentially expressed between GE and GM. Combined with other features of chicken taste buds, e.g., uniquely patterned array and short turnover cycle, our data suggest that chicken gustatory tissue provides an ideal system for multidisciplinary studies, including organogenesis and regenerative medicine.

Cutaneous transcriptome analysis in NIH hairless mice

Mice with spontaneous coat mutations are ideal animal models for studying skin development and tumorigenesis. In this study, skin hair growth cycle abnormalities were examined in NIH hairless mice 42 days after birth (P42) by using hematoxylin-eosin (H&E) staining. To examine the gene expression patterns in the skin of mutant mice, the dorsal skin of P42 female NIH mice and NIH hairless mice was sequenced by RNA-Seq, and 5,068 differentially expressed genes (DEGs) were identified (false discovery rate [FDR] ≥ 2, P < 0.05). A pathway analysis showed that basal cell carcinoma, the cell cycle and the Hippo, Hedgehog and Wnt signaling pathways were up-regulated in NIH hairless mice. Previous studies have shown that these pathways are closely associated with cell proliferation, cell cycle, organ size and cancer development. In contrast, signal transduction, bacterial and parasitic infection, and receptor-mediated pathways, including calcium signaling, were down-regulated in NIH hairless mice. A gene interaction network analysis was performed to identify genes related to hair follicle development. To verify the reliability of the RNA-Seq results, we used q-PCR to analyze 12 key genes identified from the gene interaction network analysis, including eight down-regulated and four up-regulated genes, and the results confirmed the reliability of the RNA-Seq results. Finally, we constructed the differential gene expression profiles of mutant mice by RNA-Seq. NIH hairless mice exhibited abnormalities in hair development and immune-related pathways. Pik3r1 and Pik3r3 were identified as key genes, laying the foundation for additional in-depth studies of hairless mice.

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.

A novel 29 bp insertion/deletion (indel) variant of the <i>LHX3</i> gene and its influence on growth traits in four sheep breeds of various fecundity

Abstract. Belonging to the same LIM homeobox (LHX) family, LHX3 and LHX4 are key transcription factors in animal growth and reproduction. Insertion/deletion (indel) is a relatively simple and effective DNA marker. Therefore, four sheep breeds of various fecundity were used to explore the novel indel variants within the sheep LHX3 and LHX4 gene, as well as to evaluate their effects on growth traits. Herein, only one novel 29 bp indel (NC_019460.2:g.3107494-3107522delGGCCTGGACTGTGATGGGCACCCTCCGGG) within the sheep LHX3 gene was found, and three genotypes were detected. Interestingly, the increasing trends of II (insertion/insertion) genotype frequency and I allelic frequency were the same as the growth of the fertility character. Genotypic frequency and allelic frequency distributions were significantly different between the high-fecundity breeds (HS, STHS and LFTS) and low-fecundity breed (TS) based on a χ2 test (P < 0.05). Association analyses showed that body length was significantly different in female TS and STHS and that chest width was significantly different for the female TS and male STHS (P < 0.05). These findings suggested that the 29 bp indel could extend the spectrum of genetic variations of the LHX3 gene in sheep and provide a valuable theoretical basis for the marker-assisted selection (MAS) in sheep breeding and genetics.

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.

Expression Analysis of Genes Associated with Prolificacy in FecB Carrier and Noncarrier Indian Sheep

The effect of FecB mutation on the gene expression in FecB carrier and noncarrier estrous synchronized ewes, has been analyzed. For this study the whole ovarian tissues and Graafian follicles were collected from estrus synchronized FecB carrier Garole, and non-carrier Deccani Indian sheep, showing remarkable differences in the numbers of preovulatory follicles among two groups. Eleven potential candidate genes (BMP15, GDF9, BMP4, BMP7, BMPR1B, BMPR1A, SMAD9, LHCGR, FSHR, IGF1R, and STAT5) were selected for their expression analysis by SybrGreen based real-time PCR, across ovaries and Graafian follicles of different fecundity groups, for having better insights into the effect of FecB genotypes on follicular development. Variable expression was observed for almost all the genes included in the present study among high and low fecundity groups that was most significant for the BMP7, BMP4, LHCGR, and FSHR transcripts in the ovarian follicles of high and low fecundity ewes, indicating their importance in governing the fecundity in FecB carrier, Indian Garole sheep. BMP4 expression among the genes studied was significantly higher in FecB carrier Garole sheep. This study confirms the changes in mRNA expression of the genes implicated in follicular development in FecB carrier and noncarrier Indian sheep breeds.

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.

Identification and Classification of New Transcripts in Dorper and Small-Tailed Han Sheep Skeletal Muscle Transcriptomes

High-throughput mRNA sequencing enables the discovery of new transcripts and additional parts of incompletely annotated transcripts. Compared with the human and cow genomes, the reference annotation level of the sheep genome is still low. An investigation of new transcripts in sheep skeletal muscle will improve our understanding of muscle development. Therefore, applying high-throughput sequencing, two cDNA libraries from the biceps brachii of small-tailed Han sheep and Dorper sheep were constructed, and whole-transcriptome analysis was performed to determine the unknown transcript catalogue of this tissue. In this study, 40,129 transcripts were finally mapped to the sheep genome. Among them, 3,467 transcripts were determined to be unannotated in the current reference sheep genome and were defined as new transcripts. Based on protein-coding capacity prediction and comparative analysis of sequence similarity, 246 transcripts were classified as portions of unannotated genes or incompletely annotated genes. Another 1,520 transcripts were predicted with high confidence to be long non-coding RNAs. Our analysis also revealed 334 new transcripts that displayed specific expression in ruminants and uncovered a number of new transcripts without intergenus homology but with specific expression in sheep skeletal muscle. The results confirmed a complex transcript pattern of coding and non-coding RNA in sheep skeletal muscle. This study provided important information concerning the sheep genome and transcriptome annotation, which could provide a basis for further study.

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.