Participants. This study was conducted at the Laureate Institute for Brain Research (LIBR) in Tulsa, Oklahoma (OK) between 6/30/2015 and 10/30/2015. The study was approved by the Western Institutional Review Board, and all experiments were performed in accordance with the Declaration of Helsinki; Informed consent was obtained from all participants. Participants were recruited from the general community through newspaper, flyer, radio and other media advertisements in Tulsa and the surrounding regions of OK. Subjects were screened by trained clinical interviewers to evaluate the following study exclusion criteria: (1) history of any mental health disorder such as dysthymia, simple phobia, MDD, obsessive compulsive disorder or panic disorder as a primary diagnosis currently or within 6 months prior to the screening visit; (2) history of schizophrenia, schizoaffective disorder, or a bipolar disorder; (3) current or past 6-month alcohol or drug abuse; (4) regular use (> 15 days for past 30 days) of NSAIDS; (5) history of clinically significant hepatic cardiac, renal, neurologic, cerebrovascular, metabolic, gastric, or pulmonary disease; (6) past-year use of psychotropic drugs or antidepressants; (7) history of seizure disorders (except for childhood febrile seizures); (8) serious suicidal ideation or behavior; (9) women currently pregnant or planning to become pregnant within the next 18 weeks; (10) women currently menstruating; (11) claustrophobia, or phobia for injections or blood; and (12) fMRI-related exclusion criteria (e.g., medications treating cardiovascular, respiratory, endocrine and neurological diseases likely to influence cerebral blood flow). The trial was stopped after completing target recruitment. Twenty subjects (10 females,10 males; mean age = 32 years, SD = 7, range = 27 to 51; mean body mass index [kg/m2] = 27, SD = 6, range = 20.4 to 44.7) completed this study (Table 1) (two additional subjects withdrew prior to study completion). See supplemental material for complete consort diagram. The sample size was determined based on prior ph-fMRI studies (e.g., Aupperle et al., 2012 (50)).
Table 1
Variable | Training Dataset | Ibuprofen Dataset | Cohen’s d or χ2 | p value |
N | 480 | 20 | | |
Age (Years) | 34.5 ± 10.5 | 32.4 ± 6.7 | 0.274 | < 0.01 |
Sex = Male (%) | 217 (45.2) | 10 (50.0) | 0.04 | 0.85 |
Body Mass Index (kg/m2) | 26.5 ± 5.5 | 26.6 ± 5.8 | 0.12 | 0.37 |
Ethnicity | | | 0.08 | 0.77 |
Hispanic or Latino | 27 | 2 | | |
Not Hispanic or Latino | 437 | 18 | | |
Refused/do not know | 4/14 | 0/0 | | |
Educational Status | | | 0.19 | 0.56 |
Grandes 1–11 | 1 | 0 | | |
12th grade (no diploma) | 6 | 0 | | |
Regular high school diploma | 40 | 0 | | |
GED or alternative credential | 229 | 0 | | |
Some college, no degree | 142 | 3 | | |
Associate degree | 51 | 11 | | |
Bachelor’s degree | 1 | 2 | | |
Master’s degree | 0 | 2 | | |
Missing data | 10 | 2 | | |
Note, Cohen’s d or χ2 were computed from a two-sample Kolmogorov-Simimov or a χ2 test with the p values. |
Study Design. This double-blind, randomized, repeated-measures study was registered on clinicaltrials.gov (Identifier: NCT02507219, Study of Ibuprofen Effects on Brain Function, first posted date: 23/07/2015). After the screening visit (T0), eligible subjects were tested three times (T1, T2 and T3) and at each visit received placebo, 200 mg or 600 mg dose of oral ibuprofen (dose order was counterbalanced across subjects). Ibuprofen and visually identical placebo capsules were produced by a local compounding pharmacy in Tulsa, OK. The random allocation sequence was generated using a random number generator by a statistician not involved in data collection. Drugs were labelled A/B/C and all study personnel were blinded until after data collection was complete. On the day of sessions T1, T2 and T3, subjects fasted overnight and arrived in the morning and received a snack along with either placebo, 200mg or 600mg of ibuprofen. Subjects underwent a fMRI scan approximately 1 hour after dosing and a blood draw approximately 5 hours after drug administration.
Blood samples. Venous blood was collected in BD Vacutainer Serum Blood Collection tubes with spray-coated silica as a clot activator and then transported to the University of Oklahoma Integrative Immunology Center (IIC) within two hours of collection. Blood tubes were centrifuged at 1300 x g for 10 minutes at room temperature, serum was removed, aliquoted, and then stored at − 80°C until analysis.
Total exosomes (TEs) isolation. Total exosomes (TEs) were isolated from 250 microliters (µL) serum samples using 63 µL of ExoQuick exosome precipitation solution (System Biosciences, CA, United States; Catalog #EXOQ5A-1). TE pellets were resuspended in 300 µL of 1X phosphate buffered saline (PBS) (Thermo Fisher Scientific, United States; Catalog #AM9625) with Halt protease and EDTA-free phosphatase inhibitor cocktail (Thermo Fisher Scientific, United States; Catalog #78425). TEs were used immediately or stored at − 80°C until immunochemical enrichment of exosomes from neural sources could occur.
NEEs isolation. TEs were enriched by a magnetic streptavidin bead immunocapture kit against the neural adhesion marker, L1CAM (CD171) biotinylated antibody (Fig. 4A). This technology to enrich NEEs in blood samples has been previously validated61–63. The CD171 (L1CAM, neural adhesion protein) marker was used for NEE enrichment due to its high and relatively specific expression in neurons and low levels of expression in many other cell types (neuronal marker assessments show that majority of the exosomes in NEEs have a neuronal origin61, and the level of neuronal markers neurofilament-light (NF-L) and synaptophysin (SYP) are significantly enriched in NEE by 86-fold and 951-fold compared to TEs64). Briefly, 80 µL of 9.1 µm diameter covalently cross-linked streptavidin magnetic beads (System Biosciences, CA, United States; Catalog #CSFLOWBASICA-1) and 80 µL of 100 nanograms/µL of mouse anti-human CD171 biotinylated antibody (clone 5G3, eBioscience, United States; Catalog #13-1719-82) were incubated on ice for 2 hours with gentle flicking every 30 minutes. After washing three times in 1X Bead Wash Buffer (BWB) (Systems Biosciences, CA, United States; Catalog #CSFLOWBASICA-1) using a magnetic stand, the bead/antibody complex was suspended in 400 µL of BWB. 200 µL of TE suspensions were added to the bead/antibody complex and incubated overnight at 4°C with rotation. After confirmation by flow cytometry, NEEs were eluted from the beads using 300 µL of Exosome Elution Buffer (System Biosciences, CA, United States; Catalog #CSFLOWBASICA-1), and NEEs were used immediately or stored at -80°C for miRNAs purification.
Flow Cytometry. Once the NEEs were captured and stabilized, the bead/antibody/exosome complex was coupled to a fluorescein isothiocyanate (FITC) fluorescent tag that specifically binds to exosomes (Exo-FITC, Systems Biosciences; Cat #CSFLOWBASICA-1) and subsequently analyzed by flow cytometry to confirm exosome capture (Fig. 4A). Briefly, the bead/antibody/exosome complex was washed three times with 1X BWB and then incubated in 240 µL of Exosome Stain Buffer and 10 µL of EXO-FITC for 2 hours on ice with gentle flicking every 30 minutes. The stained complex was washed three times in 1X BWB and suspended in 500 µL 1X BWB prior to flow cytometry loading. The flow cytometric data were acquired using a BD LSR II Special Order Flow Cytometer (BD Biosciences, San Jose, CA). Instrument performance was validated using BDTM Cytometer Setup and Tracking (CS&T) beads (BD Biosciences, San Jose, CA). All data were analyzed using FACS DIVA 8.0 software (BD Biosciences). Debris and small particles were excluded by gating out events with low forward scatter. Figure 4B shows an example of successful exosome capture using the beads coated with CD171 antibodies specific for NEEs. Beads without exosomes were used as a negative control.
Western Blot. Enrichment of NEEs from TEs was confirmed by western blot (Fig. 4C). Protein concentrations in TEs and NEEs and exosome-depleted serum were determined using a Pierce™ BCA protein assay (Thermo Fisher Scientific, USA, Catalog # 23225). Samples were denatured directly in a 4X Laemmli sample loading buffer and separated by SDS-PAGE using Mini PROTEAN® TGX™ precast gels (Bio-Rad, USA, Catalog # 4561044). Following electrophoresis, gels were transferred unto polyvinylidene difluoride (PVDF) membranes using a Trans-Blot® Turbo transfer system (Bio-Rad, USA, Catalog # 1704156). PVDF membranes were blocked with 5% non-fat milk powder in Tris-buffered saline containing 0.1% Tween20 (Bio-Rad, USA, Catalog # 1706435) and then incubated with primary mouse antibody against CD171 (1:1000, CD171 Monoclonal Antibody (eBio5G3 (5G3), eBioscience™, Catalog # 13-1719-82) overnight at 4°C. This was followed by incubation with a horseradish peroxidase-conjugated anti-mouse antibody (1:2000, Cell Signaling, Catalog # 7076S) for 1 hour at room temperature. PVDF membranes were visualized by Clarity Max Western ECL Substrate (Bio-Rad, USA, Catalog # 1705062) and imaged using ImageQuant LAS 4000 (GE Healthcare Bio-Science, Sweden). Image displayed in Fig. 4C was from the same gel, full-length gel is included in Supplemental Figure S1.
Nanoparticle analysis. Size and concentration of NEEs were determined using Nanoparticle Tracking Analysis system (NanoSight NS300, Malvern Panalytical Inc., Malvern, United Kingdom). Figure 4D showed that the majority of captured NEEs were in the exosome size range, with mean size of 175 nm, and a standard deviation of 53 nm; the average concentration of NEEs was approximately 1.6 x 108 particles per mL.
NEEs MiR Purification. Purification was conducted using a Qiagen miRNeasy Micro Kit (QIAGEN, United States) according to the manufacturer’s protocol. Small RNA concentration was measured using an Agilent Small RNA kit (Agilent, United States) on a Bioanalyzer 2100 instrument (Agilent, United States). MiR samples were stored at -80°C until sequencing.
NEEs MiR sequencing and data processing. MiR samples were sent to the Oklahoma Medical Research Foundation (OMRF) Clinical Genomics Center for Next Generation Sequencing (NGS). MiR libraries were generated with a Qiagen QIAseq MiR library preparation kit and NGS was performed on an Illumina NextSeq HO SR75. Raw sequence FASTQ files received from OMRF were imported to Partek Flow software for data analysis. Adapters from 3’ end were trimmed from the raw read after a quality check, and then aligned to the human genome hg38 using Bowtie alignment. Next, the aligned reads were quantified against the human miRbase mature microRNAs version 22 and reads from miR genes were normalized and scaled to reads per million for statistical data analysis.
MID task. This task contained trials where participants saw a cue then a target, and the objective was to press a button as quickly as possible while the target was on the screen. Cues indicated the possible outcomes of a trial, with circle cues indicating a gain for hitting the target and square cues indicating a loss for missing the target. The magnitude of potential gain/loss was indicated by the position of a line on the cue and text showing the trial type (-5/-1/-0/+0/+1/+5). There were 90 trials (15 of each condition) split across two 568s runs. Target duration was calibrated based on a practice session completed before the scan and adjusted during scanning, so that on average participants hit on 60% of trials and earned $30 for the task.
fMRI data were acquired during the MID task using two identical GE MR750 3T scanners using echo-planar imaging and the following parameters: 39 axial slices, TR/TE = 2000/27ms, FOV/slice = 240/2.9mm, 128 x 128 matrix. High-resolution structural images were obtained through a 3D axial T1-weighted magnetization-prepared rapid acquisition with gradient echo (MP-RAGE) sequence (TR/TE = 5/2.0 12 ms, FOV/slice = 240 × 192/0.9mm, 186 axial slices). fMRI preprocessing was done using Analysis of Functional Neuro Imaging (AFNI)65 and consisted of: removal of the first three EPI volumes for signal stabilization, despiking, slice timing correction, co-registration to the anatomical image, motion correction via rigid-body alignment, and normalization to the Montreal Neurological Institute (MNI) standard space while resampling to 2x2x2-mm voxels, and smoothing with a 4-mm full-width at half-maximum filter. A general linear model was used to model the BOLD response during the anticipation phase of the MID with regressors for each of the 6 conditions (-5/-1/-0/+0/+1/+5) as well as the six motion parameters and four polynomial terms. Voxelwise beta coefficients representing percent signal change were taken to the group level.
Statistical analysis on NEEs MiR. Normalized MiR genes miR-27b-3p and miR-320b were log-transformed due to non-normality and used as the dependent variable in a repeated measure analysis of variance (ANOVA) with dose (placebo, ibuprofen 200 mg, ibuprofen 600 mg) as the within-subjects variable; paired t-tests were employed to test mean differences between doses. Similar repeated ANOVA tests and paired t-tests were also estimated for other MiRs with enough reads for statistical analysis.
MID analyses. 3dLME66 was used to fit models with beta ~ miR*drug + visit (T1/T2/T3) + condition for gains and losses separately. Random effects of subject and visit nested within subject were included. Four models were run, for gains/losses and miR-27b-3p/miR-320b. After fitting each model, smoothness of the residuals was estimated using 3dFWHMx. Then 3dClustSim was used to estimate the family-wise error rate (FWER) given voxel-wise and cluster-size thresholds. Results are reported with a voxel-wise threshold of p < 0.001 and a FWER of α < 0.01. Effect sizes are reported based on the same LME models run post-hoc on average percent signal change in significant clusters.