Trajectories of clinical characteristics, complications and treatment choices in data-driven subgroups of type 2 diabetes.

AIMS/HYPOTHESIS: This study aimed to explore the added value of subgroups that categorise individuals with type 2 diabetes by k-means clustering for two primary care registries (the Netherlands and Scotland), inspired by Ahlqvist's novel diabetes subgroups and previously analysed by Slieker et al. METHODS: We used two Dutch and Scottish diabetes cohorts (N=3054 and 6145; median follow-up=11.2 and 12.3 years, respectively) and defined five subgroups by k-means clustering with age at baseline, BMI, HbA1c, HDL-cholesterol and C-peptide. We investigated differences between subgroups by trajectories of risk factor values (random intercept models), time to diabetes-related complications (logrank tests and Cox models) and medication patterns (multinomial logistic models). We also compared directly using the clustering indicators as predictors of progression vs the k-means discrete subgroups. Cluster consistency over follow-up was assessed. RESULTS: Subgroups' risk factors were significantly different, and these differences remained generally consistent over follow-up. Among all subgroups, individuals with severe insulin resistance faced a significantly higher risk of myocardial infarction both before (HR 1.65; 95% CI 1.40, 1.94) and after adjusting for age effect (HR 1.72; 95% CI 1.46, 2.02) compared with mild diabetes with high HDL-cholesterol. Individuals with severe insulin-deficient diabetes were most intensively treated, with more than 25% prescribed insulin at 10 years of diagnosis. For severe insulin-deficient diabetes relative to mild diabetes, the relative risks for using insulin relative to no common treatment would be expected to increase by a factor of 3.07 (95% CI 2.73, 3.44), holding other factors constant. Clustering indicators were better predictors of progression variation relative to subgroups, but prediction accuracy may improve after combining both. Clusters were consistent over 8 years with an accuracy ranging from 59% to 72%. CONCLUSIONS/INTERPRETATION: Data-driven subgroup allocations were generally consistent over follow-up and captured significant differences in risk factor trajectories, medication patterns and complication risks. Subgroups serve better as a complement rather than as a basis for compressing clustering indicators.

An omics-based machine learning approach to predict diabetes progression: a RHAPSODY study.

AIMS/HYPOTHESIS: People with type 2 diabetes are heterogeneous in their disease trajectory, with some progressing more quickly to insulin initiation than others. Although classical biomarkers such as age, HbA1c and diabetes duration are associated with glycaemic progression, it is unclear how well such variables predict insulin initiation or requirement and whether newly identified markers have added predictive value. METHODS: In two prospective cohort studies as part of IMI-RHAPSODY, we investigated whether clinical variables and three types of molecular markers (metabolites, lipids, proteins) can predict time to insulin requirement using different machine learning approaches (lasso, ridge, GRridge, random forest). Clinical variables included age, sex, HbA1c, HDL-cholesterol and C-peptide. Models were run with unpenalised clinical variables (i.e. always included in the model without weights) or penalised clinical variables, or without clinical variables. Model development was performed in one cohort and the model was applied in a second cohort. Model performance was evaluated using Harrel's C statistic. RESULTS: Of the 585 individuals from the Hoorn Diabetes Care System (DCS) cohort, 69 required insulin during follow-up (1.0-11.4 years); of the 571 individuals in the Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS) cohort, 175 required insulin during follow-up (0.3-11.8 years). Overall, the clinical variables and proteins were selected in the different models most often, followed by the metabolites. The most frequently selected clinical variables were HbA1c (18 of the 36 models, 50%), age (15 models, 41.2%) and C-peptide (15 models, 41.2%). Base models (age, sex, BMI, HbA1c) including only clinical variables performed moderately in both the DCS discovery cohort (C statistic 0.71 [95% CI 0.64, 0.79]) and the GoDARTS replication cohort (C 0.71 [95% CI 0.69, 0.75]). A more extensive model including HDL-cholesterol and C-peptide performed better in both cohorts (DCS, C 0.74 [95% CI 0.67, 0.81]; GoDARTS, C 0.73 [95% CI 0.69, 0.77]). Two proteins, lactadherin and proto-oncogene tyrosine-protein kinase receptor, were most consistently selected and slightly improved model performance. CONCLUSIONS/INTERPRETATION: Using machine learning approaches, we show that insulin requirement risk can be modestly well predicted by predominantly clinical variables. Inclusion of molecular markers improves the prognostic performance beyond that of clinical variables by up to 5%. Such prognostic models could be useful for identifying people with diabetes at high risk of progressing quickly to treatment intensification. DATA AVAILABILITY: Summary statistics of lipidomic, proteomic and metabolomic data are available from a Shiny dashboard at https://rhapdata-app.vital-it.ch .

Branched-chain amino acid levels are inversely associated with incident and prevalent chronic kidney disease in people with type 2 diabetes.

AIM: To investigate the association of plasma metabolites with incident and prevalent chronic kidney disease (CKD) in people with type 2 diabetes and establish whether this association is causal. MATERIALS AND METHODS: The Hoorn Diabetes Care System cohort is a large prospective cohort consisting of individuals with type 2 diabetes from the northwest part of the Netherlands. In this cohort we assessed the association of baseline plasma levels of 172 metabolites with incident (Ntotal = 462/Ncase = 81) and prevalent (Ntotal = 1247/Ncase = 120) CKD using logistic regression. Additionally, replication in the UK Biobank, body mass index (BMI) mediation and causality of the association with Mendelian randomization was performed. RESULTS: Elevated levels of total and individual branched-chain amino acids (BCAAs)-valine, leucine and isoleucine-were associated with an increased risk of incident CKD, but with reduced odds of prevalent CKD, where BMI was identified as an effect modifier. The observed inverse effects were replicated in the UK Biobank. Mendelian randomization analysis did not provide evidence for a causal relationship between BCAAs and prevalent CKD. CONCLUSIONS: Our study shows the intricate relationship between plasma BCAA levels and CKD in individuals with type 2 diabetes. While an association exists, its manifestation varies based on disease status and BMI, with no definitive evidence supporting a causal link between BCAAs and prevalent CKD.

Apolipoprotein-CIII O-Glycosylation Is Associated with Micro- and Macrovascular Complications of Type 2 Diabetes.

Apolipoprotein-CIII (apo-CIII) inhibits the clearance of triglycerides from circulation and is associated with an increased risk of diabetes complications. It exists in four main proteoforms: O-glycosylated variants containing either zero, one, or two sialic acids and a non-glycosylated variant. O-glycosylation may affect the metabolic functions of apo-CIII. We investigated the associations of apo-CIII glycosylation in blood plasma, measured by mass spectrometry of the intact protein, and genetic variants with micro- and macrovascular complications (retinopathy, nephropathy, neuropathy, cardiovascular disease) of type 2 diabetes in a DiaGene study (n = 1571) and the Hoorn DCS cohort (n = 5409). Mono-sialylated apolipoprotein-CIII (apo-CIII1) was associated with a reduced risk of retinopathy (ß = -7.215, 95% CI -11.137 to -3.294) whereas disialylated apolipoprotein-CIII (apo-CIII2) was associated with an increased risk (ß = 5.309, 95% CI 2.279 to 8.339). A variant of the GALNT2-gene (rs4846913), previously linked to lower apo-CIII0a, was associated with a decreased prevalence of retinopathy (OR = 0.739, 95% CI 0.575 to 0.951). Higher apo-CIII1 levels were associated with neuropathy (ß = 7.706, 95% CI 2.317 to 13.095) and lower apo-CIII0a with macrovascular complications (ß = -9.195, 95% CI -15.847 to -2.543). In conclusion, apo-CIII glycosylation was associated with the prevalence of micro- and macrovascular complications of diabetes. Moreover, a variant in the GALNT2-gene was associated with apo-CIII glycosylation and retinopathy, suggesting a causal effect. The findings facilitate a molecular understanding of the pathophysiology of diabetes complications and warrant consideration of apo-CIII glycosylation as a potential target in the prevention of diabetes complications.

Perinatal exposure to the immune-suppressant di-n-octyltin dichloride affects brain development in rats.

Disruption of the immune system during embryonic brain development by environmental chemicals was proposed as a possible cause of neurodevelopmental disorders. We previously found adverse effects of di-n-octyltin dichloride (DOTC) on maternal and developing immune systems of rats in an extended one-generation reproductive toxicity study according to the OECD 443 test guideline. We hypothesize that the DOTC-induced changes in the immune system can affect neurodevelopment. Therefore, we used in-vivo MRI and PET imaging and genomics, in addition to behavioral testing and neuropathology as proposed in OECD test guideline 443, to investigate the effect of DOTC on structural and functional brain development. Male rats were exposed to DOTC (0, 3, 10, or 30 mg/kg of diet) from 2 weeks prior to mating of the F0-generation until sacrifice of F1-animals. The brains of rats, exposed to DOTC showed a transiently enlarged volume of specific brain regions (MRI), altered specific gravity, and transient hyper-metabolism ([18F]FDG PET). The alterations in brain development concurred with hyper-responsiveness in auditory startle response and slight hyperactivity in young adult animals. Genomics identified altered transcription of key regulators involved in neurodevelopment and neural function (e.g. Nrgrn, Shank3, Igf1r, Cck, Apba2, Foxp2); and regulators involved in cell size, cell proliferation, and organ development, especially immune system development and functioning (e.g. LOC679869, Itga11, Arhgap5, Cd47, Dlg1, Gas6, Cml5, Mef2c). The results suggest the involvement of immunotoxicity in the impairment of the nervous system by DOTC and support the hypothesis of a close connection between the immune and nervous systems in brain development.

Altered blood gene expression in the obesity-related type 2 diabetes cluster may be causally involved in lipid metabolism: a Mendelian randomisation study.

AIMS/HYPOTHESIS: The aim of this study was to identify differentially expressed long non-coding RNAs (lncRNAs) and mRNAs in whole blood of people with type 2 diabetes across five different clusters: severe insulin-deficient diabetes (SIDD), severe insulin-resistant diabetes (SIRD), mild obesity-related diabetes (MOD), mild diabetes (MD) and mild diabetes with high HDL-cholesterol (MDH). This was to increase our understanding of different molecular mechanisms underlying the five putative clusters of type 2 diabetes. METHODS: Participants in the Hoorn Diabetes Care System (DCS) cohort were clustered based on age, BMI, HbA1c, C-peptide and HDL-cholesterol. Whole blood RNA-seq was used to identify differentially expressed lncRNAs and mRNAs in a cluster compared with all others. Differentially expressed genes were validated in the Innovative Medicines Initiative DIabetes REsearCh on patient straTification (IMI DIRECT) study. Expression quantitative trait loci (eQTLs) for differentially expressed RNAs were obtained from a publicly available dataset. To estimate the causal effects of RNAs on traits, a two-sample Mendelian randomisation analysis was performed using public genome-wide association study (GWAS) data. RESULTS: Eleven lncRNAs and 175 mRNAs were differentially expressed in the MOD cluster, the lncRNA AL354696.2 was upregulated in the SIDD cluster and GPR15 mRNA was downregulated in the MDH cluster. mRNAs and lncRNAs that were differentially expressed in the MOD cluster were correlated among each other. Six lncRNAs and 120 mRNAs validated in the IMI DIRECT study. Using two-sample Mendelian randomisation, we found 52 mRNAs to have a causal effect on anthropometric traits (n=23) and lipid metabolism traits (n=10). GPR146 showed a causal effect on plasma HDL-cholesterol levels (p = 2×10-15), without evidence for reverse causality. CONCLUSIONS/INTERPRETATION: Multiple lncRNAs and mRNAs were found to be differentially expressed among clusters and particularly in the MOD cluster. mRNAs in the MOD cluster showed a possible causal effect on anthropometric traits, lipid metabolism traits and blood cell fractions. Together, our results show that individuals in the MOD cluster show aberrant RNA expression of genes that have a suggested causal role on multiple diabetes-relevant traits.

IgG N-glycans are associated with prevalent and incident complications of type 2 diabetes.

AIMS/HYPOTHESIS: Inflammation is important in the development of type 2 diabetes complications. The N-glycosylation of IgG influences its role in inflammation. To date, the association of plasma IgG N-glycosylation with type 2 diabetes complications has not been extensively investigated. We hypothesised that N-glycosylation of IgG may be related to the development of complications of type 2 diabetes. METHODS: In three independent type 2 diabetes cohorts, plasma IgG N-glycosylation was measured using ultra performance liquid chromatography (DiaGene n = 1815, GenodiabMar n = 640) and mass spectrometry (Hoorn Diabetes Care Study n = 1266). We investigated the associations of IgG N-glycosylation (fucosylation, galactosylation, sialylation and bisection) with incident and prevalent nephropathy, retinopathy and macrovascular disease using Cox- and logistic regression, followed by meta-analyses. The models were adjusted for age and sex and additionally for clinical risk factors. RESULTS: IgG galactosylation was negatively associated with prevalent and incident nephropathy and macrovascular disease after adjustment for clinical risk factors. Sialylation was negatively associated with incident diabetic nephropathy after adjustment for clinical risk factors. For incident retinopathy, similar associations were found for galactosylation, adjusted for age and sex. CONCLUSIONS: We showed that IgG N-glycosylation, particularly galactosylation and to a lesser extent sialylation, is associated with a higher prevalence and future development of macro- and microvascular complications of diabetes. These findings indicate the predictive potential of IgG N-glycosylation in diabetes complications and should be analysed further in additional large cohorts to obtain the power to solidify these conclusions.

Diabetes risk loci-associated pathways are shared across metabolic tissues.

AIMS/HYPOTHESIS: Numerous genome-wide association studies have been performed to understand the influence of genetic variation on type 2 diabetes etiology. Many identified risk variants are located in non-coding and intergenic regions, which complicates understanding of how genes and their downstream pathways are influenced. An integrative data approach will help to understand the mechanism and consequences of identified risk variants. METHODS: In the current study we use our previously developed method CONQUER to overlap 403 type 2 diabetes risk variants with regulatory, expression and protein data to identify tissue-shared disease-relevant mechanisms. RESULTS: One SNP rs474513 was found to be an expression-, protein- and metabolite QTL. Rs474513 influenced LPA mRNA and protein levels in the pancreas and plasma, respectively. On the pathway level, in investigated tissues most SNPs linked to metabolism. However, in eleven of the twelve tissues investigated nine SNPs were linked to differential expression of the ribosome pathway. Furthermore, seven SNPs were linked to altered expression of genes linked to the immune system. Among them, rs601945 was found to influence multiple HLA genes, including HLA-DQA2, in all twelve tissues investigated. CONCLUSION: Our results show that in addition to the classical metabolism pathways, other pathways may be important to type 2 diabetes that show a potential overlap with type 1 diabetes.

Prognostic models for predicting the risk of foot ulcer or amputation in people with type 2 diabetes: a systematic review and external validation study.

AIMS/HYPOTHESIS: Approximately 25% of people with type 2 diabetes experience a foot ulcer and their risk of amputation is 10-20 times higher than that of people without type 2 diabetes. Prognostic models can aid in targeted monitoring but an overview of their performance is lacking. This study aimed to systematically review prognostic models for the risk of foot ulcer or amputation and quantify their predictive performance in an independent cohort. METHODS: A systematic review identified studies developing prognostic models for foot ulcer or amputation over minimal 1 year follow-up applicable to people with type 2 diabetes. After data extraction and risk of bias assessment (both in duplicate), selected models were externally validated in a prospective cohort with a 5 year follow-up in terms of discrimination (C statistics) and calibration (calibration plots). RESULTS: We identified 21 studies with 34 models predicting polyneuropathy, foot ulcer or amputation. Eleven models were validated in 7624 participants, of whom 485 developed an ulcer and 70 underwent amputation. The models for foot ulcer showed C statistics (95% CI) ranging from 0.54 (0.54, 0.54) to 0.81 (0.75, 0.86) and models for amputation showed C statistics (95% CI) ranging from 0.63 (0.55, 0.71) to 0.86 (0.78, 0.94). Most models underestimated the ulcer or amputation risk in the highest risk quintiles. Three models performed well to predict a combined endpoint of amputation and foot ulcer (C statistics >0.75). CONCLUSIONS/INTERPRETATION: Thirty-four prognostic models for the risk of foot ulcer or amputation were identified. Although the performance of the models varied considerably, three models performed well to predict foot ulcer or amputation and may be applicable to clinical practice.

Replication and cross-validation of type 2 diabetes subtypes based on clinical variables: an IMI-RHAPSODY study.

AIMS/HYPOTHESIS: Five clusters based on clinical characteristics have been suggested as diabetes subtypes: one autoimmune and four subtypes of type 2 diabetes. In the current study we replicate and cross-validate these type 2 diabetes clusters in three large cohorts using variables readily measured in the clinic. METHODS: In three independent cohorts, in total 15,940 individuals were clustered based on age, BMI, HbA1c, random or fasting C-peptide, and HDL-cholesterol. Clusters were cross-validated against the original clusters based on HOMA measures. In addition, between cohorts, clusters were cross-validated by re-assigning people based on each cohort's cluster centres. Finally, we compared the time to insulin requirement for each cluster. RESULTS: Five distinct type 2 diabetes clusters were identified and mapped back to the original four All New Diabetics in Scania (ANDIS) clusters. Using C-peptide and HDL-cholesterol instead of HOMA2-B and HOMA2-IR, three of the clusters mapped with high sensitivity (80.6-90.7%) to the previously identified severe insulin-deficient diabetes (SIDD), severe insulin-resistant diabetes (SIRD) and mild obesity-related diabetes (MOD) clusters. The previously described ANDIS mild age-related diabetes (MARD) cluster could be mapped to the two milder groups in our study: one characterised by high HDL-cholesterol (mild diabetes with high HDL-cholesterol [MDH] cluster), and the other not having any extreme characteristic (mild diabetes [MD]). When these two milder groups were combined, they mapped well to the previously labelled MARD cluster (sensitivity 79.1%). In the cross-validation between cohorts, particularly the SIDD and MDH clusters cross-validated well, with sensitivities ranging from 73.3% to 97.1%. SIRD and MD showed a lower sensitivity, ranging from 36.1% to 92.3%, where individuals shifted from SIRD to MD and vice versa. People belonging to the SIDD cluster showed the fastest progression towards insulin requirement, while the MDH cluster showed the slowest progression. CONCLUSIONS/INTERPRETATION: Clusters based on C-peptide instead of HOMA2 measures resemble those based on HOMA2 measures, especially for SIDD, SIRD and MOD. By adding HDL-cholesterol, the MARD cluster based upon HOMA2 measures resulted in the current clustering into two clusters, with one cluster having high HDL levels. Cross-validation between cohorts showed generally a good resemblance between cohorts. Together, our results show that the clustering based on clinical variables readily measured in the clinic (age, HbA1c, HDL-cholesterol, BMI and C-peptide) results in informative clusters that are representative of the original ANDIS clusters and stable across cohorts. Adding HDL-cholesterol to the clustering resulted in the identification of a cluster with very slow glycaemic deterioration.

NACHO: an R package for quality control of NanoString nCounter data.

SUMMARY: The NanoStringTM nCounter® is a platform for the targeted quantification of expression data in biofluids and tissues. While software by the manufacturer is available in addition to third parties packages, they do not provide a complete quality control (QC) pipeline. Here, we present NACHO ('NAnostring quality Control dasHbOard'), a comprehensive QC R-package. The package consists of three subsequent steps: summarize, visualize and normalize. The summarize function collects all the relevant data and stores it in a tidy format, the visualize function initiates a dashboard with plots of the relevant QC outcomes. It contains QC metrics that are measured by default by the manufacturer, but also calculates other insightful measures, including the scaling factors that are needed in the normalization step. In this normalization step, different normalization methods can be chosen to optimally preprocess data. Together, NACHO is a comprehensive method that optimizes insight and preprocessing of nCounter® data. AVAILABILITY AND IMPLEMENTATION: NACHO is available as an R-package on CRAN and the development version on GitHub https://github.com/mcanouil/NACHO. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Prediction models for development of retinopathy in people with type 2 diabetes: systematic review and external validation in a Dutch primary care setting.

AIMS/HYPOTHESIS: The aims of this study were to identify all published prognostic models predicting retinopathy risk applicable to people with type 2 diabetes, to assess their quality and accuracy, and to validate their predictive accuracy in a head-to-head comparison using an independent type 2 diabetes cohort. METHODS: A systematic search was performed in PubMed and Embase in December 2019. Studies that met the following criteria were included: (1) the model was applicable in type 2 diabetes; (2) the outcome was retinopathy; and (3) follow-up was more than 1 year. Screening, data extraction (using the checklist for critical appraisal and data extraction for systemic reviews of prediction modelling studies [CHARMS]) and risk of bias assessment (by prediction model risk of bias assessment tool [PROBAST]) were performed independently by two reviewers. Selected models were externally validated in the large Hoorn Diabetes Care System (DCS) cohort in the Netherlands. Retinopathy risk was calculated using baseline data and compared with retinopathy incidence over 5 years. Calibration after intercept adjustment and discrimination (Harrell's C statistic) were assessed. RESULTS: Twelve studies were included in the systematic review, reporting on 16 models. Outcomes ranged from referable retinopathy to blindness. Discrimination was reported in seven studies with C statistics ranging from 0.55 (95% CI 0.54, 0.56) to 0.84 (95% CI 0.78, 0.88). Five studies reported on calibration. Eight models could be compared head-to-head in the DCS cohort (N = 10,715). Most of the models underestimated retinopathy risk. Validating the models against different severities of retinopathy, C statistics ranged from 0.51 (95% CI 0.49, 0.53) to 0.89 (95% CI 0.88, 0.91). CONCLUSIONS/INTERPRETATION: Several prognostic models can accurately predict retinopathy risk in a population-based type 2 diabetes cohort. Most of the models include easy-to-measure predictors enhancing their applicability. Tailoring retinopathy screening frequency based on accurate risk predictions may increase the efficiency and cost-effectiveness of diabetic retinopathy care. REGISTRATION: PROSPERO registration ID CRD42018089122.

Serum Matrix Metalloproteinases and Left Atrial Remodeling-The Hoorn Study.

Extracellular matrix protein turnover may play an important role in left atrial (LA) remodelling. The aim is to investigate the associations between matrix metalloproteinase (MMPs), tissue inhibitor of metalloproteinase (TIMP-1) and LA volume index (LAVI) and if these associations are independent of TIMP-1 levels. Participants from The Hoorn Study, a population-based cohort study (n = 674), underwent echocardiography. Serum MMPs (i.e., MMP-1, MMP-2, MMP-3, MMP-9, and MMP-10) and TIMP-1 levels were measured with ELISA. Multiple linear regression analyses were used. MMP-1 levels were not associated with LAVI. Higher MMP-2 levels were associated with larger LAVI (regression coefficient per SD increase in MMP (95% CI); 0.03 (0.01; 0.05). Higher MMP-3 and MMP-9 levels were associated with smaller LAVI; -0.04 (-0.07; -0.01) and -0.04 (-0.06; -0.02) respectively. Only in women were higher MMP-10 levels associated with larger LAVI; 0.04 (0.00; 0.07, p-interaction 0.04). Additionally, only in women were higher TIMP-1 levels associated with smaller LAVI; -0.05 (-0.09; -0.01, p-interaction 0.03). The associations between MMPs and LAVI were independent of TIMP-1 levels. In conclusion, serum MMPs are associated with LAVI, independent of CVD risk factors and TIMP-1 levels. In addition, these associations differ according to sex and within MMP subgroups. This shows that the role of MMPs in LA remodelling is complex.

Visit-to-visit variability of glycemia and vascular complications: the Hoorn Diabetes Care System cohort.

BACKGROUND: Glycemic variation has been suggested to be a risk factor for diabetes-related complications. Previous studies did not address confounding of diabetes duration, number of visits and length of follow-up. Here, we characterize glycemic variability over time and whether its relation to diabetes-related complications and mortality is independent from diabetes- and follow-up duration. MATERIALS AND METHODS: Individuals with type 2 diabetes (n = 6770) from the Hoorn Diabetes Care System cohort were included in this study. The coefficient of variation (CV) was calculated over 5-year sliding intervals. People divided in quintiles based on their CV. Cox proportional hazard models were used to investigate the role of glycemic CV as risk factor in diabetes-related complications and mortality. RESULTS: The coefficient of variation of glucose (FG-CV) increased with time, in contrast to HbA1c (HbA1c-CV). People with a high FG-CV were those with an early age of diabetes onset (ΔQ5-Q1 = - 2.39 years), a higher BMI (ΔQ5-Q1 = + 0.92 kg/m2), an unfavorable lipid profile, i.e. lower levels of HDL-C (ΔQ5-Q1 = - 0.06 mmol/mol) and higher triglycerides (ΔQ5-Q1 =+ 1.20 mmol/mol). People with the highest FG-CV in the first 5-year interval showed an increased risk of insulin initiation, retinopathy, macrovascular complications and mortality independent of mean glycemia, classical risk factors and medication use. For HbA1c, the associations were weaker and less consistent. CONCLUSIONS: Individuals with a higher FG-CV have an unfavorable metabolic profile and have an increased risk of developing micro- and macrovascular complications and mortality. The association of HbA1c-CV with metabolic outcomes and complications was less consistent in comparison to FG-CV.

HbA1c is associated with altered expression in blood of cell cycle- and immune response-related genes.

AIMS/HYPOTHESIS: Individuals with type 2 diabetes are heterogeneous in their glycaemic control as tracked by blood HbA1c levels. Here, we investigated the extent to which gene expression levels in blood reflect current and future HbA1c levels. METHODS: HbA1c levels at baseline and 1 and 2 year follow-up were compared with gene expression levels in 391 individuals with type 2 diabetes from the Hoorn Diabetes Care System Cohort (15,564 genes, RNA sequencing). The functions of associated baseline genes were investigated further using pathway enrichment analysis. Using publicly available data, we investigated whether the genes identified are also associated with HbA1c in the target tissues, muscle and pancreas. RESULTS: At baseline, 220 genes (1.4%) were associated with baseline HbA1c. Identified genes were enriched for cell cycle and complement system activation pathways. The association of 15 genes extended to the target tissues, muscle (n = 113) and pancreatic islets (n = 115). At follow-up, expression of 25 genes (0.16%) associated with 1 year HbA1c and nine genes (0.06%) with 2 year HbA1c. Five genes overlapped across all time points, and 18 additional genes between baseline and 1 year follow-up. After adjustment for baseline HbA1c, the number of significant genes at 1 and 2 years markedly decreased, suggesting that gene expression levels in whole blood reflect the current glycaemic state and but not necessarily the future glycaemic state. CONCLUSIONS/INTERPRETATION: HbA1c levels in individuals with type 2 diabetes are associated with expression levels of genes that link to the cell cycle and complement system activation.

A SNP panel for identification of DNA and RNA specimens.

BACKGROUND: SNP panels that uniquely identify an individual are useful for genetic and forensic research. Previously recommended SNP panels are based on DNA profiles and mostly contain intragenic SNPs. With the increasing interest in RNA expression profiles, we aimed for establishing a SNP panel for both DNA and RNA-based genotyping. RESULTS: To determine a small set of SNPs with maximally discriminative power, genotype calls were obtained from DNA and blood-derived RNA sequencing data belonging to healthy, geographically dispersed, Dutch individuals. SNPs were selected based on different criteria like genotype call rate, minor allele frequency, Hardy-Weinberg equilibrium and linkage disequilibrium. A panel of 50 SNPs was sufficient to identify an individual uniquely: the probability of identity was 6.9 × 10- 20 when assuming no family relations and 1.2 × 10- 10 when accounting for the presence of full sibs. The ability of the SNP panel to uniquely identify individuals on DNA and RNA level was validated in an independent population dataset. The panel is applicable to individuals from European descent, with slightly lower power in non-Europeans. Whereas most of the genes containing the 50 SNPs are expressed in various tissues, our SNP panel needs optimization for other tissues than blood. CONCLUSIONS: This first DNA/RNA SNP panel will be useful to identify sample mix-ups in biomedical research and for assigning DNA and RNA stains in crime scenes to unique individuals.

DNA Methylation Landscapes of Human Fetal Development.

Remodelling the methylome is a hallmark of mammalian development and cell differentiation. However, current knowledge of DNA methylation dynamics in human tissue specification and organ development largely stems from the extrapolation of studies in vitro and animal models. Here, we report on the DNA methylation landscape using the 450k array of four human tissues (amnion, muscle, adrenal and pancreas) during the first and second trimester of gestation (9,18 and 22 weeks). We show that a tissue-specific signature, constituted by tissue-specific hypomethylated CpG sites, was already present at 9 weeks of gestation (W9). Furthermore, we report large-scale remodelling of DNA methylation from W9 to W22. Gain of DNA methylation preferentially occurred near genes involved in general developmental processes, whereas loss of DNA methylation mapped to genes with tissue-specific functions. Dynamic DNA methylation was associated with enhancers, but not promoters. Comparison of our data with external fetal adrenal, brain and liver revealed striking similarities in the trajectory of DNA methylation during fetal development. The analysis of gene expression data indicated that dynamic DNA methylation was associated with the progressive repression of developmental programs and the activation of genes involved in tissue-specific processes. The DNA methylation landscape of human fetal development provides insight into regulatory elements that guide tissue specification and lead to organ functionality.

Alternative Routes to Induce Naïve Pluripotency in Human Embryonic Stem Cells.

Human embryonic stem cells (hESCs) closely resemble mouse epiblast stem cells exhibiting primed pluripotency unlike mouse ESCs (mESCs), which acquire a naïve pluripotent state. Efforts have been made to trigger naïve pluripotency in hESCs for subsequent unbiased lineage-specific differentiation, a common conundrum faced by primed pluripotent hESCs due to heterogeneity in gene expression existing within and between hESC lines. This required either ectopic expression of naïve genes such as NANOG and KLF2 or inclusion of multiple pluripotency-associated factors. We report here a novel combination of small molecules and growth factors in culture medium (2i/LIF/basic fibroblast growth factor + Ascorbic Acid + Forskolin) facilitating rapid induction of transgene-free naïve pluripotency in hESCs, as well as in mESCs, which has not been shown earlier. The converted naïve hESCs survived long-term single-cell passaging, maintained a normal karyotype, upregulated naïve pluripotency genes, and exhibited dependence on signaling pathways similar to naïve mESCs. Moreover, they undergo global DNA demethylation and show a distinctive long noncoding RNA profile. We propose that in our medium, the FGF signaling pathway via PI3K/AKT/mTORC induced the conversion of primed hESCs toward naïve pluripotency. Collectively, we demonstrate an alternate route to capture naïve pluripotency in hESCs that is fast, reproducible, supports naïve mESC derivation, and allows efficient differentiation.

MethylAid: visual and interactive quality control of large Illumina 450k datasets.

UNLABELLED: The Illumina 450k array is a frequently used platform for large-scale genome-wide DNA methylation studies, i.e. epigenome-wide association studies. Currently, quality control of 450k data can be performed with Illumina's GenomeStudio and is part of a limited number 450k analysis pipelines. However, GenomeStudio cannot handle large-scale studies, and existing pipelines provide limited options for quality control and neither support interactive exploration by the user. To aid the detection of bad-quality samples in large-scale genome-wide DNA methylation studies as flexible and transparent as possible, we have developed MethylAid; a visual and interactive Web application using RStudio's shiny package. Bad-quality samples are detected using sample-dependent and sample-independent quality control probes present on the array and user-adjustable thresholds. In-depth exploration of bad-quality samples can be performed using several interactive diagnostic plots. Furthermore, plots can be annotated with user-provided metadata, for example, to identify outlying batches. Our new tool makes quality assessment of 450k array data interactive, flexible and efficient and is, therefore, expected to be useful for both data analysts and core facilities. AVAILABILITY AND IMPLEMENTATION: MethylAid is implemented as an R/Bioconductor package (www.bioconductor.org/packages/3.0/bioc/html/MethylAid.html). A demo application is available from shiny.bioexp.nl/MethylAid.

Reassessing human MHC-I genetic diversity in T cell studies.

The Major Histocompatibility Complex class I (MHC-I) system plays a vital role in immune responses by presenting antigens to T cells. Allele specific technologies, including recombinant MHC-I technologies, have been extensively used in T cell analyses for COVID-19 patients and are currently used in the development of immunotherapies for cancer. However, the immense diversity of MHC-I alleles presents challenges. The genetic diversity serves as the foundation of personalized medicine, yet it also poses a potential risk of exacerbating healthcare disparities based on MHC-I alleles. To assess potential biases, we analysed (pre)clinical publications focusing on COVID-19 studies and T cell receptor (TCR)-based clinical trials. Our findings reveal an underrepresentation of MHC-I alleles associated with Asian, Australian, and African descent. Ensuring diverse representation is vital for advancing personalized medicine and global healthcare equity, transcending genetic diversity. Addressing this disparity is essential to unlock the full potential of T cells for enhancing diagnosis and treatment across all individuals.