Robert Culverhouse

A genome-wide association study of alcohol dependence

Excessive alcohol consumption is one of the leading causes of preventable death in the United States. Approximately 14% of those who use alcohol meet criteria during their lifetime for alcohol dependence, which is characterized by tolerance, withdrawal, inability to stop drinking, and continued drinking despite serious psychological or physiological problems. We explored genetic influences on alcohol dependence among 1,897 European-American and African-American subjects with alcohol dependence compared with 1,932 unrelated, alcohol-exposed, nondependent controls. Constitutional DNA of each subject was genotyped using the Illumina 1M beadchip. Fifteen SNPs yielded P < 10−5, but in two independent replication series, no SNP passed a replication threshold of P < 0.05. Candidate gene GABRA2, which encodes the GABA receptor α2 subunit, was evaluated independently. Five SNPs at GABRA2 yielded nominal (uncorrected) P < 0.05, with odds ratios between 1.11 and 1.16. Further dissection of the alcoholism phenotype, to disentangle the influence of comorbid substance-use disorders, will be a next step in identifying genetic variants associated with alcohol dependence.

A Perspective on Epistasis: Limits of Models Displaying No Main Effect

The completion of a draft sequence of the human genome and the promise of rapid single-nucleotide-polymorphism-genotyping technologies have resulted in a call for the abandonment of linkage studies in favor of genome scans for association. However, there exists a large class of genetic models for which this approach will fail: purely epistatic models with no additive or dominance variation at any of the susceptibility loci. As a result, traditional association methods (such as case/control, measured genotype, and transmission/disequilibrium test [TDT]) will have no power if the loci are examined individually. In this article, we examine this class of models, delimiting the range of genetic determination and recurrence risks for two-, three-, and four-locus purely epistatic models. Our study reveals that these models, although giving rise to no additive or dominance variation, do give rise to increased allele sharing between affected sibs. Thus, a genome scan for linkage could detect genomic subregions harboring susceptibility loci. We also discuss some simple multilocus extensions of single-locus analysis methods, including a conditional form of the TDT.

Multiple Independent Loci at Chromosome 15q25.1 Affect Smoking Quantity: a Meta-Analysis and Comparison with Lung Cancer and COPD

Recently, genetic association findings for nicotine dependence, smoking behavior, and smoking-related diseases converged to implicate the chromosome 15q25.1 region, which includes the CHRNA5-CHRNA3-CHRNB4 cholinergic nicotinic receptor subunit genes. In particular, association with the nonsynonymous CHRNA5 SNP rs16969968 and correlates has been replicated in several independent studies. Extensive genotyping of this region has suggested additional statistically distinct signals for nicotine dependence, tagged by rs578776 and rs588765. One goal of the Consortium for the Genetic Analysis of Smoking Phenotypes (CGASP) is to elucidate the associations among these markers and dichotomous smoking quantity (heavy versus light smoking), lung cancer, and chronic obstructive pulmonary disease (COPD). We performed a meta-analysis across 34 datasets of European-ancestry subjects, including 38,617 smokers who were assessed for cigarettes-per-day, 7,700 lung cancer cases and 5,914 lung-cancer-free controls (all smokers), and 2,614 COPD cases and 3,568 COPD-free controls (all smokers). We demonstrate statistically independent associations of rs16969968 and rs588765 with smoking (mutually adjusted p-values<10−35 and <10−8 respectively). Because the risk alleles at these loci are negatively correlated, their association with smoking is stronger in the joint model than when each SNP is analyzed alone. Rs578776 also demonstrates association with smoking after adjustment for rs16969968 (p<10−6). In models adjusting for cigarettes-per-day, we confirm the association between rs16969968 and lung cancer (p<10−20) and observe a nominally significant association with COPD (p = 0.01); the other loci are not significantly associated with either lung cancer or COPD after adjusting for rs16969968. This study provides strong evidence that multiple statistically distinct loci in this region affect smoking behavior. This study is also the first report of association between rs588765 (and correlates) and smoking that achieves genome-wide significance; these SNPs have previously been associated with mRNA levels of CHRNA5 in brain and lung tissue.

The CHRNA5-CHRNA3-CHRNB4 Nicotinic Receptor Subunit Gene Cluster Affects Risk for Nicotine Dependence in African-Americans and in European-Americans

Genetic association studies have shown the importance of variants in the CHRNA5-CHRNA3-CHRNB4 cholinergic nicotinic receptor subunit gene cluster on chromosome 15q24-25.1 for the risk of nicotine dependence, smoking, and lung cancer in populations of European descent. We have carried out a detailed study of this region using dense genotyping in both European-Americans and African-Americans. We genotyped 75 known single nucleotide polymorphisms (SNPs) and one sequencing-discovered SNP in an African-American sample (N = 710) and in a European-American sample (N = 2,062). Cases were nicotine-dependent and controls were nondependent smokers. The nonsynonymous CHRNA5 SNP rs16969968 is the most significant SNP associated with nicotine dependence in the full sample of 2,772 subjects [P = 4.49 x 10(-8); odds ratio (OR), 1.42; 95% confidence interval (CI), 1.25-1.61] as well as in African-Americans only (P = 0.015; OR, 2.04; 1.15-3.62) and in European-Americans only (P = 4.14 x 10(-7); OR, 1.40; 1.23-1.59). Other SNPs that have been shown to affect the mRNA levels of CHRNA5 in European-Americans are associated with nicotine dependence in African-Americans but not in European-Americans. The CHRNA3 SNP rs578776, which has a low correlation with rs16969968, is associated with nicotine dependence in European-Americans but not in African-Americans. Less common SNPs (frequency <or= 5%) are also associated with nicotine dependence. In summary, multiple variants in this gene cluster contribute to nicotine dependence risk, and some are also associated with functional effects on CHRNA5. The nonsynonymous SNP rs16969968, a known risk variant in populations of European-descent, is also significantly associated with risk in African-Americans. Additional SNPs contribute to risk in distinct ways in these two populations.

Genetic studies of alcoholism and substance dependence.

During the past year, there has been a surge in genetic studies of alcohol and substance dependence. This increased interest is the result both of high repeated-heritability estimates of the addictive disorders and related phenotypes (typically ∼50%–60%) and of the appreciation that methods for detecting, mapping, and characterizing oligogenic susceptibility genes are rapidly improving. Thus, the full promise of genetics may become available to increase prevention and to improve the treatment of these extremely destructive common disorders.

Smoking and genetic risk variation across populations of European, Asian, and African American ancestry--a meta-analysis of chromosome 15q25.

Recent meta‐analyses of European ancestry subjects show strong evidence for association between smoking quantity and multiple genetic variants on chromosome 15q25. This meta‐analysis extends the examination of association between distinct genes in the CHRNA5‐CHRNA3‐CHRNB4 region and smoking quantity to Asian and African American populations to confirm and refine specific reported associations. Association results for a dichotomized cigarettes smoked per day phenotype in 27 datasets (European ancestry (N = 14,786), Asian (N = 6,889), and African American (N = 10,912) for a total of 32,587 smokers) were meta‐analyzed by population and results were compared across all three populations. We demonstrate association between smoking quantity and markers in the chromosome 15q25 region across all three populations, and narrow the region of association. Of the variants tested, only rs16969968 is associated with smoking (P < 0.01) in each of these three populations (odds ratio [OR] = 1.33, 95% CI = 1.25–1.42, P = 1.1 × 10−17 in meta‐analysis across all population samples). Additional variants displayed a consistent signal in both European ancestry and Asian datasets, but not in African Americans. The observed consistent association of rs16969968 with heavy smoking across multiple populations, combined with its known biological significance, suggests rs16969968 is most likely a functional variant that alters risk for heavy smoking. We interpret additional association results that differ across populations as providing evidence for additional functional variants, but we are unable to further localize the source of this association. Using the cross‐population study paradigm provides valuable insights to narrow regions of interest and inform future biological experiments. Genet. Epidemiol. 36:340–351, 2012.

Irinotecan pharmacogenetics: influence of pharmacodynamic genes.

Purpose: Irinotecan is an important drug for the treatment of solid tumors. Although genes involved in irinotecan pharmacokinetics have been shown to influence toxicity, there are no data on pharmacodynamic genes. CDC45L, NFKB1, PARP1, TDP1, and XRCC1 have been shown to influence the cytotoxic action of camptothecins, including irinotecan. Polymorphisms in the drug target of camptothecins, topoisomerase I (TOP1), and downstream effectors may influence patient outcomes to irinotecan therapy. We undertook a retrospective candidate gene haplotype association study to investigate this hypothesis. Experimental Design: Haplotype compositions of six candidate genes were constructed in European (n = 93), East Asian (n = 94), and West African (n = 95) populations. Haplotype-tagging single nucleotide polymorphisms (htSNP) were selected based on genealogic relationships between haplotypes. DNA samples from 107 European, advanced colorectal cancer patients treated with irinotecan-based regimens were genotyped for htSNPs as well as three coding region SNPs. Associations between genetic variants and toxicity (grade 3/4 diarrhea and neutropenia) or efficacy (objective response) were assessed. Results:TOP1 and TDP1 htSNPs were related to grade 3/4 neutropenia (P = 0.04) and response (P = 0.04), respectively. Patients homozygous for an XRCC1 haplotype (GGCC-G) were more likely to show an objective response to therapy than other patients (83% versus 30%; P = 0.02). This effect was also seen in a multivariate analysis (odds ratio, 11.9; P = 0.04). No genetic variants were associated with diarrhea. Conclusions: This is the first comprehensive pharmacogenetic investigation of irinotecan pharmacodynamic factors, and our findings suggest that genetic variation in the pharmacodynamic genes may influence the efficacy of irinotecan-containing therapies in advanced colorectal cancer patients.PURPOSE Irinotecan is an important drug for the treatment of solid tumors. Although genes involved in irinotecan pharmacokinetics have been shown to influence toxicity, there are no data on pharmacodynamic genes. CDC45L, NFKB1, PARP1, TDP1, and XRCC1 have been shown to influence the cytotoxic action of camptothecins, including irinotecan. Polymorphisms in the drug target of camptothecins, topoisomerase I (TOP1), and downstream effectors may influence patient outcomes to irinotecan therapy. We undertook a retrospective candidate gene haplotype association study to investigate this hypothesis. EXPERIMENTAL DESIGN Haplotype compositions of six candidate genes were constructed in European (n = 93), East Asian (n = 94), and West African (n = 95) populations. Haplotype-tagging single nucleotide polymorphisms (htSNP) were selected based on genealogic relationships between haplotypes. DNA samples from 107 European, advanced colorectal cancer patients treated with irinotecan-based regimens were genotyped for htSNPs as well as three coding region SNPs. Associations between genetic variants and toxicity (grade 3/4 diarrhea and neutropenia) or efficacy (objective response) were assessed. RESULTS TOP1 and TDP1 htSNPs were related to grade 3/4 neutropenia (P = 0.04) and response (P = 0.04), respectively. Patients homozygous for an XRCC1 haplotype (GGCC-G) were more likely to show an objective response to therapy than other patients (83% versus 30%; P = 0.02). This effect was also seen in a multivariate analysis (odds ratio, 11.9; P = 0.04). No genetic variants were associated with diarrhea. CONCLUSIONS This is the first comprehensive pharmacogenetic investigation of irinotecan pharmacodynamic factors, and our findings suggest that genetic variation in the pharmacodynamic genes may influence the efficacy of irinotecan-containing therapies in advanced colorectal cancer patients.

Single nucleotide polymorphism profiling across the methotrexate pathway in normal subjects and patients with rheumatoid arthritis

Methotrexate (MTX) is a commonly used disease-modifying antirheumatic drug in rheumatoid arthritis (RA). Polymorphisms occur in several genes encoding key enzymes in the folic acid pathway, which is influenced by MTX, but have not been evaluated in patients with RA. The effect of race on allele frequency has also not been evaluated. In this study, the allele frequencies of polymorphisms in six key enzymes in the MTX-folate pathway in patients with RA and healthy controls, including several common racial groups were studied. European- and African-American patients with RA and European and African healthy controls were genotyped for 22 genetic loci in six genes in the MTX cellular pathway. Differences in genotype distributions between the different racial groups were evaluated using chi(2) tests. Allele frequencies were significantly different (p < 0.001) for eight single nucleotide polymorphisms between the European and African controls. The allele frequencies of two polymorphisms showed significant differences (p < 0.001) between the African- and European-American patients with RA. Thus, racial differences exist between the allele frequencies of several polymorphisms in enzymes in the MTX-folate pathway in patients with RA and healthy controls. Whether such differences contribute to a differential response to MTX in patients with RA deserves to be investigated.

Protocol for a collaborative meta-analysis of 5-HTTLPR, stress, and depression

AbstractBackgroundDebate is ongoing about what role, if any, variation in the serotonin transporter linked polymorphic region (5-HTTLPR) plays in depression. Some studies report an interaction between 5-HTTLPR variation and stressful life events affecting the risk for depression, others report a main effect of 5-HTTLPR variation on depression, while others find no evidence for either a main or interaction effect. Meta-analyses of multiple studies have also reached differing conclusions.Methods/DesignTo improve understanding of the combined roles of 5-HTTLPR variation and stress in the development of depression, we are conducting a meta-analysis of multiple independent datasets. This coordinated approach utilizes new analyses performed with centrally-developed, standardized scripts. This publication documents the protocol for this collaborative, consortium-based meta-analysis of 5-HTTLPR variation, stress, and depression. Study eligibility criteria: Our goal is to invite all datasets, published or unpublished, with 5-HTTLPR genotype and assessments of stress and depression for at least 300 subjects. This inclusive approach is to minimize potential impact from publication bias. Data sources: This project currently includes investigators from 35 independent groups, providing data on at least N = 33,761 participants.The analytic plan was determined prior to starting data analysis. Analyses of individual study datasets will be performed by the investigators who collected the data using centrally-developed standardized analysis scripts to ensure a consistent analytical approach across sites. The consortium as a group will review and interpret the meta-analysis results.DiscussionVariation in 5-HTTLPR is hypothesized to moderate the response to stress on depression. To test specific hypotheses about the role of 5-HTTLPR variation on depression, we will perform coordinated meta-analyses of de novo results obtained from all available data, using variables and analyses determined a priori. Primary analyses, based on the original 2003 report by Caspi and colleagues of a GxE interaction will be supplemented by secondary analyses to help interpret and clarify issues ranging from the mechanism of effect to heterogeneity among the contributing studies. Publication of this protocol serves to protect this project from biased reporting and to improve the ability of readers to interpret the results of this specific meta-analysis upon its completion.

The age of necrotizing enterocolitis onset: an application of Sartwell's incubation period model.

Objective:Model age of necrotizing enterocolitis (NEC) onset applying Sartwells model of incubation periods, and examine its relationship to gestational age (GA).Study Design:Retrospective chart review of St Louis Childrens Hospital neonates diagnosed with NEC (⩾Bells stage II) from 2004 to 2008, inclusive.Result:The relationship between age of NEC (N=84 cases) onset and GA best fits a non-linear model, with infants ⩽28 weeks having a disproportionately longer time to onset than older GA groups and explained 50.3% of the variability in age of NEC onset. Additional clinical variables provided no improvement in explaining age of NEC onset. Application of Sartwells model to age of NEC onset proved a good fit, when birth is used as the common exposure episode, and age is equivalent of the incubation period.Conclusion:The relationship between day of NEC diagnosis and GA is non-linear, with lower GA infants having disproportionately longer time to onset. Despite these GA differences, the fit to Sartwells model for incubation periods model is consistent with NEC being a consequence of an event that occurs at or soon after birth.