Search

Abstract Background Type 2 diabetes (T2D) has reached epidemic proportions globally, including in Africa. However, molecular studies to understand the pathophysiology of T2D remain scarce outside Europe and North America. The aims of this study are to use an untargeted metabolomics approach to identify: (a) metabolites that are differentially expressed between individuals with and without T2D and (b) a metabolic signature associated with T2D in a population of Sub-Saharan Africa (SSA). Methods A total of 580 adult Nigerians from the Africa America Diabetes Mellitus (AADM) study were studied. The discovery study included 310 individuals (210 without T2D, 100 with T2D). Metabolites in plasma were assessed by reverse phase, ultra-performance liquid chromatography and mass spectrometry (RP)/UPLC-MS/MS methods on the Metabolon Platform. Welch’s two-sample t-test was used to identify differentially expressed metabolites (DEMs), followed by the construction of a biomarker panel using a random forest (RF) algorithm. The biomarker panel was evaluated in a replication sample of 270 individuals (110 without T2D and 160 with T2D) from the same study. Results Untargeted metabolomic analyses revealed 280 DEMs between individuals with and without T2D. The DEMs predominantly belonged to the lipid (51%, 142/280), amino acid (21%, 59/280), xenobiotics (13%, 35/280), carbohydrate (4%, 10/280) and nucleotide (4%, 10/280) super pathways. At the sub-pathway level, glycolysis, free fatty acid, bile metabolism, and branched chain amino acid catabolism were altered in T2D individuals. A 10-metabolite biomarker panel including glucose, gluconate, mannose, mannonate, 1,5-anhydroglucitol, fructose, fructosyl-lysine, 1-carboxylethylleucine, metformin, and methyl-glucopyranoside predicted T2D with an area under the curve (AUC) of 0.924 (95% CI: 0.845–0.966) and a predicted accuracy of 89.3%. The panel was validated with a similar AUC (0.935, 95% CI 0.906–0.958) in the replication cohort. The 10 metabolites in the biomarker panel correlated significantly with several T2D-related glycemic indices, including Hba1C, insulin resistance (HOMA-IR), and diabetes duration. Conclusions We demonstrate that metabolomic dysregulation associated with T2D in Nigerians affects multiple processes, including glycolysis, free fatty acid and bile metabolism, and branched chain amino acid catabolism. Our study replicated previous findings in other populations and identified a metabolic signature that could be used as a biomarker panel of T2D risk and glycemic control thus enhancing our knowledge of molecular pathophysiologic changes in T2D. The metabolomics dataset generated in this study represents an invaluable addition to publicly available multi-omics data on understudied African ancestry populations.

Abstract Background In vitro and in vivo studies have shown that certain cytokines and hormones may play a role in the development and progression of type 2 diabetes (T2D). However, studies on their role in T2D in humans are scarce. We evaluated associations between 11 circulating cytokines and hormones with T2D among a population of sub-Saharan Africans and tested for causal relationships using Mendelian randomization (MR) analyses. Methods We used logistic regression analysis adjusted for age, sex, body mass index, and recruitment country to regress levels of 11 cytokines and hormones (adipsin, leptin, visfatin, PAI-1, GIP, GLP-1, ghrelin, resistin, IL-6, IL-10, IL-1RA) on T2D among Ghanaians, Nigerians, and Kenyans from the Africa America Diabetes Mellitus study including 2276 individuals with T2D and 2790 non-T2D individuals. Similar linear regression models were fitted with homeostatic modelling assessments of insulin sensitivity (HOMA-S) and β-cell function (HOMA-B) as dependent variables among non-T2D individuals (n = 2790). We used 35 genetic variants previously associated with at least one of these 11 cytokines and hormones among non-T2D individuals as instrumental variables in univariable and multivariable MR analyses. Statistical significance was set at 0.0045 (0.05/11 cytokines and hormones). Results Circulating GIP and IL-1RA levels were associated with T2D. Nine of the 11 cytokines and hormones (exceptions GLP-1 and IL-6) were associated with HOMA-S, HOMA-B, or both among non-T2D individuals. Two-stage least squares MR analysis provided evidence for a causal effect of GIP and IL-RA on HOMA-S and HOMA-B in multivariable analyses (GIP ~ HOMA-S β = − 0.67, P-value = 1.88 × 10−6 and HOMA-B β = 0.59, P-value = 1.88 × 10−5; IL-1RA ~ HOMA-S β = − 0.51, P-value = 8.49 × 10−5 and HOMA-B β = 0.48, P-value = 5.71 × 10−4). IL-RA was partly mediated via BMI (30-34%), but GIP was not. Inverse variance weighted MR analysis provided evidence for a causal effect of adipsin on T2D (multivariable OR = 1.83, P-value = 9.79 × 10−6), though these associations were not consistent in all sensitivity analyses. Conclusions The findings of this comprehensive MR analysis indicate that circulating GIP and IL-1RA levels are causal for reduced insulin sensitivity and increased β-cell function. GIP’s effect being independent of BMI suggests that circulating levels of GIP could be a promising early biomarker for T2D risk. Our MR analyses do not provide conclusive evidence for a causal role of other circulating cytokines in T2D among sub-Saharan Africans.

Abstract Background Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery. Results To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3–5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism. Conclusions Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Abstract Background A complex set of perturbations occur in cytokines and hormones in the etiopathogenesis of obesity and related cardiometabolic conditions such as type 2 diabetes (T2D). Evidence for the genetic regulation of these cytokines and hormones is limited, particularly in African-ancestry populations. In order to improve our understanding of the biology of cardiometabolic traits, we investigated the genetic architecture of a large panel of obesity- related cytokines and hormones among Africans with replication analyses in African Americans. Methods We performed genome-wide association studies (GWAS) in 4432 continental Africans, enrolled from Ghana, Kenya, and Nigeria as part of the Africa America Diabetes Mellitus (AADM) study, for 13 obesity-related cytokines and hormones, including adipsin, glucose-dependent insulinotropic peptide (GIP), glucagon-like peptide-1 (GLP-1), interleukin-1 receptor antagonist (IL1-RA), interleukin-6 (IL-6), interleukin-10 (IL-10), leptin, plasminogen activator inhibitor-1 (PAI-1), resistin, visfatin, insulin, glucagon, and ghrelin. Exact and local replication analyses were conducted in African Americans (n = 7990). The effects of sex, body mass index (BMI), and T2D on results were investigated through stratified analyses. Results GWAS identified 39 significant (P value

Abstract Background Models including an interaction term and performing a joint test of SNP and/or interaction effect are often used to discover Gene-Environment (GxE) interactions. When the environmental exposure is a binary variable, analyses from exposure-stratified models which consist of estimating genetic effect in unexposed and exposed individuals separately can be of interest. In large-scale consortia focusing on GxE interactions in which only the joint test has been performed, it may be challenging to get summary statistics from both exposure-stratified and marginal (i.e not accounting for interaction) models. Results In this work, we developed a simple framework to estimate summary statistics in each stratum of a binary exposure and in the marginal model using summary statistics from the “joint” model. We performed simulation studies to assess our estimators’ accuracy and examined potential sources of bias, such as correlation between genotype and exposure and differing phenotypic variances within exposure strata. Results from these simulations highlight the high theoretical accuracy of our estimators and yield insights into the impact of potential sources of bias. We then applied our methods to real data and demonstrate our estimators’ retained accuracy after filtering SNPs by sample size to mitigate potential bias. Conclusions These analyses demonstrated the accuracy of our method in estimating both stratified and marginal summary statistics from a joint model of gene-environment interaction. In addition to facilitating the interpretation of GxE screenings, this work could be used to guide further functional analyses. We provide a user-friendly Python script to apply this strategy to real datasets. The Python script and documentation are available at https://gitlab.pasteur.fr/statistical-genetics/j2s .