Summary Height is a highly heritable, classic polygenic trait with ∼700 common associated variants identified so far through genome-wide association studies. Here, we report 83 height-associated coding variants with lower minor allele frequencies (range of 0.1-4.8%) and effects of up to 2 cm/allele (e.g. in IHH, STC2, AR and CRISPLD2), >10 times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (+1-2 cm/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 mutated in monogenic growth disorders and highlight new biological candidates (e.g. ADAMTS3, IL11RA, NOX4) and pathways (e.g. proteoglycan/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.

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, non-coding variants from which pinpointing causal genes remains challenging. Here, we combined data from 718,734 individuals to discover rare and low-frequency (MAF<5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which eight in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2, ZNF169) newly implicated in human obesity, two (MC4R, KSR2) previously observed in extreme obesity, and two variants in GIPR. Effect sizes of rare variants are ~10 times larger than of common variants, with the largest effect observed in carriers of an MC4R stop-codon (p.Tyr35Ter, MAF=0.01%), weighing ~7kg more than non-carriers. Pathway analyses confirmed enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically-supported therapeutic targets to treat obesity.

Hematological traits are important clinical parameters. To test the role of rare and low-frequency coding variants on hematological traits, we analyzed hemoglobin, hematocrit, white blood cell (WBC) and platelet count in 31,340 individuals genotyped on an exome array. We identified several missense variants of CXCR2 associated with reduced WBC count (gene-based P=2.6×10−13). In a separate family-based re-sequencing study, we identified a novel loss-of-function CXCR2 frameshift mutation in a pedigree with congenital neutropenia that abolished ligand-induced CXCR2 signal transduction and chemotaxis. We also identified novel missense or splice site variants in key hematopoiesis regulators (EPO, TRF2, HBB, TUBB1, SH2B3) associated with blood cell traits. Finally, we were able to detect associations between the rare somatic JAK2 p.Val617Phe mutation and platelet count (P=3.9×10−22) as well as hemoglobin (P=0.002), hematocrit (P=9.5×10−7) and WBC (P=3.1×10−5). In conclusion, exome arrays complement GWAS in identifying new variants that contribute to complex human traits.

Smoking is a major heritable and modifiable risk factor for many diseases, including cancer, common respiratory disorders and cardiovascular diseases. Fourteen genetic loci have previously been associated with smoking behaviour-related traits. We tested up to 235,116 single nucleotide variants (SNVs) on the exome-array for association with smoking initiation, cigarettes per day, pack-years, and smoking cessation in a fixed effects meta-analysis of up to 61 studies (up to 346,813 participants). In a subset of 112,811 participants, a further one million SNVs were also genotyped and tested for association with the four smoking behaviour traits. SNV-trait associations with P < 5 × 10 −8 in either analysis were taken forward for replication in up to 275,596 independent participants from UK Biobank. Lastly, a meta-analysis of the discovery and replication studies was performed. Sixteen SNVs were associated with at least one of the smoking behaviour traits (P < 5 × 10 −8) in the discovery samples. Ten novel SNVs, including rs12616219 near TMEM182, were followed-up and five of them (rs462779 in REV3L, rs12780116 in CNNM2, rs1190736 in GPR101, rs11539157 in PJA1, and rs12616219 near TMEM182) replicated at a Bonferroni significance threshold (P < 4.5 × 10 −3) with consistent direction of effect. A further 35 SNVs were associated with smoking behaviour traits in the discovery plus replication meta-analysis (up to 622,409 participants) including a rare SNV, rs150493199, in CCDC141 and two low-frequency SNVs in CEP350 and HDGFRP2. Functional follow-up implied that decreased expression of REV3L may lower the probability of smoking initiation. The novel loci will facilitate understanding the genetic aetiology of smoking behaviour and may lead to the identification of potential drug targets for smoking prevention and/or cessation.

BackgroundEpigenetic mechanisms may be involved in the regulation of genes found to be differentially expressed in the visceral adipose tissue (VAT) of severely obese subjects with (MetS+) versus without (MetS-) metabolic syndrome (MetS). Long interspersed nuclear element 1 (LINE-1) elements DNA methylation levels (%meth) in blood, a marker of global DNA methylation, have recently been associated with fasting glucose, blood lipids, heart diseases and stroke.AimTo test whether LINE-1%meth levels in VAT are associated with MetS phenotypes and whether they can predict MetS risk in severely obese individuals.MethodsDNA was extracted from VAT of 34 men (MetS-: n = 14, MetS+: n = 20) and 152 premenopausal women (MetS-: n = 84; MetS+: n = 68) undergoing biliopancreatic diversion for the treatment of obesity. LINE-1%meth levels were assessed by pyrosequencing of sodium bisulfite-treated DNA.ResultsThe mean LINE-1%meth in VAT was of 75.8% (SD = 3.0%). Multiple linear regression analyses revealed that LINE-1%meth was negatively associated with fasting glucose levels (β = -0.04; P = 0.03), diastolic blood pressure (β =  -0.65; P = 0.03) and MetS status (β = -0.04; P = 0.004) after adjustments for the effects of age, sex, waist circumference (except for MetS status) and smoking. While dividing subjects into quartiles based on their LINE-1%meth (Q1 to Q4: lower %meth to higher %meth levels), greater risk were observed in the first (Q1: odds ratio (OR) = 4.37, P = 0.004) and the second (Q2: OR = 4.76, P = 0.002) quartiles compared to Q4 (1.00) when adjusting for age, sex and smoking.ConclusionsThese results suggest that lower global DNA methylation, assessed by LINE-1 repetitive elements methylation analysis, would be associated with a greater risk for MetS in the presence of obesity.