Dynamics of neurotransmitter and extracellular vesicle-derived microRNA landscapes during heroin and methamphetamine withdrawal

Ethics statement and clinical sample collection

In this cross-sectional study, a total of 120 male patients with SUDs were recruited from a joint program for drug detoxification and rehabilitation in the First Affiliated Hospital of Kunming Medical University and the Kunming Drug Rehabilitation Center from January 2018 to October 2019. Twenty age-matched healthy male volunteers were recruited from the Physical Examination Center in the First Affiliated Hospital of Kunming Medical University. Prior to inclusion in the study, general demographic data and complete medical, psychiatric, and substance abuse histories for all participants were collected for all participants. Written informed consent for participation was provided by all individuals. The study was approved by the First Affiliated Hospital of Kunming Medical University Ethics Committee (2018-L-42). Each patient was diagnosed as having opioid (e.g., heroin) or amphetamine (e.g., methamphetamine) dependence based on the Diagnostic and Statistical Manual of Mental Disorders, 5th edition (DSM-5) criteria provided by the American Psychiatric Association. Individuals were excluded if they fulfilled any of the following conditions: any other DSM-5 axis I or II disorders (other than amphetamine/opioid dependence); HIV antibodies; neurological disorders or serious medical conditions; polysubstance use. Of the total 140 participants recruited, there were 60 heroin-dependent patients (HDPs, 20:20:20 at the 7-day, 3-month, and 12-month withdrawal stage), 60 methamphetamine-dependent patients (MDPs, 20:20:20 at the 7-day, 3-month, and 12-month withdrawal stage), and 20 healthy controls (HCs). An initial 10 participants per each group were used as the discovery cohort, and the subsequent 10 participants per each group were used as the validation cohort (Figure 1A). Peripheral blood samples were collected in 5 ml EDTA□2Na vacuum tubes between 8:00 and 10:00 AM in order to minimize a possible rhythm variance bias and the clinical data were obtained at the same time. Within 1 hour of blood collection, plasma was separated out by centrifugation and was immediately aliquoted and frozen in liquid nitrogen. The stored plasma samples were transported and packaged on dry ice.

• Download figure
• Open in new tab

Figure 1: Pipeline showing the experimental design in this study.

(A-C) An overview of the experimental design. (D-E) Characterisation of exosomes derived from the plasma of Substance use disorder (SUD) patients and controls; (D) The shape and structure of plasma exosomes isolated by ultracentrifugation under transmission electron microscopy (TEM); (E) The size of exosomes was analysed by nanoparticle tracking analysis (NTA); (F) Western blots of exosomal membrane markers, including Alix, CD63, and CD81. The exosomes isolated from each participant underwent both TEM, NTA and Western blotting to confirm successful enrichment and ensure limited contamination with cellular debris prior to small RNA deep sequencing.

Scales

The Hamilton-Anxiety scale (HAM-A) is a semi-structured questionnaire administered by the interviewer, who assesses 14 clusters of questions, each concerning a series of symptoms. Seven elements examine psychological stress and seven examine physical stress with a satisfactory inter-rater reliability (Hamilton, 1959). The Hamilton-Depression Scale (HAM-D) is a semi-structured questionnaire administered by the interviewer, who assesses 21 clusters of questions, measuring depression symptoms (Hamilton, 1960).

Sample preparation and instrumentation

Samples were thawed on ice-bath and 25 μl of each plasma sample was transferred to a 96-well plate. Then the plate was automated by Biomek 4000 workstation (Beckman Coulter, Inc., Brea, USA) following preprogrammed protocol. Briefly, 100 μl ice cold methanol with partial internal standards was automatically added to each sample and vortexed vigorously for 5 min. The plate was centrifuged at 4000g for 30 min (Allegra X-15R, Beckman Coulter, Inc., Indianapolis, USA). Next, 30 μl of supernatant was transferred to a clean 96-well plate, and 20 μl of freshly prepared derivative reagents was added to each well. The plate was sealed, and the derivatization was carried out at 30°C for 60 min. After derivatization, 350 μl of ice-cold 50% methanol solution was added to dilute the sample. Then the plate was stored at - 20°C for 20 min and followed by 4000g centrifugation at 4°C for 30 min. An ultra-performance liquid chromatography coupled to tandem mass spectrometry (UPLC-MS/MS) system (ACQUITY UPLC-Xevo TQ-S, Waters Corp., Milford, USA) was used to quantitate the metabolite in this project (Figure 1B).

Sample processing

The rapid turnover of many intracellular metabolites makes immediate metabolism quenching necessary. The extraction solvents were stored in -20°C freezer overnight and added to the samples immediately after the samples were thawed. An ice-salt bath was used to keep samples at a low temperature to minimize degradation during sample preparation. All the prepared samples were analyzed within 48 hours after extraction and derivatization. A comprehensive set of rigorous quality control/assurance (QC/QA) procedures was employed to ensure a consistently high quality of analytical results, throughout controlling every single step from sample receipt at laboratory to final deliverables. The goal of QC/QA was to provide the reliable data for biomarker discovery study. To achieve this, three types of quality control samples (test mixtures, internal standards, and pooled biological samples) were used in the metabolomics platform. In addition to the quality controls, conditioning samples, and solvent blank samples were also used for obtaining optimal instrument performance. Each sample was accessioned into Metabo-Profile LIMS system and was assigned by the LIMS a unique identifier, which was used to track all sample handling, tasks, results and other conditions. The samples and aliquots were barcoded and analyzed by the LIMS system.

Metabonomic data analysis

Data were analyzed using a XploreMET v3.0 system (Metabo-Profile, Shanghai, China). Partial least squares projection to latent structures and discriminant analysis (PLS-DA) were used to process the acquired GC–TOF/MS data. Principal component analysis (PCA) and partial least square discriminant analysis (PLS-DA) were performed for both positive and negative models after log transformation and pareto scaling. The variable importance in the projection value of each variable in the PLS-DA model was calculated to indicate its contribution to the classification. The significance of each metabolite with the VIP value >1 was assessed by Student’s t-test at univariate level, and the results were adjusted for multiple testing using the Benjamini–Hochberg procedure with the critical false discovery rate (FDR) set to 0.05.

Enrichment of plasma derived exosomes from peripheral blood

Fresh whole blood samples from recruited patients or HCs were collected into the EDTA□2Na vacuum tubes and the blood was centrifuged at 1500 g for 15 min. The plasma was transferred into a new tube and centrifuged at 20,000 g at 4 °C for 15 min. The supernatant was collected and kept at −80 °C until use. The exosomes were isolated by size exclusion chromatography method as described previously with minor modifications (Boing et al, 2014). In brief, 2 ml of 0.8 μm-filtered plasma sample was purified using Exosupur® columns (Echobiotech, China). The exosome samples were eluted with PBS and a total of 2 mL eluate fractions were collected according to the manufacturer’s instruction. Fractions were concentrated to 200 μl by 100 kDa molecular weight cut-off Amicon® Ultra spin filters (Merck, Germany).

Transmission electron microscopy and nanoparticle tracking analysis

For transmission electron microscopy analysis (Varga et al, 2014), 10 μl of exosomal suspension was loaded on a continuous carbon grid, and incubated for 10 min at room temperature, followed by distilled water washing. After being absorbed excessively, each sample was negatively stained with 10 μl of 2% (w/v) uranyl acetate solution at room temperature for 1 min and sucked up residue by filter paper. Then, the samples were photographed at 80kv by a Hitachi H-7650 transmission electron microscope (Hitachi, Tokyo, Japan). Particle size and concentration of plasma-derived EVs were analyzed by nanoparticle tracking analysis (NTA) on a Nano Sight NS300 instrument (Malvern Panalytical, Malvern, UK) following the manufacturer’s protocol (Varga et al., 2014). In brief, the samples were diluted to 1×10⁷ ∼ 1×10⁹ particles/ml in PBS, and slowly injected into the sample pool and detected in flow mode. These flow measurements consisted of three measurements of 60 s. The captured data were analyzed using NTA 3.2 software.

Library preparation and small RNA sequencing

Total RNAs from exosomes were isolated using the miRNeasy® Mini kit (Qiagen, MD, USA) according to the manufacturer’s protocol. Sequencing libraries were generated using the QIAseq miRNA Library Kit (Qiagen, MD, USA), according to the manufacturer’s instructions. 500 ng of exosome derived total RNA for each sample was used as input for library preparation. The RNA samples were barcoded by ligation with unique adaptor sequences to allow pooling of samples. The elution containing pooled DNA library was further processed for cluster generation and sequencing using a NovaSeq 6000 platform.

Real-time PCR analysis

To validate the miRNAs identified in RNA sequencing data, we performed qPCR analysis of selected miRNA targets. Total RNA extraction from the serum exosomes was performed with the miRNeasy® Mini kit (Qiagen, MD, USA) according to the manufacturer’s instructions. The concentration and quality of the RNA was determined using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 System (Agilent Technologies, CA, USA). Single strand cDNA was synthesized with the PrimeScript™ RT reagent Kit (Takara, Dalian, China). The RT-PCRs were conducted on an ABI Prism 7500 HT using the TaqMan Fast Advanced Master Mix. Real-time PCR amplification was performed in triplicate and a negative control was included for each primer. miRNA levels were normalized to average Ct of internal control has-miR-106a-5p in the confirmation phase.

Analysis of miRNA-predicted target genes and network visualization

The miRNet (http://www.mirnet.ca) was used to predict mRNAs targets, conduct functional enrichment analysis and generate network visualization (Fan et al, 2016). The miRNA-target interaction data were based on 10 well-annotated databases including miRTarBase, TarBase, miRecords, SM2miR, Pharmaco-miR, miR2Disease, PhenomiR, StarBase, EpimiR, miRDB, in which the datasets were merged on the same targets for each species, removed the duplicate entries, updated the literature references and stored them into SQLite database for fast retrieval. The conventional hypergeometric tests were implemented in miRNet for functional annotations-based Disease, GO, or KEGG pathway analysis.

Stem cell culture and neurite outgrowth assay

Human induced pluripotent stem cells (hiPSCs) were purchased from Guidon Pharmaceutics (Beijing, China) and cultured in mTeSR™ Plus media (STEMCELL Technologies) for feeder-free maintenance and expansion. For neurite outgrowth assay, 0.7 x 10⁵ iPSCs/well were first seeded in a 24-well plate and neuronal differentiation from hiPSCs was conducted according to the protocol previous published (Zhang et al, 2013). Ectopic expression of targeted miRNA (pEZX-MR04 or pEZX-MR04-hsa-744-5p) via lentiviral transfection (250μl media + 50μl lentivirus + 4μg/ml polybrene) for 24 hours was applied on day 4 of neuronal differentiation. Three, five and seven days after lentiviral induction of targeted miRNA expression and passage to a new dish, hiPSCs derived neurons were imaged at 20x magnification using CellInsight™ CX5 High Content Screening Platform (Thermo Scientific). Individual cell measurements of mean/median/maximum process length, total neurite outgrowth (sum of the length of all processes), number of processes, number of branches, cell body area, and mean outgrowth intensity were obtained using the HCS Studio™ 2.0 Cell Analysis Software (Thermo Scientific). Total RNA from cells was isolated using miRNeasy® Mini kit (Qiagen, MD, USA) according to the manufacturer’s instructions.

DiR-labelled exosomes and in vivo biodistribution assay

Female C57BL/6 mice (aged 6 to 8 weeks old) were kept under standard laboratory conditions. Mice were housed in a temperature-controlled room with a 12-hour light/dark cycle. The University of Kunming Medical University Committee on the Use and Care of Animals approved all experiments. The filtered exosome of serum sample from healthy donors was incubated with 1 mM fluorescent lipophilic tracer DiR (1,1-dioctadecyl-3,3,3,3-tetramethylindotricarbocyanine iodide) (D12731, Invitrogen, Life Technologies) at room temperature for 15 minutes prior to exosome isolation. Lipopolysaccharide (LPS) or PBS treated C57BL/6 mice (n = 3/group) were used and freshly purified DiR-labeled exosomes were injected through the tail vein for intravenous injections. The biodistribution of labeled exosomes was examined using 2.0 x 10¹⁰ particles/gram body weight (p/g). The particle count was measured with NTA and the samples were diluted accordingly. For the analysis of DiR-exosome distribution, IVIS Spectrum (Perkin Elmer) was used. Either live (isoflurane sedated) mice were imaged, or the animals were sacrificed and then the organs harvested prior to analysis. For the perfusion experiment, the mice were sedated, and the vascular system was flushed by transcardial perfusion for 5 minutes. The left ventricle was infused with PBS (5 ml/min) and the right atrium was perforated. After perfusion, the organs were harvested and analyzed.

Statistical analysis and bioinformatics analysis

All data were analyzed with R version 3.6.3. Kruskal–Wallis test or analysis of variance (ANOVA) test followed by a post-hoc test (Bonferroni’s t test) was used to test for differences in continuous variables. Comparison between groups was performed using a Student’s t test or a Chi-square analysis as appropriate. The Pearson correlation was used for correlation analyses. For miRNA sequencing, the raw fastq files generated from next-generation sequencing were preprocessed to remove adaptor sequences using the TagCleaner. Subsequently, known and novel miRNAs were detected within the preprocessed data using the miRDeep2 (Friedlander et al, 2012). Statistical analysis and significance were calculated using various tests and adjusted for multiple testing. Adjusted statistical significance was assessed as p < 0.05.

Characteristics of the study population

A total of 120 patients with SUDs and 20 age-matched healthy participants (without SUDs) were recruited at the First Affiliated Hospital of Kunming Medical University. The characteristics of the study participants are presented in Table 1. There were no significant differences for any variables including age, BMI, substance dependence history, annual income, and education level between patients with SUDs and healthy participants in both discovery and validation cohorts (Table 1).

Substance use, depression, and anxiety frequently co-occur (Lai et al, 2015). To better understand whether the use of heroin or methamphetamine influences mood (specifically, anxiety and depression), the HAM-A and HAM-D questionnaires were administered by experienced interviewers. As expected, the HAM-A total scores were significantly increased in the HDPs (p = 0.0005) and the MDPs (p = 0.0014) when compared to HCs in the discovery cohort (Table 1). Significant increases relative to HCs were observed in both HDPs and MDPs at the 3-month and 12-month stages after the initiation of withdrawal, respectively (Table S1). The overall differences in HAM-A for HDPs and MDPs relative to HCs, as well as the differences in the subgroups at the 3-month and 12-month withdrawal stages, were validated (Table S2). In contrast, the significant differences in HAM-D scores observed in the HDPs (p = 0.0134) in the discovery cohort did not persist in the validation cohort. In the subsequent subgroup analysis, the significant increases in HAM-D scores were only validated in the subgroup of MDPs at the 3-month withdrawal stage, but not in any of the subgroups for HDPs (Tables S1 & S2).

Dynamic changes in peripheral neurotransmitters during substance withdrawal

In order to reveal peripheral neurotransmitter changes in the patients with SUDs, we analyzed the classic neurotransmitters and their metabolites, including γ-aminobutyric acid (GABA), glutamate, serotonin, tryptophan, kynurenine, choline, etc, through UPLC-MS/MS (Table S3 and Figure 1B). Plasma samples from the 70 participants in the discovery cohort were analyzed first for the 14 neurotransmitter metabolites. Of these, seven fell below the lower level of quantitation in > 20% of samples and were excluded from analysis. Detailed results for the seven neurotransmitters are presented in Table 2. Our analysis revealed that three of the seven neurotransmitters were associated with significant changes. Two neurotransmitters, GABA and choline, were significantly decreased in both HDPs (fdr.p.GABA = 0.0026, fdr.p.choline = 0.026) and MDPs (fdr.p.GABA = 0.0009, fdr.p.choline = 0.026) relative to HCs, respectively. In contrast, serotonin was significantly increased in HDPs (fdr.p.serotonin = 0.0001) and MDPs (fdr.p.serotonin = 0.0013) relative to HCs (Table 2 and Figure S1). These significant differences in plasma GABA, choline, and serotonin were further validated in the independent cohort (Table S4). Interestingly, there were no significant differences in the tested neurotransmitters between HDPs and MDPs (Table 2 & S4).

Next, to investigate potentially differential neurotransmitters in SUDs across the various withdrawal stages, we performed additional subgroup analyses. Notably, significant increases in serotonin were observed in all three subgroups of HDPs and two subgroups of MDPs when compared to HCs, but not observed in the subgroup of MDPs at the 12-month withdrawal stage (Table 3). In contrast, the significant decrease in GABA persisted in all three subgroups of MDPs relative to HCs but was only found in the HDPs at the 3-month withdrawal stage. Moreover, the significant decreases in choline were only observed in the subgroups of HDPs and MDPs at the 3-month withdrawal stage (Table 3). To validate these findings, the plasma neurotransmitter profiles across the three withdrawal stages were further tested in an independent cohort. Compared to HCs, the significant differences in serotonin at the 7-day withdrawal stage, as well as in GABA, choline, and serotonin at the 3-month withdrawal stage were validated in both HDPs and MDPs, respectively (Table S5).

Associations of neurotransmitters and psychiatric comorbidities in SUD patients

To further determine the associations between neurotransmitters and psychiatric comorbidities in SUDs, Pearson correlations between neurotransmitters (GABA, choline, and serotonin) and HAM-A/D scales were conducted. The plasma GABA in HDPs (R = -0.40, p = 0.01) and MDPs (R = -0.35, p = 0.03) was significantly correlated with HAM-A scores. However, the significant correlation between GABA and HAM-D scores was only observed in MDPs (R = -0.43, p = 0.0052) but not HDPs (R = -0.20, p = 0.21). In addition, significant correlations were found between plasma serotonin and HAM-A score in both HDPs (R = 0.44, p = 0.0054) and MDPs (R = 0.40, p = 0.0012). Nevertheless, no significant correlation was observed between serotonin and the HAM-D score, or between choline and the HAM-A/D scores (Figure S2). In summary, these results indicate some strong correlations between specific plasma neurotransmitters and withdrawal symptoms and suggest that the dynamic changes in peripheral GABA and serotonin may play a key role in influencing moods during withdrawal.

Isolation and characterization of plasma exosomes

Despite research efforts in the field since the discovery of exosomes in peripheral circulation (Raposo & Stoorvogel, 2013), the role of exosomes in the pathophysiology of psychiatric disorders, and especially SUDs, remains poorly understood. In the present study, plasma exosomes were isolated from patients with SUDs and age-matched HCs by size exclusion chromatography in the discovery cohort, and their morphology was verified by transmission electron microscopy (TEM) and nanoparticle tracking analysis (NTA), as shown in Figure 1D & E. These analyses revealed an average particle diameter of 40–140 nm, and there was no significant difference in the exosome size distribution between the patients with SUDs and HCs (data not shown). In addition, the exosomes were verified by detecting the expression of exosomal surface markers CD81 and CD63 (Ricklefs et al, 2019) using western blotting (Figure 1F), which confirmed the successful exaction of exosomes from plasma samples.

Differential expression of exosomal miRNA between SUD patients and healthy controls

To gain insight into the transcriptional dynamics of exosomes, small RNA sequencing was applied to the total RNA derived from plasma exosomes in all 70 participants in the discovery cohort. Sequencing quality statistics are shown in Table S6. Top RNA species detected by sequencing included miRNA, piRNA, tRNA, and snoRNA, and the percentage of reads assigned to each species and the number of mapped miRNAs did not significantly differ between samples (Figure S3). Length distribution analysis of the miRNAs showed that 20–24 nucleotides (nt) was the most frequent read length, with the peak at 22 nt (Table S7), indicating that mature miRNAs were suitably enriched during the sequencing library preparation process. To focus on the miRNA ubiquitously abundant in plasma and suitable as biomarker candidates, we selected the miRNA candidates with an absolute fold change in expression of ≥ 1.5, a nominal FDR value of ≤ 0.05, and the expression (RPKM ≥ 0) across ≥ 50% of samples. This yielded 43 and 26 miRNAs with significant changes in HDPs and MDPs relative to HCs, respectively (Figure 2A-D, Table S8 & S9). The Venn diagram displays a list of 20 shared differentially expressed miRNAs (DE-miRNAs) in HDPs and MDPs when compared to HCs (Figure 2E, Table 4), suggesting there might be a variety of shared, exosomal miRNA-associated, substance-related functional abnormalities between HDPs and MDPs.

• Download figure
• Open in new tab

Figure 2: Distinct exosomal miRNAs in HC, HDPs and MDPs.

Small RNA sequencing was performed to reveal the small RNA content in plasma samples collected from HC and patients with SUDs. (A-D) Statistically significant (Fold change ≤ −1.5 and ≥ 1.5, p.fdr ≤ 0.05) tRNA fragments were detected in (A, B) HDPs vs HC and (C, D) MDPs vs HC. (E) Venn diagram showing the number of the DE-miRNAs identified in HDPs and MDPs compared to HC. (F, G) Predicted functional pathways and biological functions (F) associated with mRNA targets of the 20 overlapped DE-miRNAs and (G) miRNA-mRNA interacting network. (H) ROC AUC for the diagnosis of patients with SUDs compared to HC. (I) Expression of miR-451a in HC, HDPs and MDPs. (J) ROC AUC for the diagnosis of MDPs compared to HDPs. (K) Expression of miR-21 in HC, HDPs and MDPs. The range of expression values in the groups is indicated by the error bars.

Pathway analysis and gene target and gene ontology prediction

Network visualization and functional enrichment analysis is a powerful combination in delivering important biological insights. Using miRNet, a network-based visual exploration and functional interpretation platform (Chang et al, 2020), we identified predicted functional pathways associated with mRNA targets of the 20 shared DE-miRNAs and thereafter generated an miRNA-mRNA interacting network (Figure 2E-G). Top associated pathway/disease hits included mental and motor retardation (p = 5.2e-24), cognitive delay (p = 4.2e-23), global developmental delay (p = 3.8e-22), with a strong overall theme of developmental or intellectual abnormality (Figure 2F-G). Next, the differential metabolites and miRNAs were combined and used to perform ROC analysis to identify an optimal metabolite and/or exosomal miRNA biomarker panel in plasma for SUDs. Strikingly, from the ROC curves, the calculated sensitivities and specificities of the miRNA in diagnosing SUDs and distinguishing HDPs and MDPs were identified (Figure 2H-K). Notably, the analysis of hsa-mia-451a and hsa-mir-21a resulted in an AUC of 0.966 and 0.861, suggesting that these plasma exosomal miRNA, especially has-miR-451a and has-miR-21-5p, could serve as additional indicators for monitoring substance dependence and withdrawal. In order to verify the accuracy of miRNA-seq, hsa-miR-451a and hsa-miR-21-5p were selected as mentioned above for RT-qPCR analysis. The results indicated that the expression patterns of hsa-miR-451a and hsa-miR-21-5p were similar between small RNA-Seq and RT-qPCR (Figure S4).

Identification of exosomal miRNA signatures across heroin and methamphetamine withdrawal

We further investigated the changes in exosomal miRNAs in SUDs across the three withdrawal stages. When small RNA-seq data were subjected to weighted gene co-expression network analysis (WGCNA) (Langfelder & Horvath, 2008), multiple gene-network modules were obtained. We decided to focus mainly on the modules that could represent the dynamic changes of miRNA expression throughout withdrawal stages. Analysis of the module-trait relationship identified six unique modules (blue, brown, magenta, purple, pink, and cyan) correlated with the type of substance use or the stage of withdrawal (p < 0.05). The six modules represent different miRNA expression alteration patterns (Figure 3A-B). Hub gene/miRNA network analysis of the modules revealed that a hierarchical organization of highly connected gene/miRNA in each module, through key controlling genes/miRNAs in the modular network could be identified. Notably, the blue module miRNAs (e.g., hsa-miR-140-3p) were mainly elevated at the 3-month withdrawal stage for both HDPs and MDPs (Figure 3C-D). The purple module consisting of miRNAs significantly decreased at all three withdrawal stages (e.g., hsa-miR-92a-3p and hsa-miR-451a, Figure 3E-F). Different from the miRNAs in the blue module which showed an increase at the 3-month withdrawal stage, the miRNAs in the brown module were significantly increased at both the 7-day and 3-month stages (e.g., hsa-miR-139-3p, Figure 3G-H). By contrast, the miRNAs in the pink module were significantly decreased at both the 3-month and 12-month stages (e.g., hsa-miR-375-3p, Figure 3I-J). Moreover, the miRNAs in modules magenta and cyan had other interesting but unique expression patterns across withdrawal stages (Figure 3K-N).

• Download figure
• Open in new tab

Figure 3: WGCNA of plasma exosome samples uncovered dynamic changes of miRNA signatures during substance withdrawal.

(A) Hierarchical cluster dendrogram using WGCNA analyzing miRNA sequencing data. A total of 8 modules were identified after 0.25 threshold merging. (B) Module trait correlation analysis revealed that several key modules were significantly correlated with substance use and/or different withdrawal stages. (C, E, G, I, K, M) Eigengene expression of six selected modules across 4 groups of samples. Color code of the modules is preserved. (D, F, H, J, L, N) Expression patterns of key miRNAs representing six selected modules across subgroups of HDPs and MDPs.

We then sought to investigate the relationships between the neurotransmitters (GABA, choline, and serotonin) and the hub gene/miRNA in the blue and brown modules using Spearman’s correlations. One exciting finding was that hsa-miR-140-3p from the blue module showed strong negative correlations (p < 0.01) with both GABA and choline (Figure S5A-B). The hsa-miR-744-5p in the brown module was positively correlated with serotonin in HDPs (p = 4.9e-5) and MDPs (p = 0.032) (Figure S5C). Taken together, the WGCNA and the correlation analyses revealed a dynamic expression pattern of miRNA signatures during both heroin and methamphetamine withdrawal, and a range of neurotransmitter associated exosomal miRNAs in the peripheral blood plasma of patients with SUDs.

Effects of hsa-miR-744-5p on neurite outgrowth in human iPSC derived neurons

Emerging studies have indicated an important role for specific miRNA in neural differentiation, proliferation, and maturation during embryonic development and/or disease progression using neurite outgrowth assays (Agostini et al, 2011; Le et al, 2009). To ascertain the neuronal function and molecular mechanism of hsa-miR-744-5p, we differentiated hiPSCs to neurons and infected them with lentivirus-miR-744-5p or scrambled control (Figure 4A-B). As illustrated in Figure 4C, across three replicate experiments, when iPSC derived neurons were treated with a scrambled control, there was no decrease in relative total outgrowth and cell viability. By contrast, overexpression of hsa-miR-744-5p was sufficient to cause significant decreases in neurite processes (Figure 4D-F). For further characterization, genome-wide gene expression changes upon miRNA transfection were studied using RNA-seq and comparative transcriptome analysis on hiPSCs derived neuronal culture with ectopic hsa-miR-744-5p expression versus scrambled control. One day after the lentivirus infection, a list of 46 genes were differentially expressed with a cut-off (fold change ≤ −1.5 or ≥ 1.5, p ≤ 0.05), in which the differentially expressed genes (DEGs) were mainly over-represented in axon development (p = 2.0e-3), reproductive process (p = 5.5e-3) and system development (p = 7.6e-3) (Figure 4G-H). Likewise, 556 DEGs with the same cut-off were identified on day 3 after the induction of hsa-miR-744-5p expression in hiPSCs derived neuronal culture and these DEGs were mainly enriched for cell cycle arrest (p = 1.1e-3), cell growth (p = 1.8e-3), and positive regulation of neuron apoptotic process (p = 2.4-3) (Figure 4I-J). RNA-seq data may partly demonstrate secondary effects on global mRNA expression, while gene expression changes are triggered by primary miRNA targets. Therefore, to discriminate between primary and secondary target genes, the differential expression results obtained from the RNA-seq analysis were compared to potential hsa-miR-744-5p targets predicted by bioinformatic algorithms. The overlap of statistically significant DEGs and bioinformatically predicted target genes revealed that hsa-miR-744-5p may play a significant role involving neural development and function regulating its primary miRNA targets PSENEN, FDX1L, VPS37D, SPTBN4, PPP5C, and FAM57B (Figure 4K), although their mechanisms require further investigation. Overall, the neurite outgrowth assays and RNA-seq analyses reveal an important role of hsa-miR-744-5p in maintaining neuronal development and function.

• Download figure
• Open in new tab

Figure 4: Effects of hsa-miR-744-5p on neurite outgrowth in human iPSC derived neurons.

(A) Flow diagram of human iPSC derived neuron cell generation. (B) Representative differential interference contrast (DIC), MAP2 (Red) and TUJ1 (Green) fluorescence images illustrating the generation of human iPSC derived neurons cell. (C) Representative images illustrating the effects of hsa-miR-744-5p or control in neurite outgrowth assay on Day3, Day5 and Day7. (D-F) Statistics of neurite outgrowth assay after induction of hsa-miR-744-5p expression in human iPSC derived neurons. All data are presented as mean ± SD; n ≥ 3. p < 0.05, one-way ANOVA with Bonferroni post hoc analysis.