Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A non-coding CRHR2 SNP rs255105, a cis-eQTL for a downstream lincRNA AC005154.6, is associated with heroin addiction

  • Orna Levran ,

    Roles Conceptualization, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing

    levrano@rockefeller.edu

    Affiliation The Laboratory of the Biology of Addictive Diseases, The Rockefeller University, New York, New York, United States of America

  • Joel Correa da Rosa,

    Roles Formal analysis, Writing – review & editing

    Affiliation Center for Clinical and Translational Science, The Rockefeller University, New York, New York, United States of America

  • Matthew Randesi,

    Roles Investigation, Writing – review & editing

    Affiliation The Laboratory of the Biology of Addictive Diseases, The Rockefeller University, New York, New York, United States of America

  • John Rotrosen,

    Roles Funding acquisition, Resources, Writing – review & editing

    Affiliation NYU School of Medicine, New York, New York, United States of America

  • Miriam Adelson,

    Roles Resources, Writing – review & editing

    Affiliation Dr. Miriam and Sheldon G. Adelson Clinic for Drug Abuse Treatment and Research, Las Vegas, Nevada, United States of America

  • Mary Jeanne Kreek

    Roles Funding acquisition, Supervision, Writing – review & editing

    Affiliation The Laboratory of the Biology of Addictive Diseases, The Rockefeller University, New York, New York, United States of America

Abstract

Dysregulation of the stress response is implicated in drug addiction; therefore, polymorphisms in stress-related genes may be involved in this disease. An analysis was performed to identify associations between variants in 11 stress-related genes, selected a priori, and heroin addiction. Two discovery samples of American subjects of European descent (EA, n = 601) and of African Americans (AA, n = 400) were analyzed separately. Ancestry was verified by principal component analysis. Final sets of 414 (EA) and 562 (AA) variants were analyzed after filtering of 846 high-quality variants. The main result was an association of a non-coding SNP rs255105 in the CRH (CRF) receptor 2 gene (CRHR2), in the discovery EA sample (Pnominal = .00006; OR = 2.1; 95% CI 1.4–3.1). The association signal remained significant after permutation-based multiple testing correction. The result was corroborated by an independent EA case sample (n = 364). Bioinformatics analysis revealed that SNP rs255105 is associated with the expression of a downstream long intergenic non-coding RNA (lincRNA) gene AC005154.6. AC005154.6 is highly expressed in the pituitary but its functions are unknown. LincRNAs have been previously associated with adaptive behavior, PTSD, and alcohol addiction. Further studies are warranted to corroborate the association results and to assess the potential relevance of this lincRNA to addiction and other stress-related disorders.

Introduction

Opioids are highly addictive narcotics used mainly for pain management. Opioid addiction is a significant global public health challenge that is related in part to the unprecedented increase in the use of opioid pain relievers and more recently to the widespread availability of very inexpensive high-quality heroin and fentanyl [1]. Addictions are chronic relapsing diseases characterized by the compulsive use of a drug with destructive outcomes. Both genetic and environmental factors contribute to opioid addiction based on twin and family studies [2, 3]. Effective pharmacotherapies for opioid addiction include agonist maintenance with methadone and partial-agonist maintenance with buprenorphine as well as the antagonist treatment with extended-release injectable naltrexone.

Dysregulation of the stress response is a critical factor in the development of addiction [4]. There is a high inter-individual variability in the response to stress that is determined by an interaction of genetic and non-genetic factors. One of the mechanisms of the stress response is the adrenal secretion of glucocorticoids. Stress, endogenous opioids, and drugs of abuse modulate the hypothalamic-pituitary-adrenal (HPA) axis. Glucocorticoids regulate the activity of the HPA axis through negative feedback. Corticotropin-releasing hormone/factor (CRH/CRF) mediates the HPA axis and also acts in extra-hypothalamic regions to stimulate the mesocorticolimbic dopamine system that mediates the rewarding effects associated with drug use [5].

In our previous studies of stress-related genes that included samples that overlap the samples used in the current study, we have reported associations of variants in FKBP5, GAL, CRHR2, and AVPR1A [69] with heroin addiction. GAL SNPs were associated with opioid addiction by another group [10].

The purpose of this study was to expand previous studies by analyzing a larger number of single nucleotide polymorphisms (SNPs) within specific stress-related genes. Although the samples were genotyped for variants across the whole genome, we have decided a priori to concentrate on limited regions of a relevant function related to addiction thereby reducing the multiple testing burden. The number of subjects used in this study is relatively small but they were recruited with comprehensive ascertainment and stringent inclusion criteria. In addition, the ‘case’ subjects were preferentially selected from the extreme range of the addiction phenotype to increase the power of the study.

Materials and methods

The present study focuses on two discovery samples that are part of a multi-ancestry sample (n = 1811) that included subjects with multiple addictions. This sample is shared by the Laboratory of the Biology of Addictive Diseases from the Rockefeller University with the National Institute on Drug Abuse (NIDA) Genetics Consortium. The current study was limited to subjects with heroin addiction (‘cases’), as the major addiction. Some of the subjects were also addicted or abusing cocaine and/or alcohol.

The current study includes only subjects with predominantly European (EA) or African ancestry (AA) from the USA (n = 1001). Self-described Hispanics were not included. The EA discovery sample included 601 people that were divided into 459 cases and 142 controls. The AA sample included 400 subjects that were divided into 227 cases and 173 controls.

Subjects were recruited at the Rockefeller University or at specific opiate substitution programs (e.g., Manhattan Campus of VA NY Harbor Health Care System, Weill Medical College of Cornell University, and The Dr. Miriam and Sheldon G. Adelson Clinic for Drug Abuse Treatment and Research, in Las Vegas).

Ascertainment of cases and controls was made by personal interview, performed in a similar manner at the recruiting places, using several instruments: the Addiction Severity Index [11], KMSK [12] and Diagnostic and Statistical Manual of Mental Disorders, 4th Edition (DSM-IV). All cases had a diagnosis of opioid dependence based on lifetime DSM-IV criteria, had a history of at least one year of multiple daily uses of heroin, and were in methadone maintenance treatment in the time of recruitment. The eligibility criterion for the control group was no diagnosis of illicit drug abuse. Subjects with excessive drinking or cannabis use were excluded, as described [6]. Subjects with active DSM-IV axis I disorder were excluded from the study.

The study was approved by the Institutional Review Boards of the Rockefeller University (for Rockefeller University and the Las Vegas clinic) and the VA New York Harbor Healthcare System. All subjects signed informed consent for genetic studies and sharing DNA with NIDA.

A second EA case sample

An independent EA sample (n = 364) was analyzed to corroborate the most significant result obtained in the discovery EA sample. This sample is part of a larger sample obtained from the NIDA Clinical Trials Network study (CTN-0051), a comparative effectiveness trial of 24 weeks of treatment with extended-release naltrexone versus buprenorphine-naloxone [13]. Subjects were recruited at community treatment programs affiliated with the CTN and had DSM-5 opioid use disorder. All sites obtained local Institutional Review Board approval and all participants signed informed consent for genetic studies. This sample was not genotyped with the Smokescreen® array. European ancestry contributions were assessed by Structure analysis in the majority (60%) of the sample. Nineteen samples from individuals who self-reported to have European ancestry were excluded based on Structure estimated European ancestry contribution < 70%. The rest of the sample was included based on self-report.

Genotyping

DNA was extracted from blood and quantified using standard methods. The original multi-ancestry sample was genotyped with the Smokescreen® array [14] at RUCDR Infinite Biologics at The Rutgers University, as part of the NIDA collaborative project of opioid addiction. Smokescreen® is a genome-wide custom genotyping array of biallelic SNPs and simple indels with addiction-related gene content. CEL files of the current study samples were analyzed with Axiom Analysis Suite 2.0.0.3.5 (Affymetrix, Santa Clara, CA).

The second replication EA cohort was genotyped using Taqman pre-designed assay C_2267604_10, according to the manufacturer’s protocol (Thermo Fisher Scientific, Waltham, MA, USA). The majority of the sample (60%) was also genotyped for ancestry informative markers (AIMs) with Illumina GoldenGate Custom Panel, as described [15].

Duplicates and relatives

Familial relationships and duplications were assessed in the original cohort via pairwise Identity by Descent (IBD) analyses in PLINK. The filtered genome-wide SNP set was used to calculate PI_HAT values, which represent the proportion of genome IBD sharing [16]. All known familial relationships were consistent with estimated genome sharing. Duplicates (PI_HAT > 0.99) and relatives (PI_HAT > 0.2) were excluded.

Principal component analysis (PCA)

PCA was carried out on the original multi-ancestry cohort with the filtered genome-wide SNP set to investigate population structure, using the R package SNPRelate [17]. Seven PCs explained most of the variability of this set and were used for exclusion of outliers and adjustment of the groups that were originally based on self-reported ancestry. Three PCs were included as covariates in the analyses.

Structure analysis

Structure analysis of 155 AIMs was performed for the second replication EA sample. Each subject was anchored against 1051 samples from 51 worldwide populations, as described [18]. This sample was not genotyped with the Smokescreen® array. The markers were genotyped using an Illumina GoldenGate Custom Panel, as described [15].

Stress-related SNP set

The current analysis was limited a priori to 11 stress-related genes (Table 1). SNPolisher and Ps_Classification functions (Axiom Analysis Suite) were used to select the best performing probes based on quality control metrics (e.g., call rate, hemizygosity, frequency, and cluster separation). A total of 846 high-quality informative variants (included 30 simple indels) were selected for analysis from these genes (S1 Table). Cluster plots were visually evaluated for the SNPs with the most significant results.

Statistical analysis

Pairwise linkage disequilibrium (LD) (D’ and r2) was estimated using Haploview 4.2. LD blocks were identified using the D’ confidence interval bound of 0.7–0.98 [19]. The single-SNP association analyses were conducted using PLINK 1.9 [20] by logistic regression, under dominant or recessive model assumptions. The following filters were used: a. Exact tests for deviation from Hardy-Weinberg equilibrium (HWE) with a threshold of P = 0.005; b. MAF < 0.05; c. missing genotype data (< 90%). A maximum test statistic was applied to account for the dominant and the recessive model tests, using Sumstat [21]. Correction for multiple testing was performed by permutation test (n = 100,000), for the dominant model that showed nominally significant results, using PLINK. Sex and the PCs 1–3 were included as covariates.

The SNPs set used for association analysis was further pruned based on pairwise LD estimates using PLINK (window size of 300 SNPs, advanced by 25 SNPs and an r2 threshold of 0.8) separately for EA and AA. The EA LD-pruned set of 214 SNPs was used for permutation analysis in EA, and the AA LD-pruned set of 445 SNPs was used for permutation analysis in AA.

Analysis of the second EA case sample for the most significant SNP identified in the analysis of the discovery EA sample was conducted by logistic regression under the dominant model either separately or by adding the replication sample to the discovery EA case sample.

Results

Principal component analysis (PCA) was carried out on the original multi-ancestry cohort (n = 1811) using filtered genome-wide data from a custom addiction-related array Smokescreen® [14] (Fig 1). The first PC distinguished between African and non-African ancestry. The second PC distinguished between European and Asian ancestry as well as between European and Native American ancestry. Hispanic and mixed samples showed intermediate values.

thumbnail
Fig 1. Scatter plot of the two main principal components based on filtered genome-wide genotype data of the original cohort.

Each dot represents one individual colored by self-described ancestry. Color code: red: Africans; yellow: Asians; blue: Europeans; green: Hispanics; black: Native Americans; and gray: ‘others’. The corners correspond to high proportions of European, African, Asian and Native American ancestry.

https://doi.org/10.1371/journal.pone.0199951.g001

The two discovery samples (European ancestry and African ancestry) were defined based on data from seven PCs and analyzed separately. Outliers were excluded and misclassified samples were corrected. The three main PCs were included as continuous covariates in the analyses. Family members (PI_HAT > 0.2) were excluded based on the proportion of genome IBD sharing.

The current analysis was limited a priori to 11 stress-related genes (Table 1). A total of 846 high-quality informative variants (included 30 simple indels) were selected from this array (S1 Table).

The European ancestry (EA) sample

The EA sample was divided into 459 cases (former heroin addicts in a methadone maintenance treatment program) and 142 controls. From the original set of 846 SNPs, 389 SNPs were excluded based on low frequency (minor allele frequency, MAF < 0.05), 40 additional variants were removed due to missing genotype data (< 90%), and three variants were removed due to Hardy-Weinberg disequilibrium (HWE). A final set of 414 variants (including eight simple indels) was used for the association analysis under two models of inheritance (dominant and recessive). Linkage disequilibrium (LD) analysis revealed that 110 of these SNPs can be tagged by 38 tag SNP (r2 = 1) (S1 Table). LD-pruned set (r2 < 0.8) of 200 SNPs was used for permutation analysis.

Comparison of genotype frequency distributions between cases and controls revealed significant differences in SNPs in three genes: CRHR2, CRHR1, and FKBP5. Table 2 lists the SNPs that passed the nominal threshold of P < 0.01 with sex and the three main PCs as covariates. The strongest association that remained significant after permutation analysis was detected for the non-coding CRHR2 SNP rs255105 under the dominant model (Pcorrected = 0.01).

thumbnail
Table 2. Top associations in the discovery EA sample (Pnom < 0.01).

https://doi.org/10.1371/journal.pone.0199951.t002

SNP rs255105 is highly conserved (Fig 2). The frequency of the ancestral T allele is ~30% in HapMap populations of European ancestry (26% in the EA control sample) and 90% in HapMap populations of African ancestry (83% in the AA control sample). The frequency of the carriers of the T allele (CT and TT) was 62% in the case sample compared to 43% in the control sample suggesting that the T allele is associated with increased risk for opioid addiction in this population (OR = 2.1).

thumbnail
Fig 2. Evolutionary conservation of SNP rs255015 region.

The ancestral T allele is highly conserved across representative species.

https://doi.org/10.1371/journal.pone.0199951.g002

Independent EA case sample

The most significant result obtained in discovery EA sample (CRHR2 SNP rs255105) was genotyped in an independent EA case sample (n = 364). The difference in genotype frequency distributions between the original controls sample (n = 142) and the second EA case sample remained significant in the same direction (P = 0.002). Analysis of the combined EA case samples (n = 823) was also significant (P = 0.0001; OR = 1.65; 95% CI 1.2–2.2).

The African ancestry (AA) sample

The AA sample included 400 subjects that were divided into 227 cases and 173 controls. A total of 40 variants were removed due to missing genotype data (< 90%). A total of 239 SNPs were excluded based on low frequency (MAF < 0.05). Five variants were removed due to Hardy-Weinberg exact test. A final set of 562 variants was used for the association analysis under two models of inheritance (dominant and recessive). LD analysis revealed that 36 of these SNPs are in complete LD with another SNP in this set (S1 Table). LD-pruned set (r2 < 0.8) of 445 SNPs was used for permutation analysis.

Comparison of genotype frequency distributions between cases and controls revealed several SNPs with nominally significant differences but none of them survived permutation analysis. No association was detected for CRHR2 SNP rs255105 in this sample. The T allele, the minor allele in the EA sample, is the major allele (83%) in the AA sample.

Bioinformatics analysis

The CRHR2 (CRF2) gene (Chr 7:30,651,943–30,700,129, hg38) has several transcription variants that extend the 5' end of the gene [22]. SNP rs255105 is located upstream of the main transcript 1 (alpha, NM_001883.4) and in an intron of transcript variant 2 (beta, NM_001202475.1) and variant 3 (gamma, NM_001202481.1) (Fig 3). SNP rs255105 is also located upstream of the noncoding RNA LOC105375220 that overlaps transcript variants 2 and 3 in a reverse orientation (Fig 3).

thumbnail
Fig 3. Genomic context of SNP CRHR2 SNP rs255015.

The genomic region spanning the AC005154.6 and CRHR2 genes. The main CRHR2 transcripts (V1-V3) are shown.

https://doi.org/10.1371/journal.pone.0199951.g003

Data from the Genotype-Tissue Expression (GTEx) project release v6p [23] indicate that SNP rs255105 is a cis-expression quantitative trait locus (cis-eQTL) for a long intergenic non-coding RNA (lincRNA) gene AC005154.6 (Chr 7: 30,516,309–30,594,809; hg38) (Fig 3). The ancestral T allele of rs255105 is associated with higher expression of this gene (P = 0.000008 in the pituitary) (Fig 4). Notably, there is inconsistency in the annotations of this lincRNA in different sources (AC005154.1, ENSG00000196295; LOC401320).

thumbnail
Fig 4. GTEx box plot showing relationships between SNP rs255015 (7_30732107_T_C_b37) genotype and lincRNA AC005154.6 (ENSG00000196295.7) expression in the pituitary (P = 0.000008).

The sample groups of the different genotypes are indicated on the X-axis; and the relative expression level of lincRNA AC005154.6 is shown on the Y-axis. The T allele (the reference allele and the ancestral minor allele in the EA population) leads to higher expression of lincRNA AC005154.6. Homo Ref: homozygote for the reference T allele, Het: heterozygote; Homo Alt: homozygote for the alternative allele.

https://doi.org/10.1371/journal.pone.0199951.g004

AC005154.6 is located downstream of the CRHR2 gene between the genes GGCT (gamma-glutamylcyclotransferase) and GARS (glycyl-tRNA synthetase) (Fig 3). It has numerous alternative splice variants and its first exon overlaps the GARS gene first exon in the reverse orientation. The region overlapping CRHR2, GARS, AC005154.6 and GGCT displays strong LD and several SNPs in this region including in CRHR2 are associated with AC005154.6 expression, but they are not in strong LD with rs255105.

Only a part of AC005154.6 is evolutionally conserved among primates, and only small regions are conserved in rodents (Fig 5). There are several regions enriched with histone modifications (Fig 5). RNA-Seq data (GTEx) indicate expression of AC005154.6 in various tissues with relatively higher expression in the pituitary.

thumbnail
Fig 5. Multiple sequence alignment and histone marks of the genomic region of AC005154.6 (USCS genome browser GRCh38/hg38; http://genome.ucsc.edu/).

https://doi.org/10.1371/journal.pone.0199951.g005

Discussion

The strongest association with opioid addiction in the current study was detected for the non-coding CRHR2 SNP rs255105 in subjects with European ancestry. The result was corroborated when an independent EA sample was added to the analysis. The current study supplements our previous studies of stress-related genes in EA and AA [8, 9]. The previous EA study included a larger sample that included the majority of the current EA sample but was less homogenous and was limited in the number of SNPs analyzed per each gene. The main associations reported in the previous study were of FKBP5 SNPs. The current study identified the same FKBP5 associations although they did not reach significance after correction for multiple testing. Notably, the previous study also detected nominally significant association of CRHR2 SNP rs255102 but did not detect an association with rs255105. In the current study, the association of SNP rs255102 did not reach nominal significance (P = 0.05).

SNP rs255105 is located upstream of the main CRHR2 transcript and in an intron of several longer CRHR2 transcript variants. Interestingly, it is significantly associated with expression of a downstream lincRNA gene AC005154.6 in several tissues. The ancestral T variant of SNP rs255105, associated with the risk of developing an opioid addiction, is highly conserved. The variant C allele is relatively rare (10%) in populations of African ancestry but more frequent (~30%) in populations of European ancestry. This change in allele frequency could indicate a positive selection toward the C allele that is suggested in this study to have a protective effect on the development of opioid addiction.

CRH receptor 2

The CRH (CRF) family of peptides includes CRH and urocortin (UCN) 1–3. CRH is a neuropeptide released by the paraventricular nucleus of the hypothalamus following exposure to a stressor. It triggers the release of adrenocorticotropic hormone (ACTH) from the anterior pituitary, which in turn stimulates secretion of glucocorticoids from the adrenal gland [24]. Glucocorticoids activate the stress response and exert negative feedback on the HPA axis. CRH is also expressed as a neurotransmitter in extra-hypothalamic regions of the CNS.

CRH has two receptors: CRHR1 and CRHR2 that share ~70% of their amino acid sequences, but differ in their localization and binding affinities. They modulate anxiety-like behavior in a brain region- and cell-type dependent manner [25]. There is limited data on CRFR2 signaling but its effects seem to be plastic and dependent on past experience and ligand quantity. CRHR2 is highly expressed in the pituitary and is densely populated in brain structures involved in anxiety, fear, and arousal, such as the bed nucleus of the stria terminalis (BNST), a subregion of the extended amygdala that regulates the HPA axis. Rodent studies revealed an important role for CRHR2 in PTSD-like behavior [26], and suggest that optimal levels of CRHR2 are critical for coping with a traumatic event. Other rodent studies suggested that selective CRHR2 receptor agonists could be used as a therapy for nicotine addiction [27].

Several CRHR2 SNPs have been reported to be associated with PTSD in EA women [28], and with major depressive disorder (MDD) in Japanese [29] but none of them is in LD with SNP rs255105 indicated in the current study. It remains to be determined if the identified SNP directly affects CRHR2 expression or splicing.

Long intergenic non-coding RNAs

The association of SNP rs255105 with the expression of a downstream lincRNA AC005154.6 is intriguing and may have clinical significance. Notably, lincRNA AC005154.6 is a cis-eQTL of other SNPs in the region including SNPs that are not in strong LD with SNP rs255105.

Our understanding of the complexity of the transcriptome has transformed as a result of improved RNA sequencing technologies. Non-coding RNAs (ncRNAs) constitute the majority of the transcriptome and are important transcriptional and posttranslational regulators [30, 31]. LncRNAs, a subclass of ncRNAs, are non-coding transcripts of >200 bp with exon-intron structure but no open reading frame. LncRNAs are transcribed in complex patterns relative to protein-coding genes and participate in processes like gene expression and epigenetic regulation through both cis- and trans-mediated mechanism. LincRNAs form RNA-protein complexes, modulate chromatin-regulatory proteins or transcription factors and affect the gene expression of nearby genes. They can also interact with miRNAs to modulate their expression and activity [32].

LncRNAs that are encoded between genes, like AC005154.6, are known as long intergenic non-coding RNAs (lincRNAs). LncRNAs play a role in a wide range of functions including adaptive behavior and neurological diseases and are potential new biomarkers and drug targets for stress-related disorders [33]. The majority of lncRNAs show tissue- and temporal-specific expression in the CNS. LncRNAs have been reported to be associated with more than 200 diseases including cancer, PTSD [34] depression [35], and schizophrenia [36].

The contribution of lncRNAs to the development of drug addiction is currently largely unknown and lncRNAs have not yet been characterized in animal models of addiction [37, 38]. Several lncRNAs were shown to be upregulated in the nucleus accumbens from human postmortem brain of heroin abusers compared to drug-free controls, using Affymetrix arrays [39]. Studies have identified lincRNAs that are involved in alcohol dependence [40]. Association studies revealed that alcohol dependence was associated with lncRNA LOC100507053 [41] and lncRNA LOC339975 [42].

The finding that only part of the AC005154.6 region is evolutionally conserved among primates and only small regions are conserved in rodents is important in the evaluation of its functionality. It also implies that the effects of this gene most likely cannot be studied in rodents and demonstrates that mice models may not be reliable as preclinical models for phenotypes caused by newly evolved lincRNAs.

Conclusion

Hypothesis-driven study of variants in stress-related genes and heroin addiction identified an association with non-coding SNP in the CRHR2 gene. A downstream lincRNA gene is a cis-eQTL of this SNP and additional SNPs in the region. This result is intriguing since lincRNAs play a role in adaptive behavior and neurological diseases and are potential drug targets for stress-related disorders. Further analyses are needed to corroborate the association result, to elucidate the potential biological mechanism underlying this lincRNA, and to assess the relevance of this lincRNA to heroin addiction and other stress-related disorders.

Supporting information

Acknowledgments

We are grateful to NIDA and the NIDA genetic consortium for supporting and coordinating the Smokescreen® array genotyping at RUCDR Infinite Biologics at the Rutgers University.

References

  1. 1. Kolodny A, Courtwright DT, Hwang CS, Kreiner P, Eadie JL, Clark TW, et al. The prescription opioid and heroin crisis: a public health approach to an epidemic of addiction. Annu Rev Public Health. 2015;36: 559–74. pmid:25581144
  2. 2. Tsuang MT, Lyons MJ, Eisen SA, Goldberg J, True W, Lin N, et al. Genetic influences on DSM-III-R drug abuse and dependence: a study of 3,372 twin pairs. Am J Med Genet. 1996;67: 473–7. pmid:8886164
  3. 3. Kendler KS, Karkowski LM, Neale MC, Prescott CA. Illicit psychoactive substance use, heavy use, abuse, and dependence in a US population-based sample of male twins. Arch Gen Psychiatry. 2000;57: 261–9. pmid:10711912
  4. 4. Koob G, Kreek MJ. Stress, dysregulation of drug reward pathways, and the transition to drug dependence. Am J Psychiatry. 2007;164: 1149–59. pmid:17671276
  5. 5. Koob GF. The role of CRF and CRF-related peptides in the dark side of addiction. Brain Res. 2010;1314: 3–14. pmid:19912996
  6. 6. Levran O, Londono D, O’Hara K, Nielsen DA, Peles E, Rotrosen J, et al. Genetic susceptibility to heroin addiction: a candidate gene association study. Genes Brain Behav. 2008;7: 720–9. pmid:18518925
  7. 7. Levran O, Londono D, O’Hara K, Randesi M, Rotrosen J, Casadonte P, et al. Heroin addiction in African Americans: a hypothesis-driven association study. Genes Brain Behav. 2009;8: 531–40. pmid:19500151
  8. 8. Levran O, Randesi M, Li Y, Rotrosen J, Ott J, Adelson M, et al. Drug addiction and stress-response genetic variability: association study in African Americans. Ann Hum Genet. 2014;78: 290–8. pmid:24766650
  9. 9. Levran O, Peles E, Randesi M, Li Y, Rotrosen J, Ott J, et al. Stress-related genes and heroin addiction: a role for a functional FKBP5 haplotype. Psychoneuroendocrinology. 2014;45: 67–76. pmid:24845178
  10. 10. Beer B, Erb R, Pavlic M, Ulmer H, Giacomuzzi S, Riemer Y, et al. Association of Polymorphisms in Pharmacogenetic Candidate Genes (OPRD1, GAL, ABCB1, OPRM1) with Opioid Dependence in European Population: A Case-Control Study. PLoS One. 2013;8: e75359. pmid:24086514
  11. 11. McLellan AT, Kushner H, Metzger D, Peters R, Smith I, Grissom G, et al. The Fifth Edition of the Addiction Severity Index. J Subst Abuse Treat. 1992;9: 199–213. pmid:1334156
  12. 12. Kellogg SH, McHugh PF, Bell K, Schluger JH, Schluger RP, LaForge KS, et al. The Kreek-McHugh-Schluger-Kellogg scale: a new, rapid method for quantifying substance abuse and its possible applications. Drug Alcohol Depend. 2003;69: 137–50. pmid:12609695
  13. 13. Lee JD, Nunes EV Jr., Novo P, Bachrach K, Bailey GL, Bhatt S, et al. Comparative effectiveness of extended-release naltrexone versus buprenorphine-naloxone for opioid relapse prevention (X:BOT): a multicentre, open-label, randomised controlled trial. Lancet. 2017. pmid:29150198
  14. 14. Baurley JW, Edlund CK, Pardamean CI, Conti DV, Bergen AW. Smokescreen: a targeted genotyping array for addiction research. BMC Genomics. 2016;17: 145. pmid:26921259
  15. 15. Levran O, Randesi M, da Rosa JC, Ott J, Rotrosen J, Adelson M, et al. Overlapping dopaminergic pathway genetic susceptibility to heroin and cocaine addictions in African Americans. Ann Hum Genet. 2015;79: 188–98. pmid:25875614
  16. 16. Turner S, Armstrong LL, Bradford Y, Carlson CS, Crawford DC, Crenshaw AT, et al. Quality control procedures for genome-wide association studies. Curr Protoc Hum Genet. 2011;Chapter 1: Unit1.19. pmid:21234875
  17. 17. Zheng X, Levine D, Shen J, Gogarten SM, Laurie C, Weir BS. A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics. 2012;28: 3326–8. pmid:23060615
  18. 18. Ducci F, Roy A, Shen PH, Yuan Q, Yuan NP, Hodgkinson CA, et al. Association of substance use disorders with childhood trauma but not African genetic heritage in an African American cohort. Am J Psychiatry. 2009;166: 1031–40. pmid:19605534
  19. 19. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The structure of haplotype blocks in the human genome. Science. 2002;296: 2225–9. pmid:12029063
  20. 20. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81: 559–75. pmid:17701901
  21. 21. Hoh J, Wille A, Ott J. Trimming, weighting, and grouping SNPs in human case-control association studies. Genome Res. 2001;11: 2115–9. pmid:11731502
  22. 22. Catalano RD, Kyriakou T, Chen J, Easton A, Hillhouse EW. Regulation of corticotropin-releasing hormone type 2 receptors by multiple promoters and alternative splicing: identification of multiple splice variants. Mol Endocrinol. 2003;17: 395–410. pmid:12554761
  23. 23. GTEx Consortium. Genetic effects on gene expression across human tissues. Nature. 2017;550: 204–13. pmid:29022597
  24. 24. Henckens MJ, Deussing JM, Chen A. Region-specific roles of the corticotropin-releasing factor-urocortin system in stress. Nat Rev Neurosci. 2016;17: 636–51. pmid:27586075
  25. 25. Chen A. Genetic Dissection of the Neuroendocrine and Behavioral Responses to Stressful Challenges. In: Pfaff D, Christen Y, editors. Stem Cells in Neuroendocrinology. Cham (CH): Springer; 2016. p. 69–79.
  26. 26. Lebow M, Neufeld-Cohen A, Kuperman Y, Tsoory M, Gil S, Chen A. Susceptibility to PTSD-like behavior is mediated by corticotropin-releasing factor receptor type 2 levels in the bed nucleus of the stria terminalis. J Neurosci. 2012;32: 6906–16. pmid:22593059
  27. 27. Bagosi Z, Palotai M, Simon B, Bokor P, Buzas A, Balango B, et al. Selective CRF2 receptor agonists ameliorate the anxiety- and depression-like state developed during chronic nicotine treatment and consequent acute withdrawal in mice. Brain Res. 2016;1652: 21–9. pmid:27693397
  28. 28. Wolf EJ, Mitchell KS, Logue MW, Baldwin CT, Reardon AF, Humphries DE, et al. Corticotropin releasing hormone receptor 2 (CRHR-2) gene is associated with decreased risk and severity of posttraumatic stress disorder in women. Depress Anxiety. 2013;30: 1161–9. pmid:24123648
  29. 29. Ishitobi Y, Nakayama S, Yamaguchi K, Kanehisa M, Higuma H, Maruyama Y, et al. Association of CRHR1 and CRHR2 with major depressive disorder and panic disorder in a Japanese population. Am J Med Genet B Neuropsychiatr Genet. 2012;159B: 429–36. pmid:22467522
  30. 30. Qureshi IA, Mattick JS, Mehler MF. Long non-coding RNAs in nervous system function and disease. Brain Res. 2010;1338: 20–35. pmid:20380817
  31. 31. DiStefano JK. The Emerging Role of Long Noncoding RNAs in Human Disease. Methods Mol Biol. 2018;1706: 91–110. pmid:29423795
  32. 32. Kopp F, Mendell JT. Functional Classification and Experimental Dissection of Long Noncoding RNAs. Cell. 2018;172: 393–407. pmid:29373828
  33. 33. Spadaro PA, Flavell CR, Widagdo J, Ratnu VS, Troup M, Ragan C, et al. Long Noncoding RNA-Directed Epigenetic Regulation of Gene Expression Is Associated With Anxiety-like Behavior in Mice. Biol Psychiatry. 2015;78: 848–59. pmid:25792222
  34. 34. Guffanti G, Galea S, Yan L, Roberts AL, Solovieff N, Aiello AE, et al. Genome-wide association study implicates a novel RNA gene, the lincRNA AC068718.1, as a risk factor for post-traumatic stress disorder in women. Psychoneuroendocrinology. 2013;38: 3029–38. pmid:24080187
  35. 35. Huang X, Luo YL, Mao YS, Ji JL. The link between long noncoding RNAs and depression. Prog Neuropsychopharmacol Biol Psychiatry. 2017;73: 73–8. pmid:27318257
  36. 36. Chen S, Sun X, Niu W, Kong L, He M, Li W, et al. Aberrant Expression of Long Non-Coding RNAs in Schizophrenia Patients. Med Sci Monit. 2016;22: 3340–51. pmid:27650396
  37. 37. Walker DM, Nestler EJ. Neuroepigenetics and addiction. Handb Clin Neurol. 2018;148: 747–65. pmid:29478612
  38. 38. Sartor GC, St Laurent G 3rd, Wahlestedt C. The Emerging Role of Non-Coding RNAs in Drug Addiction. Front Genet. 2012;3: 106. pmid:22737160
  39. 39. Michelhaugh SK, Lipovich L, Blythe J, Jia H, Kapatos G, Bannon MJ. Mining Affymetrix microarray data for long non-coding RNAs: altered expression in the nucleus accumbens of heroin abusers. J Neurochem. 2011;116: 459–66. pmid:21128942
  40. 40. Mayfield RD. Emerging roles for ncRNAs in alcohol use disorders. Alcohol. 2017;60: 31–9. pmid:28438526
  41. 41. Gelernter J, Kranzler HR, Sherva R, Almasy L, Koesterer R, Smith AH, et al. Genome-wide association study of alcohol dependence:significant findings in African- and European-Americans including novel risk loci. Mol Psychiatry. 2014;19: 41–9. pmid:24166409
  42. 42. Adkins AE, Hack LM, Bigdeli TB, Williamson VS, McMichael GO, Mamdani M, et al. Genomewide Association Study of Alcohol Dependence Identifies Risk Loci Altering Ethanol-Response Behaviors in Model Organisms. Alcohol Clin Exp Res. 2017;41: 911–28. pmid:28226201