miRNA-mRNA analysis of sheep adrenal glands reveals the network regulating reproduction

Abstract

Background: The adrenal gland participates in the process of sheep reproduction. MicroRNAs (miRNAs), endogenous small noncoding RNAs, regulate gene expression at the posttranscriptional level. However, the miRNA-mRNA network profile of adrenal glands relating to reproduction in sheep is still not well-studied. As sheep with FecBBB genotype show higher lambing number compare with the sheep with FecB++ genotype. This research aims to compare gene expression by small RNA-seq in adrenal tissues at follicular (F) and luteal (L) phases in FecBBB (MM) and FecB++ (ww) sheep. After analysis of gene expression, significant differentially expressed microRNAs (DEMs) and corresponding target genes were identified. Results: A total of 180 miRNAs were found in this study, of which 19 DEMs were expressed in the four comparison groups (MM_F_A vs. MM_L_A, MM_F_A vs. ww_F_A, MM_L_A vs. ww_L_A, ww_F_A vs. ww_L_A). Subsequently, 354 target genes of 19 DEMs were predicted by integrated analysis. Cluster analysis was performed by K_means_cluster, and the expression patterns of these DEMs were separated into four subclusters. Functional analysis of target genes was performed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The results indicated that the target genes were involved mainly in the Notch signaling pathway, signal transduction, cell communication, innate immune response and amino acid metabolism. Specifically, the Notch signaling pathway, biosynthetic process and metabolic process of pyrimidine nucleotide and amino acid metabolism appear to play key regulatory roles in the sheep fertility trait. Furthermore, miRNA-mRNA interaction networks were constructed by differentially expressed genes combined with our previous study of transcriptome data. The results showed that several key genes, including TDRD3, ANAPC7, CCNL2, BRD2 and MUT, were related to the transformation from the follicular phase to the luteal phase. PLAC8L1, NFAT5, DDX24 and MBD1 were related to the high fecundity of small tail Han sheep. Conclusions: In this study, the miRNA transcriptome profile was identified, and miRNA-mRNA interaction networks were constructed in adrenal gland tissue of small tail Han sheep, the interaction between miR-370-3p and its targets were considered to play a major role in the reproduction regulation process. The results enriched the number of known miRNAs in adrenal glands and provided novel ideas and further information to demonstrate posttranscriptional regulation mechanisms at follicular and luteal phases in different genotypes of small tail Han sheep.