Biologically Informed Polygenic Scores for Brain Insulin Receptor Network Are Associated with Cardiometabolic Risk Markers and Diabetes in Women
Article information
Abstract
Background
To investigate associations between variations in the co-expression-based brain insulin receptor polygenic score and cardiometabolic risk factors and diabetes mellitus.
Methods
This cross-sectional study included 1,573 participants from the Helsinki Birth Cohort Study. Biologically informed expression-based polygenic risk scores for the insulin receptor gene network were calculated for the hippocampal (hePRS-IR) and the mesocorticolimbic (mePRS-IR) regions. Cardiometabolic markers included body composition, waist circumference, circulating lipids, insulin-like growth factor 1 (IGF-1), and insulin-like growth factor-binding protein 1 and 3 (IGFBP-1 and -3). Glucose and insulin levels were measured during a standardized 2-hour 75 g oral glucose tolerance test and impaired glucose regulation status was defined by the World Health Organization 2019 criteria. Analyzes were adjusted for population stratification, age, smoking, alcohol consumption, socioeconomic status, chronic diseases, birth weight, and leisure-time physical activity.
Results
Multinomial logistic regression indicated that one standard deviation increase in hePRS-IR was associated with increased risk of diabetes mellitus in all participants (adjusted relative risk ratio, 1.17; 95% confidence interval, 1.01 to 1.35). In women, higher hePRS-IR was associated with greater waist circumference and higher body fat percentage, levels of glucose, insulin, total cholesterol, low-density lipoprotein cholesterol, triglycerides, apolipoprotein B, insulin, and IGFBP-1 (all P≤0.02). The mePRS-IR was associated with decreased IGF-1 level in women (P=0.02). No associations were detected in men and studied outcomes.
Conclusion
hePRS-IR is associated with sex-specific differences in cardiometabolic risk factor profiles including impaired glucose regulation, abnormal metabolic markers, and unfavorable body composition in women.
Highlights
• ePRSs, based on co-expression, reflect tissue-specific biological functions.
• ePRSs associated with the insulin receptor gene network are linked to type 2 diabetes.
• Corresponding links were with adverse lipid profiles and body composition.
• These associations were seen in the hippocampal area, exclusively among older women.
• The brain insulin receptor gene network appears to influence cardiometabolic risk.
INTRODUCTION
The human brain was once thought to be insensitive to insulin. Cumulative evidence from past decades has since confirmed the existence of central insulin activity and its critical role in cognition and metabolism [1].
Central insulin activity has been proposed to be an essential regulator of peripheral glucose homeostasis [2]. In animal models, insulin receptors are expressed on the endothelial cells of the brain blood barrier, suggesting transportation of peripheral insulin into the brain, as well as in a broad range of brain areas, including the limbic regions, the hypothalamus and the hippocampus [3,4]. Through insulin receptor signaling, central insulin action has been found to regulate food intake [5] and to reduce hepatic glucose production [6]. In addition, brain insulin resistance compromises the dopaminergic system which in turn increases energy intake and peripheral glucose levels [5].
Central insulin activity further affects fat metabolism. Central insulin sensitivity is associated with more favorable fat metabolism [7], while hypothalamic insulin resistance participates in the development of obesity [8]. Furthermore, while studies investigating peripheral insulin sensitivity and central insulin activity have provided mixed findings, with some studies demonstrating that insulin resistance in obese subjects reduces transportation of insulin into the brain [9,10] and other studies have shown that peripheral hyperinsulinemia and impaired glucose tolerance are associated with higher brain glucose uptake [11,12], it is clear that peripheral and central insulin sensitivity are closely linked. Sequentially, adverse cardiometabolic health, diabetes mellitus (DM), and unfavorable fat metabolism are closely associated with impairment in brain insulin signaling [13,14].
Brain insulin action has largely been investigated in animal models. In humans, central insulin resistance has been examined by applying neuroimaging techniques and administrating intranasal insulin [15], as well as by approaching it from a genetic perspective, for example by employing polygenic risk scores (PRSs). However, conventional PRSs often fail to take into account that genes operate in networks and have tissue-specific biological functions. Thus, a recent study developed expression-based genetic scores for mesocorticolimbic and hippocampal insulin receptor-related gene networks in order to enhance traditional PRSs [16]. With this approach, it was shown that brain region specific, biologically informed PRSs for the insulin receptor gene network (ePRS-IRs) were more strongly associated with Alzheimer’s disease, addiction and childhood impulsivity than the traditional PRSs [16].
The ePRS-IRs are devised using gene co-expression data (in a given tissues), which allows identification of gene co-expression networks. Single nucleotide polymorphisms (SNPs) of these genes are then used for the calculation of the ePRS-IRs. The SNPs from the network genes are functionally annotated and subjected to linkage disequilibrium clumping for removal of highly correlated SNPs. Then a count function of the number of alleles at a given SNP, weighted by the effect size of the association between the individual SNP and gene expression data in that specific tissue, is performed using Genotype-Tissue Expression (GTeX) [17]. The sum of these values from the total number of SNPs provides the ePRS-IRs. The ePRS-IR aggregates information on the relationship between the gene of interest and other genes in the genome, the levels of tissue-specific gene expression, the genetic variation of the target sample (given by the genotyping data) and the tissue-specific effect size of the association between genotyping and gene expression (given by GTeX). Therefore, variations in the ePRS-IR represent individual variations in the expression of the tissue-specific gene co-expression network [18]. In the case of our study, variations in the ePRS-IR represent individual variations in the expression of the insulin receptor gene network in the mesocorticolimbic area or in the hippocampus.
Although widely employed in a range of fields, PRSs have not previously been applied to investigate central insulin receptor mediated pathways and metabolic factors. Overall, research on the genetic predisposition for central insulin action and peripheral metabolism has been scarce. While there is evidence that insulin receptors are expressed in the hippocampus and the mesocorticolimbic area, research on whether the variation in expression of the insulin receptors in these areas is associated with peripheral glucose and energy homeostasis has not been conducted.
In this study, we aim to employ these novel ePRS-IRs to assess the association between the central insulin receptor gene network in the hippocampus and the mesocorticolimbic area, and markers of cardiometabolic health and DM among older men and women.
METHODS
Participants
The Helsinki Birth Cohort Study (HBCS) includes 13,345 individuals who were born between 1934 and 1944 at the Helsinki University Central Hospital (HUCH) or the Helsinki City Maternity Hospital [19] and attended child welfare clinics in Helsinki. Those individuals who were living in Finland in 1971 received a unique personal identification number, as did all individuals of the Finnish population. Out of 8,760 individuals born at HUCH, 2,902 were randomly selected and invited to the clinical cohort. Between 2001 and 2004, a baseline clinical examination was conducted involving 2003 cohort members. Of those, 1,573 had sufficient data after excluding individuals with missing information on ePRS-IRs (n=383), socioeconomic status (SES, n=4), population stratification (n=1), leisure-time physical activity (LTPA, n=29), chronic diseases (n=2), alcohol consumption (n=6), and smoking (n=5). Data on circulating insulin-like growth factor 1 (IGF-1), insulin-like growth factor-binding protein 1 (IGFBP-1) and 3 (IGFBP-3) levels were available only a random subsample of individuals who participated in clinical measurements in 2003 or earlier (n=454) [20]. The study was approved by the Ethics Committee of Epidemiology and Public Health of the Hospital District of Helsinki and Uusimaa (344/E3/2000) and that of the National Public Health Institute, Helsinki and follows the guidelines of the Declaration of Helsinki. All participants gave a written informed consent before participating in the study.
Expression-based polygenic risk score for brain insulin receptor network
Genotyping and ePRS-IR calculation were performed as previously described [16]. According to standard protocols, DNA was extracted from blood samples and genotyping was performed with the modified Illumina 610 k chip by the Wellcome Trust Sanger Institute (Cambridge, UK). Genomic coverage was extended by imputation using the 1000 Genomes Phase I integrated variant set (v3/April 2012; NCBI build 37/hg19) as the reference sample and IMPUTE2 software (https://mathgen.stats.ox.ac.uk/impute/impute_v2.html). Before imputing, quality control filters were applied by setting SNP clustering probability for each genotype at >95%, call rate at >95% for individuals and markers (99% for markers with minor allele frequency [MAF] <5%), MAF at >1%, and the P value for the Hardy-Weinberg Equilibrium exact test P>1×10–6. In addition, heterozygosity and gender and relatedness checks were performed and any discrepancies removed. The total number of SNPs in the imputed data was 39282668.
For the ePRS calculation, lists of genes co-expressed with the insulin receptor in the mesocorticolimbic system or hippocampus were created. SNPs from these gene networks were mapped, and the list of SNPs was submitted to linkage disequilibrium clumping. In HBCS, the clumped list of SNPs was weighted with the betas from the GTeX, a resource database and tissue bank for studying the relationship between genetic variation and gene expression in human tissues, using data from each respective brain region. The selection of the SNPs within a given clumping window was based on the lowest P value. Thus, biologically informed mesocorticolimbic (mePRS-IR) and hippocampal (hePRS-IR) specific co-expression-based polygenic scores for the insulin receptor gene network were calculated. For the analyses both hePRS-IR and mePRS-IR were standardized and reported as z-scores. For further details on calculation examples, please see the reference article [16].
Prediabetes and diabetes mellitus
At the time of the clinical examination, fasting plasma glucose was measured in all individuals. A standard 2-hour 75 g oral glucose tolerance test (OGTT) was applied for measuring glucose and insulin levels, as well as for diagnosing DM according to the World Health Organization (WHO) 2019 criteria [21]. Individuals who met the WHO 2019 criteria for impaired fasting glucose or impaired glucose tolerance during the OGTT were considered to have prediabetes. In addition, information on diabetes medication was collected from the Finnish national medication database, and individuals who received DM medication at the age of 40 or older were considered to have type 2 diabetes mellitus (T2DM). The Finnish national prescription drug reimbursement register has previously been compared with the national hospital discharge register and have shown that around 90% of those diagnosed and receiving medication for DM after the age 40 years have T2DM. Impaired glucose regulation status was available from 1,570 participants.
Body composition and anthropometrics
Height was measured with a stadiometer (KaWe), and weight with medical scales (alpha 770, SECA, Hamburg, Germany). Body mass index was calculated as weight in kilograms divided with height in meters squared. Waist circumference in centimeters was measured twice with a soft tape from the midpoint between lowest rib and iliac crest and the mean of the two measurement was reported. Body composition including body fat percentage, fat mass, and lean body mass was assessed by bioimpedance (InBody 3.0, Biospace Co. Ltd., Seoul, Korea) [22].
Blood testing and analyses
In the OGTT, we measured plasma glucose and insulin at the time of fasting, 30 minutes and 2 hours. Plasma glucose concentrations were determined with a hexokinase method and plasma insulin was measured with 2-site immunometric assays. Fasting plasma samples were used to measure serum total cholesterol, high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), triglycerides, apolipoprotein A, and apolipoprotein B concentrations with standard enzymatic methods. Circulating IGF-1 was measured by an IGF-1 enzyme-linked immunosorbent assay (ELISA) kit (DSL-10-5600, Diagnostic Systems Laboratories Inc., Webster, TX, USA). Serum IGFBP-1 was measured by sandwich-type immunofluorometry, as reported previously [20], and serum IGFBP-3 was measured by an immunofluorometric assay, using monoclonal antibodies (1B6/5C11) against IGFBP-3 [23]. Insulin resistance and β-cell function were determined by the homeostasis model assessment (HOMA), calculated by the formulas: HOMA-IR=(fasting glucose×fasting insulin)/22.5 [24] and HOMA-β=20×fasting insulin/(fasting glucose–3.5) [25]. The insulinogenic index was calculated as the ratio of the change in insulin and glucose responses from 0 to 30 minutes [26,27]. The analyses above were performed in an accredited hospital laboratory (Huslab, Helsinki, Finland).
Covariates
Smoking was coded as never, former, and current. Alcohol consumption was coded as never or having quit, less than once a week, or weekly. Highest attained SES was obtained from Statistics Finland and coded as high official, low official, self-employed, and manual workers [28]. The subjects’ past 12-month LTPA was assessed by the validated LTPA questionnaire; the Kuopio Ischemic Heart Disease Risk Factor Study (KIHD) [29]. As previously suggested [30], LTPA was measured in metabolic equivalents of task (MET) [31], which were multiplied with time (hours) and frequency to calculate MET-hours. Through validated questionnaires about chronic diseases, subjects were asked about any health conditions. Conditions included DM, cardiovascular conditions (congestive heart failure, arrhythmias, claudication, angina pectoris, previous heart attack, and stroke), lung diseases (asthma, emphysema, and chronic bronchitis), musculoskeletal disorders (rheumatoid arthritis and osteoporosis), and presence of cancer. The presence of comorbidities was coded as none, one, or two or more. Information on birth weight was retrieved from child welfare clinics as previously described [32]. Birth weight was coded as <3,000, 3,000–3,499, and ≥3,500 g. Self-reported diabetes medication was categorized as 0 (no diabetes medication) and 1 (usage of diabetes medication).
Statistical analysis
The data are reported as means (standard deviation or 95% confidence intervals [CIs]), medians (interquartile range) or counts (percentage). Analysis of variance was applied for continuous variables and chi-square test for categorical variables when analyzing the baseline characteristics. Multinomial logistic regression was applied to investigate the association between diabetes status (normoglycemia, prediabetes, DM) and the PRS-IRs and reported as relative risk ratios (RRR). Linear regression analyses were employed to investigate the association between the ePRSIRs and continuous variables. The bootstrap method with 5,000 repetitions was used to calculate 95% CIs. All analyses were performed separately for men and women as the body composition between the sexes differ. The crude analyses were adjusted for age and population stratification [33,34]. Fully adjusted models were additionally adjusted for age, smoking, alcohol consumption, SES, presence of chronic diseases, birth weight, and LTPA. Sensitivity analyses were applied for investigating the interaction between the self-reported diabetes medication and ePRSs and also diabetes status and ePRSs on continuous outcomes. A P<0.05 was considered to be statistically significant. Statistical analyses were carried out using Stata/MP version 16.1 (Stata Corporation, College Station, TX, USA).
Data availability
The data analyzed during the current study are available from the corresponding author on reasonable request.
RESULTS
This study included 1,573 participants from the HBCS, of which 889 were women. The characteristics of the study population are presented in Table 1.
Body composition and anthropometrics
In women, higher hePRS-IR was associated with greater waist circumference and fat mass as well as higher body fat percentage (Table 2). No association was detected between the hePRS-IR and lean body mass. mePRS-IR was not associated with any of the body composition or anthropometric measure etiher in men or women.
Glucose metabolism
Higher hePRS-IR was associated with higher glucose concentrations at 2-hour in the OGTT in women but not men (Table 2). hePRS-IR was associated with higher glucose concentrations at 30 minutes in the OGTT in the crude model (P=0.02), but after adjusting for confounding factors this association weakened slightly (P=0.05). No association was found for mePRS-IR in women or men.
Insulin metabolism
In women, higher hePRS-IR was linked with higher fasting insulin concentrations as well as insulin concentrations at 30 minutes and 2 hours. None of these associations were evident in men. No associations were found for mePRS-IR and markers of insulin metabolism.
Diabetes mellitus
When analyzing all participants, a one standard deviation increase in hePRS-IR increased the risk for prediabetes by 12% (RRR, 1.12; 95% CI, 1.00 to 1.26) and for DM by 17% (RRR, 1.17; 95% CI, 1.01 to 1.35) (Fig. 1). This association seemed to be stronger in women (RRR, 1.23; 95% CI, 0.99 to 1.54), although this was not statistically significant. No associations were found in men, nor for either men or women for the mePRS-IR.
Lipid metabolism
Higher hePRS-IR was associated with higher levels of total cholesterol, LDL-C, triglycerides, and apolipoprotein B, and with lower levels of HDL-C. No significant associations were evident in the mePRS-IR (Table 3).
Subsample analyses
We investigated the association between the ePRS-IRs and IGF-1, IGFBP-1, and IGFBP-3 from a subsample of participants. In this subsample of women, higher hePRS-IR was associated with lower levels of IGFBP-1 (b=–17.1; 95% CI, –28.5 to –5.8) (Supplementary Table 1). However, in men, higher mePRS-IR was linked to lower levels of IGFBP-1 (b=–11.1; 95% CI, –21.8 to –0.5).
Sensitivity analyses
No significant diabetes medication by ePRSs interaction term was found for any of the investigated continuous variables either in crude or fully adjusted models (Supplementary Table 2). We also investigated the interaction between the diabetes status and ePRSs on continuous outcomes. We found no significant interaction between the diabetes status by ePRSs on any of the variables in women (Supplementary Table 3). In men, however, there was a significant interaction effect on IGF-1 (P for interaction=0.014 for crude and 0.01 for full model), apolipoprotein A (P for interaction=0.04 for crude and 0.044 for full model), and for triglycerides (P for interaction=0.048 for crude and 0.056 for full model) (Supplementary Table 4).
DISCUSSION
There are significant biological differences between men and women that probably contribute to sex-specific differences in cardiovascular disease risk factors as well as treatment and prognosis. It is essential to understand these sex differences in order to get a better understanding of the pathophysiology and to optimize management of cardiovascular disease in both genders. In the present study we observed an association between hePRS-IR and an increased risk for T2DM which seemed to be stronger in women. We also detected an association between hePRS-IR and unfavorable cardiometabolic health, such as poorer lipid profile and body composition in women but not in men. Diabetes medication did not affect the results, and diabetes status only displayed a marginal effect on the association between the ePRS and studied outcomes in men, not in women.
The ePRS-IRs reflect biological function of the insulin receptor gene network. The association between the hePRS-IR and peripheral glucose and insulin levels points to central insulin metabolism as regulator of systemic glucose homeostasis. Previous studies have suggested that while peripheral insulin acts mostly as a metabolic regulatory hormone, central insulin has a range of effects in the brain with paramount systemic effects [35]. In turn, impaired systemic glucose homeostasis is a hallmark of the metabolic syndrome [13] and T2DM [36]. We also detected an increased risk for T2DM in individuals with higher expression of the hePRS-IR, which further supports our hypothesis that variation in the function of central insulin receptors has peripheral outcomes.
The hePRS-IR was also associated with poorer lipid profile and unfavorable body composition. This is in line with previous research in which adiposity has been found to be accompanied by alterations in brain insulin action—more so in women [36]. Intact central insulin signaling has been proposed to reduce increase in body weight [14], while disruption of normal insulin action in the brain has been found to lead to obesity through excessive food intake [37]. Moreover, obesity is often accompanied by increased levels of cardiovascular risk factors including dyslipidemia that was associated with higher hePRS in our study.
Obesity along with disturbed glucose regulation is accompanied by aberrant levels of several of the measured biomarkers in our study. Apolipoprotein A deficiency is associated with atherosclerosis and DM [38], while increased levels of apolipopeotein B contributes to hyperlipidemia [39]. IGFBP-1, which in our study was decreased in women with higher hePRS-IR, has been proposed to be a marker for insulin sensitivity [40], in addition to being inversely associated with hyperinsulinemia in postmenopausal women [41]. These biomarkers were associated with the hePRS-IR, which could reflect that the scores were intended to indeed capture biological function rather than specific disorders.
In our study, the associations between the brain insulin receptor gene network and disturbances in glucose regulation and body composition were seen only in women. Insulin-related metabolic conditions have previously been found to be overrepresented in women, especially after menopause due to the progressive loss of estrogen and its interactions with insulin [42]. Explanations for this include that the pre-menopausal protecting effect of estrogen on insulin resistance and T2DM disappears [43], resulting in a deteriorated metabolism in the brain [44]. Indeed, the female participants in our study population are postmenopausal. In addition, female over-representation of impaired glucose regulation has been attributed to sex differences in body composition [45]. In our study populations, hePRS-IR in women was associated with increased fat mass and waist circumference, which could contribute to the observed.
Consistently throughout our study, only the hePRS-IR was found to display any significant associations. A possible explanation for the associations being seen only with the hePRS-IR and not with the mePRS-IR might be due to the higher level of gene expression for IR in the hippocampus than the mesocorticolimbic system. Moreover, the co-expression gene network in the hippocampus is much larger, and therefore the signal captured by the hePRS-IR more aligned with insulin function in a broader way, as suggested by the results.
There are several strengths as well as limitations to our study. The HBCS is a well-characterized cohort which include the use of reliable measurements from the clinical examinationand. We applied both clinical and register-based data in the analyses. We were also able to measure a variety of different markers of cardiometabolic health. Information on diabetes status was collected based on OGTT performed during the study visit. Regarding the ePRSs, they have the potential to enhance disease diagnosing and treatment when applied in combination with other clinical risk factors and disease manifestations. The limitations to our study include the fact that the participants come from a homogenous, restricted area of Finland, which may affect generalizability and applicability of our results. Body composition was assessed with bioimpedance which is not the golden standard; however, bioimpedance does accurately measure a wide range of adiposity.
To conclude, variations in hippocampal insulin receptor gene network expression were associated with impaired glucose regulation and unfavorable cardiometabolic health in women. Thus, our research extends on previous work and offers further insight into the underlying pathological features associated with DM, obesity and cardiometabolic health, and by doing so, contributes to the understanding of these diseases and how to prevent them.
SUPPLEMENTARY MATERIALS
Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2023.0039.
Notes
CONFLICTS OF INTEREST
No potential conflict of interest relevant to this article was reported.
AUTHOR CONTRIBUTIONS
Conception or design: J.S.S., J.G.E., N.S.W.
Acquisition, analysis, interpretation of data: J.S.S., P.P.S., J.L., J.G.E., N.S.W.
Drafting the work or revising: all authors.
Final approval of the manuscript: J.S.S., P.P.S., J.L., H.K., R.K., M.S., J.G.E., N.S.W.
FUNDING
This study was supported by the Finnish Foundation for Cardiovascular Research, Finnish Foundation for Diabetes Research, Juho Vainio Foundation, Academy of Finland, Novo Nordisk Foundation, Signe and Ane Gyllenberg Foundation, Samfundet Folkhälsan, Finska Läkaresällskapet, Liv och Hälsa, European Commission FP7 (DORIAN) Grant Agreement No. 278603 and EU H2020-PHC-2014-DynaHealth Grant No. 633595 and EU Horizon 2020 Award 733206 LIFECYCLE. Patricia P. Silveira is supported by Canadian Institutes of Health Research (CIHR, PJT-166066, PI Patricia P. Silveira). Hannu Koistinen is supported by Sigrid Jusélius Foundation.
Acknowledgements
The authors would like to express their gratitude to the participants in the Helsinki Birth Cohort Study. We thank Annikki Löfhjelm for technical assistance.