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.

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma, a leading cause of blindness. IOP heritability has been estimated to up to 67%, and to date only 11 IOP loci have been reported, accounting for 1.5% of IOP variability. Here, we conduct a genome-wide association study of IOP in 69,756 untreated individuals of European, Latino, Asian, and African ancestry. Multiple longitudinal IOP measurements were collected through electronic health records and, in total, 356,987 measurements were included. We identify 47 genome-wide significant IOP-associated loci (P < 5 × 10-8); of the 40 novel loci, 14 replicate at Bonferroni significance in an external genome-wide association study analysis of 37,930 individuals of European and Asian descent. We further examine their effect on the risk of glaucoma within our discovery sample. Using longitudinal IOP measurements from electronic health records improves our power to identify new variants, which together explain 3.7% of IOP variation.

We sought to estimate the causal effect of low serum 25(OH)D on multiple sclerosis (MS) susceptibility that is not confounded by environmental or lifestyle factors or subject to reverse causality.We conducted mendelian randomization (MR) analyses using an instrumental variable (IV) comprising 3 single nucleotide polymorphisms found to be associated with serum 25(OH)D levels at genome-wide significance. We analyzed the effect of the IV on MS risk and both age at onset and disease severity in 2 separate populations using logistic regression models that controlled for sex, year of birth, smoking, education, genetic ancestry, body mass index at age 18-20 years or in 20s, a weighted genetic risk score for 110 known MS-associated variants, and the presence of one or more HLA-DRB1*15:01 alleles.Findings from MR analyses using the IV showed increasing levels of 25(OH)D are associated with a decreased risk of MS in both populations. In white, non-Hispanic members of Kaiser Permanente Northern California (1,056 MS cases and 9,015 controls), the odds ratio (OR) was 0.79 (p = 0.04, 95% confidence interval (CI): 0.64-0.99). In members of a Swedish population from the Epidemiological Investigation of Multiple Sclerosis and Genes and Environment in Multiple Sclerosis MS case-control studies (6,335 cases and 5,762 controls), the OR was 0.86 (p = 0.03, 95% CI: 0.76-0.98). A meta-analysis of the 2 populations gave a combined OR of 0.85 (p = 0.003, 95% CI: 0.76-0.94). No association was observed for age at onset or disease severity.These results provide strong evidence that low serum 25(OH)D concentration is a cause of MS, independent of established risk factors.

Primary open-angle glaucoma (POAG) is a leading cause of irreversible vision loss, yet much of the genetic risk remains unaccounted for, especially in African-Americans who have a higher risk for developing POAG. We conduct a multiethnic genome-wide association study (GWAS) of POAG in the GERA cohort, with replication in the UK Biobank (UKB), and vice versa, GWAS in UKB with replication in GERA. We identify 24 loci (P < 5.0 × 10-8), including 14 novel, of which 9 replicate (near FMNL2, PDE7B, TMTC2, IKZF2, CADM2, DGKG, ANKH, EXOC2, and LMX1B). Functional studies support intraocular pressure-related influences of FMNL2 and LMX1B, with certain Lmx1b mutations causing high IOP and glaucoma resembling POAG in mice. The newly identified loci increase the proportion of variance explained in each GERA race/ethnicity group, with the largest gain in African-Americans (0.5-3.1%). A meta-analysis combining GERA and UKB identifies 24 additional loci. Our study provides important insights into glaucoma pathogenesis.

Multiple sclerosis (MS) is an autoimmune disease with high prevalence among populations of northern European ancestry. Past studies have shown that exposure to ultraviolet radiation could explain the difference in MS prevalence across the globe. In this study, we investigate whether the difference in MS prevalence could be explained by European genetic risk factors. We characterized the ancestry of MS-associated alleles using RFMix, a conditional random field parameterized by random forests, to estimate their local ancestry in the largest assembled admixed population to date, with 3,692 African Americans, 4,915 Asian Americans, and 3,777 Hispanics. The majority of MS-associated human leukocyte antigen (HLA) alleles, including the prominent HLA-DRB1*15:01 risk allele, exhibited cosmopolitan ancestry. Ancestry-specific MS-associated HLA alleles were also identified. Analysis of the HLA-DRB1*15:01 risk allele in African Americans revealed that alleles on the European haplotype conferred three times the disease risk compared to those on the African haplotype. Furthermore, we found evidence that the European and African HLA-DRB1*15:01 alleles exhibit single nucleotide polymorphism (SNP) differences in regions encoding the HLA-DRB1 antigen-binding heterodimer. Additional evidence for increased risk of MS conferred by the European haplotype were found for HLA-B*07:02 and HLA-A*03:01 in African Americans. Most of the 200 non-HLA MS SNPs previously established in European populations were not significantly associated with MS in admixed populations, nor were they ancestrally more European in cases compared to controls. Lastly, a genome-wide search of association between European ancestry and MS revealed a region of interest close to the ZNF596 gene on chromosome 8 in Hispanics; cases had a significantly higher proportion of European ancestry compared to controls. In conclusion, our study established that the genetic ancestry of MS-associated alleles is complex and implicated that difference in MS prevalence could be explained by the ancestry of MS-associated alleles.

Objective:Onset of multiple sclerosis (MS) occurs in childhood for approximately 5% of cases (pediatric MS, or ped-MS). Epigenetic influences are strongly implicated in MS pathogenesis in adults, including the contribution from microRNAs (miRNAs), small noncoding RNAs that affect gene expression by binding target gene mRNAs. Few studies have specifically examined miRNAs in ped-MS, but individuals developing MS at an early age may carry a relatively high burden of genetic risk factors, and miRNA dysregulation may therefore play a larger role in the development of ped-MS than in adult-onset MS. This study aimed to look for evidence of miRNA involvement in ped-MS pathogenesis. Methods:GWAS results from 486 ped-MS cases and 1362 controls from the U.S. Pediatric MS Network and Kaiser Permanente Northern California membership were investigated for miRNA-specific signals. First, enrichment of miRNA-target gene network signals was evaluated using MIGWAS software. Second, SNPs in miRNA genes and in target gene binding sites (miR-SNPs) were tested for association with ped-MS, and pathway analysis was performed on associated target genes. Results:MIGWAS analysis showed that miRNA-target gene signals were enriched in GWAS (P = 0.038) and identified 39 candidate biomarker miRNA-target gene pairs, including immune and neuronal signaling genes. The miR-SNP analysis implicated dysregulation of miRNA binding to target genes in five pathways, mainly involved in immune signaling. Interpretation:Evidence from GWAS suggests that miRNAs play a role in ped-MS pathogenesis by affecting immune signaling and other pathways. Candidate biomarker miRNA-target gene pairs should be further studied for diagnostic, prognostic, and/or therapeutic utility.

We investigated the association between 52 risk variants identified through genome-wide association studies and disease severity in multiple sclerosis (MS).Ten unique MS case data sets were analyzed. The Multiple Sclerosis Severity Score (MSSS) was calculated using the Expanded Disability Status Scale at study entry and disease duration. MSSS was considered as a continuous variable and as 2 dichotomous variables (median and extreme ends; MSSS of ≤5 vs >5 and MSSS of <2.5 vs ≥7.5, respectively). Single nucleotide polymorphisms (SNPs) were examined individually and as both combined weighted genetic risk score (wGRS) and unweighted genetic risk score (GRS) for association with disease severity. Random-effects meta-analyses were conducted and adjusted for cohort, sex, age at onset, and HLA-DRB1*15:01.A total of 7,125 MS cases were analyzed. The wGRS and GRS were not strongly associated with disease severity after accounting for cohort, sex, age at onset, and HLA-DRB1*15:01. After restricting analyses to cases with disease duration ≥10 years, associations were null (p value ≥0.05). No SNP was associated with disease severity after adjusting for multiple testing.The largest meta-analysis of established MS genetic risk variants and disease severity, to date, was performed. Results suggest that the investigated MS genetic risk variants are not associated with MSSS, even after controlling for potential confounders. Further research in large cohorts is needed to identify genetic determinants of disease severity using sensitive clinical and MRI measures, which are critical to understanding disease mechanisms and guiding development of effective treatments.

INTRODUCTION:Sex hormones have been implicated in the etiology of a number of diseases. To better understand disease etiology and the mechanisms of disease-risk factor associations, this analysis aimed to investigate the associations of anthropometric, sociodemographic and behavioural factors with a range of circulating sex hormones and sex hormone-binding globulin. METHODS:Statistical analyses of individual participant data from 12,330 male controls aged 25-85 years from 25 studies involved in the Endogenous Hormones Nutritional Biomarkers and Prostate Cancer Collaborative Group. Analysis of variance was used to estimate geometric means adjusted for study and relevant covariates. RESULTS:Older age was associated with higher concentrations of sex hormone-binding globulin and dihydrotestosterone and lower concentrations of dehydroepiandrosterone sulfate, free testosterone, androstenedione, androstanediol glucuronide and free estradiol. Higher body mass index was associated with higher concentrations of free estradiol, androstanediol glucuronide, estradiol and estrone and lower concentrations of dihydrotestosterone, testosterone, sex hormone-binding globulin, free testosterone, androstenedione and dehydroepiandrosterone sulfate. Taller height was associated with lower concentrations of androstenedione, testosterone, free testosterone and sex hormone-binding globulin and higher concentrations of androstanediol glucuronide. Current smoking was associated with higher concentrations of androstenedione, sex hormone-binding globulin and testosterone. Alcohol consumption was associated with higher concentrations of dehydroepiandrosterone sulfate, androstenedione and androstanediol glucuronide. East Asians had lower concentrations of androstanediol glucuronide and African Americans had higher concentrations of estrogens. Education and marital status were modestly associated with a small number of hormones. CONCLUSION:Circulating sex hormones in men are strongly associated with age and body mass index, and to a lesser extent with smoking status and alcohol consumption.

Insulin-like growth factors (IGFs) and insulin-like growth factor binding proteins (IGFBPs) have been implicated in the aetiology of several cancers. To better understand whether anthropometric, behavioural and sociodemographic factors may play a role in cancer risk via IGF signalling, we examined the cross-sectional associations of these exposures with circulating concentrations of IGFs (IGF-I and IGF-II) and IGFBPs (IGFBP-1, IGFBP-2 and IGFBP-3). The Endogenous Hormones, Nutritional Biomarkers and Prostate Cancer Collaborative Group dataset includes individual participant data from 16,024 male controls (i.e. without prostate cancer) aged 22-89 years from 22 prospective studies. Geometric means of protein concentrations were estimated using analysis of variance, adjusted for relevant covariates. Older age was associated with higher concentrations of IGFBP-1 and IGFBP-2 and lower concentrations of IGF-I, IGF-II and IGFBP-3. Higher body mass index was associated with lower concentrations of IGFBP-1 and IGFBP-2. Taller height was associated with higher concentrations of IGF-I and IGFBP-3 and lower concentrations of IGFBP-1. Smokers had higher concentrations of IGFBP-1 and IGFBP-2 and lower concentrations of IGFBP-3 than nonsmokers. Higher alcohol consumption was associated with higher concentrations of IGF-II and lower concentrations of IGF-I and IGFBP-2. African Americans had lower concentrations of IGF-II, IGFBP-1, IGFBP-2 and IGFBP-3 and Hispanics had lower IGF-I, IGF-II and IGFBP-3 than non-Hispanic whites. These findings indicate that a range of anthropometric, behavioural and sociodemographic factors are associated with circulating concentrations of IGFs and IGFBPs in men, which will lead to a greater understanding of the mechanisms through which these factors influence cancer risk.

Using the ImmunoChip custom genotyping array, we analyzed 14,498 subjects with multiple sclerosis and 24,091 healthy controls for 161,311 autosomal variants and identified 135 potentially associated regions (P < 1.0 × 10(-4)). In a replication phase, we combined these data with previous genome-wide association study (GWAS) data from an independent 14,802 subjects with multiple sclerosis and 26,703 healthy controls. In these 80,094 individuals of European ancestry, we identified 48 new susceptibility variants (P < 5.0 × 10(-8)), 3 of which we found after conditioning on previously identified variants. Thus, there are now 110 established multiple sclerosis risk variants at 103 discrete loci outside of the major histocompatibility complex. With high-resolution Bayesian fine mapping, we identified five regions where one variant accounted for more than 50% of the posterior probability of association. This study enhances the catalog of multiple sclerosis risk variants and illustrates the value of fine mapping in the resolution of GWAS signals.