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

Meta-analysis of genome-wide association studies with correlated individuals: application to the Hispanic Community Health Study/Study of Latinos (HCHS/SOL)

Investigators often meta-analyze multiple genome-wide association studies (GWASs) to increase the power to detect associations of single nucleotide polymorphisms (SNPs) with a trait. Meta-analysis is also performed within a single cohort that is stratified by, e.g., sex or ancestry group. Having correlated individuals among the strata may complicate meta-analyses, limit power, and inflate Type 1 error. For example, in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL), sources of correlation include genetic relatedness, shared household, and shared community. We propose a novel mixed-effect model for meta-analysis, “MetaCor , which accounts for correlation between stratum-specific effect estimates. Simulations show that MetaCor controls inflation better than alternatives such as ignoring the correlation between the strata or analyzing all strata together in a “pooled GWAS, especially with different minor allele frequencies (MAF) between strata. We illustrate the benefits of MetaCor on two GWASs in the HCHS/SOL. Analysis of dental caries (tooth decay) stratified by ancestry group detected a genome-wide significant SNP (rs7791001, p-value=3.66x10-8, compared to 4.67x10-7 in pooled), with different MAF between strata. Stratified analysis of BMI by ancestry group and sex reduced over-all inflation from λGC=1.050 (pooled) to λGC=1.028 (MetaCor). Furthermore, even after removing close relatives to obtain nearly uncorrelated strata, a naïve stratified analysis resulted in λGC=1.058 compare to λGC=1.027 for MetaCor

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

Diverse impacts of the rs58542926 E167K variant in TM6SF2 on viral and metabolic liver disease phenotypes

A genome-wide exome association study has identified the transmembrane 6 superfamily member 2 (TM6SF2) rs58542926 variant encoding an E167K substitution as a genetic determinant of hepatic steatosis in nonalcoholic fatty liver disease (NAFLD). The roles of this variant across a spectrum of liver diseases and pathologies and on serum lipids comparing viral hepatitis to NAFLD and viral load in chronic viral hepatitis, as well as its intrahepatic molecular signature, have not been well characterized. We undertook detailed analyses in 3260 subjects with viral and nonviral liver diseases and in healthy controls. Serum inflammatory markers and hepatic expression of TM6SF2 and genes regulating lipid metabolism were assessed in a subset with chronic hepatitis C (CHC). The rs58542926 T allele was more prevalent in 502 NAFLD patients than controls (P = 0.02) but not different in cohorts with CHC (n = 2023) and chronic hepatitis B (n = 507). The T allele was associated with alterations in serum lipids and hepatic steatosis in all diseases and with reduced hepatic TM6SF2 and microsomal triglyceride transfer protein expression. Interestingly, the substitution was associated with reduced CHC viral load but increased hepatitis B virus DNA. The rs58542926 T allele had no effect on inflammation, impacted >= F2 fibrosis in CHC and NAFLD assessed cross-sectionally (odds ratio = 1.39, 95% confidence interval 1.04-1.87, and odds ratio = 1.62, 95% confidence interval 1.03-2.52, respectively; P < 0.03 for both), but had no effect on fibrosis progression in 1174 patients with CHC and a known duration of infection. Conclusion: The TM6SF2 E167K substitution promotes steatosis and lipid abnormalities in part by altering TM6SF2 and microsomal triglyceride transfer protein expression and differentially impacts CHC and chronic hepatitis B viral load, while effects on fibrosis are marginal

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition

A genome-wide association study of sleep habits and insomnia

Several aspects of sleep behavior such as timing, duration and quality have been demonstrated to be heritable. To identify common variants that influence sleep traits in the population, we conducted a genome-wide association study of six sleep phenotypes assessed by questionnaire in a sample of 2,323 individuals from the Australian Twin Registry. Genotyping was performed on the Illumina 317, 370, and 610K arrays and the SNPs in common between platforms were used to impute non-genotyped SNPs. We tested for association with more than 2,000,000 common polymorphisms across the genome. While no SNPs reached the genome-wide significance threshold, we identified a number of associations in plausible candidate genes. Most notably, a group of SNPs in the third intron of the CACNA1C gene ranked as most significant in the analysis of sleep latency (P=1.3×10-6). We attempted to replicate this association in an independent sample from the Chronogen Consortium (n=2,034), but found no evidence of association (P=0.73). We have identified several other suggestive associations that await replication in an independent sample. We did not replicate the results from previous genome-wide analyses of self-reported sleep phenotypes after correction for multiple testing

Variation at the Melanocortin 4 Receptor gene and response to weight-loss interventions in the Diabetes Prevention Program

Objective: To assess associations and genotype × treatment interactions for melanocortin 4 receptor (MC4R) locus variants and obesity-related traits. Design and Methods Diabetes Prevention Program (DPP) participants (N=3,819, of whom 3,356 were genotyped for baseline and 3,234 for longitudinal analyses) were randomized into intensive lifestyle modification (diet, exercise, weight loss), metformin or placebo control. Adiposity was assessed in a subgroup (n=909) using computed tomography. All analyses were adjusted for age, sex, ethnicity and treatment. Results: The rs1943218 minor allele was nominally associated with short-term (6 month; P=0.032) and long-term (2 year; P=0.038) weight change. Eight SNPs modified response to treatment on short-term (rs17066856, rs9966412, rs17066859, rs8091237, rs17066866, rs7240064) or long-term (rs12970134, rs17066866) reduction in body weight, or diabetes incidence (rs17066829) (all Pinteraction <0.05). Conclusion: This is the first study to comprehensively assess the role of MC4R variants and weight regulation in a weight loss intervention trial. One MC4R variant was directly associated with obesity-related traits or diabetes; numerous other variants appear to influence body weight and diabetes risk by modifying the protective effects of the DPP interventions

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

Multi-scale modeling of gene-behavior associations in an artificial neural network model of cognitive development

In the multi-disciplinary field of developmental cognitive neuroscience, statistical associations between levels of description play an increasingly important role. One example of such associations is the observation of correlations between relatively common gene variants and individual differences in behavior. It is perhaps surprising that such associations can be detected despite the remoteness of these levels of description, and the fact that behavior is the outcome of an extended developmental process involving interaction with a variable environment. Given that they have been detected, how do such associations inform cognitive-level theories? To investigate this question, we employed a multi-scale computational model of development, using a sample domain drawn from the field of language acquisition. The model comprised an artificial neural network model of past-tense acquisition trained using the backpropagation learning algorithm, extended to incorporate population modeling and genetic algorithms. It included five levels of description, four internal: genetic, network, neurocomputation, behavior; and one external: environment. Since the mechanistic assumptions of the model were known and its operation was relatively transparent, we could evaluate whether cross-level associations gave an accurate picture of causal processes. We established that associations could be detected between artificial genes and behavioral variation, even under polygenic assumptions of a many-to-one relationship between genes and neurocomputational parameters, and when an experience-dependent developmental process interceded between the action of genes and the emergence of behavior. We evaluated these associations with respect to their specificity (to different behaviors, to function versus structure), to their developmental stability, and to their replicability, as well as considering issues of missing heritability and gene-environment interactions. We argue that gene-behavior associations can inform cognitive theory with respect to effect size, specificity, and timing. The model demonstrates a means by which researchers can undertake modeling multi-scale modeling with respect to cognition, and develop highly specific and complex hypotheses across multiple levels of description

Genetic evidence for a normal-weight "metabolically obese" phenotype linking insulin resistance, hypertension, coronary artery disease, and type 2 diabetes

PublishedJournal ArticleResearch Support, Non-U.S. Gov'tThe mechanisms that predispose to hypertension, coronary artery disease (CAD), and type 2 diabetes (T2D) in individuals of normal weight are poorly understood. In contrast, in monogenic primary lipodystrophy-a reduction in subcutaneous adipose tissue-it is clear that it is adipose dysfunction that causes severe insulin resistance (IR), hypertension, CAD, and T2D. We aimed to test the hypothesis that common alleles associated with IR also influence the wider clinical and biochemical profile of monogenic IR. We selected 19 common genetic variants associated with fasting insulin-based measures of IR. We used hierarchical clustering and results from genome-wide association studies of eight nondisease outcomes of monogenic IR to group these variants. We analyzed genetic risk scores against disease outcomes, including 12,171 T2D cases, 40,365 CAD cases, and 69,828 individuals with blood pressure measurements. Hierarchical clustering identified 11 variants associated with a metabolic profile consistent with a common, subtle form of lipodystrophy. A genetic risk score consisting of these 11 IR risk alleles was associated with higher triglycerides (β = 0.018; P = 4 × 10(-29)), lower HDL cholesterol (β = -0.020; P = 7 × 10(-37)), greater hepatic steatosis (β = 0.021; P = 3 × 10(-4)), higher alanine transaminase (β = 0.002; P = 3 × 10(-5)), lower sex-hormone-binding globulin (β = -0.010; P = 9 × 10(-13)), and lower adiponectin (β = -0.015; P = 2 × 10(-26)). The same risk alleles were associated with lower BMI (per-allele β = -0.008; P = 7 × 10(-8)) and increased visceral-to-subcutaneous adipose tissue ratio (β = -0.015; P = 6 × 10(-7)). Individuals carrying ≥17 fasting insulin-raising alleles (5.5% population) were slimmer (0.30 kg/m(2)) but at increased risk of T2D (odds ratio [OR] 1.46; per-allele P = 5 × 10(-13)), CAD (OR 1.12; per-allele P = 1 × 10(-5)), and increased blood pressure (systolic and diastolic blood pressure of 1.21 mmHg [per-allele P = 2 × 10(-5)] and 0.67 mmHg [per-allele P = 2 × 10(-4)], respectively) compared with individuals carrying ≤9 risk alleles (5.5% population). Our results provide genetic evidence for a link between the three diseases of the "metabolic syndrome" and point to reduced subcutaneous adiposity as a central mechanism