Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 07 December 2023
Sec. Molecular Diagnostics and Therapeutics

Predicting type 2 diabetes risk before and after solid organ transplantation using polygenic scores in a Danish cohort

Quenia dos Santos&#x;Quenia dos Santos1Preston Leung&#x;Preston Leung1Christian W. Thorball,Christian W. Thorball2,3Bruno LedergerberBruno Ledergerber1Jacques Fellay,Jacques Fellay2,3Cameron R. MacPherson,Cameron R. MacPherson1,4Mads HornumMads Hornum5Cynthia Terrones-CamposCynthia Terrones-Campos1Allan RasmussenAllan Rasmussen6Finn Gustafsson,Finn Gustafsson7,8Michael Perch,Michael Perch7,8Sren S. SrensenSøren S. Sørensen5Christina EkenbergChristina Ekenberg1Jens D. LundgrenJens D. Lundgren1Bo Feldt-RasmussenBo Feldt-Rasmussen5Joanne Reekie
Joanne Reekie1*
  • 1Centre of Excellence for Health, Institut Roche, Immunity and Infections (CHIP), Rigshospitalet, University of Copenhagen, Copenhagen, Denmark
  • 2Precision Medicine Unit, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland
  • 3School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
  • 4Institut Roche, Boulogne-Billancourt, France
  • 5Department of Nephrology, Copenhagen University Hospital Rigshospitalet, Copenhagen, Denmark
  • 6Department of Surgical Gastroenterology, Rigshospitalet, Copenhagen, Denmark
  • 7Department of Cardiology, Rigshospitalet, Copenhagen, Denmark
  • 8Department of Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark

Type 2 diabetes mellitus (T2DM) can be multifactorial where both genetics and environmental factors play a role. We aimed to investigate the use of polygenic risk scores (PRS) in the prediction of pre-transplant T2DM and post-transplant diabetes mellitus (PTDM) among solid organ transplant (SOT) patients. Using non-genetic risk scores alone; and the combination with PRS, separate logistic regression models were built and compared using receiver operator curves. Patients were assessed pre-transplant and in three post-transplant periods: 0–45, 46–365 and >365 days. A higher PRS was significantly associated with increased odds of pre-transplant T2DM. However, no improvement was observed for pre-transplant T2DM prediction when comparing PRS combined with non-genetic risk scores to using non-genetic risk scores alone. This was also true for predictions of PTDM in all three post-transplant periods. This study demonstrated that polygenic risk was only associated with the risk of T2DM among SOT recipients prior to transplant and not for PTDM. Combining PRS with a clinical model of non-genetic risk scores did not significantly improve the predictive ability, indicating its limited clinical utility in identifying patients at high risk for T2DM before transplantation, suggesting that non-genetic or different genetic factors may contribute to PTDM.

Introduction

Post-transplant diabetes mellitus (PTDM) is a complication that can occur after a solid organ transplant (SOT). It refers to newly diagnosed diabetes mellitus following SOT, irrespective of diagnostic timing or whether T2DM was present but undetected before transplantation (Sharif et al., 2014). The prevalence of PTDM ranges from 10%–40% (Kasiske et al., 2003; Jenssen and Hartmann, 2019) and is associated with higher risk of death and complications after transplantation when compared to non-diabetes patients (Martínez-Dolz et al., 2005; Cho et al., 2012; Kim et al., 2017; Dos Santos et al., 2022).

Polygenic risk scores (PRS) are scoring profiles that extend information provided by genome wide association studies and focus on the collective contribution of individual genetic mutations to the phenotype of interest. While genome wide association studies provide p-values and effect sizes to assess the association between individual single nucleotide polymorphisms (SNPs) with the phenotype (Marees et al., 2018), these associations are treated as single events and those that have weak associations with the phenotype are often removed through multiple testing correction. PRS in contrast calculates a score to describe the effect of each SNP, thus it is possible to assess the risk for a given patient where multiple SNPs have been identified and can be treated collectively as a polygenic risk (Choi et al., 2020).

Previous studies assessing the use of PRS to predict T2DM have reported differing results. While some studies found PRS contributed very little to the prediction of T2DM in addition to a model with non-genetic risk scores for T2DM such as age, body mass index (BMI) and sex (Lango et al., 2008; Lyssenko et al., 2008; Chikowore et al., 2016), others demonstrated an improvement in the predictive performance of the model (Meigs et al., 2008; Chatterjee et al., 2013). As the knowledge in this area remains limited, this study aimed to assess the prediction of T2DM prior to transplant and PTDM among SOT recipients using non-genetic risk scores alone and the combination of non-genetic and PRS to investigate whether adding PRS could be useful for understanding the risk of T2DM and PTDM among SOT recipients.

Patients and methods

Study population

All patients ≥18 years old who underwent a SOT (heart, liver, lung, or kidney) at Rigshospitalet, Copenhagen University Hospital between January 2010 and December 2015 were eligible for inclusion. The study included SOT recipients who had an ethylenediaminetetraacetic acid stored blood sample in either the Region Hovedstadens Biobank (Rigshospitalet) or PERSIMUNE biobank (Rigshospitalet) and were part of the Management of post-Transplant infections in Collaborating Hospitals (MATCH) cohort (Lodding et al., 2015). For patients with more than one transplantation, only data related to the first transplant after 2010 were assessed.

Data sources

Clinical characteristics, sociodemographic and biochemical data were extracted from the MATCH database and the Centre of Excellence for Personalized Medicine for Infectious Complications in Immune Deficiency (PERSIMUNE) data warehouse, (https://www.persimune.dk/), which includes both regional and nationwide data collected prospectively as part of routine care.

Information on prescribed medications including insulin and oral anti-diabetic medication were obtained from the Electronic Prescription Medication (EPM) database that had records of hospital prescriptions from 2006 to 2016, and the Danish Prescription Database, a database with outpatient prescription records from 2004 onwards (Dos Santos et al., 2022). There was a gap in the data from EPM from May 2011 to December 2011 due to a change in systems. Data on specific immunosuppressive therapies for individual patients were not available. However, detailed information on the immunosuppressive schemes per transplant type have been previously published by Ekenberg C et al. (Ekenberg et al., 2020).

Data on diagnoses were retrieved from the National Patient Registry (Lynge et al., 2011) and Sundhedsdatabanken. The National Patient Registry was established in 1977 and contains national data on all hospital admissions up to 2016 while Sundhedsdatabanken holds data records for patients in the capital region of Denmark from 2008 until 2016 (Dos Santos et al., 2022). Mortality data from the Danish Civil Registration System were used for death dates (Pedersen, 2011).

Diabetes definition

Assessment of T2DM was consistent with a previously published study (Dos Santos et al., 2022). T2DM was assessed at pre-transplant and at three different PTDM time-periods.

1. 0–45 days post-transplant- “Early Likely PTDM” (EL-PTDM).

2. 46–365 days post-transplant

3. >365 days post-transplant

The fulfillment of at least one of the following criteria during the time-period of interest (for all time-periods, except before transplant), would classify the patients as “having developed diabetes mellitus”.

• A Hemoglobin A1C test ≥6.5 mmol/L or (Sharif et al., 2014).

• A prescription of antidiabetic medication from either EPM or Danish Prescription Database (use of insulin- Anatomical Therapeutic Chemical (ATC) code A10A, or use of oral antidiabetic medication -ATC code A10B) (Methodology WCCfDS, 2023).

• A diagnosis of diabetes (International Classification of Disease (ICD)-10 codes: E11, E13) (ICD10Data.com, 2023).

T2DM prior to transplant follows the above criteria with the exception of insulin treatment used during hospitalization (from EPM database) (Dos Santos et al., 2022), since these patients present a high incidence of corticoid-induced hyperglycaemia before transplant (Dos Santos et al., 2022). Additionally, prescription for antidiabetic medication during the first 15 days post-transplantation was not included in the definition of EL-PTDM due to a high prevalence of glucose intolerance and hyperglycaemia (Chakkera et al., 2009; Hecking et al., 2012). Patients classified with pre-transplant T2DM were classified as having T2DM in the entire follow-up period. Patients who were not classified as having pre-transplant T2DM were classified as PTDM if they met the T2DM definition in one of the time-periods of interest after-transplant, but could subsequently return to non-diabetes status in the following time-period if they did not meet the T2DM definition in the new time-period (Dos Santos et al., 2022).

SOT recipients with a diagnosis code for type-1 diabetes (E10) prior to transplantation were excluded as this study focuses on polygenic scores for T2DM.

Genotyping, quality checks and imputation

Genotyping was performed using Infinium Global Screening Array-24 v1.0 BeadChip (Illumina). SNP array data were provided in plink file formats with an initial count of 673,642 SNPs and were lifted to GRCh37 using CrossMap (Zhao et al., 2014). Initial quality check was performed using PLINK software (v2.00a3LM) (Chang et al., 2015), filtering out individuals and SNPs with less than 90% genotyping and Hardy-Weinberg equilibrium p-value less than 1 × 10−6. Higher or lower than expected genotype heterozygosity was also performed through PLINK. Individuals with F-values more than three standard deviations above or below the F-value mean were removed from the data set. Subsequent quality checks such as strand flipping, position, frequency and reference/alternate allele checks were performed using PLINK through a helper script (HRC-1000G-check-bim-v4.3.0) developed by McCarthy Group tools (Centre for Human Genetics, 2023) with default recommended parameters while specifying the input population to be European. Reference files for this process used 1,000 Genomes Phase 3 combined data set (https://mathgen.stats.ox.ac.uk/impute/1000GP_Phase3.html) (Genomes Project et al., 2015). Phasing of the data was performed using SHAPEIT (v4.2) (Delaneau et al., 2019) and the imputation of SNPs was performed using IMPUTE5 (Rubinacci et al., 2020) using the 1,000 Genomes SNP set as reference panel and generated a total of 48,864,655 SNPs. Imputed SNPs with INFO score less than 0.8 were removed. An additional set of quality checks using the same parameters were performed on the imputed data to ensure high quality imputed SNPs. SNPs with minor allele frequency less than 1% were removed from the data set followed by an additional check on the heterozygosity count using same F-value filtering criterion. The final number of SNPs used for downstream analysis was 8,771,317.

Genome-wide PRS

PRS selected for the MATCH cohort were obtained from a genome wide PRS study containing T2DM scores performed by Khera et al. (Khera et al., 2018) and downloaded from PGS Catalog (ID: PGS000014) (Lambert et al., 2021). This data set consisted of 6,917,436 SNPs generated from a European cohort. The scores for the MATCH cohort were calculated with PLINK (Chang et al., 2015) and a total of 3,134,520 SNPs were processed.

Statistical analyses

Patient characteristics at the time of transplantation were described and compared for those with and without pre-transplant T2DM. Continuous variables were analysed using the Wilcoxon test (nonparametric data) and χ2 test was used for categorical variables.

PRS was split into quintiles for better visualization of their potentially non-linear associations with T2DM.

To construct the non-genetic risk score clinical risk factors were included in a multivariable logistic regression model predicting T2DM prior to transplant. The individual event probabilities were then used to create a ‘non-genetic score’ for each patient. The clinical risk factors were selected a-priori based on those previously identified in the literature (Lango et al., 2008; Lyssenko et al., 2008; Meigs et al., 2008), and available in our database. These included age at transplant (in years), sex, transplant type, BMI and Charlson Comorbidity Index (CCI) (per point) (Quan et al., 2005). For the CCI (Quan et al., 2005), the two dimensions related to T2DM (presence of T2DM with and without chronic complications) were excluded from calculation of the index to avoid collinearity issues with our outcome (Dos Santos et al., 2022). All non-genetic risk factors were calculated immediately prior to transplant and were not recalculated during follow-up.

This single measure of “non-genetic risk score” was then split into quintiles similar to the PRS. These non-genetic risk quintiles were then used to check and visualize interactions with the quintiles of the PRS. T2DM and PTDM event variation explained by the different univariable logistic regression models (non-genetic risk scores alone or PRS alone) and multivariable models (non-genetic risk scores + PRS) were documented with receiver operating characteristic values and the area under the curve (AUC) for each model was compared using χ2 tests. The area under the precision-recall curve (AUPRC) for each model is also presented. The quintiles for both the non-genetic risk scores and PRS were recalculated during each time-period to account for changes in the population. Sensitivity analyses including only patients that had a diagnosis code for T2DM prior to transplant and excluding those with a medication prescription and/or a haemoglobin A1C ≥ 6.5 mmol/L were performed to assess model consistency. For all PTDM analyses, patients with T2DM prior to transplantation were excluded.

An additional sensitivity analysis was performed, using Poisson regression analysis to account for the varying follow-up time of the SOT recipients. In this analysis patients were included from day 46 post-transplant and followed until the first time they met our definition of PTDM. Those not developing PTDM were censored at their date of death, new transplant date, or the end of follow-up (31.12.2016), whichever occurred first. All data analyses were performed using SAS 9.4.

Results

A total of 959 SOT recipients had a transplant between January 2010 and December 2015. Due to the absence of available genetic information, 133 patients were excluded from this analysis. A further 121 patients with a diagnosis code for type-1 diabetes before transplant were also excluded. Thus, a total of 705 SOT recipients were included in the final data set.

Table 1 shows the patient characteristics prior to SOT. Prior to transplant, 521 patients were categorised as non-diabetes and 184 as having pre-transplant T2DM. Within the four SOT types, the percentage of patients with pre-transplant T2DM varied (p = 0.001). Heart transplant was the only transplant procedure to have more pre-transplant T2DM patients than non-diabetes patients (50.9% vs 49.1%). Additionally significant differences in the distributions were also identified between the groups when comparing BMI (p = 0.002) and age (p = 0.001). For BMI, 26.5% of the patients with BMI <25.0 had T2DM prior to transplant compared to 32.9% of the patients with BMI ≥25.0. In the age category, patients with pre-transplant T2DM were older (54.9 years, IQR: 46.7–62.1 years) than non-diabetic patients (49.4 years, IQR: 40.1–58.3 years).

TABLE 1
www.frontiersin.org

TABLE 1. Characteristics of non-diabetes and diabetes patients at baseline.

Pre-transplant diabetes

The results of the multivariable logistic regression model used to generate the ‘non-genetic risk score’ are given in Supplementary Table S1 while results of the univariate logistic regression model are provided in Supplementary Table S2. Figures 1A,B show the distribution of pre-transplant T2DM and non-diabetes patients across the quintiles of both the ‘non-genetic risk score’ and the PRS respectively. There was a significant association with pre-transplant T2DM and the quintiles for the non-genetic risk score (p < 0.0001, Figure 1A), where a higher proportion with pre-transplant T2DM were found in the highest quintiles. This was also confirmed in the univariable logistic regression modelling in Table 2, where patients in the third quintile (Odds Ratio (OR): 2.87, 95% confidence interval (CI): 1.49–5.54), fourth (OR: 4.33, 95% CI: 2.28–8.20) and fifth quintiles (OR: 6.59, 95% CI: 3.50–12.37) had significantly higher odds of having T2DM than patients in the first quintile.

FIGURE 1
www.frontiersin.org

FIGURE 1. Distribution of patients prior to transplant, split into quintiles (x-axis) and categorized into those that are classified as having diabetes (orange) and non-diabetes (blue). The y-axis illustrates percentage distribution with each set of quintiles being based on (A) non-genetic risk factors and on (B) PRS.

TABLE 2
www.frontiersin.org

TABLE 2. Univariable logistic regression models with the non-genetic risk scores split into quintiles for the development of pre-transplant diabetes, EL-PTDM and PTDM in each study period.

The same trend was identified for the PRS, when assessing the proportion with pre-transplant T2DM across the quintiles (Figure 1B) despite a non-significant p-value (p = 0.09). Patients in the highest quintile of the PRS category were found to have higher odds of pre-transplant T2DM (OR: 2.08, 95% CI: 1.21–3.57, p = 0.007) when compared to those in the first quintile using univariable logistic regression (Table 3).

TABLE 3
www.frontiersin.org

TABLE 3. Univariable logistic regression models with the PRS split into quintiles for the development of pre-transplant diabetes, EL-PTDM and PTDM in each study period.

Figure 2 illustrates the distribution of pre-transplant T2DM and non-diabetes patients in quintiles for the combined model of non-genetic risk factor and PRS, where a higher proportion of individuals with T2DM pre-transplant are observed in the highest quintile. Table 4 shows the AUC and the AUPRC for models including the non-genetic risk scores and models with both PRS and non-genetic risks scores for the prediction of pre-transplant T2DM. The addition of the PRS to the non-genetic risk scores model did not significantly improve the fit of the model (p = 0.10), although higher odds of pre-transplant T2DM were still observed among participants in the fifth quintile of PRS (adjusted OR: 1.94, 95% CI: 1.10–3.40, p = 0.02) when compared to the first quintile. Sensitivity analysis only defining T2DM using an ICD10 diagnosis code for T2DM prior transplant showed consistent results (Supplementary Table S3).

FIGURE 2
www.frontiersin.org

FIGURE 2. Distribution of patients prior to transplant, split into quintiles (x-axis) and categorized into those that are classified as having diabetes (orange) and non-diabetes (blue). The y-axis illustrates percentage distribution with each set of quintiles being based on the combined non-genetic risk factors and PRS.

TABLE 4
www.frontiersin.org

TABLE 4. Area under the curve (AUC) for multivariable models including the non-genetic risk scores*, and combined PRS and non-genetic risks factors for the prediction of pre-transplant diabetes, EL-PTDM and PTDM.

Post-transplant diabetes (PTDM)

At 0–45 days after transplant, 60 (11.5%) of the 521 non-diabetes patients had EL-PTDM.

For the period of 46–365 days post-transplant, a total of eight patients were excluded from the analysis (four died and four censored at the end of follow-up). None of the eight were classified as EL-PTDM. This resulted in a total of 513 patients being assessed in the 46–365 days post-transplant period. A total of 127 (24.8%) patients were defined as having PTDM in this period (48 that also had EL-PTDM, and 79 newly diagnosed).

In the last period (>365 days post-transplant), a further 26 patients were excluded from the analysis (25 who had died in the year following transplant and one who had less than a year’s follow-up). Thirteen of the 26 patients excluded had PTDM in the earlier period. Thus 487 patients were assessed for PDTM >365 days after transplant of which 106 (21.8%) met our definition of PTDM (65 who also had PTDM in the early period and 41 newly classified).

During the 46–365 and >365 days post-transplant period, the highest quintile of the non-genetic risk scores was associated with an increased odds of PTDM (OR: 1.92, 95% CI: 1.04–3.53, p = 0.03, and OR: 2.78, 95% CI: 1.44–5.34, p = 0.002, respectively) compared to the lowest quintile. An increased odds of EL-PTDM (OR: 1.86, 95% CI: 0.84–4.15, p = 0.12) was also observed in the highest quintile however this difference was not statistically significant (Table 2). No significant associations were observed between the odds of developing EL-PTDM or PTDM and the PRS, in any of the time-periods in univariable analysis (Table 3).

Comparison between the AUC of the non-genetic risk scores model and the AUC with the addition of the PRS is shown in Table 4. Overall, the AUC was lower for all the models in the post-transplant periods when compared to the models for pre-transplant T2DM. Additionally, no significant improvement in the AUC was observed in any of the time-periods for predicting EL-PTDM or PTDM when the PRS was added to the model.

Sensitivity analysis using Poisson regression

513 patients were included in the sensitivity analysis using Poisson regression analysis to account for the different follow-up times. Individuals in this analysis were included from day 46-post transplant and followed until the first time a patient met our definition for PTDM, at a minimum of 46 days after transplant. Patients in the fifth quintile of non-genetic risk score had higher rate of PTDM (incidence rate ratio (IRR) 1.68, 95% CI: 0.98–2.89, p = 0.05) during follow-up (from day 46 post-transplant), compared to patients in the first quintile. However, similar to the main analysis, there was no significant association identified between the PRS and PTDM when it was included in the model with the non-genetic risk score (global p-value = 0.23). However, patients in the fifth quintile of PRS were observed to have a two times higher rate of PTDM (adjusted IRR 2.00, 95% CI: 1.07–3.75, p = 0.02) when compared to patients in the first quintile (see Supplementary Table S4).

Discussion

The aim of this study was to assess whether a PRS for T2DM was associated with pre-transplant T2DM and whether the additional information from a PRS could improve the prediction of PTDM among SOT recipients compared to one with non-genetic scores alone. While we found that a high PRS was associated with increased odds of pre-transplant T2DM, including the PRS in a model with non-genetic risk scores, did not significantly improve prediction of T2DM compared to the model with non-genetic risk scores alone. Further, when considering EL-PTDM and PTDM in various time-periods, neither the non-genetic risk scores nor the PRS were strongly associated with the outcome.

For the pre-transplant period, patients in the highest quintile of the PRS and patients in the third, fourth and fifth quintiles of the non-genetic risk were found to have had higher odds of having T2DM. The addition of the PRS did not significantly improve the model compared to the non-genetic information alone. Sensitivity analysis including only patients with ICD-10 diagnosis codes for T2DM prior to transplant also had similar results. Our results are consistent with some of the previous studies assessing PRS to predict T2DM which also did not show a significant improvement of T2DM prediction when the PRS was added to a model with non-genetic risk scores (Lango et al., 2008; Lyssenko et al., 2008). A systematic review including studies that assessed PRS to predict T2DM (Padilla-Martínez et al., 2020) through comparing AUC of different models with and/or without the genetic information showed in most of the cases, the addition of the PRS to the non-genetic risk scores had a modest effect on the ability to predict T2DM (Lango et al., 2008; Lyssenko et al., 2008; Chikowore et al., 2016) at the best case scenario.

Some studies recorded better results when incorporating the genetic information with the clinical risk factors, but this improvement on T2DM prediction was small. While this could be due to the small number of SNPs when generating the PRS, it also suggests that clinical risk factors play a larger role than genetic factors alone. For example, in a study conducted by Meigs et al., (Meigs et al., 2008), using 18 SNP PRS, the AUC for T2DM prediction adjusted for age, sex and family history overall showed minor improvements. In our study, over six million SNPs from a PRS generated to predict T2DM (Lambert et al., 2021) were mapped to over three million SNPs of the MATCH cohort. Yet similarly, minor improvements were observed on only a certain portion of the data. Thus, our findings further extend the idea that PRS has a small effect when used in addition to non-genetic factors in the prediction of T2DM.

PTDM was considered at three different periods, EL-PTDM which was assessed from day 0 to day 45 post-transplant, then PTDM assessed at both 46–365 days post-transplant with the last period at >365 days post-transplant. No significant association was observed between PRS and a diagnosis of PTDM in any of the time-periods considered. Shaked et al. performed the only study which examined PTDM and PRS in liver and kidney recipients. Their study demonstrated that recipient T2DM PRS were independently associated with PTDM risk between 6 and 12 months after transplantation in both liver and kidney transplant recipients. This significantly improved PTDM prediction compared with a model that included only non-genetic risk scores for PTDM (Shaked et al., 2022). They also reported that T2DM PRS in liver donors, but not in kidney donors, was an independent risk factor for PTDM development. In our cohort, we did not have genetic information for the donors, thus it was not possible to generate a PRS for them and assess if the addition of this information would improve the performance of the model in terms of PTDM prediction. Furthermore, we had insufficient power to consider the transplant types individually.

Interestingly, the clinical model alone did not appear to be as good at predicting EL-PTDM or PTDM when compared to predicting pre-transplant T2DM. Individuals who met our definition of PTDM, particularly in the earlier periods post-transplant may contain a mixture of those who were already at an increased risk for T2DM and those where it was a consequence of the transplant regimen. As we only considered clinical factors measured prior to transplant there may be other treatment related clinical factors such as immunosuppressive regimen including steroid use that are more important for the development of PTDM following a SOT. Additionally, it is important to note this may also be a potential explanation for the lack of association between the PRS and PTDM, as individuals with a high genetic risk for T2DM may have already been diagnosed prior to their transplant and thus were excluded from this part of our analysis. After the removal of the group with T2DM prior to transplant, the interplay between genetic risk, and clinical risk factors following transplant maybe very different.

The limitations of this study should be highlighted. Firstly, this study did not include patient self-report of diabetes status but relied on a definition to identify T2DM patients used in a previous study (Dos Santos et al., 2022). However, several studies (Miller et al., 2004; Dos Santos et al., 2022) have adopted definitions to identify patients with T2DM based on similar criteria and they all performed quite well. Additionally, validation can also be derived from our finding that the highest quintile of the PRS score was significantly associated with pre-transplant T2DM. Secondly, the risk scores used in this study was generated based on a PRS for T2DM in a general European population. This is different to post-transplant patients where some other factors, such as stress from surgery and medication, can have an impact in the genetic information (Shaked et al., 2022). One additional limitation is the lack of information about the immunosuppressive medication. It is well known that immunosuppressive medication has great influence on glucose metabolism (Boloori et al., 2015), so this would be an important adjustment variable for our analyses if it had been available. Additionally, our results were based on a cohort of SOT recipients from one hospital in Denmark, validation in an independent cohort would add strength to the generalisability of our findings.

Strengths of this study are also worth mentioning. This was the first study that proposed the use of T2DM PRS to predict PTDM in a large and more varied cohort, including heart, liver, kidney, and lung recipients and in different time-periods after transplant (0–45 days after transplant, 46–365 days and >365 days after transplant). To our knowledge, there is only one study assessing the use of T2DM polygenic score to predict PTDM in liver and kidney recipients at 6–12 months after transplant (Shaked et al., 2022).

In summary, this study demonstrated that among SOT recipients, PRS could, help in the identification of patients at risk for T2DM prior to transplant. However, traditional non-genetic risk scores were just as good at predicting T2DM prior to transplant and the PRS did not provide any improvement in identification of high-risk individuals indicating its limited potential clinical utility Furthermore, the same predictive ability of PRS was not observed for PTDM, suggesting that non-genetic or different genetic factors, possibly related to the transplant itself, may contribute to the development of PTDM. Future studies assessing prediction of PTDM with the use of polygenic scores for T2DM are needed.

Data availability statement

The data analyzed in this study is subject to the following licenses/restrictions: Data requests must be reasonable and approved by PERSIMUNE Scientific Committee prior to data sharing. Requests to access these datasets should be directed to persimune.rigshospitalet@regionh.dk.

Ethics statement

The studies involving humans were approved by Danish National Committee on Health Research Ethics and Danish Data Protection Agency. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and institutional requirements. No potentially identifiable images or data are presented in this study.

Author contributions

QdS: Conceptualization, Data curation, Investigation, Methodology, Writing–original draft, Writing–review and editing, Formal Analysis, Software, Validation, Visualization. PL: Data curation, Formal Analysis, Methodology, Validation, Writing–original draft, Writing–review and editing, Visualization. CWT: Writing–review and editing, Methodology. BL: Methodology, Supervision, Writing–review and editing. JF: Writing–review and editing, Data curation. CRM: Writing–review and editing, Methodology, Supervision. MH: Writing–review and editing, Data curation. CT-C: Writing–review and editing, Data curation, Investigation. AR: Writing–review and editing, Data curation. FG: Writing–review and editing, Data curation. MP: Writing–review and editing, Data curation. SSS: Writing–review and editing, Data curation. CE: Writing–review and editing, Data curation, Investigation. JDL: Writing–review and editing, Data curation, Supervision. BFR: Writing–review and editing, Conceptualization, Data curation, Supervision. JR: Writing–original draft, Writing–review and editing, Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Supervision, Validation, Visualization.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by the Danish National Research Foundation (DNRF126). Funder(s) supporting this project had no participation of any form including study design, data collection and analysis, as well as the preparation of the manuscript and decision to publish.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2023.1282412/full#supplementary-material

Abbreviations

Area under the curve (AUC), area under the precision-recall curve (AUPRC), body mass index (BMI), Charlson Comorbidity Index (CCI), confidence interval (CI), Early Likely PTDM (EL-PTDM), Electronic Prescription Medication (EPM), Management of post-Transplant infections in Collaborating Hospitals (MATCH), odds ratio (OR), polygenic risk scores (PRS), post-transplant type 2 diabetes mellitus (PTDM), single nucleotide polymorphisms (SNPs), solid organ transplant (SOT), Type 2 diabetes mellitus (T2DM).

References

Boloori, A., Saghafian, S., Chakkera, H. A., and Cook, C. B. (2015). Characterization of remitting and relapsing hyperglycemia in post-renal-transplant recipients. PLoS One 10 (11), e0142363. doi:10.1371/journal.pone.0142363

PubMed Abstract | CrossRef Full Text | Google Scholar

Centre for Human Genetics (2023). Nuffield department of medicine. United Kingdom: Roosevelt Drive, Oxford OX3 7BN. Available from: https://www.well.ox.ac.uk/wrayner/tools/.

Google Scholar

Chakkera, H. A., Weil, E. J., Castro, J., Heilman, R. L., Reddy, K. S., Mazur, M. J., et al. (2009). Hyperglycemia during the immediate period after kidney transplantation. Clin. J. Am. Soc. Nephrol. 4 (4), 853–859. doi:10.2215/CJN.05471008

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, C. C., Chow, C. C., Tellier, L. C., Vattikuti, S., Purcell, S. M., and Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4, 7. doi:10.1186/s13742-015-0047-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chatterjee, N., Wheeler, B., Sampson, J., Hartge, P., Chanock, S. J., and Park, J. H. (2013). Projecting the performance of risk prediction based on polygenic analyses of genome-wide association studies. Nat. Genet. 45 (4), 400–405. doi:10.1038/ng.2579

PubMed Abstract | CrossRef Full Text | Google Scholar

Chikowore, T., van Zyl, T., Feskens, E. J., and Conradie, K. R. (2016). Predictive utility of a genetic risk score of common variants associated with type 2 diabetes in a black South African population. Diabetes Res. Clin. Pract. 122, 1–8. doi:10.1016/j.diabres.2016.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, M. S., Choi, H. I., Kim, I. O., Jung, S. H., Yun, T. J., Lee, J. W., et al. (2012). The clinical course and outcomes of post-transplantation diabetes mellitus after heart transplantation. J. Korean Med. Sci. 27 (12), 1460–1467. doi:10.3346/jkms.2012.27.12.1460

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, S. W., Mak, T. S., and O'Reilly, P. F. (2020). Tutorial: a guide to performing polygenic risk score analyses. Nat. Protoc. 15 (9), 2759–2772. doi:10.1038/s41596-020-0353-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Delaneau, O., Zagury, J. F., Robinson, M. R., Marchini, J. L., and Dermitzakis, E. T. (2019). Accurate, scalable and integrative haplotype estimation. Nat. Commun. 10 (1), 5436. doi:10.1038/s41467-019-13225-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Dos Santos, Q., Hornum, M., Terrones-Campos, C., Crone, C. G., Wareham, N. E., Soeborg, A., et al. (2022). Posttransplantation diabetes mellitus among solid organ recipients in a Danish cohort. Transpl. Int. 35, 10352. doi:10.3389/ti.2022.10352

PubMed Abstract | CrossRef Full Text | Google Scholar

Ekenberg, C., da Cunha-Bang, C., Lodding, I. P., Sørensen, S. S., Sengeløv, H., Perch, M., et al. (2020). Evaluation of an electronic, patient-focused management system aimed at preventing cytomegalovirus disease following solid organ transplantation. Transpl. Infect. Dis. 22 (2), e13252. doi:10.1111/tid.13252

PubMed Abstract | CrossRef Full Text | Google Scholar

Genomes Project, C., Auton, A., Brooks, L. D., Durbin, R. M., Garrison, E. P., Kang, H. M., et al. (2015). A global reference for human genetic variation. Nature 526 (7571), 68–74. doi:10.1038/nature15393

PubMed Abstract | CrossRef Full Text | Google Scholar

Hecking, M., Haidinger, M., Döller, D., Werzowa, J., Tura, A., Zhang, J., et al. (2012). Early basal insulin therapy decreases new-onset diabetes after renal transplantation. J. Am. Soc. Nephrol. 23 (4), 739–749. doi:10.1681/ASN.2011080835

PubMed Abstract | CrossRef Full Text | Google Scholar

ICD10Data.com (2023). Diabetes mellitus E08-E13.

Google Scholar

Jenssen, T., and Hartmann, A. (2019). Post-transplant diabetes mellitus in patients with solid organ transplants. Nat. Rev. Endocrinol. 15 (3), 172–188. doi:10.1038/s41574-018-0137-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasiske, B. L., Snyder, J. J., Gilbertson, D., and Matas, A. J. (2003). Diabetes mellitus after kidney transplantation in the United States. Am. J. Transpl. 3 (2), 178–185. doi:10.1034/j.1600-6143.2003.00010.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Khera, A. V., Chaffin, M., Aragam, K. G., Haas, M. E., Roselli, C., Choi, S. H., et al. (2018). Genome-wide polygenic scores for common diseases identify individuals with risk equivalent to monogenic mutations. Nat. Genet. 50 (9), 1219–1224. doi:10.1038/s41588-018-0183-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, H. J., Jung, S. H., Kim, J. J., Yun, T. J., Kim, J. B., Choo, S. J., et al. (2017). New-Onset diabetes mellitus after heart transplantation - incidence, risk factors and impact on clinical outcome. Circ. J. 81 (6), 806–814. doi:10.1253/circj.CJ-16-0963

PubMed Abstract | CrossRef Full Text | Google Scholar

Lambert, S. A., Gil, L., Jupp, S., Ritchie, S. C., Xu, Y., Buniello, A., et al. (2021). The Polygenic Score Catalog as an open database for reproducibility and systematic evaluation. Nat. Genet. 53 (4), 420–425. doi:10.1038/s41588-021-00783-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lango, H., Palmer, C. N., Morris, A. D., Zeggini, E., Hattersley, A. T., McCarthy, M. I., et al. (2008). Assessing the combined impact of 18 common genetic variants of modest effect sizes on type 2 diabetes risk. Diabetes 57 (11), 3129–3135. doi:10.2337/db08-0504

PubMed Abstract | CrossRef Full Text | Google Scholar

Lodding, I. P., Sengeløv, H., da Cunha-Bang, C., Iversen, M., Rasmussen, A., Gustafsson, F., et al. (2015). Clinical application of variation in replication kinetics during episodes of post-transplant cytomegalovirus infections. EBioMedicine 2 (7), 699–705. doi:10.1016/j.ebiom.2015.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Lynge, E., Sandegaard, J. L., and Rebolj, M. (2011). The Danish national patient register. Scand. J. Public Health 39 (7_Suppl. l), 30–33. doi:10.1177/1403494811401482

PubMed Abstract | CrossRef Full Text | Google Scholar

Lyssenko, V., Jonsson, A., Almgren, P., Pulizzi, N., Isomaa, B., Tuomi, T., et al. (2008). Clinical risk factors, DNA variants, and the development of type 2 diabetes. N. Engl. J. Med. 359 (21), 2220–2232. doi:10.1056/NEJMoa0801869

PubMed Abstract | CrossRef Full Text | Google Scholar

Marees, A. T., de Kluiver, H., Stringer, S., Vorspan, F., Curis, E., Marie-Claire, C., et al. (2018). A tutorial on conducting genome-wide association studies: quality control and statistical analysis. Int. J. Methods Psychiatr. Res. 27 (2), e1608. doi:10.1002/mpr.1608

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez-Dolz, L., Almenar, L., Martínez-Ortiz, L., Arnau, M. A., Chamorro, C., Moro, J., et al. (2005). Predictive factors for development of diabetes mellitus post-heart transplant. Transpl. Proc. 37 (9), 4064–4066. doi:10.1016/j.transproceed.2005.09.161

CrossRef Full Text | Google Scholar

Meigs, J. B., Shrader, P., Sullivan, L. M., McAteer, J. B., Fox, C. S., Dupuis, J., et al. (2008). Genotype score in addition to common risk factors for prediction of type 2 diabetes. N. Engl. J. Med. 359 (21), 2208–2219. doi:10.1056/NEJMoa0804742

PubMed Abstract | CrossRef Full Text | Google Scholar

Methodology WCCfDS (2023). Drugs used in diabetesOslo Norway: Norwegian Institute of Public Health.

Google Scholar

Miller, D. R., Safford, M. M., and Pogach, L. M. (2004). Who has diabetes? Best estimates of diabetes prevalence in the Department of Veterans Affairs based on computerized patient data. Diabetes Care 27 (Suppl. 2), B10–B21. doi:10.2337/diacare.27.suppl_2.b10

PubMed Abstract | CrossRef Full Text | Google Scholar

Padilla-Martínez, F., Collin, F., Kwasniewski, M., and Kretowski, A. (2020). Systematic review of polygenic risk scores for type 1 and type 2 diabetes. Int. J. Mol. Sci. 21 (5), 1703. doi:10.3390/ijms21051703

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedersen, C. B. (2011). The Danish Civil registration system. Scand. J. Public Health 39 (7 Suppl. l), 22–25. doi:10.1177/1403494810387965

PubMed Abstract | CrossRef Full Text | Google Scholar

Quan, H., Sundararajan, V., Halfon, P., Fong, A., Burnand, B., Luthi, J. C., et al. (2005). Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data. Med. Care 43 (11), 1130–1139. doi:10.1097/01.mlr.0000182534.19832.83

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubinacci, S., Delaneau, O., and Marchini, J. (2020). Genotype imputation using the positional burrows wheeler transform. PLoS Genet. 16 (11), e1009049. doi:10.1371/journal.pgen.1009049

PubMed Abstract | CrossRef Full Text | Google Scholar

Shaked, A., Loza, B. L., Van Loon, E., Olthoff, K. M., Guan, W., Jacobson, P. A., et al. (2022). Donor and recipient polygenic risk scores influence the risk of post-transplant diabetes. Nat. Med. 28 (5), 999–1005. doi:10.1038/s41591-022-01758-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharif, A., Hecking, M., de Vries, A. P., Porrini, E., Hornum, M., Rasoul-Rockenschaub, S., et al. (2014). Proceedings from an international consensus meeting on posttransplantation diabetes mellitus: recommendations and future directions. Am. J. Transpl. 14 (9), 1992–2000. doi:10.1111/ajt.12850

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H., Sun, Z., Wang, J., Huang, H., Kocher, J. P., and Wang, L. (2014). CrossMap: a versatile tool for coordinate conversion between genome assemblies. Bioinformatics 30 (7), 1006–1007. doi:10.1093/bioinformatics/btt730

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: type 2 diabetes mellitus, transplant, post-transplant diabetes mellitus, solid organ transplant recipient, polygenic risk score

Citation: dos Santos Q, Leung P, Thorball CW, Ledergerber B, Fellay J, MacPherson CR, Hornum M, Terrones-Campos C, Rasmussen A, Gustafsson F, Perch M, Sørensen SS, Ekenberg C, Lundgren JD, Feldt-Rasmussen B and Reekie J (2023) Predicting type 2 diabetes risk before and after solid organ transplantation using polygenic scores in a Danish cohort. Front. Mol. Biosci. 10:1282412. doi: 10.3389/fmolb.2023.1282412

Received: 24 August 2023; Accepted: 16 November 2023;
Published: 07 December 2023.

Edited by:

Hem Chandra Jha, Indian Institute of Technology Indore, India

Reviewed by:

Yize Li, Washington University in St. Louis, United States
Chachrit Khunsriraksakul, The Pennsylvania State University, United States

Copyright © 2023 dos Santos, Leung, Thorball, Ledergerber, Fellay, MacPherson, Hornum, Terrones-Campos, Rasmussen, Gustafsson, Perch, Sørensen, Ekenberg, Lundgren, Feldt-Rasmussen and Reekie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Joanne Reekie, joanne.reekie@regionh.dk

These authors share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.