Effects of Cardiac Sympathetic Neurodegeneration and PPAR<i>γ</i> Activation on Rhesus Macaque Whole Blood miRNA and mRNA Expression Profiles

Metzger, Jeanette M.; Lopez, Mary S.; Schmidt, Jenna K.; Murphy, Megan E.; Vemuganti, Raghu; Emborg, Marina E.

doi:https://doi.org/10.1155/2020/9426204

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results and Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 9426204 | https://doi.org/10.1155/2020/9426204

Effects of Cardiac Sympathetic Neurodegeneration and PPARγ Activation on Rhesus Macaque Whole Blood miRNA and mRNA Expression Profiles

Jeanette M. Metzger,^1,2Mary S. Lopez,^2,3Jenna K. Schmidt,¹Megan E. Murphy,¹Raghu Vemuganti,^2,3and Marina E. Emborg^1,2,4

Academic Editor: Sushil K. Jha

Received05 Nov 2019

Revised13 Mar 2020

Accepted03 Apr 2020

Published02 May 2020

Abstract

Degeneration of sympathetic innervation of the heart occurs in numerous diseases, including diabetes, idiopathic REM sleep disorder, and Parkinson’s disease (PD). In PD, cardiac sympathetic denervation occurs in 80-90% of patients and can begin before the onset of motor symptoms. Today, there are no disease-modifying therapies for cardiac sympathetic neurodegeneration, and biomarkers are limited to radioimaging techniques. Analysis of expression levels of coding mRNA and noncoding RNAs, such as microRNAs (miRNAs), can uncover pathways involved in disease, leading to the discovery of biomarkers, pathological mechanisms, and potential drug targets. Whole blood in particular is a clinically relevant source of biomarkers, as blood sampling is inexpensive and simple to perform. Our research group has previously developed a nonhuman primate model of cardiac sympathetic denervation by intravenous administration of the catecholaminergic neurotoxin 6-hydroxydopamine (6-OHDA). In this rhesus macaque (Macaca mulatta) model, imaging with positron emission tomography showed that oral administration of the peroxisome proliferator-activated receptor gamma (PPARγ) agonist pioglitazone (; 5 mg/kg daily) significantly decreased cardiac inflammation and oxidative stress compared to placebo (). Here, we report our analysis of miRNA and mRNA expression levels over time in the whole blood of these monkeys. Differential expression of three miRNAs was induced by 6-OHDA (mml-miR-16-2-3p, mml-miR-133d-3p, and mml-miR-1262-5p) and two miRNAs by pioglitazone (mml-miR-204-5p and mml-miR-146b-5p) at 12 weeks posttoxin, while expression of mRNAs involved in inflammatory cytokines and receptors was not significantly affected. Overall, this study contributes to the characterization of rhesus coding and noncoding RNA profiles in normal and disease-like conditions, which may facilitate the identification and clinical translation of biomarkers of cardiac neurodegeneration and neuroprotection.

1. Introduction

Cardiac function is modulated by the sympathetic and parasympathetic input of the autonomic nervous system. Sympathetic innervation to the heart is critical to regulate cardiac activity during exercise, changes in posture, and other nonpathological activities that alter blood pressure or heart rate [1–3]. Patients with diabetes [4, 5], REM sleep disorder [6], and Parkinson’s disease (PD) frequently show loss of cardiac sympathetic innervation [7]. Sympathetic cardiac neurodegeneration occurs in 80-90% of PD patients, such that clinical guidelines now include in vivo imaging evidence of cardiac sympathetic denervation as a supportive diagnostic criterion [8–11].

Nonhuman primates are critical for modeling and studying diseases, due to their genetic and physiological similarity to humans [12]. Our research group has developed a rhesus macaque (Macaca mulatta) model of cardiac sympathetic denervation using intravenous administration of the neurotoxin 6-hydroxydopamine (6-OHDA) [13]. After entering neurons via catecholamine transporters, 6-OHDA autoxidizes and interferes with mitochondrial complex I, leading to increased oxidative stress and inflammation [14], which are neurodegenerative mechanisms shared with PD [15, 16]. We have shown that 6-OHDA-induced cardiac sympathetic loss, reactive oxygen species accumulation, and inflammatory cell recruitment can be detected by in vivo positron emission tomography (PET) [17]. Furthermore, administration of pioglitazone, a peroxisome proliferator-activated receptor gamma (PPARγ) agonist, significantly reduced cardiac oxidative stress and inflammation [17]. PPARγ is a type II nuclear receptor that classically works as a transcription factor by binding to PPARγ response elements (PPREs) in the promoters of target genes [18], including coding and noncoding RNAs [19, 20]. PPARγ activation inhibits the actions of nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), a transcription factor which increases immune cell expression of proinflammatory gene products including IL1β, TNFα, IL8, and IL6 [21, 22]. The cardiac PET findings and the transcriptional regulatory activity of PPARγ suggest that the effects of systemic 6-OHDA and pioglitazone on the heart may be, at least in part, due to alterations of the transcriptome of inflammatory cells.

Transcriptome analysis has evolved from focusing on protein coding mRNA to include noncoding RNAs, as the latter are recognized as key regulators of protein expression and cell function. Noncoding RNAs include microRNAs (miRNAs), which are short, approximately 22-nucleotide RNAs. By binding to the 3 untranslated region of target mRNAs, miRNAs inhibit translation of mRNAs into proteins. Each miRNA can target multiple mRNA sequences, affecting several functional pathways simultaneously [23]. miRNAs are key regulators of cellular health and disease progression [24] and can provide insight into disease mechanisms and serve as biomarkers. For example, analysis of whole blood samples has shown miRNA expression changes associated with renal sympathetic denervation to treat arterial hypertension [25], and circulating RNAs are an area of active research as potential biomarkers for PD [26–29]. However, alterations in mRNA and miRNA levels have not been studied in the context of cardiac sympathetic nerve loss in PD or animal models of PD. Additionally, our group is the first to evaluate miRNA expression in rhesus whole blood [30]. Here, we report our analysis of miRNA and mRNA expression levels over time in the whole blood of rhesus macaques with ongoing cardiac sympathetic neurodegeneration and PPARγ activation.

2. Materials and Methods

2.1. Ethics Statement

The present study was performed in strict accordance with the recommendations in the National Research Council Guide for the Care and Use of Laboratory Animals (2011) in an AAALAC-accredited facility (Wisconsin National Primate Research Center, University of Wisconsin-Madison). Experimental procedures were approved by the Institutional Animal Care and Use Committee (IACUC) of the University of Wisconsin-Madison (permit number: G00705). All efforts were made to minimize the number of animals used and to ameliorate any distress.

2.2. Subjects and Whole Blood Collection

Blood samples for this project were obtained from ten adult, male rhesus macaques treated with systemic 6-OHDA (50 mg/kg IV, [13]) that were part of a previously published study [17]. Twenty-four hours post-6-OHDA, animals were randomly assigned to receive daily oral dosing of placebo (; 6.2–13.0 years old; 9.8–12.3 kg) or PPARγ agonist pioglitazone (5 mg/kg; ; 5.6–11.4 years old; 9.4–10.6 kg). Blood samples were collected in PAXgene blood RNA tubes (Qiagen, Venlo, Netherlands) under ketamine anesthesia (15 mg/kg IM) before, 1 week after, and 12 weeks after 6-OHDA and stored at -80°C until RNA extraction. All 30 samples were randomized and arbitrarily numbered to achieve blinding.

2.3. RNA Isolation

Total RNA was extracted following 24 hr incubation at room temperature. Immediately after isolation, RNA quantity and quality were assessed using NanoDrop 2000 (Thermo Fisher Scientific, Halethorpe, MD, USA) and Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA) (Supplementary Table 1). Aliquots of RNA were stored at -80°C.

2.4. miRNA Library Preparation and Next-Generation Sequencing (NGS)

Library preparation, NGS, and data analysis were performed by Qiagen Genomic Services (Frederick, Maryland). Additionally, Qiagen Genomic Services assessed nucleic acid quality before and after library preparation using Bioanalyzer 2100 and TapeStation 4200 (Agilent, Santa Clara, CA, USA). The miRNA sequencing libraries were prepared with the QIAseq miRNA Library Kit (Qiagen, Venlo, Netherlands) using 100 ng total RNA. Adapters containing unique molecular indices (UMIs) were ligated to the RNA. miRNAs and small RNAs (e.g., small nucleolar RNAs (snoRNAs) and small nuclear RNAs (snRNAs)) were converted to cDNA amplified using PCR (16 cycles). Amplified RNAs were sequenced using the NextSeq 500 (Illumina, San Diego, CA, USA) generating single-end reads with a sequencing read length of 76 nucleotides. Raw data was converted to FASTQ files using the bcl2fastq software and checked using the FastQC tool. Reads with identical UMIs were collapsed into a single read and aligned to miRBase using Bowtie2 to identify known miRNAs in rhesus macaques to generate miRNA counts. The reference genome used was Mmul_1 and annotation reference with miRBase_20. Small RNAs are reads that align to the Qiagen small RNA database. Unmapped reads that did not align with rhesus or other species in miRBase were analyzed by the miRPara tool to identify putative miRNAs [31]. RNA read counts were normalized for sample sequencing depth by calculating tags per million (TPM), the number of reads for a particular miRNA in an individual subject divided by the total number of mapped reads in that subject and multiplied by 1 million.

2.5. mRNA RT² Profiler PCR Array and Additional Gene Reverse Transcription-Quantitative PCR (RT-qPCR)

RT-qPCR was used to assess genes included on Inflammatory Cytokines and Receptors RT² Profiler PCR Array (Qiagen, Venlo, Netherlands) (Supplementary Table 2) and additional genes of interest (NFKBIA, CD36, STAT1, and MAFB) (Supplementary Table 3). 5 μg of total RNA was synthesized into cDNA using the RT² First Strand Kit (Qiagen, Venlo, Netherlands) according to the manufacturer’s instructions and then either used immediately or stored at -80°C. All qPCR was carried out on the Applied Biosystems QuantStudio 3 (Thermo Fisher, Waltham, MA, USA) using RT² SYBR Green ROX qPCR Mastermix (Qiagen, Venlo, Netherlands). All control wells in the PCR arrays were within expected cycles to threshold (Cts). Additional gene qPCR was performed on 96-well plates with each gene of interest on the same plate as the reference gene, each sample in duplicate, and the inclusion of no template control wells for each gene and the reference gene. Duplicate wells for each sample were averaged for data analysis. The Ct of all no template control wells was >35. NFKBIA, CD36, STAT1, and MAFB primers were designed using the PrimerQuest online tool (Integrated DNA Technologies); qPCR products were loaded on 2% agarose gels, and amplification of a single band of the predicted size was verified for each primer pair (Supplementary Figure 1). Fold differences in gene expression were calculated using the 2^−ΔΔCt fold change method and presented as fold change [32]. Cts were normalized to the reference gene RPL13A; primers for RPL13A were included on the PCR arrays, and the same primers were purchased for additional gene RT-qPCR (Product no. 33001; Cat. no. PPQ00210B-200). This gene was selected based on having the most stable expression across all samples and previous results illustrating RPL13A as a suitable reference gene specifically in rhesus macaques [33].

2.6. Data Analysis

miRNA or other small RNA differential expression analysis was carried out using edgeR (Bioconductor). For normalization, the trimmed mean of value method based on log-fold and absolute gene-wise changes in expression levels between samples (TMM normalization) was used. All values were corrected for multiple corrections using the false discovery rate (FDR) correction method of Benjamini and Hochberg. RNAs were considered significantly differentially expressed when the FDR-corrected value was <0.05 and fold change was >|2|. Principal components analysis (PCA) was performed with R using TMM-normalized quantifications. PCA plots identified one placebo group animal (Placebo 3) as a possible outlier (Supplementary Figure 2). Differential expression analysis was performed with and without the outlier (Supplementary Table 4; Supplementary Table 5; Supplementary Table 6); results from the analysis without the outlier are reported and discussed. miRBase version 22.1 was used for comparison between the rhesus macaque (mml-), human (hsa-), mouse (mmu-), and rat (rno-) mature miRNA (miR) and miRNA precursor (mir) sequences.

mRNA gene targets of miRNAs were identified using TargetScan (v7.2) with species selected as “Rhesus” (http://www.targetscan.org). Enrichr (http://amp.pharm.mssm.edu/Enrichr/) [34, 35] was used to perform Gene Ontology (GO) function annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) biological pathway analyses. The Enrichr ontology functions analyzed were GO Biological Process 2018, GO Molecular Function 2018, and GO Cellular Component 2018. The Enrichr pathway function analyzed was KEGG 2016. GO and KEGG term values were calculated using the Fisher exact and corrected for multiple comparisons using FDR. Functional enrichment analysis was limited to mml-miR-16-2-3p as it was the most abundant differentially expressed miRNA (showed high expression levels (>10 TMM) in all groups) and exhibited the highest fold change between groups.

mRNA data was analyzed for differential expression using two-way repeated-measures analysis of variance (ANOVA) in SPSS (version 24). were replaced with 35 prior to calculating ΔCt for data analysis. All values were corrected for multiple comparisons using FDR (http://www.sdmproject.com/utilities/?show=FDR); uncorrected and FDR-corrected values are reported. Volcano plots were constructed in R (version 3.5.1).

3. Results and Discussion

3.1. NGS Overview

NGS detected an average of 10.3 million UMI-corrected reads in each library (Supplementary Table 7) at a high base quality (Supplementary Figure 3). Averaging the data from all 30 samples, 85.86% of UMI-corrected reads were mapped to known rhesus macaque miRNAs, 0.29% were mapped to small RNAs other than miRNAs, and 10.04% were aligned to the rhesus reference genome, but not to miRNA or small RNAs (Supplementary Table 7). Additionally, 0.14% of the reads represented miRNAs not before documented in rhesus macaques or putative miRNAs that did not map to any species in miRBase; a description of these novel miRNA findings has been recently published by our group [30].

3.2. miRNA

3.2.1. miRNA Expression in Rhesus Macaque Whole Blood

Of the 912 known rhesus mature miRNAs, 419 were observed in the whole blood in at least one sample of this study above a threshold of 1 TPM (Supplementary Table 8). A total of 266 and 158 miRNAs were detected at 1 TPM or 10 TPM, respectively, in all 30 of the samples (Figure 1). miRNA expression is known to show large variability between tissues [36, 37], likely explaining why this study detected approximately 50% of the known rhesus miRNome.

In naïve, baseline samples representative of the whole blood miRNA profile in normal rhesus, mml-miR-486-5p, mml-miR-16-5p, mml-miR-92a-3p, mml-miR-191-5p, and mml-miR-25 were the top 5 most abundant mature miRNAs and represented over 50% of total miRNA reads (Figure 2(a)). These findings are similar to other species, including humans [38] and rats [39], and support the concept that these miRNAs are conserved between species (Table 1). In humans, hsa-miR-16-5p, hsa-miR-486-5p, and hsa-miR-92a constitute 20.6%, 12.6%, and 12.6% of the mature miRNAs in red blood cells, respectively [40]. Our results in rhesus suggest that, like humans [41, 42], the whole blood miRNA expression pattern reflects red blood cells’ transcriptome more than leukocytes.

(a)

(b)

3.2.2. miRNA Differential Expression

Unsupervised PCA analysis indicated that each individual animal exhibited a unique miRNA profile that was stable over time, as demonstrated by samples clustering closely based on the individual (Figure 3(b)) but not treatment (Figure 3(a)). Relative to inbred rodent strains, outbred research species such as rhesus macaques exhibit high individual genetic variability [43]. This genetic variability produces individual differences in RNA expression levels [44], likely contributing to the sample clustering pattern observed in PCA (Figure 3). Despite this, five miRNAs were detected which had statistically significant changes in expression levels following systemic 6-OHDA or PPARγ activation (Table 2), including mml-miR-16-2-3p, mml-miR-204-5p, mml-miR-133d-3p, mml-miR-1262-5p, and mml-miR-146b-5p.

(a)

(b)

Systemic 6-OHDA increased the expression of mml-miR-16-2-3p and decreased mml-miR-133d-3p and mml-miR-1262-5p 12 weeks posttoxin, identifying these miRNAs as potential biomarkers of the long term (12 week) effects of sympathetic denervation. No miRNAs were found to be differentially expressed at the 1 week timepoint. The impact of increased mml-miR-16-2-3p levels 12 weeks posttoxin is unclear. Although mml-miR-16-2-3p was predicted to have over 3000 mRNA targets (Supplementary Table 9), this list of targets was not significantly enriched for mRNAs involved in any specific biological pathway defined in the KEGG database (Supplementary Table 10). The 16 GO terms significantly enriched for miR-16-2-3p gene targets were categorized into 12 biological processes and 4 molecular functions, including ubiquitin ligase activity and regulation gene transcription and expression (Table 3; Supplementary Table 11). These functions suggest that mml-miR-16-2-3p affects a broad range of cellular pathways and could have wide-reaching effects on multiple cell types. The few studies available on miR-16-2-3p have linked it to cancers [45, 46], craniofacial malformations [47, 48], treatment of PD with Levodopa/Carbidopa [49], and multiple sclerosis [50], supporting diverse roles for this miRNA. Future studies investigating the blood cell type(s), or noncellular sources such as exosomes [51, 52], in which mml-miR-16-2-3p expression is altered by 6-OHDA will give insight into the relationship between sympathetic denervation and mml-miR-16-2-3p.

Information about mml-miR-133d-3p and mml-miR-1262-5p is also limited. Mml-miR-133d-3p is now recognized to be mml-miR-133a in the current miRBase version 22.1 following the release of the rhesus macaque Mmul_8.0.1 genome assembly. Mml-miR-133a contains a different mRNA binding seed sequence, the nucleotides in the 2-7 positions at the 5 end of the miRNA, than the most similar human sequence, hsa-miR-133a-3p (Table 4). The rhesus mml-miR-1262-5p sequence is dissimilar from the human hsa-miR-1262, and a search using miRBase does not indicate an equivalent miRNA sequence in humans. Publications are not currently available on the function of mml-miR-133d-3p and mml-miR-1262-5p in rhesus, and the dissimilarity with human miRNA sequences cautions against comparison with reports in humans. The function of these RNAs will hopefully be elucidated as new studies on miRNAs in rhesus macaques emerge.

PPARγ activation by pioglitazone led to increased expression of mml-miR-146b-5p from baseline to 12 weeks post-6-OHDA, in addition to decreased levels of mml-miR-16-2-3p and mml-miR-204-5p in pioglitazone-treated animals relative to placebo-treated at the twelve-week time point. In this context, these miRNAs emerge as potential biomarkers of the long term (12 week) effects of PPARγ activation. PPARγ activation in immune cells is well documented to downregulate activity of the transcription factor NF-κB and downstream pro-inflammatory cytokines and chemokines [21, 22]. Fitting this function, miR-204-5p and miR-146b-5p are involved in inflammation. In mice, mmu-miR-204-5p is upregulated in pro-inflammatory type macrophages relative to anti-inflammatory macrophages [53]. Downregulation of mml-miR-204-5p in the current study may reflect a PPARγ-associated skewing of circulating monocytes to an anti-inflammatory phenotype. Interestingly, the opposite effect was found in human pulmonary arterial smooth muscle cells, in which hsa-miR-204-5p expression decreased following repression of PPARγ [54], suggesting a differential effect of PPARγ by cell type. Future work investigating the function of miR-204-5p in blood cells and pulmonary arterial smooth muscle with be important to understanding the impact of altered miR-204-5p expression in these cell types.

In humans, hsa-miR-146-5p has been identified as a critical regulator of the immune response, as well as in the pathophysiology of several cancers, renal, and cardiac diseases [55, 56]. Hsa-miR-146-5p exerts an anti-inflammatory effect by decreasing expression of proteins such as TNF receptor-associated factor 6 (TRAF6), a signal transducer upstream of NF-κB activation [57–59]. Our findings of differential expression of mml-miR-146b-5p and mml-miR-204-5p suggest a relationship between these miRNAs and the anti-inflammatory actions of pioglitazone observed in the heart of these animals [17]. In contrast to our findings, pioglitazone administration led to decreased levels of rno-miR-146b-5p in the rat heart in a model of heart failure; this result followed an increase rno-miR-146b-5p associated with the induction of heart failure in the model [60]. Additional research will be needed to shed light on whether these differential effects of pioglitazone on miR-146-5p expression levels are due to differences between species, tissues examined (cardiac vs. whole blood), or another factor. To the authors’ knowledge, the present study is the first to report alterations of miRNAs in whole blood following pioglitazone administration. One publication describing miRNA alterations in serum found that hsa-miR-24 levels were increased in diabetic patients following 9 months of pioglitazone [61]; miR-24 expression was not impacted in the present study. Additional work is needed to assess replicability of differentially expressed miRNAs across studies in order to identify them as reliable biomarkers.

Although in this study 6-OHDA was used to model cardiac sympathetic neurodegeneration, it is important to recognize that the neurotoxin produces sympathetic denervation outside of the heart, including spleen and bone marrow. Peripheral sympathectomy with systemic 6-OHDA is reported to alter immune cell development and phenotype, leading to complex effects on the immune response [62–65]. Whole blood miRNA changes in our monkeys may be related to the direct effects of sympathetic denervation of lymphoid organs, rather than inflammation induced by cardiac neurodegeneration. As qPCR validation of differentially expressed miRNAs was not performed in this study due to limited sample availability, candidate biomarkers should be further validated in future studies.

3.3. Other Small RNAs

3.3.1. Other Small RNA Expression in Rhesus Macaque Whole Blood

Small RNAs other than miRNAs identified by NGS included snoRNAs and snRNAs (Supplementary Table 5; Supplementary Table 6; Supplementary Table 8). A total of 85 small RNAs were present in all 30 samples at at least 1 TPM and 16 were present in all samples at at least 10 TPM (Supplementary Table 8). Of the 251 small RNAs detected at at least 1 TPM in any sample, 203 were snoRNAs and 48 were snRNAs (Supplementary Table 8). The QIAseq miRNA library kit used for the NGS portion of this study was developed primarily for miRNAs [66], thus the results may not reflect the full profile of snoRNAs and snRNAs in the whole blood of rhesus macaques.

In naïve, baseline samples, the top 5 most abundant snoRNAs or snRNAs detected in this study based on average TPM were the snoRNAs SNORND104 (ENSMMUG00000025102), SNORD83 (ENSMMUG00000026474), SCARNA3 (ENSMMUG00000036503), SNORD69 (ENSMMUG00000033050), and SNORD100 (ENSMMUG00000033610), which accounted for 40% of reads (Figure 2(b); Supplementary Table 8). The cellular function of snoRNAs is an area of active research [67] with limited information currently available [68–70]. Canonically, snoRNAs contain sequences complementary to ribosomal RNAs (rRNAs), which guide ribonucleoproteins necessary for rRNA modification such as 2-O’-methylation and pseudouridylation [67, 71].

3.3.2. Other Small RNA Differential Expression

At 12 weeks post-6-OHDA, the snoRNA SNORA46 was significantly downregulated in the 6-OHDA + Pioglitazone group relative to 6-OHDA + Placebo (Table 2). 6-OHDA + Pioglitazone animals exhibited a significant increase in expression of SNORD15 from baseline to 12 weeks posttoxin (Table 2). No snoRNAs or snRNAs were found to be differentially expressed 1 week post-6-OHDA. SNORA46 is predicted to guide the modification of uridine (U469) in ribosomal 18S rRNA [72]. SNORA15, also known as U15, is predicted to guide the modification of adenosine (A3764) in 28S rRNA [70]. Currently, only one report exists of SNORA46 or SNORD15 involvement in human disease, showing auto-antibodies to SNORD15 in patients with systemic sclerosis [73]. The functional impact of altered snoRNA expression following systemic 6-OHDA or PPARγ activation is unclear due to limited research on the effect of the above rRNA modifications on cellular activity.

3.4. mRNA

Analysis of mRNA expression levels in whole blood suggested some differences between groups (Figure 4). However, correction of values for multiple comparisons using the FDR correction revealed that no genes assessed in this study were significantly differentially expressed at the 1 week or 12 week time points (Supplementary Table 12; Supplementary Table 13). This lack of change in expression levels of inflammatory cytokine and receptor genes in whole blood is in contrast to the significant 6-OHDA-induced increase in cardiac inflammation and oxidative stress and the pioglitazone-associated amelioration of these neurodegenerative mechanisms observed 1 week post-6-OHDA by PET [17]. Pioglitazone also led to increased CD36 protein expression in the heart in these animals as detected by immunohistochemistry [74], but CD36 transcript levels in the blood were not found to be altered. These findings likely reflect differences in how the treatments affected cells in the heart compared to cells in the blood and also the unique transcriptional profiles of these tissues. Additionally, three differentially expressed miRNAs in this study, mml-miR-16-2-3p, mml-miR-133d-3p, and mml-miR-1262-5p, were predicted to target over 30 mRNAs assayed by PCR array and additional RT-qPCR (Supplementary Table 14), yet our results did not identify altered expression of these target mRNAs.

(a)

(b)

(c)

(d)

Figure 4

Volcano plots of differentially expressed mRNAs using uncorrected values. Volcano plots illustrating differential expression between the 6-OHDA + Placebo and 6-OHDA + Pioglitazone groups at (a) baseline, (b) 1 week post-6-OHDA, (c) 12 weeks post-6-OHDA and (d) between the baseline and 12 weeks post-6-OHDA timepoints in the 6-OHDA + Pioglitazone group. Each circle represents one gene. Green filled circles represent genes for which the uncorrected value is <0.05 and Log2 fold change is >1. Note that positive Log2 fold change indicates higher expression in the group to the right in the header (e.g., CCR8 showed higher expression in the 6-OHDA + Pioglitazone group at baseline) and negative Log2 fold change indicates higher expression in the group to the left. CCR8, C-C motif chemokine receptor 8; VEGFA, vascular endothelial growth factor A; CXCL9, C-X-C motif chemokine ligand 9; CD70, CD70 molecule; CCR6, C-C motif chemokine receptor 6; CCR4, C-C motif chemokine receptor 4.

The lack of significant treatment effect on mRNA expression levels in our study may be related to the use of whole blood as the source of the RNA samples. Whole blood sampling was selected to maximize clinical relevancy. Whole blood samples are easy and inexpensive to collect, and the PAXgene tubes used in this study are FDA-cleared and CE-marked for in vitro diagnostic use. The use of whole blood RNA, however, can present technical challenges. Globin mRNA can represent 50-80% of mRNA transcripts within the total RNA extracted from whole blood, lowering sensitivity [75, 76]. Although methods are available to reduce globin mRNA content, these additional steps may increase sample variation, alter RNA profiles, and reduce RNA yield [77–79]. Previous reports have shown that pioglitazone significantly altered mRNA expression in human peripheral blood mononuclear cells (PBMCs) [80, 81]. Future work in this rhesus 6-OHDA model investigating mRNA changes in globin reduced whole blood RNA or in RNA extracted from PBMCs might have increased sensitivity to detect alterations in the inflammatory cytokine and receptor genes assessed in the present study.

4. Conclusions

At twelve weeks postneurotoxin, differential expression of three miRNAs was induced by 6-OHDA (mml-miR-16-2-3p, mml-miR-133d-3p, and mml-miR-1262-5p) and two miRNAs by pioglitazone (mml-miR-204-5p and mml-miR-146b-5p), providing initial evidence that these RNAs are circulating markers of the long term (12 week) effects of sympathetic denervation and PPARγ activation, respectively. Expression of inflammatory cytokine and receptor mRNAs in rhesus macaque whole blood was unaffected by systemic 6-OHDA and PPARγ activation.

Overall, this study contributes to the characterization of rhesus coding and noncoding RNA profiles in normal and disease-like conditions, which may facilitate the identification and clinical translation of biomarkers of cardiac neurodegeneration and neuroprotection.

Data Availability

The data used to support the findings of this study are available in the supplementary materials or from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge Dr. Mark Berres, Dane Shank, Carissa Boettcher, Dr. Kevin Brunner, and the dedicated animal care and veterinary staff at the Wisconsin National Primate Research Center for their technical support. Research was supported by the National Institutes of Health (P51OD011106 and F31HL136047) and the University of Wisconsin-Madison Office of Vice Chancellor for Research and Graduate Education, Cellular and Molecular Pathology Graduate Program, and Department of Medical Physics.

Supplementary Materials

Supplementary 1. Supplementary Figure 1: gel electrophoresis of PCR products. Supplementary Figure 2: principal component analysis plots showing unsupervised assessment of placebo group samples. Supplementary Figure 3: average base quality of the UMI-corrected reads in each of the 30 RNA samples.

Supplementary 2. Supplementary Table 1: quality and concentration of all 30 extracted total RNA samples.

Supplementary 3. Supplementary Table 2: list of genes on the RT² 475 rhesus inflammatory cytokine and receptor PCR array.

Supplementary 4. Supplementary Table 3: sequences of primers used for RT-qPCR.

Supplementary 5. Supplementary Table 4: differentially expressed small RNAs.

Supplementary 6. Supplementary Table 5: results of small RNA differential expression analysis between groups with all animals included ().

Supplementary 7. Supplementary Table 6: results of small RNA differential expression analysis between groups without the outlier animal ().

Supplementary 8. Supplementary Table 7: summary of sequencing reads and mapping for each sample.

Supplementary 9. Supplementary Table 8: abundance of miRNAs and other small RNAs.

Supplementary 10. Supplementary Table 9: mRNA targets of mml-miR-16-2-3p predicted by TargetScan.

Supplementary 11. Supplementary Table 10: results of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis.

Supplementary 12. Supplementary Table 11: results of gene ontology (GO) enrichment analysis.

Supplementary 13. Supplementary Table 12: list of uncorrected values, FDR-corrected values, and fold change for all genes on the rhesus inflammatory cytokine and receptor RT² profiler array.

Supplementary 14. Supplementary Table 13: list of uncorrected values, FDR-corrected values, and fold change for additional gene qPCR.

Supplementary 15. Supplementary Table 14: predicted mRNA targets of mml-miR-16-2-3p, mml-miR-204-5p, mml-miR-133d-3p, mml-miR-1262-5p, and mml-miR-146b-5p.

References

A. Robles-Cabrera, A. Michel-Chavez, R. C. Callejas-Rojas, C. Malamud-Kessler, G. Delgado et al., “The cardiovagal, cardiosympathetic and vasosympathetic arterial baroreflexes and the neural control of short-term blood pressure,” Revista de Neurologia, vol. 59, pp. 508–516, 2014.
View at: Google Scholar
J. M. Stewart, “Mechanisms of sympathetic regulation in orthostatic intolerance,” Journal of Applied physiology, vol. 113, no. 10, pp. 1659–1668, 2012.
View at: Publisher Site | Google Scholar
J. L. Elghozi and C. Julien, “Sympathetic control of short-term heart rate variability and its pharmacological modulation,” Fundamental & Clinical Pharmacology, vol. 21, no. 4, pp. 337–347, 2007.
View at: Publisher Site | Google Scholar
M. C. Gerson, J. H. Caldwell, K. Ananthasubramaniam et al., “Influence of diabetes mellitus on prognostic utility of imaging of myocardial sympathetic innervation in heart failure patients,” Circulation Cardiovascular Imaging, vol. 4, no. 2, pp. 87–93, 2011.
View at: Publisher Site | Google Scholar
M. J. Stevens, D. M. Raffel, K. C. Allman, M. Schwaiger, and D. M. Wieland, “Regression and progression of cardiac sympathetic dysinnervation complicating diabetes: an assessment by c-11 hydroxyephedrine and positron emission tomography,” Metabolism, vol. 48, no. 1, pp. 92–101, 1999.
View at: Publisher Site | Google Scholar
T. Miyamoto, M. Miyamoto, Y. Inoue, Y. Usui, K. Suzuki, and K. Hirata, “Reduced cardiac 123i-mibg scintigraphy in idiopathic rem sleep behavior disorder,” Neurology, vol. 67, no. 12, pp. 2236–2238, 2006.
View at: Publisher Site | Google Scholar
F. Courbon, C. Brefel-Courbon, C. Thalamas et al., “Cardiac mibg scintigraphy is a sensitive tool for detecting cardiac sympathetic denervation in parkinson's disease,” Movement Disorders, vol. 18, no. 8, pp. 890–897, 2003.
View at: Publisher Site | Google Scholar
R. B. Postuma, D. Berg, M. Stern et al., “Mds clinical diagnostic criteria for parkinson's disease,” Movement Disorders, vol. 30, no. 12, pp. 1591–1601, 2015.
View at: Publisher Site | Google Scholar
D. S. Goldstein, “Cardiac denervation in patients with parkinson disease,” Cleveland Clinic Journal of Medicine, vol. 74, Supplement 1, pp. S91–S94, 2007.
View at: Publisher Site | Google Scholar
S. Jain and D. S. Goldstein, “Cardiovascular dysautonomia in parkinson disease: from pathophysiology to pathogenesis,” Neurobiology of Disease, vol. 46, no. 3, pp. 572–580, 2012.
View at: Publisher Site | Google Scholar
S. Orimo, E. Ozawa, S. Nakade, T. Sugimoto, and H. Mizusawa, “¹²³I-metaiodobenzylguanidine myocardial scintigraphy in Parkinson's disease,” Journal of Neurology, Neurosurgery, and Psychiatry, vol. 67, no. 2, pp. 189–194, 1999.
View at: Publisher Site | Google Scholar
J. P. Capitanio and M. E. Emborg, “Contributions of non-human primates to neuroscience research,” The Lancet, vol. 371, no. 9618, pp. 1126–1135, 2008.
View at: Publisher Site | Google Scholar
V. Joers, K. Seneczko, N. C. Goecks et al., “Nonuniform cardiac denervation observed by 11c-meta-hydroxyephedrine pet in 6-ohda-treated monkeys,” PLoS One, vol. 7, no. 4, article e35371, 2012.
View at: Publisher Site | Google Scholar
J. Rodriguez-Pallares, J. A. Parga, A. Muñoz, P. Rey, M. J. Guerra, and J. L. Labandeira-Garcia, “Mechanism of 6-hydroxydopamine neurotoxicity: the role of nadph oxidase and microglial activation in 6-hydroxydopamine-induced degeneration of dopaminergic neurons,” Journal of Neurochemistry, vol. 103, no. 1, pp. 145–156, 2007.
View at: Publisher Site | Google Scholar
E. C. Hirsch, S. Vyas, and S. Hunot, “Neuroinflammation in parkinson's disease,” Parkinsonism & Related Disorders, vol. 18, Supplement 1, pp. S210–S212, 2012.
View at: Publisher Site | Google Scholar
J. Blesa, I. Trigo-Damas, A. Quiroga-Varela, and V. R. Jackson-Lewis, “Oxidative stress and parkinson's disease,” Frontiers in Neuroanatomy, vol. 9, p. 91, 2015.
View at: Publisher Site | Google Scholar
J. M. Metzger, C. F. Moore, C. A. Boettcher et al., “In vivo imaging of inflammation and oxidative stress in a nonhuman primate model of cardiac sympathetic neurodegeneration,” NPJ Parkinson's Disease, vol. 4, no. 1, article 22, 2018.
View at: Publisher Site | Google Scholar
R. Kapadia, J. H. Yi, and R. Vemuganti, “Mechanisms of anti-inflammatory and neuroprotective actions of ppar-gamma agonists,” Frontiers in Bioscience, vol. 13, pp. 1813–1826, 2008.
View at: Publisher Site | Google Scholar
D. Wang, L. Shi, W. Xin et al., “Activation of PPARγ inhibits pro-inflammatory cytokines production by upregulation of miR-124 in vitro and in vivo,” Biochemical and Biophysical Research Communications, vol. 486, no. 3, pp. 726–731, 2017.
View at: Publisher Site | Google Scholar
I. S. Bae, P. J. Park, J. H. Lee, E. G. Cho, T. R. Lee, and S. H. Kim, “PPARγ-mediated G-protein coupled receptor 120 signaling pathway promotes transcriptional activation of miR-143 in adipocytes,” Gene, vol. 626, pp. 64–69, 2017.
View at: Publisher Site | Google Scholar
P. M. Flood, L. Qian, L. J. Peterson et al., “Transcriptional Factor NF-κB as a Target for Therapy in Parkinson's Disease,” Parkinson’s Disease, vol. 2011, Article ID 216298, 8 pages, 2011.
View at: Publisher Site | Google Scholar
V. H. Perry and J. Teeling, “Microglia and macrophages of the central nervous system: the contribution of microglia priming and systemic inflammation to chronic neurodegeneration,” Seminars in Immunopathology, vol. 35, no. 5, pp. 601–612, 2013.
View at: Publisher Site | Google Scholar
N. Bushati and S. M. Cohen, “Microrna functions,” Annual Review of Cell and Developmental Biology, vol. 23, pp. 175–205, 2007.
View at: Publisher Site | Google Scholar
S. Kumar, M. Vijayan, J. S. Bhatti, and P. H. Reddy, “Micrornas as peripheral biomarkers in aging and age-related diseases,” Progress in Molecular Biology and Translational Science, vol. 146, pp. 47–94, 2017.
View at: Publisher Site | Google Scholar
O. Dörr, C. Liebetrau, H. Möllmann et al., “Effect of renal sympathetic denervation on specific micrornas as an indicator of reverse remodeling processes in hypertensive heart disease,” Journal of Clinical Hypertension, vol. 18, no. 6, pp. 497–502, 2016.
View at: Publisher Site | Google Scholar
S. K. Khoo, D. Petillo, U. J. Kang et al., “Plasma-based circulating microrna biomarkers for parkinson's disease,” Journal of Parkinson's Disease, vol. 2, no. 4, pp. 321–331, 2012.
View at: Publisher Site | Google Scholar
A. E. Roser, L. Caldi Gomes, J. Schunemann, F. Maass, and P. Lingor, “Circulating mirnas as diagnostic biomarkers for parkinson's disease,” Frontiers in Neuroscience, vol. 12, p. 625, 2018.
View at: Publisher Site | Google Scholar
M. Behbahanipour, M. Peymani, M. Salari, M. S. Hashemi, M. H. Nasr-Esfahani, and K. Ghaedi, “Expression profiling of blood micrornas 885, 361, and 17 in the patients with the parkinson's disease: integrating interaction data to uncover the possible triggering age-related mechanisms,” Scientific Reports, vol. 9, no. 1, article 13759, 2019.
View at: Publisher Site | Google Scholar
M. Martins, A. Rosa, L. C. Guedes et al., “Convergence of mirna expression profiling, α-Synuclein interacton and gwas in parkinson's disease,” PLoS One, vol. 6, no. 10, article e25443, 2011.
View at: Publisher Site | Google Scholar
M. S. Lopez, J. M. Metzger, and M. E. Emborg, “Identification of novel rhesus macaque micrornas from naïve whole blood,” Molecular Biology Reports, vol. 46, no. 5, pp. 5511–5516, 2019.
View at: Publisher Site | Google Scholar
Y. Wu, B. Wei, H. Liu, T. Li, and S. Rayner, “Mirpara: a svm-based software tool for prediction of most probable microrna coding regions in genome scale sequences,” BMC Bioinformatics, vol. 12, no. 1, p. 107, 2011.
View at: Publisher Site | Google Scholar
K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative pcr and the 2^-ΔΔC_T method,” Methods, vol. 25, pp. 402–408, 2001.
View at: Publisher Site | Google Scholar
K. Ahn, J. W. Huh, S. J. Park et al., “Selection of internal reference genes for sybr green qrt-pcr studies of rhesus monkey (Macaca Mulatta) tissues,” BMC Molecular Biology, vol. 9, no. 1, p. 78, 2008.
View at: Publisher Site | Google Scholar
E. Y. Chen, C. M. Tan, Y. Kou et al., “Enrichr: interactive and collaborative html5 gene list enrichment analysis tool,” BMC Bioinformatics, vol. 14, no. 1, p. 128, 2013.
View at: Publisher Site | Google Scholar
M. V. Kuleshov, M. R. Jones, A. D. Rouillard et al., “Enrichr: a comprehensive gene set enrichment analysis web server 2016 update,” Nucleic Acids Research, vol. 44, no. W1, pp. W90–W97, 2016.
View at: Publisher Site | Google Scholar
M. Dannemann, K. Prüfer, E. Lizano, B. Nickel, H. A. Burbano, and J. Kelso, “Transcription factors are targeted by differentially expressed mirnas in primates,” Genome Biology and Evolution, vol. 4, no. 4, pp. 552–564, 2012.
View at: Publisher Site | Google Scholar
J. Yue, Y. Sheng, and K. E. Orwig, “Identification of novel homologous microrna genes in the rhesus macaque genome,” BMC Genomics, vol. 9, no. 1, p. 8, 2008.
View at: Publisher Site | Google Scholar
S. Juzenas, G. Venkatesh, M. Hübenthal et al., “A comprehensive, cell specific microrna catalogue of human peripheral blood,” Nucleic Acids Research, vol. 45, no. 16, pp. 9290–9301, 2017.
View at: Publisher Site | Google Scholar
F. Gao, P. Liu, J. Narayanan et al., “Changes in mirna in the lung and whole blood after whole thorax irradiation in rats,” Scientific Reports, vol. 7, no. 1, pp. 44132–44132, 2017.
View at: Publisher Site | Google Scholar
L. Vu, V. Ragupathy, S. Kulkarni, and C. Atreya, “Analysis of argonaute 2-microrna complexes in ex vivo stored red blood cells,” Transfusion, vol. 57, no. 12, pp. 2995–3000, 2017.
View at: Publisher Site | Google Scholar
S. Y. Chen, Y. Wang, M. J. Telen, and J. T. Chi, “The genomic analysis of erythrocyte microrna expression in sickle cell diseases,” PLoS One, vol. 3, no. 6, article e2360, 2008.
View at: Publisher Site | Google Scholar
C. W. Wu, X. Cao, C. K. Berger et al., “Novel approach to fecal occult blood testing by assay of erythrocyte-specific microrna markers,” Digestive Diseases and Sciences, vol. 62, no. 8, pp. 1985–1994, 2017.
View at: Publisher Site | Google Scholar
E. J. Vallender and G. M. Miller, “Nonhuman primate models in the genomic era: a paradigm shift,” ILAR Journal, vol. 54, no. 2, pp. 154–165, 2013.
View at: Publisher Site | Google Scholar
A. Battle, S. Mostafavi, X. Zhu et al., “Characterizing the genetic basis of transcriptome diversity through RNA-sequencing of 922 individuals,” Genome Research, vol. 24, no. 1, pp. 14–24, 2014.
View at: Publisher Site | Google Scholar
Q. Su, E. C. Zhu, Y.-L. Qu et al., “Serum level of co-expressed hub mirnas as diagnostic and prognostic biomarkers for pancreatic ductal adenocarcinoma,” Journal of Cancer, vol. 9, no. 21, pp. 3991–3999, 2018.
View at: Publisher Site | Google Scholar
X. Liu, H. Cai, W. Sheng, H. Huang, Z. Long, and Y. Wang, “Micrornas expression profile related with response to preoperative radiochemotherapy in patients with locally advanced gastric cancer,” BMC Cancer, vol. 18, no. 1, article 1048, 2018.
View at: Publisher Site | Google Scholar
J. Zou, J. Li, C. Ji, Q. Li, and X. Guo, “Expression profile of plasma micrornas in nonsyndromic cleft lip and their clinical significance as biomarkers,” Biomedicine & Pharmacotherapy, vol. 82, pp. 459–466, 2016.
View at: Publisher Site | Google Scholar
Y. Gu, Y. Zhang, C. Zhao et al., “Serum micrornas as potential biomarkers of mandibular prognathism,” Oral Diseases, vol. 20, no. 1, pp. 55–61, 2014.
View at: Publisher Site | Google Scholar
R. Margis, R. Margis, and C. R. Rieder, “Identification of blood microRNAs associated to Parkinsońs disease,” Journal of Biotechnology, vol. 152, no. 3, pp. 96–101, 2011.
View at: Publisher Site | Google Scholar
A. Keller, P. Leidinger, F. Steinmeyer et al., “Comprehensive analysis of microrna profiles in multiple sclerosis including next-generation sequencing,” Multiple Sclerosis, vol. 20, no. 3, pp. 295–303, 2014.
View at: Publisher Site | Google Scholar
D. Sanz-Rubio, I. Martin-Burriel, A. Gil et al., “Stability of circulating exosomal mirnas in healthy subjects,” Scientific Reports, vol. 8, no. 1, article 10306, 2018.
View at: Publisher Site | Google Scholar
W. Stoorvogel, “Functional transfer of microrna by exosomes,” Blood, vol. 119, no. 3, pp. 646–648, 2012.
View at: Publisher Site | Google Scholar
Y. Zhang, M. Zhang, M. Zhong, Q. Suo, and K. Lv, “Expression profiles of mirnas in polarized macrophages,” International Journal of Molecular Medicine, vol. 31, no. 4, pp. 797–802, 2013.
View at: Publisher Site | Google Scholar
T. Bertero, Y. Lu, S. Annis et al., “Systems-level regulation of microrna networks by mir-130/301 promotes pulmonary hypertension,” The Journal of Clinical Investigation, vol. 124, no. 8, pp. 3514–3528, 2014.
View at: Publisher Site | Google Scholar
M. R. Paterson and A. J. Kriegel, “Mir-146a/b: a family with shared seeds and different roots,” Physiological Genomics, vol. 49, no. 4, pp. 243–252, 2017.
View at: Publisher Site | Google Scholar
U. Testa, E. Pelosi, G. Castelli, and C. Labbaye, “Mir-146 and mir-155: two key modulators of immune response and tumor development,” Non-Coding RNA, vol. 3, no. 3, p. 22, 2017.
View at: Publisher Site | Google Scholar
Y. Lu, K. L. Hippen, A. L. Lemire et al., “Mir-146b antagomir-treated human tregs acquire increased gvhd inhibitory potency,” Blood, vol. 128, no. 10, pp. 1424–1435, 2016.
View at: Publisher Site | Google Scholar
G. Curtale, M. Mirolo, T. A. Renzi, M. Rossato, F. Bazzoni, and M. Locati, “Negative regulation of toll-like receptor 4 signaling by il-10-dependent microrna-146b,” Proceedings of the National Academy of Sciences of the United States of America, vol. 110, no. 28, pp. 11499–11504, 2013.
View at: Publisher Site | Google Scholar
H. Park, X. Huang, C. Lu, M. S. Cairo, and X. Zhou, “Microrna-146a and microrna-146b regulate human dendritic cell apoptosis and cytokine production by targeting traf6 and irak1 proteins,” The Journal of Biological Chemistry, vol. 290, no. 5, pp. 2831–2841, 2015.
View at: Publisher Site | Google Scholar
E. Legchenko, P. Chouvarine, P. Borchert et al., “PPARγ agonist pioglitazone reverses pulmonary hypertension and prevents right heart failure via fatty acid oxidation,” Science Translational Medicine, vol. 10, no. 438, article eaao0303, 2018.
View at: Publisher Site | Google Scholar
S. J. Hong, S. C. Choi, J. Y. Cho et al., “Pioglitazone increases circulating microrna-24 with decrease in coronary neointimal hyperplasia in type 2 diabetic patients- optical coherence tomography analysis,” Circulation Journal, vol. 79, no. 4, pp. 880–888, 2015.
View at: Publisher Site | Google Scholar
K. M. Grebe, K. Takeda, H. D. Hickman et al., “Cutting edge: sympathetic nervous system increases proinflammatory cytokines and exacerbates influenza a virus pathogenesis,” The Journal of Immunology, vol. 184, no. 2, pp. 540–544, 2010.
View at: Publisher Site | Google Scholar
K. M. Grebe, H. D. Hickman, K. R. Irvine, K. Takeda, J. R. Bennink, and J. W. Yewdell, “Sympathetic nervous system control of anti-influenza cd8+ t cell responses,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 13, pp. 5300–5305, 2009.
View at: Publisher Site | Google Scholar
B. Kruszewska, S. Y. Felten, and J. A. Moynihan, “Alterations in cytokine and antibody production following chemical sympathectomy in two strains of mice,” The Journal of Immunology, vol. 155, pp. 4613–4620, 1995.
View at: Google Scholar
E. J. Seeley, S. S. Barry, S. Narala, M. A. Matthay, and P. J. Wolters, “Noradrenergic neurons regulate monocyte trafficking and mortality during gram-negative peritonitis in mice,” The Journal of Immunology, vol. 190, no. 9, pp. 4717–4724, 2013.
View at: Publisher Site | Google Scholar
A. M. L. Coenen-Stass, I. Magen, T. Brooks et al., “Evaluation of methodologies for microrna biomarker detection by next generation sequencing,” RNA Biology, vol. 15, no. 8, pp. 1133–1145, 2018.
View at: Publisher Site | Google Scholar
J. Kufel and P. Grzechnik, “Small nucleolar rnas tell a different tale,” Trends in Genetics, vol. 35, no. 2, pp. 104–117, 2019.
View at: Publisher Site | Google Scholar
B. E. Jady and T. Kiss, “Characterisation of the u83 and u84 small nucleolar rnas: two novel 2'-o-ribose methylation guide rnas that lack complementarities to ribosomal rnas,” Nucleic Acids Research, vol. 28, no. 6, pp. 1348–1354, 2000.
View at: Publisher Site | Google Scholar
A. Huttenhofer, M. Kiefmann, S. Meier-Ewert et al., “Rnomics: an experimental approach that identifies 201 candidates for novel, small, non-messenger rnas in mouse,” The EMBO Journal, vol. 20, no. 11, pp. 2943–2953, 2001.
View at: Publisher Site | Google Scholar
L. Lestrade and M. J. Weber, “Snorna-lbme-db, a comprehensive database of human h/aca and c/d box snornas,” Nucleic Acids Research, vol. 34, pp. D158–D162, 2006.
View at: Publisher Site | Google Scholar
G. A. Stepanov, J. A. Filippova, A. B. Komissarov, E. V. Kuligina, V. A. Richter, and D. V. Semenov, “Regulatory role of small nucleolar rnas in human diseases,” BioMed Research International, vol. 2015, Article ID 206849, 10 pages, 2015.
View at: Publisher Site | Google Scholar
A. M. Kiss, B. E. Jady, E. Bertrand, and T. Kiss, “Human box h/aca pseudouridylation guide rna machinery,” Molecular and Cellular Biology, vol. 24, no. 13, pp. 5797–5807, 2004.
View at: Publisher Site | Google Scholar
J. H. García, M. D. Osuna, F. M. Castrejon, L. G. Enriquez, P. A. Reyes, and J. J. C. Hermosillo, “Methods to detect antifibrillarin antibodies in patients with systemic sclerosis (SSc): a comparison,” Journal of Clinical Laboratory Analysis, vol. 18, no. 1, pp. 19–26, 2004.
View at: Publisher Site | Google Scholar
J. M. Metzger, H. N. Matsoff, A. D. Zinnen et al., “Post mortem evaluation of inflammation, oxidative stress, and PPARγ activation in a nonhuman primate model of cardiac sympathetic neurodegeneration,” PLoS One, vol. 15, no. 1, article e0226999, 2020.
View at: Publisher Site | Google Scholar
H. Shin, C. P. Shannon, N. Fishbane et al., “Variation in rna-seq transcriptome profiles of peripheral whole blood from healthy individuals with and without globin depletion,” PLoS One, vol. 9, no. 3, article e91041, 2014.
View at: Publisher Site | Google Scholar
K. Krjutškov, M. Koel, A. M. Roost et al., “Globin mrna reduction for whole-blood transcriptome sequencing,” Scientific Reports, vol. 6, no. 1, article 31584, 2016.
View at: Publisher Site | Google Scholar
N. Bayatti, J. Cooper-Knock, J. J. Bury et al., “Comparison of blood rna extraction methods used for gene expression profiling in amyotrophic lateral sclerosis,” PLoS One, vol. 9, no. 1, article e87508, 2014.
View at: Publisher Site | Google Scholar
S. H. Kam, A. Singh, J. Q. He et al., “Peripheral blood gene expression changes during allergen inhalation challenge in atopic asthmatic individuals,” The Journal of Asthma, vol. 49, no. 3, pp. 219–226, 2012.
View at: Publisher Site | Google Scholar
T. Nikula, J. Mykkänen, O. Simell, and R. Lahesmaa, “Genome-wide comparison of two rna-stabilizing reagents for transcriptional profiling of peripheral blood,” Translational research : the journal of laboratory and clinical medicine, vol. 161, no. 3, pp. 181–188, 2013.
View at: Publisher Site | Google Scholar
H. Koyama, S. Tanaka, M. Monden et al., “Comparison of effects of pioglitazone and glimepiride on plasma soluble rage and rage expression in peripheral mononuclear cells in type 2 diabetes: randomized controlled trial (piorage),” Atherosclerosis, vol. 234, no. 2, pp. 329–334, 2014.
View at: Publisher Site | Google Scholar
W. Zhao, C. C. Berthier, E. E. Lewis, W. J. McCune, M. Kretzler, and M. J. Kaplan, “The peroxisome-proliferator activated receptor-γ agonist pioglitazone modulates aberrant t cell responses in systemic lupus erythematosus,” Clinical immunology, vol. 149, no. 1, pp. 119–132, 2013.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Jeanette M. Metzger et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

463

Downloads

1303

Citations

BioMed Research International

Effects of Cardiac Sympathetic Neurodegeneration and PPARγ Activation on Rhesus Macaque Whole Blood miRNA and mRNA Expression Profiles

Abstract

1. Introduction

2. Materials and Methods

2.1. Ethics Statement

2.2. Subjects and Whole Blood Collection

2.3. RNA Isolation

2.4. miRNA Library Preparation and Next-Generation Sequencing (NGS)

2.5. mRNA RT2 Profiler PCR Array and Additional Gene Reverse Transcription-Quantitative PCR (RT-qPCR)

2.6. Data Analysis

3. Results and Discussion

3.1. NGS Overview

3.2. miRNA

3.2.1. miRNA Expression in Rhesus Macaque Whole Blood

3.2.2. miRNA Differential Expression

3.3. Other Small RNAs

3.3.1. Other Small RNA Expression in Rhesus Macaque Whole Blood

3.3.2. Other Small RNA Differential Expression

3.4. mRNA

4. Conclusions

Data Availability

Conflicts of Interest

Acknowledgments

Supplementary Materials

References

Copyright

2.5. mRNA RT² Profiler PCR Array and Additional Gene Reverse Transcription-Quantitative PCR (RT-qPCR)