Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits.

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

Non-Alcoholic Fatty Liver Disease and Vitamin D in the UK Biobank: A Two-Sample Bidirectional Mendelian Randomisation Study

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

Gene set of nuclear-encoded mitochondrial regulators is enriched for common inherited variation in obesity

There are hints of an altered mitochondrial function in obesity. Nuclear-encoded genes are relevant for mitochondrial function (3 gene sets of known relevant pathways: (1) 16 nuclear regulators of mitochondrial genes, (2) 91 genes for oxidative phosphorylation and (3) 966 nuclear-encoded mitochondrial genes). Gene set enrichment analysis (GSEA) showed no association with type 2 diabetes mellitus in these gene sets. Here we performed a GSEA for the same gene sets for obesity. Genome wide association study (GWAS) data from a case-control approach on 453 extremely obese children and adolescents and 435 lean adult controls were used for GSEA. For independent confirmation, we analyzed 705 obesity GWAS trios (extremely obese child and both biological parents) and a population-based GWAS sample (KORA F4, n = 1,743). A meta-analysis was performed on all three samples. In each sample, the distribution of significance levels between the respective gene set and those of all genes was compared using the leading-edge-fraction-comparison test (cut-offs between the 50(th) and 95(th) percentile of the set of all gene-wise corrected p-values) as implemented in the MAGENTA software. In the case-control sample, significant enrichment of associations with obesity was observed above the 50(th) percentile for the set of the 16 nuclear regulators of mitochondrial genes (p(GSEA,50) = 0.0103). This finding was not confirmed in the trios (p(GSEA,50) = 0.5991), but in KORA (p(GSEA,50) = 0.0398). The meta-analysis again indicated a trend for enrichment (p(MAGENTA,50) = 0.1052, p(MAGENTA,75) = 0.0251). The GSEA revealed that weak association signals for obesity might be enriched in the gene set of 16 nuclear regulators of mitochondrial genes

Impact of FTO genotypes on BMI and weight in polycystic ovary syndrome : a systematic review and meta-analysis

Aims/hypothesis FTO gene single nucleotide polymorphisms (SNPs) have been shown to be associated with obesity-related traits and type 2 diabetes. Several small studies have suggested a greater than expected effect of the FTO rs9939609 SNP on weight in polycystic ovary syndrome (PCOS). We therefore aimed to examine the impact of FTO genotype on BMI and weight in PCOS. Methods A systematic search of medical databases (PubMed, EMBASE and Cochrane CENTRAL) was conducted up to the end of April 2011. Seven studies describing eight distinct PCOS cohorts were retrieved; seven were genotyped for SNP rs9939609 and one for SNP rs1421085. The per allele effect on BMI and body weight increase was calculated and subjected to meta-analysis. Results A total of 2,548 women with PCOS were included in the study; 762 were TT homozygotes, 1,253 had an AT/CT genotype, and 533 were AA/CC homozygotes. Each additional copy of the effect allele (A/C) increased the BMI by a mean of 0.19 z score units (95% CI 0.13, 0.24; p = 2.26 × 10−11) and body weight by a mean of 0.20 z score units (95% CI 0.14, 0.26; p = 1.02 × 10−10). This translated into an approximately 3.3 kg/m2 increase in BMI and an approximately 9.6 kg gain in body weight between TT and AA/CC homozygotes. The association between FTO genotypes and BMI was stronger in the cohorts with PCOS than in the general female populations from large genome-wide association studies. Deviation from an additive genetic model was observed in heavier populations. Conclusions/interpretation The effect of FTO SNPs on obesity-related traits in PCOS seems to be more than two times greater than the effect found in large population-based studies. This suggests an interaction between FTO and the metabolic context or polygenic background of PCOS

A Systematic Mapping Approach of 16q12.2/FTO and BMI in More Than 20,000 African Americans Narrows in on the Underlying Functional Variation: Results from the Population Architecture using Genomics and Epidemiology (PAGE) Study

Genetic variants in intron 1 of the fat mass- and obesity-associated (FTO) gene have been consistently associated with body mass index (BMI) in Europeans. However, follow-up studies in African Americans (AA) have shown no support for some of the most consistently BMI-associated FTO index single nucleotide polymorphisms (SNPs). This is most likely explained by different race-specific linkage disequilibrium (LD) patterns and lower correlation overall in AA, which provides the opportunity to fine-map this region and narrow in on the functional variant. To comprehensively explore the 16q12.2/FTO locus and to search for second independent signals in the broader region, we fine-mapped a 646-kb region, encompassing the large FTO gene and the flanking gene RPGRIP1L by investigating a total of 3,756 variants (1,529 genotyped and 2,227 imputed variants) in 20,488 AAs across five studies. We observed associations between BMI and variants in the known FTO intron 1 locus: the SNP with the most significant p-value, rs56137030 (8.3×10-6) had not been highlighted in previous studies. While rs56137030was correlated at r2>0.5 with 103 SNPs in Europeans (including the GWAS index SNPs), this number was reduced to 28 SNPs in AA. Among rs56137030 and the 28 correlated SNPs, six were located within candidate intronic regulatory elements, including rs1421085, for which we predicted allele-specific binding affinity for the transcription factor CUX1, which has recently been implicated in the regulation of FTO. We did not find strong evidence for a second independent signal in the broader region. In summary, this large fine-mapping study in AA has substantially reduced the number of common alleles that are likely to be functional candidates of the known FTO locus. Importantly our study demonstrated that comprehensive fine-mapping in AA provides a powerful approach to narrow in on the functional candidate(s) underlying the initial GWAS findings in European populations

MultiPhen: Joint Model of Multiple Phenotypes Can Increase Discovery in GWAS

The genome-wide association study (GWAS) approach has discovered hundreds of genetic variants associated with diseases and quantitative traits. However, despite clinical overlap and statistical correlation between many phenotypes, GWAS are generally performed one-phenotype-at-a-time. Here we compare the performance of modelling multiple phenotypes jointly with that of the standard univariate approach. We introduce a new method and software, MultiPhen, that models multiple phenotypes simultaneously in a fast and interpretable way. By performing ordinal regression, MultiPhen tests the linear combination of phenotypes most associated with the genotypes at each SNP, and thus potentially captures effects hidden to single phenotype GWAS. We demonstrate via simulation that this approach provides a dramatic increase in power in many scenarios. There is a boost in power for variants that affect multiple phenotypes and for those that affect only one phenotype. While other multivariate methods have similar power gains, we describe several benefits of MultiPhen over these. In particular, we demonstrate that other multivariate methods that assume the genotypes are normally distributed, such as canonical correlation analysis (CCA) and MANOVA, can have highly inflated type-1 error rates when testing case-control or non-normal continuous phenotypes, while MultiPhen produces no such inflation. To test the performance of MultiPhen on real data we applied it to lipid traits in the Northern Finland Birth Cohort 1966 (NFBC1966). In these data MultiPhen discovers 21% more independent SNPs with known associations than the standard univariate GWAS approach, while applying MultiPhen in addition to the standard approach provides 37% increased discovery. The most associated linear combinations of the lipids estimated by MultiPhen at the leading SNPs accurately reflect the Friedewald Formula, suggesting that MultiPhen could be used to refine the definition of existing phenotypes or uncover novel heritable phenotypes

Quantile-Specific Penetrance of Genes Affecting Lipoproteins, Adiposity and Height

Quantile-dependent penetrance is proposed to occur when the phenotypic expression of a SNP depends upon the population percentile of the phenotype. To illustrate the phenomenon, quantiles of height, body mass index (BMI), and plasma lipids and lipoproteins were compared to genetic risk scores (GRS) derived from single nucleotide polymorphisms (SNP)s having established genome-wide significance: 180 SNPs for height, 32 for BMI, 37 for low-density lipoprotein (LDL)-cholesterol, 47 for high-density lipoprotein (HDL)-cholesterol, 52 for total cholesterol, and 31 for triglycerides in 1930 subjects. Both phenotypes and GRSs were adjusted for sex, age, study, and smoking status. Quantile regression showed that the slope of the genotype-phenotype relationships increased with the percentile of BMI (P = 0.002), LDL-cholesterol (P = 3×10−8), HDL-cholesterol (P = 5×10−6), total cholesterol (P = 2.5×10−6), and triglyceride distribution (P = 7.5×10−6), but not height (P = 0.09). Compared to a GRS's phenotypic effect at the 10th population percentile, its effect at the 90th percentile was 4.2-fold greater for BMI, 4.9-fold greater for LDL-cholesterol, 1.9-fold greater for HDL-cholesterol, 3.1-fold greater for total cholesterol, and 3.3-fold greater for triglycerides. Moreover, the effect of the rs1558902 (FTO) risk allele was 6.7-fold greater at the 90th than the 10th percentile of the BMI distribution, and that of the rs3764261 (CETP) risk allele was 2.4-fold greater at the 90th than the 10th percentile of the HDL-cholesterol distribution. Conceptually, it maybe useful to distinguish environmental effects on the phenotype that in turn alters a gene's phenotypic expression (quantile-dependent penetrance) from environmental effects affecting the gene's phenotypic expression directly (gene-environment interaction)

A Genome-Wide Association Study on Obesity and Obesity-Related Traits

Large-scale genome-wide association studies (GWAS) have identified many loci associated with body mass index (BMI), but few studies focused on obesity as a binary trait. Here we report the results of a GWAS and candidate SNP genotyping study of obesity, including extremely obese cases and never overweight controls as well as families segregating extreme obesity and thinness. We first performed a GWAS on 520 cases (BMI>35 kg/m2) and 540 control subjects (BMI<25 kg/m2), on measures of obesity and obesity-related traits. We subsequently followed up obesity-associated signals by genotyping the top ∼500 SNPs from GWAS in the combined sample of cases, controls and family members totaling 2,256 individuals. For the binary trait of obesity, we found 16 genome-wide significant signals within the FTO gene (strongest signal at rs17817449, P = 2.5×10−12). We next examined obesity-related quantitative traits (such as total body weight, waist circumference and waist to hip ratio), and detected genome-wide significant signals between waist to hip ratio and NRXN3 (rs11624704, P = 2.67×10−9), previously associated with body weight and fat distribution. Our study demonstrated how a relatively small sample ascertained through extreme phenotypes can detect genuine associations in a GWAS

Expression of cytoplasmic and nuclear Survivin in primary and secondary human glioblastoma

Clinically, human glioblastoma (GBM) may develop de novo or from a low-grade glioma (secondary GBM), and molecular alterations in the two pathways may differ. This study examined the status of Survivin expression and apoptosis in 30 primary and 26 secondary GBMs. Our results show that cytoplasmic Survivin positivity was significantly (P<0.001) more frequent in primary GBMs (83%) than that in secondary GBMs (46%). In addition, an inverse correlation of cytoplasmc Survivin positivity with GBM apoptotic index, and a positive association between cytoplasmic Survivin and size of the tumours were observed. These results suggest that cytoplasmic Survivin, via its antiapoptotic function, may be involved in the tumorigenesis of many primary GBMs, but only in a small fraction of secondary GBMs. Furthermore, the overall progression times from low-grade precursor lesions to secondary GBMs were significantly shorter (P<0.05) in cytoplasmic Survivin-positive cases (mean, 15.6 months) than those in Survivin-negative cases (mean, 23.8 moths), and the positive expression level of Survivin in cytoplasm was upregulated in most secondary GBMs when compared to matched pre-existing low-graded lesions. These results suggest that the increased accumulation of Survivin in the cytoplasm of more malignant glioma cells may prove to be a selective advantage, thus accelerating progression to a more aggressive phenotype