Characterization of microRNAs from goat (Capra hircus) by Solexa deep-sequencing technology

Deep Small RNA Sequencing Reveals Important miRNAs Related to Muscle Development and Intramuscular Fat Deposition in Longissimus dorsi Muscle From Different Goat Breeds

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have been shown to play important post-transcriptional regulatory roles in the growth and development of skeletal muscle tissues. However, limited research into the effect of miRNAs on muscle development in goats has been reported. In this study, Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with significant phenotype difference in meat production performance were selected and the difference in Longissimus dorsi muscle tissue expression profile of miRNAs between the two goat breeds was then compared using small RNA sequencing. A total of 1,623 miRNAs were identified in Longissimus dorsi muscle tissues of the two goat breeds, including 410 known caprine miRNAs, 928 known species-conserved miRNAs and 285 novel miRNAs. Of these, 1,142 were co-expressed in both breeds, while 230 and 251 miRNAs were only expressed in LC and ZB goats, respectively. Compared with ZB goats, 24 up-regulated miRNAs and 135 miRNAs down-regulated were screened in LC goats. A miRNA-mRNA interaction network showed that the differentially expressed miRNAs would target important functional genes associated with muscle development and intramuscular fat deposition. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the target genes of differentially expressed miRNAs were significantly enriched in Ras, Rap 1, FoxO, and Hippo signaling pathways. This study suggested that these differentially expressed miRNAs may be responsible for the phenotype differences in meat production performance between the two goat breeds, thereby providing an improved understanding of the roles of miRNAs in muscle tissue of goats.

The expression and regulation of miR-1 in goat skeletal muscle and satellite cell during muscle growth and development

MicroRNA-1 (miR-1) has been shown to play an important role in muscle growth and development, however, it was mainly discovered in model animals. To explore the function and mechanism of miR-1 in goat, we firstly explored the expression profile of miR-1 in goat tissues and cells. Furthermore, the target gene of miR-1 was predicted, and the relationship between miR-1 and one of its target genes, histone deacetylase 4 (HDAC4), was analyzed through double luciferase reporter assay, real-time PCR, and western blot. It was found that the miR-1 is most abundantly expressed in goat heart and skeletal muscle tissue. Meanwhile, the expression of miR-1 showed an increasing tendency from new-born goats to the 7-month-old goats, and then its expression decreases as the goats mature further. In addition, the expression levels of miR-1 decreased in goat skeletal muscle satellite cells with the algebraic increasing of cells. At last, the results showed that HDAC4 is a target gene of miR-1 in goat, and miR-1 can inhibit the post-transcriptional expression of HDAC4, but had no significant influence on the mRNA level of HDAC4. It was hypothesized that miR-1 promotes muscle development by inhibiting the post-transcriptional expression of HDAC4 in goat.

Distribution of ncRNAs expression across hypothalamic-pituitary-gonadal axis in Capra hircus

Background

Molecular regulation of the hypothalamic-pituitary-gonadal (HPG) axis plays an essential role in the fine tuning of seasonal estrus in Capra hircus. Noncoding RNAs (ncRNAs) are emerging as key regulators in sexual development and mammalian reproduction. In order to identify ncRNAs and to assess their expression patterns, along the HPG axis, we sequenced ncRNA libraries from hypothalamus, pituitary and ovary of three goats.

Results

Among the medium length noncoding RNAs (mncRNAs) identified, small nucleolar RNAs (snoRNAs) and transfer RNAs (tRNAs) were found to be more abundant in ovary and hypothalamus, respectively. The observed GC content was representative for different classes of ncRNAs, allowing the identification of a tRNA-derived RNA fragments (tRFs) subclass, which had a peak distribution around 32–38% GC content in the hypothalamus. Differences observed among organs confirmed the specificity of microRNA (miRNA) profiles for each organ system.

Conclusions

Data on ncRNAs in organs constituting the HPG axis will contribute to understanding their role in the physiological regulation of reproduction in goats.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4767-x) contains supplementary material, which is available to authorized users.

The Binding Sites of miR-619-5p in the mRNAs of Human and Orthologous Genes

Background

Normally, one miRNA interacts with the mRNA of one gene. However, there are miRNAs that can bind to many mRNAs, and one mRNA can be the target of many miRNAs. This significantly complicates the study of the properties of miRNAs and their diagnostic and medical applications.

Results

The search of 2,750 human microRNAs (miRNAs) binding sites in 12,175 mRNAs of human genes using the MirTarget program has been completed. For the binding sites of the miR-619-5p the hybridization free energy of the bonds was equal to 100% of the maximum potential free energy. The mRNAs of 201 human genes have complete complementary binding sites of miR-619-5p in the 3’UTR (214 sites), CDS (3 sites), and 5’UTR (4 sites). The mRNAs of CATAD1, ICA1L, GK5, POLH, and PRR11 genes have six miR-619-5p binding sites, and the mRNAs of OPA3 and CYP20A1 genes have eight and ten binding sites, respectively. All of these miR-619-5p binding sites are located in the 3’UTRs. The miR-619-5p binding site in the 5’UTR of mRNA of human USP29 gene is found in the mRNAs of orthologous genes of primates. Binding sites of miR-619-5p in the coding regions of mRNAs of C8H8orf44, C8orf44, and ISY1 genes encode the WLMPVIP oligopeptide, which is present in the orthologous proteins. Binding sites of miR-619-5p in the mRNAs of transcription factor genes ZNF429 and ZNF429 encode the AHACNP oligopeptide in another reading frame. Binding sites of miR-619-5p in the 3’UTRs of all human target genes are also present in the 3’UTRs of orthologous genes of mammals. The completely complementary binding sites for miR-619-5p are conservative in the orthologous mammalian genes.

Conclusions

The majority of miR-619-5p binding sites are located in the 3’UTRs but some genes have miRNA binding sites in the 5’UTRs of mRNAs. Several genes have binding sites for miRNAs in the CDSs that are read in different open reading frames. Identical nucleotide sequences of binding sites encode different amino acids in different proteins. The binding sites of miR-619-5p in 3’UTRs, 5’UTRs and CDSs are conservative in the orthologous mammalian genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3811-6) contains supplementary material, which is available to authorized users.

Computational identification and characterization of novel microRNA in the mammary gland of dairy goat (Capra hircus)

Many studies have indicated that microRNAs (miRNAs) influence the development of the mammary gland by posttranscriptionally affecting their target genes. The objective of this research was to identify novel miRNAs in the mammary gland of dairy goats with a bioinformatics approach that was based on expressed sequence tag (EST) and genome survey sequence (GSS) analyses. We applied all known major mammals, miRNAs to search against the goat EST and GSS databases for the first time to identify new miRNAs. We, then, validated these newly predicted miRNAs with stem-loop reverse transcription followed by a SYBR Green polymerase chain reaction assay. Finally, 29 mature miRNAs were identified and verified, and of these, 14 were grouped into 13 families based on seed sequence identity and 85 potential target genes of newly verified miRNAs were subsequently predicted, most of which seemed to encode the proteins participating in regulation of metabolism, signal transduction, growth and development. The predicting accuracy of the new miRNAs was 70.37%, which confirmed that the methods used in this study were efficient and reliable. Detailed analyses of the sequence characteristics of the novel miRNAs of the goat mammary gland were performed. In conclusion, these results provide a reference for further identification of miRNAs in animals without a complete genome and thus improve the understanding of miRNAs in the caprine mammary gland.

Comparative profiling of differentially expressed microRNAs between the follicular and luteal phases ovaries of goats

To explore if the regulation at post-transcriptional level of follicular phase (Fols) to luteal phase (Luts) transition occurs in the ovaries of Anhuai goats, the differentially expressed microRNAs (miRNAs) of ovaries in the Fols and Luts were analyzed using Solexa sequencing in the study. In total, 320 known miRNAs were co-expressed in the two phases, 339 and 353 known miRNAs were expressed in the ovary in the Fols and Luts, respectively. In addition, 45 novel miRNAs were co-expressed in the two phases, 70 and 94 novel miRNAs were expressed in the ovary in the Fols and Luts, respectively. Let-7f was the highest expressed significantly different known miRNA in the two phases, and mir-159 was the highest expressed significantly different novel miRNA in the two phases, which may participate in the follicular-luteal transition of Anhuai goats. GO annotation and KEGG pathway analysis were applied to analyze the target genes of differentially expressed miRNAs detected in the two phases. The results will help to further understand the role of miRNAs in the regulation of follicular to luteal transition in goat ovaries.

The properties of binding sites of miR-619-5p, miR-5095, miR-5096, and miR-5585-3p in the mRNAs of human genes

The binding of 2,578 human miRNAs with the mRNAs of 12,175 human genes was studied. It was established that miR-619-5p, miR-5095, miR-5096, and miR-5585-3p bind with high affinity to the mRNAs of the 1215, 832, 725, and 655 genes, respectively. These unique miRNAs have binding sites in the coding sequences and untranslated regions of mRNAs. The mRNAs of many genes have multiple miR-619-5p, miR-5095, miR-5096, and miR-5585-3p binding sites. Groups of mRNAs in which the ordering of the miR-619-5p, miR-5095, miR-5096, and miR-5585-3p binding sites differ were established. The possible functional and evolutional properties of unique miRNAs are discussed.

Identification and profiling of microRNAs and their target genes from developing caprine skeletal Muscle

Goat is an important agricultural animal for meat production. Functional studies have demonstrated that microRNAs (miRNAs) regulate gene expression at the post-transcriptional level and play an important role in various biological processes. Although studies on miRNAs expression profiles have been performed in various animals, relatively limited information about goat muscle miRNAs has been reported. To investigate the miRNAs involved in regulating different periods of skeletal muscle development, we herein performed a comprehensive research for expression profiles of caprine miRNAs during two developmental stages of skeletal muscles: fetal stage and six month-old stage. As a result, 15,627,457 and 15,593,721 clean reads were obtained from the fetal goat library (FC) and the six month old goat library (SMC), respectively. 464 known miRNAs and 83 novel miRNA candidates were identified. Furthermore, by comparing the miRNA profile, 336 differentially expressed miRNAs were identified and then the potential targets of the differentially expressed miRNAs were predicted. To understand the regulatory network of miRNAs during muscle development, the mRNA expression profiles for the two development stages were characterized and 7322 differentially expressed genes (DEGs) were identified. Then the potential targets of miRNAs were compared to the DEGs, the intersection of the two gene sets were screened out and called differentially expressed targets (DE-targets), which were involved in 231 pathways. Ten of the 231 pathways that have smallest P-value were shown as network figures. Based on the analysis of pathways and networks, we found that miR-424-5p and miR-29a might have important regulatory effect on muscle development, which needed to be further studied. This study provided the first global view of the miRNAs in caprine muscle tissues. Our results help elucidation of complex regulatory networks between miRNAs and mRNAs and for the study of muscle development.

Identification and characterization of microRNAs in the ovaries of multiple and uniparous goats (Capra hircus) during follicular phase

Background

Superior kidding rate is an important economic trait in production of meat goat, and ovulation rate is the precondition of kidding rate. MicroRNAs (miRNAs) play critical roles in almost all ovarian biological processes, including folliculogenesis, follicle development, follicle atresia, luteal development and regression. To find out the different ovarian activity and follicle recruitment with miRNA-mediated posttranscriptional regulation, the small RNAs expressed pattern in the ovarian tissues of multiple and uniparous Anhui White goats during follicular phase was analyzed using Solexa sequencing data.

Results

1008 miRNAs co-expressed, 309 and 433 miRNAs specifically expressed in the ovaries of multiple and uniparous goats during follicular phase were identified. The 10 most highly expressed miRNAs in the multiple library were also the highest expressed in the uniparous library, and there were no significantly different between each other. The highest specific expressed miRNA in the multiple library was miR-29c, and the one in the uniparous library was miR-6406. 35 novel miRNAs were predicted in total. GO annotation and KEGG Pathway analyses were implemented on target genes of all miRNA in two libraries. RT-PCR was applied to detect the expression level of 5 randomly selected miRNAs in multiple and uniparous hircine ovaries, and the results were consistent with the Solexa sequencing data.

Conclusions

In the present study, the different expression of miRNAs in the ovaries of multiple and uniparous goats during follicular phase were characterized and investigated using deep sequencing technology. The result will help to further understand the role of miRNAs in kidding rate regulation and also may help to identify miRNAs which could be potentially used to increase hircine ovulation rate and kidding rate in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-339) contains supplementary material, which is available to authorized users.