Increased Plasma Branched Short-Chain Fatty Acids and Improved Glucose Homeostasis: The Microbiome and Insulin Longitudinal Evaluation Study (MILES).

Short-chain fatty acids (SCFAs) have been extensively studied for potential beneficial roles in glucose homeostasis and risk of diabetes; however, most of this research has focused on butyrate, acetate, and propionate. The effect on metabolism of branched SCFAs (BSCFAs; isobutyrate, isovalerate, and methylbutyrate) is largely unknown. In a cohort of 219 non-Hispanic White participants and 126 African American participants, we examined the association of BSCFA with dysglycemia (prediabetes and diabetes) and oral glucose tolerance test-based measures of glucose and insulin homeostasis, as well as with demographic, anthropometric, lifestyle, and lipid traits, and other SCFAs. We observed a bimodal distribution of BSCFAs, with 25 individuals having high levels (H-BSCFA group) and 320 individuals having lower levels (L-BSCFA group). The prevalence of dysglycemia was lower in the H-BSCFA group compared with the L-BSCFA group (16% vs. 49%; P = 0.0014). This association remained significant after adjustment for age, sex, race, BMI, and levels of other SCFAs. Consistent with the lower rate of dysglycemia, fasting and postprandial glucose levels were lower and the disposition index was higher in the H-BSCFA group. Additional findings in H-BSCFA versus L-BSCFA included lower fasting and postprandial C-peptide levels and lower insulin clearance without differences in insulin levels, insulin sensitivity, insulin secretion, or other variables examined, including diet and physical activity. As one of the first human studies associating higher BSCFA levels with lower odds of dysglycemia and improved glucose homeostasis, this study sets the stage for further investigation of BSCFA as a novel target for prevention or treatment of diabetes.

Genome-wide association meta-analysis identifies 17 loci associated with nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is common and partially heritable and has no effective treatments. We carried out a genome-wide association study (GWAS) meta-analysis of imaging (n = 66,814) and diagnostic code (3,584 cases versus 621,081 controls) measured NAFLD across diverse ancestries. We identified NAFLD-associated variants at torsin family 1 member B (TOR1B), fat mass and obesity associated (FTO), cordon-bleu WH2 repeat protein like 1 (COBLL1)/growth factor receptor-bound protein 14 (GRB14), insulin receptor (INSR), sterol regulatory element-binding transcription factor 1 (SREBF1) and patatin-like phospholipase domain-containing protein 2 (PNPLA2), as well as validated NAFLD-associated variants at patatin-like phospholipase domain-containing protein 3 (PNPLA3), transmembrane 6 superfamily 2 (TM6SF2), apolipoprotein E (APOE), glucokinase regulator (GCKR), tribbles homolog 1 (TRIB1), glycerol-3-phosphate acyltransferase (GPAM), mitochondrial amidoxime-reducing component 1 (MARC1), microsomal triglyceride transfer protein large subunit (MTTP), alcohol dehydrogenase 1B (ADH1B), transmembrane channel like 4 (TMC4)/membrane-bound O-acyltransferase domain containing 7 (MBOAT7) and receptor-type tyrosine-protein phosphatase δ (PTPRD). Implicated genes highlight mitochondrial, cholesterol and de novo lipogenesis as causally contributing to NAFLD predisposition. Phenome-wide association study (PheWAS) analyses suggest at least seven subtypes of NAFLD. Individuals in the top 10% and 1% of genetic risk have a 2.5-fold to 6-fold increased risk of NAFLD, cirrhosis and hepatocellular carcinoma. These genetic variants identify subtypes of NAFLD, improve estimates of disease risk and can guide the development of targeted therapeutics.

Participant-derived cell line transcriptomic analyses and mouse studies reveal a role for ZNF335 in plasma cholesterol statin response.

Background: Statins lower circulating low-density lipoprotein cholesterol (LDLC) levels and reduce cardiovascular disease risk. Though highly efficacious in general, there is considerable inter-individual variation in statin efficacy that remains largely unexplained. Methods: To identify novel genes that may modulate statin-induced LDLC lowering, we used RNA-sequencing data from 426 control- and 2 µM simvastatin-treated lymphoblastoid cell lines (LCLs) derived from European and African American ancestry participants of the Cholesterol and Pharmacogenetics (CAP) 40 mg/day 6-week simvastatin clinical trial (ClinicalTrials.gov Identifier: NCT00451828). We correlated statin-induced changes in LCL gene expression with plasma LDLC statin response in the corresponding CAP participants. For the most correlated gene identified (ZNF335), we followed up in vivo by comparing plasma cholesterol levels, lipoprotein profiles, and lipid statin response between wild-type mice and carriers of a hypomorphic (partial loss of function) missense mutation in Zfp335 (the mouse homolog of ZNF335). Results: The statin-induced expression changes of 147 human LCL genes were significantly correlated to the plasma LDLC statin responses of the corresponding CAP participants in vivo (FDR=5%). The two genes with the strongest correlations were zinc finger protein 335 (ZNF335 aka NIF-1, rho=0.237, FDR-adj p=0.0085) and CCR4-NOT transcription complex subunit 3 (CNOT3, rho=0.233, FDR-adj p=0.0085). Chow-fed mice carrying a hypomorphic missense (R1092W; aka bloto) mutation in Zfp335 had significantly lower non-HDL cholesterol levels than wild type C57BL/6J mice in a sex combined model (p=0.04). Furthermore, male (but not female) mice carrying the Zfp335R1092W allele had significantly lower total and HDL cholesterol levels than wild-type mice. In a separate experiment, wild-type mice fed a control diet for 4 weeks and a matched simvastatin diet for an additional 4 weeks had significant statin-induced reductions in non-HDLC (-43±18% and -23±19% for males and females, respectively). Wild-type male (but not female) mice experienced significant reductions in plasma LDL particle concentrations, while male mice carrying Zfp335R1092W allele(s) exhibited a significantly blunted LDL statin response. Conclusions: Our in vitro and in vivo studies identified ZNF335 as a novel modulator of plasma cholesterol levels and statin response, suggesting that variation in ZNF335 activity could contribute to inter-individual differences in statin clinical efficacy.

Multi-ancestry transcriptome-wide association analyses yield insights into tobacco use biology and drug repurposing.

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

Butyrate-Producing Bacteria and Insulin Homeostasis: The Microbiome and Insulin Longitudinal Evaluation Study (MILES).

Gut microbiome studies have documented depletion of butyrate-producing taxa in type 2 diabetes. We analyzed associations between butyrate-producing taxa and detailed measures of insulin homeostasis, whose dysfunction underlies diabetes in 224 non-Hispanic Whites and 129 African Americans, all of whom completed an oral glucose tolerance test. Stool microbiome was assessed by whole-metagenome shotgun sequencing with taxonomic profiling. We examined associations among 36 butyrate-producing taxa (n = 7 genera and 29 species) and insulin sensitivity, insulin secretion, disposition index, insulin clearance, and prevalence of dysglycemia (prediabetes plus diabetes, 46% of cohort), adjusting for age, sex, BMI, and race. The genus Coprococcus was associated with higher insulin sensitivity (ß = 0.14; P = 0.002) and disposition index (ß = 0.12; P = 0.012) and a lower rate of dysglycemia (odds ratio [OR] 0.91; 95% CI 0.85-0.97; P = 0.0025). In contrast, Flavonifractor was associated with lower insulin sensitivity (ß = -0.13; P = 0.004) and disposition index (ß = -0.11; P = 0.04) and higher prevalence of dysglycemia (OR 1.22; 95% CI 1.08-1.38; P = 0.0013). Species-level analyses found 10 bacteria associated with beneficial directions of effects and two bacteria with adverse associations on insulin homeostasis and dysglycemia. Although most butyrate producers analyzed appear to be metabolically beneficial, this is not the case for all such bacteria, suggesting that microbiome-directed therapeutic measures to prevent or treat diabetes should be targeted to specific butyrate-producing taxa rather than all butyrate producers.

Cesarean Delivery and Insulin Sensitivity in the Older Adult: The Microbiome and Insulin Longitudinal Evaluation Study.

The present study was designed to evaluate if mode of delivery at birth is associated with body mass index (BMI) and glucose homeostasis traits in later life, controlling for possible confounders, including maternal history of diabetes. Data were obtained through a racially diverse, prospective cohort study of nondiabetic, older adults, the Microbiome and Insulin Longitudinal Evaluation Study (MILES). We used generalized linear models to estimate the association between mode of delivery and glycemic status, BMI (kg/m2), waist circumference (cm), fasting glucose, fasting insulin, insulin secretion, insulin sensitivity, and insulin clearance. Further, we estimated the direct and indirect effects of cesarean delivery on glucose and insulin-related traits, as mediated by BMI status. Relative to vaginal delivery, cesarean delivery was associated with a significantly higher BMI (adjusted beta [aß] 3.53 kg/m2; 95% CI 0.15, 6.91) and fasting glucose (aß 5.12; 95% CI 0.01, 10.23), a 14% decrease in insulin sensitivity (aß -0.14; 95% CI -0.28, -0.01), and a 58% increased risk (adjusted relative risk [aRR] 1.58; 95% CI 1.08, 2.31) for prediabetes/diabetes. Associations were mediated in part by BMI, with the strongest evidence observed for glycemic status (proportion mediated 22.6%; P = .03), fasting insulin (proportion mediated 58.0%; P = .05), and insulin sensitivity index (proportion mediated 45.9%; P = .05). Independent of mediation, a significant direct effect of cesarean delivery on glycemic status was observed (aRR 1.88; 95% CI 1.16, 2.60). Cesarean delivery may lead to reduced insulin sensitivity and, ultimately, increased risk for developing prediabetes and diabetes.

Associations between adherence to the dietary approaches to stop hypertension (DASH) diet and six glucose homeostasis traits in the Microbiome and Insulin Longitudinal Evaluation Study (MILES).

BACKGROUND AND AIMS: The DASH diet conveys protection against type 2 diabetes mellitus (T2D) Via plant-based and non-plant-based recommendations. Research has not identified which glucose homeostasis pathways are improved. We examined associations between adherence to a DASH diet and six glucose homeostasis traits, probing whether associations could be attributed to the plant-based (DASH-P) and/or non-plant based (DASH-NP) components. METHODS AND RESULTS: We included data from 295 adults without T2D (age 59.3 ± 9.00 years; 63.46% non-Hispanic White and 36.54% African American, self-reported race ancestry) participating in the Microbiome and Insulin Longitudinal Evaluation Study (MILES). An oral glucose tolerance test (OGTT) yielded fasting plasma glucose, insulin, C-peptide, and insulin secretion, sensitivity, and disposition index. Habitual dietary intake was assessed by food frequency questionnaire (FFQ). Associations between DASH components and glucose homeostasis traits were examined, controlling for demographics, body mass index (BMI), physical activity, and energy intake. For significant associations, the models were repeated with scores for DASH-P and DASH-NP as predictors in the same model. DASH and DASH-P scores were inversely associated with fasting plasma glucose (DASH:ß = -0.036 ± 0.012,P = 0.005; DASH-P: ß = -0.04 ± 0.017,P = 0.002), and positively associated with insulin sensitivity (DASH:ß = 0.022 ± 0.012,P = 0.042; DASH-P: = 0.036 ± 0.015,P = 0.014). The DASH score was also associated with disposition index (ß = 0.026 ± 0.013,P = 0.038), but this association did not reach significance with DASH-P (ß = 0.035 ± 0.018,P = 0.051). No associations were observed with DASH-NP score (all P > 0.05). CONCLUSIONS: DASH diet is associated with improvement in specific glucose homeostasis traits, likely arising from increased plant-based foods. Such research may help tailor future dietary advice to specific metabolic risk, and to food groups most effective at improving these.

Type 2 Diabetes Partitioned Polygenic Scores Associate With Disease Outcomes in 454,193 Individuals Across 13 Cohorts.

OBJECTIVE: Type 2 diabetes (T2D) has heterogeneous patient clinical characteristics and outcomes. In previous work, we investigated the genetic basis of this heterogeneity by clustering 94 T2D genetic loci using their associations with 47 diabetes-related traits and identified five clusters, termed ß-cell, proinsulin, obesity, lipodystrophy, and liver/lipid. The relationship between these clusters and individual-level metabolic disease outcomes has not been assessed. RESEARCH DESIGN AND METHODS: Here we constructed individual-level partitioned polygenic scores (pPS) for these five clusters in 12 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (n = 454,193) and tested for cross-sectional association with T2D-related outcomes, including blood pressure, renal function, insulin use, age at T2D diagnosis, and coronary artery disease (CAD). RESULTS: Despite all clusters containing T2D risk-increasing alleles, they had differential associations with metabolic outcomes. Increased obesity and lipodystrophy cluster pPS, which had opposite directions of association with measures of adiposity, were both significantly associated with increased blood pressure and hypertension. The lipodystrophy and liver/lipid cluster pPS were each associated with CAD, with increasing and decreasing effects, respectively. An increased liver/lipid cluster pPS was also significantly associated with reduced renal function. The liver/lipid cluster includes known loci linked to liver lipid metabolism (e.g., GCKR, PNPLA3, and TM6SF2), and these findings suggest that cardiovascular disease risk and renal function may be impacted by these loci through their shared disease pathway. CONCLUSIONS: Our findings support that genetically driven pathways leading to T2D also predispose differentially to clinical outcomes.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes.

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

Metabolomic architecture of obesity implicates metabolonic lactone sulfate in cardiometabolic disease.

OBJECTIVE: Identify and characterize circulating metabolite profiles associated with adiposity to inform precision medicine. METHODS: Untargeted plasma metabolomic profiles in the Insulin Resistance Atherosclerosis Family Study (IRASFS) Mexican American cohort (n = 1108) were analyzed for association with anthropometric (body mass index, BMI; waist circumference, WC; waist-to-hip ratio, WHR) and computed tomography measures (visceral adipose tissue, VAT; subcutaneous adipose tissue, SAT; visceral-to-subcutaneous ratio, VSR) of adiposity. Genetic data, inclusive of genome-wide array-based genotyping, whole exome sequencing (WES) and whole genome sequencing (WGS), were evaluated to identify the genetic contributors. Phenotypic and genetic association signals were replicated across ancestries. Transcriptomic data were analyzed to explore the relationship between genetic and metabolomic data. RESULTS: A partially characterized metabolite, tentatively named metabolonic lactone sulfate (X-12063), was consistently associated with BMI, WC, WHR, VAT, and SAT in IRASFS Mexican Americans (PMA <2.02 × 10-27). Trait associations were replicated in IRASFS African Americans (PAA < 1.12 × 10-07). Expanded analyses revealed associations with multiple phenotypic measures of cardiometabolic health, e.g. insulin sensitivity (SI), triglycerides (TG), diastolic blood pressure (DBP) and plasminogen activator inhibitor-1 (PAI-1) in both ancestries. Metabolonic lactone sulfate levels were heritable (h2 > 0.47), and a significant genetic signal at the ZSCAN25/CYP3A5 locus (PMA = 9.00 × 10-41, PAA = 2.31 × 10-10) was observed, highlighting a putative functional variant (rs776746, CYP3A5∗3). Transcriptomic analysis in the African American Genetics of Metabolism and Expression (AAGMEx) cohort supported the association of CYP3A5 with metabolonic lactone sulfate levels (PFDR = 6.64 × 10-07). CONCLUSIONS: Variant rs776746 is associated with a decrease in the transcript levels of CYP3A5, which in turn is associated with increased metabolonic lactone sulfate levels and poor cardiometabolic health.

Sequencing of 640,000 exomes identifies GPR75 variants associated with protection from obesity.

Large-scale human exome sequencing can identify rare protein-coding variants with a large impact on complex traits such as body adiposity. We sequenced the exomes of 645,626 individuals from the United Kingdom, the United States, and Mexico and estimated associations of rare coding variants with body mass index (BMI). We identified 16 genes with an exome-wide significant association with BMI, including those encoding five brain-expressed G protein-coupled receptors (CALCR, MC4R, GIPR, GPR151, and GPR75). Protein-truncating variants in GPR75 were observed in ~4/10,000 sequenced individuals and were associated with 1.8 kilograms per square meter lower BMI and 54% lower odds of obesity in the heterozygous state. Knock out of Gpr75 in mice resulted in resistance to weight gain and improved glycemic control in a high-fat diet model. Inhibition of GPR75 may provide a therapeutic strategy for obesity.

Defining the Relative Role of Insulin Clearance in Early Dysglycemia in Relation to Insulin Sensitivity and Insulin Secretion: The Microbiome and Insulin Longitudinal Evaluation Study (MILES).

Insulin resistance and insufficient insulin secretion are well-recognized contributors to type 2 diabetes. A potential role of reduced insulin clearance has been suggested, but few studies have investigated the contribution of insulin clearance while simultaneously examining decreased insulin sensitivity and secretion. The goal of this study was to conduct such an investigation in a cohort of 353 non-Hispanic White and African American individuals recruited in the Microbiome and Insulin Longitudinal Evaluation Study (MILES). Participants underwent oral glucose tolerance tests from which insulin sensitivity, insulin secretion, insulin clearance, and disposition index were calculated. Regression models examined the individual and joint contributions of these traits to early dysglycemia (prediabetes or newly diagnosed diabetes). In separate models, reduced insulin sensitivity, reduced disposition index, and reduced insulin clearance were associated with dysglycemia. In a joint model, only insulin resistance and reduced insulin secretion were associated with dysglycemia. Models with insulin sensitivity, disposition index, or three insulin traits had the highest discriminative value for dysglycemia (area under the receiver operating characteristics curve of 0.82 to 0.89). These results suggest that in the race groups studied, insulin resistance and compromised insulin secretion are the main independent underlying defects leading to early dysglycemia.

Genome-wide association study of vitamin D concentrations and bone mineral density in the African American-Diabetes Heart Study.

Relative to European Americans, African Americans have lower 25-hydroxyvitamin D (25OHD) and vitamin D binding protein (VDBP) concentrations, higher 1,25-dihydroxyvitamin D (1,25(OH)2D3) concentrations and bone mineral density (BMD), and paradoxically reduced burdens of calcified atherosclerotic plaque (subclinical atherosclerosis). To identify genetic factors contributing to vitamin D and BMD measures, association analysis of >14M variants was conducted in a maximum of 697 African American-Diabetes Heart Study participants with type 2 diabetes (T2D). The most significant association signals were detected for VDBP on chromosome 4; variants rs7041 (ß = 0.44, SE = 0.019, P = 9.4x10-86) and rs4588 (ß = 0.17, SE = 0.021, P = 3.5x10-08) in the group-specific component (vitamin D binding protein) gene (GC). These variants were found to be independently associated. In addition, rs7041 was also associated with bioavailable vitamin D (BAVD; ß = 0.16, SE = 0.02, P = 3.3x10-19). Six rare variants were significantly associated with 25OHD, including a non-synonymous variant in HSPG2 (rs116788687; ß = -1.07, SE = 0.17, P = 2.2x10-10) and an intronic variant in TNIK (rs143555701; ß = -1.01, SE = 0.18, P = 9.0x10-10), both biologically related to bone development. Variants associated with 25OHD failed to replicate in African Americans from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of vitamin D metabolism and bone mineral density phenotypes in an African American population enriched for T2D could provide insight into ethnic specific differences in vitamin D metabolism and bone mineral density.

Multi-omics analysis identifies CpGs near G6PC2 mediating the effects of genetic variants on fasting glucose.

AIMS/HYPOTHESIS: An elevated fasting glucose level in non-diabetic individuals is a key predictor of type 2 diabetes. Genome-wide association studies (GWAS) have identified hundreds of SNPs for fasting glucose but most of their functional roles in influencing the trait are unclear. This study aimed to identify the mediation effects of DNA methylation between SNPs identified as significant from GWAS and fasting glucose using Mendelian randomisation (MR) analyses. METHODS: We first performed GWAS analyses for three cohorts (Taiwan Biobank with 18,122 individuals, the Healthy Aging Longitudinal Study in Taiwan with 1989 individuals and the Stanford Asia-Pacific Program for Hypertension and Insulin Resistance with 416 individuals) with individuals of Han Chinese ancestry in Taiwan, followed by a meta-analysis for combining the three GWAS analysis results to identify significant and independent SNPs for fasting glucose. We determined whether these SNPs were methylation quantitative trait loci (meQTLs) by testing their associations with DNA methylation levels at nearby CpG sites using a subsample of 1775 individuals from the Taiwan Biobank. The MR analysis was performed to identify DNA methylation with causal effects on fasting glucose using meQTLs as instrumental variables based on the 1775 individuals. We also used a two-sample MR strategy to perform replication analysis for CpG sites with significant MR effects based on literature data. RESULTS: Our meta-analysis identified 18 significant (p < 5 × 10-8) and independent SNPs for fasting glucose. Interestingly, all 18 SNPs were meQTLs. The MR analysis identified seven CpGs near the G6PC2 gene that mediated the effects of a significant SNP (rs2232326) in the gene on fasting glucose. The MR effects for two CpGs were replicated using summary data based on the European population, using an exonic SNP rs2232328 in G6PC2 as the instrument. CONCLUSIONS/INTERPRETATION: Our analysis results suggest that rs2232326 and rs2232328 in G6PC2 may affect DNA methylation at CpGs near the gene and that the methylation may have downstream effects on fasting glucose. Therefore, SNPs in G6PC2 and CpGs near G6PC2 may reside along the pathway that influences fasting glucose levels. This is the first study to report CpGs near G6PC2, an important gene for regulating insulin secretion, mediating the effects of GWAS-significant SNPs on fasting glucose.

Allele-specific variation at APOE increases nonalcoholic fatty liver disease and obesity but decreases risk of Alzheimer's disease and myocardial infarction.

Nonalcoholic fatty liver disease (NAFLD) is a leading cause of chronic liver disease and is highly correlated with metabolic disease. NAFLD results from environmental exposures acting on a susceptible polygenic background. This study performed the largest multiethnic investigation of exonic variation associated with NAFLD and correlated metabolic traits and diseases. An exome array meta-analysis was carried out among eight multiethnic population-based cohorts (n = 16 492) with computed tomography (CT) measured hepatic steatosis. A fixed effects meta-analysis identified five exome-wide significant loci (P < 5.30 × 10-7); including a novel signal near TOMM40/APOE. Joint analysis of TOMM40/APOE variants revealed the TOMM40 signal was attributed to APOE rs429358-T; APOE rs7412 was not associated with liver attenuation. Moreover, rs429358-T was associated with higher serum alanine aminotransferase, liver steatosis, cirrhosis, triglycerides and obesity; as well as, lower cholesterol and decreased risk of myocardial infarction and Alzheimer's disease (AD) in phenome-wide association analyses in the Michigan Genomics Initiative, United Kingdom Biobank and/or public datasets. These results implicate APOE in imaging-based identification of NAFLD. This association may or may not translate to nonalcoholic steatohepatitis; however, these results indicate a significant association with advanced liver disease and hepatic cirrhosis. These findings highlight allelic heterogeneity at the APOE locus and demonstrate an inverse link between NAFLD and AD at the exome level in the largest analysis to date.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

A Noncoding Variant Near PPP1R3B Promotes Liver Glycogen Storage and MetS, but Protects Against Myocardial Infarction.

CONTEXT: Glycogen storage diseases are rare. Increased glycogen in the liver results in increased attenuation. OBJECTIVE: Investigate the association and function of a noncoding region associated with liver attenuation but not histologic nonalcoholic fatty liver disease. DESIGN: Genetics of Obesity-associated Liver Disease Consortium. SETTING: Population-based. MAIN OUTCOME: Computed tomography measured liver attenuation. RESULTS: Carriers of rs4841132-A (frequency 2%-19%) do not show increased hepatic steatosis; they have increased liver attenuation indicative of increased glycogen deposition. rs4841132 falls in a noncoding RNA LOC157273 ~190 kb upstream of PPP1R3B. We demonstrate that rs4841132-A increases PPP1R3B through a cis genetic effect. Using CRISPR/Cas9 we engineered a 105-bp deletion including rs4841132-A in human hepatocarcinoma cells that increases PPP1R3B, decreases LOC157273, and increases glycogen perfectly mirroring the human disease. Overexpression of PPP1R3B or knockdown of LOC157273 increased glycogen but did not result in decreased LOC157273 or increased PPP1R3B, respectively, suggesting that the effects may not all occur via affecting RNA levels. Based on electronic health record (EHR) data, rs4841132-A associates with all components of the metabolic syndrome (MetS). However, rs4841132-A associated with decreased low-density lipoprotein (LDL) cholesterol and risk for myocardial infarction (MI). A metabolic signature for rs4841132-A includes increased glycine, lactate, triglycerides, and decreased acetoacetate and beta-hydroxybutyrate. CONCLUSIONS: These results show that rs4841132-A promotes a hepatic glycogen storage disease by increasing PPP1R3B and decreasing LOC157273. rs4841132-A promotes glycogen accumulation and development of MetS but lowers LDL cholesterol and risk for MI. These results suggest that elevated hepatic glycogen is one cause of MetS that does not invariably promote MI.

Genetic variants modulate gene expression statin response in human lymphoblastoid cell lines.

BACKGROUND: Statins are widely prescribed to lower plasma low-density lipoprotein cholesterol levels. Though statins reduce cardiovascular disease risk overall, statin efficacy varies, and some people experience adverse side effects while on statin treatment. Statins also have pleiotropic effects not directly related to their cholesterol-lowering properties, but the mechanisms are not well understood. To identify potential genetic modulators of clinical statin response, we looked for genetic variants associated with statin-induced changes in gene expression (differential eQTLs or deQTLs) in lymphoblastoid cell lines (LCLs) derived from participants of the Cholesterol and Pharmacogenetics (CAP) 40 mg/day 6-week simvastatin clinical trial. We exposed CAP LCLs to 2 µM simvastatin or control buffer for 24 h and performed polyA-selected, strand-specific RNA-seq. Statin-induced changes in gene expression from 259 European ancestry or 153 African American ancestry LCLs were adjusted for potential confounders prior to association with genotyped and imputed genetic variants within 1 Mb of each gene's transcription start site. RESULTS: From the deQTL meta-analysis of the two ancestral populations, we identified significant cis-deQTLs for 15 genes (TBC1D4, MDGA1, CHI3L2, OAS1, GATM, ASNSD1, GLUL, TDRD12, PPIP5K2, OAS3, SERPINB1, ANKDD1A, DTD1, CYFIP2, and GSDME), eight of which were significant in at least one of the ancestry subsets alone. We also conducted eQTL analyses of the endogenous (control-treated), statin-treated, and average of endogenous and statin-treated LCL gene expression levels. We identified eQTLs for approximately 6000 genes in each of the three (endogenous, statin-treated, and average) eQTL meta-analyses, with smaller numbers identified in the ancestral subsets alone. CONCLUSIONS: Several of the genes in which we identified deQTLs have functions in human health and disease, such as defense from viruses, glucose regulation, and response to chemotherapy drugs. This suggests that DNA variation may play a role in statin effects on various health outcomes. These findings could prove useful to future studies aiming to assess benefit versus risk of statin treatment using individual genetic profiles.

Physical Activity Associations with Bone Mineral Density and Modification by Metabolic Traits.

OBJECTIVE: To assess the relationship of physical activity with bone mineral density (BMD) at various sites and examine potential modifying metabolic factors. METHODS: Responses from physical activity questionnaires were used to determine total physical activity (PA), moderate physical activity (mod-PA), and sedentary time. Regression analyses were performed to evaluate association of activity traits with insulin sensitivity by euglycemic clamp, adiponectin, C-reactive protein (CRP), and plasminogen activator inhibitor-1 (PAI-1) in 741 healthy subjects. RESULTS: The cohort was relatively sedentary. Activity level was associated with arm, pelvis, and leg BMD in univariate analyses. In multivariate association analyses of arm BMD, only female sex (ß = -0.73, P < 0.0001) and adiponectin (ß = -0.076, P = 0.0091) were significant. Multivariate analyses of pelvis BMD found independent associations with body mass index (BMI) (ß = 0.33, P < 0.0001), adiponectin (ß = -0.10, P = 0.013), female sex (ß = -0.18, P < 0.0001), sedentary time (ß = -0.088, P = 0.034), PA (ß = 0.11, P = 0.01), and mod-PA (ß = 0.11, P = 0.014). Age (ß = -0.10, P = 0.0087), female sex (ß = -0.63, P < 0.0001), BMI (ß = 0.24, P < 0.0001), and mod-PA (ß = 0.10, P = 0.0024) were independently associated with leg BMD. CONCLUSIONS: These results suggest that BMD increases with physical activity in the arms, legs, and pelvis and is inversely related to sedentary time in the pelvis and legs; these associations may be modified by age, sex, BMI, and adiponectin, depending on the site, with physical activity being more important to pelvis and leg BMD than arm BMD and sedentary time being important for pelvis BMD. Moreover, we demonstrated that CRP, PAI-1, and insulin sensitivity play a minor role in BMD.

Rationale, design and baseline characteristics of the Microbiome and Insulin Longitudinal Evaluation Study (MILES).

AIM: To investigate the role of the gut microbiome in regulating key insulin homeostasis traits (insulin sensitivity, insulin secretion and insulin clearance) whose dysfunction leads to type 2 diabetes (T2D). MATERIALS AND METHODS: The Microbiome and Insulin Longitudinal Evaluation Study (MILES) focuses on African American and non-Hispanic white participants aged 40-80 years without diabetes. Three study visits are planned (at baseline, 15 and 30 months). Baseline measurements include assessment of the stool microbiome and administration of an oral glucose tolerance test, which will yield indexes of insulin sensitivity, insulin secretion and insulin clearance. The gut microbiome profile (composition and function) will be determined using whole metagenome shotgun sequencing along with analyses of plasma short chain fatty acids. Additional data collected include dietary history, sociodemographic factors, health habits, anthropometry, medical history, medications and family history. Most assessments are repeated 15 and 30 months following baseline. RESULTS: After screening 875 individuals, 129 African American and 224 non-Hispanic white participants were enrolled. At baseline, African American participants have higher blood pressure, weight, body mass index, waist and hip circumferences but similar waist-hip ratio compared with the non-Hispanic white participants. On average, African American participants are less insulin-sensitive and have higher acute insulin secretion and lower insulin clearance. CONCLUSIONS: The longitudinal design and robust characterization of potential mediators will allow for the assessment of glucose and insulin homeostasis and gut microbiota as they change over time, improving our ability to discern causal relationships between the microbiome and the insulin homeostasis traits whose deterioration determines T2D, setting the stage for future microbiome-directed therapies to prevent and treat T2D.