The impact of adjusting for baseline in pharmacogenomic genome-wide association studies of quantitative change.

In pharmacogenomic studies of quantitative change, any association between genetic variants and the pretreatment (baseline) measurement can bias the estimate of effect between those variants and drug response. A putative solution is to adjust for baseline. We conducted a series of genome-wide association studies (GWASs) for low-density lipoprotein cholesterol (LDL-C) response to statin therapy in 34,874 participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort as a case study to investigate the impact of baseline adjustment on results generated from pharmacogenomic studies of quantitative change. Across phenotypes of statin-induced LDL-C change, baseline adjustment identified variants from six loci meeting genome-wide significance (SORT/CELSR2/PSRC1, LPA, SLCO1B1, APOE, APOB, and SMARCA4/LDLR). In contrast, baseline-unadjusted analyses yielded variants from three loci meeting the criteria for genome-wide significance (LPA, APOE, and SLCO1B1). A genome-wide heterogeneity test of baseline versus statin on-treatment LDL-C levels was performed as the definitive test for the true effect of genetic variants on statin-induced LDL-C change. These findings were generally consistent with the models not adjusting for baseline signifying that genome-wide significant hits generated only from baseline-adjusted analyses (SORT/CELSR2/PSRC1, APOB, SMARCA4/LDLR) were likely biased. We then comprehensively reviewed published GWASs of drug-induced quantitative change and discovered that more than half (59%) inappropriately adjusted for baseline. Altogether, we demonstrate that (1) baseline adjustment introduces bias in pharmacogenomic studies of quantitative change and (2) this erroneous methodology is highly prevalent. We conclude that it is critical to avoid this common statistical approach in future pharmacogenomic studies of quantitative change

Association of FADS1/2 Locus Variants and Polyunsaturated Fatty Acids With Aortic Stenosis.

IMPORTANCE: Aortic stenosis (AS) has no approved medical treatment. Identifying etiological pathways for AS could identify pharmacological targets. OBJECTIVE: To identify novel genetic loci and pathways associated with AS. DESIGN, SETTING, AND PARTICIPANTS: This genome-wide association study used a case-control design to evaluate 44 703 participants (3469 cases of AS) of self-reported European ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (from January 1, 1996, to December 31, 2015). Replication was performed in 7 other cohorts totaling 256 926 participants (5926 cases of AS), with additional analyses performed in 6942 participants from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. Follow-up biomarker analyses with aortic valve calcium (AVC) were also performed. Data were analyzed from May 1, 2017, to December 5, 2019. EXPOSURES: Genetic variants (615 643 variants) and polyunsaturated fatty acids (ω-6 and ω-3) measured in blood samples. MAIN OUTCOMES AND MEASURES: Aortic stenosis and aortic valve replacement defined by electronic health records, surgical records, or echocardiography and the presence of AVC measured by computed tomography. RESULTS: The mean (SD) age of the 44 703 GERA participants was 69.7 (8.4) years, and 22 019 (49.3%) were men. The rs174547 variant at the FADS1/2 locus was associated with AS (odds ratio [OR] per C allele, 0.88; 95% CI, 0.83-0.93; P = 3.0 × 10-6), with genome-wide significance after meta-analysis with 7 replication cohorts totaling 312 118 individuals (9395 cases of AS) (OR, 0.91; 95% CI, 0.88-0.94; P = 2.5 × 10-8). A consistent association with AVC was also observed (OR, 0.91; 95% CI, 0.83-0.99; P = .03). A higher ratio of arachidonic acid to linoleic acid was associated with AVC (OR per SD of the natural logarithm, 1.19; 95% CI, 1.09-1.30; P = 6.6 × 10-5). In mendelian randomization, increased FADS1 liver expression and arachidonic acid were associated with AS (OR per unit of normalized expression, 1.31 [95% CI, 1.17-1.48; P = 7.4 × 10-6]; OR per 5-percentage point increase in arachidonic acid for AVC, 1.23 [95% CI, 1.01-1.49; P = .04]; OR per 5-percentage point increase in arachidonic acid for AS, 1.08 [95% CI, 1.04-1.13; P = 4.1 × 10-4]). CONCLUSIONS AND RELEVANCE: Variation at the FADS1/2 locus was associated with AS and AVC. Findings from biomarker measurements and mendelian randomization appear to link ω-6 fatty acid biosynthesis to AS, which may represent a therapeutic target

The impact of adjusting for baseline in pharmacogenomic genome-wide association studies of quantitative change.

In pharmacogenomic studies of quantitative change, any association between genetic variants and the pretreatment (baseline) measurement can bias the estimate of effect between those variants and drug response. A putative solution is to adjust for baseline. We conducted a series of genome-wide association studies (GWASs) for low-density lipoprotein cholesterol (LDL-C) response to statin therapy in 34,874 participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort as a case study to investigate the impact of baseline adjustment on results generated from pharmacogenomic studies of quantitative change. Across phenotypes of statin-induced LDL-C change, baseline adjustment identified variants from six loci meeting genome-wide significance (SORT/CELSR2/PSRC1, LPA, SLCO1B1, APOE, APOB, and SMARCA4/LDLR). In contrast, baseline-unadjusted analyses yielded variants from three loci meeting the criteria for genome-wide significance (LPA, APOE, and SLCO1B1). A genome-wide heterogeneity test of baseline versus statin on-treatment LDL-C levels was performed as the definitive test for the true effect of genetic variants on statin-induced LDL-C change. These findings were generally consistent with the models not adjusting for baseline signifying that genome-wide significant hits generated only from baseline-adjusted analyses (SORT/CELSR2/PSRC1, APOB, SMARCA4/LDLR) were likely biased. We then comprehensively reviewed published GWASs of drug-induced quantitative change and discovered that more than half (59%) inappropriately adjusted for baseline. Altogether, we demonstrate that (1) baseline adjustment introduces bias in pharmacogenomic studies of quantitative change and (2) this erroneous methodology is highly prevalent. We conclude that it is critical to avoid this common statistical approach in future pharmacogenomic studies of quantitative change

Assessing the clinical impact of CYP2C9 pharmacogenetic variation on phenytoin prescribing practice and patient response in an integrated health system.

ObjectiveTo assess the impact of CYP2C9 variation on phenytoin patient response and clinician prescribing practice where genotype was unknown during treatment.MethodsA retrospective analysis of Resource on Genetic Epidemiology Research on Adult Health and Aging cohort participants who filled a phenytoin prescription between 1996 and 2017. We used laboratory test results, medication dispensing records, and medical notes to identify associations of CYP2C9 genotype with phenytoin blood concentration, neurologic side effects, and medication dispensing patterns reflecting clinician prescribing practice and patient response.ResultsAmong 993 participants, we identified 69% extensive, 20% high-intermediate, 10% low-intermediate, and 2% poor metabolizers based on CYP2C9 genotypes. Compared with extensive metabolizer genotype, low-intermediate/poor metabolizer genotype was associated with increased dose-adjusted phenytoin blood concentration [21.3 pg/mL, 95% confidence interval (CI): 13.6-29.0 pg/mL; P < 0.01] and increased risk of neurologic side effects (hazard ratio: 2.40, 95% CI: 1.24-4.64; P < 0.01). Decreased function CYP2C9 genotypes were associated with medication dispensing patterns indicating dose decrease, use of alternative anticonvulsants, and worse adherence, although these associations varied by treatment indication for phenytoin.ConclusionCYP2C9 variation was associated with clinically meaningful differences in clinician prescribing practice and patient response, with potential implications for healthcare utilization and treatment efficacy

A large electronic-health-record-based genome-wide study of serum lipids

A genome-wide association study (GWAS) of 94,674 ancestrally diverse Kaiser Permanente members using 478,866 longitudinal electronic health record (EHR)-derived measurements for untreated serum lipid levels empowered multiple new findings: 121 new SNP associations (46 primary, 15 conditional, and 60 in meta-analysis with Global Lipids Genetic Consortium data); an increase of 33-42% in variance explained with multiple measurements; sex differences in genetic impact (greater impact in females for LDL, HDL, and total cholesterol and the opposite for triglycerides); differences in variance explained among non-Hispanic whites, Latinos, African Americans, and East Asians; genetic dominance and epistatic interaction, with strong evidence for both at the ABO and FUT2 genes for LDL; and tissue-specific enrichment of GWAS-associated SNPs among liver, adipose, and pancreas eQTLs. Using EHR pharmacy data, both LDL and triglyceride genetic risk scores (477 SNPs) were strongly predictive of age at initiation of lipid-lowering treatment. These findings highlight the value of longitudinal EHRs for identifying new genetic features of cholesterol and lipoprotein metabolism with implications for lipid treatment and risk of coronary heart disease

A Large-Scale Association Study Detects Novel Rare Variants, Risk Genes, Functional Elements, and Polygenic Architecture of Prostate Cancer Susceptibility

To identify rare variants associated with prostate cancer susceptibility and better characterize the mechanisms and cumulative disease risk associated with common risk variants, we conducted an integrated study of prostate cancer genetic etiology in two cohorts using custom genotyping microarrays, large imputation reference panels, and functional annotation approaches. Specifically, 11,984 men (6,196 prostate cancer cases and 5,788 controls) of European ancestry from Northern California Kaiser Permanente were genotyped and meta-analyzed with 196,269 men of European ancestry (7,917 prostate cancer cases and 188,352 controls) from the UK Biobank. Three novel loci, including two rare variants (European ancestry minor allele frequency < 0.01, at 3p21.31 and 8p12), were significant genome wide in a meta-analysis. Gene-based rare variant tests implicated a known prostate cancer gene (HOXB13), as well as a novel candidate gene (ILDR1), which encodes a receptor highly expressed in prostate tissue and is related to the B7/CD28 family of T-cell immune checkpoint markers. Haplotypic patterns of long-range linkage disequilibrium were observed for rare genetic variants at HOXB13 and other loci, reflecting their evolutionary history. In addition, a polygenic risk score (PRS) of 188 prostate cancer variants was strongly associated with risk (90th vs. 40th-60th percentile OR = 2.62, P = 2.55 × 10-191). Many of the 188 variants exhibited functional signatures of gene expression regulation or transcription factor binding, including a 6-fold difference in log-probability of androgen receptor binding at the variant rs2680708 (17q22). Rare variant and PRS associations, with concomitant functional interpretation of risk mechanisms, can help clarify the full genetic architecture of prostate cancer and other complex traits. SIGNIFICANCE: This study maps the biological relationships between diverse risk factors for prostate cancer, integrating different functional datasets to interpret and model genome-wide data from over 200,000 men with and without prostate cancer.See related commentary by Lachance, p. 1637

Association of FADS1/2 Locus Variants and Polyunsaturated Fatty Acids With Aortic Stenosis.

Importance: Aortic stenosis (AS) has no approved medical treatment. Identifying etiological pathways for AS could identify pharmacological targets. Objective: To identify novel genetic loci and pathways associated with AS. Design, Setting, and Participants: This genome-wide association study used a case-control design to evaluate 44 703 participants (3469 cases of AS) of self-reported European ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (from January 1, 1996, to December 31, 2015). Replication was performed in 7 other cohorts totaling 256 926 participants (5926 cases of AS), with additional analyses performed in 6942 participants from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. Follow-up biomarker analyses with aortic valve calcium (AVC) were also performed. Data were analyzed from May 1, 2017, to December 5, 2019. Exposures: Genetic variants (615 643 variants) and polyunsaturated fatty acids (ω-6 and ω-3) measured in blood samples. Main Outcomes and Measures: Aortic stenosis and aortic valve replacement defined by electronic health records, surgical records, or echocardiography and the presence of AVC measured by computed tomography. Results: The mean (SD) age of the 44 703 GERA participants was 69.7 (8.4) years, and 22 019 (49.3%) were men. The rs174547 variant at the FADS1/2 locus was associated with AS (odds ratio [OR] per C allele, 0.88; 95% CI, 0.83-0.93; P = 3.0 × 10-6), with genome-wide significance after meta-analysis with 7 replication cohorts totaling 312 118 individuals (9395 cases of AS) (OR, 0.91; 95% CI, 0.88-0.94; P = 2.5 × 10-8). A consistent association with AVC was also observed (OR, 0.91; 95% CI, 0.83-0.99; P = .03). A higher ratio of arachidonic acid to linoleic acid was associated with AVC (OR per SD of the natural logarithm, 1.19; 95% CI, 1.09-1.30; P = 6.6 × 10-5). In mendelian randomization, increased FADS1 liver expression and arachidonic acid were associated with AS (OR per unit of normalized expression, 1.31 [95% CI, 1.17-1.48; P = 7.4 × 10-6]; OR per 5-percentage point increase in arachidonic acid for AVC, 1.23 [95% CI, 1.01-1.49; P = .04]; OR per 5-percentage point increase in arachidonic acid for AS, 1.08 [95% CI, 1.04-1.13; P = 4.1 × 10-4]). Conclusions and Relevance: Variation at the FADS1/2 locus was associated with AS and AVC. Findings from biomarker measurements and mendelian randomization appear to link ω-6 fatty acid biosynthesis to AS, which may represent a therapeutic target