Volume 6 Supplement 1
Effect of genotype × alcoholism interaction on linkage analysis of an alcoholism-related quantitative phenotype
© Arya et al; licensee BioMed Central Ltd 2005
Published: 30 December 2005
Studies have shown that genetic and environmental factors and their interactions affect several alcoholism phenotypes. Genotype × alcoholism (G×A) interaction refers to the environmental (alcoholic and non-alcoholic) influences on the autosomal genes contributing to variation in an alcoholism-related quantitative phenotype. The purpose of this study was to examine the effects of G×A interaction on the detection of linkage for alcoholism-related phenotypes.
We used phenotypic and genotypic data from the Collaborative Study on the Genetics of Alcoholism relating to 1,388 subjects as part of Genetic Analysis Workshop 14 problem 1. We analyzed the MXDRNK phenotype to detect G×A interaction using SOLAR. Upon detecting significant interaction, we conducted variance-component linkage analyses using microsatellite marker data. For maximum number of drinks per a 24 hour period, the highest LODs were observed on chromosomes 1, 4, and 13 without G×A interaction. Interaction analysis yielded four regions on chromosomes 1, 4, 13, and 15. On chromosome 4, a maximum LOD of 1.5 at the same location as the initial analysis was obtained after incorporating G×A interaction effects. However, after correcting for extra parameters, the LOD score was reduced to a corrected LOD of 1.1, which is similar to the LOD observed in the non-interaction analysis. Thus, we see little differences in LOD scores, while some linkage regions showed large differences in the magnitudes of estimated quantitative trait loci heritabilities between the alcoholic and non-alcoholic groups. These potential hints of differences in genetic effect may influence future analyses of variants under these linkage peaks.
Family, twin, and adoption studies have indicated that genetic and environmental factors and their interactions contribute to the development of alcoholism [1–3]. Several studies have demonstrated the importance of considering environment-specific major gene effects on different phenotypes. Genotype × alcoholism (G×A) interaction refers to the environmental (alcoholic and non-alcoholic) influences on the autosomal genes contributing to variation in an alcoholism-related quantitative phenotype. Alcoholic environment refers to chronic alcohol ingestion. This may interact with gene expression in a number of different ways; affected individuals carry a different spectrum of genetic variants; or developmental differences between individuals, who are or are not at risk for alcoholism, affect gene expression; or ingestion of alcohol affects gene expression. Although it is well documented that alcoholism-related traits have strong genetic determinants, few susceptibility genes influencing these complex disease phenotypes have been identified. Because alcoholism is a complex phenotype influenced by several genes with small effects, it is difficult to detect such genes. Hence it may be helpful to examine potentially simpler endophenotypes related to disease risk . On the other hand, it may be easier to detect such susceptibility genes if they have major effect on related quantitative phenotypes [5–7]. In addition, gene × environment (G×E) interaction has been detected in quantitative genetic analyses of a variety of traits such as serum lipid concentrations and event-related evoked potentials (ERPs) [8, 9]. Furthermore, G×E interactions (e.g., genotype × age, genotype × sex, and genotype × diet) in a given quantitative trait (e.g., body composition and ERP phenotypes) under the assumption of polygenic inheritance has been considered an important component in modeling environment-specific effects for polygenic variance components and major genes [9–12]. Therefore, in this study, we examined the effects of G×A interactions on the linkage analysis of a quantitative phenotype from Collaborative Study on the Genetics of Alcoholism (COGA) data, maximum number of drinks per a 24 hour period (MXDRNK), which is a correlate of alcoholism and is expected to reflect individual's ability to metabolize alcohol as well as the effect of social environment. By using alcoholism as an environment in G×E analyses of MXDRNK (i.e., we are referring to internal/within individual environment but not family environment), we are essentially allowing for the possibility that the magnitude or source of genetic effects on variation in alcohol consumption may differ in alcoholics and non-alcoholics.
Subjects and Methods
In this study, the Genetic Analysis Workshop 14 (GAW14) COGA data (Problem 1) consisting of 1,388 family members, have been analyzed. Prior to the analysis we recoded the affection status based on the definition of alcoholism according to COGA as well as DSM-IV criteria in two ways: diagnoses 1 and 2 correspond to COGA and DSM-IV and that a includes individuals with some symptoms as unaffected (diagnoses COGA-Aldxla and DSM-IV-Aldx2a), whereas b considers them unknowns (diagnoses COGA-Aldx1b and DSM-IV-Aldx2b). In the analysis of the GAW14 COGA data, we used a maximum likelihood variance components approach for the study of G×E interaction using related individuals in different environments . To minimize the problem of non-normality, MXDRNK values were log transformed. In this interaction model, two additional parameters are modeled: a) environment-specific genetic variances, and b) a genetic correlation between groups of individuals living in different environments. A significant G×E interaction is indicated by significantly different magnitudes of genetic variances for individuals living in different exposure groups (alcoholics vs. non-alcoholics), and/or a genetic correlation (ρG) is less than 1 between exposure groups. In an interaction model, assuming the probability of an individual having a specific polygenotype is independent of environment, the expected additive genetic covariance between a pair of alcoholics is COV(alc,alc) = 2Φ σ2Galc or the covariance between a pair of non-alcoholics would be COV(noalc,noalc) = 2Φ σ2Gnalc, i.e., 2Φ times the appropriate genetic variance.
The covariance between an alcoholic individual and a nonalcoholic one is modeled as:
COV(alc,nalc) = 2Φ σGalc σGnalc ρG,
where alc and nalc denote alcoholics and non-alcoholics, respectively, Φ is the coefficient of kinship between the two individuals, ρG is the additive genetic correlation between the expression of the phenotype in the two environments, and σGalc and σGnalc are the additive genetic standard deviations of alcoholics and nonalcoholics, respectively. We first screened for the presence of G×A interaction in several quantitative phenotypes, including MXDRNK, using a quantitative genetic method. After detecting significant G×A interaction for MXDRNK phenotype, we performed variance-component linkage analyses with a customized model to include diagnosis-specific quantitative trait loci (QTL) effects and using microsatellite marker multipoint identity by descent (MIBD) matrices estimated using LOKI. The customized variance component linkage model may be defined as
COV(alc,nalc) = Π σQalc σQnalc + 2Φ σGalc σGnalc ρG.
This model has an additional QTL variance as compared with the standard linkage model. Corrected LOD (LODc) scores assume that σQalc σQnalc are independent under the null, producing a test statistic distribution of ¼χ22, ½χ21, ¼ point mass at 0. This assumption may be overly conservative. These analytical techniques were implemented using the computer program SOLAR .
p-Values for the two tests of G×A interaction (gsdalc = gsdnalc and ρG = 1) for several quantitative phenotypes.
gsdalc = gsdnalc
ρG = 1
gsdalc = gsdnalc
ρG = 1
gsdalc = gsdnalc
ρG = 1
gsdalc = gsdnalc
ρG = 1
gsdalc = gsdnalc
ρG = 1
gsdalc = gsdnalc
ρG = 1
gsdalc = gsdnalc
ρG = 1
Descriptive statistics of MXDRNK phenotype according to diagnostic criteria and affection status
Chromosomal regions linked to MXDRNK phenotype in COGA data
Environment may influence the variation in the expression of genes influencing a variety of phenotypes including alcohol-related phenotypes. Genotype × alcoholism interaction was explored for a variety of quantitative phenotypes in the COGA dataset but was detected only for the MXDRNKs phenotype. Results of G×A analyses were consistent across the COGA and DSM-IV alcoholism diagnoses, but were affected by the categorization of individuals who were unaffected with some symptoms. Results were generally stronger when these individuals were categorized as unaffected, rather than unknown. However, this result may be strongly influenced by sample size considerations as the addition of the "unaffected with some symptoms" more than doubled the size of the nonalcoholic group.
Our analysis also shows that accounting for G×E interaction may increase the linkage signal. For MXDRNK, interaction analysis failed to show evidence at the implicated region on chromosome 1 in the non-G×A analysis but yielded a slightly increased linkage signal at a different location (238 cM), which corresponds to the previously reported linkage with factor 2, a factor analysis-derived trait defined by harm avoidance, novelty seeking, and age of onset of drinking . On chromosome 4, although the increase in LOD score is only slightly higher than the LOD obtained in the linkage model without interaction, the observed environment-specific QTL effects are interesting and are consistent with the previous observation that inclusion of unaffected individuals is crucial to detection of linkage to this chromosome 4 region. Moreover, chromosome 4 QTL appears to influence unaffecteds only and this observation is also consistent with COGA findings Furthermore, this region is important because this QTL has an impact on drinking behavior in non-alcoholics and it is consistent with that identified in other studies in the literature [5, 16]. Interestingly, interaction analysis has also yielded a new region on chromosome 15 (max LOD = 2.04, LODc = 1.6, 100 cM), which is only 25 cM away from a linkage with factor 2  and QTL effect is stronger in affecteds. On the other hand, no signal was observed at the corresponding region in non-G×A analysis. In contrast, the evidence for linkage has been reduced on chromosome 13 in G×A interaction analysis and the QTL effect appears to be stronger in unaffecteds.
In conclusion, genotype × alcoholism interaction analysis yielded interesting results. Drinking behavior appears to be influenced by environment-specific genes in both alcoholics and non-alcoholics. The implicated regions on chromosomes 1, 4, and 15 are consistent with previously reported linkage findings. These results indicate that further analyses may benefit from considering the possibility of differing genetic effects in alcoholics and non-alcoholics, for example by stratifying analysis on alcoholism diagnosis.
Collaborative Study on the Genetics of Alcoholism
Event-related evoked potentials
Genotype × alcoholism
Genotype × environment
Genetic Analysis Workshop 14
Multipoint identity by descent
Maximum number of drinks per a 24-hour period
Quantitative trait loci
This work is supported in part by NIH R01MH59490.
- Heath AC, Bucholz KK, Madden PAF, Dinwiddie SH, Slutske WS, Bierut LJ, Statham DJ, Dunne MP, Whitfield JB, Martin NG: Genetic and environmental contributions to alcohol dependence risk in a national twin sample: consistency of findings in women and men. Psychol Med. 1997, 27: 1381-1396. 10.1017/S0033291797005643.View ArticlePubMedGoogle Scholar
- McGue M: A behavioral-genetic perspective on children of alcoholics. Alcohol Health Res World. 1997, 21: 210-217.PubMedGoogle Scholar
- Schuckit MA: Genetics of the risk for alcoholism. Am J Addict. 2000, 9: 103-112. 10.1080/10550490050173172.View ArticlePubMedGoogle Scholar
- Porjesz B, Begleiter H, Wang K, Almasy L, Chorlian DB, Stimus AT, Kuperman S, O'Connor SJ, Rohrbaugh J, Bauer LO, Edenberg HJ, Goate A, Rice JP, Reich T: Linkage and linkage disequilibrium mapping of ERP and EEG phenotypes. Biol Psychol. 2000, 61: 229-248. 10.1016/S0301-0511(02)00060-1.View ArticleGoogle Scholar
- Saccone NL, Kwon JM, Corbett J, Goate A, Rochberg N, Edenberg HJ, Foroud T, Li T-K, Begleiter H, Reich T, Rice JP: A genome screen of maximum number of drinks as an alcoholism phenotype. Am J Med Genet. 2000, 96: 632-637. 10.1002/1096-8628(20001009)96:5<632::AID-AJMG8>3.0.CO;2-#.View ArticlePubMedGoogle Scholar
- Mahaney MC, Blangero J, Comuzzie AG, VandeBerg JL, Stern MP, MacCluer JW: Plasma HDL-cholesterol, triglycerides, and adiposity. A quantitative genetic test of the conjoint trait hypothesis in the San Antonio Family Heart Study. Circulation. 1995, 92: 3240-3248.View ArticlePubMedGoogle Scholar
- Arya R, Duggirala R, Williams JT, Almasy L, Blangero J: Power to localize the major gene for disease liability is increased after accounting for the effects of related quantitative phenotypes. Genet Epidemiol. 2001, 21: S774-S778.PubMedGoogle Scholar
- Towne B, Blangero J, Siervogel RM: Genotype-by-sex interaction in measures of lipids, lipoproteins, and apolipoproteins. Genet Epidemiol. 1993, 10: 107-117. 10.1002/gepi.1370100647.View ArticleGoogle Scholar
- Towne B, Almasy L, Siervogel RM, Blangero J: Effects of genotype-by-sex interaction on linkage analysis of visual event-related evoked potentials. Genet Epidemiol. 1999, 17: S355-S360.View ArticlePubMedGoogle Scholar
- Blangero J, MacCluer JW, Kammerer CM, Mott CE: Genetic analysis of apolipoprotein A-I in two dietary environments. Am J Hum Genet. 1990, 47: 414-428.PubMed CentralPubMedGoogle Scholar
- Comuzzie AG, Blangero J, Mahaney MC, Mitchell BM, Hixson JE, Samollow PB, Stern MP, MacCluer JW: Major gene with sex-specific effects influences fat mass in Mexican Americans. Genet Epidemiol. 1995, 12: 475-488. 10.1002/gepi.1370120505.View ArticlePubMedGoogle Scholar
- Towne B, Blangero J, Siervogel RM: Effects of genotype-by-sex interaction on quantitative trait linkage analysis. Genet Epidemiol. 1997, 14: 1053-1058. 10.1002/(SICI)1098-2272(1997)14:6<1053::AID-GEPI82>3.0.CO;2-G.View ArticlePubMedGoogle Scholar
- Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 1198-1211. 10.1086/301844.PubMed CentralView ArticlePubMedGoogle Scholar
- Warren DM, Dyer TD, Peterson CP, Mahaney MC, Blangero J, Almasy L: A comparison of univariate, bivariate, and trivariate whole-genome linkage screens of genetically correlated electrophysiological endophenotypes. BMC Genet. 2005, 6 (Suppl 1): S117-10.1186/1471-2156-6-S1-S117.PubMed CentralView ArticlePubMedGoogle Scholar
- Dick DM, Nurnberger J, Edenberg HJ, Goate A, Crowe R, Rice J, Bucholz KK, Kramer J, Schuckit , Smith TL, Porjesz B, Begleiter H, Hesselbrock V, Foroud T: Suggestive linkage on chromosome 1 for a quantitative alcohol-related phenotype. Alcohol Clin Exp Res. 2002, 26: 1453-1460. 10.1111/j.1530-0277.2002.tb02443.x.View ArticlePubMedGoogle Scholar
- Williams JT, Begleiter H, Porjesz B, Edenberg HJ, Foroud T, Reich T, Goate A, Van Eerdewegh PV, Almasy L, Blangero J: Joint multipoint linkage analysis of multivariate qualitative and quantitative traits. II. Alcoholism and event-related potentials. Am J Hum Genet. 1999, 65: 1148-1160. 10.1086/302571.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.