Identification and characterization of miRNAs in the ovaries of a highly prolific sheep breed

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

miR-1306 induces cell apoptosis by targeting BMPR1B gene in the ovine granulosa cells

Bone morphogenetic protein receptor type-1B (BMPR1B) is one of the major gene for sheep prolificacy. However, few studies investigated its regulatory region. In this study, we reported that miR-1306 is a direct inhibitor of BMPR1B gene in the ovine granulosa cells (ovine GCs). We detected a miRNA response element of miR-1306 in the 3’ untranslated region of the ovine BMPR1B gene. Luciferase assay showed that the ovine BMPR1B gene is a direct target of miR-1306. qPCR and western blotting revealed that miR-1306 reduces the expression of BMPR1B mRNA and protein in the ovine granulosa cells. Furthermore, miR-1306 promoted cell apoptosis by suppressing BMPR1B expression in the ovine granulosa cells. Overall, our results suggest that miR-1306 is an epigenetic regulator of BMPR1B, and may serve as a potential target to improve the fecundity of sheep.

Meta-analysis of genome-wide association studies for litter size in sheep

Litter size and ovulation rate are important reproduction traits in sheep and have important impacts on the profitability of farm animals. To investigate the genetic architecture of litter size, we report the first meta-analysis of genome-wide association studies (GWAS) using 522 ewes and 564,377 SNPs from six sheep breeds. We identified 29 significant associations for litter size which 27 of which have not been reported in individual GWAS for each population. However, we could confirm the role of BMPR1B in prolificacy. Our gene set analysis discovered biological pathways related to cell signaling, communication, and adhesion. Functional clustering and enrichment using protein databases identified epidermal growth factor-like domain affecting litter size. Through analyzing protein-protein interaction data, we could identify hub genes like CASK, PLCB4, RPTOR, GRIA2, and PLCB1 that were enriched in most of the significant pathways. These genes have a role in cell proliferation, cell adhesion, cell growth and survival, and autophagy. Notably, identified SNPs were scattered on several different chromosomes implying different genetic mechanisms underlying variation of prolificacy in each breed. Given the different layers that make up the follicles and the need for communication and transfer of hormones and nutrients through these layers to the oocyte, the significance of pathways related to cell signaling and communication seems logical. Our results provide genetic insights into the litter size variation in different sheep breeds.

MicroRNAs: Potential biomarkers for reproduction, diagnosis, prognosis, and therapeutic in domestic animals

MicroRNA (miRNAs) are small non-coding RNA molecules involved in a wide range of biological processes through the post-transcriptional regulation of gene expression. Most studies evaluated microRNA expression in human, and despite fewer studies in veterinary medicine, this topic is one of the most exciting areas of modern veterinary medicine. miRNAs showed to be part of the pathogenesis of diseases and reproduction physiology in animals, making them biomarkers candidates. This review provides an overview of the current knowledge regarding miRNAs' role in reproduction and animal diseases, diagnostic and therapy.

miRNA Regulatory Functions in Farm Animal Diseases, and Biomarker Potentials for Effective Therapies

MicroRNAs (miRNAs) are small endogenous RNAs that regulate gene expression post-transcriptionally by targeting either the 3′ untranslated or coding regions of genes. They have been reported to play key roles in a wide range of biological processes. The recent remarkable developments of transcriptomics technologies, especially next-generation sequencing technologies and advanced bioinformatics tools, allow more in-depth exploration of messenger RNAs (mRNAs) and non-coding RNAs (ncRNAs), including miRNAs. These technologies have offered great opportunities for a deeper exploration of miRNA involvement in farm animal diseases, as well as livestock productivity and welfare. In this review, we provide an overview of the current knowledge of miRNA roles in major farm animal diseases with a particular focus on diseases of economic importance. In addition, we discuss the steps and future perspectives of using miRNAs as biomarkers and molecular therapy for livestock disease management as well as the challenges and opportunities for understanding the regulatory mechanisms of miRNAs related to disease pathogenesis.

Identification and characterization of miRNAs during early pregnancy in domestic sheep

MicroRNA resources in sheep are limited compared with those in other domesticated mammalian species. By sequencing small RNAs of sheep corpus luteum and endometrium, we have generated the largest amount of miRNA‐seq data and compiled the most comprehensive list thus far of miRNAs (n = 599) in sheep. Additionally, we observed a highly conserved maternally imprinted cluster of miRNAs on chromosome 18 homologous to that found on chromosome 14 in human and several other eutherian mammals.

Single nucleotide polymorphism rs737028527 (G>A) affect miR-1b-3p biogenesis and effects on chicken egg-laying traits

The abundance of miR-1b-3p in the chicken ovary is greater after sexual maturation. In the present study there was assessment of whether a single nucleotine polymorphism (SNP) led to an alteration in expression of reproductive traits. The miR-1b-3p abundance was greatest in ovarian follicles The SNP site of rs737028527 (G > A), located in the 734 bp upstream region of pre-miR-1b-3p, was identified in three different chicken breeds. Results from an association analysis of chicken egg-laying traits indicated the SNP was associated with age at first egg production (AFE) and egg number at 32 and 48 weeks (E32, E48; P < 0.01). Hens with genotype AA had an earlier AFE and greater E32 and E48 than hens with other genotypes. The abundance of mature miR-1b-3p in the hens with the GG genotype was larger than those with the AA genotype (P < 0.01), and the luciferase activity of GG genotype promoter was also greater in birds with the AA than GG genotype (P < 0.05). There was inhibition of the production of the transcription factor bound by the specificity protein 1 (Sp1) as a result of the G-to-A mutation, and the luciferase activity of the GG, but not AA, genotype was markedly increased by Sp1. In conclusion, the SNP, rs737028527 (G> A), affected the abundance of mature miR-1b-3p by Sp1 and was associated with chicken egg-laying traits. Data from the present study allow for an increased understanding of the functions and regulation of miR-1b-3p in ovarian follicle development of hens.

Proteomic analysis of sheep uterus reveals its role in prolificacy

Small Tail Han sheep have attracted attention for their high fecundity and year-round estrus. However, the molecular mechanisms of this fecundity are unknown. Polymorphism of the FecB gene has been shown to be associated with the ovulation rate and litter size in sheep. In this study, we used tandem mass tag quantitative proteomic techniques to identify the differentially abundant proteins in polytocous and monotocous Small Tail Han sheep (FecB++) uterine tissues in the follicular and luteal phases. In total, 41 and 43 differentially abundant proteins were identified in the follicular and luteal phases, respectively. Correlation analysis between the transcriptome and proteome revealed a positive correlation at the two omics levels of prolificacy. GO and KEGG pathway analyses indicated that the mRNAs and proteins upregulated in the polytocous group relative to the monotocous group are involved in sphingolipid metabolism and amino acid metabolism, and may be important in maintaining uterine functions and increasing the embryo survival rate during the estrus cycle of polytocous sheep. In conclusion, our work provides a prospective understanding of the molecular mechanism underlying the high prolificacy of Small Tail Han sheep. SIGNIFICANCE: Fecundity critically affects the profitability of sheep production, but the genetic mechanism of high-prolificacy is still unclear in sheep. We identified potential signaling pathways and differentially abundant proteins associated with reproductive performance through a combination of sheep uterus tissues proteome and transcriptome analyses. These findings will facilitate a better revealing the mechanism and provide possible targets for molecular design breeding for the formation of polytocous traits in sheep.

Comprehensive analysis of long noncoding RNA and mRNA expression patterns in sheep testicular maturation

Long noncoding RNAs (LncRNAs) have been identified as important regulators of testis development; however, their expression patterns and roles in sheep are not yet clear. Thus, we used stranded specific RNA-seq to profile the testis transcriptome (lncRNAs and mRNAs) in premature and mature sheep. Hormone levels and the testis index were examined, and histological analyses were performed at five stages of testis development, 5-day-old (D5), 3-month-old (3M), 6-month-old (6M), 9-month-old (9M), and 2-year-old (2Y), the results of which indicate a significant difference in hormone levels and testis morphometries between the 3M and 9M stages (P < 0.05). Based on the comparison between 3M and 9M samples, we found 1,118 differentially expressed (DE) lncRNAs and 7,253 DE mRNAs in the testes, and qRT-PCR analysis showed that the results correlated well with the transcriptome data. Furthermore, we constructed lncRNA-protein-coding gene interaction networks. Forty-seven DE lncRNA-targeted genes enriched for male reproduction were obtained by cis- and trans-acting; 51 DE lncRNAs and 45 cis-targets, 2 DE lncRNAs and 2 trans-targets were involved in this network. Of these, 5 lncRNAs and their targets, PRKCD, NANOS3, SERPINA5, and CYP19A1, were enriched for spermatogenesis and male gonad development signaling pathways. We further examined the expression levels of 5 candidate lncRNAs and their target genes during testis development. Lastly, the interaction of lncRNA TCONS__00863147 and its target gene PRKCD was validated in vitro in sheep Leydig cells. This study provides a valuable resource for further study of lncRNA function in sheep testis development and spermatogenesis.

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.

Whole genome sequence analysis to detect signatures of positive selection for high fecundity in sheep

Ovulation rate and prolificacy are the most important reproductive traits that have major impact on the efficiency of lamb meat production. Here, we compared the whole genomes of the Romanov sheep, known as one of the high prolific breeds, and four other sheep breeds namely Assaf, Awassi, Cambridge and British du cher, to identify genetic mechanisms underlying prolificacy in sheep. Selection signature analysis revealed 637 and 477 protein-coding genes under positive selection from F_ST and nucleotide diversity (Pi) statistics, respectively. Further analysis showed that several candidate genes including LEPR, PDGFRL and KLF5 genes are involved in sheep prolificacy. The identified candidate genes in the selected regions are novel and provide new insights into the genetic mechanisms underlying prolificacy in sheep and can be useful in sheep breeding programmes to develop improved breeds for high reproductive efficiency.

microRNA-125b Regulates Apoptosis by Targeting Bone Morphogenetic Protein Receptor 1B in Yak Granulosa Cells

The intronic microRNA, miR-125b, plays a vital role in promyelocytic and hematopoietic stem cells, and in the development and apoptosis of cancer cells. In this study, we showed that miR-125b regulates granulosa cell (GC) apoptosis in the yak ovary. Bioinformatic analyses and luciferase reporter assays demonstrated that bone morphogenetic protein receptor type 1B (BMPR1B) is an miR-125b target. miR-125b overexpression induced apoptosis in yak GC, and affected the mRNA and protein expression of BMPR1B and the ratio of Bcl2/Bax. Silencing of miR-125b decreased the rate of yak GC apoptosis and increased the ratio of Bcl2/Bax. In addition, the effects of an miR-125b inhibitor were overturned by cotransfection with siRNA-BMPR1B2 (siRNA-299) in yak GC. Together, these results demonstrated that miR-125b regulates GC apoptosis in the yak ovary by targeting BMPR1B.

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.

Genome-Wide Association Analyses Highlight the Potential for Different Genetic Mechanisms for Litter Size Among Sheep Breeds

Reproduction is an important trait in sheep breeding as well as in other livestock. However, despite its importance the genetic mechanisms of litter size in domestic sheep (Ovis aries) are still poorly understood. To explore genetic mechanisms underlying the variation in litter size, we conducted multiple independent genome-wide association studies in five sheep breeds of high prolificacy (Wadi, Hu, Icelandic, Finnsheep, and Romanov) and one low prolificacy (Texel) using the Ovine Infinium HD BeadChip, respectively. We identified different sets of candidate genes associated with litter size in different breeds: BMPR1B, FBN1, and MMP2 in Wadi; GRIA2, SMAD1, and CTNNB1 in Hu; NCOA1 in Icelandic; INHBB, NF1, FLT1, PTGS2, and PLCB3 in Finnsheep; ESR2 in Romanov and ESR1, GHR, ETS1, MMP15, FLI1, and SPP1 in Texel. Further annotation of genes and bioinformatics analyses revealed that different biological pathways could be involved in the variation in litter size of females: hormone secretion (FSH and LH) in Wadi and Hu, placenta and embryonic lethality in Icelandic, folliculogenesis and LH signaling in Finnsheep, ovulation and preovulatory follicle maturation in Romanov, and estrogen and follicular growth in Texel. Taken together, our results provide new insights into the genetic mechanisms underlying the prolificacy trait in sheep and other mammals, suggesting targets for selection where the aim is to increase prolificacy in breeding projects.

Integrated ovarian mRNA and miRNA transcriptome profiling characterizes the genetic basis of prolificacy traits in sheep (Ovis aries)

BACKGROUND

The highly prolific breeds of domestic sheep (Ovis aries) are globally valuable genetic resources for sheep industry. Genetic, nutritional and other environmental factors affect prolificacy traits in sheep. To improve our knowledge of the sheep prolificacy traits, we conducted mRNA-miRNA integrated profiling of ovarian tissues from two pure breeds with large (Finnsheep) vs. small (Texel) litter sizes and their F1 crosses, half of which were fed a flushing diet.

RESULTS

Among the samples, 16,402 genes (60.6% known ovine genes) were expressed, 79 novel miRNAs were found, and a cluster of miRNAs on chromosome 18 was detected. The majority of the differentially expressed genes between breeds were upregulated in the Texel with low prolificacy, owing to the flushing diet effect, whereas a similar pattern was not detected in the Finnsheep. F1 ewes responded similarly to Finnsheep rather than displaying a performance intermediate between the two pure breeds.

CONCLUSIONS

The identification and characterization of differentially expressed genes and miRNAs in the ovaries of sheep provided insights into genetic and environmental factors affecting prolificacy traits. The three genes (CST6, MEPE and HBB) that were differentially expressed between the group of Finnsheep and Texel ewes kept in normal diet appeared to be candidate genes of prolificacy traits and will require further validation.

miR-26a-5p Regulates TNRC6A Expression and Facilitates Theca Cell Proliferation in Chicken Ovarian Follicles

Ovarian theca cells play an indispensable role in ovarian follicular development and hormone secretion. miR-26a-5p was reported to be differentially expressed in mature and immature chicken ovaries in our previous study; however, the role of miR-26a-5p in regulating ovarian follicle function is still unclear. In this study, we demonstrated that the expression dynamics of TNRC6A mRNA in either chicken ovaries or follicles showed an opposite trend compared with that of chicken miR-26a-5p expression. miR-26a-5p inhibited TNRC6A mRNA expression by directly targeting its 3'-untranslated region in cultured chicken theca cells. Overexpression of miR-26a-5p promoted chicken follicular theca cell proliferation in vitro. Furthermore, overexpression of miR-26a-5p and knockdown of TNRC6A significantly upregulated the antiapoptotic BCL-2 gene. Taken together, this study revealed the expression dynamics of miR-26a-5p and TNRC6A in chicken ovaries and ovarian follicles and the relationship between the expression of miR-26a-5p and TNRC6A in chicken ovarian theca cells. These results suggest that miR-26a-5p facilitates chicken ovarian theca cell proliferation by targeting the TNRC6A gene.