Identification of miRNAs in Bovine Endometrium through RNAseq and Prediction of Regulated Pathways

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.

Identification of genes associated with susceptibility to persistent breeding-induced endometritis by RNA-sequencing of uterine cytobrush samples

This study aimed to investigate the susceptibility to persistent breeding-induced endometritis (PBIE). Cytobrush samples were collected from 81 broodmares 1-3 days before artificial insemination (AI). Susceptibility to PBIE was evaluated by the presence of ≥ 2 cm of intrauterine fluid 24 h after AI, besides the fertility was determined by a sonographic pregnancy diagnosis 2 weeks after ovulation. RNA expressions were compared between susceptible non-pregnant (SNP) mares (n=9) and resistant pregnant (RP) mares (n=9) as well as between susceptible pregnant (SP) mares (n=9) and susceptible non-pregnant (SNP) mares. 66 differentially expressed genes (DEGs) were identified between SNP and RP mares and 60 DEGs between SP and SNP mares. In SNP compared to RP mares, transcript levels of genes regulating steroid hormone metabolism and neutrophil chemotaxis were lower, while higher for genes participating in uterine inflammation.Transcripts of genes related to extracellular matrix degradation, tissue adhesions, and fibrosis were lower in SP mares than in SNP mares, while higher for genes related to uterine cell proliferation, differentiation, and angiogenesis in SP mares than SNP mares. In conclusion, increased transcript levels of apolipoprotein E (APOE) and roundabout 2 (ROBO2), cluster domain 44 (CD44), integrin beta 3 (ITGB3), and epidermal growth factor (EGF) are possible biomarkers for susceptibility to PBIE. While higher expression of fibroblast growth factor 9 (FGF9), kinase domain receptor (KDR), and C-X-C motif chemokine ligand (CXCL) 16, collagen type V alpha 2 (COL5A2) and fibronectin (FN1) are suggested indicators of fertility in susceptible mares if they receive proper breeding management.

A comparative analysis of gene expression induced by the embryo in the caprine endometrium

Transcriptomics is an established powerful tool to identify potential mRNAs and ncRNAs (non‐coding RNAs) for endometrial receptivity. In this study, the goat endometrium at estrus day 5 (ED5) and estrus day 15 (ED15) were selected to systematically analyse the differential expressed genes (DEGs) what were induced by the embryo. There were 1,847 genes which were significantly differential expressed in endometrium induced by the embryo at ED5, and 1,346 at ED15 (p‐value < .05). Secreted phosphoprotein 1 (SPP) was the responsive genes for embryo in the goat endometrium during estrus cycle, neurotensis (NTS) and pleiotrophin (PTN) were the responsive genes for embryo in the goat endometrium at ED5, Testin (TES) and Phosphate and Tension Homology Deleted on Chromsome ten (PTEN) at ED15. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes (KEGG) analysis revealed cytoplasm and Endocytosis were indispensable for the endometrium development in dairy goat. In a word, this resulting view of the transcriptome greatly uncovered the global trends in mRNAs expression induced by the embryo in the endometrium of dairy goats.

The mRNA and lncRNA landscape of the non-pregnant endometrium during the oestrus cycle in dairy goat

Cyclic changes in the endometrium are essential for embryo implantation in mammals; many studies report that such changes constitute a complex process involving numerous molecular mediators. In the present study, goat endometria at oestrus Day 5 and oestrus Day 15 were selected to systematically analyse the transcriptome using strand-specific Ribo-Zero RNA sequencing. Over 120 million high-quality paired-end reads were generated and 440400 transcripts were identified in the endometrial tissue of dairy goats. In total, 489 differentially expressed mRNAs and 854 differentially expressed long non-coding RNAs were identified when comparing the endometrium at goat endometria at oestrus Day 5 and oestrus Day 15. Neurotensin was found to play a potentially important role in the non-pregnant goat endometrium during the oestrus cycle. Furthermore, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes analyses of the cis-target genes of the differentially expressed long non-coding RNAs showed that GO:0005198 (structural molecule activity) and ko04510 (focal adhesion) might be involved in cyclic endometrial changes. Taken together, the resulting transcriptomic profiles elucidate global trends in mRNA and lncRNA expression in non-pregnant endometria during the oestrus cycle in dairy goats.

Protective effects of microRNA-330 on amyloid β-protein production, oxidative stress, and mitochondrial dysfunction in Alzheimer's disease by targeting VAV1 via the MAPK signaling pathway

This study aims to explore the effect of miR-330 targeting VAV1 on amyloid β-protein (Aβ) production, oxidative stress (OS), and mitochondrial dysfunction in Alzheimer's disease (AD) mice through the MAPK signaling pathway. Putative targeted gene of miR-330 was performed by a miRNA target prediction website and dual-luciferase reporter gene assay. AD mouse model was successfully established. Fourteen C57 mice were randomized into AD and control groups. The positive protein expression rate of VAV1 was measured by immunohistochemistry. Neuron cells were assigned into control, blank, negative control (NC), miR-330 mimics, miR-330 inhibitors, siRNA-VAV1, and miR-330 inhibitors + siRNA-VAV1 groups. Expression of miR-330, VAV1, ERK1, JNK1, P38MAPK, Aβ, COX, and lipoprotein receptor-related protein-1 (LRP-1) were determined using RT-qPCR and Western blotting. Colorimetry was applied to measure the levels of OS parameters of superoxide dismutase (SOD) and malondialdehyde (MDA). Aβ production in brain tissue was detected using ELISA, while that in neuron cell was measured by radioimmunoassay. MiR-330 was down-regulated in neuron cells of AD mice and VAV1 was negatively regulated by miR-330. Compared with the control group, the positive protein expression rate of VAV1 was significantly elevated in the AD group. Overexpression of miR-330 decreased the expression of VAV1, ERK1, JNK1, P38MAPK, and Aβ, but increased the expression of COX and LRP-1. AD mice revealed elevated Aβ production and MDA with decreased SOD level. The result indicates that overexpressed miR-330 targeting VAV1 through the MAPK signaling pathway reduces Aβ production and alleviates OS and mitochondrial dysfunction in AD.

The developmental transcriptome landscape of receptive endometrium during embryo implantation in dairy goats

Under natural conditions, some embryos cannot implant successfully because of the dysfunction of receptive endometrium (RE). Thus, it is imperative for us to study the molecular mechanisms involved in the formation of the RE from pre-receptive endometrium (PE). In this study, the endometrium from gestational day 5 (D5, PE) and gestational day 15 (D15, RE) dairy goats were selected to systematically analyze the transcriptome using strand-specific Ribo-Zero RNA-Seq, >120 million high-quality paired-end reads were generated and 47,616 transcripts were identified in the endometrium of dairy goats. A total of 810 mRNAs were differentially expressed genes (DEGs) between the RE and PE meeting the criteria of P-values<0.05. Bioinformatics analysis of the DEGs revealed that a number of biological processes and pathways were potentially involved in the establishment of the RE, notably energy metabolism and amino acid metabolism. Furthermore, we speculated that CXCL14, IGFBP3, and LGALS15 potentially participated in the development of endometrium. What's more, putative SNPs, InDels and AS events were identified and analyzed in the endometrium. In a word, this resulting view of the transcriptome greatly enhances the comprehensive transcript catalog and uncovers the global trends in gene expression during the formation of receptive endometrium in dairy goats.

Cytobrush: A tool for sequential evaluation of gene expression in bovine endometrium

Aims were to (i) compare specific transcript abundance between endometrial samples collected by transcervical biopsy and cytobrush and (ii) measure the abundance of endometrial transcripts involved in PGF2α synthesis in samples collected by cytobrush. In Experiment 1, endometrial samples were taken transcervically by cytobrush and biopsy 10 days after ovulation. Compared to biopsy samples, abundance of transcripts for MSTN, AKR1C4 and PGR was similar, VIM, FLT1 and PTGES was lower (p < .05) and KRT18 and CD3D was greater in cytobrush samples (p < .05). Thus, there was an enrichment of epithelial and immune cells in the cytobrush samples. In Experiment 2, endometrial samples were collected by cytobrush on days 10, 13, 16 and 19 after ovulation. Abundance of PGR2 mRNA was maximum on day 10 then decreased (p < .05). Abundance of ESR1 decreased gradually from day 10 to day 16 then increased again on day 19. The greatest abundance of OXTR was noted on day 19. The sequential alterations in abundance of these transcripts are consistent with the release of PGF2α associated with luteolysis. In summary, cytobrush sampling provides representative, physiologically relevant samples of the luminal epithelium in cattle.

MicroRNA Mediating Networks in Granulosa Cells Associated with Ovarian Follicular Development

Ovaries, which provide a place for follicular development and oocyte maturation, are important organs in female mammals. Follicular development is complicated physiological progress mediated by various regulatory factors including microRNAs (miRNAs). To demonstrate the role of miRNAs in follicular development, this study analyzed the expression patterns of miRNAs in granulosa cells through investigating three previous datasets generated by Illumina miRNA deep sequencing. Furthermore, via bioinformatic analyses, we dissected the associated functional networks of the observed significant miRNAs, in terms of interacting with signal pathways and transcription factors. During the growth and selection of dominant follicles, 15 dysregulated miRNAs and 139 associated pathways were screened out. In comparison of different styles of follicles, 7 commonly abundant miRNAs and 195 pathways, as well as 10 differentially expressed miRNAs and 117 pathways in dominant follicles in comparison with subordinate follicles, were collected. Furthermore, SMAD2 was identified as a hub factor in regulating follicular development. The regulation of miR-26a/b on smad2 messenger RNA has been further testified by real time PCR. In conclusion, we established functional networks which play critical roles in follicular development including pivotal miRNAs, pathways, and transcription factors, which contributed to the further investigation about miRNAs associated with mammalian follicular development.

Circulating microRNA Profiles during the Bovine Oestrous Cycle

Up to 50% of ovulations go undetected in modern dairy herds due to attenuated oestrus behavior and a lack of high-accuracy methods for detection of fertile oestrus. This significantly reduces overall herd productivity and constitutes a high economic burden to the dairy industry. MicroRNAs (miRNAs) are ubiquitous regulators of gene expression during both health and disease and they have been shown to regulate different reproductive processes. Extracellular miRNAs are stable and can provide useful biomarkers of tissue function; changes in circulating miRNA profiles have been reported during menstrual cycles. This study sought to establish the potential of circulating miRNAs as biomarkers of oestrus in cattle. We collected plasma samples from 8 Holstein-Friesian heifers on days Days 0, 8 and 16 of an oestrous cycle and analysed small RNA populations on each Day using two independent high-throughput approaches, namely, Illumina sequencing (n = 24 samples) and Qiagen PCR arrays (n = 9 sample pools, 3–4 samples / pool). Subsequently, we used RT-qPCR (n = 24 samples) to validate the results of high-throughput analyses, as well as to establish the expression profiles of additional miRNAs previously reported to be differentially expressed during reproductive cycles. Overall, we identified four miRNAs (let-7f, miR-125b, miR-145 and miR-99a-5p), the plasma levels of which distinctly increased (up to 2.2-fold, P < 0.05) during oestrus (Day 0) relative to other stages of the cycle (Days 8 and 16). Moreover, we identified several hundred different isomiRs and established their relative abundance in bovine plasma. In summary, our results reveal the dynamic nature of plasma miRNAs during the oestrous cycle and provide evidence of the feasibility of using circulating miRNAs as biomarkers of reproductive function in livestock in the future.

Experimental and bioinformatic analysis of cultured Bovine Endometrial Cells (BEND) responding to interferon tau (IFNT)

BACKGROUND

In ruminants, embryo implantation depends on progesterone (P4) and interferon tau (IFNT) controlling endometrial function. IFNT antagonizes bovine endometrial cells (BEND) response to phorbol 12,13-dibutyrate (PDBU) through posttranscriptional regulation of gene expression. We have previously described microRNAs (miRNAs) profiles in bovine endometrium, detecting miR-106a, relevant for embryo maternal communication. In this study, we investigated the expression miR-106a and genes for prostaglandin-endoperoxide synthase 2 (PTGS2), phospholipase A2, group IVA (PLA2G4A), estrogen receptor 1 (ESR1) and progesterone receptor (PR) in response to IFNT in BEND cells and searched for interferon responsive factors (IRFs) binding sites in their promoter genomic regions. The aim of this study was to unravel the molecular mechanisms involved in IFNT signalling and its regulation of miR-106a.

FINDINGS

PTGS2 showed increased expression under PDBU, which was antagonized by IFNT. IFNT induced expression of PR and miR-106a and downregulation of ESR1 and PR. Bioinformatic analyses detected that PLA2G4A was associated to IRF-1 and IRF-6, while ESR1, PR and PTGS2 were associated to only IRF-6. All genes exhibit one motif per IRF, except miR-106a that had three binding sites for IRF-6.

CONCLUSIONS

We report the IFNT regulatory effect on miR-106a expression through IRF-6 in bovine endometrial cells. We identified a set of potential binding sites for IRF-1 and IRF-6 within the bovine genome. A set of candidate gene regions could be characterized where IFNT can act via IRFs to regulate the expression of proteins and miRNAs. Future studies will use these data to detect new IFNT regulatory mechanisms in the endometrium.