Protocol for analysis of RNA-sequencing and proteome profiling data for subgroup identification and comparison.

RNA-sequencing and quantitative proteomic profiling simultaneously measure thousands of molecules and provide opportunities to decipher the transcriptomic and proteomic landscapes of cohort specimens for basic and health research. We present a protocol for the analysis of paired transcriptome and proteome data to identify and compare molecular subgroups among cohort specimens. We demonstrate a streamlined analysis workflow, applicable for both transcriptome and proteome data, which allows the comparison of two data types for RNA-protein variations and for derivation of biological implications. For complete details on the use and execution of this protocol, please refer to Yang et al. (2021).

Chloroquine treatment induces secretion of autophagy-related proteins and inclusion of Atg8-family proteins in distinct extracellular vesicle populations.

Chloroquine (CQ), a lysosomotropic agent, is commonly used to inhibit lysosomal degradation and macroautophagy/autophagy. Here we investigated the cell-extrinsic effects of CQ on secretion. We showed that lysosomal and autophagy inhibition by CQ altered the secretome, and induced the release of Atg8 orthologs and autophagy receptors. Atg8-family proteins, in particular, were secreted inside small extracellular vesicles (sEVs) in a lipidation-dependent manner. CQ treatment enhanced the release of Atg8-family proteins inside sEVs. Using full-length ATG16L1 and an ATG16L1 mutant that enables Atg8-family protein lipidation on double but not on single membranes, we demonstrated that LC3B is released in two distinct sEV populations: one enriched with SDCBP/Syntenin-1, CD63, and endosomal lipidated LC3B, and another that contains LC3B but is not enriched with SDCBP/Syntenin-1 or CD63, and which our data supports as originating from a double-membrane source. Our findings underscore the context-dependency of sEV heterogeneity and composition, and illustrate the integration of autophagy and sEV composition in response to lysosomal inhibition.Abbreviations: ACTB: actin beta; ANOVA: analysis of variance; ATG4B: autophagy related 4B cysteine peptidase; Atg8: autophagy related 8; ATG16L1: autophagy related 16 like 1; ATP5F1A/ATP5a: ATP synthase F1 subunit alpha; CALCOCO2: calcium binding and coiled-coil domain 2; CASP3: caspase 3; CASP7: caspase 7; CQ: chloroquine; CD9: CD9 molecule; CD63: CD63 molecule; DAPI: 4',6-diamidino-2-phenylindole; DQ-BSA: dye quenched-bovine serum albumin; ER: endoplasmic reticulum; ERN1/IRE1a: endoplasmic reticulum to nucleus signaling 1; EV: extracellular vesicles; FBS: fetal bovine serum; FDR: false discovery rate; GABARAP: GABA type A receptor-associated protein; GABARAPL2: GABA type A receptor associated protein like 2; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GFP: green fluorescent protein; GO: gene ontology; HCQ: hydroxychloroquine; HSP90AA1: heat shock protein 90 alpha family class A member 1; IP: immunoprecipitation; KO: knockout; LAMP2: lysosomal associated membrane protein 2; LIR: LC3-interacting region; LMNA: lamin A/C; MAP1LC3B/LC3B: microtubule associated protein 1 light chain 3 beta; MS: mass spectrometry; NBR1: NBR1 autophagy cargo receptor; NCOA4: nuclear receptor coactivator 4; NTA: nanoparticle tracking analysis; PE: phosphatidylethanolamine; PECA: probe-level expression change averaging; SDCBP/syntenin-1: syndecan binding protein; SD: standard deviation; SE: secreted; sEV: small extracellular vesicles; SQSTM1/p62: sequestosome 1; TAX1BP1: Tax1 binding protein 1; TEM: transmission electron microscopy; TMT: tandem-mass tag; TSG101: tumor susceptibility 101; ULK1: unc-51 like autophagy activating kinase 1; WC: whole cell.

Proteotranscriptomic classification and characterization of pancreatic neuroendocrine neoplasms.

Pancreatic neuroendocrine neoplasms (PNENs) are biologically and clinically heterogeneous. Here, we use a multi-omics approach to uncover the molecular factors underlying this heterogeneity. Transcriptomic analysis of 84 PNEN specimens, drawn from two cohorts, is substantiated with proteomic profiling and identifies four subgroups: Proliferative, PDX1-high, Alpha cell-like and Stromal/Mesenchymal. The Proliferative subgroup, consisting of both well- and poorly differentiated specimens, is associated with inferior overall survival probability. The PDX1-high and Alpha cell-like subgroups partially resemble previously described subtypes, and we further uncover distinctive metabolism-related features in the Alpha cell-like subgroup. The Stromal/Mesenchymal subgroup exhibits molecular characteristics of YAP1/WWTR1(TAZ) activation suggestive of Hippo signaling pathway involvement in PNENs. Whole-exome sequencing reveals subgroup-enriched mutational differences, supported by activity inference analysis, and identifies hypermorphic proto-oncogene variants in 14.3% of sequenced PNENs. Our study reveals differences in cellular signaling axes that provide potential directions for PNEN patient stratification and treatment strategies.

Puncta intended: connecting the dots between autophagy and cell stress networks.

Proteome profiling and global protein-interaction approaches have significantly improved our knowledge of the protein interactomes of autophagy and other cellular stress-response pathways. New discoveries regarding protein complexes, interaction partners, interaction domains, and biological roles of players that are part of these pathways are emerging. The fourth Vancouver Autophagy Symposium showcased research that expands our understanding of the protein interaction networks and molecular mechanisms underlying autophagy and other cellular stress responses in the context of distinct stressors. In the keynote presentation, Dr. Wade Harper described his team's recent discovery of a novel reticulophagy receptor for selective autophagic degradation of the endoplasmic reticulum, and discussed molecular mechanisms involved in ribophagy and non-autophagic ribosomal turnover. In other presentations, both omic and targeted approaches were used to reveal molecular players of other cellular stress responses including amyloid body and stress granule formation, anastasis, and extracellular vesicle biogenesis. Additional topics included the roles of autophagy in disease pathogenesis, autophagy regulatory mechanisms, and crosstalk between autophagy and cellular metabolism in anti-tumor immunity. The relationship between autophagy and other cell stress responses remains a relatively unexplored area in the field, with future investigations required to understand how the various processes are coordinated and connected in cells and tissues.Abbreviations: A-bodies: amyloid bodies; ACM: amyloid-converting motif; AMFR/gp78: autocrine motility factor receptor; ATG: autophagy-related; ATG4B: autophagy related 4B cysteine peptidase; CALCOCO2/NDP52: calcium binding and coiled-coil domain 2; CAR T: chimeric antigen receptor T; CASP3: caspase 3; CCPG1: cell cycle progression 1; CAR: chimeric antigen receptor; CML: chronic myeloid leukemia; CCOCs: clear cell ovarian cancers; CVB3: coxsackievirus B3; CRISPR-Cas9: clustered regularly interspaced short palindromic repeats-CRISPR associated protein 9; DDXs: DEAD-box helicases; EIF2S1/EIF-2alpha: eukaryotic translation initiation factor 2 subunit alpha; EIF2AK3: eukaryotic translation initiation factor 2 alpha kinase 3; ER: endoplasmic reticulum; EV: extracellular vesicle; FAO: fatty acid oxidation; GABARAP: GABA type A receptor-associated protein; ILK: integrin linked kinase; ISR: integrated stress response; MTOR: mechanistic target of rapamycin kinase; MPECs: memory precursory effector T cells; MAVS: mitochondrial antiviral signaling protein; NBR1: NBR1 autophagy cargo receptor; PI4KB/PI4KIIIß: phosphatidylinositol 4-kinase beta; PLEKHM1: pleckstrin homology and RUN domain containing M1; RB1CC1: RB1 inducible coiled-coil 1; RTN3: reticulon 3; rIGSRNAs: ribosomal intergenic noncoding RNAs; RPL29: ribosomal protein L29; RPS3: ribosomal protein S3; S. cerevisiae: Saccharomyces cerevisiae; sEV: small extracellular vesicles; S. pombe: Schizosaccharomyces pombe; SQSTM1: sequestosome 1; SF3B1: splicing factor 3b subunit 1; SILAC-MS: stable isotope labeling with amino acids in cell culture-mass spectrometry; SNAP29: synaptosome associated protein 29; TEX264: testis expressed 264, ER-phagy receptor; TNBC: triple-negative breast cancer; ULK1: unc-51 like autophagy activating kinase 1; VAS: Vancouver Autophagy Symposium.

Discovery of a Covalent Inhibitor of KRASG12C (AMG 510) for the Treatment of Solid Tumors.

KRASG12C has emerged as a promising target in the treatment of solid tumors. Covalent inhibitors targeting the mutant cysteine-12 residue have been shown to disrupt signaling by this long-"undruggable" target; however clinically viable inhibitors have yet to be identified. Here, we report efforts to exploit a cryptic pocket (H95/Y96/Q99) we identified in KRASG12C to identify inhibitors suitable for clinical development. Structure-based design efforts leading to the identification of a novel quinazolinone scaffold are described, along with optimization efforts that overcame a configurational stability issue arising from restricted rotation about an axially chiral biaryl bond. Biopharmaceutical optimization of the resulting leads culminated in the identification of AMG 510, a highly potent, selective, and well-tolerated KRASG12C inhibitor currently in phase I clinical trials (NCT03600883).

Genomic characterization of a well-differentiated grade 3 pancreatic neuroendocrine tumor.

Pancreatic neuroendocrine neoplasms (PanNENs) represent a minority of pancreatic neoplasms that exhibit variability in prognosis. Ongoing mutational analyses of PanNENs have found recurrent abnormalities in chromatin remodeling genes (e.g., DAXX and ATRX), and mTOR pathway genes (e.g., TSC2, PTEN PIK3CA, and MEN1), some of which have relevance to patients with related familial syndromes. Most recently, grade 3 PanNENs have been divided into two groups based on differentiation, creating a new group of well-differentiated grade 3 neuroendocrine tumors (PanNETs) that have had a limited whole-genome level characterization to date. In a patient with a metastatic well-differentiated grade 3 PanNET, our study utilized whole-genome sequencing of liver metastases for the comparative analysis and detection of single-nucleotide variants, insertions and deletions, structural variants, and copy-number variants, with their biologic relevance confirmed by RNA sequencing. We found that this tumor most notably exhibited a TSC1-disrupting fusion, showed a novel CHD7-BEND2 fusion, and lacked any somatic variants in ATRX, DAXX, and MEN1.

Diverse mechanisms of autophagy dysregulation and their therapeutic implications: does the shoe fit?

In its third edition, the Vancouver Autophagy Symposium presented a platform for vibrant discussion on the differential roles of macroautophagy/autophagy in disease. This one-day symposium was held at the BC Cancer Research Centre in Vancouver, BC, bringing together experts in cell biology, protein biochemistry and medicinal chemistry across several different disease models and model organisms. The Vancouver Autophagy Symposium featured 2 keynote speakers that are well known for their seminal contributions to autophagy research, Dr. David Rubinsztein (Cambridge Institute for Medical Research) and Dr. Kay F. Macleod (University of Chicago). Key discussions included the context-dependent roles and mechanisms of dysregulation of autophagy in diseases and the corresponding need to consider context-dependent autophagy modulation strategies. Additional highlights included the differential roles of bulk autophagy versus selective autophagy, novel autophagy regulators, and emerging chemical tools to study autophagy inhibition. Interdisciplinary discussions focused on addressing questions such as which stage of disease to target, which type of autophagy to target and which component to target for autophagy modulation. Abbreviations: AD: Alzheimer disease; AMFR/Gp78: autocrine motility factor receptor; CCCP: carbonyl cyanide m-chlorophenylhydrazone; CML: chronic myeloid leukemia; CVB3: coxsackievirus B3; DRPLA: dentatorubral-pallidoluysian atrophy; ER: endoplasmic reticulum; ERAD: ER-associated degradation; FA: focal adhesion; HCQ: hydroxychloroquine; HD: Huntingtin disease; HIF1A/Hif1α: hypoxia inducible factor 1 subunit alpha; HTT: huntingtin; IM: imatinib mesylate; MAP1LC3B: microtubule associated protein 1 light chain 3 beta; NBR1: neighbour of BRCA1; OGA: O-GlcNAcase; PDAC: pancreatic ductal adenocarcinoma; PLEKHM1: pleckstrin homology and RUN domain containing M1; polyQ: poly-glutamine; ROS: reactive oxygen species; RP: retinitis pigmentosa; SNAP29: synaptosome associated protein 29; SPCA3: spinocerebellar ataxia type 3; TNBC: triple-negative breast cancer.

Evolution of tools and methods for monitoring autophagic flux in mammalian cells.

Autophagy is an evolutionarily conserved lysosome-mediated degradation and recycling process, which functions in cellular homeostasis and stress adaptation. The process is highly dynamic and involves autophagosome synthesis, cargo recognition and transport, autophagosome-lysosome fusion, and cargo degradation. The multistep nature of autophagy makes it challenging to quantify, and it is important to consider not only the number of autophagosomes within a cell but also the autophagic degradative activity. The rate at which cargos are recognized, segregated, and degraded through the autophagy pathway is defined as autophagic flux. In practice, methods to measure autophagic flux typically evaluate the lysosome-mediated cargo degradation step by leveraging known autophagy markers such as MAP1LC3B (microtubule-associated proteins 1A/1B light chain 3 beta) or lysosome-dependent fluorescent agents. In this review, we summarize the tools and methods used in mammalian cultured cells pertaining to these two approaches, and highlight innovations that have led to their evolution in recent years. We also discuss the potential limitations of these approaches and recommend using a combination of strategies and multiple different autophagy markers to reliably evaluate autophagic flux in mammalian cells.

Molecular characterization of metastatic pancreatic neuroendocrine tumors (PNETs) using whole-genome and transcriptome sequencing.

Pancreatic neuroendocrine tumors (PNETs) are a genomically and clinically heterogeneous group of pancreatic neoplasms often diagnosed with distant metastases. Recurrent somatic mutations, chromosomal aberrations, and gene expression signatures in PNETs have been described, but the clinical significance of these molecular changes is still poorly understood, and the clinical outcomes of PNET patients remain highly variable. To help identify the molecular factors that contribute to PNET progression and metastasis, and as part of an ongoing clinical trial at the BC Cancer Agency (clinicaltrials.gov ID: NCT02155621), the genomic and transcriptomic profiles of liver metastases from five patients (four PNETs and one neuroendocrine carcinoma) were analyzed. In four of the five cases, we identified biallelic loss of MEN1 and DAXX as well as recurrent regions with loss of heterozygosity. Several novel findings were observed, including focal amplification of MYCN concomitant with loss of APC and TP53 in one sample with wild-type MEN1 and DAXX Transcriptome analyses revealed up-regulation of MYCN target genes in this sample, confirming a MYCN-driven gene expression signature. We also identified a germline NTHL1 fusion event in one sample that resulted in a striking C>T mutation signature profile not previously reported in PNETs. These varying molecular alterations suggest different cellular pathways may contribute to PNET progression, consistent with the heterogeneous clinical nature of this disease. Furthermore, genomic profiles appeared to correlate well with treatment response, lending support to the role of prospective genotyping efforts to guide therapy in PNETs.

Discovery and Structure-Guided Optimization of Diarylmethanesulfonamide Disrupters of Glucokinase-Glucokinase Regulatory Protein (GK-GKRP) Binding: Strategic Use of a N → S (nN → σ*S-X) Interaction for Conformational Constraint.

The HTS-based discovery and structure-guided optimization of a novel series of GKRP-selective GK-GKRP disrupters are revealed. Diarylmethanesulfonamide hit 6 (hGK-hGKRP IC50 = 1.2 µM) was optimized to lead compound 32 (AMG-0696; hGK-hGKRP IC50 = 0.0038 µM). A stabilizing interaction between a nitrogen atom lone pair and an aromatic sulfur system (nN → σ*S-X) in 32 was exploited to conformationally constrain a biaryl linkage and allow contact with key residues in GKRP. Lead compound 32 was shown to induce GK translocation from the nucleus to the cytoplasm in rats (IHC score = 0; 10 mg/kg po, 6 h) and blood glucose reduction in mice (POC = -45%; 100 mg/kg po, 3 h). X-ray analyses of 32 and several precursors bound to GKRP were also obtained. This novel disrupter of GK-GKRP binding enables further exploration of GKRP as a potential therapeutic target for type II diabetes and highlights the value of exploiting unconventional nonbonded interactions in drug design.

Small Molecule Disruptors of the Glucokinase-Glucokinase Regulatory Protein Interaction: 5. A Novel Aryl Sulfone Series, Optimization Through Conformational Analysis.

The glucokinase-glucokinase regulatory protein (GK-GKRP) complex plays an important role in controlling glucose homeostasis in the liver. We have recently disclosed a series of arylpiperazines as in vitro and in vivo disruptors of the GK-GKRP complex with efficacy in rodent models of type 2 diabetes mellitus (T2DM). Herein, we describe a new class of aryl sulfones as disruptors of the GK-GKRP complex, where the central piperazine scaffold has been replaced by an aromatic group. Conformational analysis and exploration of the structure-activity relationships of this new class of compounds led to the identification of potent GK-GKRP disruptors. Further optimization of this novel series delivered thiazole sulfone 93, which was able to disrupt the GK-GKRP interaction in vitro and in vivo and, by doing so, increases cytoplasmic levels of unbound GK.

Mutations in CIC and IDH1 cooperatively regulate 2-hydroxyglutarate levels and cell clonogenicity.

The majority of oligodendrogliomas (ODGs) exhibit combined losses of chromosomes 1p and 19q and mutations of isocitrate dehydrogenase (IDH1-R132H or IDH2-R172K). Approximately 70% of ODGs with 1p19q co-deletions harbor somatic mutations in the Capicua Transcriptional Repressor (CIC) gene on chromosome 19q13.2. Here we show that endogenous long (CIC-L) and short (CIC-S) CIC proteins are predominantly localized to the nucleus or cytoplasm, respectively. Cytoplasmic CIC-S is found in close proximity to the mitochondria. To study wild type and mutant CIC function and motivated by the paucity of 1p19q co-deleted ODG lines, we created HEK293 and HOG stable cell lines ectopically co-expressing CIC and IDH1. Non-mutant lines displayed increased clonogenicity, but cells co-expressing the mutant IDH1-R132H with either CIC-S-R201W or -R1515H showed reduced clonogenicity in an additive manner, demonstrating cooperative effects in our assays. Expression of mutant CIC-R1515H increased cellular 2-Hydroxyglutarate (2HG) levels compared to wild type CIC in IDH1-R132H background. Levels of phosphorylated ATP-citrate Lyase (ACLY) were lower in cell lines expressing mutant CIC-S proteins compared to cells expressing wild type CIC-S, supporting a cytosolic citrate metabolism-related mechanism bof reduced clonogenicity in our in vitro model systems. ACLY or phospho-ACLY were similarly reduced in CIC-mutant 1p19q co-deleted oligodendroglioma patient samples.

Small molecule disruptors of the glucokinase-glucokinase regulatory protein interaction: 3. Structure-activity relationships within the aryl carbinol region of the N-arylsulfonamido-N'-arylpiperazine series.

We have recently reported a novel approach to increase cytosolic glucokinase (GK) levels through the binding of a small molecule to its endogenous inhibitor, glucokinase regulatory protein (GKRP). These initial investigations culminated in the identification of 2-(4-((2S)-4-((6-amino-3-pyridinyl)sulfonyl)-2-(1-propyn-1-yl)-1-piperazinyl)phenyl)-1,1,1,3,3,3-hexafluoro-2-propanol (1, AMG-3969), a compound that effectively enhanced GK translocation and reduced blood glucose levels in diabetic animals. Herein we report the results of our expanded SAR investigations that focused on modifications to the aryl carbinol group of this series. Guided by the X-ray cocrystal structure of compound 1 bound to hGKRP, we identified several potent GK-GKRP disruptors bearing a diverse set of functionalities in the aryl carbinol region. Among them, sulfoximine and pyridinyl derivatives 24 and 29 possessed excellent potency as well as favorable PK properties. When dosed orally in db/db mice, both compounds significantly lowered fed blood glucose levels (up to 58%).

Outcomes of pancreaticoduodenectomy for pancreatic malignancy in octogenarians: an American College of Surgeons National Surgical Quality Improvement Program analysis.

BACKGROUND: Most series analyzing outcomes of pancreaticoduodenectomy in octogenarians are limited by a small sample size. The investigators used the American College of Surgeons National Surgical Quality Improvement Program database for an analysis of the impact of advanced age on outcomes after pancreatic cancer surgery. METHODS: The National Surgical Quality Improvement Program database from 2005 to 2010 was accessed to study the outcomes of 475 pancreaticoduodenectomies performed in patients ≥80 years of age compared with 4,102 patients <80 years of age using chi-square and Student's t tests. A multivariate logistic regression was used to analyze factors associated with 30-day mortality and the occurrence of major complications. RESULTS: Octogenarians had significantly more preoperative comorbidities compared with patients <80 years of age. On multivariate analysis, age ≥80 years was associated with an increased likelihood of experiencing 30-day mortality and major complications compared with patients <80 years of age. On subgroup analysis, septuagenarians had a similar odds ratio of experiencing mortality or complications compared with octogenarians, whereas patients <70 years of age were at lower risk. CONCLUSIONS: Although octogenarians have an increased risk for mortality and major complications compared with patients <80 years of age, on subgroup analysis, they do not differ from septuagenarians.

Small molecule disruptors of the glucokinase-glucokinase regulatory protein interaction: 2. Leveraging structure-based drug design to identify analogues with improved pharmacokinetic profiles.

In the previous report , we described the discovery and optimization of novel small molecule disruptors of the GK-GKRP interaction culminating in the identification of 1 (AMG-1694). Although this analogue possessed excellent in vitro potency and was a useful tool compound in initial proof-of-concept experiments, high metabolic turnover limited its advancement. Guided by a combination of metabolite identification and structure-based design, we have successfully discovered a potent and metabolically stable GK-GKRP disruptor (27, AMG-3969). When administered to db/db mice, this compound demonstrated a robust pharmacodynamic response (GK translocation) as well as statistically significant dose-dependent reductions in fed blood glucose levels.

Depression like characteristics of 5HTTLPR polymorphism and temperament in excessive internet users.

INTRODUCTION: Excessive internet use (EIU) has been reported to be comorbid with depression and the manifestation of its symptoms. This study examines the characteristics of excessive internet users that are similar to those of patients with depressive disorders in terms of serotonin transporter gene expression and harm avoidance. METHODS: 91 male adolescents with EIU and 75 healthy comparison subjects were recruited. Between group comparisons were made on genetic polymorphisms of the serotonin transport gene and with respect to novelty seeking and harm avoidance (HA) of Cloninger's Temperament Character Inventory. RESULTS: The homozygous short allelic variant of the serotonin transporter gene (SS-5HTTLPR) is more frequent in the EIU group (chi(2)=4.38, df=1, p<0.05). The HA and Beck Depression Inventory (BDI) scores were significantly higher in the EIU group than in the healthy comparison group (t=7.03, df=164, p<0.01; t=2.12, df=164, p=0.04). EIU subjects expressing SS-5HTTLPR also showed higher HA (HA1, HA2, HA4, and total HA) and Young's internet addiction scale scores than EIU subjects expressing the other serotonin transporter gene allele variants (t=2.47, df=89, p=0.01; t=2.33, df=89, p=0.02; t=2.17, df=89, p=0.03; t=2.25, df=89, p=0.03; t=2.93, df=89, p<0.01 respectively). CONCLUSIONS: The EIU group had higher SS-5HTTLPR frequencies, harm avoidance, and BDI scores. SS-5HTTLPR expression was closely related to harm avoidance in EIU. The results of this study suggest that EIU subjects may have genetic and personality traits similar to depressed patients.

Dopamine genes and reward dependence in adolescents with excessive internet video game play.

Excessive internet video game play (EIGP) has emerged as a leading cause of behavioral and developmental problems in adolescents. Recent research has implicated the role of striatal dopaminergic system in the behavioral maladaptations associated with EIGP. This study investigates the reward-dependence characteristics in EIGP adolescents as it potentially relates to genetic polymorphisms of the dopaminergic system and temperament. Seventy-nine male EIGP adolescents and 75 age- and gender-matched healthy comparison adolescents were recruited. Associations were tested with respect to the reward-dependence (RD) scale in Cloninger's Temperament and Character Inventory and the frequencies of 3 dopamine polymorphisms: Taq1A1 allele of the dopamine D2 receptor (DRD2 Taq1A1) and Val158Met in the Catecholamine-O-Methyltransferase (COMT) genes. The Taq1A1 and low activity (COMT) alleles were significantly more prevalent in the EIGP group relative to the comparison group. The present EIGP group had significantly higher RD scores than controls. Within the EIGP group, the presence of the Taq1A1 allele correlated with higher RD scores. Our findings suggest that EIGP subjects have higher reward dependency and an increased prevalence of the DRD2 Taq1A1 and COMT alleles. In particular, the DRD2 Taq1A1 allele seems to be associated with reward dependence in EIGP adolescents.