- Research article
- Open Access
Comprehensive analysis of blood cells and plasma identifies tissue-specific miRNAs as potential novel circulating biomarkers in cattle
BMC Genomicsvolume 19, Article number: 243 (2018)
The potential of circulating miRNAs as biomarkers of tissue function, both in health and disease, has been extensively demonstrated in humans. In addition, circulating miRNA biomarkers offer significant potential towards improving the productivity of livestock species, however, such potential has been hampered by the absence of information on the nature and source of circulating miRNA populations in these species. In addition, many miRNAs originally proposed as robust biomarkers of a particular tissue or disease in humans have been later shown not to be tissue specific and thus to actually have limited biomarker utility. In this study, we comprehensively analysed miRNA profiles in plasma and cell fractions of blood from cattle with the aim to identify tissue-derived miRNAs which may be useful as biomarkers of tissue function in this important food animal species.
Using small RNA sequencing, we identified 92 miRNAs with significantly higher expression in plasma compared to paired blood cell samples (n = 4 cows). Differences in miRNA levels between plasma and cell fractions were validated for eight out of 10 miRNAs using RT-qPCR (n = 10 cows). Among miRNAs found to be enriched in plasma, we confirmed miR-122 (liver), miR-133a (muscle) and miR-215 (intestine) to be tissue-enriched, as reported for other species. Profiling of additional miRNAs across different tissues identified the human homologue, miR-802, as highly enriched specifically in liver.
These results provide novel information on the source of bovine circulating miRNAs and could significantly facilitate the identification of production-relevant tissue biomarkers in livestock. In particular, miR-802, a circulating miRNA not previously identified in cattle, can reportedly regulate insulin sensitivity and lipid metabolism, and thus could potentially provide a specific biomarker of liver function, a key parameter in the context of post-partum negative energy balance in dairy cows.
MicroRNAs (miRNAs) are short, non-coding RNA molecules which are primarily involved in the post-transcriptional fine-tuning of gene expression . Since the discovery of lin-4 in Caenorhabditis elegans in 1993, thousands of miRNAs have been identified and characterised [2, 3]. It is now known that miRNAs are involved in the regulation of virtually all developmental processes, from the embryo to the adult [4,5,6,7], as well as in disease .
The discovery of cancer-associated miRNAs together with the detection of miRNAs in bio fluids  sparked a wave of research into the potential use of extracellular miRNAs as biomarkers of disease. A miRNA biomarker can be defined as a miRNA that is specifically produced or enriched in a given tissue, the circulating levels of which may reflect pathological or physiological changes in said tissue. In recent years, the value of circulating miRNAs as diagnostic biomarkers has been shown in relation to cancer (e.g. miR-21, miR-20, miR-221 [10, 11]), cardiovascular disease (miR-1, miR-133a ), liver disease (miR-122 ) and diabetes (miR-375 and miR-34 ), among many other human pathologies and physiological processes including pregnancy .
However, recent work has urged caution in interpreting the results of miRNA biomarker studies [16, 17]. Specifically, a large number of miRNAs initially proposed as cancer biomarkers were later shown not to be tissue-specific, contrary to what would be expected from a robust tissue biomarker. Many of those miRNAs were found to be present in blood cells, which are a major contributor of miRNAs in circulation , and their circulating levels were shown to be affected by changes in blood cell numbers. In addition, a review of non-cancer, disease-related miRNA biomarkers found that the majority of proposed biomarkers were not tissue-specific and, moreover, were not expressed in the disease-relevant tissues . These reports highlight the need for further studies to identify the source of proposed miRNA biomarkers and determine their enrichment in the tissue(s) of interest.
The potential of miRNA biomarkers in food-producing species has also been reported, although in a much smaller number of studies than in humans. The identification of biomarkers associated with specific production traits and disease conditions in cattle would be of significant agricultural importance. A handful of studies have already linked circulating miRNA profiles with pregnancy and the oestrous cycle [19, 20], infection [21, 22] and grazing  in cattle, as well as with bacterial disease and feed deprivation in other domestic species . However, the source and tissue-specificity of circulating miRNA populations in these species have not been determined which, in view of the above, limits the practical significance of those results.
To address these limitations and facilitate the discovery in cattle of circulating bona-fide miRNA biomarkers, we sought to profile miRNA populations in different blood fractions in order to distinguish miRNAs expressed in blood cells from those expressed in other body tissues which could provide novel specific biomarkers. Using next-generation sequencing, we provide a list of potential tissue-enriched miRNAs in circulation including miR-802, a previously uncharacterised bovine liver-specific miRNA which may prove useful as a disease biomarker in high-producing cows.
Blood samples were collected in EDTA-coated tubes (Becton Dickinson, USA) from the jugular vein of six non-pregnant Holstein-Friesian cross cows, aged between 20 and 25 months and held at Langhill Farm (University of Edinburgh), with the relevant ethics approval (see Declarations). From each blood sample, a plasma fraction was obtained using a two-step centrifugation protocol  to separate blood cells  followed by filtration through 0.2 μm syringe filters (Sartorius, Germany) to ensure the removal of any residual cellular debris. Subsequently, a cell fraction was obtained from a different volume (100 μL) of the same blood sample by centrifugation at 1900 x g for 10 min at 4 °C. The resulting supernatant was discarded and the cell pellet (consisting of red blood cell and buffy coat fractions) was re-suspended in RNase-free water up to a volume of 250 μL. The blood cell solution was immediately used for RNA extraction as described below.
Bovine tissue samples including brain, heart, intestine, kidney (cortex and medulla combined), liver, lung, skeletal muscle, ovary, placenta (placentome and interplacentome combined), skin, spleen and uterus were collected from a local abattoir in ice-cold PBS. Upon arrival in the laboratory tissues were dissected, labelled, weighed and snap-frozen in dry ice. All samples were stored at − 80 °C until further use.
All plasma, cell and tissue samples were extracted using TRIzol LS (Life Technologies, USA), as described , following the manufacturer’s protocol. Blood cells (250 μL) were homogenised using three volumes of TRIzol LS. Other tissues (each 50 mg) were thawed in 1 mL of TRIzol LS and homogenised with Lysing Matrix D (MP Biomedicals, UK) and a FastPrep FP120 Tissue Disruptor (Thermo Electron, USA). Before extraction, samples of plasma (700 μL) and blood cells (250 μL) were spiked with 3.5 μL of exogenous cel-miR-39-3p (5.6 × 108 copies per sample, Qiagen), and glycogen (40 μg; Sigma-Aldrich, USA) was added in the presence of 1/10 per volume of 5 M ammonium acetate salt (Sigma-Aldrich), as recommended by the manufacturer, to facilitate RNA precipitation. RNA pellets from plasma, blood cells and other tissues were re-suspended, respectively, in 20, 10 and 40 μL of RNase-free water (Qiagen) and used immediately or frozen at − 80 °C until use. RNA content and quality from cell and tissue samples were determined using the Nanodrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, USA).
Small RNA sequencing
Small RNA libraries were prepared from four matched samples of blood plasma and cells using the TruSeq Small RNA Library Preparation Kit (Illumina, USA) and were submitted to 36-base single-end sequencing on the HiSeq 2000 Sequencing System (Illumina). Small RNA libraries were prepared using 5 μL of RNA extract. Raw sequencing data, available on the GEO database (Accession GSE84871 ), were analysed using sRNAbench 1.0 [27, 28]. Briefly, the software was run in genome mode using default settings, with the bovine genome (bosTau4) and miRBase 21 as reference (accessed on 30 April 2015 [29, 30]) in order to identify bovine (bta) miRNAs and human (hsa) miRNA homologues. Human homologues not previously described in bovine were identified by taking reads that aligned to the bovine genome but did not match with known bovine sequences in miRBase and aligning them to human sequences in miRBase. Throughout the manuscript, these miRNAs are indicated with the suffix hsa-. A single nucleotide mismatch was allowed when mapping reads to known miRNA sequences. Before mapping, reads without sequencing adaptor or with undetermined bases, and reads below 15 nucleotides in length (after adaptor removal) were removed and not used for subsequent analyses. For additional details about the integrated analysis steps in sRNAbench please refer to the software manual (http://bioinfo2.ugr.es:8080/ceUGR/srnabench/).
Normalised reads (reads per million mapped, RPMM) obtained from sRNAbench were filtered before being passed to edgeR 3.10.2 for differential expression analysis in R language 3.2.1 [31, 32]. Prior to analysis, miRNAs which were detected with less than 25 RPMM in more than two samples of either plasma or cells were excluded in order to increase confidence in subsequent analysis and reduce technical noise. Differential expression analysis was performed using edgeR’s paired mode on 212 miRNAs (Additional file 1) using GLMfit, as described in . The statistical significance was set to false discovery rate (FDR) < 0.05.
Identification of potential novel biomarkers in cattle
Potential novel tissue biomarkers were identified using the list of miRNAs which were enriched in plasma (Additional file 1). MiRNAs that were either registered for bovine but not human, or that were human homologues not previously identified in cow were selected. To do this, Ensembl genome browser 85 (www.ensembl.org, ), NCBI (www.ncbi.nlm.nih.gov) and miRBase 21 (www.mirbase.org ) were used to identify miRNAs annotated in cow and human, and to determine their conservation between the two species. BLAST was also used (accessed via Ensembl at www.ensembl.org/Multi/Tools/Blast?db=core, ) with default settings to identify the genomic location of selected miRNAs.
MiRNA levels were quantified in matched blood plasma and cell samples. For plasma samples, 2 μL RNA were reverse-transcribed in a 10 μL reaction using the miScript II RT Kit (Qiagen) in a Whatman-Biometra Thermocycler (Biometra, USA). For blood cell and other tissue samples, 500 ng of RNA were used in each 10 μL reaction. The cDNA template was diluted 40-fold and added to 10 μL qPCR reactions which were prepared in 96-well format using the miScript SYBR Green PCR Kit (Qiagen). Template amplification was carried out in an Agilent MX3000P qPCR system (Agilent Technologies, USA). Raw fluorescence data were collected using MxPro software (Agilent Technologies). The amplification efficiency ranged between 84.5% and 117.4%, with R2 > 0.86. Expression levels were determined relative to a freshly-made standard curve and data were analysed using Microsoft Excel (Microsoft Corporation, USA). Expression levels were normalised to the mean expression of spiked-in cel-miR-39-3p for plasma and blood cells. For other tissues, expression data were normalised to RnU6–2.
Statistical analyses were carried out using GraphPad Prism 7 (GraphPad Software, USA). Differences in miRNA expression between blood plasma and cell samples were assessed using paired t-tests. In all cases, normality was tested using the Shapiro-Wilk normality test, outliers were tested with the ROUT test, and data were log2(x + 1) transformed to meet the tests’ normality criteria. Statistical significance was set to P < 0.05.
Results and discussion
Small RNA sequencing of different blood fractions
To determine the relative abundance of miRNAs in different blood compartments, we sequenced miRNAs in paired blood plasma and cell fractions from four cows. Analysis of read length distributions showed a peak at 20–23 nucleotides (characteristic of mature miRNAs) in blood cells and, although of lesser magnitude, in plasma (Fig. .1a). Unlike blood cells, plasma samples displayed an additional peak at 7–9 nucleotides. We attribute the higher proportion of short RNA fragments in plasma relative to cell samples to 1) the higher RNase content in plasma  resulting in higher levels of RNA degradation and 2) the relatively lower abundance of RNAs in plasma resulting in libraries with a higher proportion of adapter/adapter complexes. The difficulty associated with visualising low abundance cDNA bands on a gel during preparation of the plasma library likely resulted in many of the shorter sequences being included for sequencing.
Averaged for both plasma and cells, we obtained 14.9 million raw reads per sample (Table 1). After removing reads with low quality (reads without adapter, short reads and reads with multiple undetermined base calls), we mapped on average 5.9 million and 13.3 million reads to the bovine genome (see Methods) which correspond to 38.1% and 94.3% of the total reads for plasma and cells, respectively (Table 1, Fig. .1b). On average, plasma generated 2.7 million miRNA reads compared to 12.3 million reads for cells, the majority of which were, in both cases, bovine miRNAs (Table 1, Fig. .1b). After applying a detection threshold of 10 reads per sample, we detected up to 315 unique miRNAs in plasma samples and 305 in cell samples, of which only 17 and 16, respectively, were human homologues (Table 1).
Differentially expressed miRNAs in blood plasma and cells
There was incomplete overlap between the lists of most abundant miRNAs in plasma and cells, with only 5 of the 15 most abundant miRNAs being common to the two lists, i.e. miR-486, miR-142-5p, miR-191, miR-92a and miR-30e-5p (Table 2).
Principal component analysis (PCA) of a total of 212 miRNAs that were detected with ≥25 RPMM in more than two samples of either plasma or cells (see Methods, Additional file 1) produced well-separated plasma and cell populations, with a single component (PC1) accounting for 83.9% of the variation across samples (Fig. 2a). Moreover, 169 miRNAs were differentially expressed between plasma and cells (FDR < 0.05, Fig. 2b, Additional file 1), with as many as 131 miRNAs changing by more than two-fold. Of these, 92 miRNAs were enriched in plasma and 39 miRNAs were enriched in blood cells (Table 3, Fig. 3 and Additional file 1). Among miRNAs enriched in plasma were several known to be expressed exclusively or predominantly in specific tissues in humans including miR-122 (liver), miR-133a (muscle), miR-127 (adrenal gland), miR-141 (adrenal gland and reproductive system) and miR-182 (thymus) [37, 38]. Moreover, miRNAs enriched in the bovine blood cell fraction included many previously reported to be blood cell-derived in humans, such as miR-144, miR-451, let-7f, miR-26a, miR-15b, miR-20a, miR-16a, miR-16b and miR-486 [18, 25, 37, 39,40,41].
Validation of sequencing results using RT-qPCR
We used RT-qPCR to validate differences involving a total of 10 miRNAs (indicated in bold in Additional file 1) using samples from a total of six animals (including the four animals used for sequencing). Levels of one of the miRNAs selected for validation (miR-429) could not be accurately quantified in most samples because of its low abundance. For the remaining nine miRNAs (Fig. 4a) we confirmed differences (P < 0.05) in the levels of all but miR-375 between plasma and cell samples (Fig. 4b). We then sought to validate the results for these eight miRNAs (miR-122, miR-215, miR-133a, miR-144, miR-451, and miR-6119-5p, miR-26a and let-7f) using samples from an independent group of animals (four Holstein-Friesian cross cows, aged between 24 and 48 months, in late pregnancy or post-partum) and we confirmed differences in plasma and cell levels of seven miRNAs (Additional file 2); levels of miR-133a were very low (Ct 34–37) in the original group of animals and were undetectable in this second group.
MiRNA profiling across bovine tissues
To identify novel tissue-specific miRNAs in cattle we investigated the distribution across 14 different body tissues of several miRNAs which expression we had found to be enriched in plasma (Additional file 1). We first profiled three miRNAs, miR-122, miR-133a and miR-215, known to be tissue-specific in humans [37, 38]. Accordingly, the three miRNAs (Fig. 5a) were expressed predominantly in liver (miR-122, 88-fold higher than in any other tissue), muscle/heart (miR-133, 254-fold higher) and intestine (miR-215, 150-fold higher).
We then sought to identify novel tissue biomarkers. First, from the list of miRNAs enriched in plasma (Additional file 1), we looked for miRNAs that were either registered for bovine but not human, or that were human homologues not previously identified in cow, as described in the Materials and Methods. These analyses revealed 10 miRNAs which either were present in bovine but not in human (miR-1388-5p, miR-1388-3p, miR-2376, miR-2285 t, miR-3431 and miR-6529a, Fig. 5b) or were homologues of human miRNAs not identified previously in cow (miR-210-3p, miR-802, miR-4532 and miR-4792, Fig. 5c).
RT-qPCR analyses showed that the majority of these miRNAs were expressed in most of the tissues profiled, i.e., they were not enriched in any particular tissue and therefore do not provide potential biomarkers of tissue function (Fig. 5b-c). Surprisingly, miR-1388-3p, miR-2376 and miR-6529 showed high expression in blood cells despite having been identified as plasma-enriched using sequencing. A possible explanation for this result is that, despite blood cells having the highest levels of these miRNAs, the combined contribution from all other tissues to plasma levels is still greater in comparison. Interestingly, in addition to blood cells, all three miRNAs were relatively abundant in spleen (Fig. 5b), which suggests an involvement in hematopoietic regulation. However, since no functional data is available for any of the three miRNAs in any species, further study will be needed to confirm this.
Among the remaining miRNAs, miR-3431 was distinctly high in lung followed by placenta and uterus. This is a ruminant-specific miRNA which has been reported to be expressed in adipose and reproductive tissues of cattle [42, 43], although no information is available on its potential functions. Whether miR-3431 could provide a useful circulating biomarker of lung function in cows should be investigated in future studies.
Of note, the human homologue of miR-802 was distinctly enriched in liver, where its expression was 87-fold higher than the mean expression across all tissues, and 5.2-fold higher than in intestine, the tissue with the second highest mean expression levels (Fig. 5c). Using BLAT (UMD3.1, accessed 27/07/16) we mapped the mature miR-802 sequence to bovine chromosome 1 (chr1: 149720302–149,720,392 [+]). Indeed, this location in the bovine genome has been annotated as a novel miRNA with a secondary structure that matches miR-802. Although not previously reported in cow, data in humans and mice have shown miR-802 to be a key a regulator of liver function. Specifically, hepatic levels of miR-802 were elevated in obese individuals, contributing to insulin insensitivity, glucose intolerance and increased hepatic gluconeogenesis by downregulation of HNF1B . Interestingly, serum miR-802 levels were reportedly increased in type 2 diabetes patients providing a potential circulating biomarker for this disease . Another study identified miR-802 as a target of constitutive androstane receptor (CAR) signalling, which is involved in the regulation of xenobiotic detoxification, lipid homeostasis and energy metabolism . Furthermore, circulating miR-802 has been proposed as a biomarker of drug-induced liver damage in rats, with an observed increase in miRNA levels in response to tissue damage being comparable to that of the liver-enriched miR-122 [47, 48]. Overall, these results highlight the importance of miR-802 in the regulation of glucose and lipid metabolism as well as xenobiotic responses in the liver.
In dairy cattle, negative energy balance (NEB) often occurs during the post-partum period as a consequence of the high energy requirements for lactation. The liver has a central role in counteracting NEB. Liver dysfunction often ensues during the post-partum period, particularly in high-producing cows, in association with metabolic dysfunction and low fertility, constituting a major problem in modern dairy herds . Changes in the expression of some miRNAs in liver have been shown in association with NEB in cattle . In this context, future studies should investigate whether miR-802, perhaps in combination with liver-enriched miR-122, may provide a potential biomarker of negative energy balance in dairy cows.
Using next-generation sequencing, we compared for the first time miRNA levels in paired plasma and cell fractions of bovine blood and identified a total of 131 miRNAs that were differentially expressed (greater than two-fold) between the two fractions. Of those, 92 miRNAs were expressed at higher levels in plasma than in cells. Expression profiling of selected plasma-enriched miRNAs across different bovine tissues revealed miR-802, a previously unidentified bovine miRNA, to be highly enriched in liver. The miRNA database generated in this study provides a useful resource for future investigation of miRNA biology in bovine and identification of additional potential biomarkers of tissue function in livestock.
Principal component analysis
Quantitative polymerase chain reaction
Hammond SM. An overview of microRNAs. Adv Drug Deliv Rev. 2015;87:3–14.
Lee RC, Feinbaum RL, Ambros V. The C. Elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993;75:843–54.
Casey MC, Kerin MJ, Brown JA, Sweeney KJ. Evolution of a research field-a micro (RNA) example. PeerJ. 2015;3:e829.
Motti D, Bixby JL, Lemmon VP. MicroRNAs and neuronal development. Semin Fetal Neonatal Med. 2012;17:347–52.
Finch ML, Marquardt JU, Yeoh GC, Callus BA. Regulation of microRNAs and their role in liver development, regeneration and disease. Int J Biochem Cell Biol. 2014;54:288–303.
Vienberg S, Geiger J, Madsen S, Dalgaard LT. MicroRNAs in metabolism. Acta Physiol (Oxf). 2017;219(2):346–61.
Sontakke SD, Mohammed BT, McNeilly AS, Donadeu FX. Characterization of microRNAs differentially expressed during bovine follicle development. Reproduction. 2014;148:271–83.
Farazi TA, Hoell JI, Morozov P, Tuschl T. MicroRNAs in human cancer. Adv Exp Med Biol. 2013;774:1–20.
Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, Peterson A, Noteboom J, O'Briant KC, Allen A, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A. 2008;105:10513–8.
He Y, Lin J, Kong D, Huang M, Xu C, Kim TK, Etheridge A, Luo Y, Ding Y, Wang K. Current state of circulating MicroRNAs as Cancer biomarkers. Clin Chem. 2015;61:1138–55.
Schwarzenbach H, Nishida N, Calin GA, Pantel K. Clinical relevance of circulating cell-free microRNAs in cancer. Nat Rev Clin Oncol. 2014;11:145–56.
Creemers EE, Tijsen AJ, Pinto YM. Circulating microRNAs: novel biomarkers and extracellular communicators in cardiovascular disease? Circ Res. 2012;110:483–95.
Afonso MB, Rodrigues PM, Simao AL, Castro RE. Circulating microRNAs as potential biomarkers in non-alcoholic fatty liver disease and hepatocellular carcinoma. J Clin Med. 2016;5(3):30.
Guay C, Regazzi R. Circulating microRNAs as novel biomarkers for diabetes mellitus. Nat Rev Endocrinol. 2013;9:513–21.
Miura K, Miura S, Yamasaki K, Higashijima A, Kinoshita A, Yoshiura K, Masuzaki H. Identification of pregnancy-associated microRNAs in maternal plasma. Clin Chem. 2010;56:1767–71.
Pritchard CC, Kroh E, Wood B, Arroyo JD, Dougherty KJ, Miyaji MM, Tait JF, Tewari M. Blood cell origin of circulating microRNAs: a cautionary note for cancer biomarker studies. Cancer Prev Res (Phila). 2012;5:492–7.
Leidner RS, Li L, Thompson CL. Dampening enthusiasm for circulating microRNA in breast cancer. PLoS One. 2013;8:e57841.
Haider BA, Baras AS, McCall MN, Hertel JA, Cornish TC, Halushka MK. A critical evaluation of microRNA biomarkers in non-neoplastic disease. PLoS One. 2014;9:e89565.
Ioannidis J, Donadeu FX. Circulating miRNA signatures of early pregnancy in cattle. BMC Genomics. 2016;17:184.
Ioannidis J, Donadeu FX. Circulating microRNA profiles during the bovine Oestrous cycle. PLoS One. 2016;11:e0158160.
Farrell D, Shaughnessy RG, Britton L, MacHugh DE, Markey B, Gordon SV. The identification of circulating MiRNA in bovine serum and their potential as novel biomarkers of early Mycobacterium avium subsp paratuberculosis infection. PLoS One. 2015;10:e0134310.
Hansen EP, Kringel H, Thamsborg SM, Jex A, Nejsum P. Profiling circulating miRNAs in serum from pigs infected with the porcine whipworm, Trichuris suis. Vet Parasitol. 2016;223:30–3.
Muroya S, Ogasawara H, Hojito M. Grazing affects Exosomal circulating MicroRNAs in cattle. PLoS One. 2015;10:e0136475.
Ahanda ML, Zerjal T, Dhorne-Pollet S, Rau A, Cooksey A, Giuffra E. Impact of the genetic background on the composition of the chicken plasma MiRNome in response to a stress. PLoS One. 2014;9:e114598.
Cheng HH, Yi HS, Kim Y, Kroh EM, Chien JW, Eaton KD, Goodman MT, Tait JF, Tewari M, Pritchard CC. Plasma processing conditions substantially influence circulating microRNA biomarker levels. PLoS One. 2013;8:e64795.
Ioannidis J, Donadeu FX. MicroRNA biomarkers of tissue function. GEO. 2016. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE84871.
Barturen G, Rueda A, Hamberg M, Alganza A, Lebron R, Kotsyfakis M, Shi BJ, Koppers-Lalic D, Hackenberg M. sRNAbench: profiling of small RNAs and its sequence variants in single or multi-species high-throughput experiments. Methods in Next Generation Sequencing. 2014; https://doi.org/10.2478/mngs-2014-0001.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, et al. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res. 2013;41:D991–5.
Elsik CG, Unni DR, Diesh CM, Tayal A, Emery ML, Nguyen HN, Hagen DE. Bovine genome database: new tools for gleaning function from the Bos taurus genome. Nucleic Acids Res. 2016;44:D834–9.
Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ. miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008;36:D154–8.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2013.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40:4288–97.
Aken BL, Ayling S, Barrell D, Clarke L, Curwen V, Fairley S, Fernandez Banet J, Billis K, Garcia Giron C, Hourlier T, et al. The Ensembl gene annotation system. Database (Oxford). 2016; https://doi.org/10.1093/database/baw093.
Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011;39:D152–7.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinf. 2009;10:421.
Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007;129:1401–14.
Guo Z, Maki M, Ding R, Yang Y, Zhang B, Xiong L. Genome-wide survey of tissue-specific microRNA and transcription factor regulatory networks in 12 tissues. Sci Rep. 2014;4:5150.
Wang K, Yuan Y, Cho J-H, McClarty S, Baxter D, Galas DJ. Comparing the MicroRNA Spectrum between serum and plasma. PLoS One. 2012;7:e41561.
Hunter MP, Ismail N, Zhang X, Aguda BD, Lee EJ, Yu L, Xiao T, Schafer J, Lee ML, Schmittgen TD, et al. Detection of microRNA expression in human peripheral blood microvesicles. PLoS One. 2008;3:e3694.
Willeit P, Zampetaki A, Dudek K, Kaudewitz D, King A, Kirkby NS, Crosby-Nwaobi R, Prokopi M, Drozdov I, Langley SR, et al. Circulating MicroRNAs as novel biomarkers for platelet activation. Circ Res. 2013;112:595-+.
Jin W, Dodson MV, Moore SS, Basarab JA, Guan LL. Characterization of microRNA expression in bovine adipose tissues: a potential regulatory mechanism of subcutaneous adipose tissue development. BMC Mol Biol. 2010;11:29.
Huang J, Ju Z, Li Q, Hou Q, Wang C, Li J, Li R, Wang L, Sun T, Hang S, et al. Solexa sequencing of novel and differentially expressed microRNAs in testicular and ovarian tissues in Holstein cattle. Int J Biol Sci. 2011;7:1016–26.
Kornfeld JW, Baitzel C, Konner AC, Nicholls HT, Vogt MC, Herrmanns K, Scheja L, Haumaitre C, Wolf AM, Knippschild U, et al. Obesity-induced overexpression of miR-802 impairs glucose metabolism through silencing of Hnf1b. Nature. 2013;494:111–5.
Higuchi C, Nakatsuka A, Eguchi J, Teshigawara S, Kanzaki M, Katayama A, Yamaguchi S, Takahashi N, Murakami K, Ogawa D, et al. Identification of circulating miR-101, miR-375 and miR-802 as biomarkers for type 2 diabetes. Metabolism. 2015;64:489–97.
Hao R, Su S, Wan Y, Shen F, Niu B, Coslo DM, Albert I, Han X, Omiecinski CJ. Bioinformatic analysis of microRNA networks following the activation of the constitutive androstane receptor (CAR) in mouse liver. Biochim Biophys Acta. 1859;2016:1228–37.
Church RJ, Otieno M, McDuffie JE, Singh B, Sonee M, Hall L, Watkins PB, Ellinger-Ziegelbauer H, Harrill AH. Beyond miR-122: identification of MicroRNA alterations in blood during a time course of hepatobiliary injury and biliary hyperplasia in rats. Toxicol Sci. 2016;150:3–14.
Wolenski FS, Shah P, Sano T, Shinozawa T, Bernard H, Gallacher MJ, Wyllie SD, Varrone G, Cicia LA, Carsillo ME, et al. Identification of microRNA biomarker candidates in urine and plasma from rats with kidney or liver damage. J Appl Toxicol. 2016;37(3):278–86.
Fatima A, Lynn DJ, O'Boyle P, Seoighe C, Morris D. The miRNAome of the postpartum dairy cow liver in negative energy balance. BMC Genomics. 2014;15:279.
We thank Judith Risse, Helen Gunter and Mick Watson for assistance in the interpretation of sequencing data. We also thank Cheryl Ashworth for providing valuable feedback throughout the study and Helen Brown for helping with statistical analyses. Finally, we are very thankful to Tim King and Wilson Lee for advice in planning the animal study and carrying out the animal work, and to Bushra Mohammed and Sadanand Sontakke for assisting with the collection of tissue samples.
This study was funded by Zoetis, Inc. and the Biotechnology and Biological Sciences Research Council (BBSRC). The funding bodies did not participate in study design, data collection, analysis, interpretation or writing of the manuscript.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the GEO repository, Accession GSE84871, www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=olgtysuirzqfpmr&acc=GSE84871 .
All animal procedures were carried out under the UK Home Office Animals (Scientific Procedures) Act 1986, license 60/4604, with approval by the Ethical Review Committee, University of Edinburgh.
This study was partially financed by Zoetis, Inc.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sequencing miRNA analysis. Normalised miRNA expression levels in bovine blood plasma and cell samples, and the results of differential expression analysis between the two sample types. (XLSX 47 kb)
Validation of sequencing data. Validation of sequencing results for selected miRNAs using RT-qPCR on an independent group of four animals. (TIFF 1669 kb)