Exome sequencing and characterization of 49,960 individuals in the UK Biobank.

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes.

Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10-8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10-17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10-15), HERC2-OCA2 (P = 4.2 × 10-12), SLC45A2 (P = 1.7 × 10-10), IL13 (P = 2.8 × 10-9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10-9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10-8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10-7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.

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.

Factors influencing longitudinal changes of circulating liver enzyme concentrations in subjects randomized to placebo in four clinical trials.

Liver enzyme concentrations are measured as safety end points in clinical trials to detect drug-related hepatotoxicity, but little is known about the epidemiology of these biomarkers in subjects without hepatic dysfunction who are enrolled in drug trials. We studied alanine and aspartate aminotransferase (ALT and AST) in subjects randomized to placebo who completed assessments over 36 mo in a cardiovascular outcome trial [the Stabilisation of Atherosclerotic Plaque by Initiation of Darapladib Therapy ("STABILITY") trial; n = 4,264; mean age: 64.2 yr] or over 12 mo in three trials that enrolled only subjects with type 2 diabetes (T2D) [the DIA trials; n = 308; mean age: 62.4 yr] to investigate time-dependent relationships and the factors that might affect ALT and AST, including body mass index (BMI), T2D, and renal function. Multivariate linear mixed models examined time-dependent relationships between liver enzyme concentrations as response variables and BMI, baseline T2D status, hemoglobin A1c levels, and renal function, as explanatory variables. At baseline, ALT was higher in individuals who were men, <65 yr old, and obese and who had glomerular filtration rate (GFR) >60 ml·min-1·1.73 m-2. ALT was not significantly associated with T2D at baseline, although it was positively associated with HbA1c. GFR had a greater impact on ALT than T2D. ALT concentrations decreased over time in subjects who lost weight but remained stable in individuals with increasing BMI. Weight change did not alter AST concentrations. We provide new insights on the influence of time, GFR, and HbA1c on ALT and AST concentrations and confirm the effect of sex, age, T2D, BMI, and BMI change in subjects receiving placebo in clinical trials. NEW & NOTEWORTHY Clinical trials provide high-quality data on liver enzyme concentrations from subjects randomized to placebo that can be used to investigate the epidemiology of these biomarkers. The adjusted models show the influence of sex, age, time, renal function, type 2 diabetes, HbA1c, and body mass index on alanine aminotransferase and aspartate aminotransferase concentrations and their relative importance. These factors need to be considered when assessing potential signals of hepatotoxicity in trials of new drugs and in clinical trials investigating subjects with nonalcoholic fatty liver disease.

Exome-chip meta-analysis identifies association between variation in ANKRD26 and platelet aggregation.

Previous genome-wide association studies (GWAS) have identified several variants associated with platelet function phenotypes; however, the proportion of variance explained by the identified variants is mostly small. Rare coding variants, particularly those with high potential for impact on protein structure/function, may have substantial impact on phenotype but are difficult to detect by GWAS. The main purpose of this study was to identify low frequency or rare variants associated with platelet function using genotype data from the Illumina HumanExome Bead Chip. Three family-based cohorts of European ancestry, including ~4,000 total subjects, comprised the discovery cohort and two independent cohorts, one of European and one of African American ancestry, were used for replication. Optical aggregometry in platelet-rich plasma was performed in all the discovery cohorts in response to adenosine diphosphate (ADP), epinephrine, and collagen. Meta-analyses were performed using both gene-based and single nucleotide variant association methods. The gene-based meta-analysis identified a significant association (P = 7.13 × 10-7) between rare genetic variants in ANKRD26 and ADP-induced platelet aggregation. One of the ANKRD26 SNVs - rs191015656, encoding a threonine to isoleucine substitution predicted to alter protein structure/function, was replicated in Europeans. Aggregation increases of ~20-50% were observed in heterozygotes in all cohorts. Novel genetic signals in ABCG1 and HCP5 were also associated with platelet aggregation to ADP in meta-analyses, although only results for HCP5 could be replicated. The SNV in HCP5 intersects epigenetic signatures in CD41+ megakaryocytes suggesting a new functional role in platelet biology for HCP5. This is the first study to use gene-based association methods from SNV array genotypes to identify rare variants related to platelet function. The molecular mechanisms and pathophysiological relevance for the identified genetic associations requires further study.

TM6SF2 rs58542926 impacts lipid processing in liver and small intestine.

The transmembrane 6 superfamily member 2 (TM6SF2) loss-of-function variant rs58542926 is a genetic risk factor for nonalcoholic fatty liver disease and progression to fibrosis but is paradoxically associated with lower levels of hepatically derived triglyceride-rich lipoproteins. TM6SF2 is expressed predominantly in liver and small intestine, sites for triglyceride-rich lipoprotein biogenesis and export. In light of this, we hypothesized that TM6SF2 may exhibit analogous effects on both liver and intestine lipid homeostasis. To test this, we genotyped rs58542926 in 983 bariatric surgery patients from the Geisinger Medical Center for Nutrition and Weight Management, Geisinger Health System, in Pennsylvania and from 3,556 study participants enrolled in the Amish Complex Disease Research Program. Although these two cohorts have different metabolic profiles, carriers in both cohorts had improved fasting lipid profiles. Importantly, following a high-fat challenge, carriers in the Amish Complex Disease Research Program cohort exhibited significantly lower postprandial serum triglycerides, suggestive of a role for TM6SF2 in the small intestine. To gain further insight into this putative role, effects of TM6SF2 deficiency were studied in a zebrafish model and in cultured human Caco-2 enterocytes. In both systems TM6SF2 deficiency resulted in defects in small intestine metabolism in response to dietary lipids, including significantly increased lipid accumulation, decreased lipid clearance, and increased endoplasmic reticulum stress. CONCLUSIONS: These data strongly support a role of TM6SF2 in the regulation of postprandial lipemia, potentially through a similar function for TM6SF2 in the lipidation and/or export of both hepatically and intestinally derived triglyceride-rich lipoproteins. (Hepatology 2017;65:1526-1542).

FTO genotype is associated with phenotypic variability of body mass index.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Genome-wide analysis of clopidogrel active metabolite levels identifies novel variants that influence antiplatelet response.

Clopidogrel is one of the most commonly used therapeutics for the secondary prevention of cardiovascular events in patients with acute coronary syndromes. However, considerable interindividual variation in clopidogrel response has been documented, resulting in suboptimal therapy and an increased risk of recurrent events for some patients. In this investigation, we carried out the first genome-wide association study of circulating clopidogrel active metabolite levels in 513 healthy participants to directly measure clopidogrel pharmacokinetics. We observed that the CYP2C19 locus was the strongest genetic determinant of active metabolite formation (P=9.5×10). In addition, we identified novel genome-wide significant variants on chromosomes 3p25 (rs187941554, P=3.3×10) and 17q11 (rs80343429, P=1.3×10), as well as six additional loci that showed suggestive evidence of association (P≤1.0×10). Four of these loci showed nominal associations with on-clopidogrel ADP-stimulated platelet aggregation (P≤0.05). Evaluation of clopidogrel active metabolite concentration may help identify novel genetic determinants of clopidogrel response, which has implications for the development of novel therapeutics and improved antiplatelet treatment for at-risk patients in the future.

Null mutation in hormone-sensitive lipase gene and risk of type 2 diabetes.

BACKGROUND: Lipolysis regulates energy homeostasis through the hydrolysis of intracellular triglycerides and the release of fatty acids for use as energy substrates or lipid mediators in cellular processes. Genes encoding proteins that regulate energy homeostasis through lipolysis are thus likely to play an important role in determining susceptibility to metabolic disorders. METHODS: We sequenced 12 lipolytic-pathway genes in Old Order Amish participants whose fasting serum triglyceride levels were at the extremes of the distribution and identified a novel 19-bp frameshift deletion in exon 9 of LIPE, encoding hormone-sensitive lipase (HSL), a key enzyme for lipolysis. We genotyped the deletion in DNA from 2738 Amish participants and performed association analyses to determine the effects of the deletion on metabolic traits. We also obtained biopsy specimens of abdominal subcutaneous adipose tissue from 2 study participants who were homozygous for the deletion (DD genotype), 10 who were heterozygous (ID genotype), and 7 who were noncarriers (II genotype) for assessment of adipose histologic characteristics, lipolysis, enzyme activity, cytokine release, and messenger RNA (mRNA) and protein levels. RESULTS: Carriers of the mutation had dyslipidemia, hepatic steatosis, systemic insulin resistance, and diabetes. In adipose tissue from study participants with the DD genotype, the mutation resulted in the absence of HSL protein, small adipocytes, impaired lipolysis, insulin resistance, and inflammation. Transcription factors responsive to peroxisome-proliferator-activated receptor γ (PPAR-γ) and downstream target genes were down-regulated in adipose tissue from participants with the DD genotype, altering the regulation of pathways influencing adipogenesis, insulin sensitivity, and lipid metabolism. CONCLUSIONS: These findings indicate the physiological significance of HSL in adipocyte function and the regulation of systemic lipid and glucose homeostasis and underscore the severe metabolic consequences of impaired lipolysis. (Funded by the National Institutes of Health and others).

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density.

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10(-8)) level: 14q24.2 (rs227425, P-value 3.98 × 10(-13), SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10(-9), CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n = 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis.

AP2 suppresses osteoblast differentiation and mineralization through down-regulation of Frizzled-1.

Transcription factor activating protein 2 (AP2) plays an important role in cellular differentiation. Although profound craniofacial and long bone developmental abnormalities have been observed in AP2-knockout mice, the molecular effects of AP2 on osteoblasts are poorly defined. We demonstrated that AP2 regulates the expression of human Frizzled 1 (FZD1), a co-receptor for the Wnt signalling pathway, in human osteoblast cell lines and primary bone marrow stromal cells (BMSCs). We also identified a putative AP2-binding site in the FZD1 proximal promoter in silico and characterized this binding element further in Saos2 in vitro by ChIP, electrophoretic mobility shift and promoter reporter assays. The transcriptional repression of the FZD1 promoter by AP2 was confirmed in normal human fetal osteoblasts (hFOB). Furthermore, overexpression of AP2 resulted in a significant reduction in both differentiation and mineralization of Saos2 cells. Knockdown of FZD1 expression before AP2 up-regulation diminished the AP2-dependent inhibition of Saos2 cell differentiation and mineralization. Similarly, overexpressing FZD1 before AP2 treatment in both Saos2 and BMSCs diminished the inhibitory effect of AP2 on osteoblast differentiation and mineralization. Taken together, these results demonstrate that AP2 is a negative regulator of osteoblast differentiation and mineralization, and its inhibitory effect may be mediated in part through down-regulation of FZD1 expression.

Pain, Genes, and Function in the Post-Hip Fracture Period.

Post-hip fracture generalized pain can lead to a progressive decline in function and greater disability. The purpose of this study was to explore the factors that influence pain among older adults post-hip fracture, including genetic variability, and evaluate whether pain directly or indirectly influenced upper and lower extremity function. This was a secondary data analysis using data from the first 200 participants in a Baltimore Hip Study (BHS), BHS-7. Assessments were done at 2 months post-hip fracture and included age, sex, marital status, education, cognitive status, comorbidities, body mass index (BMI), upper and lower extremity function, single nucleotide polymorphisms (SNPs) from 10 candidate genes, and total areas of pain and pain intensity. Model testing was done using the AMOS statistical program. The full sample included 172 participants with an average age of 81. Fifty percent were female and the majority was Caucasian (93%). Model testing was done on 144 individuals who completed 2 month surveys. Across all models, age, cognition, and BMI were significantly associated with total areas of pain. Thirty SNPs from five genes (BDNF, FKBP5, NTRK2, NTRK3, and OXTR) were associated with areas of pain and/or pain intensity. Together, age, cognition, BMI, and the SNP from one of the five genes explained 25% of total areas of pain and 15% of pain intensity. Only age and cognition were significantly associated with lower extremity function, and only cognition was significantly associated with upper extremity function. The full model was partially supported in this study. Our genetic findings related to pain expand prior reports related to BDNF and NTRK2.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk.

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)

Characterization of European ancestry nonalcoholic fatty liver disease-associated variants in individuals of African and Hispanic descent.

UNLABELLED: Nonalcoholic fatty liver disease (NAFLD) is an obesity-related condition affecting over 50% of individuals in some populations and is expected to become the number one cause of liver disease worldwide by 2020. Common, robustly associated genetic variants in/near five genes were identified for hepatic steatosis, a quantifiable component of NAFLD, in European ancestry individuals. Here we tested whether these variants were associated with hepatic steatosis in African- and/or Hispanic-Americans and fine-mapped the observed association signals. We measured hepatic steatosis using computed tomography in five African American (n = 3,124) and one Hispanic American (n = 849) cohorts. All analyses controlled for variation in age, age(2) , gender, alcoholic drinks, and population substructure. Heritability of hepatic steatosis was estimated in three cohorts. Variants in/near PNPLA3, NCAN, LYPLAL1, GCKR, and PPP1R3B were tested for association with hepatic steatosis using a regression framework in each cohort and meta-analyzed. Fine-mapping across African American cohorts was conducted using meta-analysis. African- and Hispanic-American cohorts were 33.9/37.5% male, with average age of 58.6/42.6 years and body mass index of 31.8/28.9 kg/m(2) , respectively. Hepatic steatosis was 0.20-0.34 heritable in African- and Hispanic-American families (P < 0.02 in each cohort). Variants in or near PNPLA3, NCAN, GCKR, PPP1R3B in African Americans and PNPLA3 and PPP1R3B in Hispanic Americans were significantly associated with hepatic steatosis; however, allele frequency and effect size varied across ancestries. Fine-mapping in African Americans highlighted missense variants at PNPLA3 and GCKR and redefined the association region at LYPLAL1. CONCLUSION: Multiple genetic variants are associated with hepatic steatosis across ancestries. This explains a substantial proportion of the genetic predisposition in African- and Hispanic-Americans. Missense variants in PNPLA3 and GCKR are likely functional across multiple ancestries.

Meta-analysis of genome-wide studies identifies MEF2C SNPs associated with bone mineral density at forearm.

BACKGROUND: Forearm fractures affect 1.7 million individuals worldwide each year and most occur earlier in life than hip fractures. While the heritability of forearm bone mineral density (BMD) and fracture is high, their genetic determinants are largely unknown. AIM: To identify genetic variants associated with forearm BMD and forearm fractures. METHODS: BMD at distal radius, measured by dual-energy x-ray absorptiometry, was tested for association with common genetic variants. We conducted a meta-analysis of genome-wide association studies for BMD in 5866 subjects of European descent and then selected the variants for replication in 715 Mexican American samples. Gene-based association was carried out to supplement the single-nucleotide polymorphism (SNP) association test. We then tested the BMD-associated SNPs for association with forearm fracture in 2023 cases and 3740 controls. RESULTS: We found that five SNPs in the introns of MEF2C were associated with forearm BMD at a genome-wide significance level (p<5×10(-8)) in meta-analysis (lead SNP, rs11951031[T] -0.20 SDs per allele, p=9.01×10(-9)). The gene-based association test suggested an association between MEF2C and forearm BMD (p=0.003). The association between MEF2C variants and risk of fracture did not achieve statistical significance (SNP rs12521522[A]: OR=1.14 (95% CI 0.92 to 1.35), p=0.14). Meta-analysis also revealed two genome-wide suggestive loci at CTNNA2 and 6q23.2. CONCLUSIONS: These findings demonstrate that variants at MEF2C were associated with forearm BMD, implicating this gene in the determination of BMD at forearm.

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.

Genetic variation at NCAN locus is associated with inflammation and fibrosis in non-alcoholic fatty liver disease in morbid obesity.

OBJECTIVE: Obesity-associated non-alcoholic fatty liver disease (NAFLD) may cause liver dysfunction and failure. In a previously reported genome-wide association meta-analysis, single nucleotide polymorphisms (SNPs) near PNPLA3, NCAN, GCKR, LYPLAL1 and PPP1R3B were associated with NAFLD and with distinctive serum lipid profiles. The present study examined the relevance of these variants to NAFLD in extreme obesity. METHODS: In 1,092 bariatric surgery patients, the candidate SNPs were genotyped and association analyses with liver histology and serum lipids were performed. RESULTS: We replicated the association of hepatosteatosis with PNPLA3 rs738409[G] and with NCAN rs2228603[T]. We also replicated the association of rs2228603[T] with hepatic inflammation and fibrosis. rs2228603[T] was associated with lower serum low-density lipoprotein, total cholesterol and triglycerides. After stratification by the presence or absence of NAFLD, these associations were present predominantly in the subgroup with NAFLD. CONCLUSION: NCAN rs2228603[T] is a risk factor for liver inflammation and fibrosis, suggesting that this locus is responsible for progression from steatosis to steatohepatitis. In this bariatric cohort, rs2228603[T] was associated with low serum lipids only in patients with NAFLD. This supports a NAFLD model in which the liver may sequester triglycerides as a result of either increased triglyceride uptake and/or decreased lipolysis.

The gout epidemic in French Polynesia: a modelling study of data from the Ma'i u'u epidemiological survey.

BACKGROUND: Gout is the most common cause of inflammatory arthritis worldwide, particularly in Pacific regions. We aimed to establish the prevalence of gout and hyperuricaemia in French Polynesia, their associations with dietary habits, their comorbidities, the prevalence of the HLA-B*58:01 allele, and current management of the disease. METHODS: The Ma'i u'u survey was epidemiological, prospective, cross-sectional, and gout-focused and included a random sample of adults from the general adult population of French Polynesia. It was conducted and data were collected between April 13 and Aug 16, 2021. Participants were randomly selected to represent the general adult population of French Polynesia on the basis of housing data collected during the 2017 territorial census. Each selected household was visited by a research nurse from the Ma'i u'u survey who collected data via guided, 1-h interviews with participants. In each household, the participant was the individual older than 18 years with the closest upcoming birthday. To estimate the frequency of HLA-B*58:01, we estimated HLA-B haplotypes on individuals who had whole-genome sequencing to approximately 5× average coverage (mid-pass sequencing). A subset of individuals who self-reported Polynesian ancestry and not European, Chinese, or other ancestry were used to estimate Polynesian-ancestry specific allele frequencies. Bivariate associations were reported for weighted participants; effect sizes were estimated through the odds ratio (OR) of the association calculated on the basis of a logistic model fitted with weighted observations. FINDINGS: Among the random sample of 2000 households, 896 participants were included, 140 individuals declined, and 964 households could not be contacted. 22 participants could not be weighted due to missing data, so the final weighted analysis included 874 participants (449 [51·4%] were female and 425 [48·6%] were male) representing the 196 630 adults living in French Polynesia. The estimated prevalence of gout was 14·5% (95% CI 9·9-19·2), representing 28 561 French Polynesian adults, that is 25·5% (18·2-32·8) of male individuals and 3·5% (1·0-6·0) of female individuals. The prevalence of hyperuricaemia was estimated at 71·6% (66·7-76·6), representing 128 687 French Polynesian adults. In multivariable analysis, age (OR 1·5, 95% CI 1·2-1·8 per year), male sex (10·3, 1·8-60·7), serum urate (1·6, 1·3-2·0 per 1 mg/dL), uraturia (0·8, 0·8-0·8 per 100 mg/L), type 2 diabetes (2·1, 1·4-3·1), BMI more than 30 kg/m2 (1·1, 1·0-1·2 per unit), and percentage of visceral fat (1·7, 1·1-2·7 per 1% increase) were associated with gout. There were seven heterozygous HLA-B*58:01 carriers in the full cohort of 833 individuals (seven [0·4%] of 1666 total alleles) and two heterozygous carriers in a subset of 696 individuals of Polynesian ancestry (two [0·1%]). INTERPRETATION: French Polynesia has an estimated high prevalence of gout and hyperuricaemia, with gout affecting almost 15% of adults. Territorial measures that focus on increasing access to effective urate-lowering therapies are warranted to control this major public health problem. FUNDING: Variant Bio, the French Polynesian Health Administration, Lille Catholic University Hospitals, French Society of Rheumatology, and Novartis.

Causal relationship between obesity and vitamin D status: bi-directional Mendelian randomization analysis of multiple cohorts.

BACKGROUND: Obesity is associated with vitamin D deficiency, and both are areas of active public health concern. We explored the causality and direction of the relationship between body mass index (BMI) and 25-hydroxyvitamin D [25(OH)D] using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis. METHODS AND FINDINGS: We used information from 21 adult cohorts (up to 42,024 participants) with 12 BMI-related SNPs (combined in an allelic score) to produce an instrument for BMI and four SNPs associated with 25(OH)D (combined in two allelic scores, separately for genes encoding its synthesis or metabolism) as an instrument for vitamin D. Regression estimates for the IVs (allele scores) were generated within-study and pooled by meta-analysis to generate summary effects. Associations between vitamin D scores and BMI were confirmed in the Genetic Investigation of Anthropometric Traits (GIANT) consortium (n = 123,864). Each 1 kg/m(2) higher BMI was associated with 1.15% lower 25(OH)D (p = 6.52×10⁻²7). The BMI allele score was associated both with BMI (p = 6.30×10⁻6²) and 25(OH)D (-0.06% [95% CI -0.10 to -0.02], p = 0.004) in the cohorts that underwent meta-analysis. The two vitamin D allele scores were strongly associated with 25(OH)D (p≤8.07×10⁻57 for both scores) but not with BMI (synthesis score, p = 0.88; metabolism score, p = 0.08) in the meta-analysis. A 10% higher genetically instrumented BMI was associated with 4.2% lower 25(OH)D concentrations (IV ratio: -4.2 [95% CI -7.1 to -1.3], p = 0.005). No association was seen for genetically instrumented 25(OH)D with BMI, a finding that was confirmed using data from the GIANT consortium (p≥0.57 for both vitamin D scores). CONCLUSIONS: On the basis of a bi-directional genetic approach that limits confounding, our study suggests that a higher BMI leads to lower 25(OH)D, while any effects of lower 25(OH)D increasing BMI are likely to be small. Population level interventions to reduce BMI are expected to decrease the prevalence of vitamin D deficiency.

Association of adiposity genetic variants with menarche timing in 92,105 women of European descent.

Obesity is of global health concern. There are well-described inverse relationships between female pubertal timing and obesity. Recent genome-wide association studies of age at menarche identified several obesity-related variants. Using data from the ReproGen Consortium, we employed meta-analytical techniques to estimate the associations of 95 a priori and recently identified obesity-related (body mass index (weight (kg)/height (m)(2)), waist circumference, and waist:hip ratio) single-nucleotide polymorphisms (SNPs) with age at menarche in 92,116 women of European descent from 38 studies (1970-2010), in order to estimate associations between genetic variants associated with central or overall adiposity and pubertal timing in girls. Investigators in each study performed a separate analysis of associations between the selected SNPs and age at menarche (ages 9-17 years) using linear regression models and adjusting for birth year, site (as appropriate), and population stratification. Heterogeneity of effect-measure estimates was investigated using meta-regression. Six novel associations of body mass index loci with age at menarche were identified, and 11 adiposity loci previously reported to be associated with age at menarche were confirmed, but none of the central adiposity variants individually showed significant associations. These findings suggest complex genetic relationships between menarche and overall obesity, and to a lesser extent central obesity, in normal processes of growth and development.