Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

The genomics of heart failure: design and rationale of the HERMES consortium

Aims

The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure.

Methods and results

The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome‐wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow‐up following heart failure diagnosis ranged from 2 to 116 months. Forty‐nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34–90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low‐frequency variants (allele frequency 0.01–0.05) at P < 5 × 10⁻⁸ under an additive genetic model.

Conclusions

HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction.

Pleiotropic associations of heterozygosity for the SERPINA1 Z allele in the UK Biobank

Homozygosity for the SERPINA1 Z allele causes α1-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (β=1.02 cm, p=3.91×10-68), and with other nonrespiratory phenotypes including increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher height-adjusted forced expiratory volume in 1 s (FEV1) (β=19.36 mL, p=9.21×10-4) and FEV1/forced vital capacity (β=0.0031, p=1.22×10-5) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.

Meta-analysis uncovers genome-wide significant variants for rapid kidney function decline

Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m2/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m2 at follow-up among those with eGFRcrea 60 mL/min/1.73m2 or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or LARP4B. Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.

A Noncoding Variant Near PPP1R3B Promotes Liver Glycogen Storage and MetS, but Protects Against Myocardial Infarction

CONTEXT

Glycogen storage diseases are rare. Increased glycogen in the liver results in increased attenuation.

OBJECTIVE

Investigate the association and function of a noncoding region associated with liver attenuation but not histologic nonalcoholic fatty liver disease.

DESIGN

Genetics of Obesity-associated Liver Disease Consortium.

SETTING

Population-based.

MAIN OUTCOME

Computed tomography measured liver attenuation.

RESULTS

Carriers of rs4841132-A (frequency 2%-19%) do not show increased hepatic steatosis; they have increased liver attenuation indicative of increased glycogen deposition. rs4841132 falls in a noncoding RNA LOC157273 ~190 kb upstream of PPP1R3B. We demonstrate that rs4841132-A increases PPP1R3B through a cis genetic effect. Using CRISPR/Cas9 we engineered a 105-bp deletion including rs4841132-A in human hepatocarcinoma cells that increases PPP1R3B, decreases LOC157273, and increases glycogen perfectly mirroring the human disease. Overexpression of PPP1R3B or knockdown of LOC157273 increased glycogen but did not result in decreased LOC157273 or increased PPP1R3B, respectively, suggesting that the effects may not all occur via affecting RNA levels. Based on electronic health record (EHR) data, rs4841132-A associates with all components of the metabolic syndrome (MetS). However, rs4841132-A associated with decreased low-density lipoprotein (LDL) cholesterol and risk for myocardial infarction (MI). A metabolic signature for rs4841132-A includes increased glycine, lactate, triglycerides, and decreased acetoacetate and beta-hydroxybutyrate.

CONCLUSIONS

These results show that rs4841132-A promotes a hepatic glycogen storage disease by increasing PPP1R3B and decreasing LOC157273. rs4841132-A promotes glycogen accumulation and development of MetS but lowers LDL cholesterol and risk for MI. These results suggest that elevated hepatic glycogen is one cause of MetS that does not invariably promote MI.

Differences in geometric strength at the contralateral hip between men with hip fracture and non-fractured comparators

Older men sustain excess bone mineral density (BMD) declines after hip fracture; however, BMD provides no information on mechanical structure and strength. The aim was to assess whether changes in hip bone geometry in older men after hip fracture differ than that expected with aging. Two cohorts were used: Baltimore Hip Studies 7th cohort (BHS-7) and Baltimore Men's Osteoporosis Study (MOST). The sample (N = 170) included older Caucasian men with hip fracture that were propensity score matched (1:1) to community-dwelling non-fractured comparators. Hip Structural Analysis (HSA) calculated aerial BMD and metrics of bone structural strength: cross-sectional bone area (CSA), cortical outer diameter (OD), section modulus (SM), and centroid position (CP). Mixed-effect models estimated changes in HSA parameters and adjusted robust regression models evaluated between-cohort differences in annual percent change at the narrow neck (NN), intertrochanteric (IT), and femoral shaft (FS). Hip fracture was associated with statistically greater declines in NN CSA (β = -2.818; 95% CI: -3.300%, -2.336%), SM (β = -1.896%; 95% CI: -2.711%, -1.080%) and CP (β = -0.884%; 95% CI: -0.889%, -0.880%) and significantly larger increases in NN OD (β = 0.187%; 95% CI: 0.185%, 0.190%). Differences in IT HSA parameters were like the NN but larger in magnitude, while there were favorable changes in FS geometry where fragility fractures are rare. Findings indicate there are declines in bone structure and strength at the NN and IT regions of the proximal femur in older men during hip fracture recovery that far exceed what occurs during normal aging.

Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure

Heart failure (HF) is a leading cause of morbidity and mortality worldwide. A small proportion of HF cases are attributable to monogenic cardiomyopathies and existing genome-wide association studies (GWAS) have yielded only limited insights, leaving the observed heritability of HF largely unexplained. We report results from a GWAS meta-analysis of HF comprising 47,309 cases and 930,014 controls. Twelve independent variants at 11 genomic loci are associated with HF, all of which demonstrate one or more associations with coronary artery disease (CAD), atrial fibrillation, or reduced left ventricular function, suggesting shared genetic aetiology. Functional analysis of non-CAD-associated loci implicate genes involved in cardiac development (MYOZ1, SYNPO2L), protein homoeostasis (BAG3), and cellular senescence (CDKN1A). Mendelian randomisation analysis supports causal roles for several HF risk factors, and demonstrates CAD-independent effects for atrial fibrillation, body mass index, and hypertension. These findings extend our knowledge of the pathways underlying HF and may inform new therapeutic strategies.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Author Correction: New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A catalog of genetic loci associated with kidney function from analyses of a million individuals

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries

Reduced lung function predicts mortality and is key to the diagnosis of chronic obstructive pulmonary disease (COPD). In a genome-wide association study in 400,102 individuals of European ancestry, we define 279 lung function signals, 139 of which are new. In combination, these variants strongly predict COPD in independent populations. Furthermore, the combined effect of these variants showed generalizability across smokers and never smokers, and across ancestral groups. We highlight biological pathways, known and potential drug targets for COPD and, in phenome-wide association studies, autoimmune-related and other pleiotropic effects of lung function-associated variants. This new genetic evidence has potential to improve future preventive and therapeutic strategies for COPD.

Disentangling the genetics of lean mass

Background

Lean body mass (LM) plays an important role in mobility and metabolic function. We previously identified five loci associated with LM adjusted for fat mass in kilograms. Such an adjustment may reduce the power to identify genetic signals having an association with both lean mass and fat mass.

Objectives

To determine the impact of different fat mass adjustments on genetic architecture of LM and identify additional LM loci.

Methods

We performed genome-wide association analyses for whole-body LM (20 cohorts of European ancestry with n = 38,292) measured using dual-energy X-ray absorptiometry) or bioelectrical impedance analysis, adjusted for sex, age, age2, and height with or without fat mass adjustments (Model 1 no fat adjustment; Model 2 adjustment for fat mass as a percentage of body mass; Model 3 adjustment for fat mass in kilograms).

Results

Seven single-nucleotide polymorphisms (SNPs) in separate loci, including one novel LM locus (TNRC6B), were successfully replicated in an additional 47,227 individuals from 29 cohorts. Based on the strengths of the associations in Model 1 vs Model 3, we divided the LM loci into those with an effect on both lean mass and fat mass in the same direction and refer to those as "sumo wrestler" loci (FTO and MC4R). In contrast, loci with an impact specifically on LM were termed "body builder" loci (VCAN and ADAMTSL3). Using existing available genome-wide association study databases, LM increasing alleles of SNPs in sumo wrestler loci were associated with an adverse metabolic profile, whereas LM increasing alleles of SNPs in "body builder" loci were associated with metabolic protection.

Conclusions

In conclusion, we identified one novel LM locus (TNRC6B). Our results suggest that a genetically determined increase in lean mass might exert either harmful or protective effects on metabolic traits, depending on its relation to fat mass.

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes

Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10-8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10-17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10-15), HERC2-OCA2 (P = 4.2 × 10-12), SLC45A2 (P = 1.7 × 10-10), IL13 (P = 2.8 × 10-9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10-9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10-8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10-7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

Older men who sustain a hip fracture experience greater declines in bone mineral density at the contralateral hip than non-fractured comparators

Men experience declining bone mineral density (BMD) after hip fracture; however, changes attributable to fracture are unknown. This study evaluated the excess BMD decline attributable to hip fracture among older men. Older men with hip fracture experienced accelerated BMD declines and are at an increased risk of secondary fractures.

INTRODUCTION

The objective was to determine whether bone mineral density (BMD) changes in men after hip fracture exceed that expected with aging.

METHODS

Two cohorts were used: Baltimore Hip Studies 7th cohort (BHS-7) and Baltimore Men's Osteoporosis Study (MOST). BHS-7 recruited older adults (N = 339) hospitalized for hip fracture; assessments occurred within 22 days of admission and at 2, 6, and 12 months follow-up. MOST enrolled age-eligible men (N = 694) from population-based listings; data were collected at a baseline visit and a second visit that occurred between 10 and 31 months later. The combined sample (n = 452) consisted of Caucasian men from BHS-7 (n = 89) and MOST (n = 363) with ≥ 2 dual-energy X-ray absorptiometry scans and overlapping ranges of age, height, and weight. Mixed-effect models estimated rates of BMD change, and generalized linear models evaluated differences in annual bone loss at the total hip and femoral neck between cohorts.

RESULTS

Adjusted changes in total hip and femoral neck BMD were - 4.16% (95% CI, - 4.87 to - 3.46%) and - 4.90% (95% CI, - 5.88 to - 3.92%) in BHS-7 participants; - 1.57% (95% CI, - 2.19 to - 0.96%) and - 0.99% (95% CI, - 1.88 to - 0.10%) in MOST participants; and statistically significant (P < 0.001) between-group differences in change were - 2.59% (95% CI, - 3.26 to - 1.91%) and - 3.91% (95% CI, - 4.83 to - 2.98%), respectively.

CONCLUSION

Hip fracture in older men is associated with accelerated BMD declines at the non-fractured hip that are greater than those expected during aging, and pharmacological interventions in this population to prevent secondary fractures may be warranted.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

Exome-chip meta-analysis identifies association between variation in ANKRD26 and platelet aggregation

Previous genome-wide association studies (GWAS) have identified several variants associated with platelet function phenotypes; however, the proportion of variance explained by the identified variants is mostly small. Rare coding variants, particularly those with high potential for impact on protein structure/function, may have substantial impact on phenotype but are difficult to detect by GWAS. The main purpose of this study was to identify low frequency or rare variants associated with platelet function using genotype data from the Illumina HumanExome Bead Chip. Three family-based cohorts of European ancestry, including ~4,000 total subjects, comprised the discovery cohort and two independent cohorts, one of European and one of African American ancestry, were used for replication. Optical aggregometry in platelet-rich plasma was performed in all the discovery cohorts in response to adenosine diphosphate (ADP), epinephrine, and collagen. Meta-analyses were performed using both gene-based and single nucleotide variant association methods. The gene-based meta-analysis identified a significant association (P = 7.13 × 10^-7) between rare genetic variants in ANKRD26 and ADP-induced platelet aggregation. One of the ANKRD26 SNVs - rs191015656, encoding a threonine to isoleucine substitution predicted to alter protein structure/function, was replicated in Europeans. Aggregation increases of ~20-50% were observed in heterozygotes in all cohorts. Novel genetic signals in ABCG1 and HCP5 were also associated with platelet aggregation to ADP in meta-analyses, although only results for HCP5 could be replicated. The SNV in HCP5 intersects epigenetic signatures in CD41+ megakaryocytes suggesting a new functional role in platelet biology for HCP5. This is the first study to use gene-based association methods from SNV array genotypes to identify rare variants related to platelet function. The molecular mechanisms and pathophysiological relevance for the identified genetic associations requires further study.

Genetic Determinants of Radiographic Knee Osteoarthritis in African Americans

OBJECTIVE

The etiology of knee osteoarthritis (OA), the most common form of arthritis, is complex and may differ by race or ethnicity. In recent years, genetic studies have identified many genetic variants associated with OA, but nearly all the studies were conducted in European whites and Asian Americans. Few studies have focused on the genetics of knee OA in African Americans.

METHODS

We performed a genome-wide association study of radiographic knee OA in 1217 African Americans from 2 North American cohort studies: 590 subjects from the Johnston County Osteoarthritis Project and 627 subjects from the Osteoarthritis Initiative. Analyses were conducted in each cohort separately and combined in an inverse variance fixed effects metaanalysis, which were then included in pathway analyses. We additionally tested 12 single-nucleotide polymorphisms robustly associated with OA in European white populations for association in African Americans.

RESULTS

We identified a genome-wide significant variant in LINC01006 (minor allele frequency 12%; p = 4.11 × 10-9) that is less common in European white populations (minor allele frequency < 3%). Five other independent loci reached suggestive significance (p < 1 × 10-6). In pathway analyses, dorsal/ventral neural tube patterning and iron ion transport pathways were significantly associated with knee OA in African Americans (false discovery rate < 0.05). We found no evidence that previously reported OA susceptibility variants in European whites were associated with knee OA in African Americans.

CONCLUSION

These results highlight differences in the genetic architecture of knee OA between African American and European whites. This finding underscores the need to include more diverse populations in OA genetics studies.

Oxaliplatin-Induced Peripheral Neuropathy and Identification of Unique Severity Groups in Colorectal Cancer

CONTEXT

Oxaliplatin-induced peripheral neuropathy (OIPN) is a dose-limiting toxicity of oxaliplatin and affects most colorectal cancer patients. OIPN is commonly evaluated by patient symptom report, using scales to reflect impairment. They do not discriminate between unique grouping of symptoms and signs, which impedes prompt identification of OIPN.

OBJECTIVE

The objective of this study was to identify clusters of symptoms and signs that differentiated underlying clinical severity and segregated patients within our population into OIPN subgroups.

METHODS

Chemotherapy-naive colorectal cancer patients (N = 148) receiving oxaliplatin were administered the Total Neuropathy Score clinical (TNSc^©), which includes symptom report (sensory, motor, autonomic) and sensory examination (pin sense, vibration, reflexes). The TNSc was administered before chemotherapy initiation (T0) and after cumulative doses of oxaliplatin 510-520 mg/m² (T1) and 1020-1040 mg/m² of oxaliplatin (T2). Using mean T2 TNSc scores, latent class analysis grouped patients into OIPN severity cohorts.

RESULTS

Latent class analysis categorized patients into four distinct OIPN groups: low symptoms and low signs (n = 54); low symptoms and intermediate signs (n = 44); low symptoms and high signs (n = 21); and high symptoms and high signs (n = 29). No differences were noted among OIPN groups on age, sex, chemotherapy regimen, or cumulative oxaliplatin dose.

CONCLUSION

We identified OIPN patient groups with distinct symptoms/signs, demonstrating variability of OIPN presentation regardless of cumulative oxaliplatin dose. Over half of the sample had positive findings on OIPN examination despite little or no symptoms. Sensory examination of all patients receiving oxaliplatin is indicated for timely identification of OIPN, which will allow earlier symptom management.

Large meta-analysis of genome-wide association studies identifies five loci for lean body mass

Lean body mass, consisting mostly of skeletal muscle, is important for healthy aging. We performed a genome-wide association study for whole body (20 cohorts of European ancestry with n = 38,292) and appendicular (arms and legs) lean body mass (n = 28,330) measured using dual energy X-ray absorptiometry or bioelectrical impedance analysis, adjusted for sex, age, height, and fat mass. Twenty-one single-nucleotide polymorphisms were significantly associated with lean body mass either genome wide (p < 5 × 10⁻⁸) or suggestively genome wide (p < 2.3 × 10⁻⁶). Replication in 63,475 (47,227 of European ancestry) individuals from 33 cohorts for whole body lean body mass and in 45,090 (42,360 of European ancestry) subjects from 25 cohorts for appendicular lean body mass was successful for five single-nucleotide polymorphisms in/near HSD17B11, VCAN, ADAMTSL3, IRS1, and FTO for total lean body mass and for three single-nucleotide polymorphisms in/near VCAN, ADAMTSL3, and IRS1 for appendicular lean body mass. Our findings provide new insight into the genetics of lean body mass.

Lean body mass is a highly heritable trait and is associated with various health conditions. Here, Kiel and colleagues perform a meta-analysis of genome-wide association studies for whole body lean body mass and find five novel genetic loci to be significantly associated.

Genome-Wide Association Study of Radiographic Knee Osteoarthritis in North American Caucasians

OBJECTIVE

A major barrier to genetic studies of osteoarthritis (OA) is the need to obtain large numbers of individuals with standardized radiographic evaluations for OA. To address this gap, we performed a genome-wide association study (GWAS) of radiographically defined tibiofemoral knee OA in 3,898 cases and 3,168 controls from 4 well-characterized North American cohorts, and we performed replication analysis of previously reported OA loci.

METHODS

We performed meta-analysis using a 2-stage design. Stage 1 (discovery) consisted of a GWAS meta-analysis of radiographic knee OA carried out in the Osteoarthritis Initiative and the Johnston County Osteoarthritis Project. Knee OA was defined as definite osteophytes and possible joint space narrowing or total joint replacement in one or both knees. Stage 2 (validation) was performed in the Multicenter Osteoarthritis Study and the Genetics of Osteoarthritis study. We genotyped lead meta-analysis variants (P ≤ 1 × 10^-4 ) from stage 1 and tested the association between these variants and knee OA. We then combined results from all cohorts in a meta-analysis.

RESULTS

Lead variants from stage 1, representing 49 unique loci, were analyzed in stage 2; none met genome-wide significance in the combined analysis of stage 1 and stage 2. We validated 1 locus (rs4867568 near LSP1P3) with nominal significance (P < 0.05), which was also our top finding in the combined meta-analysis (odds ratio [OR] 0.84 [95% confidence interval (95% CI) 0.79-0.91], P = 3.02 × 10^-6 ). We observed nominally significant associations (P < 0.05) with 3 previously reported OA loci: rs143383 in GDF5 (OR 1.12 [95% CI 1.04-1.21], P = 2.13 × 10^-3 ), rs835487 in CHST11 (OR 0.93 [95% CI 0.85-0.99], P = 0.03), and rs8044769 in FTO (OR 1.10 [95% CI 1.03-1.19], P = 6.13 × 10^-3 ).

CONCLUSION

These findings provide suggestive evidence of a novel knee OA locus and confirm previously reported associations in GDF5, CHST11, and FTO.

Association Between Telomere Length and Risk of Cancer and Non-Neoplastic Diseases: A Mendelian Randomization Study

IMPORTANCE

The causal direction and magnitude of the association between telomere length and incidence of cancer and non-neoplastic diseases is uncertain owing to the susceptibility of observational studies to confounding and reverse causation.

OBJECTIVE

To conduct a Mendelian randomization study, using germline genetic variants as instrumental variables, to appraise the causal relevance of telomere length for risk of cancer and non-neoplastic diseases.

DATA SOURCES

Genomewide association studies (GWAS) published up to January 15, 2015.

STUDY SELECTION

GWAS of noncommunicable diseases that assayed germline genetic variation and did not select cohort or control participants on the basis of preexisting diseases. Of 163 GWAS of noncommunicable diseases identified, summary data from 103 were available.

DATA EXTRACTION AND SYNTHESIS

Summary association statistics for single nucleotide polymorphisms (SNPs) that are strongly associated with telomere length in the general population.

MAIN OUTCOMES AND MEASURES

Odds ratios (ORs) and 95% confidence intervals (CIs) for disease per standard deviation (SD) higher telomere length due to germline genetic variation.

RESULTS

Summary data were available for 35 cancers and 48 non-neoplastic diseases, corresponding to 420 081 cases (median cases, 2526 per disease) and 1 093 105 controls (median, 6789 per disease). Increased telomere length due to germline genetic variation was generally associated with increased risk for site-specific cancers. The strongest associations (ORs [95% CIs] per 1-SD change in genetically increased telomere length) were observed for glioma, 5.27 (3.15-8.81); serous low-malignant-potential ovarian cancer, 4.35 (2.39-7.94); lung adenocarcinoma, 3.19 (2.40-4.22); neuroblastoma, 2.98 (1.92-4.62); bladder cancer, 2.19 (1.32-3.66); melanoma, 1.87 (1.55-2.26); testicular cancer, 1.76 (1.02-3.04); kidney cancer, 1.55 (1.08-2.23); and endometrial cancer, 1.31 (1.07-1.61). Associations were stronger for rarer cancers and at tissue sites with lower rates of stem cell division. There was generally little evidence of association between genetically increased telomere length and risk of psychiatric, autoimmune, inflammatory, diabetic, and other non-neoplastic diseases, except for coronary heart disease (OR, 0.78 [95% CI, 0.67-0.90]), abdominal aortic aneurysm (OR, 0.63 [95% CI, 0.49-0.81]), celiac disease (OR, 0.42 [95% CI, 0.28-0.61]) and interstitial lung disease (OR, 0.09 [95% CI, 0.05-0.15]).

CONCLUSIONS AND RELEVANCE

It is likely that longer telomeres increase risk for several cancers but reduce risk for some non-neoplastic diseases, including cardiovascular diseases.

Genome-wide analysis of clopidogrel active metabolite levels identifies novel variants that influence antiplatelet response

Clopidogrel is one of the most commonly used therapeutics for the secondary prevention of cardiovascular events in patients with acute coronary syndromes. However, considerable interindividual variation in clopidogrel response has been documented, resulting in suboptimal therapy and an increased risk of recurrent events for some patients. In this investigation, we carried out the first genome-wide association study of circulating clopidogrel active metabolite levels in 513 healthy participants to directly measure clopidogrel pharmacokinetics. We observed that the CYP2C19 locus was the strongest genetic determinant of active metabolite formation (P=9.5×10). In addition, we identified novel genome-wide significant variants on chromosomes 3p25 (rs187941554, P=3.3×10) and 17q11 (rs80343429, P=1.3×10), as well as six additional loci that showed suggestive evidence of association (P≤1.0×10). Four of these loci showed nominal associations with on-clopidogrel ADP-stimulated platelet aggregation (P≤0.05). Evaluation of clopidogrel active metabolite concentration may help identify novel genetic determinants of clopidogrel response, which has implications for the development of novel therapeutics and improved antiplatelet treatment for at-risk patients in the future.

Rare and low-frequency coding variants alter human adult height

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

Prospective Evaluation of Genetic Variation in Platelet Endothelial Aggregation Receptor 1 Reveals Aspirin-Dependent Effects on Platelet Aggregation Pathways

Genetic variation in the platelet endothelial aggregation receptor 1 (PEAR1) gene, most notably rs12041331, is implicated in altered on‐aspirin platelet aggregation and increased cardiovascular event risk. We prospectively tested the effects of aspirin administration at commonly prescribed doses (81, 162, and 324 mg/day) on agonist‐induced platelet aggregation by rs12041331 genotype in 67 healthy individuals. Prior to aspirin administration, rs12041331 minor allele carriers had significantly reduced adenosine diphosphate (ADP)‐induced platelet aggregation compared with noncarriers (P = 0.03) but was not associated with other platelet pathways. In contrast, rs12041331 was significantly associated with on‐aspirin platelet aggregation when collagen and epinephrine were used to stimulate platelet aggregation (P < 0.05 for all associations), but not ADP. The influence of PEAR1 rs12041331 on platelet aggregation is pathway‐specific and is altered by aspirin at therapeutic doses, but not in a dose‐dependent manner. Additional studies are needed to determine the impact of PEAR1 on cardiovascular events in aspirin‐treated patients.

The CAPN2/CAPN8 Locus on Chromosome 1q Is Associated with Variation in Serum Alpha-Carotene Concentrations

BACKGROUND/AIMS

Alpha-carotene is a provitamin A carotenoid present in fruits and vegetables. Higher serum concentrations of α-carotene have been associated with lower risk of cancer and all-cause mortality. Previous studies have suggested that genetic variants influence serum concentrations of provitamin A carotenoids, but to date no variants have been robustly associated with serum α-carotene concentrations. The aim of this study was to identify genetic associations with serum α-carotene concentrations using the genome-wide association study (GWAS) approach.

METHODS

A GWAS of serum α-carotene concentrations was conducted in 433 Old Order Amish adults who had consumed a 6-day controlled diet. Linear regression models adjusting for age, gender, and family structure were utilized to evaluate associations between genetic variants and serum α-carotene concentrations.

RESULTS

Genome-wide significant associations with α-carotene concentrations were observed for loci on chromosome 1q41 between the genes CAPN2 and CAPN8 (rs12137025, p = 3.55 × 10-8), chromosome 2p21 in PRKCE (rs2594495, p = 1.01 × 10-8), and chromosome 4q34 (rs17830069, p = 2.89 × 10-8).

CONCLUSIONS

We identified 3 novel loci associated with serum α-carotene concentrations among a population that consumed a controlled diet. While replication is necessary, the CAPN2/CAPN8 locus provides compelling evidence for an association with serum α-carotene concentrations and may suggest a relationship with the development and progression of cancers.

SOS2 and ACP1 Loci Identified through Large-Scale Exome Chip Analysis Regulate Kidney Development and Function

Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium (n_Stage1: 111,666; n_Stage2: 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea (PPM1J, EDEM3, ACP1, SPEG, EYA4, CYP1A1, and ATXN2L; P_Stage1<3.7×10^-7), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 (P=5.4×10^-8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2-knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.

A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Multiethnic Exome-Wide Association Study of Subclinical Atherosclerosis

BACKGROUND

The burden of subclinical atherosclerosis in asymptomatic individuals is heritable and associated with elevated risk of developing clinical coronary heart disease. We sought to identify genetic variants in protein-coding regions associated with subclinical atherosclerosis and the risk of subsequent coronary heart disease.

METHODS AND RESULTS

We studied a total of 25 109 European ancestry and African ancestry participants with coronary artery calcification (CAC) measured by cardiac computed tomography and 52 869 participants with common carotid intima-media thickness measured by ultrasonography within the CHARGE Consortium (Cohorts for Heart and Aging Research in Genomic Epidemiology). Participants were genotyped for 247 870 DNA sequence variants (231 539 in exons) across the genome. A meta-analysis of exome-wide association studies was performed across cohorts for CAC and carotid intima-media thickness. APOB p.Arg3527Gln was associated with 4-fold excess CAC (P=3×10^-10). The APOE ε2 allele (p.Arg176Cys) was associated with both 22.3% reduced CAC (P=1×10^-12) and 1.4% reduced carotid intima-media thickness (P=4×10^-14) in carriers compared with noncarriers. In secondary analyses conditioning on low-density lipoprotein cholesterol concentration, the ε2 protective association with CAC, although attenuated, remained strongly significant. Additionally, the presence of ε2 was associated with reduced risk for coronary heart disease (odds ratio 0.77; P=1×10^-11).

CONCLUSIONS

Exome-wide association meta-analysis demonstrates that protein-coding variants in APOB and APOE associate with subclinical atherosclerosis. APOE ε2 represents the first significant association for multiple subclinical atherosclerosis traits across multiple ethnicities, as well as clinical coronary heart disease.

Clopidogrel Improves Skin Microcirculatory Endothelial Function in Persons With Heightened Platelet Aggregation

Background

Platelet activation can lead to enhanced oxidative stress, inflammatory response, and endothelial dysfunction. To quantify the effects of platelet inhibition on endothelial function, we assessed platelet activity of healthy persons before and after clopidogrel administration and evaluated its effects on endothelial function. We hypothesized that clopidogrel, by attenuating platelet activity, would result in enhanced endothelial function.

Methods and Results

Microcirculatory endothelial function was quantified by laser Doppler flowmetry (LDF) mediated by thermal hyperemia (TH) and postocclusive reactive hyperemia, respectively, in 287 and 241 relatively healthy and homogenous Old Order Amish persons. LDF and platelet aggregation measures were obtained at baseline and after 7 days of clopidogrel administration. Our primary outcome was percentage change in post‐ versus preclopidogrel LDF measures. Preclopidogrel TH‐LDF and platelet aggregation were higher in women than in men (P<0.001). Clopidogrel administration was associated with ≈2‐fold higher percentage change in TH‐LDF in participants with high versus low baseline platelet aggregation (39.4±10.1% versus 17.4±5.6%, P=0.03). Clopidogrel also increased absolute TH‐LDF measures in persons with high platelet aggregation (1757±766 to 2154±1055, P=0.03), with a more prominent effect in women (1909±846 to 2518±1048, P=0.001). There was no evidence that clopidogrel influenced postocclusive reactive hyperemia LDF measures.

Conclusions

The administration of clopidogrel in healthy persons with high baseline platelet aggregation results in improved TH‐induced microcirculatory endothelial function. These data suggest that clopidogrel may have a beneficial effect on microcirculatory endothelial function, presumably through antiplatelet activity, and may confer additional vascular benefits.

Clinical Trial Registration

URL: https://www.clinicaltrials.gov. Unique identifier: NCT00799396.

Transcriptional Regulation of Frizzled-1 in Human Osteoblasts by Sp1

The wingless pathway has a powerful influence on bone metabolism and is a therapeutic target in skeletal disorders. Wingless signaling is mediated in part through the Frizzled (FZD) receptor family. FZD transcriptional regulation is poorly understood. Herein we tested the hypothesis that Sp1 plays an important role in the transcriptional regulation of FZD1 expression in osteoblasts and osteoblast mineralization. To test this hypothesis, we conducted FZD1 promoter assays in Saos2 cells with and without Sp1 overexpression. We found that Sp1 significantly up-regulates FZD1 promoter activity in Saos2 cells. Chromatin immunoprecipitation (ChIP) and electrophoretic mobility shift (EMSA) assays identified a novel and functional Sp1 binding site at -44 to -40 from the translation start site in the FZD1 promoter. The Sp1-dependent activation of the FZD1 promoter was abolished by mithramycin A (MMA), an antibiotic affecting both Sp1 binding and Sp1 protein levels in Saos2 cells. Similarly, down-regulation of Sp1 in hFOB cells resulted in less FZD1 expression and lower alkaline phosphatase activity. Moreover, over-expression of Sp1 increased FZD1 expression and Saos2 cell mineralization while MMA decreased Sp1 and FZD1 expression and Saos2 cell mineralization. Knockdown of FZD1 prior to Sp1 overexpression partially abolished Sp1 stimulation of osteoblast differentiation markers. Taken together, our results suggest that Sp1 plays a role in human osteoblast differentiation and mineralization, which is at least partially mediated by Sp1-dependent transactivation of FZD1.

The Interplay of Genetics, Behavior, and Pain with Depressive Symptoms in the Elderly

PURPOSE OF STUDY

About 25% of older adults suffer from depressive symptoms. Commonly studied candidate genes associated with depression include those that influence serotonin (SLC6A4), dopamine (COMT), or neuroplasticity (BDNF, NTRK3). However, the majority of candidate gene studies do not consider the interplay of genetics, demographic, clinical, and behavioral factors and how they jointly contribute to depressive symptoms among older adults. The purpose of this study was to gain a more comprehensive understanding of depressive symptoms among older adults.

DESIGN AND METHODS

In this descriptive study, demographic, behavioral, and clinical characteristics (age, gender, comorbidities, volunteering, physical activity, pain, and fear of falling) were obtained via interview of 114 residents in a continuing care retirement community. Peripheral whole blood was collected for DNA extraction. We examined common single nucleotide polymorphisms (SNPs) in the aforementioned genes using path analyses.

RESULTS

SNPs in the NTRK3 gene, pain, physical activity, and fear of falling were directly associated with depressive symptoms in older adults. Those who had polymorphisms in the NTRK3 gene, pain, fear of falling, and were less physically active were more likely to exhibit depressive symptoms. None of the SNPs in SLC6A4, COMT, or BDNF genes were significantly associated with depressive symptoms.

IMPLICATIONS

Our use of a path analysis to examine a biopsychosocial model of depressive symptoms provided the opportunity to describe a comprehensive clinical picture of older adults at risk for depressive symptoms. Thus, interventions could be implemented to identify older adults at risk for depressive symptoms.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Abstract 32: Analysis of Serum Clopidogrel Active Metabolite Concentration Identifies Novel Genetic Variants Associated With Clopidogrel Pharmacokinetics

Clopidogrel with aspirin is standard of care for the secondary prevention of ischemic events in patients with acute coronary syndromes. However, substantial inter-individual variation in clopidogrel response results in sub-optimal therapy in some patients and an increased risk of recurrent events. Up to 73% of the variation in clopidogrel efficacy is found to be heritable, and while key genetic variants such as CYP2C19 *2 have been identified, a large portion of the genetic heritability remains unaccounted for. To date, few studies have used clopidogrel active metabolite to evaluate clopidogrel pharmacogenomics. The evaluation of this biomarker may be particularly informative as it enables us to focus almost exclusively on clopidogrel pharmacokinetics and may allow new discoveries leading to a more thorough understanding of clopidogrel activation.

Serum clopidogrel active metabolite concentration was measured by HPLC-MS/MS in 506 healthy, Amish participants of the Pharmacogenomics of Anti-Platelet Intervention (PAPI) Study after a 7-day 75-mg/day clopidogrel intervention. Genome-wide genotyping was performed using Affymetrix 5.0 and 6.0 arrays, and SNPs were imputed to the 1000 Genomes cosmopolitan panel. Mixed models analysis was conducted under an additive model and adjusted for age, sex, BMI, and relatedness among study participants with MMAP (http://edn.som.umaryland.edu/mmap).

The most significant signal of our GWAS was the previously identified CYP2C19 / CYP2C9 locus (P= 9.2x10 -15 ), and was fully accounted for by the CYP2C19 *2 variant. In addition to this signal, novel genome-wide significant signals were identified in loci on Chr. 3 (rs187941554, P=4.7x10 -11 ), Chr. 7 (rs73407739, p=3.6x10 -8 ), and Chr. 17 (rs80343429, P=3.7x10 -8 ).

In this study we conducted the first ever GWAS of serum clopidogrel active metabolite concentration. Novel associations for clopidogrel efficacy were identified, and further study is required to determine the impact of these variants on clinical care Continual analysis of clopidogrel active metabolite concentration will contribute to a better understanding of the clopidogrel metabolic pathway, and potentially lay the groundwork for improved anti-platelet treatment for CVD patients in the future.

Assignment of Functional Relevance to Genes at Type 2 Diabetes-Associated Loci Through Investigation of β-Cell Mass Deficits

Type 2 diabetes (T2D) has been associated with a large number of genomic loci, many of which encompass multiple genes without a definitive causal gene. This complexity has hindered efforts to clearly identify functional candidate genes and interpret their role in mediating susceptibility to disease. Here we examined the relevance of individual genes found at T2D-associated loci by assessing their potential contribution to a phenotype relevant to the disease state: production and maintenance of β-cell mass. Using transgenic zebrafish in which β-cell mass could be rapidly visualized in vivo, we systematically suppressed the expression of orthologs of genes found at T2D-associated genomic loci. Overall, we tested 67 orthologs, many of which had no known relevance to β-cell mass, at 62 human T2D-associated loci, including eight loci with multiple candidate genes. In total we identified 25 genes that were necessary for proper β-cell mass, providing functional evidence for their role in a physiological phenotype directly related to T2D. Of these, 16 had not previously been implicated in the regulation of β-cell mass. Strikingly, we identified single functional candidate genes at the majority of the loci for which multiple genes were analyzed. Further investigation into the contribution of the 25 genes to the adaptive capacity of β-cells suggested that the majority of genes were not required for glucose-induced expansion of β-cell mass but were significantly necessary for the regeneration of β-cells. These findings suggest that genetically programmed deficiencies in β-cell mass may be related to impaired maintenance. Finally, we investigated the relevance of our findings to human T2D onset in diabetic individuals from the Old Order Amish and found that risk alleles in β-cell mass genes were associated with significantly younger age of onset and lower body mass index. Taken together, our study offers a functional approach to assign relevance to genes at T2D-associated loci and offers experimental evidence for the defining role of β-cell mass maintenance in genetic susceptibility to T2D onset.

Gene Expression Differences Between Offspring of Long-Lived Individuals and Controls in Candidate Longevity Regions: Evidence for PAPSS2 as a Longevity Gene

Although there is compelling evidence for a genetic contribution to longevity, identification of specific genes that robustly associate with longevity has been a challenge. In order to identify longevity-enhancing genes, we measured differential gene expression between offspring of long-lived Amish (older than 90 years; cases, n = 128) and spouses of these offspring (controls, n = 121) and correlated differentially expressed transcripts with locations of longevity-associated variants detected in a prior genome-wide association study (GWAS) of survival to age 90. Expression of one of these transcripts, 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (PAPSS2), was significantly higher in offspring versus controls (4×10(-4)) and this association was replicated using quantitative real-time polymerase chain reaction. PAPSS2, a sulfation enzyme located on chromosome 10, is ~80kb upstream of the PAPSS2 transcription start site. We found evidence of cis-expression for the originally reported GWAS SNP and PAPSS2 Monogenic conditions linked to PAPSS2 include andrenocortical androgen excess resulting in premature pubarche and skeletal dysplasias, both of which have premature aging features. In summary, these findings provide novel evidence for PAPSS2 as a longevity locus and illustrate the value of harnessing multiple "-omic" approaches to identify longevity candidates.

A Common Variant in the SETD7 Gene Predicts Serum Lycopene Concentrations

Dietary intake and higher serum concentrations of lycopene have been associated with lower incidence of prostate cancer and other chronic diseases. Identifying determinants of serum lycopene concentrations may thus have important public health implications. Prior studies have suggested that serum lycopene concentrations are under partial genetic control. The goal of this research was to identify genetic predictors of serum lycopene concentrations using the genome-wide association study (GWAS) approach among a sample of 441 Old Order Amish adults that consumed a controlled diet. Linear regression models were utilized to evaluate associations between genetic variants and serum concentrations of lycopene. Variant rs7680948 on chromosome 4, located in the intron region of the SETD7 gene, was significantly associated with serum lycopene concentrations (p = 3.41 × 10⁻⁹). Our findings also provided nominal support for the association previously noted between SCARB1 and serum lycopene concentrations, although with a different SNP (rs11057841) in the region. This study identified a novel locus associated with serum lycopene concentrations and our results raise a number of intriguing possibilities regarding the nature of the relationship between SETD7 and lycopene, both of which have been independently associated with prostate cancer. Further investigation into this relationship might help provide greater mechanistic understanding of these associations.

Corrigendum: Rare coding variants and X-linked loci associated with age at menarche

Nature Communications 6, Article number: 7756 (2015); Published 4 August 2015; Updated 17 December 2015 In the Results section and in the legend of Table 1 of this Article, the company deCODE genetics, Inc. is incorrectly referred to as ‘Diabetes Epidemiology: collaborative analysis of Diagnostic criteria in Europe’.

Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Development of a physiology-directed population pharmacokinetic and pharmacodynamic model for characterizing the impact of genetic and demographic factors on clopidogrel response in healthy adults

Clopidogrel (Plavix®), is a widely used antiplatelet agent, which shows high inter-individual variability in treatment response in patients following the standard dosing regimen. In this study, a physiology-directed population pharmacokinetic/pharmacodynamic (PK/PD) model was developed based on clopidogrel and clopidogrel active metabolite (clop-AM) data from the PAPI and the PGXB2B studies using a step-wise approach in NONMEM (version 7.2). The developed model characterized the in vivo disposition of clopidogrel, its bioactivation into clop-AM in the liver and subsequent platelet aggregation inhibition in the systemic circulation reasonably well. It further allowed the identification of covariates that significantly impact clopidogrel's dose-concentration-response relationship. In particular, CYP2C19 intermediate and poor metabolizers converted 26.2% and 39.5% less clopidogrel to clop-AM, respectively, compared to extensive metabolizers. In addition, CES1 G143E mutation carriers have a reduced CES1 activity (82.9%) compared to wild-type subjects, which results in a significant increase in clop-AM formation. An increase in BMI was found to significantly decrease clopidogrel's bioactivation, whereas increased age was associated with increased platelet reactivity. Our PK/PD model analysis suggests that, in order to optimize clopidogrel dosing on a patient-by-patient basis, all of these factors have to be considered simultaneously, e.g. by using quantitative clinical pharmacology tools.

Oxylipid Profile of Low-Dose Aspirin Exposure: A Pharmacometabolomics Study

Background

While aspirin is a well‐established and generally effective anti‐platelet agent, considerable inter‐individual variation in drug response exists, for which mechanisms are not completely understood. Metabolomics allows for extensive measurement of small molecules in biological samples, enabling detailed mapping of pathways involved in drug response.

Methods and Results

We used a mass‐spectrometry‐based metabolomics platform to investigate the changes in the serum oxylipid metabolome induced by an aspirin intervention (14 days, 81 mg/day) in healthy subjects (n=156). We observed a global decrease in serum oxylipids in response to aspirin (25 metabolites decreased out of 30 measured) regardless of sex. This decrease was concomitant with a significant decrease in serum linoleic acid levels (−19%, P=1.3×10⁻⁵), one of the main precursors for oxylipid synthesis. Interestingly, several linoleic acid‐derived oxylipids were not significantly associated with arachidonic‐induced ex vivo platelet aggregation, a widely accepted marker of aspirin response, but were significantly correlated with platelet reactivity in response to collagen.

Conclusions

Together, these results suggest that linoleic acid‐derived oxylipids may contribute to the non‐COX1 mediated variability in response to aspirin. Pharmacometabolomics allowed for more comprehensive interrogation of mechanisms of action of low dose aspirin and of variation in aspirin response.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHR_adjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHR_adjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHR_adjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHR_adjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHR_adjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Adult body size and body shape differ substantially between men and women and change over time. More than 100 genetic variants that influence body mass index (measure of body size) or waist-to-hip ratio (measure of body shape) have been identified. While there is evidence that some genetic loci affect body shape differently in men than in women, little is known about whether genetic effects differ in older compared to younger adults, and whether such changes differ between men and women. Therefore, we conducted a systematic genome-wide search, including 114 studies (>320,000 individuals), to specifically identify genetic loci with age- and or sex-dependent effects on body size and shape. We identified 15 loci of which the effect on BMI was different in older compared to younger adults, whereas we found no evidence for loci with different effects in men compared to women. The opposite was seen for body shape as we identified 44 loci of which the effect on waist-to-hip ratio differed between men and women, but no difference between younger and older adults were observed. Our observations may provide new insights into the biology that underlies weight change with age or the sexual dimorphism of body shape.

Genetic Variation in the Platelet Endothelial Aggregation Receptor 1 Gene Results in Endothelial Dysfunction

Platelet Endothelial Aggregation Receptor 1 (PEAR1) is a newly identified membrane protein reported to be involved in multiple vascular and thrombotic processes. While most studies to date have focused on the effects of this receptor in platelets, PEAR1 is located in multiple tissues including the endothelium, where it is most highly expressed. Our first objective was to evaluate the role of PEAR1 in endothelial function by examining flow-mediated dilation of the brachial artery in 641 participants from the Heredity and Phenotype Intervention Heart Study. Our second objective was to further define the impact of PEAR1 on cardiovascular disease computationally through meta-analysis of 75,000 microarrays, yielding insights regarding PEAR1 function, and predictions of phenotypes and diseases affected by PEAR1 dysregulation. Based on the results of this meta-analysis we examined whether genetic variation in PEAR1 influences endothelial function using an ex vivo assay of endothelial cell migration. We observed a significant association between rs12041331 and flow-mediated dilation in participants of the Heredity and Phenotype Intervention Heart Study (P = 0.02). Meta-analysis results revealed that PEAR1 expression is highly correlated with several genes (e.g. ANG2, ACVRL1, ENG) and phenotypes (e.g. endothelial cell migration, angiogenesis) that are integral to endothelial function. Functional validation of these results revealed that PEAR1 rs12041331 is significantly associated with endothelial migration (P = 0.04). Our results suggest for the first time that genetic variation of PEAR1 is a significant determinant of endothelial function through pathways implicated in cardiovascular disease.

Rare coding variants and X-linked loci associated with age at menarche

More than 100 loci have been identified for age at menarche by genome-wide association studies; however, collectively these explain only ∼3% of the trait variance. Here we test two overlooked sources of variation in 192,974 European ancestry women: low-frequency protein-coding variants and X-chromosome variants. Five missense/nonsense variants (in ALMS1/LAMB2/TNRC6A/TACR3/PRKAG1) are associated with age at menarche (minor allele frequencies 0.08-4.6%; effect sizes 0.08-1.25 years per allele; P<5 × 10(-8)). In addition, we identify common X-chromosome loci at IGSF1 (rs762080, P=9.4 × 10(-13)) and FAAH2 (rs5914101, P=4.9 × 10(-10)). Highlighted genes implicate cellular energy homeostasis, post-transcriptional gene silencing and fatty-acid amide signalling. A frequently reported mutation in TACR3 for idiopathic hypogonatrophic hypogonadism (p.W275X) is associated with 1.25-year-later menarche (P=2.8 × 10(-11)), illustrating the utility of population studies to estimate the penetrance of reportedly pathogenic mutations. Collectively, these novel variants explain ∼0.5% variance, indicating that these overlooked sources of variation do not substantially explain the 'missing heritability' of this complex trait.