- Open Access
Association of common variation in the PPARAgene with incident myocardial infarction in individuals with type 2 diabetes: A Go-DARTS study
Nuclear Receptor volume 3, Article number: 4 (2005)
Common variants of the PPARA gene have been found to associate with ischaemic heart disease in non diabetic men. The L162V variant was found to be protective while the C2528G variant increased risk. L162V has also been associated with altered lipid measures. We therefore sought to determine the effect of PPARA gene variation on susceptibility to myocardial infarction in patients with type 2 diabetes. 1810 subjects with type 2 diabetes from the prospective Go-DARTS study were genotyped for the L162V and C2528G variants in the PPARA gene and the association of the variants with incident non-fatal myocardial infarction was examined. Cox's proportional hazards was used to interrogate time to event from recruitment, and linear regression for analysing association of genotype with quantitative clinical traits.
The V162 allele was associated with decreased risk of non-fatal myocardial infarction (HR = 0.31, 95%CI 0.10–0.93 p = 0.037) whereas the C2528 allele was associated with increased risk (HR = 2.77 95%CI 1.34–5.75 p = 0.006). Similarly V162 was associated with a later mean age of diagnosis with type 2 diabetes and C2582 an earlier age of diagnosis. C2528 was also associated with increased total cholesterol and LDL cholesterol, which did not account for the observed increased risk. Haplotype analysis demonstrated that when both rare variants occurred on the same haplotype the effect of each was abrogated.
Genetic variation at the PPARA locus is important in determining cardiovascular risk in both male and female patients with diabetes. This genotype associated risk appears to be independent of the effect of these genotypes on lipid profiles and age of diagnosis with diabetes.
Dysregulation of fatty acid metabolism plays a pivotal role in the aetiology of type 2 diabetes , and explains, at least in part, the association between obesity, type 2 diabetes and cardiovascular disease (CVD). PPARα is a member of the nuclear receptor super-family of ligand-activated transcription factors. Ligands for PPARα include polyunsaturated fatty acids and the fibrate class of lipid-lowering drugs . It is expressed at high levels in several cell types involved in the atherosclerotic process , and its activation has beneficial effects on plasma lipids, endothelial function and markers of inflammation . Thus, the PPARA gene is a strong candidate for a genetic determinant of CVD risk in people with type 2 diabetes .
The PPARA gene has been screened for common variation [6–8]. The most studied variant is the leucine 162 valine (L162V) polymorphism, present at allele frequencies between 5 and 10%, and situated in the DNA binding domain. Functional studies have demonstrated that the V162 allele is more active in vitro [7, 9], and the V162 allele has been associated with altered plasma lipid levels [6, 8, 10], delayed progression of angiographically determined CV disease in the Lopid Coronary Angiography Trial (LOCAT), and reduced risk of ischemic heart disease in the Second Northwick Park Heart Study (NPHS2) . A second, more common, G→C variant situated in intron 7 (G2528C) is in partial allelic association with the L162V variant and shows opposing effects on cardiovascular risk and cardiac growth [9, 11]. Recently it has been demonstrated that haplotypes of these variants in association with a further A→C variant in intron 1 influence age of onset of type 2 and time to requiring insulin .
PPARα activators improve the dyslipidemia associated with type 2 diabetes and may be particularly beneficial in lowering risk of CVD in subjects with type 2 diabetes or metabolic syndrome . We therefore investigated the association between PPARA gene variation with risk of CVD and diabetes related traits in Caucasian subjects with type 2 diabetes participating in the prospective population-based Genetics of Diabetes Audit and Research in Tayside Scotland (Go-DARTS) study [14–16].
The clinical characteristics of the genotyped cohort are shown in table 1. The allele frequencies of both polymorphisms were consistent with those previously published for European non-diabetic populations (table 2). The two polymorphisms were both in Hardy-Weinberg equilibrium and were in significant linkage disequilibrium D' = 0.204 p =< 0.00001. Estimated haplotype frequencies indicated very similar values with those previously published (table 3). There was little evidence that the genotypes either singly, or when included in the model together, were associated with blood pressure or lipid measurements (table 4). V162 was associated with a generally more cardioprotective profile with V/V homozygotes having lower systolic and diastolic blood pressure, lower LDL cholesterol and higher HDL cholesterol than L/L homozygotes, although none of these differences were significant. Conversely C2528 homozygotes had a small but significantly higher total cholesterol and calculated LDL cholesterol compared to G/G homozygotes. We also found an association of genotype with age diagnosed with type 2 diabetes, with the V162 allele being associated with a significantly later age of diagnosis and the C2528 allele with a significantly earlier age of diagnosis (Table 5). When we considered haplotypes we found that V162-G2528 was associated with almost a 4 year delay in diagnosis with diabetes compared to the common L162-G2582 haplotype (p = 0.004). This association was completely abrogated when C2528 occurred together with V162 as a haplotype.
During a median follow up time 49.6 months there were 108 non-fatal myocardial infarction events and 355 deaths from all causes. In a fully adjusted Cox's proportional hazards model (table 6) we found that V162 was significantly protective against non-fatal myocardial infarction (HR 0.31, 95%CI 0.10–0.93, p = 0.037), while the C2528 variant was associated with a significantly higher risk of non-fatal myocardial infarction (HR 2.77, 95%CI 1.34–5.75, p = 0.006). This association was found to be similar in both sexes. Neither variant demonstrated any evidence of an association with risk of myocardial infarction when considered in isolation. Again, when we considered haplotypes, we found that compared to the haplotype with both common alleles, the haplotype L162-C2528 was associated with a significantly increased cardiovascular risk (HR 1.68 95%CI 1.16–2.43 p = 0.006) and the V162-G2528 a decreased risk although in this case this was not significant (HR 0.54, 95%CI 0.20–1.48, p = 0.23). Again the relative associations of each variant were completely abrogated when they occurred together on the same haplotype. The inclusion of total cholesterol in the model did not attenuate these observed associations but rather further strengthened them (V162: HR 0.28, 95% CI 0.09–0.89, p = 0.031 and C2528: HR 2.87, 95%CI 1.38–5.95, p = 0.005) demonstrating that the increased risk associated with the C2528 was not linked to raised cholesterol levels. When a combined endpoint of death from all cause, and non-fatal myocardial infarction was considered in the same model it was found that the V162 continued to demonstrate a reduced risk of an event although the association was attenuated (HR 0.52, 95%CI 0.28–0.98, p = 0.044). C2528 again demonstrated an increased risk although this was now weak and borderline non-significant (HR 1.52, 95% CI 0.99–2.31, p = 0.052).
It has been previously demonstrated that two common variants at the PPARA locus are associated with opposing risks of development of atherosclerotic vascular disease and myocardial infarction in two separate populations of non-diabetic male subjects taking part in the LOCAT and NPHS2 studies . Individuals with type 2 diabetes are however particularly susceptible to atherosclerotic macrovascular disease, and PPARα activators such as the fibrate group of drugs appear to be particularly beneficial in reducing cardiovascular events in this group of patients . In this study we have confirmed the observation that V162 is associated with a decreased risk and the C2528 variant is associated with an increased risk of cardiovascular disease and that this observation can be extended to individuals with type 2 diabetes. We also found that the association is similar in both male and female individuals. Finally we confirm a recent finding that these variants are associated with opposing influences on age of diagnosis with type 2 diabetes , and that the C2528 variant is also associated with significantly higher total cholesterol and calculated LDL cholesterol levels.
Several studies have considered the potential clinical importance of genetic variation at the PPARA locus although most have concentrated on lipid levels and have considered the L162V variant in isolation. These studies have been inconsistent indicating that L162V may influence levels of cholesterol or other lipoproteins, depending on the population analysed [6, 8, 10, 12, 17, 18], while other studies have found no evidence for such an association [19, 20]. These inconsistencies may be due in large part to differing environments, genetic background and diseased status (including medications prescribed) between the populations considered. For instance the relative concentration in the diet of saturated to polyunsaturated fat has recently been demonstrated to significantly affect association of L162Vgenotype with lipoprotein profile [21, 22]. Furthermore, it is likely that there will be differential usage of fibrate (as well as other lipid modifying) drugs between individuals with type 2 diabetes and non-diabetic populations which may also influence the lipid levels differentially by PPARα haplotype [23, 24]. Gene/gene interactions may also be important as evidenced by the observation that variants in the PPARD and APOE genes can influence the observed association [10, 20]. The inevitability of gene/environment interactions, and the observed inconsistency between studies, illustrate the difficulties of considering single measures of quantitative traits. Such measures are likely to vary considerably over an individual's lifetime, depending on health status and diet. Importantly, the clinical measures in this study were a mean of multiple measures taken over up to a three year period and therefore represent a limited integration of such temporal fluctuations.
In this study we found that the G2582C variant, but not the L162V variant, influenced lipid levels and this association was not influenced by L162V. The difference between the mean values of LDL cholesterol for GG individuals compared to CC was rather small (0.36 mmol/L) and even in this high risk population did not account for the increased cardiovascular risk associated with C2582. This was not unexpected, as in keeping with the previous studies in non-diabetic men, inclusion of total cholesterol or LDL cholesterol in the model did not affect the association with cardiovascular outcomes, indicating that the increased risk associated with the C2528 variant is not likely to be through its influence on lipid levels.
Few studies have considered cardiovascular disease or considered variants other than the L162V. One recent study also demonstrated a non significant trend towards a protective effect of V162 in individuals with diabetes . This study did not consider the G2825C variant. The present study however confirms that the V162 variant is protective against nonfatal myocardial infarction, while the C2528 variant is associated with an increased risk in this population with type 2 diabetes. These observations also appear to affect overall mortality in this population. This apparent consistency across studies with respect to cardiovascular events probably reflects the small, though global, modulation in phenotype acting across an individual's lifetime due to the slight changes in activity of PPARα associated with each variant. Unlike single measures of lipids or lipoproteins, this is less likely to be effected by temporal environmental changes. This is also likely to be true for clinical events such as age of diagnosis. In the previous study that considered age of diagnosis a further variant in intron 1 of PPARA was used to construct haplotypes with L162V and G2825C .
In this study, as in ours, the C2825 was associated with an earlier age of diagnosis, with the V162 allele being associated with protection from early diagnosis, in a manner consistent with the modulation of CVD risk. This is the first study to present the association with age-of-diagnosis of type 2 diabetes in the same population with the association with cardiovascular risk, and we can state that inclusion of age of diagnosis does not modulate the observed associations with cardiovascular risk and vice-versa, demonstrating that these are independent observations. This is not surprising, as PPARA variation is associated with CVD risk regardless of diabetic status.
The G2825C variant is in a non-coding region (intron 7) and therefore unlikely to be a directly causal variant, however several other studies have indicated its potential biological importance [9, 11, 12] and also a reduced response to fibrate therapy . These observations together with that of PPARα activation through fibrate therapy have given rise to the suggestion that the C2528 variant is associated with reduced RNA transcription and hence lower PPARα levels. Although this suggestion together with a mechanism remains to be demonstrated, it is probably due to a further (as yet) unidentified variant in or near the PPARA gene.
We have previously demonstrated a similar situation in the PPARG locus in which the biological effect of a variant with probable functional consequences is consistently modified in an opposing manner by the presence of a variant for which a function is difficult to ascribe [14, 15]. Given the role of the PPARs as master controllers of energy metabolism it is likely, as the accumulating evidence appears to suggest, that they will manifest epistasis and balanced polymorphism allowing for rapid adaptive evolutionary response to widely differing environmental challenges. Indeed, the existence of widespread, balanced variation has recently been suggested for genes involved in complex traits .
We have confirmed the potential importance of genetic variation at the PPARA locus in modulating susceptibility to cardiovascular disease, and have shown that this association is relevant to individuals with type 2 diabetes.
In the Tayside region of Scotland detailed clinical information on all individuals with diabetes mellitus is recorded on a continuously updated electronic clinical information system known as DARTS (Diabetes Audit and Research in Tayside Scotland) . Validated, region-wide electronic record-linkage techniques facilitate the identification of individuals with diabetes in the Tayside population with a sensitivity of 97% . Relevant clinical data is linked to databases containing all inpatient hospital admissions in Tayside from 1980 with diagnostic codes from ICD-9/10 (International Classification of Diseases, ninth and tenth revisions), and records of death certificates from the registrar general. This automated electronic follow-up is manually validated through a continuous cycle of review by dedicated study clinicians. Incident cardiovascular events in this population have been described previously .
Following written informed consent from individuals registered on DARTS, blood samples for genetic studies have been collected, thereby creating a genetic sub-study known as Go-DARTS. Rigorous compliance with NHS data protection and encryption standards is maintained and the study was approved by the local research ethics committee.
The PPARα L162V and G2528C genotypes were determined in 1,810 individuals all of whom were Caucasian with type 2 diabetes diagnosed between the age of 35 and 70 years. Taqman (Applied Biosystems) allelic discrimination assays were used. The primers and probes used for the allelic discrimination assays were as follows: L162V Forward primer-CAGAAACAAATGCCAGTATTGTCG, Reverse primer-GGCCACCTTACCTACCGTTGTG, L162 probe (FAM labelled) – TTCACAAGTGCCTTTCTGTCGGGATGT, V162 probe (TET labelled) – TTCACAAGTGCGTTTCTGTCGGGATGT.
G2528C Forward primer-TCCTTAAATATGGTGGAACACTTGAAG, Reverse primer-TCACAACCACCAGTTTTGCAT, G2528 probe (FAM labelled) – ATATCTAGTTTGGATTCAAAAGCTTCATTTCCCA, C2528 probe (TET labelled) – ATATCTAGTTTCGATTCAAAAGCTTCATTTCCCA.
For each clinical measure the mean was determined from multiple measures obtained up to a maximum period of three years (up to two years prior to enrolment, and up to one year following enrolment). LDL cholesterol was estimated through the use of the Freidwald equation. Linear regression was used to determine the association of genotype with each measure corrected for age at genotyping. For determining the association of genotype with age of diagnosis, this was corrected for gender and presence of a history of smoking by determining residuals and adding these to the overall mean age of diagnosis. Cox's proportional hazards was used to model time to first event. All individuals were followed from the point of genotyping until a non-fatal myocardial infarction occurred, or a composite of non-fatal myocardial infarct or all cause death. Censoring occurred either at the end of the study or death from any cause. Both variants were included in the model and it was found that a co-dominant model for each variant produced the best fit. The following variables were also included in the model; age at genotyping, gender, history of smoking, treatment with insulin, a previous history of a myocardial infarction, prevalent angina and prevalent cerebrovascular disease. Haplotype frequency estimates together with haplotype effects were determined using the THESIAS Program [28, 29]. In the case of survival analysis by haplotype using THESIAS only age at genotyping, prevalent angina and a previous history of a myocardial infarction were included in the model.
Diabetes Audit and Research in Tayside Study
Genetics of DARTS
Peroxisome Proliferator Activated Receptor alpha
- PPARA :
Gene for PPARα
- PPARD :
Gene for PPARδ
Lopoid Coronary Angiography Trial
Second Northwick Park Heart Study
McGarry JD: Banting lecture 2001: dysregulation of fatty acid metabolism in the etiology of type 2 diabetes. Diabetes. 2002, 51: 7-18.
Krey G, Braissant O, L'Horset F, Kalkhoven E, Perroud M, Parker MG, Wahli W: Fatty acids, eicosanoids, and hypolipidemic agents identified as ligands of peroxisome proliferator-activated receptors by coactivator-dependent receptor ligand assay. Mol Endocrinol. 1997, 11: 779-791. 10.1210/me.11.6.779.
Braissant O, Foufelle F, Scotto C, Dauca M, Wahli W: Differential expression of peroxisome proliferator-activated receptors (PPARs): tissue distribution of PPAR-alpha, -beta, and -gamma in the adult rat. Endocrinology. 1996, 137: 354-366. 10.1210/en.137.1.354.
Marx N, Duez H, Fruchart JC, Staels B: Peroxisome proliferator-activated receptors and atherogenesis: regulators of gene expression in vascular cells. Circ Res. 2004, 94: 1168-1178. 10.1161/01.RES.0000127122.22685.0A.
Vosper H, Khoudoli GA, Graham TL, Palmer CN: Peroxisome proliferator-activated receptor agonists, hyperlipidaemia, and atherosclerosis. Pharmacol Ther. 2002, 95: 47-62. 10.1016/S0163-7258(02)00232-2.
Flavell DM, Pineda Torra I, Jamshidi Y, Evans D, Diamond JR, Elkeles RS, Bujac SR, Miller G, Talmud PJ, Staels B, Humphries SE: Variation in the PPARalpha gene is associated with altered function in vitro and plasma lipid concentrations in Type II diabetic subjects. Diabetologia. 2000, 43: 673-680. 10.1007/s001250051357.
Sapone A, Peters JM, Sakai S, Tomita S, Papiha SS, Dai R, Friedman FK, Gonzalez FJ: The human peroxisome proliferator-activated receptor alpha gene: identification and functional characterization of two natural allelic variants. Pharmacogenetics. 2000, 10: 321-333. 10.1097/00008571-200006000-00006.
Vohl MC, Lepage P, Gaudet D, Brewer CG, Betard C, Perron P, Houde G, Cellier C, Faith JM, Despres JP, Morgan K, Hudson TJ: Molecular scanning of the human PPARa gene: association of the L162v mutation with hyperapobetalipoproteinemia. J Lipid Res. 2000, 41: 945-952.
Jamshidi Y, Montgomery HE, Hense HW, Myerson SG, Torra IP, Staels B, World MJ, Doering A, Erdmann J, Hengstenberg C, Humphries SE, Schunkert H, Flavell DM: Peroxisome proliferator--activated receptor alpha gene regulates left ventricular growth in response to exercise and hypertension. Circulation. 2002, 105: 950-955. 10.1161/hc0802.104535.
Tai ES, Demissie S, Cupples LA, Corella D, Wilson PW, Schaefer EJ, Ordovas JM: Association between the PPARA L162V polymorphism and plasma lipid levels: the Framingham Offspring Study. Arterioscler Thromb Vasc Biol. 2002, 22: 805-810. 10.1161/01.ATV.0000012302.11991.42.
Flavell DM, Jamshidi Y, Hawe E, Pineda Torra I, Taskinen MR, Frick MH, Nieminen MS, Kesaniemi YA, Pasternack A, Staels B, Miller G, Humphries SE, Talmud PJ, Syvanne M: Peroxisome proliferator-activated receptor alpha gene variants influence progression of coronary atherosclerosis and risk of coronary artery disease. Circulation. 2002, 105: 1440-1445. 10.1161/01.CIR.0000012145.80593.25.
Flavell DM, Ireland H, Stephens JW, Hawe E, Acharya J, Mather H, Hurel SJ, Humphries SE: Peroxisome proliferator-activated receptor alpha gene variation influences age of onset and progression of type 2 diabetes. Diabetes. 2005, 54: 582-586.
Robins SJ: Cardiovascular disease with diabetes or the metabolic syndrome: should statins or fibrates be first line lipid therapy?. Curr Opin Lipidol. 2003, 14: 575-583. 10.1097/00041433-200312000-00005.
Doney AS, Fischer B, Leese G, Morris AD, Palmer CN: Cardiovascular risk in type 2 diabetes is associated with variation at the PPARG locus: a Go-DARTS study. Arterioscler Thromb Vasc Biol. 2004, 24: 2403-2407. 10.1161/01.ATV.0000147897.57527.e4.
Doney AS, Fischer B, Cecil JE, Boylan K, McGuigan FE, Ralston SH, Morris AD, Palmer CN: Association of the Pro12Ala and C1431T variants of PPARG and their haplotypes with susceptibility to Type 2 diabetes. Diabetologia. 2004, 47: 555-558. 10.1007/s00125-003-1323-1.
Morris AD, Boyle DI, MacAlpine R, Emslie-Smith A, Jung RT, Newton RW, MacDonald TM: The diabetes audit and research in Tayside Scotland (DARTS) study: electronic record linkage to create a diabetes register. DARTS/MEMO Collaboration.[see comment]. BMJ. 1997, 315: 524-528.
Lacquemant C, Lepretre F, Pineda Torra I, Manraj M, Charpentier G, Ruiz J, Staels B, Froguel P: Mutation screening of the PPARalpha gene in type 2 diabetes associated with coronary heart disease. Diabetes & Metabolism. 2000, 26: 393-401.
Jamshidi Y, Flavell DM, Hawe E, MacCallum PK, Meade TW, Humphries SE: Genetic determinants of the response to bezafibrate treatment in the lower extremity arterial disease event reduction (LEADER) trial. Atherosclerosis. 2002, 163: 183-192. 10.1016/S0021-9150(02)00002-3.
Gouni-Berthold I, Giannakidou E, Muller-Wieland D, Faust M, Kotzka J, Berthold HK, Krone W: Association between the PPARalpha L162V polymorphism, plasma lipoprotein levels, and atherosclerotic disease in patients with diabetes mellitus type 2 and in nondiabetic controls. American Heart Journal. 2004, 147: 1117-1124. 10.1016/j.ahj.2003.12.005.
Skogsberg J, Kannisto K, Cassel TN, Hamsten A, Eriksson P, Ehrenborg E: Evidence that peroxisome proliferator-activated receptor delta influences cholesterol metabolism in men. Arterioscler Thromb Vasc Biol. 2003, 23 (4): 637-643. 10.1161/01.ATV.0000064383.88696.24.
Paradis AM, Fontaine-Bisson B, Bosse Y, Robitaille J, Lemieux S, Jacques H, Lamarche B, Tchernof A, Couture P, Vohl MC: The peroxisome proliferator-activated receptor alpha Leu162Val polymorphism influences the metabolic response to a dietary intervention altering fatty acid proportions in healthy men. American Journal of Clinical Nutrition. 2005, 81: 523-530.
Tai ES, Corella D, Demissie S, Cupples LA, Coltell O, Schaefer EJ, Tucker KL, Ordovas JM: Polyunsaturated fatty acids interact with the PPARA-L162V polymorphism to affect plasma triglyceride and apolipoprotein C-III concentrations in the Framingham Heart Study. J Nutr. 2005, 135: 397-403.
Bosse Y, Pascot A, Dumont M, Brochu M, Prud'homme D, Bergeron J, Despres JP, Vohl MC: Influences of the PPAR alpha-L162V polymorphism on plasma HDL(2)-cholesterol response of abdominally obese men treated with gemfibrozil. Genetics in Medicine. 2002, 4: 311-315. 10.1097/00125817-200207000-00010.
Brisson D, Ledoux K, Bosse Y, St-Pierre J, Julien P, Perron P, Hudson TJ, Vohl MC, Gaudet D: Effect of apolipoprotein E, peroxisome proliferator-activated receptor alpha and lipoprotein lipase gene mutations on the ability of fenofibrate to improve lipid profiles and reach clinical guideline targets among hypertriglyceridemic patients. Pharmacogenetics. 2002, 12: 313-320. 10.1097/00008571-200206000-00007.
Foucher C, Rattier S, Flavell DM, Talmud PJ, Humphries SE, Kastelein JJ, Ayyobi A, Pimstone S, Frohlich J, Ansquer JC, Steiner G, investigators D: Response to micronized fenofibrate treatment is associated with the peroxisome-proliferator-activated receptors alpha G/C intron7 polymorphism in subjects with type 2 diabetes. Pharmacogenetics. 2004, 14: 823-829. 10.1097/00008571-200412000-00005.
Kroymann J, Mitchell-Olds T: Epistasis and balanced polymorphism influencing complex trait variation. Nature. 2005, 435: 95-98. 10.1038/nature03480.
Evans JM, Wang J, Morris AD: Comparison of cardiovascular risk between patients with type 2 diabetes and those who had had a myocardial infarction: cross sectional and cohort studies.[erratum appears in BMJ 2002 Jun 8;324(7350):1357]. BMJ. 2002, 324: 939-942. 10.1136/bmj.324.7343.939.
Tregouet DA, Escolano S, Tiret L, Mallet A, Golmard JL: A new algorithm for haplotype-based association analysis: the Stochastic-EM algorithm. Ann Hum Genet. 2004, 68: 165-177. 10.1046/j.1529-8817.2003.00085.x.
Tregouet DA, Tiret L: Cox proportional hazards survival regression in haplotype-based association analysis using the Stochastic-EM algorithm. Eur J Hum Genet. 2004, 12: 971-974. 10.1038/sj.ejhg.5201238.
The Go-DARTS recruitment was funded by an anonymous trust donations to Tenovus Tayside. Colin Palmer and Andrew Morris are supported by the Scottish Higher Education Funding Council (SHEFC) and the Scottish Executive Genetic Health Initiative. Go-DARTS is currently supported by the Wellcome Trust.
The author(s) declare that they have no competing interests.
ASFD and CNAP wrote the manuscript and performed the analysis, CNAP and ADM conceived the study and participated in its design and coordination. GL helped draft the manuscript and contributed to the data analysis. BF and SL performed the genotyping. All authors read and approved the final manuscript.