A whole-blood transcriptome meta-analysis identifies gene expression signatures of cigarette smoking

Cigarette smoking is a leading modifiable cause of death worldwide. We hypothesized that cigarette smoking induces extensive transcriptomic changes that lead to target-organ damage and smoking-related diseases. We performed a meta-analysis of transcriptome-wide gene expression using whole blood-derived RNA from 10,233 participants of European ancestry in six cohorts (including 1421 current and 3955 former smokers) to identify associations between smoking and altered gene expression levels. At a false discovery rate (FDR) <0.1, we identified 1270 differentially expressed genes in current vs. never smokers, and 39 genes in former vs. never smokers. Expression levels of 12 genes remained elevated up to 30 years after smoking cessation, suggesting that the molecular consequence of smoking may persist for decades. Gene ontology analysis revealed enrichment of smoking-related genes for activation of platelets and lymphocytes, immune response, and apoptosis. Many of the top smoking-related differentially expressed genes, including LRRN3 and GPR15, have DNA methylation loci in promoter regions that were recently reported to be hypomethylated among smokers. By linking differential gene expression with smoking-related disease phenotypes, we demonstrated that stroke and pulmonary function show enrichment for smoking-related gene expression signatures. Mediation analysis revealed the expression of several genes (e.g. ALAS2) to be putative mediators of the associations between smoking and inflammatory biomarkers (IL6 and C-reactive protein levels). Our transcriptomic study provides potential insights into the effects of cigarette smoking on gene expression in whole blood and their relations to smoking-related diseases. The results of such analyses may highlight attractive targets for treating or preventing smoking-related health effects.

Fifteen Genetic Loci Associated With the Electrocardiographic P Wave

BACKGROUND

The P wave on an ECG is a measure of atrial electric function, and its characteristics may serve as predictors for atrial arrhythmias. Increased mean P-wave duration and P-wave terminal force traditionally have been used as markers for left atrial enlargement, and both have been associated with increased risk of atrial fibrillation. Here, we explore the genetic basis of P-wave morphology through meta-analysis of genome-wide association study results for P-wave duration and P-wave terminal force from 12 cohort studies.

METHODS AND RESULTS

We included 44 456 individuals, of which 6778 (16%) were of African ancestry. Genotyping, imputation, and genome-wide association study were performed at each study site. Summary-level results were meta-analyzed centrally using inverse-variance weighting. In meta-analyses of P-wave duration, we identified 6 significant (P<5×10^-8) novel loci and replicated a prior association with SCN10A. We identified 3 loci at SCN5A, TBX5, and CAV1/CAV2 that were jointly associated with the PR interval, PR segment, and P-wave duration. We identified 6 novel loci in meta-analysis of P-wave terminal force. Four of the identified genetic loci were significantly associated with gene expression in 329 left atrial samples. Finally, we observed that some of the loci associated with the P wave were linked to overall atrial conduction, whereas others identified distinct phases of atrial conduction.

CONCLUSIONS

We have identified 6 novel genetic loci associated with P-wave duration and 6 novel loci associated with P-wave terminal force. Future studies of these loci may aid in identifying new targets for drugs that may modify atrial conduction or treat atrial arrhythmias.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV ₁), forced vital capacity (FVC) and the ratio of FEV ₁ to FVC (FEV ₁/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 ^-7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Circulating metabolites and general cognitive ability and dementia: Evidence from 11 cohort studies

INTRODUCTION

Identifying circulating metabolites that are associated with cognition and dementia may improve our understanding of the pathogenesis of dementia and provide crucial readouts for preventive and therapeutic interventions.

METHODS

We studied 299 metabolites in relation to cognition (general cognitive ability) in two discovery cohorts (N total = 5658). Metabolites significantly associated with cognition after adjusting for multiple testing were replicated in four independent cohorts (N total = 6652), and the associations with dementia and Alzheimer's disease (N = 25,872) and lifestyle factors (N = 5168) were examined.

RESULTS

We discovered and replicated 15 metabolites associated with cognition including subfractions of high-density lipoprotein, docosahexaenoic acid, ornithine, glutamine, and glycoprotein acetyls. These associations were independent of classical risk factors including high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, triglycerides, glucose, and apolipoprotein E (APOE) genotypes. Six of the cognition-associated metabolites were related to the risk of dementia and lifestyle factors.

DISCUSSION

Circulating metabolites were consistently associated with cognition, dementia, and lifestyle factors, opening new avenues for prevention of cognitive decline and dementia.

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.

Genome-wide gene-environment interaction in depression: A systematic evaluation of candidate genes: The childhood trauma working-group of PGC-MDD

Gene by environment (GxE) interaction studies have investigated the influence of a number of candidate genes and variants for major depressive disorder (MDD) on the association between childhood trauma and MDD. Most of these studies are hypothesis driven and investigate only a limited number of SNPs in relevant pathways using differing methodological approaches. Here (1) we identified 27 genes and 268 SNPs previously associated with MDD or with GxE interaction in MDD and (2) analyzed their impact on GxE in MDD using a common approach in 3944 subjects of European ancestry from the Psychiatric Genomics Consortium who had completed the Childhood Trauma Questionnaire. (3) We subsequently used the genome-wide SNP data for a genome-wide case-control GxE model and GxE case-only analyses testing for an enrichment of associated SNPs. No genome-wide significant hits and no consistency among the signals of the different analytic approaches could be observed. This is the largest study for systematic GxE interaction analysis in MDD in subjects of European ancestry to date. Most of the known candidate genes/variants could not be supported. Thus, their impact on GxE interaction in MDD may be questionable. Our results underscore the need for larger samples, more extensive assessment of environmental exposures, and greater efforts to investigate new methodological approaches in GxE models for MDD.

Genetic Association of Major Depression With Atypical Features and Obesity-Related Immunometabolic Dysregulations

IMPORTANCE

The association between major depressive disorder (MDD) and obesity may stem from shared immunometabolic mechanisms particularly evident in MDD with atypical features, characterized by increased appetite and/or weight (A/W) during an active episode.

OBJECTIVE

To determine whether subgroups of patients with MDD stratified according to the A/W criterion had a different degree of genetic overlap with obesity-related traits (body mass index [BMI] and levels of C-reactive protein [CRP] and leptin).

DESIGN, SETTING, AND PATIENTS

This multicenter study assembled genome-wide genotypic and phenotypic measures from 14 data sets of the Psychiatric Genomics Consortium. Data sets were drawn from case-control, cohort, and population-based studies, including 26 628 participants with established psychiatric diagnoses and genome-wide genotype data. Data on BMI were available for 15 237 participants. Data were retrieved and analyzed from September 28, 2015, through May 20, 2017.

MAIN OUTCOMES AND MEASURES

Lifetime DSM-IV MDD was diagnosed using structured diagnostic instruments. Patients with MDD were stratified into subgroups according to change in the DSM-IV A/W symptoms as decreased or increased.

RESULTS

Data included 11 837 participants with MDD and 14 791 control individuals, for a total of 26 628 participants (59.1% female and 40.9% male). Among participants with MDD, 5347 (45.2%) were classified in the decreased A/W and 1871 (15.8%) in the increased A/W subgroups. Common genetic variants explained approximately 10% of the heritability in the 2 subgroups. The increased A/W subgroup showed a strong and positive genetic correlation (SE) with BMI (0.53 [0.15]; P = 6.3 × 10-4), whereas the decreased A/W subgroup showed an inverse correlation (-0.28 [0.14]; P = .06). Furthermore, the decreased A/W subgroup had a higher polygenic risk for increased BMI (odds ratio [OR], 1.18; 95% CI, 1.12-1.25; P = 1.6 × 10-10) and levels of CRP (OR, 1.08; 95% CI, 1.02-1.13; P = 7.3 × 10-3) and leptin (OR, 1.09; 95% CI, 1.06-1.12; P = 1.7 × 10-3).

CONCLUSIONS AND RELEVANCE

The phenotypic associations between atypical depressive symptoms and obesity-related traits may arise from shared pathophysiologic mechanisms in patients with MDD. Development of treatments effectively targeting immunometabolic dysregulations may benefit patients with depression and obesity, both syndromes with important disability.

Investigation of Naturally Occurring Single-Nucleotide Variants in Human TAAR1

Activation of trace amine-associated receptor 1 (TAAR1) in endocrine pancreas is involved in weight regulation and glucose homeostasis. The purpose of this study was the identification and characterization of potential TAAR1 variants in patients with overweight/obesity and disturbed glucose homeostasis. Screening for TAAR1 variants was performed in 314 obese or overweight patients with impaired insulin secretion. The detected variants were functionally characterized concerning TAAR1 cell surface expression and signaling properties and their allele frequencies were determined in the population-based Study of Health in Pomerania (SHIP). Three heterozygous carriers of the single nucleotide missense variants p.Arg23Cys (R23C, rs8192618), p.Ser49Leu (S49L, rs140960896), and p.Ille171Leu (I171L, rs200795344) were detected in the patient cohort. While p.Ser49Leu and p.Ille171Leu were found in obese/overweight patients with slightly impaired glucose homeostasis, p.Arg23Cys was identified in a patient with a complete loss of insulin production. Functional in vitro characterization revealed a like wild-type function for I171L, partial loss of function for S49L and a complete loss of function for R23C. The frequency of the R23C variant in 2018 non-diabetic control individuals aged 60 years and older in the general population-based SHIP cohort was lower than in the analyzed patient sample. Both variants are rare in the general population indicating a recent origin in the general gene pool and/or the consequence of pronounced purifying selection, in line with the obvious detrimental effect of the mutations. In conclusion, our study provides hints for the existence of naturally occurring TAAR1 variants with potential relevance for weight regulation and glucose homeostasis.

A genome-wide association study identifies nucleotide variants at SIGLEC5 and DEFA1A3 as risk loci for periodontitis

Periodontitis is one of the most common inflammatory diseases, with a prevalence of 11% worldwide for the severe forms and an estimated heritability of 50%. The disease is characterized by destruction of the alveolar bone due to an aberrant host inflammatory response to a dysbiotic oral microbiome. Previous genome-wide association studies (GWAS) have reported several suggestive susceptibility loci. Here, we conducted a GWAS using a German and Dutch case-control sample of aggressive periodontitis (AgP, 896 cases, 7,104 controls), a rare but highly severe and early-onset form of periodontitis, validated the associations in a German sample of severe forms of the more moderate phenotype chronic periodontitis (CP) (993 cases, 1,419 controls). Positive findings were replicated in a Turkish sample of AgP (223 cases, 564 controls). A locus at SIGLEC5 (sialic acid binding Ig-like lectin 5) and a chromosomal region downstream of the DEFA1A3 locus (defensin alpha 1-3) showed association with both disease phenotypes and were associated with periodontitis at a genome-wide significance level in the pooled samples, with P = 1.09E-08 (rs4284742,-G; OR = 1.34, 95% CI = 1.21-1.48) and P = 5.48E-10 (rs2738058,-T; OR = 1.28, 95% CI = 1.18-1.38), respectively. SIGLEC5 is expressed in various myeloid immune cells and classified as an inhibitory receptor with the potential to mediate tyrosine phosphatases SHP-1/-2 dependent signaling. Alpha defensins are antimicrobial peptides with expression in neutrophils and mucosal surfaces and a role in phagocyte-mediated host defense. This study identifies the first shared genetic risk loci of AgP and CP with genome-wide significance and highlights the role of innate and adaptive immunity in the etiology of periodontitis.

An Analysis of Two Genome-wide Association Meta-analyses Identifies a New Locus for Broad Depression Phenotype

BACKGROUND

The genetics of depression has been explored in genome-wide association studies that focused on either major depressive disorder or depressive symptoms with mostly negative findings. A broad depression phenotype including both phenotypes has not been tested previously using a genome-wide association approach. We aimed to identify genetic polymorphisms significantly associated with a broad phenotype from depressive symptoms to major depressive disorder.

METHODS

We analyzed two prior studies of 70,017 participants of European ancestry from general and clinical populations in the discovery stage. We performed a replication meta-analysis of 28,328 participants. Single nucleotide polymorphism (SNP)-based heritability and genetic correlations were calculated using linkage disequilibrium score regression. Discovery and replication analyses were performed using a p-value-based meta-analysis. Lifetime major depressive disorder and depressive symptom scores were used as the outcome measures.

RESULTS

The SNP-based heritability of major depressive disorder was 0.21 (SE = 0.02), the SNP-based heritability of depressive symptoms was 0.04 (SE = 0.01), and their genetic correlation was 1.001 (SE = 0.2). We found one genome-wide significant locus related to the broad depression phenotype (rs9825823, chromosome 3: 61,082,153, p = 8.2 × 10^-9) located in an intron of the FHIT gene. We replicated this SNP in independent samples (p = .02) and the overall meta-analysis of the discovery and replication cohorts (1.0 × 10^-9).

CONCLUSIONS

This large study identified a new locus for depression. Our results support a continuum between depressive symptoms and major depressive disorder. A phenotypically more inclusive approach may help to achieve the large sample sizes needed to detect susceptibility loci for depression.

Novel Blood Pressure Locus and Gene Discovery Using Genome-Wide Association Study and Expression Data Sets From Blood and the Kidney

Elevated blood pressure is a major risk factor for cardiovascular disease and has a substantial genetic contribution. Genetic variation influencing blood pressure has the potential to identify new pharmacological targets for the treatment of hypertension. To discover additional novel blood pressure loci, we used 1000 Genomes Project-based imputation in 150 134 European ancestry individuals and sought significant evidence for independent replication in a further 228 245 individuals. We report 6 new signals of association in or near HSPB7, TNXB, LRP12, LOC283335, SEPT9, and AKT2, and provide new replication evidence for a further 2 signals in EBF2 and NFKBIA. Combining large whole-blood gene expression resources totaling 12 607 individuals, we investigated all novel and previously reported signals and identified 48 genes with evidence for involvement in blood pressure regulation that are significant in multiple resources. Three novel kidney-specific signals were also detected. These robustly implicated genes may provide new leads for therapeutic innovation.

Transcriptome-Wide Analysis Identifies Novel Associations With Blood Pressure

Hypertension represents a major cardiovascular risk factor. The pathophysiology of increased blood pressure (BP) is not yet completely understood. Transcriptome profiling offers possibilities to uncover genetics effects on BP. Based on 2 populations including 2549 individuals, a meta-analyses of monocytic transcriptome-wide profiles were performed to identify transcripts associated with BP. Replication was performed in 2 independent studies of whole-blood transcriptome data including 1990 individuals. For identified candidate genes, a direct link between long-term changes in BP and gene expression over time and by treatment with BP-lowering therapy was assessed. The predictive value of protein levels encoded by candidate genes for subsequent cardiovascular disease was investigated. Eight transcripts (CRIP1, MYADM, TIPARP, TSC22D3, CEBPA, F12, LMNA, and TPPP3) were identified jointly accounting for up to 13% (95% confidence interval, 8.7-16.2) of BP variability. Changes in CRIP1, MYADM, TIPARP, LMNA, TSC22D3, CEBPA, and TPPP3 expression associated with BP changes-among these, CRIP1 gene expression was additionally correlated to measures of cardiac hypertrophy. Assessment of circulating CRIP1 (cystein-rich protein 1) levels as biomarkers showed a strong association with increased risk for incident stroke (hazard ratio, 1.06; 95% confidence interval, 1.03-1.09; P=5.0×10^-5). Our comprehensive analysis of global gene expression highlights 8 novel transcripts significantly associated with BP, providing a link between gene expression and BP. Translational approaches further established evidence for the potential use of CRIP1 as emerging disease-related biomarker.

Cystatin C and Cardiovascular Disease: A Mendelian Randomization Study

BACKGROUND

Epidemiological studies show that high circulating cystatin C is associated with risk of cardiovascular disease (CVD), independent of creatinine-based renal function measurements. It is unclear whether this relationship is causal, arises from residual confounding, and/or is a consequence of reverse causation.

OBJECTIVES

The aim of this study was to use Mendelian randomization to investigate whether cystatin C is causally related to CVD in the general population.

METHODS

We incorporated participant data from 16 prospective cohorts (n = 76,481) with 37,126 measures of cystatin C and added genetic data from 43 studies (n = 252,216) with 63,292 CVD events. We used the common variant rs911119 in CST3 as an instrumental variable to investigate the causal role of cystatin C in CVD, including coronary heart disease, ischemic stroke, and heart failure.

RESULTS

Cystatin C concentrations were associated with CVD risk after adjusting for age, sex, and traditional risk factors (relative risk: 1.82 per doubling of cystatin C; 95% confidence interval [CI]: 1.56 to 2.13; p = 2.12 × 10⁻¹⁴). The minor allele of rs911119 was associated with decreased serum cystatin C (6.13% per allele; 95% CI: 5.75 to 6.50; p = 5.95 × 10⁻²¹¹), explaining 2.8% of the observed variation in cystatin C. Mendelian randomization analysis did not provide evidence for a causal role of cystatin C, with a causal relative risk for CVD of 1.00 per doubling cystatin C (95% CI: 0.82 to 1.22; p = 0.994), which was statistically different from the observational estimate (p = 1.6 × 10⁻⁵). A causal effect of cystatin C was not detected for any individual component of CVD.

CONCLUSIONS

Mendelian randomization analyses did not support a causal role of cystatin C in the etiology of CVD. As such, therapeutics targeted at lowering circulating cystatin C are unlikely to be effective in preventing CVD.

Genetic loci associated with heart rate variability and their effects on cardiac disease risk

Reduced cardiac vagal control reflected in low heart rate variability (HRV) is associated with greater risks for cardiac morbidity and mortality. In two-stage meta-analyses of genome-wide association studies for three HRV traits in up to 53,174 individuals of European ancestry, we detect 17 genome-wide significant SNPs in eight loci. HRV SNPs tag non-synonymous SNPs (in NDUFA11 and KIAA1755), expression quantitative trait loci (eQTLs) (influencing GNG11, RGS6 and NEO1), or are located in genes preferentially expressed in the sinoatrial node (GNG11, RGS6 and HCN4). Genetic risk scores account for 0.9 to 2.6% of the HRV variance. Significant genetic correlation is found for HRV with heart rate (-0.74<r_g<-0.55) and blood pressure (-0.35<r_g<-0.20). These findings provide clinically relevant biological insight into heritable variation in vagal heart rhythm regulation, with a key role for genetic variants (GNG11, RGS6) that influence G-protein heterotrimer action in GIRK-channel induced pacemaker membrane hyperpolarization.

Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk

The timing of puberty is a highly polygenic childhood trait that is epidemiologically associated with various adult diseases. Using 1000 Genomes Project-imputed genotype data in up to ∼370,000 women, we identify 389 independent signals (P < 5 × 10-8) for age at menarche, a milestone in female pubertal development. In Icelandic data, these signals explain ∼7.4% of the population variance in age at menarche, corresponding to ∼25% of the estimated heritability. We implicate ∼250 genes via coding variation or associated expression, demonstrating significant enrichment in neural tissues. Rare variants near the imprinted genes MKRN3 and DLK1 were identified, exhibiting large effects when paternally inherited. Mendelian randomization analyses suggest causal inverse associations, independent of body mass index (BMI), between puberty timing and risks for breast and endometrial cancers in women and prostate cancer in men. In aggregate, our findings highlight the complexity of the genetic regulation of puberty timing and support causal links with cancer susceptibility.

Large-scale genome-wide analysis identifies genetic variants associated with cardiac structure and function

BACKGROUND

Understanding the genetic architecture of cardiac structure and function may help to prevent and treat heart disease. This investigation sought to identify common genetic variations associated with inter-individual variability in cardiac structure and function.

METHODS

A GWAS meta-analysis of echocardiographic traits was performed, including 46,533 individuals from 30 studies (EchoGen consortium). The analysis included 16 traits of left ventricular (LV) structure, and systolic and diastolic function.

RESULTS

The discovery analysis included 21 cohorts for structural and systolic function traits (n = 32,212) and 17 cohorts for diastolic function traits (n = 21,852). Replication was performed in 5 cohorts (n = 14,321) and 6 cohorts (n = 16,308), respectively. Besides 5 previously reported loci, the combined meta-analysis identified 10 additional genome-wide significant SNPs: rs12541595 near MTSS1 and rs10774625 in ATXN2 for LV end-diastolic internal dimension; rs806322 near KCNRG, rs4765663 in CACNA1C, rs6702619 near PALMD, rs7127129 in TMEM16A, rs11207426 near FGGY, rs17608766 in GOSR2, and rs17696696 in CFDP1 for aortic root diameter; and rs12440869 in IQCH for Doppler transmitral A-wave peak velocity. Findings were in part validated in other cohorts and in GWAS of related disease traits. The genetic loci showed associations with putative signaling pathways, and with gene expression in whole blood, monocytes, and myocardial tissue.

CONCLUSION

The additional genetic loci identified in this large meta-analysis of cardiac structure and function provide insights into the underlying genetic architecture of cardiac structure and warrant follow-up in future functional studies.

FUNDING

For detailed information per study, see Acknowledgments.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Genome-wide physical activity interactions in adiposity - A meta-analysis of 200,452 adults

Physical activity (PA) may modify the genetic effects that give rise to increased risk of obesity. To identify adiposity loci whose effects are modified by PA, we performed genome-wide interaction meta-analyses of BMI and BMI-adjusted waist circumference and waist-hip ratio from up to 200,452 adults of European (n = 180,423) or other ancestry (n = 20,029). We standardized PA by categorizing it into a dichotomous variable where, on average, 23% of participants were categorized as inactive and 77% as physically active. While we replicate the interaction with PA for the strongest known obesity-risk locus in the FTO gene, of which the effect is attenuated by ~30% in physically active individuals compared to inactive individuals, we do not identify additional loci that are sensitive to PA. In additional genome-wide meta-analyses adjusting for PA and interaction with PA, we identify 11 novel adiposity loci, suggesting that accounting for PA or other environmental factors that contribute to variation in adiposity may facilitate gene discovery.

NFAT5 and SLC4A10 Loci Associate with Plasma Osmolality

Disorders of water balance, an excess or deficit of total body water relative to body electrolyte content, are common and ascertained by plasma hypo- or hypernatremia, respectively. We performed a two-stage genome-wide association study meta-analysis on plasma sodium concentration in 45,889 individuals of European descent (stage 1 discovery) and 17,637 additional individuals of European descent (stage 2 replication), and a transethnic meta-analysis of replicated single-nucleotide polymorphisms in 79,506 individuals (63,526 individuals of European descent, 8765 individuals of Asian Indian descent, and 7215 individuals of African descent). In stage 1, we identified eight loci associated with plasma sodium concentration at P<5.0 × 10^-6 Of these, rs9980 at NFAT5 replicated in stage 2 meta-analysis (P=3.1 × 10^-5), with combined stages 1 and 2 genome-wide significance of P=5.6 × 10^-10 Transethnic meta-analysis further supported the association at rs9980 (P=5.9 × 10^-12). Additionally, rs16846053 at SLC4A10 showed nominally, but not genome-wide, significant association in combined stages 1 and 2 meta-analysis (P=6.7 × 10^-8). NFAT5 encodes a ubiquitously expressed transcription factor that coordinates the intracellular response to hypertonic stress but was not previously implicated in the regulation of systemic water balance. SLC4A10 encodes a sodium bicarbonate transporter with a brain-restricted expression pattern, and variant rs16846053 affects a putative intronic NFAT5 DNA binding motif. The lead variants for NFAT5 and SLC4A10 are cis expression quantitative trait loci in tissues of the central nervous system and relevant to transcriptional regulation. Thus, genetic variation in NFAT5 and SLC4A10 expression and function in the central nervous system may affect the regulation of systemic water balance.

Genetic effects influencing risk for major depressive disorder in China and Europe

Major depressive disorder (MDD) is a common, complex psychiatric disorder and a leading cause of disability worldwide. Despite twin studies indicating its modest heritability (~30-40%), extensive heterogeneity and a complex genetic architecture have complicated efforts to detect associated genetic risk variants. We combined single-nucleotide polymorphism (SNP) summary statistics from the CONVERGE and PGC studies of MDD, representing 10 502 Chinese (5282 cases and 5220 controls) and 18 663 European (9447 cases and 9215 controls) subjects. We determined the fraction of SNPs displaying consistent directions of effect, assessed the significance of polygenic risk scores and estimated the genetic correlation of MDD across ancestries. Subsequent trans-ancestry meta-analyses combined SNP-level evidence of association. Sign tests and polygenic score profiling weakly support an overlap of SNP effects between East Asian and European populations. We estimated the trans-ancestry genetic correlation of lifetime MDD as 0.33; female-only and recurrent MDD yielded estimates of 0.40 and 0.41, respectively. Common variants downstream of GPHN achieved genome-wide significance by Bayesian trans-ancestry meta-analysis (rs9323497; log10 Bayes Factor=8.08) but failed to replicate in an independent European sample (P=0.911). Gene-set enrichment analyses indicate enrichment of genes involved in neuronal development and axonal trafficking. We successfully demonstrate a partially shared polygenic basis of MDD in East Asian and European populations. Taken together, these findings support a complex etiology for MDD and possible population differences in predisposing genetic factors, with important implications for future genetic studies.

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.

PCSK9 genetic variants and risk of type 2 diabetes: a mendelian randomisation study

BACKGROUND

Statin treatment and variants in the gene encoding HMG-CoA reductase are associated with reductions in both the concentration of LDL cholesterol and the risk of coronary heart disease, but also with modest hyperglycaemia, increased bodyweight, and modestly increased risk of type 2 diabetes, which in no way offsets their substantial benefits. We sought to investigate the associations of LDL cholesterol-lowering PCSK9 variants with type 2 diabetes and related biomarkers to gauge the likely effects of PCSK9 inhibitors on diabetes risk.

METHODS

In this mendelian randomisation study, we used data from cohort studies, randomised controlled trials, case control studies, and genetic consortia to estimate associations of PCSK9 genetic variants with LDL cholesterol, fasting blood glucose, HbA1c, fasting insulin, bodyweight, waist-to-hip ratio, BMI, and risk of type 2 diabetes, using a standardised analysis plan, meta-analyses, and weighted gene-centric scores.

FINDINGS

Data were available for more than 550 000 individuals and 51 623 cases of type 2 diabetes. Combined analyses of four independent PCSK9 variants (rs11583680, rs11591147, rs2479409, and rs11206510) scaled to 1 mmol/L lower LDL cholesterol showed associations with increased fasting glucose (0·09 mmol/L, 95% CI 0·02 to 0·15), bodyweight (1·03 kg, 0·24 to 1·82), waist-to-hip ratio (0·006, 0·003 to 0·010), and an odds ratio for type diabetes of 1·29 (1·11 to 1·50). Based on the collected data, we did not identify associations with HbA1c (0·03%, -0·01 to 0·08), fasting insulin (0·00%, -0·06 to 0·07), and BMI (0·11 kg/m2, -0·09 to 0·30).

INTERPRETATION

PCSK9 variants associated with lower LDL cholesterol were also associated with circulating higher fasting glucose concentration, bodyweight, and waist-to-hip ratio, and an increased risk of type 2 diabetes. In trials of PCSK9 inhibitor drugs, investigators should carefully assess these safety outcomes and quantify the risks and benefits of PCSK9 inhibitor treatment, as was previously done for statins.

FUNDING

British Heart Foundation, and University College London Hospitals NHS Foundation Trust (UCLH) National Institute for Health Research (NIHR) Biomedical Research Centre.