Discovery of porcine microRNAs in multiple tissues by a Solexa deep sequencing approach

Improvement, identification, and target prediction for miRNAs in the porcine genome by using massive, public high-throughput sequencing data

Despite the broad variety of available microRNA (miRNA) research tools and methods, their application to the identification, annotation, and target prediction of miRNAs in nonmodel organisms is still limited. In this study, we collected nearly all public sRNA-seq data to improve the annotation for known miRNAs and identify novel miRNAs that have not been annotated in pigs (Sus scrofa). We newly annotated 210 mature sequences in known miRNAs and found that 43 of the known miRNA precursors were problematic due to redundant/missing annotations or incorrect sequences. We also predicted 811 novel miRNAs with high confidence, which was twice the current number of known miRNAs for pigs in miRBase. In addition, we proposed a correlation-based strategy to predict target genes for miRNAs by using a large amount of sRNA-seq and RNA-seq data. We found that the correlation-based strategy provided additional evidence of expression compared with traditional target prediction methods. The correlation-based strategy also identified the regulatory pairs that were controlled by nonbinding sites with a particular pattern, which provided abundant complementarity for studying the mechanism of miRNAs that regulate gene expression. In summary, our study improved the annotation of known miRNAs, identified a large number of novel miRNAs, and predicted target genes for all pig miRNAs by using massive public data. This large data-based strategy is also applicable for other nonmodel organisms with incomplete annotation information.

Differential Regulatory Roles of MicroRNAs in Porcine Intramuscular and Subcutaneous Adipocytes

The deposition of intramuscular (IM) and subcutaneous (SC) fat is an important trait influencing pork quality. Understanding the genetic differences between these two types of adipose tissues is consequently of great importance for pig breeding. Here, we established primary cultures of IM and SC adipocytes from Jiaxing black pigs. The microRNA (miRNA) expression profiles of the two types of adipocytes were obtained by RNA-seq. A total of 741 miRNAs were identified in IM and SC adipocytes, including 155 significant differentially expressed (SDE) miRNAs. According to gene ontology and Kyoto Encyclopedia of Genes analysis, the target genes of the SDE miRNAs were enriched in categories and pathways related to transcriptional regulation, fatty acid biosynthesis, as well as the MAPK and PI3K/Akt pathways. Notably, miR-206 expression was 36-fold higher in IM adipocytes than in SC adipocytes. The overexpression of miR-206 in IM and SC adipocytes decreased cell proliferation and triglyceride accumulation. Luciferase activity assays and quantitative polymerase chain reaction confirmed that miR-206 regulates adipocyte proliferation by targeting STARD7 and inhibits adipogenesis by repressing Krüppel-like factor 4 (KLF4) expression. Accordingly, the effect of miR-206 mimics was attenuated by the overexpression of KLF4 in adipocytes. Taken together, we identified the expression profiles of miRNAs in adipocytes, which revealed that miR-206 acts as a suppressor of adipogenesis.

The Effect of MicroRNA-331-3p on Preadipocytes Proliferation and Differentiation and Fatty Acid Accumulation in Laiwu Pigs

Objective

The proliferation and differentiation of preadipocytes are regulated by microRNAs (miRNAs), hormones, and other factors. This study aimed to investigate the effects of miR-331-3p on the proliferation and differentiation of preadipocytes in addition to fatty acid metabolism.

Methods

Preadipocytes were transfected with miR-331-3p mimics, miR-NC, or miR-331-3p inhibitor to explore its effect on cell proliferation and fatty acid accumulation. Furthermore, preadipocytes were transfected with pre-miR-331-3p, pcDNA3.1(+), or miR-331-3p inhibitor to explore its effect on differentiation.

Results

It was observed that miR-331-3p could inhibit preadipocytes proliferation. Furthermore, miR-331-3p was highly expressed during cellular differentiation and appeared to promote the process. In addition, dual fluorescein analysis showed that dihydrolipoamide S-succinyltransferase (DLST) is a target gene of miR-331-3p, and overexpression of miR-331-3p could regulate the metabolism of fatty acids in the citrate pyruvate cycle by targeting DLST expression.

Conclusion

In summary, these findings indicated that miR-331-3p exerts contrasting effects on the processes of fat deposition.

Overexpression of MicroRNA-10a in Germ Cells Causes Male Infertility by Targeting Rad51 in Mouse and Human

Spermatogenesis is a complicated process including spermatogonial stem cells self-renewal and differentiates into mature spermatozoa. MicroRNAs (miRNAs) as a class of small non-coding RNAs play a crucial role during the process of spermatogenesis. However, the function of a plenty of miRNAs on spermatogenesis and the potential mechanisms remain largely unknown. Here, we show that genetically conditional overexpressed miR-10a in germ cells caused complete male sterility, characterized by meiotic arrested in germ cells. Analysis of miR-10a overexpression mouse testes reveals that failure of double strand break (DSB) repairs and aberrant spermatogonial differentiation. Furthermore, we identified Rad51 as a key target of miR-10a in germ cell by bioinformatics prediction and luciferase assay, which may be responsible for the infertility of the miR-10a overexpressed mice and germ cell arrested patients. Our data show that miR-10a dependent genetic regulation of meiotic process is crucial for male germ cell development and spermatogenesis in both mouse and human. These findings facilitate our understanding of the roles of miRNA-10a in spermatogenesis and male fertility.

MicroRNAs regulate the main events in rice drought stress response by manipulating the water supply to shoots

MicroRNAs (miRNAs) are small endogenous regulatory RNAs that are involved in a variety of biological processes related to proliferation, development, and response to biotic and abiotic stresses. miRNA profiles of rice (Oryza sativa L. cv. IR64.) leaves in a partial root zone drying (PRD) system were analysed using a high-throughput sequencing approach to identify miRNAs associated with drought signalling. The treatments performed in this study were as follows: well-watered ("wet" roots, WW), wherein both halves of the pot were watered daily; drought ("dry" roots, DD), wherein water was withheld from both halves of the pot; and well-watered/drought ("wet" and "dry" roots, WD), wherein one half of each pot was watered daily, the same as in WW, and water was withheld from the other part, the same as in DD. High-throughput sequencing enabled us to detect novel miRNAs and study the differential expression of known miRNAs. A total of 209 novel miRNAs were detected in this study. Differential miRNA profiling of the DD, WD and WW conditions showed differential expression of 159 miRNAs, among which 83, 44 and 32 miRNAs showed differential expression under both DD and WD conditions. The detection of putative targets of the differentially expressed miRNAs and investigation of their functions showed that most of these genes encode transcription factors involved in growth and development, leaf morphology, regulation of hormonal homeostasis, and stress response. The most important differences between the DD and WD conditions involved regulation of the levels of hormones such as auxin, cytokinin, abscisic acid, and jasmonic acid and also regulation of phosphor homeostasis. Overall, differentially expressed miRNAs under WD conditions were found to differ from those under DD conditions, with such differences playing a role in adaptation and inducing the normal condition. The mechanisms involved in regulating hormonal homeostasis and involved in energy production and consumption were found to be the most important regulatory pathways distinguishing the DD and WD conditions.

Identification and characterization of microRNAs in the pituitary of pubescent goats

BACKGROUND

Puberty is the period during a female mammal's life when it enters estrus and ovulates for the first time; this indicates that a mammal is capable of reproduction. The onset of puberty is a complex and tightly coordinated biological event; it has been reported that microRNAs (miRNAs) are involved in regulating the initiation of puberty.

METHODS

We performed miRNA sequencing on pituitary tissue from prepubescent and pubescent goats to investigate differences in miRNA expression during the onset of puberty in female goats. The target genes of these miRNAs were evaluated by GO enrichment and KEGG pathway analysis to identify critical pathways regulated by these miRNAs during puberty in goats. Finally, we selected four known miRNA and one novel miRNAs to evaluate expression patterns in two samples via qRT-PCR to validate the RNA-seq data.

RESULTS

In this study, 476 miRNAs were detected in goat pituitary tissue; 13 of these were specifically expressed in the pituitary of prepubescent goats, and 17 were unique to the pituitary of pubescent goats. Additionally, 73 novel miRNAs were predicted in these two libraries. 20 differentially expressed miRNAs were identified in this study. KEGG pathway enrichment analysis revealed that the differentially expressed miRNA target genes were enriched in pathways related to ovary development during puberty, including the GABAergic synapse, oxytocin signaling pathway, the cAMP signaling pathway, progesterone-mediated oocyte maturation. In this study, differential miRNA expression in the pituitary tissue of prepubescent and pubescent goats were identified and characterized.

CONCLUSION

These results provide important information regarding the potential regulation of the onset of goat puberty by miRNAs, and contribute to the elucidation of miRNA regulated processes during maturation and reproduction.

Expression and Regulation Profile of Mature MicroRNA in the Pig: Relevance to Xenotransplantation

The pig is an important source of meat production and provides a valuable model for certain human diseases. MicroRNA (miRNA), which is noncoding RNA and regulates gene expression at the posttranscriptional level, plays a critical role in various biological processes. Studies on identification and function of mature miRNAs in multiple pig tissues are increasing, yet the literature is limited. Therefore, we reviewed current research to determine the miRNAs expressed in specific pig tissues that are involved in carcass values (including muscle and adipocytes), reproduction (including pituitary, testis, and ovary), and development of some solid organs (e.g., brain, lung, kidney, and liver). We also discuss the possible regulating mechanisms of miRNA. Finally, as pig organs are suitable candidates for xenotransplantation, biomarkers of their miRNA in xenotransplantation were evaluated.

High-throughput sequencing analysis of differentially expressed miRNAs and target genes in ischemia/reperfusion injury and apelin-13 neuroprotection

miRNAs regulate a variety of biological processes through pairing-based regulation of gene expression at the 3′ end of the noncoding region of the target miRNA. miRNAs were found to be abnormally expressed in ischemia/reperfusion injury models. High-throughput sequencing is a recently developed method for sequencing miRNAs and has been widely used in the analysis of miRNAs. In this study, ischemia/reperfusion injury models were intracerebroventricularly injected with 50 μg/kg apelin-13. High-throughput sequencing showed that 357 known miRNAs were differentially expressed among rat models, among which 78 changed to > 2-fold or < 0.5-fold. Quantitative real-time polymerase chain reaction was selected to confirm the expression levels of four miRNAs that were differentially expressed, the results of which were consistent with the results of high-throughput sequencing. Gene Ontology analysis revealed that the predicted targets of the different miRNAs are particularly associated with cellular process, metabolic process, single-organism process, cell, and binding. Kyoto Encyclopedia of Gene and Genome analysis showed that the target genes are involved in metabolic pathways, mitogen-activated protein kinase signaling pathway, calcium signaling pathway, and nuclear factor-κB signaling pathway. Our findings suggest that differentially expressed miRNAs and their target genes play an important role in ischemia/reperfusion injury and neuroprotection by apelin-13.

Profiling and characterization of a longissimus dorsi muscle microRNA dataset from an F2 Duroc × Pietrain pig resource population

To elucidate the effects of microRNA (miRNA) regulation in skeletal muscle of adult pigs, miRNA expression profiling was performed with RNA extracted from longissimus dorsi (LD) muscle samples from 174 F₂ pigs (~ 5.5 months of age) from a Duroc × Pietrain resource population. Total RNA was extracted from LD samples, and libraries were sequenced on an Illumina HiSeq 2500 platform in 1 × 50 bp format. After processing, 232,826,977 total reads were aligned to the Sus scrofa reference genome (v10.2.79), with 74.8% of total reads mapping successfully. The miRDeep2 software package was utilized to quantify annotated Sus scrofa mature miRNAs from miRBase (Release 21) and to predict candidate novel miRNA precursors. Among the retained 295 normalized mature miRNA expression profiles ssc­miR­1, ssc­miR­133a­3p, ssc­miR­378, ssc­miR­206, and ssc­miR­10b were the most abundant, all of which have previously been shown to be expressed in pig skeletal muscle. Additionally, 27 unique candidate novel miRNA precursors were identified exhibiting homologous sequence to annotated human miRNAs. The composition of classes of small RNA present in this dataset was also characterized; while the majority of unique expressed sequence tags were not annotated in any of the queried databases, the most abundantly expressed class of small RNA in this dataset was miRNAs. This data provides a resource to evaluate miRNA regulation of gene expression and effects on complex trait phenotypes in adult pig skeletal muscle. The raw sequencing data were deposited in the Sequence Read Archive, BioProject PRJNA363073.

Identification and analysis of the expression of microRNA from lactating and nonlactating mammary glands of the Chinese swamp buffalo

To study the role of microRNA (miR) in the lactation physiology of water buffalo, 2 multiparous dairy buffaloes (including an 8-yr-old buffalo that had been lactating for 3 mo, as well as a 10-yr-old nonlactating, nonpregnant buffalo) were used for miR library construction. The profile of differentially expressed miR in lactating and nonlactating mammary gland tissues of these water buffalo were investigated using Illumina-Solexa high-throughput sequencing technology (Illumina, San Diego, CA). The data identified 259 miR families, 359 mature miR, 363 pre-miR, 230 novel buffalo miR, and 5 buffalo-specific miR that were expressed in mammary tissues. Some highly significantly differentially expressed miR were explored, including bbu-miR-497, bbu-miR-30a-5p, bbu-miR-148a, bbu-miR-29a, bbu-miR-125a, bbu-miR-125b, and bbu-miR-103. The expression patterns of 18 miR were confirmed by quantitative real-time PCR in both tissues, and the expression of bbu miR-103 and novel miR-57 constituted the largest differences between lactating and nonlactating tissues. Further functional analysis indicated that the overexpression or suppression of miR-103 in buffalo mammary epithelial cells downregulated or upregulated the expression of pantothenate kinase 3, and also significantly increased the transcription factor steroid regulatory element binding protein, followed by the acceleration of de novo synthesis of fatty acids by upregulation of acetyl coenzyme A carboxylase α expression. The expression levels of 34 predicted target genes of novel-miR-57 in lactating and nonlactating mammary gland tissues were all analyzed by quantitative real-time PCR. Finally, only the expression of docking protein 4 could be upregulated or downregulated selectively by bbu-novel-miR-57 in buffalo mammary epithelial cells and the Bcap-37 cell line. This study provides an overview of the miR expression profile of water buffalo and the interaction between some key miR and their target genes, which may improve understanding of the important roles of miR in buffalo milk fat synthesis.

Transcriptome-Wide Comparative Analysis of microRNA Profiles in the Telogen Skins of Liaoning Cashmere Goats (Capra hircus) and Fine-Wool Sheep (Ovis aries) by Solexa Deep Sequencing

Compare the microRNA (miRNA) trancriptomes of goat and sheep skin using Solexa sequencing to understand the development of skin and hair follicles (HFs). miRNA expression patterns vary in the two small RNA libraries from goat (G library) and sheep (S library) telogen skin samples. Analysis of the size distribution of 25.32 million clean reads revealed that most are 21-23 nucleotides. A total of 1910 known miRNAs and 2261 novel mature miRNAs were identified in this study. Among them, 107 novel miRNAs and 1246 known miRNAs were differentially expressed in the two libraries; 10 of the known miRNAs were identified using stem-loop quantitative real-time PCR. Furthermore, GO and KEGG pathway analysis of predicted miRNA targets illustrated the roles of these differentially expressed miRNAs in telogen HF development and growth. This study provides important information about the role of miRNAs in the regulation of HF development and their function in the telogen phase. This observation may help future investigations of the regulation of miRNAs during wool quality improvement.

Circulating Organ-Specific MicroRNAs Serve as Biomarkers in Organ-Specific Diseases: Implications for Organ Allo- and Xeno-Transplantation

Different cell types possess different miRNA expression profiles, and cell/tissue/organ-specific miRNAs (or profiles) indicate different diseases. Circulating miRNA is either actively secreted by living cells or passively released during cell death. Circulating cell/tissue/organ-specific miRNA may serve as a non-invasive biomarker for allo- or xeno-transplantation to monitor organ survival and immune rejection. In this review, we summarize the proof of concept that circulating organ-specific miRNAs serve as non-invasive biomarkers for a wide spectrum of clinical organ-specific manifestations such as liver-related disease, heart-related disease, kidney-related disease, and lung-related disease. Furthermore, we summarize how circulating organ-specific miRNAs may have advantages over conventional methods for monitoring immune rejection in organ transplantation. Finally, we discuss the implications and challenges of applying miRNA to monitor organ survival and immune rejection in allo- or xeno-transplantation.

Identification and differential expression of microRNAs in the ovaries of pigs (Sus scrofa) with high and low litter sizes

Litter size affects profitability in the swine industry. Mammalian ovaries play important roles during reproduction, including ovulation and hormone secretion, which are tightly regulated by specific microRNAs (miRNAs). In this study, we investigated the effects of specific miRNAs on porcine litter size. We compared the ovarian miRNAs of Yorkshire pigs with high (YH) and low (YL) litter sizes using Solexa sequencing technology. We identified 327 and 320 miRNAs in the ovaries of YH and YL pigs respectively. A total of 297 miRNAs were co-expressed; 30 and 23 miRNAs respectively were specifically expressed in the two libraries. A total of 83 novel miRNAs were predicted; 37 specific miRNAs were obtained, of which 21 miRNAs were upregulated and 16 miRNAs were downregulated in YH compared with YL. Additionally, 19 628 and 19 250 target genes were predicted in the two libraries respectively. The results revealed that specific miRNAs (i.e., miR-224, miR-99a, let-7c, miR-181c, miR-214 and miR-21) may affect porcine litter size. The results of this study will help in gaining understanding of the role of miRNAs in porcine litter size regulation.

Feeding of Enterococcus faecium NCIMB 10415 Leads to Intestinal miRNA-423-5p-Induced Regulation of Immune-Relevant Genes

Probiotics are widely used in human and animal health, but little is known about the mode of action of probiotics. One possible mechanism at the molecular level could be an influence on microRNAs (miRNAs) and the related immune-relevant target genes. Here, we analyzed differential expression of miRNA and potential target genes of ileal and jejunal lymphatic tissues from Enterococcus faeciumNCIMB 10415-fed piglets versus untreated controls by using next-generation sequencing. We identified miR-423-5p as being greatly affected by the treatment group (2.32-fold;P= 0.014). Validation by reverse transcription-quantitative PCR (RT-qPCR) confirmed a significant upregulation of miR-423-5p (2.11-fold;P= 0.03) and, additionally, downregulation of the important immune-relevant immunoglobulin lambda light C region (IGLC) (0.61-fold;P= 0.03) and immunoglobulin kappa constant (IGKC) (0.69-fold;P= 0.04) target genes. Expression analysis of miR-423-5p and IGLC at different age points shows a clear anti correlated relationship. Luciferase reporter assays with a HeLa cell line verified IGLC as a target of miR-423-5p. The results provided evidence for an effect of feeding of E. faeciumon the expression of miR-423-5p and on the regulation of the IGLC gene through miR-423-5p. This might be a possible mode of action of E. faeciumon immune cell regulation in the small intestine.

Identification of Novel and Conserved miRNAs in Leaves of In vitro Grown Citrus reticulata "Lugan" Plantlets by Solexa Sequencing

MicroRNAs (miRNAs) play essential roles in plant development, but the roles in the in vitro plant development are unknown. Leaves of ponkan plantlets derived from mature embryos at in vitro culture conditions were used to sequence small RNA fraction via Solexa sequencing, and the miRNAs expression was analyzed. The results showed that there were 3,065,625 unique sequences in ponkan, of which 0.79% were miRNAs. The RNA sequences with lengths of 18-25 nt derived from the library were analyzed, leading to the identification of 224 known miRNAs, of which the most abundant were miR157, miR156, and miR166. Three hundred and fifty-eight novel miRNA candidates were also identified, and the number of reads of ponkan novel miRNAs varied from 5 to 168,273. The expression of the most known miRNAs obtained was at low levels, which varied from 5 to 4,946,356. To better understand the role of miRNAs during the preservation of ponkan in vitro plantlet, the expression patterns of cre-miR156a/159b/160a/166a/167a/168a/171/398b were validated by quantitative real-time PCR (qPCR). The results showed that not only the development-associated miRNAs, e.g., cre-miR156/159/166/396, expressed highly at the early preservation period in the in vitro ponkan plantlet leaves but also the stress-related miRNAs, e.g., cre-miR171 and cre-miR398b, expressed highly at the same time. The expression levels of most tested miRNAs were found to decrease after 6 months and the amounts of these miRNAs were kept at low levels at 18 months. After analyzing the expression level of their targets during the reservation of the ponkan in vitro plantlet, development-associated cre-ARF6 and stress-related cre-CSD modules exhibited negative correlation with miR167 and miR398, respectively, indicating an involvement of the miRNAs in the in vitro development of ponkan and function in the conservation of ponkan germplasm.

Deciphering the microRNA transcriptome of skeletal muscle during porcine development

MicroRNAs (miRNAs) play critical roles in many important biological processes, such as growth and development in mammals. Various studies of porcine muscle development have mainly focused on identifying miRNAs that are important for fetal and adult muscle development; however, little is known about the role of miRNAs in middle-aged muscle development. Here, we present a comprehensive investigation of miRNA transcriptomes across five porcine muscle development stages, including one prenatal and four postnatal stages. We identified 404 known porcine miRNAs, 118 novel miRNAs, and 101 miRNAs that are conserved in other mammals. A set of universally abundant miRNAs was found across the distinct muscle development stages. This set of miRNAs may play important housekeeping roles that are involved in myogenesis. A short time-series expression miner analysis indicated significant variations in miRNA expression across distinct muscle development stages. We also found enhanced differentiation- and morphogenesis-related miRNA levels in the embryonic stage; conversely, apoptosis-related miRNA levels increased relatively later in muscle development. These results provide integral insight into miRNA function throughout pig muscle development stages. Our findings will promote further development of the pig as a model organism for human age-related muscle disease research.

Genome-wide analysis of microRNAs identifies the lipid metabolism pathway to be a defining factor in adipose tissue from different sheep

MicroRNAs are short (17–24 nt) non-coding RNAs that are involved in post-transcriptional regulation of gene expression in multicellular organisms by affecting both the stability and translation of mRNAs. In recent years, deep sequencing of the transcriptome is increasingly being utilized with the promise of higher sensitivity for the identification of differential expression patterns as well as the opportunity to discover new transcripts, including new alternative isoforms and miRNAs. Here, we utilized RNA-seq technology to perform a genome-wide analysis of miRNAs from the adipose tissue of the two species of sheep to look for clues that might explain the fat deposition differences between the sheep. The RNA-seq analysis detected 3132 miRNAs from the adipose tissue of the Small-tail Han and Dorset sheep, of which 2893 were defined as potential new miRNAs. In addition, 54 miRNAs were differentially expressed between the two breeds of sheep. Gene ontology and pathway analyses of the predicted target genes that negatively associated with the differentially expressed miRNAs revealed that there was less active lipid metabolism in the adipose tissue of Small Tail Han sheep. This study can help understand the underling mechanisms responsible for the morphological differences related to fat deposition between two breeds of sheep.

Expression Patterns of miRNA-423-5p in the Serum and Pericardial Fluid in Patients Undergoing Cardiac Surgery

BACKGROUND

Recently, it has been reported that specific microRNA (miRNA) levels are elevated in serum and can be used as biomarkers in patients with cardiovascular diseases. However, miRNAs expression profiles and their sources in pericardial fluid (PF) are unclear.

METHODS AND RESULTS

The purpose of this study was to identify the levels of miRNAs in PF in relation to those in the serum in patients undergoing cardiac surgery. Serum (S) and PF from patients undergoing coronary artery bypass graft (CABG) due to stable angina pectoris (sAP) and unstable AP (uAP) and aortic valve replacement due to aortic stenosis (AS) were analyzed for the detection of miRNAs. We named these samples S-sAP, S-uAP, S-AS, PF-sAP, PF-uAP, and PF-AS, respectively. We first measured the levels of miR-423-5p, which was recognized previously as a biomarker for heart failure. miR-423-5p levels were significantly higher in PF than serum. Although there was no difference in miR-423-5p levels among the PF-AS, PF-sAP, and PF-uAP, its levels were significantly elevated in S-uAP compared with those in S-AS and S-sAP. In order to clarify the source of miR-423-5p in PF, we measured the levels of muscle-enriched miR-133a and vascular-enriched miR-126 and miR-92a in the same samples. miR-133a levels were significantly higher in serum than in PF, and it was elevated in S-uAP compared with S-AS. miR-126 level was significantly increased in serum compared with PF, and the level of miR-92a the similar tendency. miR-423-5p is located in the first intron of NSRP1. There is another miRNA, miR-3184, encoded in the opposite direction in the same region. In vitro experiments indicated that the duplex of miR-423-5p and miR-3184-3p was more resistant to RNase than the duplex of miR-423-5p and miR-133-3p, which may help to stabilize miR-423-5p in the PF.

CONCLUSIONS

Our results suggested that miR-423-5p is enriched in PF, and serum miR-423-5p may be associate with uAP. Its expression pattern was different to that of muscle- and vascular-enriched miRNAs, miR-133a, miR-126, and miR-92a.

Identification and profiling of microRNAs in goat endometrium during embryo implantation

Background

MicroRNAs (miRNAs) are short, highly conserved small noncoding RNAs that had fundamental roles in post-transcriptional gene expression, and they are crucial for proper control of biological processes and known to participate in embryo implantation. However, miRNA expression profiles in the pre-receptive and receptive phases of the goat endometrium during embryo implantation are unknown.

Results

A total of 1,069 and 847 miRNAs were expressed in receptive (R) and pre-receptive (P) goat endometrium, and 632 miRNAs were co-expressed in both phases. We identified 545 (50.98%) known miRNAs in the R library and 522 (61.63%) in the P library. There were 110 up-expressed miRNAs and 33 down-expressed miRNAs in receptive endometrium compared with the pre-receptive endometrium meeting the criteria of P-values< 0.05. Moreover, GO and KEGG analysis of the target genes of the differentially expressed miRNAs revealed some candidate miRNAs, genes and pathways that may involve in the formation of the receptive endometrium. Based on stem-loop RT-qPCR, 15 miRNAs were detected and the results suggested that the majority of the miRNA expression data measured by Solexa deep sequencing could represent actual miRNA expression levels.

Conclusions

Our data revealed the first miRNA profile related to the biology of the goat receptive endometrium during embryo implantation, and the results suggested that a subset of miRNAs might play important roles in the formation of endometrial receptivity. Thus, elucidating the physiological roles of endometrial miRNAs will help us better understand the genetic control of embryo implantation in goats.

Identification of differentially expressed microRNAs in placentas of cloned and normally produced calves by Solexa sequencing

Aberrant gene expression caused by aberrant nuclear reprogramming results in developmental abnormalities in cloned animals and ultimately their death. MicroRNAs (miRNAs), a family of ∼ 22 nucleotide, non-coding, single-stranded RNA molecules, are considered as key regulators of gene expression. Numerous miRNAs and their expression patterns have been identified in various species. However, the significance of miRNAs in developmental abnormalities in cloned animals is unclear. Small RNA libraries were generated from the placentas of cloned (somatic cell nuclear transfer, SCNT) and normally produced (control) calves. A total of 18,815,541 clean reads were obtained from the SCNT library and 19,329,352 from the control library. In total, 430 conserved bovine miRNAs were identified in bovine placenta. Furthermore, the family, expression predominance, and base substitution of the conserved miRNAs were also analyzed. We found 135 conserved miRNAs that were differentially expressed significantly between the two samples, which suggest that these miRNAs may affect developmental abnormalities in cloned cattle and ultimately their death. The miRNA target prediction, gene ontology, and pathway analysis for these target genes were also carried out. The present study expands the collection of bovine miRNAs and could initiate further studies on the functions of miRNAs in developmental abnormalities and death in cloned animals.

Interactome mapping reveals important pathways in skeletal muscle development of pigs

The regulatory relationship and connectivity among genes involved in myogenesis and hypertrophy of skeletal muscle in pigs still remain large challenges. Presentation of gene interactions is a potential way to understand the mechanisms of developmental events in skeletal muscle. In this study, genome-wide transcripts and miRNA profiling was determined for Landrace pigs at four time points using microarray chips. A comprehensive method integrating gene ontology annotation and interactome network mapping was conducted to analyze the biological patterns and interaction modules of muscle development events based on differentially expressed genes and miRNAs. Our results showed that in total 484 genes and 34 miRNAs were detected for the duration from embryonic stage to adult in pigs, which composed two linear expression patterns with consensus changes. Moreover, the gene ontology analysis also disclosed that there were three typical biological events i.e., microstructure assembly of sarcomere at early embryonic stage, myofibril formation at later embryonic stage and function establishments of myoblast cells at postnatal stage. The interactome mappings of different time points also found the down-regulated trend of gene expression existed across the whole duration, which brought a possibility to introduce the myogenesis related miRNAs into the interactome regulatory networks of skeletal muscle in pigs.

Current transcriptomics in pig immunity research

Swine performance in the face of disease challenge is becoming progressively more important. To improve the pig’s robustness and resilience against pathogens through selection, a better understanding of the genetic and epigenetic factors in the immune response is required. This review highlights results from the most recent transcriptome research, and the meta-analyses performed, in the context of pig immunity. A technological overview is given including wholegenome microarrays, immune-specific arrays, small-scale high-throughput expression methods, high-density tiling arrays, and next generation sequencing (NGS). Although whole genome microarray techniques will remain complementary to NGS for some time in domestic species, research will transition to sequencing-based methods due to cost-effectiveness and the extra information that such methods provide. Furthermore, upcoming high-throughput epigenomic studies, which will add greatly to our knowledge concerning the impact of epigenetic modifications on pig immune response, are listed in this review. With emphasis on the insights obtained from transcriptomic analyses for porcine immunity, we also discuss the experimental design in pig immunity research and the value of the newly published porcine genome assembly in using the pig as a model for human immune response. We conclude by discussing the importance of establishing community standards to maximize the possibility of integrative computational analyses, such as was clearly beneficial for the human ENCODE project.

Growth hormone-regulated mRNAs and miRNAs in chicken hepatocytes

Growth hormone (GH) is a key regulatory factor in animal growth, development and metabolism. Based on the expression level of the GH receptor, the chicken liver is a major target organ of GH, but the biological effects of GH on the chicken liver are not fully understood. In this work we identified mRNAs and miRNAs that are regulated by GH in primary hepatocytes from female chickens through RNA-seq, and analyzed the functional relevance of these mRNAs and miRNAs through GO enrichment analysis and miRNA target prediction. A total of 164 mRNAs were found to be differentially expressed between GH-treated and control chicken hepatocytes, of which 112 were up-regulated and 52 were down-regulated by GH. A total of 225 chicken miRNAs were identified by the RNA-Seq analysis. Among these miRNAs 16 were up-regulated and 1 miRNA was down-regulated by GH. The GH-regulated mRNAs were mainly involved in growth and metabolism. Most of the GH-upregulated or GH-downregulated miRNAs were predicted to target the GH-downregulated or GH-upregulated mRNAs, respectively, involved in lipid metabolism. This study reveals that GH regulates the expression of many mRNAs involved in metabolism in female chicken hepatocytes, which suggests that GH plays an important role in regulating liver metabolism in female chickens. The results of this study also support the hypothesis that GH regulates lipid metabolism in chicken liver in part by regulating the expression of miRNAs that target the mRNAs involved in lipid metabolism.

Expression profiling of exosomal miRNAs derived from human esophageal cancer cells by Solexa high-throughput sequencing

Cellular genetic materials, such as microRNAs (miRNAs), mRNAs and proteins, are packaged inside exosomes, small membrane vesicles of endocytic origin that are released into the extracellular environment. These cellular genetic materials can be delivered into recipient cells, where they exert their respective biological effects. However, the miRNA profiles and biological functions of exosomes secreted by cancer cells remain unknown. The present study explored the miRNA expression profile and distribution characteristics of exosomes derived from human esophageal cancer cells through Solexa high-throughput sequencing. Results showed that 56,421 (2.94%) unique sequences in cells and 7727 (0.63%) in exosomes matched known miRNAs. A total of 342 and 48 known miRNAs were identified in cells and exosomes, respectively. Moreover, 64 and 32 novel miRNAs were predicted in cells and exosomes, respectively. Significant differences in miRNA expression profiles were found between human esophageal cancer cells and exosomes. These findings provided new insights into the characteristics of miRNAs in exosomes derived from human esophageal cancer cells and the specific roles of miRNAs in intercellular communication mediated by exosomes in esophageal cancer.

Identification of conserved and novel microRNAs in cerebral ischemia-reperfusion injury of rat using deep sequencing

MicroRNAs are a class of noncoding small RNAs that regulate gene expression by inhibiting target genes at post-transcriptional levels. MicroRNAs have been highlighted in many organs and tissues, including the brain. To identify special microRNAs involved in ischemia-reperfusion injury, we performed a comprehensive small RNA profiling in rat model and the control using Illumina high-throughput sequencing. A total of 9,444,562 and 10,290,391 clean reads were sequenced from two small RNA libraries constructed, respectively. Three hundred fifty-eight known microRNAs were identified, in which 78 microRNAs exhibited significantly differential expression between model and control. In addition, 62 and 68 novel miRNAs were found in model and control, respectively. Comparative analysis showed that 24 novel microRNAs were differentially expressed with greater than six-fold change. The GO annotation suggested that predicted targets of microRNAs were enriched into the category of metabolic process, cell part, cell-extracellular communications, and so on. KEGG pathway analysis suggested that these genes were involved in many important pathways, mainly including signaling transduction, MAPK signaling pathway, NF-κB signaling pathway, and neurotrophin signaling pathway. Our findings provided a deeper understanding to the regulatory mechanism of microRNAs underlying cerebral ischemia, therefore benefitting the improvement of the protection and treatment strategies of this disease.

MicroRNA expression profiling of the fifth-instar posterior silk gland of Bombyx mori

BACKGROUND

The growth and development of the posterior silk gland and the biosynthesis of the silk core protein at the fifth larval instar stage of Bombyx mori are of paramount importance for silk production.

RESULTS

Here, aided by next-generation sequencing and microarry assay, we profile 1,229 microRNAs (miRNAs), including 728 novel miRNAs and 110 miRNA/miRNA* duplexes, of the posterior silk gland at the fifth larval instar. Target gene prediction yields 14,222 unique target genes from 1,195 miRNAs. Functional categorization classifies the targets into complex pathways that include both cellular and metabolic processes, especially protein synthesis and processing.

CONCLUSION

The enrichment of target genes in the ribosome-related pathway indicates that miRNAs may directly regulate translation. Our findings pave a way for further functional elucidation of these miRNAs and their targets in silk production.

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.

Tissue-specific expression and regulatory networks of pig microRNAome

Background

Despite the economic and medical importance of the pig, knowledge about its genome organization, gene expression regulation, and molecular mechanisms involved in physiological processes is far from that achieved for mouse and rat, the two most used model organisms in biomedical research. MicroRNAs (miRNAs) are a wide class of molecules that exert a recognized role in gene expression modulation, but only 280 miRNAs in pig have been characterized to date.

Results

We applied a novel computational approach to predict species-specific and conserved miRNAs in the pig genome, which were then subjected to experimental validation. We experimentally identified candidate miRNAs sequences grouped in high-confidence (424) and medium-confidence (353) miRNAs according to RNA-seq results. A group of miRNAs was also validated by PCR experiments. We established the subtle variability in expression of isomiRs and miRNA-miRNA star couples supporting a biological function for these molecules. Finally, miRNA and mRNA expression profiles produced from the same sample of 20 different tissue of the animal were combined, using a correlation threshold to filter miRNA-target predictions, to identify tissue-specific regulatory networks.

Conclusions

Our data represent a significant progress in the current understanding of miRNAome in pig. The identification of miRNAs, their target mRNAs, and the construction of regulatory circuits will provide new insights into the complex biological networks in several tissues of this important animal model.

A feedback circuit between miR-133 and the ERK1/2 pathway involving an exquisite mechanism for regulating myoblast proliferation and differentiation

MiR-133 was found to be specifically expressed in cardiac and skeletal muscle in previous studies. There are two members in the miR-133 family: miR-133a and miR-133b. Although previous studies indicated that miR-133a was related to myogenesis, the signaling pathways regulated by miR-133 were still not very clear. In this study, we showed that both miR-133a and miR-133b were upregulated during myogenesis through Solexa sequencing. We confirmed that miR-133 could promote myoblast differentiation and inhibit cell proliferation through the regulation of the extracellular signal-regulated kinase (ERK) signaling pathway in C2C12 cells. FGFR1 and PP2AC, which both participate in signal transduction of the ERK1/2 pathway, were found to be negatively regulated by miR-133a and miR-133b at the post-transcriptional level. Also, downregulation of ERK1/2 phosphorylation by miR-133 was detected. FGFR1 and PP2AC were also found to repress C2C12 differentiation by specific siRNAs. In addition, we found that inhibition of ERK1/2 pathway activity can inhibit C2C12 cell proliferation and promote the initiation of differentiation but form short and small myotubes. Furthermore, we found that the expression of miR-133 was negatively regulated by ERK1/2 signaling pathway. In summary, we demonstrated the role of miR-133 in myoblast and further revealed a new feedback loop between miR-133 and the ERK1/2 signaling pathway involving an exquisite mechanism for regulating myogenesis.

MicroRNA expression profiling of the porcine developing hypothalamus and pituitary tissue

MicroRNAs (miRNAs), a class of small non-coding RNA molecules, play important roles in gene expressions at transcriptional and post-transcriptional stages in mammalian brain. So far, a growing number of porcine miRNAs and their function have been identified, but little is known regarding the porcine developing hypothalamus and pituitary. In the present study, Solexa sequencing analysis showed 14,129,397 yielded reads, 6,680,678 of which were related to 674 unique miRNAs. After a microarray assay, we detected 175 unique miRNAs in the hypothalamus, including 136 previously known miRNAs and 39 novel candidates, while a total of 140 miRNAs, including 104 known and 36 new candidate miRNAs, were discovered in pituitary. More importantly, 37 and 30 differentially expressed miRNAs from several developmental stages of hypothalamus and pituitary were revealed, respectively. The 37 differentially expressed miRNAs in hypothalamus represented 6 different expression patterns, while the 30 differentially expressed miRNAs in pituitary represented 7 different expression patterns. To clarify potential target genes and specific functions of these differentially expressed miRNAs in hypothalamus and pituitary, TargetScan and Gorilla prediction tools were then applied. The current functional analysis showed that the differentially expressed miRNAs in hypothalamus and pituitary shared many biological processes, with the main differences being found in tissue-specific processes including: CDP-diacylglycerol biosynthetic/metabolic process; phosphatidic acid biosynthetic/metabolic process; energy reserve metabolic process for hypothalamus; adult behavior; sterol transport/homeostasis; and cholesterol/reverse cholesterol transport for pituitary. Overall, this study identified miRNA profiles and differentially expressed miRNAs among various developmental stages in hypothalamus and pituitary and indicated miRNA profiles change with age and brain location, enhancing our knowledge about spatial and temporal expressions of miRNAs in the porcine developing brain.

Selected isomiR expression profiles via arm switching?

A mature miRNA may be generated from 5p or 3p arm of a hairpin precursor. The selection may be flexible via "arm switching". However, accumulating evidences suggest that both arms of many pre-miRNAs can yield mature functional miRNAs. Herein, we attempted to compare the isomiR expression profiles between the two arms through analyzing in-house and published small RNA deep sequencing datasets. Although many miR-#-5p and miR-#-3p have been reported as functional miRNAs, fewer miRNA pairs (11 and 6 pairs are collected in tumor and normal cells, respectively) are simultaneously identified as abundant miRNA species. According to isomiR types and dominant isomiR species, miR-#-5p and miR-#-3p show various isomiR expression profiles as well as diverse enrichment levels. IsomiR profiles of non-dominant arm are not well-conserved in 5' ends as well as isomiR profiles of dominant arm. If both the miR-#-5p and miR-#-3p are abundantly expressed, their isomiR expression profiles are always stable across different samples. Similar to diverse enrichment levels of miR-#-5p and miR-#-3p, the isomiR expression patterns may also be influenced by the phenomenon of "arm switching". The diverged isomiR expression profiles further enrich the complexity of multiple isomiRs, and complicate the coding-non-coding RNA regulatory network.

PHDcleav: a SVM based method for predicting human Dicer cleavage sites using sequence and secondary structure of miRNA precursors

BACKGROUND

Dicer, an RNase III enzyme, plays a vital role in the processing of pre-miRNAs for generating the miRNAs. The structural and sequence features on pre-miRNA which can facilitate position and efficiency of cleavage are not well known. A precise cleavage by Dicer is crucial because an inaccurate processing can produce miRNA with different seed regions which can alter the repertoire of target genes.

RESULTS

In this study, a novel method has been developed to predict Dicer cleavage sites on pre-miRNAs using Support Vector Machine. We used the dataset of experimentally validated human miRNA hairpins from miRBase, and extracted fourteen nucleotides around Dicer cleavage sites. We developed number of models using various types of features and achieved maximum accuracy of 66% using binary profile of nucleotide sequence taken from 5p arm of hairpin. The prediction performance of Dicer cleavage site improved significantly from 66% to 86% when we integrated secondary structure information. This indicates that secondary structure plays an important role in the selection of cleavage site. All models were trained and tested on 555 experimentally validated cleavage sites and evaluated using 5-fold cross validation technique. In addition, the performance was also evaluated on an independent testing dataset that achieved an accuracy of ~82%.

CONCLUSION

Based on this study, we developed a webserver PHDcleav (http://www.imtech.res.in/raghava/phdcleav/) to predict Dicer cleavage sites in pre-miRNA. This tool can be used to investigate functional consequences of genetic variations/SNPs in miRNA on Dicer cleavage site, and gene silencing. Moreover, it would also be useful in the discovery of miRNAs in human genome and design of Dicer specific pre-miRNAs for potent gene silencing.

Characterization of microRNAs in mud crab Scylla paramamosain under Vibrio parahaemolyticus infection

BACKGROUND

Infection of bacterial Vibrio parahaemolyticus is common in mud crab farms. However, the mechanisms of the crab's response to pathogenic V. parahaemolyticus infection are not fully understood. MicroRNAs (miRNAs) are a class of small noncoding RNAs that function as regulators of gene expression and play essential roles in various biological processes. To understand the underlying mechanisms of the molecular immune response of the crab to the pathogens, high-throughput Illumina/Solexa deep sequencing technology was used to investigate the expression profiles of miRNAs in S. paramamosain under V. parahaemolyticus infection.

METHODOLOGY/PRINCIPAL FINDINGS

Two mixed RNA pools of 7 tissues (intestine, heart, liver, gill, brain, muscle and blood) were obtained from V. parahaemolyticus infected crabs and the control groups, respectively. By aligning the sequencing data with known miRNAs, we characterized 421 miRNA families, and 133 conserved miRNA families in mud crab S. paramamosain were either identical or very similar to existing miRNAs in miRBase. Stem-loop qRT-PCRs were used to scan the expression levels of four randomly chosen differentially expressed miRNAs and tissue distribution. Eight novel potential miRNAs were confirmed by qRT-PCR analysis and the precursors of these novel miRNAs were verified by PCR amplification, cloning and sequencing in S. paramamosain. 161 miRNAs (106 of which up-regulated and 55 down-regulated) were significantly differentially expressed during the challenge and the potential targets of these differentially expressed miRNAs were predicted. Furthermore, we demonstrated evolutionary conservation of mud crab miRNAs in the animal evolution process.

CONCLUSIONS/SIGNIFICANCE

In this study, a large number of miRNAs were identified in S. paramamosain when challenged with V. parahaemolyticus, some of which were differentially expressed. The results show that miRNAs might play some important roles in regulating gene expression in mud crab under V. parahaemolyticus infection, providing a basis for further investigation of miRNA-modulating networks in innate immunity of mud crab.

Close association between paralogous multiple isomiRs and paralogous/orthologues miRNA sequences implicates dominant sequence selection across various animal species

MicroRNAs (miRNAs) are crucial negative regulators of gene expression at the post-transcriptional level. Next-generation sequencing technologies have identified a series of miRNA variants (named isomiRs). In this study, paralogous isomiR assemblies (from the miRNA locus) were systematically analyzed based on data acquired from deep sequencing data sets. Evolutionary analysis of paralogous (members in miRNA gene family in a specific species) and orthologues (across different animal species) miRNAs was also performed. The sequence diversity of paralogous isomiRs was found to be similar to the diversity of paralogous and orthologues miRNAs. Additionally, both isomiRs and paralogous/orthologues miRNAs were implicated in 5' and 3' ends (especially 3' ends), nucleotide substitutions, and insertions and deletions. Generally, multiple isomiRs can be produced from a single miRNA locus, but most of them had lower enrichment levels, and only several dominant isomiR sequences were detected. These dominant isomiR groups were always stable, and one of them would be selected as the most abundant miRNA sequence in specific animal species. Some isomiRs might be consistent to miRNA sequences in some species but not the other. Homologous miRNAs were often detected in similar isomiR repertoires, and showed similar expression patterns, while dominant isomiRs showed complex evolutionary patterns from miRNA sequences across the animal kingdom. These results indicate that the phenomenon of multiple isomiRs is not a random event, but rather the result of evolutionary pressures. The existence of multiple isomiRs enables different species to express advantageous sequences in different environments. Thus, dominant sequences emerge in response to functional and evolutionary pressures, allowing an organism to adapt to complex intra- and extra-cellular events.

Phytoplasma-Responsive microRNAs Modulate Hormonal, Nutritional, and Stress Signalling Pathways in Mexican Lime Trees

BACKGROUND

Witches' broom disease of Mexican lime (Citrus aurantifolia L.), which is associated to the phytoplasma 'Candidatus Phytoplasma aurantifolia', is a devastating disease that results in significant economic losses. Plants adapt to biotic stresses by regulating gene expression at the transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) are a recently identified family of molecules that regulate plant responses to environmental stresses through post-transcriptional gene silencing.

METHODS

Using a high-throughput approach to sequence small RNAs, we compared the expression profiles of miRNAs in healthy Mexican lime trees and in plants infected with 'Ca. P. aurantifolia'.

RESULTS

Our results demonstrated the involvement of different miRNAs in the response of Mexican lime trees to infection by 'Ca. P. aurantifolia'. We identified miRNA families that are expressed differentially upon infection with phytoplasmas. Most of the miRNAs had variants with small sequence variations (isomiRs), which are expressed differentially in response to pathogen infection.

CONCLUSIONS

It is likely that the miRNAs that are expressed differentially in healthy and phytoplasma-infected Mexican lime trees are involved in coordinating the regulation of hormonal, nutritional, and stress signalling pathways, and the complex interactions between them. Future research to elucidate the roles of these miRNAs should improve our understanding of the level of diversity of specific plant responses to phytoplasmas.

A Review of Computational Tools in microRNA Discovery

Since microRNAs (miRNAs) were discovered, their impact on regulating various biological activities has been a surprising and exciting field. Knowing the entire repertoire of these small molecules is the first step to gain a better understanding of their function. High throughput discovery tools such as next-generation sequencing significantly increased the number of known miRNAs in different organisms in recent years. However, the process of being able to accurately identify miRNAs is still a complex and difficult task, requiring the integration of experimental approaches with computational methods. A number of prediction algorithms based on characteristics of miRNA molecules have been developed to identify new miRNA species. Different approaches have certain strengths and weaknesses and in this review, we aim to summarize several commonly used tools in metazoan miRNA discovery.

MicroRNA transcriptomes relate intermuscular adipose tissue to metabolic risk

Intermuscular adipose tissue is located between the muscle fiber bundles in skeletal muscles, and has similar metabolic features to visceral adipose tissue, which has been found to be related to a number of obesity-related diseases. Although various miRNAs are known to play crucial roles in adipose deposition and adipogenesis, the microRNA transcriptome of intermuscular adipose tissue has not, until now, been studied. Here, we sequenced the miRNA transcriptomes of porcine intermuscular adipose tissue by small RNA-sequencing and compared it to a representative subcutaneous adipose tissue. We found that the inflammation- and diabetes-related miRNAs were significantly enriched in the intermuscular rather than in the subcutaneous adipose tissue. A functional enrichment analysis of the genes predicted to be targeted by the enriched miRNAs also indicated that intermuscular adipose tissue was associated mainly with immune and inflammation responses. Our results suggest that the intermuscular adipose tissue should be recognized as a potential metabolic risk factor of obesity.

Characterization and differential expression of microRNAs in the ovaries of pregnant and non-pregnant goats (Capra hircus)

Background

Ovarian follicular development and hormone secretion are complex and coordinated biological processes which will usually be altered during pregnancy. Ovarian function is tightly regulated by a multitude of genes, and also by some specific miRNAs. It is necessary to identify the differentially expressed miRNAs in the ovaries of pregnant and non-pregnant mammals, in order to further understand the role of miRNA-mediated post-transcriptional regulation in mammalian reproduction. Here, we performed a comprehensive search for hircine miRNAs using two small RNA sequencing libraries prepared from the ovaries of pregnant and non-pregnant goats.

Results

617 conserved and 7 putative novel miRNAs were identified in the hircine ovaries. A total of 471 conserved miRNAs (76.34%) were co-expressed in both pregnant and non-pregnant libraries, and 90 pregnancy-specific and 56 non-pregnancy-specific conserved miRNAs were identified. Additionally, 407 unique miRNAs (65.96%) were significantly differentially expressed in the pregnant and non-pregnant libraries, of which 294 were upregulated and 113 were downregulated in the pregnant library compared to the non-pregnant library. Further analysis showed that miR-143 was predicted to bind to the target sequences of Frizzled-6 and -3 receptor genes in the Wnt/beta-catenin signaling pathway, and let-7b may target the Activin receptor I and Smad 2/3 genes in the TGF-beta signaling pathway. The expression level of 5 randomly selected miRNAs were analyzed by quantitative real-time PCR (q-PCR), and the results demonstrated that the expression patterns were consistent with the Solexa sequencing results.

Conclusions

The identification and characterization of differentially expressed miRNAs in the ovaries of pregnant and non-pregnant goats provides important information on the role of miRNA in the regulation of the ovarian development and function. This data will be helpful to facilitate studies on the regulation of miRNAs during mammalian reproduction.

Identification of differences in microRNA transcriptomes between porcine oxidative and glycolytic skeletal muscles

Background

MicroRNAs (miRNAs) are a type of non-coding small RNA ~22 nucleotides in length that regulate the expression of protein coding genes at the post-transcriptional level. Glycolytic and oxidative myofibers, the two main types of skeletal muscles, play important roles in metabolic health as well as in meat quality and production in the pig industry. Previous expression profile studies of different skeletal muscle types have focused on these aspects of mRNA and proteins; nonetheless, an explanation of the miRNA transcriptome differences between these two distinct muscles types is long overdue.

Results

Herein, we present a comprehensive analysis of miRNA expression profiling between the porcine longissimus doris muscle (LDM) and psoas major muscle (PMM) using a deep sequencing approach. We generated a total of 16.62 M (LDM) and 18.46 M (PMM) counts, which produced 15.22 M and 17.52 M mappable sequences, respectively, and identified 114 conserved miRNAs and 89 novel miRNA*s. Of 668 unique miRNAs, 349 (52.25%) were co-expressed, of which 173 showed significant differences (P < 0.01) between the two muscle types. Muscle-specific miR-1-3p showed high expression levels in both libraries (LDM, 32.01%; PMM, 20.15%), and miRNAs that potentially affect metabolic pathways (such as the miR-133 and -23) showed significant differences between the two libraries, indicating that the two skeletal muscle types shared mainly muscle-specific miRNAs but expressed at distinct levels according to their metabolic needs. In addition, an analysis of the Gene Ontology (GO) terms and KEGG pathway associated with the predicted target genes of the differentially expressed miRNAs revealed that the target protein coding genes of highly expressed miRNAs are mainly involved in skeletal muscle structural development, regeneration, cell cycle progression, and the regulation of cell motility.

Conclusion

Our study indicates that miRNAs play essential roles in the phenotypic variations observed in different muscle fiber types.

A highly expressed miR-101 isomiR is a functional silencing small RNA

Background

MicroRNAs (miRNAs) are short non-coding regulatory RNAs that control gene expression usually producing translational repression and gene silencing. High-throughput sequencing technologies have revealed heterogeneity at length and sequence level for the majority of mature miRNAs (IsomiRs). Most isomiRs can be explained by variability in either Dicer1 or Drosha cleavage during miRNA biogenesis at 5’ or 3’ of the miRNA (trimming variants). Although isomiRs have been described in different tissues and organisms, their functional validation as modulators of gene expression remains elusive. Here we have characterized the expression and function of a highly abundant miR-101 5’-trimming variant (5’-isomiR-101).

Results

The analysis of small RNA sequencing data in several human tissues and cell lines indicates that 5’-isomiR-101 is ubiquitously detected and a highly abundant, especially in the brain. 5’-isomiR-101 was found in Ago-2 immunocomplexes and complementary approaches showed that 5’-isomiR-101 interacted with different members of the silencing (RISC) complex. In addition, 5’-isomiR-101 decreased the expression of five validated miR-101 targets, suggesting that it is a functional variant. Both the binding to RISC members and the degree of silencing were less efficient for 5’-isomiR-101 compared with miR-101. For some targets, both miR-101 and 5’-isomiR-101 significantly decreased protein expression with no changes in the respective mRNA levels. Although a high number of overlapping predicted targets suggest similar targeted biological pathways, a correlation analysis of the expression profiles of miR-101 variants and predicted mRNA targets in human brains at different ages, suggest specific functions for miR-101- and 5’-isomiR-101.

Conclusions

These results suggest that isomiRs are functional variants and further indicate that for a given miRNA, the different isomiRs may contribute to the overall effect as quantitative and qualitative fine-tuners of gene expression.

miRNA expression profile analysis in kidney of different porcine breeds

microRNAs (miRNAs) are important post-transcriptional regulators in eukaryotes that target mRNAs repressing their expression. The uncertain process of pig domestication, with different origin focuses, and the selection process that commercial breeds suffered, have generated a wide spectrum of breeds with clear genetic and phenotypic variability. The aim of this work was to define the miRNAs expression profile in kidney of several porcine breeds. Small RNA libraries from kidney were elaborated and high-throughput sequenced with the 454 Genome Sequencer FLX (Roche). Pigs used were classified into three groups: the European origin group (Iberian breed and European Wild Boar ancestor), European commercial breeds (Landrace, Large White and Piétrain breeds) and breeds with Asian origin (Meishan and Vietnamese breeds). A total of 229 miRNAs were described in the pig kidney miRNA profile, including 110 miRNAs out of the 257 previously described pig miRNAs and 119 orthologous miRNAs. The most expressed miRNAs in pig kidney microRNAome were Hsa-miR-200b-3p, Ssc-miR-125b and Ssc-miR-23b. Moreover, 5 novel porcine miRNAs and 3 orthologous miRNAs could be validated through RT-qPCR. miRNA sequence variation was determined in 116 miRNAs, evidencing the presence of isomiRs. 125 miRNAs were differentially expressed between breed groups. The identification of breed-specific miRNAs, which could be potentially associated to certain phenotypes, is becoming a new tool for the study of the genetic variability underlying complex traits and furthermore, it adds a new layer of complexity to the interesting process of pig evolution.

Isolation of total RNA and detection procedures for miRNA present in bovine-cultured adipocytes and adipose tissues

Micro-ribonucleic acids (miRNA) regulate gene expression posttranscriptionally by altering translation of protein(s) encoded by specific messenger RNA. Therefore the ability to detect and quantify the expression levels of specific miRNA present within a cell or tissue is necessary to thoroughly examine cellular physiology and gene expression. Here we describe procedures that allow for the isolation and quantification of miRNA in bovine adipocytes and adipose tissue.

Comparative transcriptome profiling of dairy goat microRNAs from dry period and peak lactation mammary gland tissues

BACKGROUND

MicroRNAs (miRNAs) are small noncoding RNA molecules that serve as important post-transcriptional gene expression regulators by targeting messenger RNAs for post-transcriptional endonucleolytic cleavage or translational inhibition. miRNAs play important roles in many biological processes. Extensive high-throughput sequencing studies of miRNAs have been performed in several animal models. However, little is known about the diversity of these regulatory RNAs in goat (Capra hircus), which is one of the most important agricultural animals and the oldest domesticated species raised worldwide. Goats have long been used for their milk, meat, hair (including cashmere), and skins throughout much of the world.

RESULTS

In this study, two small RNA libraries were constructed based on dry period and peak lactation dairy goat mammary gland tissues and sequenced using the Illumina-Solexa high-throughput sequencing technology. A total of 346 conserved and 95 novel miRNAs were identified in the dairy goat. miRNAs expression was confirmed by qRT-PCR in nine tissues and in the mammary gland during different stages of lactation. In addition, several candidate miRNAs that may be involved in mammary gland development and lactation were found by comparing the miRNA expression profiles in different tissues and developmental stages of the mammary gland.

CONCLUSIONS

This study reveals the first miRNAs profile related to the biology of the mammary gland in the dairy goat. The characterization of these miRNAs could contribute to a better understanding of the molecular mechanisms of lactation physiology and mammary gland development in the dairy goat.

Identification of novel and differentially expressed MicroRNAs of dairy goat mammary gland tissues using solexa sequencing and bioinformatics

MicroRNAs are small, noncoding RNA molecules that regulate gene expression at the post-transcriptional level and play an important role in various biological processes. Although most microRNAs expression profiles studies have been performed in humans or rodents, relatively limited knowledge also exists in other mammalian species. The identification of the full repertoire of microRNAs expressed in the lactating mammary gland of Capra hircus would significantly increase our understanding of the physiology of lactating mammary glands. In this study, two libraries were constructed using the lactating mammary gland tissues of Laoshan dairy goats (Capra hircus) during peak and late lactation. Solexa high-throughput sequencing technique and bioinformatics were used to determine the abundance and differential expression of the microRNAs between peak and late lactation. As a result, 19,044,002 and 7,385,833 clean reads were obtained, respectively, and 1,113 conserved known microRNAs and 31 potential novel microRNA candidates were identified. A total of 697 conserved microRNAs were significantly differentially expressed with a P-value<0.01, 272 microRNAs were up-regulated and 425 microRNAs were down-regulated during peak lactation. The results were validated using real-time quantitative RT-PCR. 762,557 annotated mRNA transcripts were predicted as putative target gene candidates. The GO annotation and KEGG pathway analysis suggested that differentially expressed microRNAs were involved in mammary gland physiology, including signal transduction, and cell-cell and cell-extracellular communications. This study provided the first global of the microRNA in Capra hircus and expanded the repertoire of microRNAs. Our results have great significance and value for the elucidation of complex regulatory networks between microRNAs and mRNAs and for the study of mammary gland physiology and lactation.

Circulating microRNAs as diagnostic biomarkers for cardiovascular diseases

One of the major challenges in cardiovascular disease is the identification of reliable clinical biomarkers that can be routinely measured in plasma. MicroRNAs (miRNAs) were recently discovered to circulate in the bloodstream in a remarkably stable form. Because of their stability and often tissue- and disease-specific expression and the possibility to measure them with high sensitivity and specificity, miRNAs are emerging as new diagnostic biomarkers. In this review we will provide an overview of the potential of circulating miRNAs as biomarkers for a wide range of cardiovascular diseases such as coronary artery disease, myocardial infarction, hypertension, heart failure, viral myocarditis, and type-2 diabetes mellitus. Furthermore, we will discuss the challenges with regard to further validation in large patient cohorts, and we will discuss how the measurement of multiple miRNAs simultaneously might improve the accuracy of the diagnostic test.