Obesity is globally prevalent and highly heritable, but the underlying genetic factors remain largely elusive. To identify genetic loci for obesity-susceptibility, we examined associations between body mass index (BMI) and ~2.8 million SNPs in up to 123,865 individuals, with targeted follow-up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity-susceptibility loci and identified 18 new loci associated with BMI (P<5×10−8), one of which includes a copy number variant near GPRC5B. Some loci (MC4R, POMC, SH2B1, BDNF) map near key hypothalamic regulators of energy balance, and one is near GIPR, an incretin receptor. Furthermore, genes in other newly-associated loci may provide novel insights into human body weight regulation.

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence phenotype. Genome-wide association (GWA) studies have identified >600 variants associated with human traits1, but these typically explain small fractions of phenotypic variation, raising questions about the utility of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait2,3. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P=0.016), and that underlie skeletal growth defects (P<0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants, and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented amongst variants that alter amino acid structure of proteins and expression levels of nearby genes. Our data explain ∼10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to ∼16% of phenotypic variation (∼20% of heritable variation). Although additional approaches are needed to fully dissect the genetic architecture of polygenic human traits, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

We performed a meta-analysis of 14 genome-wide association studies of coronary artery disease (CAD) comprising 22,233 cases and 64,762 controls of European descent, followed by genotyping of top association signals in 60,738 additional individuals. This genomic analysis identified 13 novel loci harboring one or more SNPs that were associated with CAD at P<5×10−8 and confirmed the association of 10 of 12 previously reported CAD loci. The 13 novel loci displayed risk allele frequencies ranging from 0.13 to 0.91 and were associated with a 6 to 17 percent increase in the risk of CAD per allele. Notably, only three of the novel loci displayed significant association with traditional CAD risk factors, while the majority lie in gene regions not previously implicated in the pathogenesis of CAD. Finally, five of the novel CAD risk loci appear to have pleiotropic effects, showing strong association with various other human diseases or traits.

Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body-mass-index (up to 77,167 participants), following up 16 loci in an additional 29 studies (up to 113,636 subjects). We identified 13 novel loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1, and CPEB4 (P 1.9 × 10−9 to 1.8 × 10−40), and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex-difference 1.9 × 10−3 to 1.2 × 10−13). These findings provide evidence for multiple loci that modulate body fat distribution, independent of overall adiposity, and reveal powerful gene-by-sex interactions.

Consistent but indirect evidence has implicated genetic factors in smoking behavior1,2. We report meta-analyses of several smoking phenotypes within cohorts of the Tobacco and Genetics Consortium (n = 74,053). We also partnered with the European Network of Genetic and Genomic Epidemiology (ENGAGE) and Oxford-GlaxoSmithKline (Ox-GSK) consortia to follow up the 15 most significant regions (n > 140,000). We identified three loci associated with number of cigarettes smoked per day. The strongest association was a synonymous 15q25 SNP in the nicotinic receptor gene CHRNA3 (rs1051730[A], β = 1.03, standard error (s.e.) = 0.053, P = 2.8 × 10−73). Two 10q25 SNPs (rs1329650[G], β = 0.367, s.e. = 0.059, P = 5.7 × 10−10; and rs1028936[A], β = 0.446, s.e. = 0.074, P = 1.3 × 10−9) and one 9q13 SNP in EGLN2 (rs3733829[G], β = 0.333, s.e. = 0.058, P = 1.0 × 10−8) also exceeded genome-wide significance for cigarettes per day. For smoking initiation, eight SNPs exceeded genome-wide significance, with the strongest association at a nonsynonymous SNP in BDNF on chromosome 11 (rs6265[C], odds ratio (OR) = 1.06, 95% confidence interval (Cl) 1.04–1.08, P = 1.8 × 10−8). One SNP located near DBH on chromosome 9 (rs3025343[G], OR = 1.12, 95% Cl 1.08–1.18, P = 3.6 × 10−8) was significantly associated with smoking cessation.

Reduced glomerular filtration rate defines chronic kidney disease and is associated with cardiovascular and all-cause mortality. We conducted a meta-analysis of genome-wide association studies for estimated glomerular filtration rate (eGFR), combining data across 133,413 individuals with replication in up to 42,166 individuals. We identify 24 new and confirm 29 previously identified loci. Of these 53 loci, nineteen associate with eGFR among individuals with diabetes. Using bioinformatics, we show that identified genes at eGFR loci are enriched for expression in kidney tissues and in pathways relevant for kidney development and transmembrane transporter activity, kidney structure, and regulation of glucose metabolism. Chromatin state mapping and DNase I hypersensitivity analyses across adult tissues demonstrate preferential mapping of associated variants to regulatory regions in kidney but not extra-renal tissues. These findings suggest that genetic determinants of eGFR are mediated largely through direct effects within the kidney and highlight important cell types and biologic pathways.

SUMMARY Variability in induced pluripotent stem cell (iPSC) lines remains a concern for disease modeling and regenerative medicine. We have used RNA sequencing analysis and linear mixed models to examine the sources of gene expression variability in 317 human iPSC lines from 101 individuals. We found that ~50% of genome-wide expression variability is explained by variation across individuals and identified a set of expression quantitative trait loci that contribute to this variation. These analyses coupled with allele specific expression show that iPSCs retain a donor specific gene expression pattern. Network, pathway and key driver analyses showed that Polycomb targets contribute significantly to the non-genetic variability seen within and across individuals, highlighting this chromatin regulator as a likely source of reprogramming-based variability. Our findings therefore shed light on variation between iPSC lines and illustrate the potential for our dataset and other similar large-scale analyses to identify underlying drivers relevant to iPSC applications.

IntroductionThe natriuretic peptide system plays a crucial role in blood pressure (BP) and blood volume (BV) homeostasis. Increases in atrial stretch trigger the release of atrial natriuretic peptide (ANP) and brain natriuretic peptide (BNP), which act through natriuretic peptide receptor A (NPRA) in the kidneys and vasculature to increase natriuresis, diuresis, and vasorelaxation (1-4).Elevations in circulating ANP and BNP are cardinal features of cardiac hypertrophy and heart failure in both humans and animal models (5-8), and induction of the ANP gene is one of the most robust responses to hypertrophic stimuli in the heart (9-12). Recent results suggest that these natriuretic peptides are not simply involved in BP and BV homeostasis but that they are also directly involved in moderating the cardiac growth response to hypertrophic stimuli. For example, the NPRA system has intrinsic growth inhibitory properties in noncardiac and cardiac cells in vitro (13). In endothelial and vascular smooth muscle cells (VSMCs), ANP is antimitogenic (14-17) and in cultured neonatal rat cardiomyocytes, ANP both attenuates the growth response to adrenergic stimuli (18) through a pathway that requires extracellular signal-related protein kinase (ERK) activation (19) and induces apoptosis (20). Furthermore, NPRA knockout mice (which have greatly attenuated responses to ANP and BNP) not only have increased BP, but also display marked cardiac hypertrophy and chamber dilatation by 3 months of age (21), whereas other mutant mice with similar increases in BP do not have this hypertrophic phenotype (22). These observations have led to the hypothesis that the NPRA signaling pathway plays an important autocrine role in the prevention of cardiac hypertrophy and heart failure beyond its effects on BP and BV regulation.To test the hypothesis that the NPRA pathway directly modulates the hypertrophic response we compare the response of Npr1 +/+ and Npr1 -/-mice to a pressure overload induced by transverse aortic constriction (TAC). We have also investigated the natural progression to cardiac hypertrophy that occurs in Npr1 -/-mice in the absence of TAC and the effects of chronic (3 months) BP normalization by enalapril, furosemide, hydralazine, or losartan. We conclude that the NPRA system indeed plays a primary role in moderating cardiac hypertrophy in vivo independently of its effects on BP regulation. Mice lacking natriuretic peptide receptor A (NPRA) have marked cardiac hypertrophy and chamber dilatation disproportionate to their increased blood pressure (BP), suggesting, in support of previous in vitro data, that the NPRA system moderates the cardiac response to hypertrophic stimuli. Here, we have followed the changes in cardiac function in response to altered mechanical load on the heart of NPRA-null mice (Npr1 -/-). Chronic treatment with either enalapril, furosemide, hydralazine, or losartan were all effective in reducing and maintaining BP at normal levels without affecting heart weight/body weight. In the reverse direction, we used t...