The association of cardiometabolic, diet and lifestyle parameters with plasma glucagon-like peptide-1: An IMI DIRECT study.

CONTEXT: The role of glucagon-like peptide-1(GLP-1) in Type 2 diabetes (T2D) and obesity is not fully understood. OBJECTIVE: We investigate the association of cardiometabolic, diet and lifestyle parameters on fasting and postprandial GLP-1 in people at risk of, or living with, T2D. METHOD: We analysed cross-sectional data from the two Innovative Medicines Initiative (IMI) Diabetes Research on Patient Stratification (DIRECT) cohorts, cohort 1(n=2127) individuals at risk of diabetes; cohort 2 (n=789) individuals with new-onset of T2D. RESULTS: Our multiple regression analysis reveals that fasting total GLP-1 is associated with an insulin resistant phenotype and observe a strong independent relationship with male sex, increased adiposity and liver fat particularly in the prediabetes population. In contrast, we showed that incremental GLP-1 decreases with worsening glycaemia, higher adiposity, liver fat, male sex and reduced insulin sensitivity in the prediabetes cohort. Higher fasting total GLP-1 was associated with a low intake of wholegrain, fruit and vegetables inpeople with prediabetes, and with a high intake of red meat and alcohol in people with diabetes. CONCLUSION: These studies provide novel insights into the association between fasting and incremental GLP-1, metabolic traits of diabetes and obesity, and dietary intake and raise intriguing questions regarding the relevance of fasting GLP-1 in the pathophysiology T2D.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits.

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

A meta-analysis and a functional study support the influence of mtDNA variant m.16519C on the risk of rapid progression of knee osteoarthritis.

OBJECTIVES: To identify mitochondrial DNA (mtDNA) genetic variants associated with the risk of rapid progression of knee osteoarthritis (OA) and to characterise their functional significance using a cellular model of transmitochondrial cybrids. METHODS: Three prospective cohorts contributed participants. The osteoarthritis initiative (OAI) included 1095 subjects, the Cohort Hip and Cohort Knee included 373 and 326 came from the PROspective Cohort of Osteoarthritis from A Coruña. mtDNA variants were screened in an initial subset of 450 subjects from the OAI by in-depth sequencing of mtDNA. A meta-analysis of the three cohorts was performed. A model of cybrids was constructed to study the functional consequences of harbouring the risk mtDNA variant by assessing: mtDNA copy number, mitochondrial biosynthesis, mitochondrial fission and fusion, mitochondrial reactive oxygen species (ROS), oxidative stress, autophagy and a whole transcriptome analysis by RNA-sequencing. RESULTS: mtDNA variant m.16519C is over-represented in rapid progressors (combined OR 1.546; 95% CI 1.163 to 2.054; p=0.0027). Cybrids with this variant show increased mtDNA copy number and decreased mitochondrial biosynthesis; they produce higher amounts of mitochondrial ROS, are less resistant to oxidative stress, show a lower expression of the mitochondrial fission-related gene fission mitochondrial 1 and an impairment of autophagic flux. In addition, its presence modulates the transcriptome of cybrids, especially in terms of inflammation, where interleukin 6 emerges as one of the most differentially expressed genes. CONCLUSIONS: The presence of the mtDNA variant m.16519C increases the risk of rapid progression of knee OA. Among the most modulated biological processes associated with this variant, inflammation and negative regulation of cellular process stand out. The design of therapies based on the maintenance of mitochondrial function is recommended.

A phenome-wide comparative analysis of genetic discordance between obesity and type 2 diabetes.

Obesity and type 2 diabetes are causally related, yet there is considerable heterogeneity in the consequences of both conditions and the mechanisms of action are poorly defined. Here we show a genetic-driven approach defining two obesity profiles that convey highly concordant and discordant diabetogenic effects. We annotate and then compare association signals for these profiles across clinical and molecular phenotypic layers. Key differences are identified in a wide range of traits, including cardiovascular mortality, fat distribution, liver metabolism, blood pressure, specific lipid fractions and blood levels of proteins involved in extracellular matrix remodelling. We find marginal differences in abundance of Bacteroidetes and Firmicutes bacteria in the gut. Instrumental analyses reveal prominent causal roles for waist-to-hip ratio, blood pressure and cholesterol content of high-density lipoprotein particles in the development of diabetes in obesity. We prioritize 17 genes from the discordant signature that convey protection against type 2 diabetes in obesity, which may represent logical targets for precision medicine approaches.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries.

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

Impact of insufficient sleep on dysregulated blood glucose control under standardised meal conditions.

AIMS/HYPOTHESIS: Sleep, diet and exercise are fundamental to metabolic homeostasis. In this secondary analysis of a repeated measures, nutritional intervention study, we tested whether an individual's sleep quality, duration and timing impact glycaemic response to a breakfast meal the following morning. METHODS: Healthy adults' data (N = 953 [41% twins]) were analysed from the PREDICT dietary intervention trial. Participants consumed isoenergetic standardised meals over 2 weeks in the clinic and at home. Actigraphy was used to assess sleep variables (duration, efficiency, timing) and continuous glucose monitors were used to measure glycaemic variation (>8000 meals). RESULTS: Sleep variables were significantly associated with postprandial glycaemic control (2 h incremental AUC), at both between- and within-person levels. Sleep period time interacted with meal type, with a smaller effect of poor sleep on postprandial blood glucose levels when high-carbohydrate (low fat/protein) (pinteraction = 0.02) and high-fat (pinteraction = 0.03) breakfasts were consumed compared with a reference 75 g OGTT. Within-person sleep period time had a similar interaction (high carbohydrate: pinteraction = 0.001, high fat: pinteraction = 0.02). Within- and between-person sleep efficiency were significantly associated with lower postprandial blood glucose levels irrespective of meal type (both p < 0.03). Later sleep midpoint (time deviation from midnight) was found to be significantly associated with higher postprandial glucose, in both between-person and within-person comparisons (p = 0.035 and p = 0.051, respectively). CONCLUSIONS/INTERPRETATION: Poor sleep efficiency and later bedtime routines are associated with more pronounced postprandial glycaemic responses to breakfast the following morning. A person's deviation from their usual sleep pattern was also associated with poorer postprandial glycaemic control. These findings underscore sleep as a modifiable, non-pharmacological therapeutic target for the optimal regulation of human metabolic health. Trial registration ClinicalTrials.gov NCT03479866.

A prospective study of the relationships between movement and glycemic control during day and night in pregnancy.

Both disturbed sleep and lack of exercise can disrupt metabolism in pregnancy. Accelerometery was used to objectively assess movement during waking (physical activity) and movement during sleeping (sleep disturbance) periods and evaluated relationships with continuous blood glucose variation during pregnancy. Data was analysed prospectively. 15-women without pre-existing diabetes mellitus wore continuous glucose monitors and triaxial accelerometers from February through June 2018 in Sweden. The relationships between physical activity and sleep disturbance with blood glucose rate of change were assessed. An interaction term was fitted to determine difference in the relationship between movement and glucose variation, conditional on waking/sleeping. Total movement was inversely related to glucose rate of change (p < 0.001, 95% CI (- 0.037, - 0.026)). Stratified analyses showed total physical activity was inversely related to glucose rate of change (p < 0.001, 95% CI (- 0.040, - 0.028)), whereas sleep disturbance was not related to glucose rate of change (p = 0.07, 95% CI (< - 0.001, 0.013)). The interaction term was positively related to glucose rate of change (p < 0.001, 95% CI (0.029, 0.047)). This study provides temporal evidence of a relationship between total movement and glycemic control in pregnancy, which is conditional on time of day. Movement is beneficially related with glycemic control while awake, but not during sleep.

Validation of Polymorphisms Associated with the Risk of Radiation-Induced Oesophagitis in an Independent Cohort of Non-Small-Cell Lung Cancer Patients.

Several studies have identified single-nucleotide polymorphisms (SNPs) associated with adverse effects in non-small-cell lung cancer (NSCLC) patients treated with radiation therapy. Here, using an independent cohort, we aimed to validate the reported associations. We selected 23 SNPs in 17 genes previously associated with radiation-induced oesophagitis for validation in a cohort of 178 Spanish NSCLC patients. Of them, 18 SNPs were finally analysed, following the methods described in the original published studies. Two SNPs replicated their association with radiation-induced oesophagitis (rs7165790 located in the BLM gene: odds ratio (OR) = 0.16, 95% CI = 0.04-0.65, p-value = 0.010; rs4772468 at FGF14: OR = 4.36, 95% CI = 1.15-16.46, p-value = 0.029). The SNP rs2868371 at HSPB1 was also validated but displayed an opposite effect to the formerly described (OR = 3.72; 95% CI = 1.49-9.25; p-value = 0.004). Additionally, we tested a meta-analytic approach including our results and the previous datasets reported in the referenced publications. Twelve SNPs (including the two previously validated) retained their statistically significant association with radiation-induced oesophagitis. This study strengthens the role of inflammation and DNA double-strand break repair pathways in the risk prediction of developing radiation-induced oesophagitis in NSCLC patients. The validated variants are good candidates to be evaluated in risk prediction models for patient stratification based on their radiation susceptibility.

Genetic Predisposition to Coronary Artery Disease in Type 2 Diabetes Mellitus.

BACKGROUND: Coronary artery disease (CAD) is accelerated in subjects with type 2 diabetes mellitus (T2D). METHODS: To test whether this reflects differential genetic influences on CAD risk in subjects with T2D, we performed a systematic assessment of genetic overlap between CAD and T2D in 66 643 subjects (27 708 with CAD and 24 259 with T2D). Variants showing apparent association with CAD in stratified analyses or evidence of interaction were evaluated in a further 117 787 subjects (16 694 with CAD and 11 537 with T2D). RESULTS: None of the previously characterized CAD loci was found to have specific effects on CAD in T2D individuals, and a genome-wide interaction analysis found no new variants for CAD that could be considered T2D specific. When we considered the overall genetic correlations between CAD and its risk factors, we found no substantial differences in these relationships by T2D background. CONCLUSIONS: This study found no evidence that the genetic architecture of CAD differs in those with T2D compared with those without T2D.

A Multi-omic Integrative Scheme Characterizes Tissues of Action at Loci Associated with Type 2 Diabetes.

Resolving the molecular processes that mediate genetic risk remains a challenge because most disease-associated variants are non-coding and functional characterization of these signals requires knowledge of the specific tissues and cell-types in which they operate. To address this challenge, we developed a framework for integrating tissue-specific gene expression and epigenomic maps to obtain "tissue-of-action" (TOA) scores for each association signal by systematically partitioning posterior probabilities from Bayesian fine-mapping. We applied this scheme to credible set variants for 380 association signals from a recent GWAS meta-analysis of type 2 diabetes (T2D) in Europeans. The resulting tissue profiles underscored a predominant role for pancreatic islets and, to a lesser extent, adipose and liver, particularly among signals with greater fine-mapping resolution. We incorporated resulting TOA scores into a rule-based classifier and validated the tissue assignments through comparison with data from cis-eQTL enrichment, functional fine-mapping, RNA co-expression, and patterns of physiological association. In addition to implicating signals with a single TOA, we found evidence for signals with shared effects in multiple tissues as well as distinct tissue profiles between independent signals within heterogeneous loci. Lastly, we demonstrated that TOA scores can be directly coupled with eQTL colocalization to further resolve effector transcripts at T2D signals. This framework guides mechanistic inference by directing functional validation studies to the most relevant tissues and can gain power as fine-mapping resolution and cell-specific annotations become richer. This method is generalizable to all complex traits with relevant annotation data and is made available as an R package.

Low Vitamin D Levels and Frailty Status in Older Adults: A Systematic Review and Meta-Analysis.

Serum vitamin D deficiency is widespread among older adults and is a potential modifiable risk factor for frailty. Moreover, frailty has been suggested as an intermediate step in the association between low levels of vitamin D and mortality. Hence, we conducted a systematic review of the literature and meta-analysis to test the possible association of low concentrations of serum 25-hydroxyvitamin D (25(OH)D), a marker of vitamin D status, with frailty in later life. We reviewed cross-sectional or longitudinal studies evaluating populations of older adults and identifying frailty by a currently validated scale. Meta-analyses were restricted to cross-sectional data from studies using Fried's phenotype to identify frailty. Twenty-six studies were considered in the qualitative synthesis, and thirteen studies were included in the meta-analyses. Quantitative analyses showed significant differences in the comparisons of frail (standardized mean difference (SMD)-1.31, 95% confidence interval (CI) (-2.47, -0.15), p = 0.0271) and pre-frail (SMD-0.79, 95% CI (-1.58, -0.003), p = 0.0491) subjects vs. non-frail subjects. Sensitivity analyses reduced heterogeneity, resulting in a smaller but still highly significant between-groups difference. Results obtained indicate that lower 25(OH)D levels are significantly associated with increasing frailty severity. Future challenges include interventional studies testing the possible benefits of vitamin D supplementation in older adults to prevent/palliate frailty and its associated outcomes.

Association of inflammatory mediators with frailty status in older adults: results from a systematic review and meta-analysis.

Frailty is a geriatric syndrome defined as a status of extreme vulnerability to stressors, leading to a higher risk of negative health-related outcomes. "Inflammaging", an age-related state of low-grade chronic inflammation, is characterized by an increased concentration of pro-inflammatory cytokines and acute phase proteins. Inflammaging has been postulated as an underlying mechanism of frailty, and several studies tested the relationship between frailty and concentration of inflammatory mediators. The aim of this systematic review and meta-analysis was to test whether inflammatory mediators are overproduced in frail older adults. Among the 758 articles identified in the literature search, 50 were included in the systematic review, and 39 in the three meta-analyses, i.e., C-reactive protein (CRP), interleukin 6 (IL6), and tumor necrosis factor α. To reduce heterogeneity, meta-analyses were restricted to studies identifying frailty by the Fried et al. [1] [J. Gerontol. A. Biol. Sci. Med. Sci. 56, M146-56] phenotypic criteria. Quantitative analyses measuring the association between frailty and biomarker concentrations showed significant differences when frail subjects were compared to non-frail and pre-frail subjects for CRP and IL6. This work established strong association between inflammatory biomarkers and frailty, confirming the role of age-related chronic inflammation in frailty development.

Correction: Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

[This corrects the article DOI: 10.1371/journal.pgen.1007813.].

Association analyses identify 31 new risk loci for colorectal cancer susceptibility.

Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide, and has a strong heritable basis. We report a genome-wide association analysis of 34,627 CRC cases and 71,379 controls of European ancestry that identifies SNPs at 31 new CRC risk loci. We also identify eight independent risk SNPs at the new and previously reported European CRC loci, and a further nine CRC SNPs at loci previously only identified in Asian populations. We use in situ promoter capture Hi-C (CHi-C), gene expression, and in silico annotation methods to identify likely target genes of CRC SNPs. Whilst these new SNP associations implicate target genes that are enriched for known CRC pathways such as Wnt and BMP, they also highlight novel pathways with no prior links to colorectal tumourigenesis. These findings provide further insight into CRC susceptibility and enhance the prospects of applying genetic risk scores to personalised screening and prevention.

Developing a network view of type 2 diabetes risk pathways through integration of genetic, genomic and functional data.

BACKGROUND: Genome-wide association studies (GWAS) have identified several hundred susceptibility loci for type 2 diabetes (T2D). One critical, but unresolved, issue concerns the extent to which the mechanisms through which these diverse signals influencing T2D predisposition converge on a limited set of biological processes. However, the causal variants identified by GWAS mostly fall into a non-coding sequence, complicating the task of defining the effector transcripts through which they operate. METHODS: Here, we describe implementation of an analytical pipeline to address this question. First, we integrate multiple sources of genetic, genomic and biological data to assign positional candidacy scores to the genes that map to T2D GWAS signals. Second, we introduce genes with high scores as seeds within a network optimization algorithm (the asymmetric prize-collecting Steiner tree approach) which uses external, experimentally confirmed protein-protein interaction (PPI) data to generate high-confidence sub-networks. Third, we use GWAS data to test the T2D association enrichment of the "non-seed" proteins introduced into the network, as a measure of the overall functional connectivity of the network. RESULTS: We find (a) non-seed proteins in the T2D protein-interaction network so generated (comprising 705 nodes) are enriched for association to T2D (p = 0.0014) but not control traits, (b) stronger T2D-enrichment for islets than other tissues when we use RNA expression data to generate tissue-specific PPI networks and (c) enhanced enrichment (p = 3.9 × 10- 5) when we combine the analysis of the islet-specific PPI network with a focus on the subset of T2D GWAS loci which act through defective insulin secretion. CONCLUSIONS: These analyses reveal a pattern of non-random functional connectivity between candidate causal genes at T2D GWAS loci and highlight the products of genes including YWHAG, SMAD4 or CDK2 as potential contributors to T2D-relevant islet dysfunction. The approach we describe can be applied to other complex genetic and genomic datasets, facilitating integration of diverse data types into disease-associated networks.

Differential Association of Mitochondrial DNA Haplogroups J and H With the Methylation Status of Articular Cartilage: Potential Role in Apoptosis and Metabolic and Developmental Processes.

OBJECTIVE: To analyze the influence of mitochondrial genome variation on the DNA methylome of articular cartilage. METHODS: DNA methylation profiling was performed using data deposited in the NCBI Gene Expression Omnibus database (accession no. GSE43269). Data were obtained for 14 cartilage samples from subjects with haplogroup J and 20 cartilage samples from subjects with haplogroup H. Subsequent validation was performed in an independent subset of 7 subjects with haplogroup J and 9 with haplogroup H by RNA-seq. Correlated genes were validated by real-time polymerase chain reaction in an independent cohort of 12 subjects with haplogroup J and 12 with haplogroup H. Appropriate analyses were performed using R Bioconductor and qBasePlus software, and gene ontology analysis was conducted using DAVID version 6.8. RESULTS: DNA methylation profiling revealed 538 differentially methylated loci, while whole-transcriptome profiling identified 2,384 differentially expressed genes, between cartilage samples from subjects with haplogroup H and those with haplogroup J. Seventeen genes showed an inverse correlation between methylation and expression. In terms of gene ontology, differences in correlations between methylation and expression were also detected between cartilage from subjects with haplogroup H and those with haplogroup J, highlighting a significantly enhanced apoptotic process in cartilage from subjects with haplogroup H (P = 0.007 for methylation and P = 0.019 for expression) and repressed apoptotic process in cartilage from subjects with haplogroup J (P = 0.021 for methylation), as well as a significant enrichment of genes related to metabolic processes (P = 1.93 × 10-4 for methylation and P = 6.79 x 10-4 for expression) and regulation of gene expression (P = 0.012 for methylation) in cartilage from subjects with haplogroup H, and to developmental processes (P = 0.015 for methylation and P = 8.25 x 10-12 for expression) in cartilage from subjects with haplogroup J. CONCLUSION: Mitochondrial DNA variation differentially associates with the methylation status of articular cartilage by acting on key mechanisms involved in osteoarthritis, such as apoptosis and metabolic and developmental processes.

Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

Polycystic ovary syndrome (PCOS) is a disorder characterized by hyperandrogenism, ovulatory dysfunction and polycystic ovarian morphology. Affected women frequently have metabolic disturbances including insulin resistance and dysregulation of glucose homeostasis. PCOS is diagnosed with two different sets of diagnostic criteria, resulting in a phenotypic spectrum of PCOS cases. The genetic similarities between cases diagnosed based on the two criteria have been largely unknown. Previous studies in Chinese and European subjects have identified 16 loci associated with risk of PCOS. We report a fixed-effect, inverse-weighted-variance meta-analysis from 10,074 PCOS cases and 103,164 controls of European ancestry and characterisation of PCOS related traits. We identified 3 novel loci (near PLGRKT, ZBTB16 and MAPRE1), and provide replication of 11 previously reported loci. Only one locus differed significantly in its association by diagnostic criteria; otherwise the genetic architecture was similar between PCOS diagnosed by self-report and PCOS diagnosed by NIH or non-NIH Rotterdam criteria across common variants at 13 loci. Identified variants were associated with hyperandrogenism, gonadotropin regulation and testosterone levels in affected women. Linkage disequilibrium score regression analysis revealed genetic correlations with obesity, fasting insulin, type 2 diabetes, lipid levels and coronary artery disease, indicating shared genetic architecture between metabolic traits and PCOS. Mendelian randomization analyses suggested variants associated with body mass index, fasting insulin, menopause timing, depression and male-pattern balding play a causal role in PCOS. The data thus demonstrate 3 novel loci associated with PCOS and similar genetic architecture for all diagnostic criteria. The data also provide the first genetic evidence for a male phenotype for PCOS and a causal link to depression, a previously hypothesized comorbid disease. Thus, the genetics provide a comprehensive view of PCOS that encompasses multiple diagnostic criteria, gender, reproductive potential and mental health.

Author Correction: Regulatory variants at KLF14 influence type 2 diabetes risk via a female-specific effect on adipocyte size and body composition.

In the version of this article originally published, minus signs were missing from the three ß-values for BMI given in Table 1. The errors have been corrected in the HTML and PDF versions of the article.