GET THE APP

Selection and Optimization of Circulating Exosomal Reference microRNAs for Quantitative Real-time PCR Studies of Porcine Body Fluids
..

Journal of Blood & Lymph

ISSN: 2165-7831

Open Access

Research Article - (2022) Volume 12, Issue 3

Selection and Optimization of Circulating Exosomal Reference microRNAs for Quantitative Real-time PCR Studies of Porcine Body Fluids

Qiang Pu1, Lin Wang2, Guojun Kang3, Changbao Liu4, Changfeng Yang4, Jia Luo4* and Yongfu Huang4*
*Correspondence: Dr. Jia Luo, Department of Animal Science and Technology, Southwest University, Chongqing, China, Tel: +8615283524028, Email: Yongfu Huang, Department of Animal Science and Technology, Southwest University, Chongqing, China, Tel: +8615283524028, Email:
1Chongqing Liu Jiu Animal Husbandry Technology Co., LTD, Chongqing, China
2Chongqing Academy of Animal Sciences, Chongqing, China
3Agriculture and Rural Affairs Committee, Xiushan, Chongqing, China
4Department of Animal Science and Technology, Southwest University, Chongqing, China

Received: 02-Jun-2022, Manuscript No. JBL-22-65694; Editor assigned: 06-Jun-2022, Pre QC No. JBL-22-65694; Reviewed: 20-Jun-2022, QC No. JBL-22-65694; Revised: 27-Jun-2022, Manuscript No. JBL-22-65694; Published: 04-Jul-2022 , DOI: 10.37421/ 2165-7831.2022.12.287
Citation: Pu, Qiang , Lin Wang, Guojun Kang and Changbao Lie, et al. “Selection and Optimization of Circulating Exosomal Reference microRNAs for Quantitative Real-time PCR Studies of Porcine Body Fluids.” J Blood Lymph 12 (2022): 287.
Copyright: © 2022 Pu Q, et al. This is an open-access article distributed under the terms of the creative commons attribution license which permits unrestricted use, distribution and reproduction in any medium, provided the original author and source are credited.

Abstract

Circulating exosomal miRNAs released into all body fluids have incredible functionality and stability. Their expression is associated with multiple pathological conditions, can be used as informative biomarkers when assessing and monitoring the body’s physiopathological status. However, there is no consensus on reference miRNAs for circulating exosomal reference and abundance normalization. The present study aimed to quantify 16 potential reference miRNAs in ten porcine body fluids using qRT-PCR. Further, their stability was quantified by combining multiple gold-standard statistical tools, including BestKeeper, GeNorm, and NormFinder. The identified miRNAs were comprehensively ranked. The top-ranked miRNA was recommended as the optimal reference miRNAs for data normalization. To identify more stable genes, the body fluids were assigned into three groups based on the collection point, they are in vivo (bile, bladder fluid, and gastric juice), in vitro (colostrum, ordinary milk, semen, and urine) and in the blood (UVBP, UABP and PBS). The most stable optimal circulating exosomal reference miRNAs in the body fluids were let-7b-5p (miR-93) in bile, miR-92a in bladder fluid, miR-93 in gastric juice, let-7b-5p in colostrum, miR-92a in ordinary milk and urine, miR-25 in semen, let-7b-5p in UVBP, miR-25 in UABP and U6 in PBS. Overall, miR-93, miR-451 (miR-92a), and miR-25 are the bona fide reference miRNA for qRT-PCR data normalization of body fluids in vivo, in vitro, and blood, respectively. Across all body fluids, miR-451 was the most stable when determining the miRNA abundance in the circulating exosomes

Keywords

Circulating miRNA • Exosomes • Reference genes • Porcine body fluids

Introduction

AmicroRNAs (miRNAs) found in exosomes in various body fluids have been considered promising diagnostic and prognostic biomarkers for diseases. This is due to their relatively high stability and role as post-transcriptional regulators of gene expression [1-4]. Exosomes are membrane-bound vesicles acting as biological cargo carriers of various molecules, including miRNAs, released from cells into the extracellular environment [5]. They are present in different body fluids with incredible stability [6]. The circulating exosomal miRNAs in various body fluids, including the amniotic fluid, breast milk, and blood have been linked to multiple pathological conditions [7,8]. Moreover, fluid-type-specific miRNAs have functional roles on the surrounding tissues [6]. Besides, there is a higher correlation between the clinical-pathological variables with miRNA level in exosomes than with the vesicle-free miRNAs in biofluids. Also, exosomal miRNAs are important biomarkers easily detected non-invasive from all body fluids [9]. However, their detection requires a reliable and quantitative assessment of the miRNA expression.

MiRNAs are short RNA molecules abundant and stable in biofluids, which regulate gene expression in eukaryotic cells. There has been heightened interest in the utility of circulating exosomal miRNAs as minimally invasive diagnostic and monitoring biomarkers for a wide range of human diseases.

However, the quantification of circulating exosomal miRNAs is faced with specific challenges, including the relatively low RNA content in biofluids and lack well-established reference miRNAs [10]. Therefore, there is a need to develop protocols for detecting and quantifying miRNAs in biofluids. Reverse Transcription-quantitative Polymerase Chain Reaction (RT-qPCR) routinely used in detecting the levels of circulating miRNA remains the gold standard for specific detection of selected sets of miRNAs [11]. It easily and with high accuracy quantifies the miRNA content in biological fluids [12]. Unfortunately, the RT-qPCR normalization procedures are sometimes disputable. Generally, the housekeeping molecules with stable expressions, including SNORD44 (RNU44), SNORD48 (RNU48), and nuclear RNA RNU6-1 (U6), are extensively used as reference elements during miRNA quantification in cell and tissue samples. However, there is considerable variation in their expressions in body fluids, such as plasma and serum; thus, it may not be suitable for normalizing the circulating miRNAs [13-16]. Therefore, the selection and optimization of circulating exosomal reference miRNAs in body fluids using RT-qPCR is extremely urgent.

Pigs( Susscrofa ) are emerging as an attractive biomedical model for studying human diseases due to their similar anatomy, genetics, and physiological characteristics, unlike other commonly used mammals and mice [17,18]. Therefore, the present study aimed to determine the suitability of RT-qPCR in quantifying 16 potential reference miRNAs in ten pig body fluids. In addition, we aimed to analyze the differences in miRNAs content in vivo obtained in body fluids (bile, bladder fluid, and gastric juice) following a surgical procedure, in vitro obtained nonsurgically (colostrum, ordinary milk, semen, and urine) and in the blood (Umbilical Venous Blood Plasma (UVBP), Umbilical Arterial Blood Plasma (UABP) and Peripheral Blood Serum (PBS)). The findings in this study will be used in recommending the top-ranked miRNAs as reference miRNAs for data normalization in similar experimental systems.

Materials and Methods

All animal procedures were conducted in compliance with the institutional guidelines for the care and use of laboratory animals. The procedures were reviewed and approved by the Animal Care and Ethics Committee of Southwest University, Chong Qing, China (Permit number 20190723).

Sample collection

Animal experiments were organized by Southwest University and conducted at Chongqing Liu Jiu Animal Husbandry Technology Co., LTD. Ten fluid samples were collected from ten Yorkshire pigs, of which semen were come from breeding boar, colostrum and ordinary milk come from multiparous sows and the rest come from female commercial pigs. They were assigned into three groups based on the location of collection: body fluids in vivo (bile, bladder fluid, and gastric juice), body fluids in vitro (colostrum, ordinary milk, semen, and urine) and the blood (UVBP, UABP and PBS) (n=6). For UVBP and UABP collection, a 15% EDTA anticoagulant solution was prepared and 0.1 ml was added in a 5 ml syringe before their collection at birth. Peripheral blood serum was collected using an EDTA- containing anticoagulant tube. Colostrum and ordinary milk were collected from multiparous sows on the first and seventh days after birth, respectively. Semen was artificially collected from adult boars in the boars' station, while urine entailed the morning urine of adult commercial pigs. Body fluids in vivo (bile, bladder fluid, and gastric juice) were collected during adult pigs' dissection, immediately frozen in liquid nitrogen then transferred to a freezer at-80°C for long-term preservation.

Exosome preparation, candidate Reference miRNAs Se- lectionand and RT-qPCR

The exosomes in the body fluids were isolated as previously described by Luo, et al. [19]. Multiple studies have been conducted to validate the stable expression of reference miRNA genes across different porcine tissues and organs [20,21]. Since circulating exosomal miRNAs in body fluids are mainly related to surrounding tissues and organs, the relatively stable miRNAs in tissues screened previously were used included 16 miRNAs: let-7b-5p, let-7a, miR-103, miR-17-3p, miR-106a, miR-191, miR- 9, miR-25, miR-26a, miR-17-5p, miR-335, miR-148a, miR-451, miR-92a, miR-16 and U6 [22-24].

According to the manufacturer's instructions, the total exosomal RNA was extracted using the TRIzol LS reagent (Invitrogen). The purity of the extracted RNA was determined using a NanoDrop ND2000 spectrophotometer (hermo Scientific, Wilmington, DE, USA). Only RNA samples with a 260/280 ratio between 1.8 and 2.0 purity were selected for subsequent analyses. The obtained cDNA was diluted tenfold with ultrapure water to generate cDNA templates for qRT-PCR. According to the manufacturer's instructions, the reverse transcription of miRNAs was performed using a commercial kit (TaKaRa, China). Serially diluted cDNA templates were analyzed to generate standard curves for each reference gene. All analyses across all the samples were run in triplicate. Quantitative Real-time-PCR (qRT-qPCR) was performed using an SYBR Premix Ex Taq kit (TaKaRa, China) on a Bio-rad IQTM5 system (Bio-Rad, Hercules, CA, USA). The amplification conditions were as follows: initial denaturation at 95°C for 30 s, followed by 40 cycles of denaturation at 95°C for 30s, annealing at 60°C for 40s, and extension at 72°C for 30s. The specificity of each primer pair was verified by melting curves analysis using a corresponding cDNA sample.

Determination of reference gene expression stability

The initial qRT-PCR data were obtained from the StepOnePlus™ software (Applied BioSystems, Foster City, CA, USA) and put into a Microsoft Excel 2010 datasheet. The mean quantification cycle (Cq) value was converted into the relative expression level using the 2−ΔΔCt method during miRNA expression stability analysis. Next, GeNorm, NormFinder and BestKeeper were used to evaluate the gene expression stability [25-27]. The Cq values averages run in triplicate were used to compare the stability of each candidate house-keeping miRNA gene.

Statistical analysis

The statistical analysis of the expression levels of the reference miRNA genes was performed using Duncan multiple comparisons using the SPSS11.0 software. P<0.05 was regarded as a level of statistical significance. The stability of miRNAs in the body fluids was quantified by combining multiple gold-standard statistical tools, including GeNorm, NormFinder and BestKeeper [25-27].

Results

Expression levels of the candidate reference miRNA

The Cycle threshold (Ct) value for each reference miRNA is presented in Supplementary Table 1. There were significant differences across the 16 reference miRNAs in the ten body fluids. The mean Ct value of the candidate genes ranged from 22 to 40 cycles, with miR-148a having the highest transcript levels (22.73 ± 0.02) in ordinary milk and miR-16 with the lowest (39.49 ± 0.14) in urine. Based on the geNorm calculations of the preliminary Ct values and stability sequencing (Supplementary Table 2), ten miRNAs were selected for subsequent stability analysis including, miR- 92a, miR-93, let-7a, miR-451, miR-106a, U6, miR-17-5P, miR-26a, miR-25, let-7b-5p.

GeNorm analysis of candidate reference miRNAs in body fluids

The geNorm analysis revealed that all candidate genes had an M value below 1.5, which is the maximum value for acceptable gene expression stability [25]. The gene with the lowest M value has the most stable expression. However, the optimal candidate miRNA stability expressed in each body fluid differed (Supplementary Table 2). The most stable miRNA were miR-25 and let-7b-5p in URBP, let-7b-5p and miR-93 in UVBP, miR- 26a and miR-335 in PBS, miR-17-5P and miR-148a in ordinary milk, miR- 451 and miR-16 in colostrum, miR-103 and miR-451 in semen, miR-191 and miR-92a in urine, let-7b-5p and miR-93 in bile, miR-191 and miR-16 in bladder fluid, and miR-17-5p and miR-93 in gastric juice. Further analyses of the M value ranking in all body fluids in vivo, in vitro, and the blood revealed the first five ranking of stability in all body fluids as miR-451>miR- 92a (let-7b-5p)>miR-106a (U6). The average stability rank in the blood, in vitro and in vivo was let-7b-5p>miR-25>U6, miR-451>miR-92a>miR-106a, and U6>miR-92a>let-7a, respectively.

Analysis of candidate reference miRNAs in blood

To explore the stability of miRNA expression, the stability of miRNA in URBP, UVBP, and PBS were analyzed (Table 1). A comprehensive ranking of the miRNA expression stability revealed that let-7b-5p expressed in URBP was the most stable, ranking 1, 2, and 1 with GeNorm, NormFinder, and BestKeeper, respectively. Moreover, miR-25 expression was the most stable in UVBP and PBS. The U6 expression was also stable in PBS (Table 1).

Table 1: Ranking and best combination of candidate reference genes in blood.

URBP Genormer (M value) Normfinder stability value Bestkeeper CV (%CP) UVBP Genormer (M value) Normfinder stability value Bestkeeper CV (%CP) PBS Genormer (M value) Normfinder stability value Bestkeeper CV (%CP)
miR-92a 0.12 (10) 0.11 (10) 0.27 (8) miR-17-5P 0.10 (10) 0.09 (10) 0.03 (2) miR-451 0.11 (10) 0.11 (10) 0.33 (10)
miR-93 0.11 (9) 0.10 (9) 0.32 (9) miR-92a 0.09 (9) 0.08 (9) 0.31 (8) miR-106a 0.09 (9) 0.10 (9) 0.29 (9)
let-7a 0.10 (8) 0.09 (8) 0.34 (10) miR-26a 0.08 (8) 0.08 (8) 0.31 (7) miR-93 0.08 (8) 0.08 (8) 0.24 (8)
miR-451 0.08 (7) 0.07 (7) 0.23 (7) miR-106a 0.07 (7) 0.04 (5) 0.19 (5) miR-92a 0.06 (7) 0.04 (6) 0.12 (6)
miR-106a 0.06 (6) 0.05 (6) 0.12 (5) miR-451 0.07 (6) 0.04 (6) 0.31 (9) let-7b-5p 0.06 (6) 0.05 (7) 0.13 (7)
U6 0.06 (5) 0.04 (5) 0.11 (4) miR-93 0.06 (5) 0.02 (1)* 0.20 (6) miR-26a 0.04 (5) 0.01 (3) 0.03 (2)
miR-17-5P 0.04 (4) 0.01 (1)* 0.19 (6) let-7b-5p 0.06 (4) 0.02 (2) 0.16 (4) let-7a 0.03 (4) 0.01 (2) 0.10 (3)
miR-26a 0.03 (3) 0.04 (4) 0.10 (3) U6 0.04 (3) 0.05 (7) 0.15 (3) U6* 0.03 (3) 0.01 (1)* 0.02 (1)*
miR-25 0.02 (1)* 0.02 (3) 0.09 (2) miR-25* 0.03 (1)* 0.03 (3) 0.01 (1)* miR-25* 0.01 (1)* 0.04 (5) 0.11 (4)
let-7b-5p* 0.02 (1)* 0.01 (2) 0.08 (1)* let-7a 0.03 (1)* 0.03 (4) 0.38 (10) miR-17-5P 0.01 (1)* 0.04 (4) 0.12 (5)

Analysis of the candidate reference miRNAs in body flu- ids in vitro

Except for blood, commonly used to detect disease, body fluids in vitro are more readily available than in vivo. Here, the stability of miRNA expression in body fluids in vitro (ordinary milk, colostrum, semen, and urine) was analyzed. The findings are presented in Table 2. A comprehensive ranking presented miR-106a, let-7b-5p, miR-25, and miR-92a as the most stably expressed in ordinary milk (ranks 1, 1 and 3), colostrum (ranks 1, 1 and 1), semen (ranks 5, 1 and 3), and urine (ranks 1, 2 and 1) using GeNorm, NormFinder, and BestKeeper, respectively (Table 2).

Table 2: Ranking and best combination of candidate reference genes of body fluids in vitro.

Ordinary milk Genormer (M value) Normfinder stability value Bestkeeper CV (% CP) Colostrum Genormer (M value) Normfinder stability value Bestkeeper CV (%CP) Semen Genormer (M value) Normfinder stability value Bestkeeper CV (%CP) Urine Genormer (M value) Normfinder stability value Bestkeeper CV (%CP)
let-7a 0.06 (10) 0.07 (10) 0.36 (9) let-7a 0.09 (10) 0.08 (9) 0.37 (10) let-7b-5p 0.12 (10) 0.12 (10) 0.27 (8) U6 0.11 (10) 0.15 (10) 0.38 (9)
miR-25 0.05 (9) 0.06 (9) 0.25 (8) miR-25 0.08 (9) 0.08 (10) 0.22 (8) let-7a 0.10 (9) 0.08 (9) 0.25 (7) miR-25 0.08 (9) 0.10 (9) 0.29 (7)
U6 0.05 (8) 0.05 (8) 0.2 (7) miR-93 0.07 (8) 0.08 (8) 0.22 (9) miR-93 0.09 (8) 0.08 (8) 0.24 (5) let-7a 0.07 (8) 0.08 (8) 0.30 (8)
miR-26a 0.04 (7) 0.03 (6) 0.07 (1)* miR-26a 0.06 (7) 0.06 (7) 0.13 (7) miR-92a 0.08 (7) 0.04 (4) 0.10 (2) miR-106a 0.06 (7) 0.05 (7) 0.10 (4)
miR-93 0.04 (6) 0.03 (7) 0.12 (6) miR-17-5P 0.05 (6) 0.03 (5) 0.08 (4) miR-451 0.08 (6) 0.03 (3) 0.08 (1)* miR-26a 0.05 (6) 0.03 (6) 0.05 (2)
let-7b-5p 0.03 (5) 0.02 (5) 0.10 (4) miR-106a 0.04 (5) 0.01 (2) 0.07 (3) miR-25* 0.07 (5) 0.01 (1)* 0.11 (3) let-7b-5p 0.04 (5) 0.02 (5) 0.11 (6)
miR-17-5P 0.02 (4) 0.01 (3) 0.37 (10) miR-92a 0.04 (4) 0.05 (6) 0.11 (6) miR-26a 0.06 (4) 0.01 (1)* 0.14 (4) miR-93 0.04 (4) 0.02 (4) 0.11 (5)
miR-451 0.02 (3) 0.02 (4) 0.11 (5) U6 0.03 (3) 0.03 (4) 0.09 (5) miR-17-5P 0.03 (3) 0.08 (7) 0.32 (10) miR-451 0.03 (3) 0.01 (1)* 0.09 (3)
miR-106a* 0.01 (1)* 0.01 (1)* 0.08 (3) miR-451 0.02 (2) 0.02 (3) 0.02 (2) miR-106a 0.02 (1)* 0.05 (5) 0.24 (6) miR-92a* 0.02 (1)* 0.01 (2) 0.01 (1)*
miR-92a 0.01 (1)* 0.01 (2) 0.07 (2) let-7b-5p* 0.02 (1)* 0.01 (1)* 0.02 (1)* U6 0.02 (1)* 0.06 (6) 0.32 (9) miR-17-5P 0.02 (1)* 0.01 (3) 0.38 (10)

Analysis of the candidate reference miRNAs in body flu- ids in vivo

The stability of miRNA expressed in body fluids in vitro, including in bile, bladder fluid, and gastric juice, is presented in Table 3. A comprehensive ranking of the miRNA expression stability revealed that miR-93 (let-7b-5p), miR-92a, and miR-93 were most stably expressed in bile (ranks 1, 2, and 1) , bladder fluid (ranks 5, 2, and 2), and gastric juice (ranks 1, 2 and 2) based on the GeNorm, NormFinder and BestKeeper analyses, respectively (Table 3).

Table 3: Ranking and best combination of candidate reference genes of body fluids in vivo.

Bile Genormer (M value) Normfinder stability value Bestkeeper CV (%CP) Bladder fluid Genormer (M value) Normfinder stability value Bestkeeper CV (%CP) Gastric juice Genormer (M value) Normfinder stability value Bestkeeper CV (%CP)
miR-451 0.12 (10) 0.11 (10) 0.40 (10) miR-93 0.11 (10) 0.11 (10) 0.28 (10) miR-26a 0.07 (10) 0.09 (10) 0.27 (10)
miR-25 0.11 (9) 0.10 (9) 0.28 (9) miR-17-5P 0.10 (9) 0.10 (9) 0.13 (3) let-7b-5p 0.05 (9) 0.08 (9) 0.20 (9)
miR-17-5P 0.09 (8) 0.07 (8) 0.10 (3) let-7b-5p 0.08 (8) 0.08 (7) 0.17 (6) let-7a 0.04 (8) 0.03 (7) 0.10 (6)
miR-92a 0.09 (7) 0.03 (3) 0.15 (6) miR-25 0.08 (7) 0.08 (8) 0.25 (9) miR-25 0.04 (7) 0.03 (8) 0.13 (8)
miR-26a 0.08 (6) 0.06 (4) 0.15 (7) U6 0.07 (6) 0.03 (3) 0.04 (1)* U6 0.03 (6) 0.02 (6) 0.07 (4)
miR-106a 0.07 (5) 0.07 (7) 0.15 (5) miR-92a* 0.06 (5) 0.02 (2) 0.06 (2) miR-92a 0.03 (5) 0.02 (5) 0.08 (5)
U6 0.05 (4) 0.06 (6) 0.15 (4) let-7a 0.04 (4) 0.01 (1)* 0.19 (8) miR-451 0.02 (4) 0.01 (1)* 0.02 (1)*
let-7a 0.04 (3) 0.06 (5) 0.22 (8) miR-451 0.01 (3) 0.03 (4) 0.15 (4) miR-106a 0.02 (3) 0.01 (3) 0.04 (3)
let-7b-5p* 0.04 (1)* 0.03 (2) 0.04 (1)* miR-106a 0.01 (1)* 0.04 (5) 0.17 (5) miR-93* 0.02 (1)* 0.01 (2) 0.04 (2)
miR-93* 0.04 (1)* 0.02 (1)* 0.05 (2) miR-26a 0.01 (1)* 0.04 (6) 0.18 (7) miR-17-5P 0.02 (1)* 0.01 (4) 0.12 (7)

Comprehensive ranking of the reference miRNAs

Given that each body fluid has its distinct characteristics, it is necessary to identify the optimal reference miRNAs in each fluid for a better practical application. The comprehensive comparative analysis revealed that miR- 25, let-7b-5p, miR-25, and U6 were the optimal reference miRNA in blood, URBP, UVBP, and PBS, respectively (Table 4). The optimal reference miRNA in body fluids obtained in vitro was miR-451 and miR-92a, with miR- 92a, let-7b-5p, miR-25, and miR-92a as the optimal reference miRNA in ordinary milk, colostrum, semen, and urine, respectively. In addition, the optimal reference miRNA in body fluids in vivo was miR-93, with let-7b-5p (miR-93), miR-92a, and miR-93 as the optimal reference miRNA in bile, bladder fluid, and gastric juice, respectively. Overall, miR-25, miR-92, and miR-93 were the most common in blood, in vivo, and in vitro, respectively (Table 4).

Table 4: The Optimal reference gene in each body fluids calculated by geNorm, NormFinder and Bestkeeper algorithms.

Type of fluid URBP UVBP PBS  Blood  Ordinary milk Colostrum Semen Urine  Body fluids in vitro Bile Bladder fluid Gastric juice Body fluids in vivo
Optimal reference gene let-7b-5p miR-25 U6 miR-25 miR-92a let-7b-5p miR-25 miR-92a miR-451 let-7b-5p miR-92a miR-93 miR-93
                  miR-92a miR-93      

Table S1: Optimal reference gene and expression stability values (M value) of the candidate reference miRNAs calculated by geNorm (ranking in parentheses).

Sample source miR-103 miR-106a miR-148a miR-17-3P miR-191 miR-451 miR-92a miR-16 miR-25 miR-26a miR-335 let-7b-5p miR-93 U6 miR-17-5P let-7a
Mean ± SD
URBP 36.35 ± 0.13 34.59 ± 0.06 34.27 ± 0.23 37.51 ± 0.07 34.83 ± 0.02 31.18 ± 0.10 35.70 ± 0.13 32.81 ± 0.04 37.65 ± 0.05 37.16 ± 0.06 38.25 ± 0.13 38.37 ± 0.04 35.41 ± 0.15 29.19 ± 0.05 33.28 ± 0.10 33.39 ± 0.17
UVBP 30.28 ± 0.14 28.52 ± 0.08 28.99 ± 0.07 35.21 ± 0.09 32.19 ± 0.03 25.56 ± 0.11 30.53 ± 0.13 33.17 ± 0.11 31.03 ± 0.01 33.73 ± 0.16 37.59 ± 0.15 35.33 ± 0.08 30.14 ± 0.08 29.64 ± 0.07 33.41 ± 0.02 33.54 ± 0.17
PBS 37.45 ± 0.12 33.81 ± 0.15 35.43 ± 0.19 34.12 ± 0.09 34.36 ± 0.18 28.85 ± 0.13 36.85 ± 0.07 36.20 ± 0.06 34.32 ± 0.06 37.05 ± 0.02 36.05 ± 0.03 38.07 ± 0.07 35.51 ± 0.12 29.37 ± 0.01 33.64 ± 0.06 34.21 ± 0.05
Ordinary milk 25.52 ± 0.03 24.62 ± 0.03 22.73 ± 0.02 30.46 ± 0.15 23.27 ± 0.05 29.18 ± 0.05 24.56 ± 0.03 27.94 ± 0.07 24.84 ± 0.09 26.13 ± 0.02 36.11 ± 0.08 24.41 ± 0.03 25.10 ± 0.04 22.97 ± 0.07 24.38 ± 0.12 23.24 ± 0.11
Colostrum 33.07 ± 0.13 31.69 ± 0.04 31.79 ± 0.12 39.55 ± 0.16 31.69 ± 0.13 31.14 ± 0.01 32.89 ± 0.06 36.31 ± 0.01 33.29 ± 0.11 33.98 ± 0.07 37.58 ± 0.14 32.99 ± 0.01 33.23 ± 0.11 28.39 ± 0.04 30.10 ± 0.04 27.26 ± 0.14
Semen 35.10 ± 0.02 37.27 ± 0.13 33.03 ± 0.04 39.57 ± 0.09 34.53 ± 0.14 31.76 ± 0.04 34.15 ± 0.05 38.49 ± 0.14 34.42 ± 0.06 35.28 ± 0.07 36.11 ± 0.09 32.47 ± 0.13 34.91 ± 0.12 32.43 ± 0.14 32.43 ± 0.14 29.50 ± 0.11
Urine 37.45 ± 0.15 36.20 ± 0.05 36.36 ± 0.13 33.47 ± 0.18 37.31 ± 0.01 33.2 ± 0.05 38.54 ± 0.01 39.49 ± 0.14 38.21 ± 0.15 36.19 ± 0.03 36.92 ± 0.12 37.77 ± 0.06 38.48 ± 0.06 32.47 ± 0.18 32.47 ± 0.18 33.41 ± 0.15
Bile 33.38 ± 0.05 34.14 ± 0.07 36.19 ± 0.18 38.43 ± 0.13 33.65 ± 0.13 33.49 ± 0.17 36.09 ± 0.08 38.38 ± 0.10 36.27 ± 0.15 35.99 ± 0.08 38.80 ± 0.10 34.16 ± 0.02 34.35 ± 0.03 32.01 ± 0.07 34.05 ± 0.05 31.25 ± 0.09
Bladder fluid 36.03 ± 0.04 33.58 ± 0.08 36.26 ± 0.09 36.12 ± 0.10 33.46 ± 0.02 32.44 ± 0.07 35.80 ± 0.03 36.47 ± 0.04 36.89 ± 0.12 35.05 ± 0.09 35.43 ± 0.09 36.21 ± 0.09 35.27 ± 0.15 29.14 ± 0.02 30.66 ± 0.06 32.12 ± 0.09
Gastric juice 32.79 ± 0.14 31.15 ± 0.02 32.62 ± 0.05 38.31 ± 0.10 33.54 ± 0.16 24.13 ± 0.01 32.26 ± 0.04 32.93 ± 0.08 30.22 ± 0.06 34.18 ± 0.13 38.75 ± 0.15 30.10 ± 0.09 32.39 ± 0.02 28.83 ± 0.03 30.98 ± 0.06 29.08 ± 0.04

Table S2. Optimal reference gene and expression stability values (M value) of the candidate reference miRNAs calculated by geNorm (ranking in parentheses).

  URBP UVBP PBS Ordinary milk Colostrum Semen Urine Bile Bladder fluid Gastric juice The average rank in all Average rank in blood Average rank of body fluids in vitro Average rank of body fluids in vivo
Optimal reference gene miR-25 let-7b-5p miR-26a miR-148a miR-451 miR-103 miR-191 let-7b-5p miR-191 miR-17-5P  
  let-7b-5p miR-93 miR-335 miR-17-5P miR-16 miR-451 miR-92a miR-93 miR-16 miR-93
miR-148a 0.13 (15) 0.08 (9) 0.13 (15) 0.01 (1)* 0.06 (8) 0.06 (5) 0.08 (9) 0.13 (15) 0.09 (13) 0.04 (8) 9.8 13 5.75 12
miR-92a 0.11 (14) 0.09 (11) 0.06 (7) 0.02 (3) 0.04 (5) 0.05 (4) 0.01 (1) 0.08 (6) 0.03 (3) 0.03 (4) 5.8 10.67 3.25 4.33
miR-93 0.10 (13) 0.02 (1)* 0.10 (12) 0.03 (7) 0.06 (7) 0.07 (7) 0.04 (4) 0.04 (1)* 0.11 (15) 0.02 (1)* 6.8 8.67 6.25 5.67
miR-103 0.10 (12) 1.12 (15) 0.08 (10) 0.03 (5) 0.07 (9) 0.02 (1)* 0.14 (15) 0.06 (4) 0.03 (4) 0.08 (12) 8.7 12.33 7.5 6.67
miR-335 0.09 (11) 0.10 (12) 0.01 (1)* 0.07 (14) 0.12 (15) 0.09 (10) 0.07 (8) 0.07 (5) 0.08 (9) 0.10 (14) 9.9 8 11.75 9.33
let-7a 0.08 (10) 0.07 (6) 0.03 (3) 0.06 (13) 0.09 (11) 0.11 (14) 0.09 (10) 0.04 (2) 0.04 (5) 0.04 (7) 8.5 6.33 13 4.67
miR-451 0.07 (9) 0.03 (2) 0.09 (11) 0.02 (4) 0.01 (1)* 0.02 (1)* 0.03 (3) 0.12 (14) 0.06 (6) 0.02 (3) 5.40* 7.33 2.25* 7.67
miR-17-3P 0.06 (8) 0.05 (3) 0.07 (9) 0.08 (15) 0.11 (14) 0.06 (6) 0.12 (13) 0.11 (12) 0.08 (10) 0.07 (11) 10.1 6.67 12 11
miR-106a 0.06 (7) 0.06 (4) 0.11 (13) 0.01 (2) 0.03 (4) 0.08 (8) 0.05 (7) 0.09 (7) 0.07 (7) 0.02 (2) 6.1 8 5.25 5.33
U6 0.05 (6) 0.08 (7) 0.02 (2) 0.05 (10) 0.02 (3) 0.09 (9) 0.13 (14) 0.05 (3) 0.02 (2) 0.03 (5) 6.1 5 9 3.33*
miR-16 0.05 (5) 0.09 (10) 0.06 (8) 0.06 (12) 0.01 (1)* 0.10 (11) 0.11 (12) 0.10 (9) 0.02 (1)* 0.05 (9) 7.8 7.67 9 6.33
miR-17-5P 0.04 (4) 0.10 (13) 0.04 (5) 0.01 (1)* 0.04 (6) 0.11 (13) 0.01 (2) 0.10 (10) 0.10 (14) 0.02 (1)* 6.9 7.33 5.5 8.33
miR-26a 0.04 (3) 0.11 (14) 0.01 (1)* 0.04 (8) 0.08 (10) 0.04 (3) 0.04 (6) 0.09 (8) 0.07 (8) 0.09 (13) 7.4 6 6.75 9.67
miR-191 0.03 (2) 0.08 (8) 0.12 (14) 0.04 (9) 0.10 (13) 0.10 (12) 0.01 (1)* 0.10 (11) 0.02 (1)* 0.11 (15) 8.6 8 8.75 9
miR-25 0.02 (1)* 0.07 (5) 0.04 (4) 0.05 (11) 0.09 (12) 0.04 (2) 0.10 (11) 0.11 (13) 0.09 (12) 0.04 (6) 7.7 3.33 9 10.33
let-7b-5p 0.02 (1)* 0.02 (1)* 0.05 (6) 0.03 (6) 0.01 (2) 0.12 (15) 0.04 (5) 0.04 (1)* 0.08 (11) 0.06 (10) 5.8 2.67* 7 7.33

Discussion

Exosomes are nanovesicles released from almost all cell types under physiological and pathological conditions. Due to their unusually high stability in biofluids, circulating exosomal miRNAs have emerged as promising non- invasive biomarkers for human disease [28]. MicroRNA expression profiling in body fluids, specifically the plasma and the serum, is key to identifying miRNA biomarkers informing early disease diagnosis and predicting response to therapy [29]. However, ribonucleases in body fluids degrade RNA, including the extracellular miRNA; thus, obtaining miRNAs expression profiles in biofluids is technically challenging [30]. Several reference genes for miRNA expression profiles have been identified in tissues [31] and blood in the umbilical cord [32] using different analysis strategies [33]. Besides, the reference genes in exosomes have been screened [34]. However, the miRNA stability in the specific body fluids exosomes remains unknown. In addition, there are no reference genes identified in the circulating exosomal miRNA in body fluids. Given the protective effect of exosomes on RNA, the circulating exosomal RNA can be used as biomarkers for biological humoral diseases [30,35,36]. This has increased the interest in extracellular RNA, specifically in circulating exosomal miRNA. As a result, there is a growing interest in identifying these circulating exosomal miRNAs as non-invasive biomarkers for the early detection of diseases or selection of treatment options.

Circulating exosome miRNA studies using qRT-PCR depends on an endogenous reference miRNA to normalize the abundance of the miRNAs [37]. In the present study, we analyzed the stability of miRNA in ten common body fluids, including those obtained in vivo, in vitro, and the blood. To identify the most suitable reference exosomal miRNA for stable expression in the different body fluids, a sequencing analysis of miRNA present in all body fluids was conducted to screen out the optimal circulating exosomal miRNAs. We performed a preliminary analysis on all the 16 candidate miRNAs using GeNorm. We selected the top ten stably expressed in the body fluids for subsequent stability analysis. The most suitable candidate reference miRNAs in the different body fluids were revealed by geNorm analysis. There were huge differences in miRNAs stability between the different body fluids. For example, the most stable miRNA in all body fluids was miR-451. However, let-7b-5p was the most stable in the blood, miR- 451 in body fluids obtained in vitro and U6 in body fluids obtained in vivo.

This is consistent with previous studies, which revealed that miR-451had strong stability; thus, it could be detected at all detection points [38]. In contrast, U6 has been recommended as the reference gene for the relative quantification of miRNA expression levels [39].

Based on the GeNorm analysis, the most stable miRNA in blood, body fluids obtained in vitro and in vivo were let-7b-5p>miR-25>U6, miR- 451>miR-92a>miR-106a, and U6 >miR-92a>let-7a, respectively. Moreover, let-7b-5p was the most stably expressed in URBP and miR-25 in UVBP and PBS in blood. This finding is consistent with previous studies where let-7b- 5p and miR-25 were identified as stably expressed genes in blood [40,41]. Moreover, miR-106a, let-7b-5p, miR-25, miR-92a, MiR-93 (let-7b-5p), miR-92a, miR-93 were the most stably expressed miRNA in ordinary milk, colostrum, semen, urine, bile, bladder fluid, and gastric juice, respectively. This provides a reference for selecting marker genes during specific physiological conditions.

A comprehensive analysis revealed that the optimal reference miR-25, miR-92, and miR-93 were common in body fluids obtained in the blood, in vivo, and in vitro, respectively. This result is consistent with previous findings where miR-92 was identified as an internal reference for stable expression in body fluids and miR-93 as a suitable reference gene for serum miRNA analysis in gastric cancer patients [33,42]. The preliminary findings in this study proved a huge potential in body fluid-based miRNAs used as references in diagnosing diseases and predicting results following a therapeutic intervention.

Conclusion

In the present study, we analyzed the stability of 16 miRNAs obtained from exosomes in ten body fluids and optimized the reference miRNAs. Based on the comprehensive ranking, the optimal reference genes in vivo were let-7b-5p (miR-93) in bile, miR-92a in bladder fluid, and miR-93 in gastric juice, while let-7b-5p in colostrum, miR-92a in ordinary milk and urine, and miR-25 in semen were the optimal reference genes in body fluids in vitro and let-7b-5p in UVBP, miR-25 in UABP and U6 in PBS in the blood. Across all body fluids, miR-451 was the most stable circulating exosomal reference miRNA.

Statements and Declarations

Funding

The project was supported by Chongqing Natural Science Foundation (No. cstc2020jcyj-msxmX0961) and China Agriculture Research System of MOF and MARA.

Competing interests

The authors have no relevant financial or non-financial interests to disclose.

Author contribution

The authors contribute are as follows: the conception and design of the study (QiangPu, Lin Wang), acquisition of data (QiangPu, Lin Wang, Guojun Kang, Changbao Liu), analysis and interpretation of data (QiangPu, Lin Wang, Changfeng Yang), drafting the article (QiangPu, Lin Wang), final approval of the version to be submitted (JiaLuo, Yongfu Huang).

Ethics approval

This study was granted by the Ethics Committee of Southwest University (No.20190723).

References

Google Scholar citation report
Citations: 443

Journal of Blood & Lymph received 443 citations as per Google Scholar report

Journal of Blood & Lymph peer review process verified at publons

Indexed In

 
arrow_upward arrow_upward