We conducted genome-wide association studies (GWAS) of relative intake from the macronutrients fat, protein, carbohydrates, and sugar in over 235,000 individuals of European ancestries. We identified 21 unique, approximately independent lead SNPs. Fourteen lead SNPs are uniquely associated with one macronutrient at genome-wide significance (P < 5 × 10-8), while five of the 21 lead SNPs reach suggestive significance (P < 1 × 10-5) for at least one other macronutrient. While the phenotypes are genetically correlated, each phenotype carries a partially unique genetic architecture. Relative protein intake exhibits the strongest relationships with poor health, including positive genetic associations with obesity, type 2 diabetes, and heart disease (rg ≈ 0.15-0.5). In contrast, relative carbohydrate and sugar intake have negative genetic correlations with waist circumference, waist-hip ratio, and neighborhood deprivation (|rg| ≈ 0.1-0.3) and positive genetic correlations with physical activity (rg ≈ 0.1 and 0.2). Relative fat intake has no consistent pattern of genetic correlations with poor health but has a negative genetic correlation with educational attainment (rg ≈-0.1). Although our analyses do not allow us to draw causal conclusions, we find no evidence of negative health consequences associated with relative carbohydrate, sugar, or fat intake. However, our results are consistent with the hypothesis that relative protein intake plays a role in the etiology of metabolic dysfunction.

C-reactive protein (CRP) is a circulating biomarker indicative of systemic inflammation. We aimed to evaluate genetic associations with CRP levels among non-European-ancestry populations through discovery, fine-mapping and conditional analyses. A total of 30 503 non-European-ancestry participants from 6 studies participating in the Population Architecture using Genomics and Epidemiology study had serum high-sensitivity CRP measurements and ∼200 000 single nucleotide polymorphisms (SNPs) genotyped on the Metabochip. We evaluated the association between each SNP and log-transformed CRP levels using multivariate linear regression, with additive genetic models adjusted for age, sex, the first four principal components of genetic ancestry, and study-specific factors. Differential linkage disequilibrium patterns between race/ethnicity groups were used to fine-map regions associated with CRP levels. Conditional analyses evaluated for multiple independent signals within genetic regions. One hundred and sixty-three unique variants in 12 loci in overall or race/ethnicity-stratified Metabochip-wide scans reached a Bonferroni-corrected P-value <2.5E-7. Three loci have no (HACL1, OLFML2B) or only limited (PLA2G6) previous associations with CRP levels. Six loci had different top hits in race/ethnicity-specific versus overall analyses. Fine-mapping refined the signal in six loci, particularly in HNF1A. Conditional analyses provided evidence for secondary signals in LEPR, IL1RN and HNF1A, and for multiple independent signals in CRP and APOE. We identified novel variants and loci associated with CRP levels, generalized known CRP associations to a multiethnic study population, refined association signals at several loci and found evidence for multiple independent signals at several well-known loci. This study demonstrates the benefit of conducting inclusive genetic association studies in large multiethnic populations.

Aromatase inhibitors (AIs) reduce breast cancer recurrence and prolong survival, but up to 30% of patients exhibit recurrence. Using a genome-wide association study of patients entered on MA.27, a phase III randomized trial of anastrozole versus exemestane, we identified a single nucleotide polymorphism (SNP) in CUB And Sushi multiple domains 1 (CSMD1) associated with breast cancer-free interval, with the variant allele associated with fewer distant recurrences. Mechanistically, CSMD1 regulates CYP19 expression in an SNP- and drug-dependent fashion, and this regulation is different among 3 AIs: anastrozole, exemestane, and letrozole. Overexpression of CSMD1 sensitized AI-resistant cells to anastrozole but not to the other 2 AIs. The SNP in CSMD1 that was associated with increased CSMD1 and CYP19 expression levels increased anastrozole sensitivity, but not letrozole or exemestane sensitivity. Anastrozole degrades estrogen receptor α (ERα), especially in the presence of estradiol (E2). ER+ breast cancer organoids and AI- or fulvestrant-resistant breast cancer cells were more sensitive to anastrozole plus E2 than to AI alone. Our findings suggest that the CSMD1 SNP might help to predict AI response, and anastrozole plus E2 serves as a potential new therapeutic strategy for patients with AI- or fulvestrant-resistant breast cancers.

BACKGROUND: Cetuximab, an EGFR inhibitor used to treat multiple cancer types, including colon cancer, causes severe skin toxicity in 5%-20% of patients, leading to decreased quality of life and treatment delays. Our understanding of which patients have an increased risk of severe toxicities is limited. We conducted a genome-wide association study to identify germline variants predictive of cetuximab-induced severe skin toxicity. METHODS: Our study included 1,209 patients with stage III colon cancer randomized to receive cetuximab plus 5-fluorouracil and oxaliplatin as part of the NCCTG N0147 (Alliance) clinical trial. Skin toxicity outcomes were collected using the Common Toxicity Criteria for Adverse Events version 3.0. We performed genotyping, evaluating approximately 10 million genetic variants. We used logistic regression to evaluate the association of each genetic variant and severe (grade ≥ 3) skin toxicity, adjusting for age, sex, and genetic ancestry. Genome-wide significance was defined as P < 5 × 10-8. RESULTS: Participants were predominantly middle-aged white men; 20% (n = 243) experienced severe skin toxicity. Two genetic variants in the retinoic acid receptor alpha (RARA) gene were significantly associated with severe skin toxicity [OR, 3.93; 95% confidence interval (CI), 2.47-6.25; P < 7.8 × 10-9]. Functional annotations indicate these variants are in the RARA promoter. Additional significantly associated variants were identified in chromosome 2 intergenic regions. CONCLUSIONS: Identified variants could represent a potential target for risk stratification of patients with colon cancer receiving cetuximab. IMPACT: Retinoids have shown promise in the treatment of cetuximab-induced skin toxicity, so follow-up work could evaluate whether individuals with the RARA variant would benefit from retinoid therapy.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is the third leading cause of death globally. COPD patients with cachexia or weight loss have increased risk of death independent of body mass index (BMI) and lung function. We tested the hypothesis genetic variation is associated with weight loss in COPD using a genome-wide association study approach. METHODS: Participants with COPD (N = 4308) from three studies (COPDGene, ECLIPSE, and SPIROMICS) were analysed. Discovery analyses were performed in COPDGene with replication in SPIROMICS and ECLIPSE. In COPDGene, weight loss was defined as self-reported unintentional weight loss > 5% in the past year or low BMI (BMI < 20 kg/m2 ). In ECLIPSE and SPIROMICS, weight loss was calculated using available longitudinal visits. Stratified analyses were performed among African American (AA) and Non-Hispanic White (NHW) participants with COPD. Single variant and gene-based analyses were performed adjusting for confounders. Fine mapping was performed using a Bayesian approach integrating genetic association results with linkage disequilibrium and functional annotation. Significant gene networks were identified by integrating genetic regions associated with weight loss with skeletal muscle protein-protein interaction (PPI) data. RESULTS: At the single variant level, only the rs35368512 variant, intergenic to GRXCR1 and LINC02383, was associated with weight loss (odds ratio = 3.6, 95% confidence interval = 2.3-5.6, P = 3.2 × 10-8 ) among AA COPD participants in COPDGene. At the gene level in COPDGene, EFNA2 and BAIAP2 were significantly associated with weight loss in AA and NHW COPD participants, respectively. The EFNA2 association replicated among AA from SPIROMICS (P = 0.0014), whereas the BAIAP2 association replicated in NHW from ECLIPSE (P = 0.025). The EFNA2 gene encodes the membrane-bound protein ephrin-A2 involved in the regulation of developmental processes and adult tissue homeostasis such as skeletal muscle. The BAIAP2 gene encodes the insulin-responsive protein of mass 53 kD (IRSp53), a negative regulator of myogenic differentiation. Integration of the gene-based findings participants with PPI data revealed networks of genes involved in pathways such as Rho and synapse signalling. CONCLUSIONS: The EFNA2 and BAIAP2 genes were significantly associated with weight loss in COPD participants. Collectively, the integrative network analyses indicated genetic variation associated with weight loss in COPD may influence skeletal muscle regeneration and tissue remodelling.

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Walking is a simple form of exercise, widely promoted for its health benefits. Self-reported walking pace has been associated with a range of cardiorespiratory and cancer outcomes, and is a strong predictor of mortality. Here we perform a genome-wide association study of self-reported walking pace in 450,967 European ancestry UK Biobank participants. We identify 70 independent associated loci (P < 5 × 10-8), 11 of which are novel. We estimate the SNP-based heritability as 13.2% (s.e. = 0.21%), reducing to 8.9% (s.e. = 0.17%) with adjustment for body mass index. Significant genetic correlations are observed with cardiometabolic, respiratory and psychiatric traits, educational attainment and all-cause mortality. Mendelian randomization analyses suggest a potential causal link of increasing walking pace with a lower cardiometabolic risk profile. Given its low heritability and simple measurement, these findings suggest that self-reported walking pace is a pragmatic target for interventions aiming for general benefits on health.

Signalling lipids of the N-acyl ethanolamine (NAE) and ceramide (CER) classes have emerged as potential biomarkers of cardiovascular disease (CVD). We sought to establish the heritability of plasma NAEs (including the endocannabinoid anandamide) and CERs, to identify common DNA variants influencing the circulating concentrations of the heritable lipids, and assess causality of these lipids in CVD using 2-sample Mendelian randomization (2SMR). Nine NAEs and 16 CERs were analyzed in plasma samples from 999 members of 196 British Caucasian families, using targeted ultra-performance liquid chromatography with tandem mass spectrometry. All lipids were significantly heritable (h2 = 36-62%). A missense variant (rs324420) in the gene encoding the enzyme fatty acid amide hydrolase (FAAH), which degrades NAEs, associated at genome-wide association study (GWAS) significance (P < 5 × 10-8) with four NAEs (DHEA, PEA, LEA and VEA). For CERs, rs680379 in the SPTLC3 gene, which encodes a subunit of the rate-limiting enzyme in CER biosynthesis, associated with a range of species (e.g. CER[N(24)S(19)]; P = 4.82 × 10-27). We observed three novel associations between SNPs at the CD83, SGPP1 and DEGS1 loci, and plasma CER traits (P < 5 × 10-8). 2SMR in the CARDIoGRAMplusC4D cohorts (60 801 cases; 123 504 controls) and in the DIAGRAM cohort (26 488 cases; 83 964 controls), using the genetic instruments from our family-based GWAS, did not reveal association between genetically determined differences in CER levels and CVD or diabetes. Two of the novel GWAS loci, SGPP1 and DEGS1, suggested a casual association between CERs and a range of haematological phenotypes, through 2SMR in the UK Biobank, INTERVAL and UKBiLEVE cohorts (n = 110 000-350 000).

Glycosuria is a condition where glucose is detected in urine at higher concentrations than normal (i.e. not detectable). Glycosuria at some point during pregnancy has an estimated prevalence of 50% and is associated with adverse outcomes in both mothers and offspring. Little is currently known about the genetic contribution to this trait or the extent to which it overlaps with other seemingly related traits, e.g. diabetes. We performed a genome-wide association study (GWAS) for self-reported glycosuria in pregnant mothers from the Avon Longitudinal Study of Parents and Children (cases/controls = 1249/5140). We identified two loci, one of which (lead SNP = rs13337037; chromosome 16; odds ratio of glycosuria per effect allele: 1.42; 95% CI: 1.30, 1.56; P = 1.97 × 10-13) was then validated using an obstetric measure of glycosuria measured in the same cohort (227/6639). We performed a secondary GWAS in the 1986 Northern Finland Birth Cohort (NFBC1986; 747/2991) using midwife-reported glycosuria and offspring genotype as a proxy for maternal genotype. The combined results revealed evidence for a consistent effect on glycosuria at the chromosome 16 locus. In follow-up analyses, we saw little evidence of shared genetic underpinnings with the exception of urinary albumin-to-creatinine ratio (Rg = 0.64; SE = 0.22; P = 0.0042), a biomarker of kidney disease. In conclusion, we identified a genetic association with self-reported glycosuria during pregnancy, with the lead SNP located 15kB upstream of SLC5A2, a target of antidiabetic drugs. The lack of strong genetic correlation with seemingly related traits such as type 2 diabetes suggests different genetic risk factors exist for glycosuria during pregnancy.