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.

Serum liver enzyme concentrations are the most frequently-used laboratory markers of liver disease, a major cause of mortality. We conduct a meta-analysis of genome-wide association studies of liver enzymes from UK BioBank and BioBank Japan. We identified 160 previously-unreported independent alanine aminotransferase, 190 aspartate aminotransferase, and 199 alkaline phosphatase genome-wide significant associations, with some affecting multiple different enzymes. Associated variants implicate genes that demonstrate diverse liver cell type expression and promote a range of metabolic and liver diseases. These findings provide insight into the pathophysiology of liver and other metabolic diseases that are associated with serum liver enzyme concentrations.

Evidence for a role for vitamin D in non-alcoholic fatty liver disease (NAFLD) pathogenesis is conflicting. As Mendelian randomisation (MR) avoids many limitations of conventional observational studies, this two-sample bidirectional MR analysis was conducted to determine the following: (i) whether genetically predicted 25-hydroxyvitamin D [25(OH)D] levels are a risk factor for NAFLD, and (ii) whether genetic risk for NAFLD influences 25(OH)D levels. Single-nucleotide polymorphisms (SNPs) associated with serum 25(OH)D levels were obtained from the European ancestry-derived SUNLIGHT consortium. SNPs associated with NAFLD or NASH (p-value < 1 × 10-5) were extracted from previous studies and supplemented by genome-wide association studies (GWASs) performed in the UK Biobank. These GWASs were done both without (primary analysis) and with (sensitivity analysis) the population-level exclusion of other liver diseases (e.g., alcoholic liver diseases, toxic liver diseases, viral hepatitis, etc.). Subsequently, MR analyses were performed to obtain effect estimates using inverse variance weighted (IVW) random effect models. Cochran's Q statistic, MR-Egger regression intercept, MR pleiotropy residual sum and outlier (MR-PRESSO) analyses were used to assess pleiotropy. No causal association of genetically predicted serum 25(OH)D (per standard deviation increase) with risk of NAFLD was identified in either the primary analysis: n = 2757 cases, n = 460,161 controls, odds ratio (95% confidence interval): 0.95 (0.76, -1.18), p = 0.614; or the sensitivity analysis. Reciprocally, no causal association was identified between the genetic risk of NAFLD and serum 25(OH)D levels, OR = 1.00 (0.99, 1.02, p = 0.665). In conclusion, this MR analysis found no evidence of an association between serum 25(OH)D levels and NAFLD in a large European cohort.

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.

Nonalcoholic fatty liver disease (NAFLD) clusters in families, but the only known common genetic variants influencing risk are near PNPLA3. We sought to identify additional genetic variants influencing NAFLD using genome-wide association (GWA) analysis of computed tomography (CT) measured hepatic steatosis, a non-invasive measure of NAFLD, in large population based samples. Using variance components methods, we show that CT hepatic steatosis is heritable (∼26%-27%) in family-based Amish, Family Heart, and Framingham Heart Studies (n = 880 to 3,070). By carrying out a fixed-effects meta-analysis of genome-wide association (GWA) results between CT hepatic steatosis and ∼2.4 million imputed or genotyped SNPs in 7,176 individuals from the Old Order Amish, Age, Gene/Environment Susceptibility-Reykjavik study (AGES), Family Heart, and Framingham Heart Studies, we identify variants associated at genome-wide significant levels (p<5×10(-8)) in or near PNPLA3, NCAN, and PPP1R3B. We genotype these and 42 other top CT hepatic steatosis-associated SNPs in 592 subjects with biopsy-proven NAFLD from the NASH Clinical Research Network (NASH CRN). In comparisons with 1,405 healthy controls from the Myocardial Genetics Consortium (MIGen), we observe significant associations with histologic NAFLD at variants in or near NCAN, GCKR, LYPLAL1, and PNPLA3, but not PPP1R3B. Variants at these five loci exhibit distinct patterns of association with serum lipids, as well as glycemic and anthropometric traits. We identify common genetic variants influencing CT-assessed steatosis and risk of NAFLD. Hepatic steatosis associated variants are not uniformly associated with NASH/fibrosis or result in abnormalities in serum lipids or glycemic and anthropometric traits, suggesting genetic heterogeneity in the pathways influencing these traits.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Diverticular disease is common and has a high morbidity. Treatments are limited owing to the poor understanding of its pathophysiology. Here, to elucidate its etiology, we performed a genome-wide association study of diverticular disease (27,444 cases; 382,284 controls) from the UK Biobank and tested for replication in the Michigan Genomics Initiative (2,572 cases; 28,649 controls). We identified 42 loci associated with diverticular disease; 39 of these loci are novel. Using data-driven expression-prioritized integration for complex traits (DEPICT), we show that genes in these associated regions are significantly enriched for expression in mesenchymal stem cells and multiple connective tissue cell types and are co-expressed with genes that have a role in vascular and mesenchymal biology. Genes in these associated loci have roles in immunity, extracellular matrix biology, cell adhesion, membrane transport and intestinal motility. Phenome-wide association analysis of the 42 variants shows a common etiology of diverticular disease with obesity and hernia. These analyses shed light on the genomic landscape of diverticular disease.

Genome-wide association studies (GWAS) have identified 36 loci associated with body mass index (BMI), predominantly in populations of European ancestry. We conducted a meta-analysis to examine the association of >3.2 million SNPs with BMI in 39,144 men and women of African ancestry and followed up the most significant associations in an additional 32,268 individuals of African ancestry. We identified one new locus at 5q33 (GALNT10, rs7708584, P = 3.4 × 10(-11)) and another at 7p15 when we included data from the GIANT consortium (MIR148A-NFE2L3, rs10261878, P = 1.2 × 10(-10)). We also found suggestive evidence of an association at a third locus at 6q16 in the African-ancestry sample (KLHL32, rs974417, P = 6.9 × 10(-8)). Thirty-two of the 36 previously established BMI variants showed directionally consistent effect estimates in our GWAS (binomial P = 9.7 × 10(-7)), five of which reached genome-wide significance. These findings provide strong support for shared BMI loci across populations, as well as for the utility of studying ancestrally diverse populations.

Nonalcoholic fatty liver disease (NAFLD) is prevalent worldwide. NAFLD is associated with elevated serum triglycerides (TG), low-density lipoprotein cholesterol (LDL), and reduced high-density lipoprotein cholesterol (HDL). Both NAFLD and blood lipid levels are genetically influenced and may share a common genetic etiology. We used genome-wide association studies (GWAS)-ranked genes and gene-set enrichment analysis to identify pathways that affect serum lipids and NAFLD. We identified credible genes in these pathways and characterized missense variants in these for effects on serum traits. We used MAGENTA to identify 58 enriched pathways from publicly available TG, LDL, and HDL GWAS (n = 99,000). Three of these pathways were also enriched for associations with European-ancestry NAFLD GWAS (n = 7176). One pathway, farnesoid X receptor (FXR)/retinoid X receptor (RXR) activation, was replicated for association in an African-ancestry NAFLD GWAS (n = 3214) and plays a role in serum lipids and NAFLD. Credible genes (proteins) in FXR/RXR activation include those associated with cholesterol/bile/bilirubin transport/absorption (ABCC2 (MRP2) [ATP binding cassette subfamily C member (multidrug resistance-associated protein 2)], ABCG5, ABCG8 [ATP-binding cassette (ABC) transporters G5 and G8], APOB (APOB) [apolipoprotein B], FABP6 (ILBP) [fatty acid binding protein 6 (ileal lipid-binding protein)], MTTP (MTP) [microsomal triglyceride transfer protein], SLC4A2 (AE2) [solute carrier family 4 member 2 (anion exchange protein 2)]), nuclear hormone-mediated control of metabolism (NR0B2 (SHP) [nuclear receptor subfamily 0 group B member 2 (small heterodimer partner)], NR1H4 (FXR) [nuclear receptor subfamily 1 group H member 4 (FXR)], PPARA (PPAR) [peroxisome proliferator activated receptor alpha], FOXO1 (FOXO1A) [forkhead box O1]), or other pathways (FETUB (FETUB) [fetuin B]). Missense variants in ABCC2 (MRP2), ABCG5 (ABCG5), ABCG8 (ABCG8), APOB (APOB), MTTP (MTP), NR0B2 (SHP), NR1H4 (FXR), and PPARA (PPAR) that associate with serum LDL levels also associate with serum liver function tests in UK Biobank. Conclusion: Genetic variants in NR1H4 (FXR) that protect against liver steatosis increase serum LDL cholesterol while variants in other members of the family have congruent effects on these traits. Human genetic pathway enrichment analysis can help guide therapeutic development by identifying effective targets for NAFLD/serum lipid manipulation while minimizing side effects. In addition, missense variants could be used in companion diagnostics to determine their influence on drug effectiveness.