Glycogen metabolism links glucose homeostasis to thermogenesis in adipocytes.

Adipocytes increase energy expenditure in response to prolonged sympathetic activation via persistent expression of uncoupling protein 1 (UCP1)1,2. Here we report that the regulation of glycogen metabolism by catecholamines is critical for UCP1 expression. Chronic ß-adrenergic activation leads to increased glycogen accumulation in adipocytes expressing UCP1. Adipocyte-specific deletion of a scaffolding protein, protein targeting to glycogen (PTG), reduces glycogen levels in beige adipocytes, attenuating UCP1 expression and responsiveness to cold or ß-adrenergic receptor-stimulated weight loss in obese mice. Unexpectedly, we observed that glycogen synthesis and degradation are increased in response to catecholamines, and that glycogen turnover is required to produce reactive oxygen species leading to the activation of p38 MAPK, which drives UCP1 expression. Thus, glycogen has a key regulatory role in adipocytes, linking glucose metabolism to thermogenesis.

Role of Polyunsaturated Fat in Modifying Cardiovascular Risk Associated With Family History of Cardiovascular Disease: Pooled De Novo Results From 15 Observational Studies.

BACKGROUND: It is unknown whether dietary intake of polyunsaturated fatty acids (PUFA) modifies the cardiovascular disease (CVD) risk associated with a family history of CVD. We assessed interactions between biomarkers of low PUFA intake and a family history in relation to long-term CVD risk in a large consortium. METHODS: Blood and tissue PUFA data from 40 885 CVD-free adults were assessed. PUFA levels ≤25th percentile were considered to reflect low intake of linoleic, alpha-linolenic, and eicosapentaenoic/docosahexaenoic acids (EPA/DHA). Family history was defined as having ≥1 first-degree relative who experienced a CVD event. Relative risks with 95% CI of CVD were estimated using Cox regression and meta-analyzed. Interactions were assessed by analyzing product terms and calculating relative excess risk due to interaction. RESULTS: After multivariable adjustments, a significant interaction between low EPA/DHA and family history was observed (product term pooled RR, 1.09 [95% CI, 1.02-1.16]; P=0.01). The pooled relative risk of CVD associated with the combined exposure to low EPA/DHA, and family history was 1.41 (95% CI, 1.30-1.54), whereas it was 1.25 (95% CI, 1.16-1.33) for family history alone and 1.06 (95% CI, 0.98-1.14) for EPA/DHA alone, compared with those with neither exposure. The relative excess risk due to interaction results indicated no interactions. CONCLUSIONS: A significant interaction between biomarkers of low EPA/DHA intake, but not the other PUFA, and a family history was observed. This novel finding might suggest a need to emphasize the benefit of consuming oily fish for individuals with a family history of CVD.

Subcutaneous adipose tissue splice quantitative trait loci reveal differences in isoform usage associated with cardiometabolic traits.

Alternate splicing events can create isoforms that alter gene function, and genetic variants associated with alternate gene isoforms may reveal molecular mechanisms of disease. We used subcutaneous adipose tissue of 426 Finnish men from the METSIM study and identified splice junction quantitative trait loci (sQTLs) for 6,077 splice junctions (FDR < 1%). In the same individuals, we detected expression QTLs (eQTLs) for 59,443 exons and 15,397 genes (FDR < 1%). We identified 595 genes with an sQTL and exon eQTL but no gene eQTL, which could indicate potential isoform differences. Of the significant sQTL signals, 2,114 (39.8%) included at least one proxy variant (linkage disequilibrium r2 > 0.8) located within an intron spanned by the splice junction. We identified 203 sQTLs that colocalized with 141 genome-wide association study (GWAS) signals for cardiometabolic traits, including 25 signals for lipid traits, 24 signals for body mass index (BMI), and 12 signals for waist-hip ratio adjusted for BMI. Among all 141 GWAS signals colocalized with an sQTL, we detected 26 that also colocalized with an exon eQTL for an exon skipped by the sQTL splice junction. At a GWAS signal for high-density lipoprotein cholesterol colocalized with an NR1H3 sQTL splice junction, we show that the alternative splice product encodes an NR1H3 transcription factor that lacks a DNA binding domain and fails to activate transcription. Together, these results detect splicing events and candidate mechanisms that may contribute to gene function at GWAS loci.

Extent to which array genotyping and imputation with large reference panels approximate deep whole-genome sequencing.

Understanding the genetic basis of human diseases and traits is dependent on the identification and accurate genotyping of genetic variants. Deep whole-genome sequencing (WGS), the gold standard technology for SNP and indel identification and genotyping, remains very expensive for most large studies. Here, we quantify the extent to which array genotyping followed by genotype imputation can approximate WGS in studies of individuals of African, Hispanic/Latino, and European ancestry in the US and of Finnish ancestry in Finland (a population isolate). For each study, we performed genotype imputation by using the genetic variants present on the Illumina Core, OmniExpress, MEGA, and Omni 2.5M arrays with the 1000G, HRC, and TOPMed imputation reference panels. Using the Omni 2.5M array and the TOPMed panel, ≥90% of bi-allelic single-nucleotide variants (SNVs) are well imputed (r2 > 0.8) down to minor-allele frequencies (MAFs) of 0.14% in African, 0.11% in Hispanic/Latino, 0.35% in European, and 0.85% in Finnish ancestries. There was little difference in TOPMed-based imputation quality among the arrays with >700k variants. Individual-level imputation quality varied widely between and within the three US studies. Imputation quality also varied across genomic regions, producing regions where even common (MAF > 5%) variants were consistently not well imputed across ancestries. The extent to which array genotyping and imputation can approximate WGS therefore depends on reference panel, genotype array, sample ancestry, and genomic location. Imputation quality by variant or genomic region can be queried with our new tool, RsqBrowser, now deployed on the Michigan Imputation Server.

Integrating transcriptomics, metabolomics, and GWAS helps reveal molecular mechanisms for metabolite levels and disease risk.

Transcriptomics data have been integrated with genome-wide association studies (GWASs) to help understand disease/trait molecular mechanisms. The utility of metabolomics, integrated with transcriptomics and disease GWASs, to understand molecular mechanisms for metabolite levels or diseases has not been thoroughly evaluated. We performed probabilistic transcriptome-wide association and locus-level colocalization analyses to integrate transcriptomics results for 49 tissues in 706 individuals from the GTEx project, metabolomics results for 1,391 plasma metabolites in 6,136 Finnish men from the METSIM study, and GWAS results for 2,861 disease traits in 260,405 Finnish individuals from the FinnGen study. We found that genetic variants that regulate metabolite levels were more likely to influence gene expression and disease risk compared to the ones that do not. Integrating transcriptomics with metabolomics results prioritized 397 genes for 521 metabolites, including 496 previously identified gene-metabolite pairs with strong functional connections and suggested 33.3% of such gene-metabolite pairs shared the same causal variants with genetic associations of gene expression. Integrating transcriptomics and metabolomics individually with FinnGen GWAS results identified 1,597 genes for 790 disease traits. Integrating transcriptomics and metabolomics jointly with FinnGen GWAS results helped pinpoint metabolic pathways from genes to diseases. We identified putative causal effects of UGT1A1/UGT1A4 expression on gallbladder disorders through regulating plasma (E,E)-bilirubin levels, of SLC22A5 expression on nasal polyps and plasma carnitine levels through distinct pathways, and of LIPC expression on age-related macular degeneration through glycerophospholipid metabolic pathways. Our study highlights the power of integrating multiple sets of molecular traits and GWAS results to deepen understanding of disease pathophysiology.

Author Correction: Exome sequencing of Finnish isolates enhances rare-variant association power.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exome sequencing of Finnish isolates enhances rare-variant association power.

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power.

Omega-3 Blood Levels and Stroke Risk: A Pooled and Harmonized Analysis of 183 291 Participants From 29 Prospective Studies.

BACKGROUND: The effect of marine omega-3 PUFAs on risk of stroke remains unclear. METHODS: We investigated the associations between circulating and tissue omega-3 PUFA levels and incident stroke (total, ischemic, and hemorrhagic) in 29 international prospective cohorts. Each site conducted a de novo individual-level analysis using a prespecified analytical protocol with defined exposures, covariates, analytical methods, and outcomes; the harmonized data from the studies were then centrally pooled. Multivariable-adjusted HRs and 95% CIs across omega-3 PUFA quintiles were computed for each stroke outcome. RESULTS: Among 183 291 study participants, there were 10 561 total strokes, 8220 ischemic strokes, and 1142 hemorrhagic strokes recorded over a median of 14.3 years follow-up. For eicosapentaenoic acid, comparing quintile 5 (Q5, highest) with quintile 1 (Q1, lowest), total stroke incidence was 17% lower (HR, 0.83 [CI, 0.76-0.91]; P<0.0001), and ischemic stroke was 18% lower (HR, 0.82 [CI, 0.74-0.91]; P<0.0001). For docosahexaenoic acid, comparing Q5 with Q1, there was a 12% lower incidence of total stroke (HR, 0.88 [CI, 0.81-0.96]; P=0.0001) and a 14% lower incidence of ischemic stroke (HR, 0.86 [CI, 0.78-0.95]; P=0.0001). Neither eicosapentaenoic acid nor docosahexaenoic acid was associated with a risk for hemorrhagic stroke. These associations were not modified by either baseline history of AF or prevalent CVD. CONCLUSIONS: Higher omega-3 PUFA levels are associated with lower risks of total and ischemic stroke but have no association with hemorrhagic stroke.

Association of structural variation with cardiometabolic traits in Finns.

The contribution of genome structural variation (SV) to quantitative traits associated with cardiometabolic diseases remains largely unknown. Here, we present the results of a study examining genetic association between SVs and cardiometabolic traits in the Finnish population. We used sensitive methods to identify and genotype 129,166 high-confidence SVs from deep whole-genome sequencing (WGS) data of 4,848 individuals. We tested the 64,572 common and low-frequency SVs for association with 116 quantitative traits and tested candidate associations using exome sequencing and array genotype data from an additional 15,205 individuals. We discovered 31 genome-wide significant associations at 15 loci, including 2 loci at which SVs have strong phenotypic effects: (1) a deletion of the ALB promoter that is greatly enriched in the Finnish population and causes decreased serum albumin level in carriers (p = 1.47 × 10-54) and is also associated with increased levels of total cholesterol (p = 1.22 × 10-28) and 14 additional cholesterol-related traits, and (2) a multi-allelic copy number variant (CNV) at PDPR that is strongly associated with pyruvate (p = 4.81 × 10-21) and alanine (p = 6.14 × 10-12) levels and resides within a structurally complex genomic region that has accumulated many rearrangements over evolutionary time. We also confirmed six previously reported associations, including five led by stronger signals in single nucleotide variants (SNVs) and one linking recurrent HP gene deletion and cholesterol levels (p = 6.24 × 10-10), which was also found to be strongly associated with increased glycoprotein level (p = 3.53 × 10-35). Our study confirms that integrating SVs in trait-mapping studies will expand our knowledge of genetic factors underlying disease risk.

SARS-CoV-2 Infection of Human Neurons Is TMPRSS2 Independent, Requires Endosomal Cell Entry, and Can Be Blocked by Inhibitors of Host Phosphoinositol-5 Kinase.

2019 coronavirus disease (COVID-19) is a disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In addition to respiratory illness, COVID-19 patients exhibit neurological symptoms lasting from weeks to months (long COVID). It is unclear whether these neurological manifestations are due to an infection of brain cells. We found that a small fraction of human induced pluripotent stem cell (iPSC)-derived neurons, but not astrocytes, were naturally susceptible to SARS-CoV-2. Based on the inhibitory effect of blocking antibodies, the infection seemed to depend on the receptor angiotensin-converting enzyme 2 (ACE2), despite very low levels of its expression in neurons. The presence of double-stranded RNA in the cytoplasm (the hallmark of viral replication), abundant synthesis of viral late genes localized throughout infected cells, and an increase in the level of viral RNA in the culture medium (viral release) within the first 48 h of infection suggested that the infection was productive. Productive entry of SARS-CoV-2 requires the fusion of the viral and cellular membranes, which results in the delivery of the viral genome into the cytoplasm of the target cell. The fusion is triggered by proteolytic cleavage of the viral surface spike protein, which can occur at the plasma membrane or from endosomes or lysosomes. We found that SARS-CoV-2 infection of human neurons was insensitive to nafamostat and camostat, which inhibit cellular serine proteases, including transmembrane serine protease 2 (TMPRSS2). Inhibition of cathepsin L also did not significantly block infection. In contrast, the neuronal infection was blocked by apilimod, an inhibitor of phosphatidyl-inositol 5 kinase (PIK5K), which regulates early to late endosome maturation. IMPORTANCE COVID-19 is a disease caused by the coronavirus SARS-CoV-2. Millions of patients display neurological symptoms, including headache, impairment of memory, seizures, and encephalopathy, as well as anatomical abnormalities, such as changes in brain morphology. SARS-CoV-2 infection of the human brain has been documented, but it is unclear whether the observed neurological symptoms are linked to direct brain infection. The mechanism of virus entry into neurons has also not been characterized. Here, we investigated SARS-CoV-2 infection by using a human iPSC-derived neural cell model and found that a small fraction of cortical-like neurons was naturally susceptible to infection. The productive infection was ACE2 dependent and TMPRSS2 independent. We also found that the virus used the late endosomal and lysosomal pathway for cell entry and that the infection could be blocked by apilimod, an inhibitor of cellular PIK5K.

Healthy dietary pattern is associated with lower glycemia independently of the genetic risk of type 2 diabetes: a cross-sectional study in Finnish men.

PURPOSE: Hyperglycemia is affected by lifestyle and genetic factors. We investigated if dietary patterns associate with glycemia in individuals with high or low genetic risk for type 2 diabetes (T2D). METHODS: Men (n = 1577, 51-81 years) without T2D from the Metabolic Syndrome in Men (METSIM) cohort filled a food-frequency questionnaire and participated in a 2-hour oral glucose tolerance test. Polygenetic risk score (PRS) including 76 genetic variants was used to stratify participants into low or high T2D risk groups. We established two data-driven dietary patterns, termed healthy and unhealthy, and investigated their association with plasma glucose concentrations and hyperglycemia risk. RESULTS: Healthy dietary pattern was associated with lower fasting and 2-hour plasma glucose, glucose area under the curve, and better insulin sensitivity (Matsuda insulin sensitivity index) and insulin secretion (disposition index) in unadjusted and adjusted models, whereas the unhealthy pattern was not. No interaction was observed between the patterns and PRS on glycemic measures. Healthy dietary pattern was negatively associated with the risk for hyperglycemia in an adjusted model (OR 0.69, 95% CI 0.51-0.95, in the highest tertile), whereas unhealthy pattern was not (OR 1.08, 95% CI 0.79-1.47, in the highest tertile). No interaction was found between diet and PRS on the risk for hyperglycemia (p = 0.69 for healthy diet, p = 0.54 for unhealthy diet). CONCLUSION: Our findings suggest that healthy diet is associated with lower glucose concentrations and lower risk for hyperglycemia in men with no interaction with the genetic risk.

Chromatin accessibility and gene expression during adipocyte differentiation identify context-dependent effects at cardiometabolic GWAS loci.

Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and mechanisms at genome-wide association study (GWAS) loci. To identify regulatory elements that display differential activity across adipocyte differentiation, we performed ATAC-seq and RNA-seq in a human cell model of preadipocytes and adipocytes at days 4 and 14 of differentiation. For comparison, we created a consensus map of ATAC-seq peaks in 11 human subcutaneous adipose tissue samples. We identified 58,387 context-dependent chromatin accessibility peaks and 3,090 context-dependent genes between all timepoint comparisons (log2 fold change>1, FDR<5%) with 15,919 adipocyte- and 18,244 preadipocyte-dependent peaks. Adipocyte-dependent peaks showed increased overlap (60.1%) with Roadmap Epigenomics adipocyte nuclei enhancers compared to preadipocyte-dependent peaks (11.5%). We linked context-dependent peaks to genes based on adipocyte promoter capture Hi-C data, overlap with adipose eQTL variants, and context-dependent gene expression. Of 16,167 context-dependent peaks linked to a gene, 5,145 were linked by two or more strategies to 1,670 genes. Among GWAS loci for cardiometabolic traits, adipocyte-dependent peaks, but not preadipocyte-dependent peaks, showed significant enrichment (LD score regression P<0.005) for waist-to-hip ratio and modest enrichment (P < 0.05) for HDL-cholesterol. We identified 659 peaks linked to 503 genes by two or more approaches and overlapping a GWAS signal, suggesting a regulatory mechanism at these loci. To identify variants that may alter chromatin accessibility between timepoints, we identified 582 variants in 454 context-dependent peaks that demonstrated allelic imbalance in accessibility (FDR<5%), of which 55 peaks also overlapped GWAS variants. At one GWAS locus for palmitoleic acid, rs603424 was located in an adipocyte-dependent peak linked to SCD and exhibited allelic differences in transcriptional activity in adipocytes (P = 0.003) but not preadipocytes (P = 0.09). These results demonstrate that context-dependent peaks and genes can guide discovery of regulatory variants at GWAS loci and aid identification of regulatory mechanisms.

Genetic evidence that high BMI in childhood has a protective effect on intermediate diabetes traits, including measures of insulin sensitivity and secretion, after accounting for BMI in adulthood.

AIMS/HYPOTHESIS: Determining how high BMI at different time points influences the risk of developing type 2 diabetes and affects insulin secretion and insulin sensitivity is critical. METHODS: By estimating childhood BMI in 441,761 individuals in the UK Biobank, we identified which genetic variants had larger effects on adulthood BMI than on childhood BMI, and vice versa. All genome-wide significant genetic variants were then used to separate the independent genetic effects of high childhood BMI from those of high adulthood BMI on the risk of type 2 diabetes and insulin-related phenotypes using Mendelian randomisation. We performed two-sample MR using external studies of type 2 diabetes, and oral and intravenous measures of insulin secretion and sensitivity. RESULTS: We found that a childhood BMI that was one standard deviation (1.97 kg/m2) higher than the mean, corrected for the independent genetic liability to adulthood BMI, was associated with a protective effect for seven measures of insulin sensitivity and secretion, including increased insulin sensitivity index (ß=0.15; 95% CI 0.067, 0.225; p=2.79×10-4) and reduced fasting glucose levels (ß=-0.053; 95% CI -0.089, -0.017; p=4.31×10-3). However, there was little to no evidence of a direct protective effect on type 2 diabetes (OR 0.94; 95% CI 0.85, 1.04; p=0.228) independently of genetic liability to adulthood BMI. CONCLUSIONS/INTERPRETATION: Our results provide evidence of the protective effect of higher childhood BMI on insulin secretion and sensitivity, which are crucial intermediate diabetes traits. However, we stress that our results should not currently lead to any change in public health or clinical practice, given the uncertainty regarding the biological pathway of these effects and the limitations of this type of study.

Protective Alzheimer's disease-associated APP A673T variant predominantly decreases sAPPß levels in cerebrospinal fluid and 2D/3D cell culture models.

The rare A673T variant was the first variant found within the amyloid precursor protein (APP) gene conferring protection against Alzheimer's disease (AD). Thereafter, different studies have discovered that the carriers of the APP A673T variant show reduced levels of amyloid beta (Aß) in the plasma and better cognitive performance at high age. Here, we analyzed cerebrospinal fluid (CSF) and plasma of APP A673T carriers and control individuals using a mass spectrometry-based proteomics approach to identify differentially regulated targets in an unbiased manner. Furthermore, the APP A673T variant was introduced into 2D and 3D neuronal cell culture models together with the pathogenic APP Swedish and London mutations. Consequently, we now report for the first time the protective effects of the APP A673T variant against AD-related alterations in the CSF, plasma, and brain biopsy samples from the frontal cortex. The CSF levels of soluble APPß (sAPPß) and Aß42 were significantly decreased on average 9-26% among three APP A673T carriers as compared to three well-matched controls not carrying the protective variant. Consistent with these CSF findings, immunohistochemical assessment of cortical biopsy samples from the same APP A673T carriers did not reveal Aß, phospho-tau, or p62 pathologies. We identified differentially regulated targets involved in protein phosphorylation, inflammation, and mitochondrial function in the CSF and plasma samples of APP A673T carriers. Some of the identified targets showed inverse levels in AD brain tissue with respect to increased AD-associated neurofibrillary pathology. In 2D and 3D neuronal cell culture models expressing APP with the Swedish and London mutations, the introduction of the APP A673T variant resulted in lower sAPPß levels. Concomitantly, the levels of sAPPα were increased, while decreased levels of CTFß and Aß42 were detected in some of these models. Our findings emphasize the important role of APP-derived peptides in the pathogenesis of AD and demonstrate the effectiveness of the protective APP A673T variant to shift APP processing towards the non-amyloidogenic pathway in vitro even in the presence of two pathogenic mutations.

Downregulation of diacylglycerol kinase delta contributes to hyperglycemia-induced insulin resistance.

Type 2 (non-insulin-dependent) diabetes mellitus is a progressive metabolic disorder arising from genetic and environmental factors that impair beta cell function and insulin action in peripheral tissues. We identified reduced diacylglycerol kinase delta (DGKdelta) expression and DGK activity in skeletal muscle from type 2 diabetic patients. In diabetic animals, reduced DGKdelta protein and DGK kinase activity were restored upon correction of glycemia. DGKdelta haploinsufficiency increased diacylglycerol content, reduced peripheral insulin sensitivity, insulin signaling, and glucose transport, and led to age-dependent obesity. Metabolic flexibility, evident by the transition between lipid and carbohydrate utilization during fasted and fed conditions, was impaired in DGKdelta haploinsufficient mice. We reveal a previously unrecognized role for DGKdelta in contributing to hyperglycemia-induced peripheral insulin resistance and thereby exacerbating the severity of type 2 diabetes. DGKdelta deficiency causes peripheral insulin resistance and metabolic inflexibility. These defects in glucose and energy homeostasis contribute to mild obesity later in life.

The causal effect of obesity on prediabetes and insulin resistance reveals the important role of adipose tissue in insulin resistance.

Reverse causality has made it difficult to establish the causal directions between obesity and prediabetes and obesity and insulin resistance. To disentangle whether obesity causally drives prediabetes and insulin resistance already in non-diabetic individuals, we utilized the UK Biobank and METSIM cohort to perform a Mendelian randomization (MR) analyses in the non-diabetic individuals. Our results suggest that both prediabetes and systemic insulin resistance are caused by obesity (p = 1.2×10-3 and p = 3.1×10-24). As obesity reflects the amount of body fat, we next studied how adipose tissue affects insulin resistance. We performed both bulk RNA-sequencing and single nucleus RNA sequencing on frozen human subcutaneous adipose biopsies to assess adipose cell-type heterogeneity and mitochondrial (MT) gene expression in insulin resistance. We discovered that the adipose MT gene expression and body fat percent are both independently associated with insulin resistance (p≤0.05 for each) when adjusting for the decomposed adipose cell-type proportions. Next, we showed that these 3 factors, adipose MT gene expression, body fat percent, and adipose cell types, explain a substantial amount (44.39%) of variance in insulin resistance and can be used to predict it (p≤2.64×10-5 in 3 independent human cohorts). In summary, we demonstrated that obesity is a strong determinant of both prediabetes and insulin resistance, and discovered that individuals' adipose cell-type composition, adipose MT gene expression, and body fat percent predict their insulin resistance, emphasizing the critical role of adipose tissue in systemic insulin resistance.

Adiponectin GWAS loci harboring extensive allelic heterogeneity exhibit distinct molecular consequences.

Loci identified in genome-wide association studies (GWAS) can include multiple distinct association signals. We sought to identify the molecular basis of multiple association signals for adiponectin, a hormone involved in glucose regulation secreted almost exclusively from adipose tissue, identified in the Metabolic Syndrome in Men (METSIM) study. With GWAS data for 9,262 men, four loci were significantly associated with adiponectin: ADIPOQ, CDH13, IRS1, and PBRM1. We performed stepwise conditional analyses to identify distinct association signals, a subset of which are also nearly independent (lead variant pairwise r2<0.01). Two loci exhibited allelic heterogeneity, ADIPOQ and CDH13. Of seven association signals at the ADIPOQ locus, two signals colocalized with adipose tissue expression quantitative trait loci (eQTLs) for three transcripts: trait-increasing alleles at one signal were associated with increased ADIPOQ and LINC02043, while trait-increasing alleles at the other signal were associated with decreased ADIPOQ-AS1. In reporter assays, adiponectin-increasing alleles at two signals showed corresponding directions of effect on transcriptional activity. Putative mechanisms for the seven ADIPOQ signals include a missense variant (ADIPOQ G90S), a splice variant, a promoter variant, and four enhancer variants. Of two association signals at the CDH13 locus, the first signal consisted of promoter variants, including the lead adipose tissue eQTL variant for CDH13, while a second signal included a distal intron 1 enhancer variant that showed ~2-fold allelic differences in transcriptional reporter activity. Fine-mapping and experimental validation demonstrated that multiple, distinct association signals at these loci can influence multiple transcripts through multiple molecular mechanisms.

Adipose Tissue Gene Expression Associations Reveal Hundreds of Candidate Genes for Cardiometabolic Traits.

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with cardiometabolic traits including type 2 diabetes (T2D), lipid levels, body fat distribution, and adiposity, although most causal genes remain unknown. We used subcutaneous adipose tissue RNA-seq data from 434 Finnish men from the METSIM study to identify 9,687 primary and 2,785 secondary cis-expression quantitative trait loci (eQTL; <1 Mb from TSS, FDR < 1%). Compared to primary eQTL signals, secondary eQTL signals were located further from transcription start sites, had smaller effect sizes, and were less enriched in adipose tissue regulatory elements compared to primary signals. Among 2,843 cardiometabolic GWAS signals, 262 colocalized by LD and conditional analysis with 318 transcripts as primary and conditionally distinct secondary cis-eQTLs, including some across ancestries. Of cardiometabolic traits examined for adipose tissue eQTL colocalizations, waist-hip ratio (WHR) and circulating lipid traits had the highest percentage of colocalized eQTLs (15% and 14%, respectively). Among alleles associated with increased cardiometabolic GWAS risk, approximately half (53%) were associated with decreased gene expression level. Mediation analyses of colocalized genes and cardiometabolic traits within the 434 individuals provided further evidence that gene expression influences variant-trait associations. These results identify hundreds of candidate genes that may act in adipose tissue to influence cardiometabolic traits.

Exome-Derived Adiponectin-Associated Variants Implicate Obesity and Lipid Biology.

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.