Real-time quantitative PCR and droplet digital PCR for plant miRNAs in mammalian blood provide little evidence for general uptake of dietary miRNAs: limited evidence for general uptake of dietary plant xenomiRs

Kenneth W Witwer, Melissa A McAlexander, Suzanne E Queen, Robert J Adams, Kenneth W Witwer, Melissa A McAlexander, Suzanne E Queen, Robert J Adams

Abstract

Evidence that exogenous dietary miRNAs enter the bloodstream and tissues of ingesting animals has been accompanied by an indication that at least one plant miRNA, miR168, participates in "cross-kingdom" regulation of a mammalian transcript. If confirmed, these findings would support investigation of miRNA-based dietary interventions in disease. Here, blood was obtained pre- and post-prandially (1, 4, 12 h) from pigtailed macaques that received a miRNA-rich plant-based substance. Plant and endogenous miRNAs were measured by RT-qPCR. Although low-level amplification was observed for some plant miRNA assays, amplification was variable and possibly non-specific, as suggested by droplet digital PCR. A consistent response to dietary intake was not observed. While our results do not support general and consistent uptake of dietary plant miRNAs, additional studies are needed to establish whether or not plant or animal xenomiRs are transferred across the gut in sufficient quantity to regulate endogenous genes.

Keywords: RT-qPCR; diet; digital PCR; exosome; extracellular vesicle; microRNA; nutrition; plant; xenomiR.

Figures

https://www.ncbi.nlm.nih.gov/pmc/articles/instance/3849155/bin/rna-10-1080-g1.jpg
Figure 1. miRNA abundance and amplification efficiency in a plant-based dietary substance and in mammalian plasma. (A) Relative abundance of several highly conserved plant miRNAs in a soy- and fruit substance. RNA was purified from 100 µl “Silk” using the mirVana total RNA isolation protocol (Ambion/Life Technologies). Hydrolysis probe-based RT-qPCR measurements were made from up to three independent RNA isolations. Fold abundance was calculated relative to miR160. Similar results were obtained at different input RNA dilutions (not shown). Error bars represent standard deviation. (B and C) RT-qPCR assay linearity and amplification efficiency was assessed (shown here: miR160) using 2-fold serial dilutions of cDNA prepared with mature miRNA-specific stem-loop primers and plant RNA from the dietary substance (B) or RNA from pre- or post-prandial plasma samples (C). For plant RNA, (B) starting dilution (#1) corresponded to cDNA from approximately 6.3 nanoliters of the original substance. Final dilution (#6) was from about 200 picoliters original volume. Amplification occurred in all wells, and amplification efficiency was approximately 100%. Error bars are standard deviation of triplicate measurements (less than half a cycle for all dilutions). For plasma, dilutions corresponded to material from 2 µl plasma (#1) down to 60 nl plasma (#6). Standard deviation corresponded to half a cycle or greater for 8 of 12 samples, and several wells within technical replicate groups did not amplify. Amplification efficiencies were 94% (plasma 1) and 65% (plasma 2).
https://www.ncbi.nlm.nih.gov/pmc/articles/instance/3849155/bin/rna-10-1080-g2.jpg
Figure 2. Plant miRNAs in circulation before and following intake of a miRNA-rich food source. (A) Apparent abundance in plasma of miR166 relative to miR160, as assessed by RT-qPCR with RNA samples purified by mirVana technique, did not match relative abundance in the food source at any time point pre- or post-prandial (compare with Fig. 1A) in either of two macaque subjects designated “67X” and “HL26.” (B and C) Plasma miR160 and miR166 Cq values and standard deviations by time point and subject (67X and HL26). Missing standard deviation bars are due to insufficient numbers of amplifying wells. Note that, to the extent it could be estimated, direction of apparent regulation often differed between the two subjects.
https://www.ncbi.nlm.nih.gov/pmc/articles/instance/3849155/bin/rna-10-1080-g3.jpg
Figure 3. RT-qPCR results for RNA isolated with a high-performance biofluids-specific protocol (Exiqon). Cq and standard deviation of RT-qPCR results for five plant miRNAs and one exogenous spike-in RNA are shown for plasma from subjects 67X (left) and HL26 (right) obtained at the pre-prandial (0) time point and at 1, 4 and 12 h post-prandial. Asterisk indicates a time point for which no reactions amplified for that miRNA. In the bottom panels, highly consistent recovery is shown for a miRNA (cel-miR-39) that was spiked in to the RNA isolation reactions to assess variability in the isolation technique that could have affected results.
https://www.ncbi.nlm.nih.gov/pmc/articles/instance/3849155/bin/rna-10-1080-g4.jpg
Figure 4. Droplet digital PCR indicates specific amplification of plant miRNA from plant source but not from pre- or postprandial animal plasma. (A and B) Intensity plots of DX100 QuantaSoft results for miR167 and miR168: for RNA from the dietary plant material (left), positive droplets clustered at high and relatively consistent intensity (specific amplification). Droplets from four hour postprandial plasma samples (right) had low and varying intensity, especially for miR168, consistent with inefficient or spurious amplification. Triplicate reactions are shown. Note that the intensity scale is not necessarily equivalent between different plots. (C) miR160 showed the highest level of amplification from plasma for miRNAs we examined by RT-qPCR (see Fig. 2), but by ddPCR with pre- and post-prandial plasma samples, there was a wide range of droplet intensities. From plant material, droplet intensities clustered consistently as for the miRNAs in (A and B) (not shown). The horizontal line is an example of where a threshold might be drawn. (D) Low-level non-specific amplification for plant miRs166, 167 and 168 from one no-template control reaction each (water). (E) Copy number counts and confidence interval minimum and maximum for miR160 per microliter of reaction volume over a 1:4 dilution series of material corresponding to cDNA from input into a reverse transcription reaction that corresponded to RNA obtained from 6.25 nl (dilution 1), 1.6 nl (dilution 2) and 400 pl (dilution 3) of the “Silk” plant-based substance. (F) Counts and confidence intervals for undiluted material from plasma RNA isolations for the pre-intake and 4 h postprandial time points for miR160 for subjects 67X and HL26.

References

    1. Zhang L, Hou D, Chen X, Li D, Zhu L, Zhang Y, et al. Exogenous plant MIR168a specifically targets mammalian LDLRAP1: evidence of cross-kingdom regulation by microRNA. Cell Res. 2012;22:107–26. doi: 10.1038/cr.2011.158.
    1. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33. doi: 10.1016/j.cell.2009.01.002.
    1. Mendell JT, Olson EN. MicroRNAs in stress signaling and human disease. Cell. 2012;148:1172–87. doi: 10.1016/j.cell.2012.02.005.
    1. Valadi H, Ekström K, Bossios A, Sjöstrand M, Lee JJ, Lötvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007;9:654–9. doi: 10.1038/ncb1596.
    1. Pegtel DM, Cosmopoulos K, Thorley-Lawson DA, van Eijndhoven MA, Hopmans ES, Lindenberg JL, et al. Functional delivery of viral miRNAs via exosomes. Proc Natl Acad Sci USA. 2010;107:6328–33. doi: 10.1073/pnas.0914843107.
    1. Witwer KW, Buzas EI, Bemis LT, Bora A, Lasser C, Lotvall J, et al. Standardisation of sample collection, isolation and analysis methods in extracellular vesicle research: An ISEV position paper. Journal of Extracellular Vesicles. 2013 doi: 10.3402/jev.v2i0.20360. In press.
    1. Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci USA. 2011;108:5003–8. doi: 10.1073/pnas.1019055108.
    1. Turchinovich A, Weiz L, Langheinz A, Burwinkel B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011;39:7223–33. doi: 10.1093/nar/gkr254.
    1. Turchinovich A, Weiz L, Burwinkel B. Extracellular miRNAs: the mystery of their origin and function. Trends Biochem Sci. 2012;37:460–5. doi: 10.1016/j.tibs.2012.08.003.
    1. Witwer KW. XenomiRs and miRNA homeostasis in health and disease: evidence that diet and dietary miRNAs directly and indirectly influence circulating miRNA profiles. RNA Biol. 2012;9:1147–54. doi: 10.4161/rna.21619.
    1. Jiang M, Sang X, Hong Z. Beyond nutrients: food-derived microRNAs provide cross-kingdom regulation. Bioessays. 2012;34:280–4. doi: 10.1002/bies.201100181.
    1. Zhang Y, Wiggins BE, Lawrence C, Petrick J, Ivashuta S, Heck G. Analysis of plant-derived miRNAs in animal small RNA datasets. BMC Genomics. 2012;13:381. doi: 10.1186/1471-2164-13-381.
    1. Petrick JS, Brower-Toland B, Jackson AL, Kier LD. Safety assessment of food and feed from biotechnology-derived crops employing RNA-mediated gene regulation to achieve desired traits: A scientific review. Regul Toxicol Pharmacol. 2013;66:167–76. doi: 10.1016/j.yrtph.2013.03.008. [epub ahead of print]
    1. Wang K, Li H, Yuan Y, Etheridge A, Zhou Y, Huang D, et al. The complex exogenous RNA spectra in human plasma: an interface with human gut biota? PLoS One. 2012;7:e51009. doi: 10.1371/journal.pone.0051009.
    1. McAlexander MA, Phillips MJ, Witwer KW. Comparison of methods for miRNA extraction from plasma and quantitative recovery of RNA from plasma and cerebrospinal fluid. Front Genet. 2013;4:83. doi: 10.3389/fgene.2013.00083.
    1. Hindson BJ, Ness KD, Masquelier DA, Belgrader P, Heredia NJ, Makarewicz AJ, et al. High-throughput droplet digital PCR system for absolute quantitation of DNA copy number. Anal Chem. 2011;83:8604–10. doi: 10.1021/ac202028g.
    1. Watson AK, Witwer KW. Do platform-specific factors explain microRNA profiling disparities? Clin Chem. 2012;58:472–4, author reply 474-5. doi: 10.1373/clinchem.2011.175281.
    1. Mullokandov G, Baccarini A, Ruzo A, Jayaprakash AD, Tung N, Israelow B, et al. High-throughput assessment of microRNA activity and function using microRNA sensor and decoy libraries. Nat Methods. 2012;9:840–6. doi: 10.1038/nmeth.2078.
    1. Peng P, Chan SW, Shah GA, Jacobsen SE. Plant genetics: increased outcrossing in hothead mutants. Nature. 2006;443:E8–, discussion E8-9. doi: 10.1038/nature05251.
    1. Lolle SJ, Victor JL, Young JM, Pruitt RE. Genome-wide non-mendelian inheritance of extra-genomic information in Arabidopsis. Nature. 2005;434:505–9. doi: 10.1038/nature03380.
    1. Wei LQ, Yan LF, Wang T. Deep sequencing on genome-wide scale reveals the unique composition and expression patterns of microRNAs in developing pollen of Oryza sativa. Genome Biol. 2011;12:R53. doi: 10.1186/gb-2011-12-6-r53.
    1. Hirschi KD. New foods for thought. Trends Plant Sci. 2012;17:123–5. doi: 10.1016/j.tplants.2012.01.004.
    1. Kosaka N, Ochiya T. Unraveling the Mystery of Cancer by Secretory microRNA: Horizontal microRNA Transfer between Living Cells. Front Genet. 2011;2:97.
    1. Semenov DV, Baryakin DN, Brenner EV, Kurilshikov AM, Vasiliev GV, Bryzgalov LA, et al. Unbiased approach to profile the variety of small non-coding RNA of human blood plasma with massively parallel sequencing technology. Expert Opin Biol Ther. 2012;12(Suppl 1):S43–51. doi: 10.1517/14712598.2012.679653.
    1. Witwer KW. Data submission and quality in microarray-based microRNA profiling. Clin Chem. 2013;59:392–400. doi: 10.1373/clinchem.2012.193813.
    1. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:H0034. doi: 10.1186/gb-2002-3-7-research0034.
    1. Meyer SU, Kaiser S, Wagner C, Thirion C, Pfaffl MW. Profound effect of profiling platform and normalization strategy on detection of differentially expressed microRNAs--a comparative study. PLoS One. 2012;7:e38946. doi: 10.1371/journal.pone.0038946.
    1. Witwer KW, Sarbanes SL, Liu J, Clements JE. A plasma microRNA signature of acute lentiviral infection: biomarkers of CNS disease. AIDS. 2011;204:1104–14.
    1. Liang H, Zen K, Zhang J, Zhang CY, Chen X. New roles for microRNAs in cross-species communication. RNA Biol. 2013;10 doi: 10.4161/rna.23663. [epub ahead of print.] In press.
    1. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci USA. 2008;105:10513–8. doi: 10.1073/pnas.0804549105.
    1. Sisk JM, Clements JE, Witwer KW. miRNA profiles of monocyte-lineage cells are consistent with complicated roles in HIV-1 restriction. Viruses. 2012;4:1844–64. doi: 10.3390/v4101844.
    1. Witwer KW, Watson AK, Blankson JN, Clements JE. Relationships of PBMC microRNA expression, plasma viral load, and CD4+ T-cell count in HIV-1-infected elite suppressors and viremic patients. Retrovirology. 2012;9:5. doi: 10.1186/1742-4690-9-5.
    1. Snow JW, Hale AE, Isaccs SK, Baggish AL, Chan SY. Ineffective delivery of diet-derived microRNAs to recipient animal organisms. RNA Biology. 2013;10 doi: 10.4161/rna.25246. In press.

Source: PubMed

3
Iratkozz fel