Identification of reference miRNAs in plasma useful for the study of oestrogen-responsive miRNAs associated with acquired Protein S deficiency in pregnancy
BMC Research Notes volume 10, Article number: 312 (2017)
Accumulating evidence indicate that circulating microRNAs (miRNAs) are useful independent non-invasive biomarkers, with unique miRNA signatures defined for various pathophysiological conditions. However, there are no established universal housekeeping miRNAs for the normalisation of miRNAs in body fluids. We have previously identified an oestrogen-responsive miRNA, miR-494, in regulating the anticoagulant, Protein S, in HuH-7 liver cells. Moreover, increased thrombotic risk associated with elevated circulating oestrogen levels is frequently observed in pregnant women and oral contraceptive users. In order to identify other oestrogen-responsive miRNAs, including miR-494, that may be indicative of increased thrombotic risk in plasma, we used nanoString analysis to identify robust and stable endogenous reference miRNAs for the study of oestrogen-responsive miRNAs in plasma.
We compared the plasma miRNA expression profile of individuals with: (1) Low circulating oestrogens (healthy men and non-pregnant women not taking oral contraceptives), (2) High circulating synthetic oestrogens, (women taking oral contraceptives) and (3) High circulating natural oestrogens (pregnant females >14 weeks gestation). From the nanoString analyses, 11 candidate reference miRNAs which exhibited high counts and not significantly differentially expressed between groups were selected for validation using realtime quantitative polymerase chain reaction (RT-qPCR) and digital droplet PCR (DDPCR) in pooled plasma samples, and the stability of their expression evaluated using NormFinder and BestKeeper algorithms. Four miRNAs (miR-25-5p, miR-188-5p, miR-222-3p and miR-520f) demonstrated detectable stable expression between groups and were further analysed by RT-qPCR in individual plasma samples, where miR-188-5p and miR-222-3p expression were identified as a stable pair of reference genes. The miRNA reference panel consisting of synthetic spike-ins cel-miR-39 and ath-miR159a, and reference miRNAs, miR-188-5p and miR-222-3p was useful in evaluating fold-change of the pregnancy-associated miRNA, miR-141-3p, between groups.
The miRNA reference panel will be useful for normalising qPCR data comparing miRNA expression between men and women, non-pregnant and pregnant females, and the potential effects of endogenous and synthetic oestrogens on plasma miRNA expression.
MiRNAs are short, non-coding RNA species approximately 22 nucleotides in length that function to downregulate the expression of its target genes by binding to specific sequences in the 3′ untranslated region (UTR) of target mRNAs leading to the destabilisation of the mRNA or translational repression. Since their discovery in 1993 [1, 2], the role of miRNAs in regulating key cellular processes and their involvement in various pathophysiological conditions, particularly in solid tumours are well studied, and have been extensively reviewed [3,4,5]. It has also been demonstrated that miRNA profiles for human tumours are superior to mRNA markers, where poorly differentiated tumours could be successfully classified using miRNA expression profiles, whereas mRNA profiles were highly inaccurate when applied to the same tissue samples . A comparison of miRNA and gene expression in frozen breast cancer tissue samples and their paired formalin-fixed paraffin-embedded tissue samples obtained during primary surgical resection showed that expression of miRNAs remains robust even in samples with degraded total RNA, clearly demonstrating the stability of miRNAs in compromised samples in comparison to mRNA markers .
The value of miRNAs as useful biomarkers was elevated when plasma and serum were found to be rich in miRNAs that appeared to be protected from RNase degradation [8, 9], and remained highly stable at room temperature and under adverse conditions such as multiple freeze–thaw cycles [10, 11]. A large variety of miRNAs have been detected in various fractions of circulating blood  such as erythrocytes , anucleate platelets , apoptotic bodies  and packaged in circulating microparticles shed from platelets and all other cells in the body . As such, we hypothesise that circulating levels of tissue specific miRNAs can provide important information regarding the health status of target tissues, serving as useful non-invasive biomarkers. Due to significant protein and lipid variability between individual plasma and serum samples, a spike-in control is often introduced to correct for variations in RNA extraction efficiency and contamination from inhibitors of PCR (Kroh , Sourvinou ); although spike-ins have their own caveats. Currently there are no universal endogenous miRNA controls for the study of circulating miRNAs; instead, individual studies arbitrarily designating miRNAs for normalisation.
We demonstrated in HuH-7 liver carcinoma cells that oestrogen treatment resulted in the induction of miR-494 expression in association with downregulated levels of the anticoagulant, Protein S, indicating a mechanism for miRNA-mediated regulation of acquired Protein S deficiency, with a potential role for increased miR-494 levels in pregnant women as an indicator for high thrombotic risk . It is likely that increased circulating oestrogen levels during pregnancy alters the expression of a suite of miRNAs in oestrogen-sensitive tissues, with changes in specific miRNAs possibly detectable in circulating plasma. To this end, we set out to identify stable populations of miRNA which can act as endogenous normalisation markers in human plasma for investigating oestrogen-responsive circulating miRNAs.
All participants provided written informed consent for the use of venous blood samples in this study prior to blood collection. This project was approved by the Human Research Ethics Committee at Murdoch University (HREC Ref: 2013/131).
Participants and plasma preparation
Peripheral whole blood was collected from healthy males, non-pregnant females not taking oral contraceptives, non-pregnant females taking oral contraceptives and pregnant females >14 weeks gestation in EDTA Vacutainer® tubes (Becton–Dickinson, NJ, USA). Whole blood was centrifuged at 2000×g for 15 min within 4 h of blood collection to separate plasma from cells and platelets as previously described , and the supernatant plasma was transferred into RNase-free tubes and stored at −80 °C until required for RNA extraction.
Total RNA, including small RNAs was extracted using the miRVANA PARIS RNA extraction kit (Life Technologies, CA, USA), according to the manufacturer’s protocol. For plasma samples, each sample was spiked with 5 μL of 200 nM synthetic cel-miR-39 and ath-miR159a (Sigma-Aldrich, MO, USA) after RNases were inactivated. The extraction protocol was modified with an extended centrifugation time of 60 min at 13,000×g at the organic extraction step, and each RNA sample was eluted in 50 μL of hot elution (95 °C) buffer.
NanoString nCounter® miRNA Assay
miRNA profiling of total RNA extracted from human plasma samples was performed using the nCounter® Human miRNA Panel v2 (NanoString, Seattle, USA). Plasma total RNA samples were concentrated using a 3 kDa size exclusion filter (Sigma-Aldrich, MO, USA) and 3 μL of concentrated total RNA was hybridised to probes, purified and counted on the nCounter® Prep Station and Digital Analyser according to the manufacturer’s instructions. For each sample, 555 fields of view (FOV) were imaged and counted by the nCounter® Digital Analyser.
Availability of data and materials
All supporting data for all samples analysed by NanoString nCounter® miRNA Assay is presented in this article or as Additional file 1.
Reverse transcription and miRNA expression studies by real-time PCR (RT-PCR) and droplet digital PCR (DDPCR)
Total RNA samples were reverse transcribed using the TaqMan® microRNA reverse transcription kit (Life Technologies, CA, USA), and selected mature miRNAs were detected using TaqMan® miRNA assays (Table 1) (Life Technologies, CA, USA) on the CFX96/384 real-time PCR system (Bio-Rad, CA, USA). Absolute quantitation of miRNA concentration in plasma samples was directly determined using the same TaqMan® miRNA assays and droplet digital PCR (DDPCR) (Bio-Rad, CA, USA). The miRNA concentrations (copies/μL plasma) were calculated using the formula: (copy number of miRNA per microlitre of DDPCR reaction mixture) × (the dilution factor in the plasma with RNA extraction, reverse transcription and DDPCR reaction), and finally normalised to the concentration of cel-miR-39 and ath-miR159a.
Statistical analysis and evaluation of candidate reference miRNA stability
The raw counts obtained from the nCounter® Digital Analyser were compared overall among the four cohorts and via pairwise cohort comparisons based on quasi-likelihood regression with a log link and variance proportional to the square of the mean , which upon inspection indicated was an appropriate model. The regressions adjusted simultaneously for the sum of the logs of the internal positive controls and the synthetic spike-in miRNA, ath-miR159a. Candidate housekeeping miRNAs were identified if they had raw counts >40 and were not significantly different between samples analysed (p > 0.05). Analyses were carried out in TIBCO Spotfire S+ ver 8.2 (Boston, MA).
The stability of miRNA expression was evaluated using NormFinder  and Bestkeeper  software with averaged Cq values to select the most stable reference miRNAs. The 2−∆∆Ct method was used to calculate the expression of the target miRNAs relative to suitable reference miRNAs, including the two synthetic spike-in miRNAs, cel-miR-39 and ath-miR159a.
Identification of candidate housekeeping miRNAs by nanoString nCounter® miRNA assay
In order to identify suitable candidate housekeeping miRNAs for studying oestrogen-responsive miRNAs in plasma, the expression profiles of 800 miRNAs in human plasma samples collected from healthy males, non-pregnant females not taking oral contraceptives, non-pregnant females currently taking oral contraceptives and pregnant females (>14 weeks gestation) were obtained using the nanoString nCounter® miRNA assay v2, which comprised the screening phase (Fig. 1). The candidate housekeeping miRNAs were then validated using RT-qPCR and DDPCR for oestrogen-non-responsiveness between the cohorts in total RNA extracted from pooled plasma samples and the miRNAs identified as maintaining robust expression and concordance with nanoString were further analysed in individual plasma samples by RT-qPCR (Fig. 1).
Out of the 800 miRNAs analysed using the nanoString nCounter® miRNA assay v2, 760 miRNAs had average counts less than 46, of which ~600 miRNAs showed average nanoString counts of less than 20. As such, the initial selection criterion for classifying detectable miRNAs at the screening phase was set at 40 counts to include miRNAs that were considered “low positives”. The results showed that 44 out of 800 miRNAs analysed showed average nanoString counts greater than 40 in all 16 plasma samples, of which the counts for 11 miRNAs were not significantly different by pairwise comparisons (p > 0.05) and identified as suitable candidate reference miRNAs for further analysis (Table 1). Of the 11 candidate miRNAs, miR-25-3p was the most abundant miRNA in the plasma samples analysed, followed by miR-520f, miR-149-5p, let-7a-5p, miR-2682-5p, miR-612, miR-222-3p, miR-598, miR-188-5p, miR-489 and lastly miR-1183 exhibiting the lowest average nanoString counts (Table 1).
Validation of candidate housekeeping miRNAs by qPCR and DDPCR in pooled and individual plasma samples
To validate the stability of the selected candidate housekeeping miRNAs, equal volumes of plasma samples collected were pooled within each cohort and total RNA was extracted from the pooled plasma, reverse transcribed and analysed by qPCR and DDPCR, the whole process performed twice as technical replicates. There was poor concordance of miRNA expression between nanoString and qPCR results, miRNAs which showed high counts in nanoString, miR-149-5p, miR-489, miR-612, miR-1183 and miR-2682 were not amplified under 40 cycles using qPCR TaqMan® microRNA assays. This observation was matched with low miRNA concentrations using DDPCR methods (Fig. 2).
Based on Cq values, miR-222-3p exhibited the highest relative abundance, with an average Cq of 29, followed by miR-25-3p, miR-520f, let-7a-5p and miR-188-5p (Fig. 2a), where the order of abundance was matched to the corresponding miRNA concentrations determined by DDPCR (Fig. 2b). MiR-598 was an exception, however, where its absolute counts determined by DDPCR were markedly lower than expected (<5copies/μL plasma), in contrast to its relative expression levels determined by RT-qPCR with an average Cq of 35 (Fig. 2). The stability scores of the candidate miRNAs were determined using NormFinder and Bestkeeper which identified miR-25-3p, miR-188-5p, miR-222-3p and miR-520f as the miRNAs that exhibited low variability across the four groups analysed (Table 2), as well as high expression determined by qPCR and DDPCR. Although let-7a-5p was ranked the 4th most stable miRNA, with miR-25-3p being 5th, the expression levels of miR-25-3p were markedly higher than that of let-7a-5p expression levels in nanoString, qPCR and DDPCR analyses, in association with good stability values for miR-25-3p, as determined by NormFinder and Bestkeeper. Consistent and robust detection across platforms is the top selection criteria; therefore, miR-25-3p, miR-188-5p, miR-222-3p and miR-520f were selected for further analyses in individual plasma samples.
The expression levels of the four candidate miRNAs were analysed by qPCR in 40 individual plasma samples (10 plasma samples per group) and the variability of each miRNA across the plasma samples analysed determined by the stability scores calculated using NormFinder and BestKeeper. The relative abundance, determined by RT-qPCR, of miR-222-3p was the most highly expressed in all plasma samples with the lowest average Cq, followed by miR-25-3p, then miR-188-5p and lastly miR-520f (Fig. 3). These results were in agreement with that obtained using pooled plasma samples where the selected miRNAs exhibited the same order of abundance.
Stability analysis of the four candidate miRNAs in individual plasma samples using NormFinder generated stability scores of 0.115–0.601, indicating that all four miRNAs were good reference miRNAs, with miR-25-3p, miR-188-5p and miR-222-3p being the top three least variable miRNAs (in order of increasing variability) in the individual plasma samples tested (Table 3). Bestkeeper analysis showed that miR-188-5p was the most stable miRNA, followed by miR-222-3p and miR-25-3p, where the stability score for miR-25-3p was 1.30 and not in agreement with the level of variability determined by the NormFinder algorithm. As the stability analysis for miR-188-5p and miR-222-3p were consistent between NormFinder and Bestkeeper with stability scores in both cases indicating low variability, miR-188-5p and miR-222-3p were chosen as the two endogenous reference miRNAs for plasma miRNA quantitation.
Reference housekeeping panel is effective for normalising RT-qPCR data for circulating miRNAs
Many oestrogen-responsive miRNAs have been identified, however the studies were performed in breast cancer cell lines and other cell types [22,23,24]. The effects of high oestrogen levels on miRNAs in circulation have not been investigated, and there is also little information regarding gender associated miRNA expression profiles. As such, in order to validate the housekeeping miRNA panel, we determined the relative expression levels of a pregnancy-associated miRNA, miR-141-3p in total RNA extracted from individual plasma samples from each group using RT-qPCR. Levels of miR-141-3p have been reported in maternal serum and plasma samples to be increased with gestational age [11, 25], and significantly decreased following termination of pregnancy [25, 26]. The raw counts for miR-141-3p obtained from nanoString nCounter® analysis here ranged from 15 to 35 counts in the samples analysed and did not show significant difference between the groups by pairwise comparisons.
The relative expression of miR-141-3p in individual plasma samples between the four groups was analysed using RT-qPCR with various normalisation strategies (Fig. 4). When normalised to spike-ins alone, levels of miR-141-3p were significantly higher in pregnant samples compared to males, non-pregnant and oral contraceptive (p < 0.05) as expected (Fig. 4a). Overall there was a significant difference between the means in the four groups (p = 0.031). Individually, miR-141-3p expression in each cohort was lower than the pregnant group, namely healthy males (p = 0.031), oral contraceptives (p = 0.013) and non-pregnant women (p = 0.0091). There was no significant difference between the pairwise comparisons of miR-141-3p expression between healthy males, non-pregnant females not on oral contraceptives and non-pregnant females taking oral contraceptives with pregnant females (p = 0.87). In addition, as a group miR-141-3p expression was also significantly lower than that of the pregnant group (p = 0.0030).
When the RT-qPCR data was normalised to both the synthetic spike-ins and endogenous miRNA controls miR-188-5p and miR-222-3p, the relative levels of miR-141-3p was also significantly increased in pregnant samples compared to non-pregnant and oral contraceptive samples (p < 0.05), similar to results obtained when the data was normalised to spike-ins alone. The significance levels between cohorts, however, increased with the inclusion of the endogenous reference miRNAs in the normalisation (Fig. 4b). Overall there was a significant difference between the means in the four groups (p = 0.0166). Individually, miR-141-3p expression in each cohort is lower than that in the pregnant group, namely healthy males (p = 0.027), oral contraceptive women (p = 0.0057) and non-pregnant women (p = 0.0058). There was no significant difference in miR-141-3p expression between the healthy male, oral contraceptive women and non-pregnant women groups (p = 0.77), but as a group, their miR-141-3p levels were significantly lower than the pregnant women group (p = 0.0016).
Mounting evidence highlight the value of circulating miRNAs as a source of stable non-invasive biomarkers, with unique circulating microRNA profiles defined for various pathological conditions [27,28,29,30,31]. However, the lack of universal housekeeping controls and standardised methods for quantitating circulating miRNAs remains a major hurdle for the accurate analysis of miRNA expression in biofluids. This caveat also hinders translation of disease-associated miRNA signatures characterised using high throughput methods such as microarrays and next generation sequencing methods into diagnostic tests in the laboratory. As qPCR-based methods are relatively affordable and accessible for most laboratories, it remains the gold standard for miRNA analysis. Therefore, the identification of a robust endogenous miRNA control(s) specific for each condition with stable and consistent level of expression across platforms is fundamentally important.
Precedence exists for utilising short RNAs as endogenous gene normalisation factor(s), albeit with certain pitfalls. Small nucleolar RNAs (snoRNA) such as RNU6B, which is highly expressed in cells , and miR-16 which is abundantly expressed in blood are frequently used normalisation factors for miRNA studies [33,34,35,36], but their suitability as endogenous controls is controversial. RNU6B is present at low concentrations or shown to exhibit wide concentration variations in circulation , and was reported to be less resistant to multiple freeze–thaw cycles compared with other mature miRNA species analysed in serum samples . Red blood cells contain high levels of miR-16 , and studies have demonstrated that expression of miR-16 was significantly increased in haemolysed plasma samples, with the miR-16 levels proportional to the degree of haemolysis [40, 41], indicating that miR-16 is more useful as a haemolysis marker than a normalisation factors. Furthermore, miR-16 has been reported to be responsive to steroid hormones, downregulated by progestin and oestradiol in breast cancer cells [42, 43], and serum miR-16 is part of a diagnostic/prognostic miRNA panel for ovarian cancer patients . Taken together, there exists a need for a stable and non-hormone responsive miRNA in normalising endogenous miRNA expression levels in human plasma.
As an initial screen of 11 candidate reference miRNAs selected from the nanoString nCounter® assay, we analysed the expression of those candidate miRNAs in pooled plasma samples of (1) males, (2) non-pregnant females not on oral contraceptives, (3) non-pregnant oral contraceptive users and (4) pregnant females. Out of the 11 candidate reference miRNAs, 5 of them were detected at low levels (Cq > 38) or not amplified by qPCR or DDPCR, demonstrating marked discordance between the detection methods. This observation highlighted platform variations between nanoString and qPCR-based methods. The results reported here showed that many of the candidate miRNAs which exhibited high counts in nanoString did not correspondingly demonstrate high abundance when analysed using RT-qPCR and DDPCR, and low nanoString counts may not necessarily mean poor detectability via PCR amplification methods. nCounter® is a hybridisation based technology with high specificity without an amplification step; eliminating amplification bias introduced in qPCR methods. Low sensitivity, however, is a major limitation of nanoString, especially for the analysis of low input RNA samples such as serum/plasma RNA, which can be overcome by qPCR techniques. As such, it may not be surprising that the low concordance between nanoString and qPCR-based results observed in this present study has also been reported by other groups looking at circulating miRNAs. A study evaluating the performance of microRNA array platforms using different breast cancer cell lines found that hybridisation-based technology such as the nCounter® system demonstrated lower sensitivity compared to next generation sequencing and PCR-based platforms . Mestdagh et al.  conducted a large scale analysis of commercially available quantitative miRNA expression profiling platforms and evaluated them based on reproducibility, sensitivity, accuracy, specificity and concordance of differential expression. Their results demonstrated substantial inter-platform differences associated with each type of technology coupled with an inverse correlation between sensitivity and specificity . In this present study, the nanoString nCounter® miRNA assay v2 was useful as an initial screening tool to identify highly expressed miRNAs but users should be aware of false negatives due to the reduced assay sensitivity. Regardless, the results obtained in this study using nanoString led to the identification of miR-188-5p and miR-222-3p from human plasma samples, which exhibited consistently stable levels of expression between the cohorts analysed, across two measurement platforms.
During RNA extraction, there is significant protein and lipid variability between individual plasma or serum samples, which may affect RNA extraction efficiency and introduce potential PCR inhibitors [33, 47], therefore, a synthetic spike-in control is necessary. However, synthetic spike-in miRNAs only correct for differences in extraction efficiency and are not useful in normalising endogenous miRNA content. This may be sufficient for studies that have standardised pre-analytical conditions including collection, plasma preparation, sample storage and total RNA extraction methods. Hence it is prudent that suitable reference miRNAs are established for each study with appropriate test samples to avoid introducing bias in the analysis and skewing the final results. In this study, we used the nanoString nCounter® miRNA assay to obtain the plasma miRNA profiles of non-pregnant women, women who are currently taking oral contraceptives and pregnant females, mixed groups that represent individuals with low and high circulating oestradiol concentrations, respectively. From this data, we identified two stable endogenous reference miRNAs that are suitable for the normalisation of plasma miRNA expression in addition to synthetic spike in miRNAs, such as cel-miR-39 and ath-miR159a, which will allow for more accurate qPCR normalisation for studying oestrogen responsive miRNAs that may be involved in regulating other haemostatic factors leading to Protein S deficiency during pregnancy.
The independent plasma analyses of four candidate miRNAs, miR-25-3p, miR-188-5p, miR-222-5p and miR-520f showed that miR-520f expression exhibited high inter-individual variability, with Cq > 38 for many samples tested. The stability analyses using NormFinder and Bestkeeper showed that the remaining three miRNAs were stable between individual plasma samples tested but based on the stability scores, miR-188-5p and miR-222-3p were selected as the reference miRNAs. Both miR-188-5p and miR-222-3p have been previously evaluated as housekeeping miRNAs in plasma samples from cancer patients  and in breast tumour tissue biopsies , respectively, but were not ultimately chosen as reference miRNAs due to variability between the cancer cohorts. This indicates that miR-188-5p and miR-222-3p were not suitable for cancer studies, but it validates the data in this study that miR-188-5p and miR-222-3p exhibited and robust expression in multiple tissue types.
Furthermore, it should be noted that while miR-188-5p and miR-222-3p were found to be stable in the cohorts analysed in this study, they do have ascribed biological functions. In hepatocellular carcinoma cells, miR-188-5p was shown to inhibit growth of hepatocellular carcinoma cells by targeting the fibroblast growth factor  and was also demonstrated to be growth inhibitory in prostate cancer cells . In contrast, miR-222-3p has been reported to be important in the epithelial mesenchymal transition in breast cancer cell lines and its expression associated with ERK1/2 signalling activity . Expression of miR-222-3p has also been shown to be upregulated by selective estrogen receptor downregulator (SERD), fulvestrant in the breast cancer cell line, MCF-7 . In circulation, Vickers et al.  have found miR-188-5p to be highly abundant in high density lipoproteins and circulating levels of miR-222-3p have been demonstrated to be a useful diagnostic marker in gastric cancer; and miR-222-3p has been reported to negatively regulate the expression of ERα in endometrial carcinoma cells , but oestrogen-mediated regulation of miR-222-3p and miR-188-5p have not been previously reported.
RT-qPCR analyses of miR-141-3p expression was chosen to adjudicate the novel miRNA normalisation panel established in this study. Levels of serum and plasma miR-141-3p have been reported by multiple studies to be significantly upregulated in serum and plasma samples of pregnant women using PCR-based miRNA arrays [11, 56], and as such, serves as an ideal candidate for validating the reference normalisation panel. The results presented here showed that significant changes in miR-141-3p expression between non-pregnant (with or without oral contraceptives) and pregnant females can be observed when the RT-qPCR data was normalised to the synthetic spike-ins. This is consistent with previously published data, despite a different normalisation strategy being applied in each study. Here, we demonstrated that greater levels of significance in the differences in expression levels were observed when miR-188-5p and miR-222-3p were introduced into the normalisation panel with the spike-ins, compared to normalisation with the spike-ins alone. Importantly, this suggests that small, but significant differences in miRNA expression initially overlooked when results were normalised to synthetic spike-ins alone will be identified by including endogenous miRNA normalisation factor(s). Circulating miR-141-3p has been previously investigated as a marker for metastatic disease in prostate cancer patients  and also in colon cancer patients . A more recent study has identified miR-141-3p to be a novel biomarker for metastatic colorectal cancer in plasma samples from male and female patients [58, 59]. However, miR-141-3p levels were either only determined in men or there was not a comparison of relative miR-141-3p expression between male and female patient cohorts in these studies.
To date, there are no reports describing reference miRNAs for the analysis of hormone-responsive miRNAs in human serum and plasma samples. Furthermore, studies investigating oestrogen-responsive miRNAs have predominantly been performed in human breast cancer cell lines [22, 23, 60, 61]; whereas, hormone induced miRNAs in thrombosis remains poorly characterised. Thus the accurate identification of miRNAs regulated during pregnancy under high levels of circulating oestrogen is critical in assessing which target haemostatic genes could heighten thrombotic risks. It is vital, therefore, that suitable endogenous miRNAs are established in the appropriate study cohorts. In addition, there have been many studies which characterised miRNA expression profiles associated with pregnancy complications such as preeclampsia [62, 63], ectopic pregnancies  and gestational diabetes . Baseline levels, however, are absent as these studies compared expression profiles between pregnant women only and not against non-pregnant participants. This current study is the first report for the systematic identification of reference miRNAs for the normalisation of miRNA qPCR data in studies investigating oestrogen-responsive miRNAs.
While the miRNA reference panel reported in this study has been established for identifying novel miRNAs that contribute to increased thrombotic risk in pregnant women with high circulating oestrogen levels compared to non-pregnant women and women taking oral contraceptives, it can be applied in a number of other biological contexts. The control miRNAs were validated to be stable between healthy men, non-pregnant women, women on oral contraceptive and pregnant women. As such, it may be used to normalise miRNA expression in studies investigating gender-specific associations of miRNA expression in disease. A recent study conducted a comprehensive comparison of miRNA expression profiles in healthy and tumour tissues between men and women demonstrated marked gender associated differences in miRNA/isomiR expression . Reports of associations of gender in a population based study  and influence of gender differences on microRNA gene regulation in disease  all support the need for the identification of gender neutral miRNAs that can be used in studies that compare miRNA expression between male and female patient cohorts. The miRNA reference panel in this report will also be useful in the comparison of plasma miRNA expression regulated by synthetic and natural oestrogens in healthy individuals, with the caveat that the stability of the reference miRNAs requires validation in breast cancer patients or women with polycystic ovarian syndrome. Furthermore, the less than ideal inter-platform comparability of miRNA expression between currently available commercial platforms reported by multiple studies should be taken into account [45, 46, 68, 69]. As such, the strategy for establishing suitable reference miRNAs in any biological context should employ a different detection platform for screening and validation, i.e. reverse transcription quantitative PCR versus hybridisation.
Endogenous controls for circulating miRNAs are not universal and the choice of normalising miRNAs must be carefully selected for each study. If a previously reported endogenous control miRNA is absent, it is pertinent that suitable housekeeping miRNAs be established in the study cohorts using our reported strategy. Also, the variability in lipid content of biological fluids between individuals can affect RNA extraction efficiencies, as such, the addition of synthetic spike-ins cel-miR-39 and ath-miR159a are necessary as first line controls for normalising variations in extraction efficiency; however, synthetic spike-ins alone, are not adequate for normalisation. The results presented here identified two plasma miRNAs, miR-188-5p and miR-222-3p, whose expression levels were stable between the four cohorts analysed, representing individuals with low and high circulating endogenous oestradiol concentrations, as well as high circulating synthetic oestrogens, and should be included in the normalisation panel for the identification of oestrogen responsive miRNAs in circulation. The four reference miRNA panel established in this study may also be useful for normalising qPCR data comparing miRNA expression between men and women, non-pregnant and pregnant females, and also the potential effects of endogenous and synthetic oestrogens on plasma miRNA expression.
3 prime untranslated region
real time quantitative polymerase chain reaction
Wightman B, Ha I, Ruvkun G. Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans. Cell. 1993;75(5):855–62.
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(5):843–54.
Li Y, Kowdley KV. MicroRNAs in common human diseases. Genom Proteom Bioinform. 2012;10(5):246–53.
Nishiguchi T, Imanishi T, Akasaka T. MicroRNAs and cardiovascular diseases. Biomed Res Int. 2015;2015:682857.
Bhatt K, Kato M, Natarajan R. Mini-review: emerging roles of microRNAs in the pathophysiology of renal diseases. Am J Physiol Renal Physiol. 2016;310(2):F109–18.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al. MicroRNA expression profiles classify human cancers. Nature. 2005;435(7043):834–8.
Peiro-Chova L, Pena-Chilet M, Lopez-Guerrero JA, Garcia-Gimenez JL, Alonso-Yuste E, Burgues O, Lluch A, Ferrer-Lozano J, Ribas G. High stability of microRNAs in tissue samples of compromised quality. Virchows Archiv. 2013;463(6):765–74.
Chen X, Ba Y, Ma L, Cai X, Yin Y, Wang K, Guo J, Zhang Y, Chen J, Guo X, et al. Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008;18(10):997–1006.
Cortez MA, Bueso-Ramos C, Ferdin J, Lopez-Berestein G, Sood AK, Calin GA. MicroRNAs in body fluids—the mix of hormones and biomarkers. Nat Rev Clin Oncol. 2011;8(8):467–77.
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 USA. 2008;105(30):10513–8.
Gilad S, Meiri E, Yogev Y, Benjamin S, Lebanony D, Yerushalmi N, Benjamin H, Kushnir M, Cholakh H, Melamed N, et al. Serum microRNAs are promising novel biomarkers. PLoS ONE. 2008;3(9):e3148.
Freedman JE, Ercan B, Morin KM, Liu CT, Tamer L, Ayaz L, Kanadasi M, Cicek D, Seyhan AI, Akilli RE, et al. The distribution of circulating microRNA and their relation to coronary disease. F1000Research. 2012;1:50.
Chen SY, Wang Y, Telen MJ, Chi JT. The genomic analysis of erythrocyte microRNA expression in sickle cell diseases. PLoS ONE. 2008;3(6):e2360.
Landry P, Plante I, Ouellet DL, Perron MP, Rousseau G, Provost P. Existence of a microRNA pathway in anucleate platelets. Nat Struct Mol Biol. 2009;16(9):961–6.
Zernecke A, Bidzhekov K, Noels H, Shagdarsuren E, Gan L, Denecke B, Hristov M, Koppel T, Jahantigh MN, Lutgens E, et al. Delivery of microRNA-126 by apoptotic bodies induces CXCL12-dependent vascular protection. Sci Signal. 2009;2(100):ra81.
Diehl P, Fricke A, Sander L, Stamm J, Bassler N, Htun N, Ziemann M, Helbing T, El-Osta A, Jowett JB, et al. Microparticles: major transport vehicles for distinct microRNAs in circulation. Cardiovasc Res. 2012;93(4):633–44.
Tay JW, Romeo G, Hughes QW, Baker RI. Micro-ribonucleic Acid 494 regulation of protein S expression. J Thromb Haemost. 2013;11(8):1547–55.
Mackie I, Cooper P, Lawrie A, Kitchen S, Gray E, Laffan M. British committee for standards in H: guidelines on the laboratory aspects of assays used in haemostasis and thrombosis. Int J Lab Hematol. 2013;35(1):1–13.
Wedderburn RWM. Quasi-likelihood functions, generalized linear models, and the Gauss–Newton method. Biometrika. 1974;61(3):439–47.
Andersen CL, Jensen JL, Orntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64(15):5245–50.
Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: Bestkeeper–Excel-based tool using pair-wise correlations. Biotech Lett. 2004;26(6):509–15.
Ferraro L, Ravo M, Nassa G, Tarallo R, De Filippo MR, Giurato G, Cirillo F, Stellato C, Silvestro S, Cantarella C, et al. Effects of oestrogen on microRNA expression in hormone-responsive breast cancer cells. Horm Cancer. 2012;3(3):65–78.
Bhat-Nakshatri P, Wang G, Collins NR, Thomson MJ, Geistlinger TR, Carroll JS, Brown M, Hammond S, Srour EF, Liu Y, et al. Estradiol-regulated microRNAs control estradiol response in breast cancer cells. Nucleic Acids Res. 2009;37(14):4850–61.
He YQ, Sheng JQ, Ling XL, Fu L, Jin P, Yen L, Rao J. Estradiol regulates miR-135b and mismatch repair gene expressions via estrogen receptor-beta in colorectal cells. Exp Mol Med. 2012;44(12):723–32.
Chim SS, Shing TK, Hung EC, Leung TY, Lau TK, Chiu RW, Lo YM. Detection and characterization of placental microRNAs in maternal plasma. Clin Chem. 2008;54(3):482–90.
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(11):1767–71.
Zuo Z, Maiti S, Hu S, Loghavi S, Calin GA, Garcia-Manero G, Kantarjian HM, Medeiros LJ, Cooper LJ, Bueso-Ramos CE. Plasma circulating-microRNA profiles are useful for assessing prognosis in patients with cytogenetically normal myelodysplastic syndromes. Mod Pathol. 2015;28(3):373–82.
Polytarchou C, Oikonomopoulos A, Mahurkar S, Touroutoglou A, Koukos G, Hommes DW, Iliopoulos D. Assessment of circulating MicroRNAs for the diagnosis and disease activity evaluation in patients with ulcerative colitis by using the nanostring technology. Inflamm Bowel Dis. 2015;21(11):2533–9.
Lian F, Cui Y, Zhou C, Gao K, Wu L. Identification of a plasma four-microRNA panel as potential noninvasive biomarker for osteosarcoma. PLoS ONE. 2015;10(3):e0121499.
Fornari F, Ferracin M, Trere D, Milazzo M, Marinelli S, Galassi M, Venerandi L, Pollutri D, Patrizi C, Borghi A, et al. Circulating microRNAs, miR-939, miR-595, miR-519d and miR-494, identify cirrhotic patients with HCC. PLoS ONE. 2015;10(10):e0141448.
Faruq O, Vecchione A. microRNA: diagnostic perspective. Front Med. 2015;2:51.
Iio A, Takagi T, Miki K, Naoe T, Nakayama A, Akao Y. DDX6 post-transcriptionally down-regulates miR-143/145 expression through host gene NCR143/145 in cancer cells. Biochim Biophys Acta. 2013;1829(10):1102–10.
Kroh EM, Parkin RK, Mitchell PS, Tewari M. Analysis of circulating microRNA biomarkers in plasma and serum using quantitative reverse transcription-PCR (qRT-PCR). Methods. 2010;50(4):298–301.
Song J, Bai Z, Han W, Zhang J, Meng H, Bi J, Ma X, Han S, Zhang Z. Identification of suitable reference genes for qPCR analysis of serum microRNA in gastric cancer patients. Dig Dis Sci. 2012;57(4):897–904.
Wang J, Chen J, Chang P, LeBlanc A, Li D, Abbruzzesse JL, Frazier ML, Killary AM, Sen S. MicroRNAs in plasma of pancreatic ductal adenocarcinoma patients as novel blood-based biomarkers of disease. Cancer Prev Res. 2009;2(9):807–13.
Redova M, Poprach A, Nekvindova J, Iliev R, Radova L, Lakomy R, Svoboda M, Vyzula R, Slaby O. Circulating miR-378 and miR-451 in serum are potential biomarkers for renal cell carcinoma. J Transl Med. 2012;10:55.
Wang K, Yuan Y, Cho JH, McClarty S, Baxter D, Galas DJ. Comparing the MicroRNA spectrum between serum and plasma. PLoS ONE. 2012;7(7):e41561.
Xiang M, Zeng Y, Yang R, Xu H, Chen Z, Zhong J, Xie H, Xu Y, Zeng X. U6 is not a suitable endogenous control for the quantification of circulating microRNAs. Biochem Biophys Res Commun. 2014;454(1):210–4.
Merkerova M, Belickova M, Bruchova H. Differential expression of microRNAs in hematopoietic cell lineages. Eur J Haematol. 2008;81(4):304–10.
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. 2012;5(3):492–7.
Kirschner MB, Kao SC, Edelman JJ, Armstrong NJ, Vallely MP, van Zandwijk N, Reid G. Haemolysis during sample preparation alters microRNA content of plasma. PLoS ONE. 2011;6(9):e24145.
Rivas MA, Venturutti L, Huang YW, Schillaci R, Huang TH, Elizalde PV. Downregulation of the tumor-suppressor miR-16 via progestin-mediated oncogenic signaling contributes to breast cancer development. Breast Cancer Res. 2012;14(3):R77.
Yu X, Zhang X, Dhakal IB, Beggs M, Kadlubar S, Luo D. Induction of cell proliferation and survival genes by estradiol-repressed microRNAs in breast cancer cells. BMC Cancer. 2012;12:29.
Meng X, Joosse SA, Muller V, Trillsch F, Milde-Langosch K, Mahner S, Geffken M, Pantel K, Schwarzenbach H. Diagnostic and prognostic potential of serum miR-7, miR-16, miR-25, miR-93, miR-182, miR-376a and miR-429 in ovarian cancer patients. Br J Cancer. 2015;113(9):1358–66.
Knutsen E, Fiskaa T, Ursvik A, Jorgensen TE, Perander M, Lund E, Seternes OM, Johansen SD, Andreassen M. Performance comparison of digital microRNA profiling technologies applied on human breast cancer cell lines. PLoS ONE. 2013;8(10):e75813.
Mestdagh P, Hartmann N, Baeriswyl L, Andreasen D, Bernard N, Chen C, Cheo D, D’Andrade P, DeMayo M, Dennis L, et al. Evaluation of quantitative miRNA expression platforms in the microRNA quality control (miRQC) study. Nat Methods. 2014;11(8):809–15.
Sourvinou IS, Markou A, Lianidou ES. Quantification of circulating miRNAs in plasma: effect of preanalytical and analytical parameters on their isolation and stability. J Mol Diagn. 2013;15(6):827–34.
Hu J, Wang Z, Liao BY, Yu L, Gao X, Lu S, Wang S, Dai Z, Zhang X, Chen Q, et al. Human miR-1228 as a stable endogenous control for the quantification of circulating microRNAs in cancer patients. Int J Cancer. 2014;135(5):1187–94.
Rinnerthaler G, Hackl H, Gampenrieder SP, Hamacher F, Hufnagl C, Hauser-Kronberger C, Zehentmayr F, Fastner G, Sedlmayer F, Mlineritsch B, et al. miR-16-5p Is a stably-expressed housekeeping microrna in breast cancer tissues from primary tumors and from metastatic sites. Int J Mol Sci. 2016;17(2):156.
Fang F, Chang RM, Yu L, Lei X, Xiao S, Yang H, Yang LY. MicroRNA-188-5p suppresses tumor cell proliferation and metastasis by directly targeting FGF5 in hepatocellular carcinoma. J Hepatol. 2015;63(4):874–85.
Zhang H, Qi S, Zhang T, Wang A, Liu R, Guo J, Wang Y, Xu Y. miR-188-5p inhibits tumour growth and metastasis in prostate cancer by repressing LAPTM4B expression. Oncotarget. 2015;6(8):6092–104.
Stinson S, Lackner MR, Adai AT, Yu N, Kim HJ, O’Brien C, Spoerke J, Jhunjhunwala S, Boyd Z, Januario T, et al. TRPS1 targeting by miR-221/222 promotes the epithelial-to-mesenchymal transition in breast cancer. Sci Signal. 2011;4(177):ra41.
Rao X, Di Leva G, Li M, Fang F, Devlin C, Hartman-Frey C, Burow ME, Ivan M, Croce CM, Nephew KP. MicroRNA-221/222 confers breast cancer fulvestrant resistance by regulating multiple signaling pathways. Oncogene. 2011;30(9):1082–97.
Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011;13(4):423–33.
Liu B, Che Q, Qiu H, Bao W, Chen X, Lu W, Li B, Wan X. Elevated MiR-222-3p promotes proliferation and invasion of endometrial carcinoma via targeting ERalpha. PLoS ONE. 2014;9(1):e87563.
Miura K, Higashijima A, Mishima H, Miura S, Kitajima M, Kaneuchi M, Yoshiura K, Masuzaki H. Pregnancy-associated microRNAs in plasma as potential molecular markers of ectopic pregnancy. Fertil Steril. 2015;103(5):1202–8.
Yaman Agaoglu F, Kovancilar M, Dizdar Y, Darendeliler E, Holdenrieder S, Dalay N, Gezer U. Investigation of miR-21, miR-141, and miR-221 in blood circulation of patients with prostate cancer. Tumour Biol. 2011;32(3):583–8.
Cheng H, Zhang L, Cogdell DE, Zheng H, Schetter AJ, Nykter M, Harris CC, Chen K, Hamilton SR, Zhang W. Circulating plasma MiR-141 is a novel biomarker for metastatic colon cancer and predicts poor prognosis. PLoS ONE. 2011;6(3):e17745.
Sun Y, Liu Y, Cogdell D, Calin GA, Sun B, Kopetz S, Hamilton SR, Zhang W. Examining plasma microRNA markers for colorectal cancer at different stages. Oncotarget. 2016;7(10):11434–49.
Maillot G, Lacroix-Triki M, Pierredon S, Gratadou L, Schmidt S, Benes V, Roche H, Dalenc F, Auboeuf D, Millevoi S, et al. Widespread estrogen-dependent repression of microRNAs involved in breast tumor cell growth. Cancer Res. 2009;69(21):8332–40.
Wickramasinghe NS, Manavalan TT, Dougherty SM, Riggs KA, Li Y, Klinge CM. Estradiol downregulates miR-21 expression and increases miR-21 target gene expression in MCF-7 breast cancer cells. Nucleic Acids Res. 2009;37(8):2584–95.
Gunel T, Zeybek YG, Akcakaya P, Kalelioglu I, Benian A, Ermis H, Aydinli K. Serum microRNA expression in pregnancies with preeclampsia. Genet Mol Res. 2011;10(4):4034–40.
Yang Q, Lu J, Wang S, Li H, Ge Q, Lu Z. Application of next-generation sequencing technology to profile the circulating microRNAs in the serum of preeclampsia versus normal pregnant women. Clin Chim Acta. 2011;412(23–24):2167–73.
Zhu Y, Tian F, Li H, Zhou Y, Lu J, Ge Q. Profiling maternal plasma microRNA expression in early pregnancy to predict gestational diabetes mellitus. Int J Gynaecol Obstet. 2015;130(1):49–53.
Guo L, Liang T, Yu J, Zou Q. A comprehensive analysis of miRNA/isomiR expression with gender difference. PLoS ONE. 2016;11(5):e0154955.
Ameling S, Kacprowski T, Chilukoti RK, Malsch C, Liebscher V, Suhre K, Pietzner M, Friedrich N, Homuth G, Hammer E, et al. Associations of circulating plasma microRNAs with age, body mass index and sex in a population-based study. BMC Med Genom. 2015;8:61.
Sharma S, Eghbali M. Influence of sex differences on microRNA gene regulation in disease. Biol sex differ. 2014;5(1):3.
Pradervand S, Weber J, Lemoine F, Consales F, Paillusson A, Dupasquier M, Thomas J, Richter H, Kaessmann H, Beaudoing E, et al. Concordance among digital gene expression, microarrays, and qPCR when measuring differential expression of microRNAs. Biotechniques. 2010;48(3):219–22.
Sato F, Tsuchiya S, Terasawa K, Tsujimoto G. Intra-platform repeatability and inter-platform comparability of microRNA microarray technology. PLoS ONE. 2009;4(5):e5540.
JWT: Wrote the manuscript, conceptualised and performed the experiments, analysed the data, corresponding author. IJ: Performed statistical analyses on Nanostring data, evaluation of data and edited the manuscript. QWH: Evaluation of data and edited the manuscript. JT: Evaluation of data and edited the manuscript. RIB: Evaluation of data and edited the manuscript. All authors read and approved the final manuscript.
Many thanks to Dr Abha Chopra for technical advice regarding sample preparation for nanoString nCounter® analyses and Dr Gisela Mir Arnau at the Molecular Genomics Core Facility, Peter MacCallum Cancer Centre for running total RNA samples on the nCounter®. Acknowledgments also extend to Dr Anthony Beckhouse at Bio-Rad Laboratories for extensive technical support in the optimisation of miRNA digital droplet PCR experiments.
The authors declare that they have no competing interests.
Availability of data and materials
Nanostring raw counts and pairwise comparisons are available for download from Additional Material.
Consent to publish
All participants provided written informed consent as part of their consent for participating in this study for the use of the results obtained for scientific publication.
Ethics approval and consent to participate
All participants provided written informed consent for the use of venous blood samples in this study prior to blood collection. This project was approved by the Human Research Ethics Committee at Murdoch University (HREC Ref: 2013/131).
This work has been funded by Murdoch University, the Perth Blood Institute and the Hollywood Private Hospital Research Foundation.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Tay, J.W., James, I., Hughes, Q.W. et al. Identification of reference miRNAs in plasma useful for the study of oestrogen-responsive miRNAs associated with acquired Protein S deficiency in pregnancy. BMC Res Notes 10, 312 (2017). https://doi.org/10.1186/s13104-017-2636-3
- Reference miRNAs
- Droplet Digital PCR (DDPCR)
- Individual Plasma Samples