Search

It has been suggested that diagnostic yield (DY) from Exome Sequencing (ES) may be lower among patients with non-European ancestries than those with European ancestry. We examined the association of DY with estimated continental/subcontinental genetic ancestry in a racially/ethnically diverse pediatric and prenatal clinical cohort. Cases (N = 845) with suspected genetic disorders underwent ES for diagnosis. Continental/subcontinental genetic ancestry proportions were estimated from the ES data. We compared the distribution of genetic ancestries in positive, negative, and inconclusive cases by Kolmogorov-Smirnov tests and linear associations of ancestry with DY by Cochran-Armitage trend tests. We observed no reduction in overall DY associated with any genetic ancestry (African, Native American, East Asian, European, Middle Eastern, South Asian). However, we observed a relative increase in proportion of autosomal recessive homozygous inheritance versus other inheritance patterns associated with Middle Eastern and South Asian ancestry, due to consanguinity. In this empirical study of ES for undiagnosed pediatric and prenatal genetic conditions, genetic ancestry was not associated with the likelihood of a positive diagnosis, supporting the equitable use of ES in diagnosis of previously undiagnosed but potentially Mendelian disorders across all ancestral populations.

Narcolepsy type 1 (NT1) is caused by a loss of hypocretin/orexin transmission. Risk factors include pandemic 2009 H1N1 influenza A infection and immunization with Pandemrix®. Here, we dissect disease mechanisms and interactions with environmental triggers in a multi-ethnic sample of 6,073 cases and 84,856 controls. We fine-mapped GWAS signals within HLA (DQ0602, DQB1*03:01 and DPB1*04:02) and discovered seven novel associations (CD207, NAB1, IKZF4-ERBB3, CTSC, DENND1B, SIRPG, PRF1). Significant signals at TRA and DQB1*06:02 loci were found in 245 vaccination-related cases, who also shared polygenic risk. T cell receptor associations in NT1 modulated TRAJ*24, TRAJ*28 and TRBV*4-2 chain-usage. Partitioned heritability and immune cell enrichment analyses found genetic signals to be driven by dendritic and helper T cells. Lastly comorbidity analysis using data from FinnGen, suggests shared effects between NT1 and other autoimmune diseases. NT1 genetic variants shape autoimmunity and response to environmental triggers, including influenza A infection and immunization with Pandemrix®.

The diagnostic yield of exome sequencing (ES) has primarily been evaluated in individuals of European ancestry, with less focus on underrepresented minority (URM) and underserved (US) patients. We evaluated the diagnostic yield of ES in a cohort of predominantly US and URM pediatric and prenatal patients suspected to have a genetic disorder. Eligible pediatric patients had multiple congenital anomalies and/or neurocognitive disabilities and prenatal patients had one or more structural anomalies, disorders of fetal growth, or fetal effusions. URM and US patients were prioritized for enrollment and underwent ES at a single academic center. We identified definitive positive or probable positive results in 201/845 (23.8%) patients, with a significantly higher diagnostic rate in pediatric (26.7%) compared to prenatal patients (19.0%) (P = 0.01). For both pediatric and prenatal patients, the diagnostic yield and frequency of inconclusive findings did not differ significantly between URM and non-URM patients or between patients with US status and those without US status. Our results demonstrate a similar diagnostic yield of ES between prenatal and pediatric URM/US patients and non-URM/US patients for positive and inconclusive results. These data support the use of ES to identify clinically relevant variants in patients from diverse populations.

BackgroundPrior studies of the glycemic effect of statins have been inconsistent. Also, most studies have only considered a short duration of statin use; the effect of long-term statin use on fasting glucose (FG) has not been well examined. The aim of this work is to investigate the effect of long-term statin exposure on FG levels.MethodsUsing electronic health record (EHR) data from a large and diverse longitudinal cohort, we defined long-term statin exposure in two ways: the cumulative years of statin use (cumulative supply) and the years' supply-weighted sum of doses (cumulative dose). Simvastatin, lovastatin, atorvastatin and pravastatin were included in the analysis. The relationship between statin exposure and FG was examined using linear regression with mixed effects modeling, comparing statin users before and after initiating statins and statin never-users.ResultsWe examined 593,130 FG measurements from 87,151 individuals over a median follow up of 20 years. Of these, 42,678 were never-users and 44,473 were statin users with a total of 730,031 statin prescriptions. FG was positively associated with cumulative supply of statin but not comulative dose when both measures were in the same model. While statistically significant, the annual increase in FG attributable to statin exposure was modest at only 0.14 mg/dl, with only slight and non-significant differences among statin types.ConclusionsElevation in FG level is associated with statin exposure, but the effect is modest. The results suggest that the risk of a clinically significant increase in FG attributable to long-term statin use is small for most individuals.

Genetic substudies of randomized controlled trials demonstrate that high coronary heart disease (CHD) polygenic risk score modifies statin CHD relative risk reduction; it is unknown if the association extends to statin users undergoing routine care. We sought to determine how statin effectiveness is modified by CHD polygenic risk score in a real-world cohort of participants without previous myocardial infarction. We determined CHD polygenic risk scores in participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. Covariate-adjusted Cox regression models were used to compare the risk of cardiovascular outcomes between statin users and matched nonusers. Statin effectiveness on incident myocardial infarction showed no gradient with increasing 10-year Pooled Cohort Equations atherosclerotic cardiovascular disease (ASCVD) risk across low, borderline, intermediate, and high ASCVD risk score groups. In contrast, statin effectiveness by polygenic risk was largest in the high polygenic risk score group (hazard ratio (HR) 0.41, 95% confidence interval (CI), 0.31-0.53; P = 1.5E-11), intermediate in the intermediate polygenic risk score group (HR 0.56, 95% CI, 0.47-0.66; P = 8.4E-12), and smallest in the low polygenic risk score group (HR 0.67, 95% CI, 0.47-0.97; P = 0.03; P for high vs. low = 0.01). ASCVD risk and statin low-density lipoprotein cholesterol (LDL-C) lowering did not differ across polygenic risk score groups. In patients undergoing routine care, CHD polygenic risk modified statin relative risk reduction of incident myocardial infarction independent of LDL-C lowering. Our findings extend prior work by identifying a subset (i.e., self-identified White individuals with low CHD polygenic risk scores) with attenuated clinical benefit from statins.

We report results from the Bipolar Exome (BipEx) collaboration analysis of whole-exome sequencing of 13,933 patients with bipolar disorder (BD) matched with 14,422 controls. We find an excess of ultra-rare protein-truncating variants (PTVs) in patients with BD among genes under strong evolutionary constraint in both major BD subtypes. We find enrichment of ultra-rare PTVs within genes implicated from a recent schizophrenia exome meta-analysis (SCHEMA; 24,248 cases and 97,322 controls) and among binding targets of CHD8. Genes implicated from genome-wide association studies (GWASs) of BD, however, are not significantly enriched for ultra-rare PTVs. Combining gene-level results with SCHEMA, AKAP11 emerges as a definitive risk gene (odds ratio (OR) = 7.06, P = 2.83 × 10-9). At the protein level, AKAP-11 interacts with GSK3B, the hypothesized target of lithium, a primary treatment for BD. Our results lend support to BD's polygenicity, demonstrating a role for rare coding variation as a significant risk factor in BD etiology.

The association between the c.521T>C variant allele in SLCO1B1 (reference single nucleotide polymorphism (rs)4149056) and simvastatin-induced myotoxicity was discovered over a decade ago; however, whether this relationship represents a class effect is still not fully known. The aim of this study was to investigate the relationship between rs4149056 genotype and statin-induced myotoxicity in patients taking atorvastatin and lovastatin. Study participants were from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. A total of 233 statin-induced myopathy + rhabdomyolysis cases met the criteria for inclusion and were matched to 2,342 controls. To validate the drug response phenotype, we replicated the previously established association between rs4149056 genotype and simvastatin-induced myotoxicity. In particular, compared with homozygous T allele carriers, there was a significantly increased risk of simvastatin-induced myopathy + rhabdomyolysis in homozygous carriers of the C allele (CC vs. TT, odds ratio [OR] 4.62, 95% confidence interval [CI] 1.58-11.90, P = 0.003). For lovastatin users, homozygous carriers of the C allele were also at increased risk of statin-induced myopathy + rhabdomyolysis (CC vs. TT, OR 4.49, 95% CI 1.68-10.80, P = 0.001). In atorvastatin users, homozygous carriers of the C allele were twice as likely to experience statin-induced myopathy, though this association did not achieve statistical significance (CC vs. TT, OR 2.00, 95% CI 0.44-6.59, P = 0.30). In summary, our findings suggest that the association of rs4149056 with simvastatin-related myotoxicity may also extend to lovastatin. More data is needed to determine the extent of the association in atorvastatin users. Altogether, these data expand the evidence base for informing guidelines of pharmacogenetic-based statin prescribing practices.

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

Plasma lipids are known heritable risk factors for cardiovascular disease, but increasing evidence also supports shared genetics with diseases of other organ systems. We devised a comprehensive three-phase framework to identify new lipid-associated genes and study the relationships among lipids, genotypes, gene expression and hundreds of complex human diseases from the Electronic Medical Records and Genomics (347 traits) and the UK Biobank (549 traits). Aside from 67 new lipid-associated genes with strong replication, we found evidence for pleiotropic SNPs/genes between lipids and diseases across the phenome. These include discordant pleiotropy in the HLA region between lipids and multiple sclerosis and putative causal paths between triglycerides and gout, among several others. Our findings give insights into the genetic basis of the relationship between plasma lipids and diseases on a phenome-wide scale and can provide context for future prevention and treatment strategies.

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.

Public health newborn screening (NBS) programs provide population-scale ascertainment of rare, treatable conditions that require urgent intervention. Tandem mass spectrometry (MS/MS) is currently used to screen newborns for a panel of rare inborn errors of metabolism (IEMs)1-4. The NBSeq project evaluated whole-exome sequencing (WES) as an innovative methodology for NBS. We obtained archived residual dried blood spots and data for nearly all IEM cases from the 4.5 million infants born in California between mid-2005 and 2013 and from some infants who screened positive by MS/MS, but were unaffected upon follow-up testing. WES had an overall sensitivity of 88% and specificity of 98.4%, compared to 99.0% and 99.8%, respectively for MS/MS, although effectiveness varied among individual IEMs. Thus, WES alone was insufficiently sensitive or specific to be a primary screen for most NBS IEMs. However, as a secondary test for infants with abnormal MS/MS screens, WES could reduce false-positive results, facilitate timely case resolution and in some instances even suggest more appropriate or specific diagnosis than that initially obtained. This study represents the largest, to date, sequencing effort of an entire population of IEM-affected cases, allowing unbiased assessment of current capabilities of WES as a tool for population screening.

Hundreds of loci have been associated with blood pressure (BP) traits from many genome-wide association studies. We identified an enrichment of these loci in aorta and tibial artery expression quantitative trait loci in our previous work in ~100 000 Genetic Epidemiology Research on Aging study participants. In the present study, we sought to fine-map known loci and identify novel genes by determining putative regulatory regions for these and other tissues relevant to BP. We constructed maps of putative cis-regulatory elements (CREs) using publicly available open chromatin data for the heart, aorta and tibial arteries, and multiple kidney cell types. Variants within these regions may be evaluated quantitatively for their tissue- or cell-type-specific regulatory impact using deltaSVM functional scores, as described in our previous work. We aggregate variants within these putative CREs within 50 Kb of the start or end of 'expressed' genes in these tissues or cell types using public expression data and use deltaSVM scores as weights in the group-wise sequence kernel association test to identify candidates. We test for association with both BP traits and expression within these tissues or cell types of interest and identify the candidates MTHFR, C10orf32, CSK, NOV, ULK4, SDCCAG8, SCAMP5, RPP25, HDGFRP3, VPS37B and PPCDC. Additionally, we examined two known QT interval genes, SCN5A and NOS1AP, in the Atherosclerosis Risk in Communities Study, as a positive control, and observed the expected heart-specific effect. Thus, our method identifies variants and genes for further functional testing using tissue- or cell-type-specific putative regulatory information.

We describe an 11-year old boy with severe global developmental delays, failure to thrive and growth retardation, refractory seizures with recurrent status epilepticus, hypogammaglobulinemia, hypergonadotropic hypogonadism, and duodenal strictures. He had facial and skin findings compatible with trichothiodystrophy, including sparse and brittle hair, thin eyebrows, and dry skin. Exome sequencing showed a hemizygous, truncating variant in RNF113A, c.903_910delGCAGACCA, predicting p.(Gln302fs*12), that was inherited from his mother. Although his clinical features overlap closely with features described in the two previously reported male first cousins with RNF113A loss of function mutations, the duodenal strictures seen in this patient have not been reported. Interestingly, the patient's mother had short stature and 100% skewed X-inactivation as seen in other obligate female carriers. A second male with developmental delays, microcephaly, seizures, ambiguous genitalia, and facial anomalies that included sparse and brittle hair, thin eyebrows and dry skin was recently reported to have c.897_898delTG, predicting p.(Cys299*) in RNF113A and we provide additional clinical details for this patient. This report further supports deleterious variants in RNF113A as a cause of a novel trichothiodystrophy syndrome.

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.

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.

Mammographic density (MD) phenotypes are strongly associated with breast cancer risk and highly heritable. In this GWAS meta-analysis of 24,192 women, we identify 31 MD loci at P

ImportancePrimary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders.ObjectivesTo perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma.Design, settings, and participantsA 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma.ExposuresGenetic variants associated with primary open-angle glaucoma.Main outcomes and measuresPresence of primary open-angle glaucoma. Genome-wide significance was defined as P C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P

Importance:Primary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders. Objectives:To perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma. Design, Settings, and Participants:A 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma. Exposures:Genetic variants associated with primary open-angle glaucoma. Main Outcomes and Measures:Presence of primary open-angle glaucoma. Genome-wide significance was defined as P C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P

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.

Body mass index (BMI), a proxy measure for obesity, is determined by both environmental (including ethnicity, age, and sex) and genetic factors, with > 400 BMI-associated loci identified to date. However, the impact, interplay, and underlying biological mechanisms among BMI, environment, genetics, and ancestry are not completely understood. To further examine these relationships, we utilized 427,509 calendar year-averaged BMI measurements from 100,418 adults from the single large multiethnic Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. We observed substantial independent ancestry and nationality differences, including ancestry principal component interactions and nonlinear effects. To increase the list of BMI-associated variants before assessing other differences, we conducted a genome-wide association study (GWAS) in GERA, with replication in the Genetic Investigation of Anthropomorphic Traits (GIANT) consortium combined with the UK Biobank (UKB), followed by GWAS in GERA combined with GIANT, with replication in the UKB. We discovered 30 novel independent BMI loci (P < 5.0 × 10-8) that replicated. We then assessed the proportion of BMI variance explained by sex in the UKB using previously identified loci compared to previously and newly identified loci and found slight increases: from 3.0 to 3.3% for males and from 2.7 to 3.0% for females. Further, the variance explained by previously and newly identified variants decreased with increasing age in the GERA and UKB cohorts, echoed in the variance explained by the entire genome, which also showed gene-age interaction effects. Finally, we conducted a tissue expression QTL enrichment analysis, which revealed that GWAS BMI-associated variants were enriched in the cerebellum, consistent with prior work in humans and mice.

BackgroundLow-density lipoprotein cholesterol (LDL-C) response to statin therapy has not been fully elucidated in real-world populations. The primary objective of this study was to characterize statin LDL-C dose-response and its heritability in a large, multiethnic population of statin users.MethodsWe determined the effect of statin dosing on lipid measures utilizing electronic health records in 33 139 statin users from the Kaiser Permanente GERA cohort (Genetic Epidemiology Research on Adult Health and Aging). The relationship between statin defined daily dose and lipid parameter response (percent change) was determined.ResultsDefined daily dose and LDL-C response was associated in a log-linear relationship (β, -6.17; SE, 0.09; P

The Clinical Sequencing Evidence-Generating Research (CSER) consortium, now in its second funding cycle, is investigating the effectiveness of integrating genomic (exome or genome) sequencing into the clinical care of diverse and medically underserved individuals in a variety of healthcare settings and disease states. The consortium comprises a coordinating center, six funded extramural clinical projects, and an ongoing National Human Genome Research Institute (NHGRI) intramural project. Collectively, these projects aim to enroll and sequence over 6,100 participants in four years. At least 60% of participants will be of non-European ancestry or from underserved settings, with the goal of diversifying the populations that are providing an evidence base for genomic medicine. Five of the six clinical projects are enrolling pediatric patients with various phenotypes. One of these five projects is also enrolling couples whose fetus has a structural anomaly, and the sixth project is enrolling adults at risk for hereditary cancer. The ongoing NHGRI intramural project has enrolled primarily healthy adults. Goals of the consortium include assessing the clinical utility of genomic sequencing, exploring medical follow up and cascade testing of relatives, and evaluating patient-provider-laboratory level interactions that influence the use of this technology. The findings from the CSER consortium will offer patients, healthcare systems, and policymakers a clearer understanding of the opportunities and challenges of providing genomic medicine in diverse populations and settings, and contribute evidence toward developing best practices for the delivery of clinically useful and cost-effective genomic sequencing in diverse healthcare settings.

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

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.

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

Hypertension is a leading cause of global disease, mortality, and disability. While individuals of African descent suffer a disproportionate burden of hypertension and its complications, they have been underrepresented in genetic studies. To identify novel susceptibility loci for blood pressure and hypertension in people of African ancestry, we performed both single and multiple-trait genome-wide association analyses. We analyzed 21 genome-wide association studies comprised of 31,968 individuals of African ancestry, and validated our results with additional 54,395 individuals from multi-ethnic studies. These analyses identified nine loci with eleven independent variants which reached genome-wide significance (P < 1.25×10-8) for either systolic and diastolic blood pressure, hypertension, or for combined traits. Single-trait analyses identified two loci (TARID/TCF21 and LLPH/TMBIM4) and multiple-trait analyses identified one novel locus (FRMD3) for blood pressure. At these three loci, as well as at GRP20/CDH17, associated variants had alleles common only in African-ancestry populations. Functional annotation showed enrichment for genes expressed in immune and kidney cells, as well as in heart and vascular cells/tissues. Experiments driven by these findings and using angiotensin-II induced hypertension in mice showed altered kidney mRNA expression of six genes, suggesting their potential role in hypertension. Our study provides new evidence for genes related to hypertension susceptibility, and the need to study African-ancestry populations in order to identify biologic factors contributing to hypertension.

The rapid development of genomic sequencing technologies has decreased the cost of genetic analysis to the extent that it seems plausible that genome-scale sequencing could have widespread availability in pediatric care. Genomic sequencing provides a powerful diagnostic modality for patients who manifest symptoms of monogenic disease and an opportunity to detect health conditions before their development. However, many technical, clinical, ethical, and societal challenges should be addressed before such technology is widely deployed in pediatric practice. This article provides an overview of the Newborn Sequencing in Genomic Medicine and Public Health Consortium, which is investigating the application of genome-scale sequencing in newborns for both diagnosis and screening.

Prostate-specific antigen (PSA) levels have been used for detection and surveillance of prostate cancer (PCa). However, factors other than PCa-such as genetics-can impact PSA. Here we present findings from a genome-wide association study (GWAS) of PSA in 28,503 Kaiser Permanente whites and 17,428 men from replication cohorts. We detect 40 genome-wide significant (P

Longitudinal electronic health records on 99,785 Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort individuals provided 1,342,814 systolic and diastolic blood pressure measurements for a genome-wide association study on long-term average systolic, diastolic, and pulse pressure. We identified 39 new loci among 75 genome-wide significant loci (P ≤ 5 × 10-8), with most replicating in the combined International Consortium for Blood Pressure (ICBP; n = 69,396) and UK Biobank (UKB; n = 152,081) studies. Combining GERA with ICBP yielded 36 additional new loci, with most replicating in UKB. Combining all three studies (n = 321,262) yielded 241 additional genome-wide significant loci, although no replication sample was available for these. All associated loci explained 2.9%, 2.5%, and 3.1% of variation in systolic, diastolic, and pulse pressure, respectively, in GERA non-Hispanic whites. Using multiple blood pressure measurements in GERA doubled the variance explained. A normalized risk score was associated with time to onset of hypertension (hazards ratio = 1.18, P = 8.2 × 10-45). Expression quantitative trait locus analysis of blood pressure loci showed enrichment in aorta and tibial artery.

PurposeTo evaluate the association between refractive error and the prevalence of glaucoma by race or ethnicity.DesignCross-sectional study.ParticipantsKaiser Permanente Northern California Health Plan members with refractive error measured at 35 years of age or older between 2008 and 2014 and with no history of cataract surgery, refractive surgery, or a corneal disorder.MethodsWe identified 34 040 members with glaucoma or ocular hypertension (OHTN; cases) and 403 398 members without glaucoma (controls). Glaucoma cases were classified as primary angle-closure glaucoma (PACG); 1 of the 4 forms of open-angle glaucoma: primary open-angle glaucoma (POAG), normal-tension glaucoma (NTG), pigmentary glaucoma (PIGM), and pseudoexfoliation glaucoma (PEX); or OHTN. Refractive error, expressed as spherical equivalent (SE), was coded as a continuous trait and also as categories. Logistic regression analyses were used to estimate the association between refractive error and the prevalence of glaucoma overall and in specific racial or ethnic groups.Main outcome measuresThe association between refractive error and glaucoma subtypes evaluated as odds ratios (ORs) with 95% confidence intervals (CIs).ResultsIn controls, the mean SE was -0.59 diopters (D) (standard deviation, 2.62 D). Each 1-D reduction in SE was associated with a 22% decrease in the odds of PACG (OR, 0.78; 95% CI, 0.77-0.80) and with increases in the odds of open-angle glaucoma ranging from 1.23 (95% CI, 1.20-1.26) for PIGM, to 1.07 (95% CI, 1.03-1.11) for PEX, and to 1.05 (95% CI, 1.04-1.06) for OHTN. In addition, we observed a stronger association between myopia and POAG among non-Hispanic whites (OR, 1.12; 95% CI, 1.11-1.13) and NTG among Asians (OR, 1.17; 95% CI, 1.15-1.20) and non-Hispanic whites (OR, 1.19; 95% CI, 1.15-1.22).ConclusionsMyopia was associated with an increased prevalence of all forms of open-angle glaucoma and OHTN, whereas hyperopia was associated with a substantially increased prevalence of PACG. Although high myopia is a strong risk factor for glaucoma subtypes, low and moderate myopia also have a significant effect on glaucoma risk. Additionally, there were moderate racial differences in the association of myopia with the risk of POAG and NTG.

UnlabelledA genome-wide association study (GWAS) of prostate cancer in Kaiser Permanente health plan members (7,783 cases, 38,595 controls; 80.3% non-Hispanic white, 4.9% African-American, 7.0% East Asian, and 7.8% Latino) revealed a new independent risk indel rs4646284 at the previously identified locus 6q25.3 that replicated in PEGASUS (N = 7,539) and the Multiethnic Cohort (N = 4,679) with an overall P = 1.0 × 10(-19) (OR, 1.18). Across the 6q25.3 locus, rs4646284 exhibited the strongest association with expression of SLC22A1 (P = 1.3 × 10(-23)) and SLC22A3 (P = 3.2 × 10(-52)). At the known 19q13.33 locus, rs2659124 (P = 1.3 × 10(-13); OR, 1.18) nominally replicated in PEGASUS. A risk score of 105 known risk SNPs was strongly associated with prostate cancer (P < 1.0 × 10(-8)). Comparing the highest to lowest risk score deciles, the OR was 6.22 for non-Hispanic whites, 5.82 for Latinos, 3.77 for African-Americans, and 3.38 for East Asians. In non-Hispanic whites, the 105 risk SNPs explained approximately 7.6% of disease heritability. The entire GWAS array explained approximately 33.4% of heritability, with a 4.3-fold enrichment within DNaseI hypersensitivity sites (P = 0.004).SignificanceTaken together, our findings of independent risk variants, ethnic variation in existing SNP replication, and remaining unexplained heritability have important implications for further clarifying the genetic risk of prostate cancer. Our findings also suggest that there may be much promise in evaluating understudied variation, such as indels and ethnically diverse populations.

The Kaiser Permanente (KP) Research Program on Genes, Environment and Health (RPGEH), in collaboration with the University of California-San Francisco, undertook genome-wide genotyping of >100,000 subjects that constitute the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. The project, which generated >70 billion genotypes, represents the first large-scale use of the Affymetrix Axiom Genotyping Solution. Because genotyping took place over a short 14-month period, creating a near-real-time analysis pipeline for experimental assay quality control and final optimized analyses was critical. Because of the multi-ethnic nature of the cohort, four different ethnic-specific arrays were employed to enhance genome-wide coverage. All assays were performed on DNA extracted from saliva samples. To improve sample call rates and significantly increase genotype concordance, we partitioned the cohort into disjoint packages of plates with similar assay contexts. Using strict QC criteria, the overall genotyping success rate was 103,067 of 109,837 samples assayed (93.8%), with a range of 92.1-95.4% for the four different arrays. Similarly, the SNP genotyping success rate ranged from 98.1 to 99.4% across the four arrays, the variation depending mostly on how many SNPs were included as single copy vs. double copy on a particular array. The high quality and large scale of genotype data created on this cohort, in conjunction with comprehensive longitudinal data from the KP electronic health records of participants, will enable a broad range of highly powered genome-wide association studies on a diversity of traits and conditions.

Using genome-wide genotypes, we characterized the genetic structure of 103,006 participants in the Kaiser Permanente Northern California multi-ethnic Genetic Epidemiology Research on Adult Health and Aging Cohort and analyzed the relationship to self-reported race/ethnicity. Participants endorsed any of 23 race/ethnicity/nationality categories, which were collapsed into seven major race/ethnicity groups. By self-report the cohort is 80.8% white and 19.2% minority; 93.8% endorsed a single race/ethnicity group, while 6.2% endorsed two or more. Principal component (PC) and admixture analyses were generally consistent with prior studies. Approximately 17% of subjects had genetic ancestry from more than one continent, and 12% were genetically admixed, considering only nonadjacent geographical origins. Self-reported whites were spread on a continuum along the first two PCs, indicating extensive mixing among European nationalities. Self-identified East Asian nationalities correlated with genetic clustering, consistent with extensive endogamy. Individuals of mixed East Asian-European genetic ancestry were easily identified; we also observed a modest amount of European genetic ancestry in individuals self-identified as Filipinos. Self-reported African Americans and Latinos showed extensive European and African genetic ancestry, and Native American genetic ancestry for the latter. Among 3741 genetically identified parent-child pairs, 93% were concordant for self-reported race/ethnicity; among 2018 genetically identified full-sib pairs, 96% were concordant; the lower rate for parent-child pairs was largely due to intermarriage. The parent-child pairs revealed a trend toward increasing exogamy over time; the presence in the cohort of individuals endorsing multiple race/ethnicity categories creates interesting challenges and future opportunities for genetic epidemiologic studies.

The Kaiser Permanente Research Program on Genes, Environment, and Health (RPGEH) Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort includes DNA specimens extracted from saliva samples of 110,266 individuals. Because of its relationship to aging, telomere length measurement was considered an important biomarker to develop on these subjects. To assay relative telomere length (TL) on this large cohort over a short time period, we created a novel high throughput robotic system for TL analysis and informatics. Samples were run in triplicate, along with control samples, in a randomized design. As part of quality control, we determined the within-sample variability and employed thresholds for the elimination of outlying measurements. Of 106,902 samples assayed, 105,539 (98.7%) passed all quality control (QC) measures. As expected, TL in general showed a decline with age and a sex difference. While telomeres showed a negative correlation with age up to 75 years, in those older than 75 years, age positively correlated with longer telomeres, indicative of an association of longer telomeres with more years of survival in those older than 75. Furthermore, while females in general had longer telomeres than males, this difference was significant only for those older than age 50. An additional novel finding was that the variance of TL between individuals increased with age. This study establishes reliable assay and analysis methodologies for measurement of TL in large, population-based human studies. The GERA cohort represents the largest currently available such resource, linked to comprehensive electronic health and genotype data for analysis.

PurposeWe compared across age-related macular degeneration (AMD) subtypes the effect of AMD risk variants, their predictive power, and heritability.MethodsThe prevalence of AMD was estimated among active non-Hispanic white Kaiser Permanente Northern California members who were at least 65 years of age as of June 2013. The genetic analysis included 5,170 overall AMD cases ascertained from electronic health records (EHR), including 1,239 choroidal neovascularization (CNV) cases and 1,060 nonexudative AMD cases without CNV, and 23,130 controls of non-Hispanic white ancestry from the Kaiser Permanente Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. Imputation was based on the 1000 Genomes Project reference panel.ResultsThe narrow-sense heritability due to common autosomal single nucleotide polymorphisms (SNPs) was 0.37 for overall AMD, 0.19 for AMD unspecified, 0.20 for nonexudative AMD, and 0.60 for CNV. For the 19 previously reported AMD risk loci, the area under the receiver operating characteristic (ROC) curve was 0.675 for overall AMD, 0.640 for AMD unspecified, 0.678 for nonexudative AMD, and 0.766 for CNV. The individual effects on the risk of AMD for 18 of the 19 SNPs were in a consistent direction with those previously reported, including a protective effect of the APOE ε4 allele. Conversely, the risk of AMD was significantly increased in carriers of the ε2 allele.ConclusionsThese findings provide an independent confirmation of many of the previously identified AMD risk loci, and support a potentially greater role of genetic factors in the development of CNV. The replication of established associations validates the use of EHR in genetic studies of ophthalmologic traits.

BackgroundAdmixture mapping of lipids was followed-up by family-based association analysis to identify variants for cardiovascular disease in African Americans.Methods and resultsThe present study conducted admixture mapping analysis for total cholesterol, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and triglycerides. The analysis was performed in 1905 unrelated African American subjects from the National Heart, Lung and Blood Institute's Family Blood Pressure Program (FBPP). Regions showing admixture evidence were followed-up with family-based association analysis in 3556 African American subjects from the FBPP. The admixture mapping and family-based association analyses were adjusted for age, age(2), sex, body mass index, and genome-wide mean ancestry to minimize the confounding caused by population stratification. Regions that were suggestive of local ancestry association evidence were found on chromosomes 7 (low-density lipoprotein cholesterol), 8 (high-density lipoprotein cholesterol), 14 (triglycerides), and 19 (total cholesterol and triglycerides). In the fine-mapping analysis, 52 939 single-nucleotide polymorphisms (SNPs) were tested and 11 SNPs (8 independent SNPs) showed nominal significant association with high-density lipoprotein cholesterol (2 SNPs), low-density lipoprotein cholesterol (4 SNPs), and triglycerides (5 SNPs). The family data were used in the fine-mapping to identify SNPs that showed novel associations with lipids and regions, including genes with known associations for cardiovascular disease.ConclusionsThis study identified regions on chromosomes 7, 8, 14, and 19 and 11 SNPs from the fine-mapping analysis that were associated with high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and triglycerides for further studies of cardiovascular disease in African Americans.

© 2015 Hoffmann et al.An efficient approach to characterizing the disease burden of rare genetic variants is to impute them into large well-phenotyped cohorts with existing genome-wide genotype data using large sequenced referenced panels. The success of t

We report a genome-wide association study (GWAS) and admixture analysis of glaucoma in 12 008 African-American and Hispanic women (age 50-79 years) from the Women's Health Initiative (WHI). Although GWAS of glaucoma have been conducted on several populations, this is the first to look at glaucoma in individuals of African-American and Hispanic race/ethnicity. Prevalent and incident glaucoma was determined by self-report from study questionnaires administered at baseline (1993-1998) and annually through 2005. For African Americans, there was a total of 658 prevalent cases, 1062 incident cases and 6067 individuals who never progressed to glaucoma. For our replication cohort, we used the WHI Hispanics, including 153 prevalent cases, 336 incident cases and 2685 non-cases. We found an association of African ancestry with glaucoma incidence in African Americans (hazards ratio 1.62, 95% CI 1.023-2.56, P = 0.038) and in Hispanics (hazards ratio 3.21, 95% CI 1.32-7.80, P = 0.011). Although we found that no previously identified glaucoma SNPs replicated in either the WHI African Americans or Hispanics, a risk score combining all previously reported hits was significant in African-American prevalent cases (P = 0.0046), and was in the expected direction in the incident cases, as well as in the Hispanic incident cases. Additionally, after imputing to 1000 Genomes, two less common independent SNPs were suggestive in African Americans, but had too low of an allele frequency in Hispanics to test for replication. These results suggest the possibility of a distinct genetic architecture underlying glaucoma in individuals of African ancestry.

Exome and whole-genome sequencing studies are becoming increasingly common, but little is known about the accuracy of the genotype calls made by the commonly used platforms. Here we use replicate high-coverage sequencing of blood and saliva DNA samples from four European-American individuals to estimate lower bounds on the error rates of Complete Genomics and Illumina HiSeq whole-genome and whole-exome sequencing. Error rates for nonreference genotype calls range from 0.1% to 0.6%, depending on the platform and the depth of coverage. Additionally, we found (1) no difference in the error profiles or rates between blood and saliva samples; (2) Complete Genomics sequences had substantially higher error rates than Illumina sequences had; (3) error rates were higher (up to 6%) for rare or unique variants; (4) error rates generally declined with genotype quality (GQ) score, but in a nonlinear fashion for the Illumina data, likely due to loss of specificity of GQ scores greater than 60; and (5) error rates increased with increasing depth of coverage for the Illumina data. These findings, especially (3)-(5), suggest that caution should be taken in interpreting the results of next-generation sequencing-based association studies, and even more so in clinical application of this technology in the absence of validation by other more robust sequencing or genotyping methods.