Psikiatri Depan. 2018; 9: 81.
Diterbitkan secara online 2018 Mar 12. doi: 10.3389 / fpsyt.2018.00081
PMCID: PMC5858605
PMID: 29593587
Minho Lee,1,† Hyeyoung Cho,2,3,† Seung Hyun Jung,2,4 Seon-Hee Yim,2 Sung-Min Cho,2,3 Ji-Won Chun,5 Soo-Hyun Paik,5 Taman Yae Eun,6,7 Dong Huey Cheon,7,8 Ji Eun Lee,7 Jung-Seok Choi,9 Dai-Jin Kim,5,* dan Yeun-Jun Chung1,2,3,*
Abstrak
Latar Belakang
Penggunaan adiktif dari Internet dan permainan online adalah gangguan kejiwaan yang potensial yang disebut gangguan permainan Internet (IGD). Perubahan profil ekspresi microRNA (miRNA) telah dilaporkan dalam darah dan jaringan otak pasien dengan gangguan kejiwaan tertentu dan disarankan sebagai biomarker. Namun, belum ada laporan profil miRNA darah di IGD.
metode
Untuk menemukan miRNA terkait IGD, kami menganalisis profil ekspresi miRNA dari sampel 51 (25 IGD dan kontrol 26) menggunakan TaqMan Low Density miRNA Array. Untuk validasi, kami melakukan PCR transkripsi balik kuantitatif dengan sampel independen 36 (kontrol 20 dan 16).
Hasil
Melalui penemuan dan validasi independen, kami mengidentifikasi tiga miRNA (hsa-miR-200c-3p, hsa-miR-26b-5p, hsa-miR-652-3p) yang secara signifikan diregulasi ke bawah dalam kelompok IGD. Individu dengan ketiga perubahan miRNA memiliki risiko IGD yang jauh lebih tinggi daripada mereka yang tidak memiliki perubahan [rasio odds (OR) 22, 95% CI 2.29-211.11], dan OR meningkat dosis secara dependen dengan jumlah miRNA yang diubah. Gen target yang diprediksi dari tiga miRNA dikaitkan dengan jalur saraf. Kami mengeksplorasi ekspresi protein dari tiga gen target hilir oleh western blot dan mengkonfirmasi bahwa ekspresi GABRB2 dan DPYSL2 secara signifikan lebih tinggi pada kelompok IGD.
Kesimpulan
Kami mengamati bahwa ekspresi hsa-miR-200c-3p, hsa-miR-26b-5p, dan hsa-miR-652-3p diturunkan regulasi pada pasien IGD. Hasil kami akan sangat membantu untuk memahami patofisiologi IGD.
Pengantar
Penggunaan adiktif dari Internet dan permainan berbasis Internet bukan hanya sebuah fenomena sosial di negara-negara dengan infrastruktur akses Internet yang luas, tetapi potensi gangguan kejiwaan yang dinamakan Internet gaming disorder (IGD) (1-3). Menurut laporan epidemiologis, tingkat prevalensi IGD pada remaja bervariasi di berbagai negara, mulai dari 0.8 hingga 26.7% (4). Khususnya, penelitian menunjukkan tingkat prevalensi di atas 10% pada remaja di banyak negara Asia seperti Korea Selatan, Cina, Taiwan, Hong Kong, dan Singapura (4). IGD dikaitkan dengan gangguan dalam kognisi, hubungan psiko-sosial, dan kehidupan sehari-hari; misalnya, menurunnya kinerja akademik atau pekerjaan (4-7). IGD sekarang dimasukkan dalam Bagian III (Ketentuan untuk Studi Lebih Lanjut) dari revisi kelima Manual Diagnostik dan Statistik Gangguan Mental (DSM-V) (8). Namun, terlepas dari kepentingan klinis-sosialnya, sedikit yang diketahui tentang mekanisme genetik molekuler di balik IGD.
Studi kembar skala besar baru-baru ini menyarankan latar belakang genetik untuk IGD (9, 10). Vink et al. menyelidiki perbedaan individu dalam penggunaan Internet kompulsif dengan 5,247 kembar monozigot dan diozgotik remaja di Twin Register Belanda dan melaporkan bahwa 48% dari perbedaan dijelaskan oleh faktor genetik (9). Li et al. mengamati pasangan 825 dari kembar remaja Tionghoa dan melaporkan bahwa faktor genetik menjelaskan 58 – 66% dari perbedaan (10). Dengan demikian, polimorfisme gen yang terlibat dalam transmisi neurotransmisi, kognisi, dan perhatian seperti reseptor dopamin gen D2 (DRD2), gen katekolamin-O-metiltransferase (COMT), gen transporter serotonin (5HTTLPR), dan gen reseptor kolinergik nikotin alfa 4 (CHRNA4) telah dilaporkan terkait secara signifikan dengan kecanduan Internet (11-13). Baru-baru ini, Kim et al. varian yang disaring lebih dari gen kandidat 100 terkait dengan produksi, aksi, dan metabolisme neurotransmiter oleh analisis sequencing generasi berikutnya dan melaporkan bahwa rs2229910 dari NTRK3 gen dikaitkan dengan IGD (14).
Selain faktor genetik, juga diketahui bahwa fenotip neurobehavioral secara epigenetik dikendalikan oleh RNA non-coding termasuk microRNAs (miRNAs) (15, 16). miRNA adalah molekul RNA beruntai tunggal kecil non-coding (panjangnya sekitar 20-23 nukleotida), yang secara negatif mengatur ekspresi gen penyandi protein dengan menurunkan mRNA dan memainkan peran penting dalam proses patofisiologis berbagai penyakit (17). Garis bukti telah menunjukkan bahwa miRNA berlimpah dalam sistem saraf pusat manusia (CNS) dan bertindak untuk menyempurnakan tingkat ekspresi gen target mereka, yang terlibat dalam pengembangan dan pematangan sistem SSP (15). Memang, penelitian terbaru telah mengungkapkan bahwa profil ekspresi miRNA diubah dalam jaringan otak pasien dengan gangguan kejiwaan, menunjukkan bahwa profil ekspresi mereka dapat menjadi penanda biologis untuk gangguan kejiwaan (15, 16, 18). Sebagai contoh, melalui analisis postmortem, Lopez et al. melaporkan bahwa ekspresi miR-1202, yang mengatur ekspresi gen reseptor-glutamat metabotropik-4 dan memprediksi respons terhadap antidepresan, diregulasi dalam jaringan korteks prefrontal pada pasien gangguan depresi mayor (19). Dalam hal skrining biomarker, pendekatan ini memiliki batasan yang jelas karena melakukan biopsi jaringan SSP untuk skrining tidak mungkin. Karena miRNA dapat dideteksi dalam darah (plasma atau serum), miRNA yang bersirkulasi memiliki keuntungan yang pasti sebagai biomarker non-invasif pada gangguan neuropsikiatri. Namun, hingga saat ini, belum ada penelitian tentang sirkulasi profil miRNA di IGD. Pemahaman yang lebih baik tentang profil ekspresi miRNA yang beredar dapat membantu untuk memperjelas mekanisme pengembangan IGD dan memfasilitasi terjemahan klinis.
Dalam penelitian ini, kami bertujuan untuk mengidentifikasi penanda miRNA terkait IGD dengan mengamati miRNA plasma yang diekspresikan secara berbeda antara IGD dan kelompok kontrol dan mengeksplorasi implikasi biologis mereka.
Bahan dan Metode
Subjek studi
Kami mensurvei remaja 3,166 (berusia 12 – 18 tahun) menggunakan skor DSM-V IGD. Di antara mereka, 251 (laki-laki 168 dan perempuan 83) didiagnosis sebagai IGD sesuai dengan kriteria DSM-V (8). Total individu 91 (49 IGD dan kontrol 42) memberikan persetujuan untuk penelitian ini. Di antara mereka, empat individu dikeluarkan sesuai dengan kriteria eksklusi. Akhirnya, individu 87 (subjek 45 IGD dan individu kontrol sehat 42) terdaftar untuk penelitian ini. Di antara mereka, peserta 51 (pasien 25 IGD dan kontrol 26) direkrut sebagai penemuan yang ditetapkan dari 2014 ke 2016. Peserta 36 lainnya (pasien 20 IGD dan kontrol 16) direkrut sebagai set validasi independen dari 2016. Semua peserta adalah individu Korea, terdaftar dari Rumah Sakit St. Mary Seoul (Seoul, Korea Selatan) dan Rumah Sakit Boramae Universitas Nasional Seoul (Seoul, Korea Selatan). Semua peserta menjalani wawancara terstruktur oleh seorang psikiater berdasarkan Jadwal Kiddie Korea untuk Gangguan Afektif dan Skizofrenia (K-SADS-PL) (20). Semua peserta menyelesaikan Subjek Desain Blok dan Kosakata Skala Intelegensi Korea-Wechsler untuk Anak-anak, edisi 4th (K-WISC-IV) (21). Impulsif dinilai oleh Barratt Impulsiveness Scale (BIS) (22). Skala Sistem Penghambat Perilaku (BInS) dan Sistem Aktivasi Perilaku (BAS) diukur untuk menilai dimensi kepribadian (23). Kriteria eksklusi meliputi gangguan medis utama masa lalu atau saat ini (misalnya, diabetes mellitus), gangguan neurologis (misalnya, gangguan kejang, cedera kepala), gangguan kejiwaan (misalnya, gangguan depresi mayor, gangguan kecemasan), keterbelakangan mental, atau penyalahgunaan zat apa pun (mis. , tembakau, ganja, alkohol). Karakteristik umum dari subyek penelitian dirangkum dalam Tabel Tabel1.1. Penelitian ini telah disetujui oleh Dewan Peninjauan Institusi dari Universitas Kedokteran Universitas Katolik Korea (MC16SISI0120). Semua peserta dan orang tua mereka memberikan persetujuan tertulis.
Tabel 1
Karakteristik umum dari subyek penelitian.
penemuan | Pengesahan | Gabungan | |||||||
---|---|---|---|---|---|---|---|---|---|
kontrol | IGD | P-nilai | kontrol | IGD | P-nilai | kontrol | IGD | P-nilai | |
N | 26 | 25 | 16 | 20 | 42 | 45 | |||
Umur (tahun) | |||||||||
Median (min – maks) | 13 (12 - 17) | 13 (12 - 15) | 0.759 | 15 (13 - 18) | 14.5 (12 - 18) | 0.628 | 14 (12 - 18) | 14 (12 - 18) | 0.509 |
Jam gaming Internet mingguan (h) | |||||||||
Median (min – maks) | 5.25 (2 - 17) | 18 (6 - 46) | 1.27E − 6a | 5.5 (2 - 23) | 8 (1 - 112) | 0.374 | 5.5 (2 - 23) | 14 (1 - 112) | 1.63E − 5a |
Pendapatan rumah tangga bulanan (juta KRW) | |||||||||
Median (min – maks) | 5 (1 - 9) | 3 (1 - 9) | 0.588 | 4 (4 - 4) | 2 (2 - 2) | 1.000 | 5 (1 - 9) | 3 (1 - 9) | 0.460 |
Pendidikan (tahun) | |||||||||
Median (min – maks) | 8 (7 - 9) | 8 (7 - 9) | 0.584 | 12 (12 - 12) | 6 (6 - 13) | 0.305 | 8 (7 - 12) | 8 (6 - 13) | 0.269 |
K-WISC: desain blok | |||||||||
Median (min – maks) | 10.5 (4 - 17) | 10 (4 - 16) | 0.544 | 10 (3 - 16) | 12.5 (4 - 15) | 0.125 | 10 (3 - 17) | 11 (4 - 16) | 0.598 |
K-WISC: kosakata | |||||||||
Median (min – maks) | 9 (5 - 17) | 7 (5 - 13) | 0.174 | 9.5 (8 - 15) | 11.5 (5 - 15) | 0.595 | 9 (5 - 17) | 9 (5 - 15) | 0.527 |
KS | |||||||||
Median (min – maks) | 24 (17 - 36) | 37 (22 - 51) | 3.81E − 6a | 29 (17 - 34) | 59 (22 - 108) | 1.2E − 5a | 25 (17 - 36) | 40 (22 - 108) | 2.05E − 10a |
BIS | |||||||||
Median (min – maks) | 63 (35 - 75) | 67.5 (45 - 81) | 0.080 | 61 (45 - 79) | 63 (32 - 82) | 0.835 | 62 (35 - 79) | 65 (32 - 82) | 0.240 |
BAS | |||||||||
Median (min – maks) | 31 (15 - 40) | 31 (13 - 51) | 0.558 | 36.5 (22 - 48) | 34 (27 - 52) | 1.000 | 32 (15 - 48) | 34 (13 - 52) | 0.637 |
BIN | |||||||||
Median (min – maks) | 18 (10 - 26) | 17.5 (13 - 27) | 0.642 | 18.5 (12 - 25) | 20 (13 - 21) | 0.138 | 18 (10 - 26) | 19 (13 - 27) | 0.302 |
IGD, pasien gangguan game internet; KS, Skala Kecanduan Internet Korea; BIS, Skala Impulsif Barrat; BAS, Sistem Aktivasi Perilaku; BInS, Sistem Penghambatan Perilaku; KRW, Won Korea.
aP <0.05 (uji Mann – Whitney – Wilcoxon).
Eksperimen TaqMan Low Density miRNA Array (TLDA)
Darah tepi dikumpulkan dari masing-masing peserta dan dipindahkan ke laboratorium dalam 4 jam untuk meminimalkan lisis sel darah. Spesimen disentrifugasi pada 3,000 rpm selama 10 menit pada suhu kamar. Kemudian, supernatan (lapisan plasma) dikumpulkan tanpa mencemari sel-sel darah. MiRNA yang beredar diekstraksi menggunakan TaqMan miRNA ABC Purification Kit (Panel Manusia A; Thermo Fisher Scientific, Waltham, MA, USA) sesuai dengan instruksi pabrik. Singkatnya, 50 µL sampel plasma dan 100 µL buffer ABC dicampur. Setelah hibridisasi dengan manik-manik magnetik anti-miRNA target-spesifik, miRNA yang bersirkulasi terbatas dielusi dari manik-manik dengan 100 µL buffer elution. Pada fase penemuan, 381 miRNA diperiksa dari sampel plasma 51 (25 IGD dan kontrol 26) menggunakan Kit Pemurnian miRNA miqNA TaqMan (Panel Manusia A; Thermo Fisher Scientific, Waltham, MA, USA) sesuai dengan instruksi pabriknya. Transkripsi balik Megaplex dan reaksi pra-amplifikasi dilakukan untuk meningkatkan kuantitas cDNA untuk analisis ekspresi miRNA menggunakan MegaplexPreAmp Primers Human Pool A dan TaqManPreAmp Master Mix (Thermo Fisher Scientific). Panel TLDA A v2.0 (Thermo Fisher Scientific) dijalankan pada sistem PCR waktu-nyata ViiA7 (Thermo Fisher Scientific) untuk mengevaluasi ekspresi miRNAs. Data mentah diproses menggunakan ExpressionSuite Software v1.0.4 (Thermo Fisher Scientific) untuk menentukan nilai Ct untuk setiap miRNA.
Analisis Data untuk TLDA
Kami pertama kali mengukur siklus ambang (nilai Ct) dari setiap miRNA. miRNA dengan nilai Ct> 35 dianggap tidak terdeteksi dan dikeluarkan dari analisis selanjutnya. Semua nilai Ct dinormalisasi ke nilai Ct miR-374b (nilai ΔCt), salah satu miRNA paling stabil yang beredar dalam plasma manusia (24). Rasio lipat-perubahan log2 (nilai ΔΔCt) ekspresi dihitung menggunakan nilai rata-rata sampel kontrol sebagai kalibrator dalam paket HTqPCR dalam Bioconductor (25). Kuantifikasi relatif (RQ) dari setiap target miRNA didefinisikan sebagai 2−ΔΔCt. Untuk pengujian hipotesis dari perbedaan ekspresi antara dua kelompok, kami menerapkan analisis variabel pengganti (SVA) untuk menangkap heterogenitas seperti efek batch dalam percobaan menggunakan sva paket dalam Bioconductor (26). miRNA dengan a P-nilai <0.05 dianggap berbeda secara signifikan antara dua kelompok.
Analisis Pengayaan Set Gen
Untuk analisis pengayaan set gen, kami menggunakan ToppFun di ToppGene Suite (27) untuk menyimpulkan Gen Ontologi (GO) yang diperkaya secara signifikan (GO)28) istilah, jalur, dan istilah penyakit. Sebagai input untuk pendekatan ini, kami menggunakan 1,230 yang diprediksi gen target dari kandidat miRNAs. Analisis jalur digunakan untuk menemukan jalur signifikan dari gen target yang diprediksi menurut KEGG, BioCarta, Reactome, GeneMAPP, dan MSigDB di jalur ToppGene. Signifikansi istilah pengayaan fungsional ditentukan berdasarkan Bonferroni-disesuaikan P-nilai.
Validasi dan Replikasi Reverse Transcription PCR (qRT-PCR)
Untuk memvalidasi miRNA 10 yang secara diferensial diekspresikan dalam tahap penemuan, qRT-PCR dilakukan menggunakan TaqMan MicroRNA Assay (miR-15b-5p, #000390; miR-26b-5p, #000407; miR-29p, #3; miR-000413b; 125; miR-5b-000449p, #200; miR-3c-002300p, #337; miR-5-002156p, #411; miR-5-001610p, #423; mi-5; mi-002340; mi-483; -5p, #002338; dan miR-652-3p, #002352) dan sistem ViiA7 (Life Technologies) sesuai dengan protokol pabrikan. Sepuluh nanogram dari total RNA dikonversi menjadi cDNA untai pertama dengan primer khusus miRNA menggunakan Kit Transkripsi Terbalik MicroRNA TaqMan (#4366596, Life Technologies), diikuti oleh PCR waktu-nyata dengan Probe TaqMan. RQ setiap miRNA didefinisikan sebagai 2−ΔCt, di mana ΔCt adalah perbedaan dalam siklus ambang untuk sampel yang bersangkutan, dinormalisasi terhadap miRNA endogen (miR-374b-5p, #001319). Semua reaksi PCR dilakukan dalam rangkap tiga, dan nilai Ctnya dirata-ratakan. Kami menghitung rasio perubahan lipat log2 (ΔΔCt) dari masing-masing miRNA dengan cara yang sama seperti dalam analisis berbasis array. Tes Mann-Whitney-Wilcoxon non-parametrik dilakukan untuk menguji perbedaan tingkat ekspresi miRNAs dalam dua kelompok dengan ambang P-nilai 0.05.
Analisis Western Blot
Setiap sampel serum pertama kali habis dari protein dengan kelimpahan tinggi 14 teratas (albumin, imunoglobulin G, imunoglobulin A, serotransferrin, haptoglobin, antitripsin alfa-1, fibrinogen, alfa-2 glikoprotein, alfa-1 glikoprotein asam, imunoglobulin Mipro , apolipoprotein A-II, pelengkap C3, dan transthyretin) menggunakan kolom MARS-14 (4.6 × 50 mm, Teknologi Agilent, Santa Clara, CA, USA) sebelum analisis western blot. Fraksi tak terikat yang diperoleh dari kolom MARS-14 dikonsentrasikan menggunakan filter sentrifugal Amicon Ultracel-3 (cutoff 3 kDa), dan kemudian konsentrasi protein ditentukan menggunakan metode asam bicinchoninic. Jumlah yang sama (dari 10 ke 30 μg) dari sampel kontrol dan serum IGD dipisahkan pada 4-20% Mini-PROTEAN TGX gel pracetak (Bio-Rad, CA, USA) dan dipindahkan ke membran difluoride polivinilidena. Selanjutnya, membran diblokir dalam TBS-T (190 mM NaCl, 25 mM Tris-HCl, pH 7.5, dan 0.05% Tween 20) dengan 5% susu kering tanpa lemak pada suhu kamar selama 30 min. Membran kemudian diinkubasi dengan antibodi primer terhadap DPYSL2 (1: 500, Novus Biologicals, Littleton, CO, USA), GABRB2 (1: AbnumX, Abcam, Cambridge, MA, USA), dan CNR1000 (1: Bioteknologi 1 , Inc., Santa Cruz, CA, USA), DUSP100 (4: 1, MybioSource, San Diego, CA, USA), dan PI500 (15: 1, MybioSource, San Diego, CA, USA) di TBS-T dengan 500 % susu kering tanpa lemak pada 5 ° C semalaman, dan kemudian dengan antibodi sekunder yang sesuai, baik bovine anti-mouse (4: 1, Santa Cruz Biotechnology) atau kambing anti-kelinci (1,000: 1, Cell Signaling, Beverly, MA, USA ) dikonjugasikan ke peroksidase lobak pada suhu kamar selama 1,000 h. Deteksi sinyal dilakukan menggunakan chemiluminescence dengan reagen ECL (GE health, Piscataway, NJ, USA). Kami menghitung hasil western blot menggunakan perangkat lunak analisis TotalLab 1D (Non-linear Dynamics, Newcastle upon Tyne, UK). Kemudian, nilai rasio densitometri dihitung dengan membagi nilai densitometri setiap sampel seperti yang dijelaskan di tempat lain (29). Sebagai kontrol untuk normalisasi, sampel serum dikumpulkan dari 46 IGD dan sampel kontrol digunakan untuk setiap percobaan. Signifikansi statistik ditentukan dengan menggunakan uji Mann-Whitney-Wilcoxon non-parametrik dengan ambang batas P-nilai 0.05.
Hasil
Karakteristik Subjek Penelitian
Gambaran demografis dan klinis dari subyek penelitian ditunjukkan pada Tabel Tabel1.1. Ketika kami membandingkan IGD dan kelompok kontrol menurut Skala Kecanduan Internet Korea (K-Skala) seperti yang dijelaskan di tempat lain (20, 30), kelompok IGD menunjukkan nilai median K-Scale yang jauh lebih tinggi daripada kelompok kontrol (37 vs 24, P = 3.81 × 10-6) (Meja (Table1) .1). Waktu rata-rata mingguan yang dihabiskan untuk permainan internet di grup IGD secara signifikan lebih lama daripada kontrol (18 vs 5.25 jam, P = 1.27 × 10-6). Sedangkan tidak ada perbedaan yang signifikan antara dua kelompok dalam usia, pendapatan rumah tangga bulanan, durasi pendidikan, desain blok, dan hasil subtest kosakata dari K-WISC, BIS, BInS, dan BAS.
MiRNA yang Dinyatakan Berbeda Antara IGD dan Kontrol
Untuk menemukan miRNA terkait IGD, kami mengadopsi pendekatan dua langkah (penemuan dan validasi independen). Desain penelitian dan strategi keseluruhan diilustrasikan dalam Gambar S1 dalam Bahan Tambahan. Pada tahap penemuan, kami menganalisis profil ekspresi miRNA dari sampel 51 (25 IGDs dan kontrol 26) menggunakan array miRNA yang mengandung 384 miRNAs. Tingkat ekspresi miRNAs 10 ditemukan berbeda secara signifikan antara IGD dan kelompok kontrol (Tabel (Table2) .2). Level ekspresi relatif dari miRNA 10 ini ditunjukkan pada Gambar Figur1.1. Di antara mereka, dua (hsa-miR-423-5p dan hsa-miR-483-5p) diregulasi dan delapan (hsa-miR-15b-5p, hsa-miR-26b-5p, hsa-miR 29p, hsa-miR 3p hsa-miR-125b-5p, hsa-miR-200c-3p, hsa-miR-337c-5p, hsa-miR-411-5p, dan hsa-miR-652-3p grup di downregulated dalam IGD).
Tabel 2
MikroRNA yang dinyatakan berbeda (miRNA) dan lipat perubahan.
miRNA | penemuan | Pengesahan | Gabungan | |||
---|---|---|---|---|---|---|
P-nilai | Lipat perubahan | P-nilai | Lipat perubahan | P-nilai | Lipat perubahan | |
hsa-miR-15b-5p | 0.033 | 0.829 | 0.694 | 1.119 | 0.381 | 0.947 |
hsa-miR-26b-5pa | 0.008 | 0.871 | 0.049 | 0.841 | 0.013 | 0.857 |
hsa-miR-29b-3p | 0.005 | 0.400 | 0.560 | 1.187 | 0.089 | 0.647 |
hsa-miR-125b-5p | 0.021 | 0.582 | 0.290 | 0.950 | 0.069 | 0.723 |
hsa-miR-200c-3pa | 0.011 | 0.336 | 0.003 | 0.542 | 2.93 × 10-5 | 0.415 |
hsa-miR-337c-5p | 0.009 | 0.385 | 0.582 | 0.872 | 0.020 | 0.553 |
hsa-miR-411-5p | 0.004 | 0.322 | 0.336 | 1.282 | 0.158 | 0.595 |
hsa-miR-423-5p | 0.026 | 1.387 | 0.189 | 0.955 | 0.518 | 1.175 |
hsa-miR-483-5p | 0.018 | 1.861 | 0.765 | 1.413 | 0.211 | 1.647 |
hsa-miR-652-3pa | 0.019 | 0.715 | 0.049 | 0.877 | 0.011 | 0.782 |
amiRNA secara signifikan diubah dalam kedua penemuan dan validasi secara konsisten.
qRT-PCR Validasi miRNA Kandidat
Untuk memvalidasi miRNA kandidat 10, kami melakukan qRT-PCR dengan set validasi independen (20 IGD dan kontrol 16) (Tabel S1 dalam Bahan Tambahan). Tiga dari miRNA ini (hsa-miR-200c-3p, hsa-miR-26b-5p, dan hsa-miR-652-3p) secara signifikan diregulasi ke bawah dalam kelompok IGD dari set validasi (Tabel (Table2) .2). Tiga miRNA lainnya (hsa-miR-337c-5p, hsa-miR-125b, dan hsa-miR-423-5p) juga diregulasi ke bawah dalam kelompok IGD tetapi tidak secara signifikan. Sisa empat miRNAs (hsa-miR-15b-5p, hsa-miR-29b-3p, hsa-miR-411-5p, dan hsa-miR-423-5p) diekspresikan secara berlawanan dalam set validasi. Ketika kami menggabungkan set penemuan dan validasi (total subjek 45 IGD dan kontrol 42), tiga miRNA yang divalidasi secara signifikan signifikan (Tabel (Table2) .2). Informasi terperinci, lokasi kromosom, sekuens matang, dan level ekspresi dalam SSP dari ketiga miRNA ini tersedia di Tabel S2 dalam Bahan Tambahan.
Efek Sinergis dari Perubahan Simultan Tiga MiRNA pada Risiko IGD
Untuk mengevaluasi efek gabungan dari tiga miRNA, kami mengamati odds ratio (OR) dari empat subkelompok (dengan perubahan 0, 1, 2, atau 3 miRNA). perubahan miRNA didefinisikan oleh nilai RQ seperti yang dijelaskan dalam Bagian “Bahan dan Metode. ”Karena ketiga penanda miRNA diregulasi ke bawah dalam kelompok IGD, miRNA yang nilai RQ-nya di bawah satu ditolak seperti yang diubah. Informasi terperinci dari nilai RQ setiap mata pelajaran studi untuk tiga miRNA tersedia di Tabel S3 dalam Bahan Tambahan. Untuk setiap subkelompok, odds dihitung sebagai rasio jumlah kontrol dengan IGD, maka setiap OR dihitung dengan membagi peluang masing-masing subkelompok dengan odds subkelompok tanpa perubahan miRNA. Individu dengan tiga perubahan miRNA menunjukkan risiko 22 kali lebih tinggi daripada mereka yang tanpa perubahan miRNA (OR 22, 95% CI 2.29-211.11). OR menunjukkan tren peningkatan dengan jumlah miRNA yang diubah dari 0 ke 3 (r2 = 0.996) (Gambar (Gambar 22).
GO dan Pathway Analysis of Gen Target dari Kandidat miRNAs
Untuk mendapatkan wawasan tentang fungsi dari tiga penanda miRNA yang secara signifikan diturunkan regulasinya dalam kelompok IGD, gen target mereka diprediksi menggunakan basis data miRWalk 2.0 (31). Total gen 1,230 secara konsisten diprediksi sebagai target hilir oleh empat algoritma (miRWalk, miRanda, RNA22, dan Targetcan) menggunakan database miRWalk (32-34) (Tabel S4 dalam Bahan Tambahan). Analisis pengayaan set gen menggunakan ToppFun di ToppGene Suite menunjukkan bahwa gen target miRNA tersebut secara signifikan terkait dengan jalur pengembangan saraf seperti "panduan akson" dan istilah GO seperti "neurogenesis" (Tabel S5 dalam Bahan Tambahan).
Ekspresi Gen Target yang Diprediksi
Di antara gen target hilir dari tiga miRNA, 140 diprediksi secara bersamaan untuk dua atau lebih miRNA (Tabel S4 dalam Bahan Tambahan). Untuk mengeksplorasi apakah level ekspresi protein dari gen target hilir berbeda antara IGD dan kelompok kontrol, kami memilih gen 2 (DUSP4 dan PI15), yang diprediksi sebagai target hilir semua 3 miRNAs dan gen 3 tambahan (GABRB2, DPYSL2, dan CNR1) dari yang diprediksi untuk miRNA 2 dan melakukan analisis western blot dengan sampel plasma dari 28 IGD dan kontrol 28 yang tersedia untuk percobaan. Kami membandingkan ekspresi dari lima target antara IGD dan kelompok kontrol dengan mengukur intensitas pita dan area seperti yang dijelaskan di tempat lain (29). Di antara mereka, tingkat ekspresi kontrol DPYSL2 (28 IGDs dan 28, P = 0.0037) dan GABBR2 (27 IGD dan 28 kontrol, P = 0.0052) secara signifikan lebih tinggi pada kelompok IGD (Gambar (Figure3) .3). Namun, kami tidak dapat mengamati ekspresi diferensial dari CNR1 (P = 0.0853), DUSP4 (P = 0.5443), dan PI15 (P = 0.6346).
Diskusi
Telah dilaporkan bahwa miRNA terlibat dalam perkembangan saraf (35, 36), dan ekspresi diferensial miRNA otak diamati pada penyakit kejiwaan seperti skizofrenia (37). Oleh karena itu, masuk akal bahwa profil miRNA yang bersirkulasi dapat menjadi biomarker yang berguna untuk IGD. MIRNA yang bersirkulasi telah disarankan sebagai biomarker untuk beragam gangguan neuropsikiatri (38-40); Namun, mekanisme molekuler di balik pengembangan IGD sebagian besar masih belum diketahui meskipun penting secara klinis dan sosial. Secara khusus, belum ada penelitian tentang miRNA terkait IGD. Tujuan dari penelitian ini ada dua. Pertama, kami berusaha menemukan miRNAs plasma yang terkait dengan IGD. Kedua, kami mengevaluasi implikasi biologis dari kandidat miRNA dengan mengeksplorasi ekspresi protein dan GO dari gen target hilir. Melalui penapisan genom profil ekspresi miRNA dan validasi hilir kandidat, kami menemukan bahwa ekspresi tiga miRNA (hsa-miR-200c-3p, hsa-miR-26b-5p, dan hsa-miR-652-3p) adalah secara signifikan lebih rendah pada pasien IGD daripada kontrol. Meskipun pola ekspresi tujuh kandidat miRNA lainnya tidak direplikasi dalam validasi, itu bisa menjadi negatif palsu karena ukuran sampel yang kecil dalam penelitian ini. Sepengetahuan kami, ini adalah laporan pertama tentang kemungkinan bahwa profil ekspresi miRNA darah bisa menjadi biomarker yang berguna untuk IGD. Kombinasi dari tiga penanda miRNA dapat berfungsi sebagai alat invasif minimal untuk identifikasi awal orang yang berisiko IGD.
MiRNA yang diidentifikasi dalam penelitian ini telah dilaporkan terlibat dalam berbagai gangguan neuropsikiatri. Ekspresi hsa-miR-200c dalam darah telah dilaporkan diturunkan regulasinya pada beberapa gangguan kejiwaan seperti skizofrenia (41) dan episode depresi besar (42). miR-200c dilaporkan lebih tinggi diekspresikan dalam fraksi sinaptik daripada total otak depan (43) dan juga dikaitkan dengan kematian sel saraf (44). Berdasarkan laporan-laporan sebelumnya, miR-200c terlibat dalam pengembangan saraf dan dapat dikaitkan dengan gangguan neuropsikiatri jika ekspresinya terganggu. Beberapa penelitian telah menyarankan hubungan antara miR-652 dan risiko gangguan neuropsikiatri. Mirip dengan pendekatan kami, untuk mengidentifikasi biomarker darah untuk skizofrenia, Lai et al. melakukan analisis TLDA dengan pasien skizofrenia dan kontrol normal, dan menemukan bahwa tujuh miRNA termasuk hsa-miR-652 diekspresikan secara berbeda pada pasien skizofrenia (45). Dalam studi berikutnya, mereka merancang model prediksi menggunakan data ekspresi miRNA dan berhasil membedakan skizofrenia dari kontrol normal (46). Perubahan ekspresi hsa-miR-652 juga diamati pada pecandu alkohol (47). Hsa-miR-26b ditemukan diaktifkan selama diferensiasi sel neuron (48). Perkins et al. melaporkan bahwa hsa-miR-26b diturunkan regulasi di korteks prefrontal pasien skizofrenia (49).
Meskipun tidak ada bukti langsung untuk mendukung hubungan antara ekspresi terganggu miRNA ini dan patofisiologi IGD, kita dapat menyimpulkan bahwa disregulasi miRNA ini dapat dikaitkan dengan patofisiologi IGD berdasarkan berbagai laporan sebelumnya tentang gen hilir yang kami prediksi . Beberapa gen hilir dari tiga miRNA seperti GABRB2, CNR1, NRXN1, dan DPYSL2 dilaporkan dikaitkan dengan gangguan neuropsikiatri. Asam gamma-aminobutyric (GABA) adalah neurotransmitter penghambat utama dalam SSP. Disregulasi reseptor GABA, berimplikasi pada gangguan neuropsikiatri termasuk kecanduan, kecemasan, dan depresi (50), yang juga merupakan fitur utama IGD (8). Polimorfisme genetik pada gen reseptor GABA dilaporkan berhubungan dengan kecanduan alkohol dan skizofrenia (51, 52). Dihydropyrimidinase-like 2 (DPYSL2) adalah anggota dari keluarga mediator protein collapsin response, yang berperan dalam perakitan mikrotubulus, pensinyalan sinaptik, dan regulasi pertumbuhan aksonal. Akibatnya, molekul ini telah disarankan sebagai biomarker untuk gangguan kejiwaan (53, 54). Polimorfisme dalam DPYSL2 gen juga dilaporkan dikaitkan dengan gangguan penggunaan alkohol (55). Laporan sebelumnya dan data kami menunjukkan bahwa ekspresi berlebih dari GABRB2 dan DPYSL2, target hilir miRNA yang diregulasi, memiliki implikasi untuk patogenesis gangguan neuropsikiatri termasuk IGD. Cannabinoid receptor type 1 (CNR1) adalah heteroreseptor presinaptik yang memodulasi pelepasan neurotransmitter dan gangguan dalam pensinyalan cannabinoid yang dikaitkan dengan berbagai gangguan neuropsikiatri (56). Polimorfisme genetik CNR1 gen diketahui terkait dengan ketergantungan zat pada Kaukasia (57). Dalam model tikus, aktivasi ventral hippocampus CNR1 mengganggu perilaku sosial dan kognisi normal (58). Perubahan genetik dalam keluarga NRXN diketahui terlibat dalam beragam gangguan neuropsikiatri termasuk kecanduan (59).
Untuk menguji implikasi biologis dari ketiga kandidat miRNA dengan cara yang lebih langsung, kami mengeksplorasi ekspresi protein dari gen target hilir mereka. Karena terbatasnya ketersediaan sampel plasma, dari kandidat umum 140 (diperkirakan sebagai downstream 2 atau lebih miRNAs), kami memeriksa target 5 (GABRB2, DPYSL2, CNR1, DUSP4, dan PI15) dengan western blot dan mengkonfirmasikan ekspresi dari GAB barat dan DPYSL2 secara signifikan lebih tinggi pada kelompok IGD. Laporan sebelumnya dan data kami menunjukkan bahwa ekspresi berlebih dari GABRB2 dan DPYSL2, target hilir miRNA yang diregulasi, mungkin memiliki implikasi untuk patogenesis gangguan neuropsikiatri termasuk IGD. Hasil GO dan analisis jalur dari jalur perkembangan saraf juga mendukung implikasi neurobiologis dari penanda miRNA. Temuan lain yang menarik adalah efek sinergis dari perubahan simultan miRNA. Individu dengan downregulasi semua miRNA 2 menunjukkan risiko 3 kali lebih tinggi daripada mereka yang tanpa downregulasi, dan OR meningkat secara tergantung dosis. Meskipun CI untuk ketiga perubahan ini luas karena ukuran sampel yang terbatas, korelasi positif yang jelas (r2 = 0.996) mendukung efek sinergis dari tiga miRNA.
Meskipun kami menemukan penanda miRNA terkait IGD dan individu dengan ketiga perubahan miRNA memiliki risiko 22 kali lebih tinggi daripada yang tanpa perubahan miRNA, ada beberapa keterbatasan dalam penelitian ini. Pertama, ukuran sampel kecil meningkatkan kemungkinan hilang penanda miRNA signifikan lainnya. Kedua, karena data kami tidak cukup untuk mengklarifikasi apakah profil miRNA plasma merupakan sebab atau akibat, kami tidak dapat mengkonfirmasi peran biologis dari penanda non-invasif ini dalam pengaturan klinis. Profiling miRNA lebih lanjut dan analisis gen hilirnya menggunakan jaringan otak manusia dari bank jaringan otak dapat memberikan jawaban yang lebih langsung. Analisis jaringan otak dengan model hewan gangguan permainan juga akan sangat membantu. Ketiga, karena terbatasnya ketersediaan sampel plasma, kami memeriksa hanya lima molekul kandidat hilir. Menjelajahi lebih banyak target hilir dengan set sampel yang lebih besar akan sangat membantu untuk lebih memahami mekanisme molekuler IGD.
Singkatnya, melalui penapisan genom profil ekspresi miRNA dan validasi independen, kami menemukan tiga miRNA terkait IGD (hsa-miR-200c-3p, hsa-miR-26b-5p, dan hsa-miR-652-3p). Banyak gen hilir mereka dilaporkan terlibat dalam beragam gangguan neuropsikiatrik, dan validasi eksperimental dari perubahan ekspresi gen-gen hilir ini mendukung implikasi miRNA yang diidentifikasi dalam penelitian ini. Kami menemukan bahwa individu dengan downregulasi ketiga miRNA berisiko tinggi terkena IGD. Bersama-sama dengan faktor-faktor risiko klinis atau lingkungan yang diketahui dan kriteria diagnostik, temuan kami dapat memfasilitasi intervensi dini untuk membantu orang dengan risiko IGD yang lebih tinggi.
Pernyataan etika
Penelitian ini telah disetujui oleh Dewan Peninjauan Institusi dari Universitas Kedokteran Universitas Katolik Korea (MC16SISI0120). Semua peserta dan orang tua mereka memberikan persetujuan tertulis.
Kontribusi Penulis
ML dan HC berkontribusi secara merata pada makalah ini. ML, D-JK, dan Y-JC merancang penelitian ini. SJ, S-MC, YP, DC, dan JL melakukan percobaan dan pembuatan data. J-WC, S-HP, J-SC, dan D-JK mengumpulkan sampel darah dan informasi klinis. ML, HC, S-HY, dan Y-JC menganalisis data. ML, HC, S-HY, dan Y-JC menggambarkan naskah itu. Y-JC mengawasi proyek.
Pernyataan Benturan Kepentingan
Para penulis menyatakan bahwa penelitian ini dilakukan tanpa adanya hubungan komersial atau keuangan yang dapat ditafsirkan sebagai potensi konflik kepentingan.
Catatan kaki
Pendanaan. Pekerjaan ini didukung oleh hibah dari Program Penelitian Otak melalui National Research Foundation of Korea (NRF), yang didanai oleh Kementerian Sains dan TIK dan Perencanaan Masa Depan (NRF-2015M3C7A1064778).
Materi tambahan
Materi Tambahan untuk artikel ini dapat ditemukan online di http://www.frontiersin.org/articles/10.3389/fpsyt.2018.00081/full#supplementary-material.
Referensi