Intra-individual variation of miRNA expression levels in human plasma samples

Background: Circulating miRNAs as potential non-invasive biomarkers for disease risk assessment and cancer early diagnosis have attracted increasing interest. Little information, however, is available regarding the intra-individual variation of circulating miRNA levels.

Methods: We measured expression levels of a panel of 800 miRNAs in repeated plasma samples from 51 healthy individuals that were collected 6 to 12?months apart and evaluated the intra-individual variation by the intra-class correlation coefficient (ICC).

Results: After background correction, a total of 185 miRNAs were detected in at least 10% of the plasma samples, with 69 and 28 miRNAs being detected in 50% and 90% of samples, respectively. The median ICC was 0.46 for these 185 miRNAs. Among them, 41% (75 miRNAs) had an ICC?≥?0.5, and 23% (42 miRNAs) had an ICC?≥?0.6. The ICC is higher for miRNAs with higher expression levels or higher detection rates, when compared to those with lower expression levels or lower detection rates.

Conclusions: These results suggest that common circulating miRNAs are stable over a relatively long period and can serve as reliable biomarkers for epidemiological and clinical research.     

Crosstalk between circulating microRNAs and chronotypical features in subjects with metabolic syndrome

ABSTRACT

Circulating microRNAs (miRNAs) are valuable biomarkers that may provide important insight into the pathogenesis of metabolic syndrome (MetS). Moreover, there is an association between chronotypical characteristics and MetS predisposition. Considering that expression of some miRNAs is circadian-rhythm-dependent, the aim of this study was to investigate the circulating miRNA profile in subjects with and without MetS in association with chronotype. The expression of 86 metabolic syndrome-related miRNAs was investigated in the plasma of 21 subjects with MetS and in 82 subjects without MetS using miRCURY LNA miRNA PCR System technology. Chronotype was assessed using the Horne and Östberg Morningness-Eveningness Questionnaire. Bioinformatic analyses were performed to explore the target genes and biological pathways regulated by the selected miRNAs. Subjects with MetS were more often evening chronotype compared to non-MetS controls. Additionally, four miRNAs (miR-140-3p, miR-150-5p, miR-375, and miR-29 c-3p) demonstrated interaction with MetS and chronotype. Interestingly, the target genes of these four miRNAs participate in pathways related to the circadian clock. In conclusion, we identified four circulating miRNAs whose circulating levels could interact with MetS and chronotype.     

A pilot study of circulating miRNAs as potential biomarkers of early stage breast cancer

Background

To date, there are no highly sensitive and specific minimally invasive biomarkers for detection of breast cancer at an early stage. The occurrence of circulating microRNAs (miRNAs) in blood components (including serum and plasma) has been repeatedly observed in cancer patients as well as healthy controls. Because of the significance of miRNA in carcinogenesis, circulating miRNAs in blood may be unique biomarkers for early and minimally invasive diagnosis of human cancers. The objective of this pilot study was to discover a panel of circulating miRNAs as potential novel breast cancer biomarkers.

Methodology/Principal Findings

Using microarray-based expression profiling followed by Real-Time quantitative Polymerase Cycle Reaction (RT-qPCR) validation, we compared the levels of circulating miRNAs in plasma samples from 20 women with early stage breast cancer (10 Caucasian American (CA) and 10 African American (AA)) and 20 matched healthy controls (10 CAs and 10 AAs). Using the significance level of p<0.05 constrained by at least two-fold expression change as selection criteria, we found that 31 miRNAs were differentially expressed in CA study subjects (17 up and 14 down) and 18 miRNAs were differentially expressed in AA study subjects (9 up and 9 down). Interestingly, only 2 differentially expressed miRNAs overlapped between CA and AA study subjects. Using receiver operational curve (ROC) analysis, we show that not only up-regulated but also down-regulated miRNAs can discriminate patients with breast cancer from healthy controls with reasonable sensitivity and specificity. To further explore the potential roles of these circulating miRNAs in breast carcinogenesis, we applied pathway-based bioinformatics exploratory analysis and predicted a number of significantly enriched pathways which are predicted to be regulated by these circulating miRNAs, most of which are involved in critical cell functions, cancer development and progression.

Conclusions

Our observations from this pilot study suggest that the altered levels of circulating miRNAs might have great potential to serve as novel, noninvasive biomarkers for early detection of breast cancer.     

Influence of Plasma Processing on Recovery and Analysis of Circulating Nucleic Acids

Circulating nucleic acids (CNAs) are under investigation as a liquid biopsy in cancer. However there is wide variation in blood processing and methods for isolation of circulating free DNA (cfDNA) and microRNAs (miRNAs). Here we compare the extraction efficiency and reproducibility of 4 commercially available kits for cfDNA and 3 for miRNA using spike-in of reference templates. We also compare the effects of increasing time between venepuncture and centrifugation and differential centrifugation force on recovery of CNAs. cfDNA was quantified by TaqMan qPCR and targeted deep sequencing. miRNA profiles were assessed with TaqMan low-density arrays and assays. The QIAamp^® DNA Blood Mini and Circulating nucleic acid kits gave the highest recovery of cfDNA and efficient recovery (>90%) of a 564bp spike-in. Moreover, targeted sequencing revealed overlapping cfDNA profiles and variant depth, including detection of HER2 gene amplification, using the Ion AmpliSeq™Cancer Hotspot Panel v2. Highest yields of miRNA and the synthetic Arabidopsis thaliana miR-159a spike-in were obtained using the miRNeasy Serum/Plasma kit, with saturation above 200 µl of plasma. miRNA profiles showed significant variation with increasing time before centrifugation (p<0.001) and increasing centrifugation force, with depletion of platelet associated miRNAs, whereas cfDNA was unaffected. However, sample replicates showed excellent reproducibility on TaqMan low density arrays (ρ = 0.96, p<0.0001). We also successfully generated miRNA profiles for plasma samples stored > 12 years, highlighting the potential for analysis of stored sample biobanks. In the era of the liquid biopsy, standardisation of methods is required to minimise variation, particularly for miRNA.     

Pre-storage centrifugation conditions have significant impact on measured microRNA levels in biobanked EDTA plasma samples

BackgroundIn the past few years, an increasing number of studies have reported the potential use of microRNAs (miRNA) as circulating biomarkers for diagnosis or prognosis of a wide variety of diseases. There is, however, a lack of reproducibility between studies. Due to the high miRNA content in platelets this may partly be explained by residual platelets in the plasma samples used. When collecting fresh plasma samples, it is possible to produce cell-free/platelet-poor plasma by centrifugation. In this study, we systematically investigated whether biobanked EDTA plasma samples could be processed to be suitable for miRNA analysis.Materials and methodsBlood samples were collected from ten healthy volunteers and centrifuged to produce platelet-poor-plasma (PPP) and standard biobank plasma. After one week at ?80 °C the biobanked EDTA plasma was re-centrifuged by different steps to remove residual platelets. Using RT-qPCR the levels of 14 miRNAs in the different plasma preparations were compared to that of PPP.ResultsWe were able to remove residual platelets from biobanked EDTA plasma by re-centrifugation of the thawed samples. Nevertheless, for most of the investigated miRNAs, the miRNA level was significantly higher in the re-centrifuged biobanked plasma compared to PPP, even when the platelet count was reduced to 0–1×10⁹/L.ConclusionWe found, that pre-storage centrifugation conditions have a significant impact on the measured EDTA plasma level of miRNAs known to be present in platelets. Even for the miRNAs found to be less effected, we showed that a 1.5–3 fold change in plasma levels may possible be caused by or easily overseen due to sample preparation and/or storage.     

Identification of a Plasma Four-microRNA Panel as Potential Noninvasive Biomarker for Osteosarcoma

Background

Circulating microRNAs (miRNAs) are emerging as promising biomarkers for human cancer. Osteosarcoma is the most common human primary malignant bone tumor in children and young adults. The objective of this study was to investigate whether circulating miRNAs in plasma could be a useful biomarker for detecting osteosarcoma and monitoring tumor removal dynamics.

Methods

Plasma samples were obtained from 90 patients before surgery, 50 patients after one month of surgery, and 90 healthy individuals. The study was divided into three steps: First, initial screening of the profiles of circulating miRNAs in pooled plasma samples from healthy controls and pre-operative osteosarcoma patients using a TaqMan low density array (TLDA). Second, evaluation of miRNA concentration in individual plasma samples from 90 pre-operative osteosarcoma patients and 90 healthy controls by a quantitative real time PCR (qRT-PCR) assay. Third, evaluation of miRNA concentration in paired plasma samples from 50 pre- and post-operative osteosarcoma patients by qRT-PCR assay.

Results

Four plasma miRNAs including miR-195-5p, miR-199a-3p, miR-320a, and miR-374a-5p were significantly increased in the osteosarcoma patients. Receiver operating characteristics curve analysis of the combined populations demonstrated that the four-miRNA signature could discriminate cases from controls with an area under the curve of 0.9608 (95% CI 0.9307-0.9912). These 4 miRNAs were markedly decreased in the plasma after operation. In addition, circulating miR-195-5p and miR-199a-3p were correlated with metastasis status, while miR-199a-3p and miR-320a were correlated with histological subtype.

Conclusions

Our data suggest that altered levels of circulating miRNAs might have great potential to serve as novel, non-invasive biomarkers for osteosarcoma.     

Plasma Processing Conditions Substantially Influence Circulating microRNA Biomarker Levels

Circulating, cell-free microRNAs (miRNAs) are promising candidate biomarkers, but optimal conditions for processing blood specimens for miRNA measurement remain to be established. Our previous work showed that the majority of plasma miRNAs are likely blood cell-derived. In the course of profiling lung cancer cases versus healthy controls, we observed a broad increase in circulating miRNA levels in cases compared to controls and that higher miRNA expression correlated with higher platelet and particle counts. We therefore hypothesized that the quantity of residual platelets and microparticles remaining after plasma processing might impact miRNA measurements. To systematically investigate this, we subjected matched plasma from healthy individuals to stepwise processing with differential centrifugation and 0.22 µm filtration and performed miRNA profiling. We found a major effect on circulating miRNAs, with the majority (72%) of detectable miRNAs substantially affected by processing alone. Specifically, 10% of miRNAs showed 4–30x variation, 46% showed 30-1,000x variation, and 15% showed >1,000x variation in expression solely from processing. This was predominantly due to platelet contamination, which persisted despite using standard laboratory protocols. Importantly, we show that platelet contamination in archived samples could largely be eliminated by additional centrifugation, even in frozen samples stored for six years. To minimize confounding effects in microRNA biomarker studies, additional steps to limit platelet contamination for circulating miRNA biomarker studies are necessary. We provide specific practical recommendations to help minimize confounding variation attributable to plasma processing and platelet contamination.     

Circulating microRNAs as Specific Biomarkers for Breast Cancer Detection

Background

We previously showed microRNAs (miRNAs) in plasma are potential biomarkers for colorectal cancer detection. Here, we aimed to develop specific blood-based miRNA assay for breast cancer detection.

Methodology/Principal Findings

TaqMan-based miRNA profiling was performed in tumor, adjacent non-tumor, corresponding plasma from breast cancer patients, and plasma from matched healthy controls. All putative markers identified were verified in a training set of breast cancer patients. Selected markers were validated in a case-control cohort of 170 breast cancer patients, 100 controls, and 95 other types of cancers and then blindly validated in an independent set of 70 breast cancer patients and 50 healthy controls. Profiling results showed 8 miRNAs were concordantly up-regulated and 1 miRNA was concordantly down-regulated in both plasma and tumor tissue of breast cancer patients. Of the 8 up-regulated miRNAs, only 3 were significantly elevated (p<0.0001) before surgery and reduced after surgery in the training set. Results from the validation cohort showed that a combination of miR-145 and miR-451 was the best biomarker (p<0.0001) in discriminating breast cancer from healthy controls and all other types of cancers. In the blind validation, these plasma markers yielded Receiver Operating Characteristic (ROC) curve area of 0.931. The positive predictive value was 88% and the negative predictive value was 92%. Altered levels of these miRNAs in plasma have been detected not only in advanced stages but also early stages of tumors. The positive predictive value for ductal carcinoma in situ (DCIS) cases was 96%.

Conclusions

These results suggested that these circulating miRNAs could be a potential specific biomarker for breast cancer screening.     

Detection of microRNA Expression in Human Peripheral Blood Microvesicles

Background

MicroRNAs (miRNA) are small non-coding RNAs that regulate translation of mRNA and protein. Loss or enhanced expression of miRNAs is associated with several diseases, including cancer. However, the identification of circulating miRNA in healthy donors is not well characterized. Microvesicles, also known as exosomes or microparticles, circulate in the peripheral blood and can stimulate cellular signaling. In this study, we hypothesized that under normal healthy conditions, microvesicles contain miRNAs, contributing to biological homeostasis.

Methodology/Principal Findings

Microvesicles were isolated from the plasma of normal healthy individuals. RNA was isolated from both the microvesicles and matched mononuclear cells and profiled for 420 known mature miRNAs by real-time PCR. Hierarchical clustering of the data sets indicated significant differences in miRNA expression between peripheral blood mononuclear cells (PBMC) and plasma microvesicles. We observed 71 miRNAs co-expressed between microvesicles and PBMC. Notably, we found 33 and 4 significantly differentially expressed miRNAs in the plasma microvesicles and mononuclear cells, respectively. Prediction of the gene targets and associated biological pathways regulated by the detected miRNAs was performed. The majority of the miRNAs expressed in the microvesicles from the blood were predicted to regulate cellular differentiation of blood cells and metabolic pathways. Interestingly, a select few miRNAs were also predicted to be important modulators of immune function.

Conclusions

This study is the first to identify and define miRNA expression in circulating plasma microvesicles of normal subjects. The data generated from this study provides a basis for future studies to determine the predictive role of peripheral blood miRNA signatures in human disease and will enable the definition of the biological processes regulated by these miRNA.     

Identification of Suitable Reference Genes for qRT-PCR Analysis of Circulating microRNAs in Hepatitis B Virus-Infected Patients

Circulating microRNAs (miRNAs) were found to exist in serum/plasma in a highly stable, cell-free form, and aberrantly expressed in many human diseases. Currently, the expression levels of circulating miRNAs are estimated by quantitative real-time polymerase chain reaction. However, no study has systematically evaluated reference genes for evaluating circulating microRNA expression. This study describes the identification and characterization of an appropriate reference gene for the normalization of circulating miRNA levels in hepatitis B virus (HBV)-infected patients and healthy people. Ten miRNAs that resemble the mean expression of the TaqMan low density array together with U6, RNU6B, and miR-16 were validated with two algorithms, geNorm, and NormFinder, after ensuring their equivalent expression between the two study groups. The combination of miR-26a, miR-221, and miR-22* is recommended as the most stable set of reference genes for circulating miRNA evaluation in HBV patients and healthy people.     

Circulating miRNA-33: a potential biomarker in patients with coronary artery disease

Background: Circulating microRNAs (miRNA) are present in body fluids in stable, cell-free form. Likewise, these miRNAs can be identified in various stages of coronary artery disease (CAD) such as inflammation, endothelial dysfunction, proliferation and atherosclerosis among others. miRNA expression levels can be identified.

Aims and objectives: To determine the expression of circulating miRNAs (miR-126, miR-92, miR-33, miR-145 and miR-155) in CAD patients of Indian origin.

Material and methods: miRNA profiling analysis in blood plasma was performed by quantitative real-time-PCR (qRT-PCR) in 60 angiographically verified subjects including 30 CAD patients and 30 age- and gender-matched controls. Association between the expression of all five circulating miRNAs and clinical characteristics of patients with CAD were analysed using Medcalc statistics. The severity of CAD was assessed using SYNTAX score (SS).

Results: Expression of plasma miR-33 increased by 2.9 folds in CAD patients than in control group (p value ≥0.002) also it was found that miR-33 expression levels in mild cases (SS: ≤22) were significantly higher than CAD controls. There was a modest negative correlation between miR-33 and total cholesterol/high density lipoprotein ratio, triglycerides and very low density lipoprotein.

Conclusion: The study reports a significant association between increased levels of plasma miR-33 and CAD. Thus, plasma miR-33 appears to be a promising non-invasive biomarker, but requires further validation in a large cohort.     

Plasma microRNA profiling highlights miR-1260b and miR-720 as novel diagnostic and prognostic biomarkers of esophageal squamous cell carcinoma

This paper was designed to explore the value of miRNAs as diagnostic biomarkers that may facilitate the early detection of esophageal squamous cell carcinoma (ESCC). Plasma miRNA profiles were defined via an array-based approach using samples from ESCC patients and healthy controls (n=5 each). Differentially expressed miRNAs in these samples were validated via qPCR in ESCC patients (n=96) and healthy controls (n=51), and the relationship between ESCC patient plasma miR-1260b and miR-720 levels and clinicopathological characteristics were additionally examined. In total, 12 plasma miRNAs that were differentially expressed between ESCC patients and healthy controls were identified via miRNA. Six of these miRNAs were subsequently validated, revealing that both miR-1260b and miR-720 were significantly differentially abundant in ESCC patients and controls, with miR-1260b being significantly upregulated in ESCC patients relative to controls (2.24, 1.41 respectively, P<0.001), while the opposite was observed with respect to miR-720 (0.66, 2.27 respectively, P=0.001). The use of both miR-720 and miR-1260b as a combined diagnostic tool was highly efficacious, yielding an AUC of 0.814, a sensitivity of 86.3%, and a specificity of 73.2% as a means of detecting ESCC patients. Elevated plasma miR-1260b level was also associated with a poorer patient prognosis when compared to patients with a low plasma miRNA level (P=0.021). This study has successfully developed a plasma miRNA biomarker signature of ESCC that may offer value as a diagnostic or prognostic tool when evaluating patients with ESCC.     

U6 is not a suitable endogenous control for the quantification of circulating microRNAs

Recently, microRNAs have been detected in serum and plasma, and circulating microRNA (miRNA) profiles have now been associated with many diseases such as cancers and heart disease, as well as altered physiological states. Because of their stability and disease resistance, circulation miRNAs appear to be an ideal material for biomarkers of diseases and physiological states in blood. However, the lack of a suitable internal reference gene (internal reference miRNA) has hampered research and application of circulating miRNAs. Currently, U6 and miR-16 are the most common endogenous controls in the research of miRNAs in tissues and cells. We performed microarray-based serum miRNA profiling on the serum of 20 nasopharyngeal carcinoma patients and 20 controls to detect the expressions of U6 and miRNAs. Profiling was followed by real-time quantitative Polymerase Chain Reaction (qPCR) in 80 patients (20 each with gastric cancer, nasopharyngeal carcinoma, colorectal cancer, and breast cancer) and 30 non-cancerous controls. qPCR was also performed to detect miRNAs in serum with repeated freezing and thawing. The results of microarray showed that with the exception of U6, Ct values of miR-16, miR-24, miR-142-3p, miR-19b and miR-192 in serum samples of nasopharyngeal carcinoma were greater than control samples. The results of 110 cases showed large fluctuations in U6 expression. The difference between the greatest and the least levels of expression was 3.29 for delta Ct values, and 1.23 for miR-16. The expressions of U6, miR-16 and miR-24 in serum subjected to different freeze–thaw cycles showed that U6 expression gradually decreased after 1, 2, and 4 cycles of freezing and thawing, while the expression of miR-16 and miR-24 remained relatively stable. Collectively, our results suggested that U6 is unsuitable as an internal reference gene in the research of circulating miRNAs.     

A panel of three plasma microRNAs for colorectal cancer diagnosis

BackgroundDifferential microRNA (miRNA) expression profiles in plasma or serum were identified, providing foundation for studying their potentially diagnostic role in colorectal cancer (CRC).MethodsWe performed S-poly(T) Plus PCR assay to select and validate differentially expressed plasma miRNAs from a sample set including 101 CRC patients, 20 patients with colorectal noncancerous polyps (NCP), and 134 healthy controls. And bioinformatics methods was used to integrated predicted or validated targets of the differentially dysregulated miRNAs and analyzed their overrepresented pathways.ResultsAfter the two-phase selection and validation process, we identified a miRNA panel (miR-144-3p, miR-425-5p, and miR-1260b) with high diagnostic efficiency for CRC; the panel distinguished CRC patients from controls with 93.8% sensitivity and 91.3% specificity. Results indicated that the dysregulated miRNAs in CRC were functionally involved in several key cancer-related pathways, such as axonal guidance, PI3K, and calcium signaling pathways.ConclusionsOur study demonstrated that a plasma 3-miRNA panel may serve as a novel noninvasive biomarker to diagnose CRC. This plasma 3-miRNA panel may be related to CRC development. However, further studies are needed to highlight its theoretical strengths.     

Alterations in The Plasma Expression of mir-15b,mir-195 and the Tumor-Suppressor Gene DLEU7 in Patients with B-Cell Chronic Lymphocytic Leukemia

Background:Chronic lymphocytic leukemia (CLL) is one of the most prevalent forms of leukemia in adults. Inactivation of the DLEU7 gene is frequently observed in patients with CLL. Furthermore, microRNAs (miRNAs) have been observed to have a critical role in the pathogenesis of several cancers, including leukemia. Considering the tumor-suppressive role of DLEU7, as well as the tumor suppressor or oncogenic role of microRNAs (miRNAs), the aim of the present study was to evaluate the potential miRNAs targeting the DLEU7 gene in B-cells and explore expression changes these genes in the plasma of B-CLL patients. Methods:The miRNAs interacting with the DLEU7 gene were predicted and selected using bioinformatics tools. A total of 80 plasma samples were collected from 40 patients with B-cells and 40 healthy individuals, then subjected to RNA extraction and cDNA synthesis. The expression profiles of the predicted miRNAs and the DLEU7 gene in the plasma of B-CLL patients and healthy individuals were determined by RT-qPCR analysis. Results:The bioinformatics prediction indicated that miR-15b and miR-195 target the DLEU7 gene. The expression levels of miR-15b and miR-195 were significantly higher in the plasma of patients with B-CLL compared to the healthy individuals (91.6, p= 0.001) (169, p= 0.001). However, the expression level of the DLEU7 gene was found to be significantly lower in the patient group compared to healthy controls (0.304, p= 0.001).Conclusion:Both miR-15b and miR-195, have the potential to function as novel and non-invasive biomarkers in the diagnosis and prognosis of patients with B-CLL.Key Words: B-CLL, miRNA, Biomarker, DLEU7, RT-QPCR     

Circulating microRNAs as Biomarkers for Detection of Autologous Blood Transfusion

MicroRNAs (miRNAs) are small non-coding RNAs that regulate various biological processes. Cell-free miRNAs measured in blood plasma have emerged as specific and sensitive markers of physiological processes and disease. In this study, we investigated whether circulating miRNAs can serve as biomarkers for the detection of autologous blood transfusion, a major doping technique that is still undetectable. Plasma miRNA levels were analyzed using high-throughput quantitative real-time PCR. Plasma samples were obtained before and at several time points after autologous blood transfusion (blood bag storage time 42 days) in 10 healthy subjects and 10 controls without transfusion. Other serum markers of erythropoiesis were determined in the same samples. Our results revealed a distinct change in the pattern of circulating miRNAs. Ten miRNAs were upregulated in transfusion samples compared with control samples. Among these, miR-30b, miR-30c, and miR-26b increased significantly and showed a 3.9-, 4.0-, and 3.0-fold change, respectively. The origin of these miRNAs was related to pulmonary and liver tissues. Erythropoietin (EPO) concentration decreased after blood reinfusion. A combination of miRNAs and EPO measurement in a mathematical model enhanced the efficiency of autologous transfusion detection through miRNA analysis. Therefore, our results lay the foundation for the development of miRNAs as novel blood-based biomarkers to detect autologous transfusion.     

Identification of endogenous microRNA references in porcine serum for quantitative real-time PCR normalization

MicroRNAs (miRNAs) are evolutionarily conserved small non-coding RNAs that regulate the expression of genes, and they affect important biological and physiological states. Circulating miRNAs in blood are useful markers of metabolism and economic traits. Expression levels of circulating miRNAs have been estimated using quantitative real-time PCR (qPCR). Proper normalization is critical for accurate miRNA expression analysis. However, there is no study which systematically presented endogenous reference genes for evaluating circulating miRNA expression in pigs. In this study, ten porcine miRNAs (let-7a, miR-16, miR-17, miR-23a, miR-26a, miR-93, miR-103, miR-107, miR-127 and miR-191), based on the literature, were chosen as candidate reference miRNAs in serum. We evaluated the expression stability value of these miRNAs in Berkshire, Duroc, Landrace and Yorkshire pigs using geNorm and NormFinder. We determined the optimal combination of reference miRNAs for qPCR experiments: miR-127 and miR-17 in Berkshire pigs; miR-127 and miR-93 in Duroc and Landrace pigs; miR-127 and miR-16 in Yorkshire pigs. miR-127 was the best reference gene in pigs, regardless of the breed. Our study is crucial for the discovery of novel biomarkers in pigs. The reference miRNAs presented in this study could be used as appropriate reference genes for the measurement of circulating miRNA levels in studies of physiological blood metabolites.