Skip to main content
  • Research article
  • Open access
  • Published:

Comparative miRNAome analysis revealed different miRNA expression profiles in bovine sera and exosomes

Abstract

Background

Extra-cellular components, such as serum and exosome, have drawn great attention as a readily accessible source of biomarkers for mammalian health. However, the contribution of different blood components to the signature of respective microRNAs (miRNAs) remains unknown, especially in cattle. In this study we aimed to investigate the miRNAs from bovine sera and exosomes, and to provide insights into their future applications.

Results

Blood collected from four healthy dairy cows were used for this study. The serum and exosomal RNAs were extracted using two commonly used commercial kits (Norgen and Invitrogen total RNA isolation kit), respectively. The miRNA profiles were then generated using RNA-seq. Sera had higher complexity of miRNAome consisting of 328 ± 17 miRNAs, while less number of miRNAs (260 ± 15, P = 0.001) was detected in exosomes. The profile of total detected miRNAs in sera and exosomes was different, while exosomes represented about 78 % of total miRNAs expressed in sera, suggesting that exosomes are the major miRNAs carriers in bovine sera. A total of 24 and 3 miRNAs (RPM > 5) were exclusively expressed in sera and exosomes, respectively. In addition, 12 miRNAs were differentially expressed between sera and exosomes (FDR < 0.05), with the expression of four of them being further validated by stem-loop RT-qPCR. Moreover, functional analysis showed that uniquely and highly expressed miRNAs in sera were mainly related to diseases and disorders, while the predicted functions of those in exosomes were enriched in tissue development and lipid metabolism.

Conclusion

Our results provide evidence that bovine sera and exosomes miRNAomes are different with regarding to the miRNA numbers, types and expressions. Based on their distinct profiles, miRNAs from sera and exosomes may reflect different aspects of physiological and pathological conditions in cattle. The functional analysis suggest that sera may be preferable for the purpose of detecting inflammation in cattle, while exosomes may be a better choice for monitoring the status of muscle development and lipid metabolism.

Background

As an important source of high quality protein for human consumption, cattle are considered to be one of the most important domestic animals. At present, intensive production systems have aggravated the health issues in cattle, such as mastitis [1], acidosis [2, 3], metritis [4], as well as welfare issues which have been leading to increasing foodborne and zoonotic risks to humans and economic losses for the industry. Due to the lack of validated early diagnosis procedures, most diseases and disorders are usually diagnosed at an advanced stage, and treatments for these pathologies are often expensive and ineffective. Therefore, there is an urgent need to identify reliable methods to diagnose the health status of cattle, which can then be used for early detection and monitoring disease progression.

Accumulating evidences suggest that expression patterns of microRNAs (miRNAs) in biofluids represent the in vivo status of many physiological changes and diseases [57], revealing that they can be served as diagnostic markers for multiple human diseases including cancers [811]. Blood is a non-invasive and the easiest obtained biofluid, and miRNAs in blood hold great promise to discover biomarkers for a wide range of diseases and biological processes [1215]. Whole blood was the most frequently used biofluid in detecting miRNAs, but the outcomes could be biased due to the complexity of various cell types and components [16, 17]. Recently, the miRNAs in body fluids, such as plasma, serum, urine, saliva and sputum [1820] have been used as informative biomarkers to assess and monitor the body’s physiological and pathological status due to their stability even under extreme conditions in human [21, 22]. In addition, it has been shown that the circulating miRNAs from plasma and serum could serve as potential biomarkers for livestock health and disease, such as miR- 26a for cattle early pregnancy [23], miR-19a and miR-19b for cattle heat stress [24], miR-29c and miR-375 for chicken puberty onset [25], and miR-122 for pig cardiogenic shock [26].

These free circulating miRNAs can be protected from nucleases by various types of carriers [27], such as exosomes. Exosomes are the most studied carriers, which are small (30–90 nm) and derived from the multivesicular body-sorting pathway [2830]. Recent researches have proposed to use exosomal miRNAs for diagnositic markers in human diseases [31], since the miRNAs in exosomes have specific function and higher variability than blood cells [16]. In addition, the quantity of miRNAs in exosomes exhibited more difference between healthy individuals and cancer patients than that in sera [27]. However, the exosomal miRNA profiling needs extra steps in RNA extraction and which miRNAs (sera vs. exosomes) are more representative for physiological and health changes in cattle have not been defined.

The next-generation sequencing has made it possible to obtain highly detailed information of miRNAs on the types and abundance from various biomaterials [32]. However, consensus has not been reached with regard to the sample types used for isolation of total RNA (such as sera or exosomes), especially in cattle. To date, there is little information on miRNAomes in bovine sera and exosomes, which could be potential diagnostic biomarkers relating to cattle health. Therefore, the aim of the current study was to compare the bovine miRNAomes of sera and exosomes, and to provide insights into their future applications.

Results

MiRNA profiles of sera and exosomes

An average of 9.25 ± 1.44 and 15.70 ± 12.43 ng small RNAs were obtained from sera and exosomes (Table 1), respectively. RNA sequencing resulted in 29,692,695 reads for sera and 6,581,761 reads for exosomes, respectively. After quality and length filter, 22,745,381 reads (76.6 %) in sera and 3,960,646 reads (60.2 %) in exosomes were used for further analysis (Additional file 1: Table S1). In exosomes, higher proportion of reads failed the size trimming compared to sera (39.9 % vs. 23.8 %). After mapping to bovine genome (Baylor Btau 4.6.1), the proportion of annotated miRNAs was 6.9 % in sera, and 11.3 % in exosomes (Fig. 1a). And the reads mapped to tRNA, rRNA snoRNA and snRNA were low (< 1 %) (Fig. 1a) from both sera and exosomes with large proportion reads unidentified. Size distribution of the reads between 19 and 40 nt revealed 2 peaks at 19–25 nt and 30–33 nt for both sample types (Fig. 1b).

Table 1 RNA extraction from bovine sera and exosomes
Fig. 1
figure 1

Profiling of small RNAs in cattle sera and exosomes samples. a The relative abundance of different classes of small RNAs. b Size and frequency distribution of detected small RNAs (19-40 nt). SE Sera, EX Exosomes

Complexity and specificity of sera and exosomes miRNAomes

The miRNAs that were detected in at least two cattle with more than 1 reads per million total mapped reads (RPM) in sera or exosomes were considered as expressed miRNAs. Sera had higher number of miRNAs (328 ± 17) expressed, while significantly less miRNAs (260 ± 15, P = 0.001) were expressed in exosomes (Fig. 2a). The respective miRNAomes of sera and exosomes are available in Additional file 2: Table S2. Using principle component analysis (PCA) it showed that the miRNAs profiles between sera and exosomes were different (Fig. 2b). More uniquely expressed miRNAs were found in sera than those in exosomes, with 24 and 3 uniquely expressed miRNAs detected (with mean RPM > 5) in sera and exosomes, respectively (Table 2).

Fig. 2
figure 2

Complexity of miRNAs detected in the sera and exosomes. a Average number of miRNAs expressed in sera and exosomes, different letters (a, b) represent significant difference (P < 0.05). b Principle component analysis of the total detected miRNAs in sera and exosomes (blue dots represent sera, red dots represent exosomes). c Venn diagram showing the profile of miRNAs expressed in sera and exosomes. d Correlation analysis of the commonly expressed miRNAs in sera and exosomes

Table 2 Unique miRNAs expressed in bovine sera and exosomes

It was observed that 282 miRNAs were commonly expressed in sera and exosomes (Fig. 2c), and their expression profiles had a high correlation of 0.97 (Fig. 2d). From the commonly expressed miRNAs, the proportion of the top 10 miRNAs accounted 95.6 and 96.2 % of total reads for sera and exosomes, respectively, with the same predominant miRNAs identified. The comparison of the relative abundance of top 10 annotated miRNAs showed that sera and exosomes had similar profiles of abundant miRNAs except miR-191 in sera and miR-10b in exosmes (Fig. 3a, b).

Fig. 3
figure 3

The pie chart of ten most abundant miRNAs. a Sera. b Exosomes. MiRNAs differed in two samples were colored in red

In addition, 12 miRNAs were differentially expressed (DE) between sera and exosomes. The expression of miR-361, miR-16a, miR-146b and miR-24-3p were higher in sera, while the expression of other eight miRNAs (miR-328, miR-1249, miR-1306, miR-339b, miR-339a, miR-296-5p, miR-21-3p and miR-92b) were higher in exosomes (FDR < 0.05, Table 3).

Table 3 Differentially expressed miRNAs in sera and exosomes

Experimental validation of miRNA expression

A total of 10 uniquely expressed and DE miRNAs that identified from miRNA-seq (four uniquely expressed miRNAs in sera: miR-196b, miR-196a, miR-18b, and miR-9-3p; two uniquely expressed miRNAs in exosomes: miR-2284k and miR-9-5p; four DE miRNAs between sera and exosomes: miR-146b, miR-24-3p, miR-328, and miR-21-3p) were selected for validation using stem-loop RT-qPCR. Three out of four DE miRNAs showed similar trend as detected by miRNA-seq. Specifically, the expressions of miR-146b and miR-24-3p were highly expressed in the sera, while miR-328 was highly expressed in the exosomes by miRNA-seq and RT-qPCR (P < 0.05) (Fig. 4 a-c). Although the expression of miR-21-3p, detected by RT-qPCR showed no significant difference between sera and exosomes (P < 0.1), its expression was numerically higher in exosomes, which was similar with the results detected by miRNA-seq (Fig. 4d). However, the expression of those sera or exosomes uniquely expressed miRNAs could not be detected by RT-qPCR (CT value < 35), which may be due to their low abundance and the relative low content of miRNAs in serum.

Fig. 4
figure 4

Expression of selected DE miRNAs between sera and exosomes detected by RT-qPCR and miRNA-seq. a Expression of miR-146b; b Expression of miR-24-3p; c Expression of miR-328; d Expression of miR-21-3p. miRNA expression from RT-qPCR represented by lines on the top and values are shown on the right vertical axis as relative abundance. miRNA expression from miRNA-seq represented by bars on the bottom and values are shown on the left vertical axis as log2RPM (normalized reads number). * on the top of lines or bars indicate significant (P < 0.05 or FDR < 0.05) difference between sera and exosomes. Data are presented as Mean ± Standard deviation

Stability of the bovine sera and exosomes miRNAomes in healthy animals

It is important to determine how miRNAomes varies among individuals, especially for the healthy animals. As the sera and exosomes were taken from the same four cattle, the individual variation between sera and exosomes reflected the stability of miRNAs in two types of samples. The high correlation was observed with r = 0.96 for sera (Additional file 3: Table S3) and r = 0.95 for exosomal miRNAomes (Additional file 4: Table S4) among four animals. The hierarchical clustering analysis on the 50 most variable miRNAs revealed that sera and exosomes were clustered to different groups (Fig. 5).

Fig. 5
figure 5

Heatmap generated by clustering of the 50 most variable miRNAs in sera (SE) and exosomes (EX) from four cattle. Colors represent different normalized sequencing reads number as indicated by the color bar (Log2RPM)

Functional categories and pathways of uniquely expressed and differentially expressed miRNAs in sera and exosomes

Based on predicted targets of miRNAs using Targetscan and mirRanda as well as functional analysis using Ingenuity Pathway Analysis (IPA), the functional categories of miRNAomes in sera and exosomes were identified. Further functional analysis of uniquely expressed miRNAs in sera and exosomes revealed 3,847 and 1,520 gene targets, respectively. Among the top 20 functions of predicted gene targets of the uniquely expressed miRNAs in sera and exosomes, 14 functional categories were common in two types of samples (Table 4). Besides the common functions, the sera uniquely expressed miRNAs targeted the function of diseases and disorders such as neurological/infections diseases and skeletal, muscular, hereditary and psychological disorders, while the exosomes uniquely expressed miRNAs packaged into exosomes targeted the function of tissue development and lipid metabolism (Table 4).

Table 4 Top 20 functional categories of predicted gene targets by uniquely expressed miRNAs in sera and exosomes

There were 3,183 genes putatively targeting by the 12 DE miRNAs. Functional analysis revealed that most of the functional categories were shared by the highly expressed DE miRNAs in sera or exosomes (Table 5). Additionally, the function of infectious disease and cell-mediated immune response were targeted by sera enriched miRNAs, while the skeletal/muscular system development and lipid metabolism were targeted by exosomes enriched miRNAs. Besides, the most targeted pathways predicted based on sera enriched miRNAs were LPS-stimulated MAPK signaling, while the exosomes enriched miRNAs targeted cancer pathways (Table 6).

Table 5 Top 20 functions of predicted gene targets of differentially expressed miRNAs in sera and exosomes
Table 6 Predicted pathways targeted by miRNAs differentially expressed in sera and exosomes

Discussion

To date, the searching of blood miRNA biomarkers mainly focused on the altered expression levels of miRNAs in diseases. However, little is known about their expression levels associated with the physiological changes and composition of miRNAome, especially in cattle. Here, we conducted a comprehensive miRNAome comparison between cattle sera and exosomes that have been commonly used in miRNA studies.

While surveying the profiles of cell-free small RNAs in cattle, we observed that the annotated miRNAs were low from the sequencing datasets. Size distribution of the total reads revealed two peaks in sera and exosomes, which is similar as reported in the blood samples of mouse and human [33, 34]. The nucleic acid fragments ranging 19–25 nt were mostly annotated as miRNAs. It has been reported that, in addition to miRNAs, several classes of small RNAs mediated multiple biological functions, including siRNAs (20–24 nt), piRNAs (25–32 nt), and others that bound by Argonaute proteins also could be detected in blood [3537]. Moreover, studies on the processed small RNAs, such as 5’ tRNA halves and and 5’ YRNA fragments (30–33 nt) suggested that they were abundant in circulating small RNAs and might act as novel form of signaling molecules [38, 39]. Thus, it is postulated that the 30–33 nt peak in our study might be processed small RNAs.

The miRNAs detected in exosomes represented about 78 % of total miRNAs expressed in sera, suggesting exosomes are the major carriers for miRNAs in bovine sera. In particular, there were 282 miRNAs expressed in all individuals with high correlation between sera and exosomes, indicating those miRNAs may play important roles in extra-cellular blood fluid of cattle. MiR-486, miR-92a, miR-16, miR-144, and miR-451 have been commonly applied as erythrocyte markers for human and mice [40, 41], and they were all expressed in our bovine sera and exosomes dataset with high proportion. Besides, miR-22-3p, the secondly abundant miRNA in sera (12.0 %) and exosomes (10.4 %), was miRNA biomarker expressed by myeloid [42]. Therefore, red blood cells and myeloid might contribute significantly to circulating miRNAs and this could have important implications for circulating miRNA biomarker studies. Some of top abundance miRNAs detected in bovine sera were not abundant ones in human [17], vice versa, which were not found in bovine sera. The inconsistence in miRNA spectrums might come from the different erythropoiesis between human and bovine [42].

At present, little is known about the origin and destination of these sera/exosomes miRNAs in the circulatory system. Based on identified miRNA profiles, our results showed that spectra of miRNAomes were different between sera and exosomes, and miRNAome of sera was more complex than that of exosomes. The complexity of sera miRNAomes is reasonable, because circulating miRNAs were bounded to various carriers including exosomes, proteins, lipoprotein particles, and cell-derived extracellular vesicles [28, 29, 43]. For example, it has been reported that Argonaute2 complexes carried a population of circulating miRNAs independent of vesicles in human plasma [44]. It is noticeable that the sera miRNAomes detected in this study may also contain the exosome miRNAs. To precisely identify sera specific miRNAs, sera-depletion from exosomes using size-exclusion chromatography in the future studies will provide the detailed contribution and function of the different circulating miRNAs components in sera. In addition, the unique origin may contribute to the miRNAs patterns in sera and exosomes. The assumption above could be supported by the fact that 24 and 3 miRNAs were exclusively expressed with high level in sera and exosomes, respectively. In addition, 12 miRNAs were differentially expressed between two samples. Because of their distinct profiles, miRNAs from sera and exosomes might reflect different aspects of physiological and pathological conditions in cattle.

Many miRNAs are evolutionarily conserved [45]. Therefore, information from human studies on the highly conserved miRNAs that have important roles in regulating health may also be applied to cattle. For example, miR-196, which was uniquely expressed in sera, has recently been reported to be involved in the antiviral immunity in human [46]. And antiviral interferon response could induce the expression of miR-196, which then inhibited the gene expression and replication of hepatitis C virus [47, 48]. Mitogen activated protein kinases (MAPK) are a group of signaling molecules that play important roles in the inflammatory processes [4951]. Together with LPS-stimulated MAPK pathway that is the most enriched for up-regulated miRNAs in sera, our results suggest that sera uniquely expressed miRNAs (such as miR-196) may play important roles in regulating inflammation and immune response in cattle. MiR-9 is a brain specific miRNA [52], and we found that miR-9-3p and miR-9-5p was uniquely expressed in bovine sera and exosomes, respectively. The differential enrichment of miR-9-3p and miR-9-5p further demonstrated the specific function of sera and exosomes. In addition, miR-92b was differentially expressed between sera and exosomes, with higher expression in exosomes. Nass et al. reported that miR-92b was specifically expressed in the heart and muscle, and overexpression of miR-92b reduced Mef2 levels and caused muscle attachment defects [53]. In addition, miR-92b and miR-9 were over-expressed in developing brain and neuronal stem cells compared to the adult brain, and thus have been implicated as players in human nervous system development [54, 55]. These suggest that exosomes uniquely expressed and/or highly expressed miRNAs may be used as markers to study changes in development.

The high correlation of expression profiles among individuals was observed, demonstrating the low individual variation of miRNAs in bovine extra-cellular blood fluid, the reasonable approach of sera/exosomes sample preparation and miRNA profiling in this study. Overall, the identified expression of miRNAs in sera and exosomes enabled an insight into the basis miRNAs pattern of healthy cattle, which provides a strong basis to inform the further investigation of circulating miRNAomes as biomarkers in cattle.

Conclusion

Using RNA-seq, our results revealed that 1) exosomes showed a less complex with respect to the miRNAomes as compared to sera; 2) sera and exosomes had their specificities in profiling of miRNAs, which may reflect different aspects of physiological and pathological conditions in cattle; 3) exosomes were major miRNAs carriers in cattle sera, and both samples could be used as potential candidates in biomarker study for different goals in future. Nevertheless, this study laid a foundation for understanding the dysregulations diverging from ordinary physiology by providing information about miRNAomes of healthy individual of bovine sera and exosomes. Furthermore, our results provide the guidance in selection of sample type for particular goals in cattle miRNAs study. Based on our preliminary findings, sera may be preferable for selecting biomarkers of heath such as disease, while exosomes may be better for monitoring the muscle development and lipid metabolism, which should be validated in the future studies using disease and production deficient animal models.

Methods

Blood sample collection

The whole blood samples were collected from 4 healthy lactating Holstein cow (33 ± 7 mon, 654 ± 54 kg; mean ± SD). All cows had normal milk production and intake and did not have any incidence of health related issues (data not shown). Experimental protocols were reviewed and approved by the Livestock Animal Care Committee of the University of Alberta and all procedures were conducted following the guidelines of Canadian Council on Animal Care. Briefly, 9 mL whole blood was collected in S-Movovette Sera-Gel tubes (Sarstedt, Nümbrecht, Germany) containing clot activator, and then the samples were incubated at room temperature for 1 h. The whole blood samples were firstly centrifuged at 4 °C (1,900 × g, 10 min), followed by a second centrifugation at 4 °C (16,000 × g, 10 min). After this step, the circulating cell-free nucleic acid were in the supernatant (sera) and then the sera were transferred into a fresh 2.0 mL tube and stored at -80 °C for further analysis.

Total RNA isolation

Total RNA in sera was extracted using Total RNA Purification kit (Norgen, Thorold, CA) following manufacturer’s protocols. For exosomal RNA extraction, exosomes were firstly isolated from the same sera samples using the Total Exosomes Isolation kit (Invitrogen, Carlsbad, CA), and then the total RNA was extracted using Total Exosome RNA and Protein Isolation Kit (Invitrogen, Carlsbad, CA). The concentration of RNA was determined using the Qubit microRNA assay kit and the 2.0 Fluorometer (Invitrogen, Carlsbad, CA). The RNA samples were stored at -80 °C for further analysis.

Library construction and sequencing

Total RNA (5 μL) from each sample was used to construct miRNA library using the TruSeq Small RNA Sample Preparation Kit (Illumina, San Diego, CA) according to the manufacturer’s instruction. Small RNA libraries were then pooled together in equal volumes for gel purification. The pooled library was sequenced at Génome Québec (Montréal, Canada) using the Illumina HiSeq 2000 system (Illumina) as 50 bp single reads.

Small RNA sequence analysis

Low-quality reads were removed from raw data using CASAVA 1.8 based on chastity. High-quality reads were then subjected to sRNAbench for adaptor sequence trimming and reads length distribution analysis [56]. The sequences with read length larger than 15 nt were aligned against bovine miRNA database (miRBase, release version 21) with the default parameters to identify known miRNAs using sRNAbench. Each library was processed separately. The expression level of miRNAs in each library was estimated by sRNAbench, which normalized reads count number of each miRNA RPM by the following formula: RPM = (miRNA reads number/total mapped reads per library) × 1,000,000. Correlations were calculated using the Pearson method in R software (version 3.0). The differentially expressed (DE) miRNAs were investigated by using bioinformatics tool edgeR [57]. The DE miRNAs were determined by Logfold change (FC) >1 or < -1 and false discovery rate (FDR) < 0.05 based on Benjamini and Hochberg multiple testing correction [58].

Target gene prediction and functional analysis

The functional analysis of selected miRNAs was performed following a previous study [59]. Briefly, target genes of selected miRNAs were commonly predicted by TargetScan Relase 6.0 (http://www.targetscan.org/) [60] and miRanda (http://www.microrna.org/microrna/home.do). Then the predicted target genes were analyzed through IPA using core analysis module (www.ingenuity.com). A threshold of FDR < 0.05 and molecule number > 1 were used to enrich significant biological functions of each miRNA.

Experimental validation of miRNA expression using stem-loop RT-qPCR

The expression of uniquely expressed and DE miRNAs that were identified by miRNA-seq were validated by stem-loop RT-qPCR using TAQMAN miRNA assays following the manufacturer’s instruction (Applied Biosystems, Foster City, CA). Briefly, cDNAs were reverse transcribed from 5 μL of total RNA using 5 X specific miRNA RT primer and then were amplified using a 20 X TAQMAN miRNA assay. Fluorescence signal was detected with StepOnePlus™ Real-Time PCR System (Applied Biosystems, Foster City, CA). In this study, miR-93 was used as an internal reference to calculate the relative expression of target miRNA using the 2-ΔΔCT (cycle threshold, CT) method [61]. This miRNA was selected due to its stable expression between sera and exosomes after comparing with the most used reference U6 [62, 63], and previously reported serum miRNA reference candidates miR-127 and miR-744 in bovine and mouse [64, 65]. In addition, no difference was observed between sera and exosomes in miR-93 expression detected by RNA-seq (P > 0.05). Differences were considered statistically different at P < 0.05 and analysis were performed using t-test.

References

  1. Seegers H, Fourichon C, Beaudeau F. Production effects related to mastitis and mastitis economics in dairy cattle herds. Vet Res. 2002;34(5):475–91.

    Article  Google Scholar 

  2. Schwaiger T, Beauchemin KA, Penner GB. The duration of time that beef cattle are fed a high-grain diet affects the recovery from a bout of ruminal acidosis: dry matter intake and ruminal fermentation. J Anim Sci. 2013;91(12):5729–42.

    Article  CAS  PubMed  Google Scholar 

  3. Wierenga K, McAllister TA, Gibb DJ, Chaves AV, Okine EK, Beauchemin KA, Oba M. Evaluation of triticale dried distillers grains with solubles as a substitute for barley grain and barley silage in feedlot finishing diets. J Anim Sci. 2010;88(9):3018–29.

    Article  CAS  PubMed  Google Scholar 

  4. Vergara CF, Döpfer D, Cook NB, Nordlund KV, McArt JAA, Nydam DV, Oetzel GR. Risk factors for postpartum problems in dairy cows: explanatory and predictive modeling. J Dairy Sci. 2014;97(7):4127–40.

    Article  CAS  PubMed  Google Scholar 

  5. Gilad S, Meiri E, Yogev Y, Benjamin S, Lebanony D, Yerushalmi N, et al. Sera microRNAs are promising novel biomarkers. PLoS One. 2007;3(9):e3148.

    Article  Google Scholar 

  6. Wang K, Zhang S, Marzolf B. Circulating microRNAs, potential biomarkers for drug-induced liver injury. Proc Natl Acad Sci U S A. 2009;106(11):4402–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Qin X, Xu H, Gong W, Deng W. The Tumor Cytosol miRNAs, Fluid miRNAs, and Exosomes miRNAs in Lung Cancer. Front Oncol. 2013;4:357.

    Google Scholar 

  8. Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.

    Article  CAS  PubMed  Google Scholar 

  9. Qin X, Yan L, Zhao X, Li C, Fu Y. microRNA-21 overexpression contributes to cell proliferation by targeting PTEN in endometrioid endometrial cancer. Oncol Lett. 2012;4(6):1290–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Kim YJ, Hwang SH, Cho HH, Shin KK. MicroRNA 21 regulates the proliferation of human adipose tissue derived mesenchymal stem cells and high‐fat diet‐induced obesity alters microRNA 21 expression in white adipose tissues. J Cell Physiol. 2012;227(1):183–93.

    Article  CAS  PubMed  Google Scholar 

  11. Zhang Z, Sun J, Bai Z, Li H, He S, Chen R. MicroRNA-153 acts as a prognostic marker in gastric cancer and its role in cell migration and invasion. Onco Targets Ther. 2015;8:357–64.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Burgos KL, Javaherian A, Bomprezzi R, Ghaffari L, Rhodes S, Courtright A, et al. Identification of extracellular miRNA in human cerebrospinal fluid by next-generation sequencing. RNA. 2013;19(5):712–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wu F, Guo NJ, Tian H, Marohn M, Gearhart S, Bayless TM, et al. Peripheral blood microRNAs distinguish active ulcerative colitis and Crohn’s disease. Inflamm Bowel Dis. 2011;17(1):241–50.

    Article  PubMed  Google Scholar 

  14. Hoekstra M, van der Lans CAC, Halvorsen B, Gullestad L, Kuiper J, Aukrust P, et al. The peripheral blood mononuclear cell microRNA signature of coronary artery disease. Biochem Biophys Res Commun. 2010;394(3):792–7.

    Article  CAS  PubMed  Google Scholar 

  15. Meder B, Keller A, Vogel B, Haas J, Sedaghat-Hanmedani F, Kayvanpour E, et al. MicroRNA signatures in total peripheral blood as novel biomarkers for acute myocardial infarction. Basic Res Cardiol. 2011;106(1):13–23.

    Article  CAS  PubMed  Google Scholar 

  16. Leidinger P, Backes C, Meder B, Meese E, Keller A. The human miRNA repertoire of different blood compounds. BMC Genomics. 2014;15:474.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Cheng L, Sharples RA, Scicluna BJ, Hill AF. Exosomes provide a protective and enriched source of miRNA for biomarker profiling compared to intracellular and cell-free blood. J Extracell Vesicles. 2014;3:23743.

    Article  Google Scholar 

  18. Zubakov D, Boersma AW, Choi Y, van Kuijk PF, Wiemer EA, Kayser M. MicroRNA markers for forensic body fluid identification obtained from microarray screening and quantitative RT-PCR confirmation. Int Legal Med. 2010;124(3):217–26.

    Article  Google Scholar 

  19. Weber JA, Baxter DH, Zhang S, Huang DY, Huang KH, Lee MJ, et al. The microRNA spectrum in 12 body fluids. Clin Chem. 2010;56(11):1733–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Tzimagiorgis G, Michailidou EZ, Kritis A, Markopoulos AK, Kouidou S. Recovering circulating extracellular or cell-free RNA from bodily fluids. Cancer Epidemiol. 2011;35(6):580–9.

    Article  CAS  PubMed  Google Scholar 

  21. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A. 2008;105(30):10513–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Chen X, Ba Y, Ma L, Cai X, Yin Y, Wang K, et al. Characterization of microRNAs in sera: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008;18(10):997–1006.

    Article  CAS  PubMed  Google Scholar 

  23. Ioannidis J, Donadeu FX. Circulating miRNA signatures of early pregnancy in cattle. BMC Genomics. 2016;17:184.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Zheng Y, Chen KL, Zheng XM, Li HX, Wang GL. Identification and bioinformatics analysis of microRNAs associated with stress and immune response in serum of heat-stressed and normal Holstein cows. Cell Stress Chaperones. 2014;19(6):973–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Han W, Zhu Y, Su Y, Li G, Qu L, Zhang H, et al. High-throughput sequencing reveals circulating miRNAs as potential biomarkers for measuring puberty onset in chicken (Gallus gallus). PLoS One. 2016;11(5):e0154958.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Andersson P, Gidlöf O, Braun OO, Götberg M, van der Pals J, Olde B, et al. Plasma levels of liver-specific miR-122 is massively increased in a porcine cardiogenic shock model and attenuated by hypothermia. Shock. 2012;37(2):234–8.

    Article  CAS  PubMed  Google Scholar 

  27. Ashby J, Flack K, Jimenez LA, Duan Y, Khatib A-KK, Somlo G, et al. Distribution profiling of circulating microRNAs in sera. Anal Chem. 2014;86(18):9343–9.

    Article  CAS  PubMed  Google Scholar 

  28. Raposo G, Nijman HW, Stoorvogel W, Liejendekker R, Harding CV, Melief CJ, et al. B lymphocytes secrete antigen-presenting vesicles. J Exp Med. 1996;183(3):1161–72.

    Article  CAS  PubMed  Google Scholar 

  29. Denzer K, van Eijk M, Kleijmeer MJ, Jakobson E, de Groot C, Geuze HJ. Follicular dendritic cells carry MHC class II-expressing microvesicles at their surface. J Immunol. 2000;165(3):1259–65.

    Article  CAS  PubMed  Google Scholar 

  30. Hu G, Drescher KM, Chen XM. Exosomal miRNAs: biological properties and therapeutic potential. Front Genet. 2011;3:56.

    Google Scholar 

  31. Kourembanas S. Exosomes: vehicles of intercellular signaling, biomarkers, and vectors of cell therapy. Annu Rev Physiol. 2015;77:13–27.

    Article  CAS  PubMed  Google Scholar 

  32. Pritchard CC, Cheng HH, Tewari M. MicroRNA profiling: approaches and considerations. Nat Rev Gene. 2012;13(5):358–69.

    Article  CAS  Google Scholar 

  33. Dhahbi JM, Spindler SR, Atamna H, Boffelli D, Mote P, Martin DI. 5’-YRNA fragments derived by processing of transcripts from specific YRNA genes and pseudogenes are abundant in human sera and plasma. Physiol Genomics. 2013;45(21):990–8.

    Article  CAS  PubMed  Google Scholar 

  34. Dhahbi JM, Spindler SR, Atamna H, Yamakawa A, Boffelli D, Mote P, et al. 5’ tRNA halves are present as abundant complexes in sera, concentrated in blood cells, and modulated by aging and calorie restriction. BMC Genomics. 2012;14:298.

    Article  Google Scholar 

  35. Okamura K. Diversity of animal small RNA pathways and their biological utility. Wiley Interdiscip Rev RNA. 2012;3(3):351–68.

    Article  CAS  PubMed  Google Scholar 

  36. Zhang C. Novel functions for small RNA molecules. Curr Opin Mol Ther. 2009;11(6):641–51.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Esteller M. Non-coding RNAs in human disease. Nat Rev Genet. 2011;12(12):861–74.

    Article  CAS  PubMed  Google Scholar 

  38. Kramer BS. The science of early detection. Urol Oncol. 2004;22(4):344–7.

    Article  PubMed  Google Scholar 

  39. Mestdagh P, Hartmann N, Baeriswyl L, Andreasen D, Bernard N, Chen C, et al. Evaluation of quantitative miRNA expression platforms in the microRNA quality control (miRQC) study. Nat Methods. 2014;11(8):809–15.

    Article  CAS  Google Scholar 

  40. Rasmussen KD, Simmini S, Abreu-Goodger C, Bartonicek N, Di Giacomo M, Bilbao-Cortes D, et al. The miR-144/451 locus is required for erythroid homeostasis. J Experi Med. 2010;207(7):1351–8.

    Article  CAS  Google Scholar 

  41. Pritchard CC, Kroh E, Wood B, Arroyo JD, Dougherty KJ, Miyaji MM, et al. Blood cell origin of circulating microRNAs: A cautionary note for cancer biomarker studies. Cancer Pre Res (Phila). 2012;5(3):492–7.

    Article  CAS  Google Scholar 

  42. Spornraft M, Kirchner B, Pfaffl MW, Riedmaier I. Comparison of the miRNAome and piRNome of bovine blood and plasma by small RNA sequencing. Biotechnol Lett. 2015;37(6):1165–76.

    Article  CAS  PubMed  Google Scholar 

  43. Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011;13(4):423–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci U S A. 2011;108(12):5003–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Siomi H, Siomi MC. Posttranscriptional regulation of MicroRNA biogenesis in animals. Mol Cell. 2010;38:323–32.

    Article  CAS  PubMed  Google Scholar 

  46. Mahajan VS, Drake A, Chen J. Virus-specific host miRNAs: antiviral defenses or promoters of persistent infection? Trends Immunol. 2009;30(1):1–7.

    Article  CAS  PubMed  Google Scholar 

  47. Hou W, Tian Q, Zheng J, Bonkovsky HL. MicroRNA‐196 represses Bach1 protein and hepatitis C virus gene expression in human hepatoma cells expressing hepatitis C viral proteins. Hepatology. 2010;51(5):1494–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Pedersen IM, Cheng G, Wieland S, Volinia S, Croce CM, Chisari FV, et al. Interferon modulation of cellular microRNAs as an antiviral mechanism. Nature. 2007;449(7164):919–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Ajizian SJ, English BK, Meals EA. Specific inhibitors of p38 and extracellular signal-regulated kinase mitogen-activated protein kinase pathways block inducible nitric oxide synthase and tumor necrosis factor accumulation murine macrophages stimulated with lipopolysaccharide and interferon-gamma. J Infect Dis. 1999;179(4):939–44.

    Article  CAS  PubMed  Google Scholar 

  50. Haddad EB, Birrell M, McCluskie K, Ling A, Webber SE, Foster ML, et al. Role of p38 MAP kinase in LPS-induced airway inflammation in the rat. Br J Pharmacol. 2001;132(8):1715–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Schmeck B, Zahlten J, Moog K, van Laak V, Huber S, Hocke AC, et al. Streptococcus pneumoniae-induced p38 MAPK-dependent phosphorylation of RelA at the interleukin-8 promotor. J Biol Chem. 2004;279(51):53241–7.

    Article  CAS  PubMed  Google Scholar 

  52. Jung HJ, Coffinier C, Choe Y, Beigneux AP, Davies BS, Yang SH, et al. Regulation of prelamin A but not lamin C by miR-9, a brain-specific microRNA. Proc Natl Acad Sci U S A. 2012;109(7):E423–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Nass D, Rosenwald S, Meiri E, Gilad S, Tabibian-Keissar H, Schlosberg A, et al. MiR-92b and miR‐9/9* are specifically expressed in brain primary tumors and can be used to differentiate primary from metastatic brain tumors. Brain Pathol. 2009;19(3):375–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Zhao JJ, Hua YJ, Sun DG, Meng XX, Xiao HS, Ma X. Genome-wide microRNA profiling in human fetal nervous tissues by oligonucleotide microarray. Childs Nervs Syst. 2006;22(11):1419–25.

    Article  Google Scholar 

  55. Kapsimali M, Kloosterman WP, de Bruijn E, Rosa F, Plasterk RH, Wilson SW. MicroRNAs show a wide diversity of expression profiles in the developing and mature central nervous system. Genome Biol. 2007;8(8):R173.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Barturen G, Rueda A, Hamberg M. sRNAbench: profiling of small RNAs and its sequence variants in single or multi-species high-throughput experiments. Meth Next Generation Sequencing. 2014;1(1):2084–7173.

    Article  Google Scholar 

  57. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    Article  CAS  PubMed  Google Scholar 

  58. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57(1):289–300.

    Google Scholar 

  59. Liang G, Malmuthuge N, McFadden TB, Bao H, Griebel PJ, Stothard P, et al. Potential regulatory role of microRNAs in the development of bovine gastrointestinal tract during early life. PLoS One. 2014;9(3):e92592.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Enright AJ, John B, Gaul U, Tuschl T, Sander C. MicroRNA targets in Drosophila. Genome Biol. 2003;5(1):R1.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  62. Ng EK, Chong WW, Jin H, Lam EK, Shin VY, Yu J, et al. Differential expression of microRNAs in plasma of patients with colorectal cancer: a potential marker for colorectal cancer screening. Gut. 2009;58(10):1375–81.

    Article  CAS  PubMed  Google Scholar 

  63. Tsujiura M, Ichikawa D, Komatsu S, Shiozaki A, Takeshita H, Konishi H, et al. Circulating microRNAs in plasma of patients with gastric cancers. Br J Cancer. 2010;102(7):1174–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Mi QS, Weiland M, Qi RQ, Gao XH, Poisson LM, Zhou L. Identification of mouse serum miRNA endogenous references by global gene expression profiles. PLoS One. 2012;7(2):e31278.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Bae IS, Chung KY, Yi J, Kim TI, Choi HS, Cho YM, et al. Identification of reference genes for relative quantification of circulating microRNAs in bovine serum. PLoS One. 2015;10(3):e0122554.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgement

This study is supported by Alberta Livestock Meat Agency (ALMA 2011F129R), NSERC Discovery grant for L.L. Guan, Alberta Innovates Doctoral Graduate Student Scholarship and Chinese Council Scholarship for G. Liang and International Postdoctoral Exchange Fellowship Program (20140020) for K. Zhao.

Availability of supporting data

All the sequencing data sets supporting the results in the current study have been deposited in the publicly available NCBI’s Gene Expression Omnibus Database (GSE72643) (http://www.ncbi. nlm.nih.gov/geo/).

Authors’ contributions

LLG, GL and KZ designed the experiments and wrote the paper. XS collected the samples and constructed the libraries. KZ and GL did the data analysis. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Le Luo Guan.

Additional files

Additional file 1: Table S1.

Compilation of evaluated reads data in the process of generating miRNAs. (DOCX 70 kb)

Additional file 2: Table S2.

Excel spreadsheet containing total miRNAs expressed in sera and exosomes. Expression data for the miRNAs were normalized to RPM, and the miRNAs that were detected in at least 2 animals in each group were considered expressed in this study. (DOCX 206 kb)

Additional file 3: Table S3.

Pearson correlation of the detected miRNAs among individuals in sera. (DOCX 38 kb)

Additional file 4: Table S4.

Pearson correlation of the detected miRNAs among individuals in exosomes. (DOCX 38 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhao, K., Liang, G., Sun, X. et al. Comparative miRNAome analysis revealed different miRNA expression profiles in bovine sera and exosomes. BMC Genomics 17, 630 (2016). https://doi.org/10.1186/s12864-016-2962-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-016-2962-1

Keywords