Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

The Polygenic and Monogenic Basis of Blood Traits and Diseases.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

A saturated map of common genetic variants associated with human height.

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Primary weight loss failure after Roux-en-Y gastric bypass is characterized by impaired gut-hormone mediated regulation of food intake.

BACKGROUND/OBJECTIVES: After Roux-en-Y gastric bypass (RYGB) a subset of patients never obtain excess BMI loss (EBMIL) > 50% and are categorized as having primary weight loss (WL) failure. We hypothesized that postprandial concentrations of glucagon-like peptide 1 (GLP-1) and peptide YY (PYY) would be lower in patients with primary WL failure compared with patients with successfully maintained WL. Furthermore, that inhibition of gut hormone secretions would increase ad libitum food intake less in patients with primary WL failure. SUBJECTS/METHODS: Twenty women with primary WL failure (LowEBMIL < 50%) were individually matched to twenty women with successful WL (HighEBMIL > 60%) on age, preoperative BMI and time from RYGB. On separate days performed in a random order, patient-blinded subcutaneous injections of octreotide or saline (placebo) were followed by a fixed breakfast and an ad libitum lunch with blood sampling for appetite regulating hormones and Visual-Analogue-Scale (VAS)-scoring of hunger/satiety. Furthermore, participants underwent gene variant analysis for GLP-1, PYY and their receptors, indirect calorimetry, dual-energy X-ray absorptiometry (DXA)-scans, 4-days at-home food registration and 14-days step counting. RESULTS: On placebo days, postprandial GLP-1, PYY and cholecystokinin (CCK) concentrations were similar between groups after breakfast. Fasting ghrelin was lower in LowEBMIL, but the postprandial suppression was similar. LowEBMIL had lower satiety VAS-scores and less suppression of hunger VAS-scores. Gene variants did not differ between groups. Octreotide diminished GLP-1, PYY, CCK and ghrelin concentrations in both groups. Octreotide did not affect ad libitum food intake in LowEBMIL (-1% [-13, 12], mean [95%CI]), while food intake increased in HighEBMIL (+23% [2,44]). CONCLUSIONS: Primary WL failure after RYGB was not characterized by impaired secretions of appetite regulating gut hormones. Interestingly, inhibition of gut hormone secretions with octreotide only increased food intake in patients with successful WL post-RYGB. Thus, an impaired central anorectic response to gut hormones may contribute to primary WL failure after RYGB.

Sequencing and de novo assembly of 150 genomes from Denmark as a population reference.

Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

Genome-wide study of early and severe childhood asthma identifies interaction between CDHR3 and GSDMB.

BACKGROUND: Asthma with severe exacerbation is one of the most common causes of hospitalization among young children. Exacerbations are typically triggered by respiratory infections, but the host factors causing recurrent infections and exacerbations in some children are poorly understood. As a result, current treatment options and preventive measures are inadequate. OBJECTIVE: We sought to identify genetic interaction associated with the development of childhood asthma. METHODS: We performed an exhaustive search for pairwise interaction between genetic single nucleotide polymorphisms using 1204 cases of a specific phenotype of early childhood asthma with severe exacerbations in patients aged 2 to 6 years combined with 5328 nonasthmatic controls. Replication was attempted in 3 independent populations, and potential underlying immune mechanisms were investigated in the COPSAC2010 and COPSAC2000 birth cohorts. RESULTS: We found evidence of interaction, including replication in independent populations, between the known childhood asthma loci CDHR3 and GSDMB. The effect of CDHR3 was dependent on the GSDMB genotype, and this interaction was more pronounced for severe and early onset of disease. Blood immune analyses suggested a mechanism related to increased IL-17A production after viral stimulation. CONCLUSIONS: We found evidence of interaction between CDHR3 and GSDMB in development of early childhood asthma, possibly related to increased IL-17A response to viral infections. This study demonstrates the importance of focusing on specific disease subtypes for understanding the genetic mechanisms of asthma.

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.

A novel rare CUBN variant and three additional genes identified in Europeans with and without diabetes: results from an exome-wide association study of albuminuria.

AIMS/HYPOTHESIS: Identifying rare coding variants associated with albuminuria may open new avenues for preventing chronic kidney disease and end-stage renal disease, which are highly prevalent in individuals with diabetes. Efforts to identify genetic susceptibility variants for albuminuria have so far been limited, with the majority of studies focusing on common variants. METHODS: We performed an exome-wide association study to identify coding variants in a two-stage (discovery and replication) approach. Data from 33,985 individuals of European ancestry (15,872 with and 18,113 without diabetes) and 2605 Greenlanders were included. RESULTS: We identified a rare (minor allele frequency [MAF]: 0.8%) missense (A1690V) variant in CUBN (rs141640975, ß = 0.27, p = 1.3 × 10-11) associated with albuminuria as a continuous measure in the combined European meta-analysis. The presence of each rare allele of the variant was associated with a 6.4% increase in albuminuria. The rare CUBN variant had an effect that was three times stronger in individuals with type 2 diabetes compared with those without (pinteraction = 7.0 × 10-4, ß with diabetes = 0.69, ß without diabetes = 0.20) in the discovery meta-analysis. Gene-aggregate tests based on rare and common variants identified three additional genes associated with albuminuria (HES1, CDC73 and GRM5) after multiple testing correction (pBonferroni < 2.7 × 10-6). CONCLUSIONS/INTERPRETATION: The current study identifies a rare coding variant in the CUBN locus and other potential genes associated with albuminuria in individuals with and without diabetes. These genes have been implicated in renal and cardiovascular dysfunction. The findings provide new insights into the genetic architecture of albuminuria and highlight target genes and pathways for the prevention of diabetes-related kidney disease.

Exome Genotyping Identifies Pleiotropic Variants Associated with Red Blood Cell Traits.

Red blood cell (RBC) traits are important heritable clinical biomarkers and modifiers of disease severity. To identify coding genetic variants associated with these traits, we conducted meta-analyses of seven RBC phenotypes in 130,273 multi-ethnic individuals from studies genotyped on an exome array. After conditional analyses and replication in 27,480 independent individuals, we identified 16 new RBC variants. We found low-frequency missense variants in MAP1A (rs55707100, minor allele frequency [MAF] = 3.3%, p = 2 × 10(-10) for hemoglobin [HGB]) and HNF4A (rs1800961, MAF = 2.4%, p < 3 × 10(-8) for hematocrit [HCT] and HGB). In African Americans, we identified a nonsense variant in CD36 associated with higher RBC distribution width (rs3211938, MAF = 8.7%, p = 7 × 10(-11)) and showed that it is associated with lower CD36 expression and strong allelic imbalance in ex vivo differentiated human erythroblasts. We also identified a rare missense variant in ALAS2 (rs201062903, MAF = 0.2%) associated with lower mean corpuscular volume and mean corpuscular hemoglobin (p < 8 × 10(-9)). Mendelian mutations in ALAS2 are a cause of sideroblastic anemia and erythropoietic protoporphyria. Gene-based testing highlighted three rare missense variants in PKLR, a gene mutated in Mendelian non-spherocytic hemolytic anemia, associated with HGB and HCT (SKAT p < 8 × 10(-7)). These rare, low-frequency, and common RBC variants showed pleiotropy, being also associated with platelet, white blood cell, and lipid traits. Our association results and functional annotation suggest the involvement of new genes in human erythropoiesis. We also confirm that rare and low-frequency variants play a role in the architecture of complex human traits, although their phenotypic effect is generally smaller than originally anticipated.

Carriers of a VEGFA enhancer polymorphism selectively binding CHOP/DDIT3 are predisposed to increased circulating levels of thyroid-stimulating hormone.

BACKGROUND: Levels of serum thyroid-stimulating hormone (TSH) indicate thyroid function, because thyroid hormone negatively controls TSH release. Genetic variants in the vascular endothelial growth factor A (VEGFA) gene are associated with TSH levels. The aim of this study was to characterise the association of VEGFA variants with TSH in a Danish cohort and to identify and characterise functional variants. METHODS: We performed an association study of the VEGFA locus for circulating TSH levels in 8445 Danish individuals. Lead variants were tested for allele-specific effects in vitro using luciferase reporter and gel-shift assays. RESULTS: Four SNPs in VEGFA were associated with circulating TSH (rs9472138, rs881858, rs943080 and rs4711751). For rs881858, the presence of each G-allele was associated with a corresponding decrease in TSH levels of 2.3% (p=8.4×10-9) and an increase in circulating free T4 levels (p=0.0014). The SNP rs881858 is located in a binding site for CHOP (C/EBP homology protein) and c/EBPß (ccaat enhancer binding protein ß). Reporter-gene analysis showed increased basal enhancer activity of the rs881858 A-allele versus the G-allele (34.5±9.9% (average±SEM), p=0.0012), while co-expression of CHOP effectively suppressed the rs881858 A-allele activity. The A-allele showed stronger binding to CHOP in gel-shift assays. CONCLUSIONS: VEGF is an important angiogenic signal required for tissue expansion. We show that VEGFA variation giving allele-specific response to transcription factors with overlapping binding sites associate closely with circulating TSH levels. Because CHOP is induced by several types of intracellular stress, this indicates that cellular stress could be involved in the normal or pathophysiological response of the thyroid to TSH. TRIAL REGISTRATION NUMBER: NCT00289237, NCT00316667; Results.

Identification and functional characterization of G6PC2 coding variants influencing glycemic traits define an effector transcript at the G6PC2-ABCB11 locus.

Genome wide association studies (GWAS) for fasting glucose (FG) and insulin (FI) have identified common variant signals which explain 4.8% and 1.2% of trait variance, respectively. It is hypothesized that low-frequency and rare variants could contribute substantially to unexplained genetic variance. To test this, we analyzed exome-array data from up to 33,231 non-diabetic individuals of European ancestry. We found exome-wide significant (P<5×10-7) evidence for two loci not previously highlighted by common variant GWAS: GLP1R (p.Ala316Thr, minor allele frequency (MAF)=1.5%) influencing FG levels, and URB2 (p.Glu594Val, MAF = 0.1%) influencing FI levels. Coding variant associations can highlight potential effector genes at (non-coding) GWAS signals. At the G6PC2/ABCB11 locus, we identified multiple coding variants in G6PC2 (p.Val219Leu, p.His177Tyr, and p.Tyr207Ser) influencing FG levels, conditionally independent of each other and the non-coding GWAS signal. In vitro assays demonstrate that these associated coding alleles result in reduced protein abundance via proteasomal degradation, establishing G6PC2 as an effector gene at this locus. Reconciliation of single-variant associations and functional effects was only possible when haplotype phase was considered. In contrast to earlier reports suggesting that, paradoxically, glucose-raising alleles at this locus are protective against type 2 diabetes (T2D), the p.Val219Leu G6PC2 variant displayed a modest but directionally consistent association with T2D risk. Coding variant associations for glycemic traits in GWAS signals highlight PCSK1, RREB1, and ZHX3 as likely effector transcripts. These coding variant association signals do not have a major impact on the trait variance explained, but they do provide valuable biological insights.

Glucose tolerance is associated with differential expression of microRNAs in skeletal muscle: results from studies of twins with and without type 2 diabetes.

AIMS/HYPOTHESIS: We aimed to identify microRNAs (miRNAs) associated with type 2 diabetes and risk of developing the disease in skeletal muscle biopsies from phenotypically well-characterised twins. METHODS: We measured muscle miRNA levels in monozygotic (MZ) twins discordant for type 2 diabetes using arrays. Further investigations of selected miRNAs included target prediction, pathway analysis, silencing in cells and association analyses in a separate cohort of 164 non-diabetic MZ and dizygotic twins. The effects of elevated glucose and insulin levels on miRNA expression were examined, and the effect of low birthweight (LBW) was studied in rats. RESULTS: We identified 20 miRNAs that were downregulated in MZ twins with diabetes compared with their non-diabetic co-twins. Differences for members of the miR-15 family (miR-15b and miR-16) were the most statistically significant, and these miRNAs were predicted to influence insulin signalling. Indeed, miR-15b and miR-16 levels were associated with levels of key insulin signalling proteins, miR-15b was associated with the insulin receptor in non-diabetic twins and knockdown of miR-15b/miR-16 in myocytes changed the levels of insulin signalling proteins. LBW in twins and undernutrition during pregnancy in rats were, in contrast to overt type 2 diabetes, associated with increased expression of miR-15b and/or miR-16. Elevated glucose and insulin suppressed miR-16 expression in vitro. CONCLUSIONS: Type 2 diabetes is associated with non-genetic downregulation of several miRNAs in skeletal muscle including miR-15b and miR-16, potentially targeting insulin signalling. The paradoxical findings in twins with overt diabetes and twins at increased risk of the disease underscore the complexity of the regulation of muscle insulin signalling in glucose homeostasis.

Young men with low birthweight exhibit decreased plasticity of genome-wide muscle DNA methylation by high-fat overfeeding.

AIMS/HYPOTHESIS: The association between low birthweight (LBW) and risk of developing type 2 diabetes may involve epigenetic mechanisms, with skeletal muscle being a prime target tissue. Differential DNA methylation patterns have been observed in single genes in muscle tissue from type 2 diabetic and LBW individuals, and we recently showed multiple DNA methylation changes during short-term high-fat overfeeding in muscle of healthy people. In a randomised crossover study, we analysed genome-wide DNA promoter methylation in skeletal muscle of 17 young LBW men and 23 matched normal birthweight (NBW) men after a control and a 5 day high-fat overfeeding diet. METHODS: DNA methylation was measured using Illumina's Infinium BeadArray covering 27,578 CpG sites representing 14,475 different genes. RESULTS: After correction for multiple comparisons, DNA methylation levels were found to be similar in the LBW and NBW groups during the control diet. Whereas widespread DNA methylation changes were observed in the NBW group in response to high-fat overfeeding, only a few methylation changes were seen in the LBW group (χ(2), p < 0.001). CONCLUSIONS/INTERPRETATION: Our results indicate lower DNA methylation plasticity in skeletal muscle from LBW vs NBW men, potentially contributing to understanding the link between LBW and increased risk of type 2 diabetes.

Characterization of genetic variants of GIPR reveals a contribution of ß-arrestin to metabolic phenotypes.

Incretin-based therapies are highly successful in combatting obesity and type 2 diabetes1. Yet both activation and inhibition of the glucose-dependent insulinotropic polypeptide (GIP) receptor (GIPR) in combination with glucagon-like peptide-1 (GLP-1) receptor (GLP-1R) activation have resulted in similar clinical outcomes, as demonstrated by the GIPR-GLP-1R co-agonist tirzepatide2 and AMG-133 (ref. 3) combining GIPR antagonism with GLP-1R agonism. This underlines the importance of a better understanding of the GIP system. Here we show the necessity of ß-arrestin recruitment for GIPR function, by combining in vitro pharmacological characterization of 47 GIPR variants with burden testing of clinical phenotypes and in vivo studies. Burden testing of variants with distinct ligand-binding capacity, Gs activation (cyclic adenosine monophosphate production) and ß-arrestin 2 recruitment and internalization shows that unlike variants solely impaired in Gs signalling, variants impaired in both Gs and ß-arrestin 2 recruitment contribute to lower adiposity-related traits. Endosomal Gs-mediated signalling of the variants shows a ß-arrestin dependency and genetic ablation of ß-arrestin 2 impairs cyclic adenosine monophosphate production and decreases GIP efficacy on glucose control in male mice. This study highlights a crucial impact of ß-arrestins in regulating GIPR signalling and overall preservation of biological activity that may facilitate new developments in therapeutic targeting of the GIPR system.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Physical Activity and Insulin Sensitivity Independently Attenuate the Effect of FTO rs9939609 on Obesity.

OBJECTIVE: The association between FTO rs9939609 and obesity is modified by physical activity (PA) and/or insulin sensitivity (IS). We aimed to assess whether these modifications are independent, to assess whether PA and/or IS modify the association between rs9939609 and cardiometabolic traits, and to elucidate underlying mechanisms. RESEARCH DESIGN AND METHODS: Genetic association analyses comprised up to 19,585 individuals. PA was self-reported, and IS was defined based on inverted HOMA insulin resistance index. Functional analyses were performed in muscle biopsies from 140 men and in cultured muscle cells. RESULTS: The BMI-increasing effect of the FTO rs9939609 A allele was attenuated by 47% with high PA (ß [SE], -0.32 [0.10] kg/m2, P = 0.0013) and by 51% with high IS (-0.31 [0.09] kg/m2, P = 0.00028). Interestingly, these interactions were essentially independent (PA, -0.20 [0.09] kg/m2, P = 0.023; IS, -0.28 [0.09] kg/m2, P = 0.0011). The rs9939609 A allele was also associated with higher all-cause mortality and certain cardiometabolic outcomes (hazard ratio, 1.07-1.20, P > 0.04), and these effects tended to be weakened by greater PA and IS. Moreover, the rs9939609 A allele was associated with higher expression of FTO in skeletal muscle tissue (0.03 [0.01], P = 0.011), and in skeletal muscle cells, we identified a physical interaction between the FTO promoter and an enhancer region encompassing rs9939609. CONCLUSIONS: Greater PA and IS independently reduced the effect of rs9939609 on obesity. These effects might be mediated through altered expression of FTO in skeletal muscle. Our results indicated that PA and/or other means of increasing insulin sensitivity could counteract FTO-related genetic predisposition to obesity.

Rare Heterozygous Loss-of-Function Variants in the Human GLP-1 Receptor Are Not Associated With Cardiometabolic Phenotypes.

CONTEXT: Lost glucagon-like peptide 1 receptor (GLP-1R) function affects human physiology. OBJECTIVE: This work aimed to identify coding nonsynonymous GLP1R variants in Danish individuals to link their in vitro phenotypes and clinical phenotypic associations. METHODS: We sequenced GLP1R in 8642 Danish individuals with type 2 diabetes or normal glucose tolerance and examined the ability of nonsynonymous variants to bind GLP-1 and to signal in transfected cells via cyclic adenosine monophosphate (cAMP) formation and ß-arrestin recruitment. We performed a cross-sectional study between the burden of loss-of-signaling (LoS) variants and cardiometabolic phenotypes in 2930 patients with type 2 diabetes and 5712 participants in a population-based cohort. Furthermore, we studied the association between cardiometabolic phenotypes and the burden of the LoS variants and 60 partly overlapping predicted loss-of-function (pLoF) GLP1R variants found in 330 566 unrelated White exome-sequenced participants in the UK Biobank cohort. RESULTS: We identified 36 nonsynonymous variants in GLP1R, of which 10 had a statistically significant loss in GLP-1-induced cAMP signaling compared to wild-type. However, no association was observed between the LoS variants and type 2 diabetes, although LoS variant carriers had a minor increased fasting plasma glucose level. Moreover, pLoF variants from the UK Biobank also did not reveal substantial cardiometabolic associations, despite a small effect on glycated hemoglobin A1c. CONCLUSION: Since no homozygous LoS nor pLoF variants were identified and heterozygous carriers had similar cardiometabolic phenotype as noncarriers, we conclude that GLP-1R may be of particular importance in human physiology, due to a potential evolutionary intolerance of harmful homozygous GLP1R variants.

European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation.

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Skeletal muscle enhancer interactions identify genes controlling whole-body metabolism.

Obesity and type 2 diabetes (T2D) are metabolic disorders influenced by lifestyle and genetic factors that are characterized by insulin resistance in skeletal muscle, a prominent site of glucose disposal. Numerous genetic variants have been associated with obesity and T2D, of which the majority are located in non-coding DNA regions. This suggests that most variants mediate their effect by altering the activity of gene-regulatory elements, including enhancers. Here, we map skeletal muscle genomic enhancer elements that are dynamically regulated after exposure to the free fatty acid palmitate or the inflammatory cytokine TNFα. By overlapping enhancer positions with the location of disease-associated genetic variants, and resolving long-range chromatin interactions between enhancers and gene promoters, we identify target genes involved in metabolic dysfunction in skeletal muscle. The majority of these genes also associate with altered whole-body metabolic phenotypes in the murine BXD genetic reference population. Thus, our combined genomic investigations identified genes that are involved in skeletal muscle metabolism.