Human inherited PD-L1 deficiency is clinically and immunologically less severe than PD-1 deficiency

We previously reported two siblings with inherited PD-1 deficiency who died from autoimmune pneumonitis at 3 and 11 years of age after developing other autoimmune manifestations, including type 1 diabetes (T1D). We report here two siblings, aged 10 and 11 years, with neonatal-onset T1D (diagnosed at the ages of 1 day and 7 wk), who are homozygous for a splice-site variant of CD274 (encoding PD-L1). This variant results in the exclusive expression of an alternative, loss-of-function PD-L1 protein isoform in overexpression experiments and in the patients' primary leukocytes. Surprisingly, cytometric immunophenotyping and single-cell RNA sequencing analysis on blood leukocytes showed largely normal development and transcriptional profiles across lymphoid and myeloid subsets in the PD-L1-deficient siblings, contrasting with the extensive dysregulation of both lymphoid and myeloid leukocyte compartments in PD-1 deficiency. Our findings suggest that PD-1 and PD-L1 are essential for preventing early-onset T1D but that, unlike PD-1 deficiency, PD-L1 deficiency does not lead to fatal autoimmunity with extensive leukocytic dysregulation.

Phenotype-based targeted treatment of SGLT2 inhibitors and GLP-1 receptor agonists in type 2 diabetes

AIMS/HYPOTHESIS

A precision medicine approach in type 2 diabetes could enhance targeting specific glucose-lowering therapies to individual patients most likely to benefit. We aimed to use the recently developed Bayesian causal forest (BCF) method to develop and validate an individualised treatment selection algorithm for two major type 2 diabetes drug classes, sodium-glucose cotransporter 2 inhibitors (SGLT2i) and glucagon-like peptide-1 receptor agonists (GLP1-RA).

METHODS

We designed a predictive algorithm using BCF to estimate individual-level conditional average treatment effects for 12-month glycaemic outcome (HbA1c) between SGLT2i and GLP1-RA, based on routine clinical features of 46,394 people with type 2 diabetes in primary care in England (Clinical Practice Research Datalink; 27,319 for model development, 19,075 for hold-out validation), with additional external validation in 2252 people with type 2 diabetes from Scotland (SCI-Diabetes [Tayside & Fife]). Differences in glycaemic outcome with GLP1-RA by sex seen in clinical data were replicated in clinical trial data (HARMONY programme: liraglutide [n=389] and albiglutide [n=1682]). As secondary outcomes, we evaluated the impacts of targeting therapy based on glycaemic response on weight change, tolerability and longer-term risk of new-onset microvascular complications, macrovascular complications and adverse kidney events.

RESULTS

Model development identified marked heterogeneity in glycaemic response, with 4787 (17.5%) of the development cohort having a predicted HbA1c benefit >3 mmol/mol (>0.3%) with SGLT2i over GLP1-RA and 5551 (20.3%) having a predicted HbA1c benefit >3 mmol/mol with GLP1-RA over SGLT2i. Calibration was good in hold-back validation, and external validation in an independent Scottish dataset identified clear differences in glycaemic outcomes between those predicted to benefit from each therapy. Sex, with women markedly more responsive to GLP1-RA, was identified as a major treatment effect modifier in both the UK observational datasets and in clinical trial data: HARMONY-7 liraglutide (GLP1-RA): 4.4 mmol/mol (95% credible interval [95% CrI] 2.2, 6.3) (0.4% [95% CrI 0.2, 0.6]) greater response in women than men. Targeting the two therapies based on predicted glycaemic response was also associated with improvements in short-term tolerability and long-term risk of new-onset microvascular complications.

CONCLUSIONS/INTERPRETATION

Precision medicine approaches can facilitate effective individualised treatment choice between SGLT2i and GLP1-RA therapies, and the use of routinely collected clinical features for treatment selection could support low-cost deployment in many countries.

Rare variant associations with birth weight identify genes involved in adipose tissue regulation, placental function and insulin-like growth factor signalling

Investigating the genetic factors influencing human birth weight may lead to biological insights into fetal growth and long-term health. Genome-wide association studies of birth weight have highlighted associated variants in more than 200 regions of the genome, but the causal genes are mostly unknown. Rare genetic variants with robust evidence of association are more likely to point to causal genes, but to date, only a few rare variants are known to influence birth weight. We aimed to identify genes that harbour rare variants that impact birth weight when carried by either the fetus or the mother, by analysing whole exome sequence data in UK Biobank participants. We annotated rare (minor allele frequency <0.1%) protein-truncating or high impact missense variants on whole exome sequence data in up to 234,675 participants with data on their own birth weight (fetal variants), and up to 181,883 mothers who reported the birth weight of their first child (maternal variants). Variants within each gene were collapsed to perform gene burden tests and for each associated gene, we compared the observed fetal and maternal effects. We identified 8 genes with evidence of rare fetal variant effects on birth weight, of which 2 also showed maternal effects. One additional gene showed evidence of maternal effects only. We observed 10/11 directionally concordant associations in an independent sample of up to 45,622 individuals (sign test P=0.01). Of the genes identified, IGF1R and PAPPA2 (fetal and maternal-acting) have known roles in insulin-like growth factor bioavailability and signalling. PPARG, INHBE and ACVR1C (all fetal-acting) have known roles in adipose tissue regulation and rare variants in the latter two also showed associations with favourable adiposity patterns in adults. We highlight the dual role of PPARG in both adipocyte differentiation and placental angiogenesis. NOS3, NRK, and ADAMTS8 (fetal and maternal-acting) have been implicated in both placental function and hypertension. Analysis of rare coding variants has identified regulators of fetal adipose tissue and fetoplacental angiogenesis as determinants of birth weight, as well as further evidence for the role of insulin-like growth factors.

Genetic drivers of heterogeneity in type 2 diabetes pathophysiology

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

Association of Gestational Free and Total Triiodothyronine With Gestational Hypertension, Preeclampsia, Preterm Birth, and Birth Weight: An Individual Participant Data Meta-analysis

CONTEXT

Triiodothyronine (T3) is the bioactive form of thyroid hormone. In contrast to thyroid-stimulating hormone and free thyroxine, we lack knowledge on the association of gestational T3 with adverse obstetric outcomes.

OBJECTIVE

To investigate the associaiton of gestational free or total T3 (FT3 or TT3) with adverse obstetric outcomes.

METHODS

We collected individual participant data from prospective cohort studies on gestational FT3 or TT3, adverse obstetric outcomes (preeclampsia, gestational hypertension, preterm birth and very preterm birth, small for gestational age [SGA], and large for gestational age [LGA]), and potential confounders. We used mixed-effects regression models adjusting for potential confounders.

RESULTS

The final study population comprised 33 118 mother-child pairs of which 27 331 had data on FT3 and 16 164 on TT3. There was a U-shaped association of FT3 with preeclampsia (P = .0069) and a J-shaped association with the risk of gestational hypertension (P = .029). Higher TT3 was associated with a higher risk of gestational hypertension (OR per SD of TT3 1.20, 95% CI 1.08 to 1.33; P = .0007). A lower TT3 but not FT3 was associated with a higher risk of very preterm birth (OR 0.72, 95% CI 0.55 to 0.94; P = .018). TT3 but not FT3 was positively associated with birth weight (mean difference per 1 SD increase in TT3 12.8, 95% CI 6.5 to 19.1 g, P < .0001) but there was no association with SGA or LGA.

CONCLUSION

This study provides new insights on the association of gestational FT3 and TT3 with major adverse pregnancy outcomes that form the basis for future studies required to elucidate the effects of thyroid function on pregnancy outcomes. Based on the current study, routine FT3 or TT3 measurements for the assessment of thyroid function during pregnancy do not seem to be of added value in the risk assessment for adverse outcomes.

Penetrance and expressivity of mitochondrial variants in a large clinically unselected population

Whole genome sequencing (WGS) from large clinically unselected cohorts provides a unique opportunity to assess the penetrance and expressivity of rare and/or known pathogenic mitochondrial variants in population. Using WGS from 179 862 clinically unselected individuals from the UK Biobank, we performed extensive single and rare variant aggregation association analyses of 15 881 mtDNA variants and 73 known pathogenic variants with 15 mitochondrial disease-relevant phenotypes. We identified 12 homoplasmic and one heteroplasmic variant (m.3243A>G) with genome-wide significant associations in our clinically unselected cohort. Heteroplasmic m.3243A>G (MAF = 0.0002, a known pathogenic variant) was associated with diabetes, deafness and heart failure and 12 homoplasmic variants increased aspartate aminotransferase levels including three low-frequency variants (MAF ~0.002 and beta~0.3 SD). Most pathogenic mitochondrial disease variants (n = 66/74) were rare in the population (<1:9000). Aggregated or single variant analysis of pathogenic variants showed low penetrance in unselected settings for the relevant phenotypes, except m.3243A>G. Multi-system disease risk and penetrance of diabetes, deafness and heart failure greatly increased with m.3243A>G level ≥ 10%. The odds ratio of these traits increased from 5.61, 12.3 and 10.1 to 25.1, 55.0 and 39.5, respectively. Diabetes risk with m.3243A>G was further influenced by type 2 diabetes genetic risk. Our study of mitochondrial variation in a large-unselected population identified novel associations and demonstrated that pathogenic mitochondrial variants have lower penetrance in clinically unselected settings. m.3243A>G was an exception at higher heteroplasmy showing a significant impact on health making it a good candidate for incidental reporting.

Hyperglycaemia is a causal risk factor for upper limb pathologies

BACKGROUND

Diabetes (regardless of type) and obesity are associated with a range of musculoskeletal disorders. The causal mechanisms driving these associations are unknown for many upper limb pathologies. We used genetic techniques to test the causal link between glycemia, obesity and musculoskeletal conditions.

METHODS

In the UK Biobank's unrelated European cohort (N = 379 708) we performed mendelian randomisation (MR) analyses to test for a causal effect of long-term high glycaemia and adiposity on four musculoskeletal pathologies: frozen shoulder, Dupuytren's disease, carpal tunnel syndrome and trigger finger. We also performed single-gene MR using rare variants in the GCK gene.

RESULTS

Using MR, we found evidence that long-term high glycaemia has a causal role in the aetiology of upper limb conditions. A 10-mmol/mol increase in genetically predicted haemoglobin A1C (HbA1c) was associated with frozen shoulder: odds ratio (OR) = 1.50 [95% confidence interval (CI), 1.20-1.88], Dupuytren's disease: OR = 1.17 (95% CI, 1.01-1.35), trigger finger: OR = 1.30 (95% CI, 1.09-1.55) and carpal tunnel syndrome: OR = 1.20 (95% CI, 1.09-1.33). Carriers of GCK mutations have increased odds of frozen shoulder: OR = 7.16 (95% CI, 2.93-17.51) and carpal tunnel syndrome: OR = 2.86 (95% CI, 1.50-5.44) but not Dupuytren's disease or trigger finger. We found evidence that an increase in genetically predicted body mass index (BMI) of 5 kg/m2 was associated with carpal tunnel syndrome: OR = 1.13 (95% CI, 1.10-1.16) and associated negatively with Dupuytren's disease: OR = 0.94 (95% CI, 0.90-0.98), but no evidence of association with frozen shoulder or trigger finger. Trigger finger (OR 1.96 (95% CI, 1.42-2.69) P = 3.6e-05) and carpal tunnel syndrome [OR 1.63 (95% CI, 1.36-1.95) P = 8.5e-08] are associated with genetically predicted unfavourable adiposity increase of one standard deviation of body fat.

CONCLUSIONS

Our study consistently demonstrates a causal role of long-term high glycaemia in the aetiology of upper limb musculoskeletal conditions. Clinicians treating diabetes patients should be aware of these complications in clinic, specifically those managing the care of GCK mutation carriers. Upper limb musculoskeletal conditions should be considered diabetes complications.

Integrating multiple lines of evidence to assess the effects of maternal BMI on pregnancy and perinatal outcomes

BACKGROUND

Higher maternal pre-pregnancy body mass index (BMI) is associated with adverse pregnancy and perinatal outcomes. However, whether these associations are causal remains unclear.

METHODS

We explored the relation of maternal pre-/early-pregnancy BMI with 20 pregnancy and perinatal outcomes by integrating evidence from three different approaches (i.e. multivariable regression, Mendelian randomisation, and paternal negative control analyses), including data from over 400,000 women.

RESULTS

All three analytical approaches supported associations of higher maternal BMI with lower odds of maternal anaemia, delivering a small-for-gestational-age baby and initiating breastfeeding, but higher odds of hypertensive disorders of pregnancy, gestational hypertension, preeclampsia, gestational diabetes, pre-labour membrane rupture, induction of labour, caesarean section, large-for-gestational age, high birthweight, low Apgar score at 1 min, and neonatal intensive care unit admission. For example, higher maternal BMI was associated with higher risk of gestational hypertension in multivariable regression (OR = 1.67; 95% CI = 1.63, 1.70 per standard unit in BMI) and Mendelian randomisation (OR = 1.59; 95% CI = 1.38, 1.83), which was not seen for paternal BMI (OR = 1.01; 95% CI = 0.98, 1.04). Findings did not support a relation between maternal BMI and perinatal depression. For other outcomes, evidence was inconclusive due to inconsistencies across the applied approaches or substantial imprecision in effect estimates from Mendelian randomisation.

CONCLUSIONS

Our findings support a causal role for maternal pre-/early-pregnancy BMI on 14 out of 20 adverse pregnancy and perinatal outcomes. Pre-conception interventions to support women maintaining a healthy BMI may reduce the burden of obstetric and neonatal complications.

FUNDING

Medical Research Council, British Heart Foundation, European Research Council, National Institutes of Health, National Institute for Health Research, Research Council of Norway, Wellcome Trust.

Incretin hormone responses to carbohydrate and protein/fat are preserved in adults with sulfonylurea-treated KCNJ11 neonatal diabetes

The incretin hormones glucagon-like peptide-1 (GLP-1) and glucose-dependent insulinotropic polypeptide (GIP), are thought to be the main drivers of insulin secretion in individuals with sulfonylurea (SU)-treated KCNJ11 permanent neonatal diabetes. The aim of this study was to assess for the first time the incretin hormone response to carbohydrate and protein/fat in adults with sulfonylurea-treated KCNJ11 permanent neonatal diabetes compared with that of controls without diabetes. Participants were given a breakfast high in carbohydrate and an isocaloric breakfast high in protein/fat on two different mornings. Incremental area under the curve and total area under the curve (0-240 minutes) for total GLP-1 and GIP were compared between groups, using non-parametric statistical methods. Post-meal GLP-1 and GIP secretion were similar in cases and controls, suggesting this process is adenosine triphosphate-sensitive potassium channel-independent. Future research will investigate whether treatments targeting the incretin pathway are effective in individuals with KCNJ11 permanent neonatal diabetes who do not have good glycemic control on sulfonylurea alone.

Primate-specific ZNF808 is essential for pancreatic development in humans

Identifying genes linked to extreme phenotypes in humans has the potential to highlight biological processes not shared with all other mammals. Here, we report the identification of homozygous loss-of-function variants in the primate-specific gene ZNF808 as a cause of pancreatic agenesis. ZNF808 is a member of the KRAB zinc finger protein family, a large and rapidly evolving group of epigenetic silencers which target transposable elements. We show that loss of ZNF808 in vitro results in aberrant activation of regulatory potential contained in the primate-specific transposable elements it represses during early pancreas development. This leads to inappropriate specification of cell fate with induction of genes associated with liver identity. Our results highlight the essential role of ZNF808 in pancreatic development in humans and the contribution of primate-specific regions of the human genome to congenital developmental disease.

Bringing precision medicine to the management of pregnancy in women with glucokinase-MODY: a study of diagnostic accuracy and feasibility of non-invasive prenatal testing

AIMS/HYPOTHESIS

In pregnancies where the mother has glucokinase-MODY (GCK-MODY), fetal growth is determined by fetal genotype. When the fetus inherits a maternal pathogenic GCK variant, normal fetal growth is anticipated, and insulin treatment of maternal hyperglycaemia is not recommended. At present, fetal genotype is estimated from measurement of fetal abdominal circumference on ultrasound. Non-invasive prenatal testing of fetal GCK genotype (NIPT-GCK) using cell-free DNA in maternal blood has recently been developed. We aimed to compare the diagnostic accuracy of NIPT-GCK with that of ultrasound, and determine the feasibility of using NIPT-GCK to guide pregnancy management.

METHODS

We studied an international cohort of pregnant women with hyperglycaemia due to GCK-MODY. We compared the diagnostic accuracy of NIPT-GCK with that of measurement of fetal abdominal circumference at 28 weeks' gestation (n=38) using a directly genotyped offspring sample as the reference standard. In a feasibility study, we assessed the time to result given to clinicians in 43 consecutive pregnancies affected by GCK-MODY between July 2019 and September 2021.

RESULTS

In terms of diagnostic accuracy, NIPT-GCK was more sensitive and specific than ultrasound in predicting fetal genotype (sensitivity 100% and specificity 96% for NIPT-GCK vs sensitivity 53% and specificity 61% for fetal abdominal circumference 75th percentile). In terms of feasibility, a valid NIPT-GCK fetal genotype (≥95% probability) was reported in all 38 pregnancies with an amenable variant and repeated samples when needed. The median time to report was 5 weeks (IQR 3-8 weeks). For the 25 samples received before 20 weeks' gestation, results were reported at a median gestational age of 20 weeks (IQR 18-24), with 23/25 (92%) reported before 28 weeks.

CONCLUSIONS/INTERPRETATION

Non-invasive prenatal testing of fetal genotype in GCK-MODY pregnancies is highly accurate and is capable of providing a result before the last trimester for most patients. This means that non-invasive prenatal testing of fetal genotype is the optimal approach to management of GCK-MODY pregnancies.

Genome-wide association study of placental weight identifies distinct and shared genetic influences between placental and fetal growth

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

The Role of ONECUT1 Variants in Monogenic and Type 2 Diabetes Mellitus

ONECUT1 (also known as HNF6) is a transcription factor involved in pancreatic development and β-cell function. Recently, biallelic variants in ONECUT1 were reported as a cause of neonatal diabetes mellitus (NDM) in two subjects, and missense monoallelic variants were associated with type 2 diabetes and possibly maturity-onset diabetes of the young (MODY). Here we examine the role of ONECUT1 variants in NDM, MODY, and type 2 diabetes in large international cohorts of subjects with monogenic diabetes and >400,000 subjects from UK Biobank. We identified a biallelic frameshift ONECUT1 variant as the cause of NDM in one individual. However, we found no enrichment of missense or null ONECUT1 variants among 484 individuals clinically suspected of MODY, in whom all known genes had been excluded. Finally, using a rare variant burden test in the UK Biobank European cohort, we identified a significant association between heterozygous ONECUT1 null variants and type 2 diabetes (P = 0.006) but did not find an association between missense variants and type 2 diabetes. Our results confirm biallelic ONECUT1 variants as a cause of NDM and highlight monoallelic null variants as a risk factor for type 2 diabetes. These findings confirm the critical role of ONECUT1 in human β-cell function.

ARTICLE HIGHLIGHTS

Second international consensus report on gaps and opportunities for the clinical translation of precision diabetes medicine

Precision medicine is part of the logical evolution of contemporary evidence-based medicine that seeks to reduce errors and optimize outcomes when making medical decisions and health recommendations. Diabetes affects hundreds of millions of people worldwide, many of whom will develop life-threatening complications and die prematurely. Precision medicine can potentially address this enormous problem by accounting for heterogeneity in the etiology, clinical presentation and pathogenesis of common forms of diabetes and risks of complications. This second international consensus report on precision diabetes medicine summarizes the findings from a systematic evidence review across the key pillars of precision medicine (prevention, diagnosis, treatment, prognosis) in four recognized forms of diabetes (monogenic, gestational, type 1, type 2). These reviews address key questions about the translation of precision medicine research into practice. Although not complete, owing to the vast literature on this topic, they revealed opportunities for the immediate or near-term clinical implementation of precision diabetes medicine; furthermore, we expose important gaps in knowledge, focusing on the need to obtain new clinically relevant evidence. Gaps include the need for common standards for clinical readiness, including consideration of cost-effectiveness, health equity, predictive accuracy, liability and accessibility. Key milestones are outlined for the broad clinical implementation of precision diabetes medicine.

Waist circumference and glycaemia are strong predictors of progression to diabetes in individuals with prediabetes in sub-Saharan Africa: 4-year prospective cohort study in Malawi

Sub-Saharan Africa is projected to have the highest increase in the number of people with diabetes worldwide. However, the drivers of diabetes in this region have not been clearly elucidated. The aim of this study was to evaluate the incidence of diabetes and the predictors of progression in a population-based cohort with impaired fasting glucose (IFG) in Malawi. We used data from an extensive rural and urban non-communicable disease survey. One hundred seventy-five, of 389 individuals with impaired fasting glucose (IFG) at baseline, age 48 ±15 years and body mass index 27.5 ±5.9 kg/m2 were followed up for a median of 4.2 years (714 person-years). Incidence rates were calculated, and predictors of progression to diabetes were analysed using multivariable logistic regression models, with overall performance determined using receiver operator characteristics (ROC) curves. The median follow-up was 4.2 (IQR 3.4-4.7) years. Forty-five out of 175 (26%) progressed to diabetes. Incidence rates of diabetes were 62.9 per 1000 person-years 95% CI, 47.0-84.3. The predictors of progression were higher; age (odds ratio [OR] 1.48, P = 0.046), BMI (OR 1.98, P = 0.001), waist circumference (OR 2.50,P<0.001), waist-hip ratio (OR 1.40, P = 0.03), systolic blood pressure (OR 1.56, P = 0.01), fasting plasma glucose (OR 1.53, P = 0.01), cholesterol (OR 1.44, P = 0.05) and low-density lipoprotein cholesterol (OR 1.80, P = 0.002). A simple model combining fasting plasma glucose and waist circumference was predictive of progression to diabetes (ROC area under the curve = 0.79). The incidence of diabetes in people with IFG is high in Malawi and predictors of progression are like those seen in other populations. Our data also suggests that a simple chart with probabilities of progression to diabetes based on waist circumference and fasting plasma glucose could be used to identify those at risk of progression in clinical settings in sub-Saharan Africa.

Systematic Review of Treatment of Beta-Cell Monogenic Diabetes

Background

Beta-cell monogenic forms of diabetes are the area of diabetes care with the strongest support for precision medicine. We reviewed treatment of hyperglycemia in GCK-related hyperglycemia, HNF1A-HNF4A- and HNF1B-diabetes, Mitochondrial diabetes (MD) due to m.3243A>G variant, 6q24-transient neonatal diabetes (TND) and SLC19A2-diabetes.

Methods

Systematic reviews with data from PubMed, MEDLINE and Embase were performed for the different subtypes. Individual and group level data was extracted for glycemic outcomes in individuals with genetically confirmed monogenic diabetes.

Results

147 studies met inclusion criteria with only six experimental studies and the rest being single case reports or cohort studies. Most studies had moderate or serious risk of bias.For GCK-related hyperglycemia, six studies (N=35) showed no deterioration in HbA1c on discontinuing glucose lowering therapy. A randomized trial (n=18 per group) showed that sulfonylureas (SU) were more effective in HNF1A-diabetes than in type 2 diabetes, and cohort and case studies supported SU effectiveness in lowering HbA1c. Two crossover trials (n=15 and n=16) suggested glinides and GLP-1 receptor agonists might be used in place of SU. Evidence for HNF4A-diabetes was limited. While some patients with HNF1B-diabetes (n=301) and MD (n=250) were treated with oral agents, most were on insulin. There was some support for the use of oral agents after relapse in 6q24-TND, and for thiamine improving glycemic control and reducing insulin requirement in SLC19A2-diabetes (less than half achieved insulin-independency).

Conclusion

There is limited evidence to guide the treatment in monogenic diabetes with most studies being non-randomized and small. The data supports: no treatment in GCK-related hyperglycemia; SU for HNF1A-diabetes. Further evidence is needed to examine the optimum treatment in monogenic subtypes.

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

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

Comparison of causal forest and regression-based approaches to evaluate treatment effect heterogeneity: an application for type 2 diabetes precision medicine

OBJECTIVE

Precision medicine requires reliable identification of variation in patient-level outcomes with different available treatments, often termed treatment effect heterogeneity. We aimed to evaluate the comparative utility of individualized treatment selection strategies based on predicted individual-level treatment effects from a causal forest machine learning algorithm and a penalized regression model.

METHODS

Cohort study characterizing individual-level glucose-lowering response (6 month reduction in HbA1c) in people with type 2 diabetes initiating SGLT2-inhibitor or DPP4-inhibitor therapy. Model development set comprised 1,428 participants in the CANTATA-D and CANTATA-D2 randomised clinical trials of SGLT2-inhibitors versus DPP4-inhibitors. For external validation, calibration of observed versus predicted differences in HbA1c in patient strata defined by size of predicted HbA1c benefit was evaluated in 18,741 patients in UK primary care (Clinical Practice Research Datalink).

RESULTS

Heterogeneity in treatment effects was detected in clinical trial participants with both approaches (proportion predicted to have a benefit on SGLT2-inhibitor therapy over DPP4-inhibitor therapy: causal forest: 98.6%; penalized regression: 81.7%). In validation, calibration was good with penalized regression but sub-optimal with causal forest. A strata with an HbA1c benefit > 10 mmol/mol with SGLT2-inhibitors (3.7% of patients, observed benefit 11.0 mmol/mol [95%CI 8.0-14.0]) was identified using penalized regression but not causal forest, and a much larger strata with an HbA1c benefit 5-10 mmol with SGLT2-inhibitors was identified with penalized regression (regression: 20.9% of patients, observed benefit 7.8 mmol/mol (95%CI 6.7-8.9); causal forest 11.6%, observed benefit 8.7 mmol/mol (95%CI 7.4-10.1).

CONCLUSIONS

Consistent with recent results for outcome prediction with clinical data, when evaluating treatment effect heterogeneity researchers should not rely on causal forest or other similar machine learning algorithms alone, and must compare outputs with standard regression, which in this evaluation was superior.

Type 2 Diabetes risk alleles in Peptidyl-glycine Alpha-amidating Monooxygenase influence GLP-1 levels and response to GLP-1 Receptor Agonists

Patients with type 2 diabetes vary in their response to currently available therapeutic agents (including GLP-1 receptor agonists) leading to suboptimal glycemic control and increased risk of complications. We show that human carriers of hypomorphic T2D-risk alleles in the gene encoding peptidyl-glycine alpha-amidating monooxygenase (PAM), as well as Pam-knockout mice, display increased resistance to GLP-1 in vivo. Pam inactivation in mice leads to reduced gastric GLP-1R expression and faster gastric emptying: this persists during GLP-1R agonist treatment and is rescued when GLP-1R activity is antagonized, indicating resistance to GLP-1's gastric slowing properties. Meta-analysis of human data from studies examining GLP-1R agonist response (including RCTs) reveals a relative loss of 44% and 20% of glucose lowering (measured by glycated hemoglobin) in individuals with hypomorphic PAM alleles p.S539W and p.D536G treated with GLP-1R agonist. Genetic variation in PAM has effects on incretin signaling that alters response to medication used commonly for treatment of T2D.

Genetic effects on the timing of parturition and links to fetal birth weight

The timing of parturition is crucial for neonatal survival and infant health. Yet, its genetic basis remains largely unresolved. We present a maternal genome-wide meta-analysis of gestational duration (n = 195,555), identifying 22 associated loci (24 independent variants) and an enrichment in genes differentially expressed during labor. A meta-analysis of preterm delivery (18,797 cases, 260,246 controls) revealed six associated loci and large genetic similarities with gestational duration. Analysis of the parental transmitted and nontransmitted alleles (n = 136,833) shows that 15 of the gestational duration genetic variants act through the maternal genome, whereas 7 act both through the maternal and fetal genomes and 2 act only via the fetal genome. Finally, the maternal effects on gestational duration show signs of antagonistic pleiotropy with the fetal effects on birth weight: maternal alleles that increase gestational duration have negative fetal effects on birth weight. The present study provides insights into the genetic effects on the timing of parturition and the complex maternal-fetal relationship between gestational duration and birth weight.

FOXP3 TSDR Measurement Could Assist Variant Classification and Diagnosis of IPEX Syndrome

Pathogenic FOXP3 variants cause immune dysregulation polyendocrinopathy enteropathy X-linked (IPEX) syndrome, a progressive autoimmune disease resulting from disruption of the regulatory T cell (Treg) compartment. Assigning pathogenicity to novel variants in FOXP3 is challenging due to the heterogeneous phenotype and variable immunological abnormalities. The number of cells with demethylation at the Treg cell-specific demethylated region (TSDR) is an independent biomarker of IPEX. We aimed to investigate if diagnosing IPEX at presentation with isolated diabetes could allow for effective monitoring of disease progression and assess whether TSDR analysis can aid FOXP3 variant classification and predict disease course. We describe a large genetically diagnosed IPEX cohort (n = 65) and 13 individuals with other monogenic autoimmunity subtypes in whom we quantified the proportion of cells with FOXP3 TSDR demethylation, normalized to the number with CD4 demethylation (%TSDR/CD4) and compare them to 29 unaffected controls. IPEX patients presenting with isolated diabetes (50/65, 77%) often later developed enteropathy (20/50, 40%) with a median interval of 23.5 weeks. %TSDR/CD4 was a good discriminator of IPEX vs. unaffected controls (ROC-AUC 0.81, median 13.6% vs. 8.5%, p < 0.0001) with higher levels of demethylation associated with more severe disease. Patients with other monogenic autoimmunity had a similar %TSDR/CD4 to controls (median 8.7%, p = 1.0). Identifying increased %TSDR/CD4 in patients with novel FOXP3 mutations presenting with isolated diabetes facilitates diagnosis and could offer an opportunity to monitor patients and begin immune modulatory treatment before onset of severe enteropathy.

Multi-ancestry genome-wide study in >2.5 million individuals reveals heterogeneity in mechanistic pathways of type 2 diabetes and complications

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes. To characterise the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study (GWAS) data from 2,535,601 individuals (39.7% non-European ancestry), including 428,452 T2D cases. We identify 1,289 independent association signals at genome-wide significance (P<5×10-8) that map to 611 loci, of which 145 loci are previously unreported. We define eight non-overlapping clusters of T2D signals characterised by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial, and enteroendocrine cells. We build cluster-specific partitioned genetic risk scores (GRS) in an additional 137,559 individuals of diverse ancestry, including 10,159 T2D cases, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned GRS are more strongly associated with coronary artery disease and end-stage diabetic nephropathy than an overall T2D GRS across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings demonstrate the value of integrating multi-ancestry GWAS with single-cell epigenomics to disentangle the aetiological heterogeneity driving the development and progression of T2D, which may offer a route to optimise global access to genetically-informed diabetes care.

Infancy-onset diabetes caused by de-regulated AMPylation of the human endoplasmic reticulum chaperone BiP

Dysfunction of the endoplasmic reticulum (ER) in insulin-producing beta cells results in cell loss and diabetes mellitus. Here we report on five individuals from three different consanguineous families with infancy-onset diabetes mellitus and severe neurodevelopmental delay caused by a homozygous p.(Arg371Ser) mutation in FICD. The FICD gene encodes a bifunctional Fic domain-containing enzyme that regulates the ER Hsp70 chaperone, BiP, via catalysis of two antagonistic reactions: inhibitory AMPylation and stimulatory deAMPylation of BiP. Arg371 is a conserved residue in the Fic domain active site. The FICD^R371S mutation partially compromises BiP AMPylation in vitro but eliminates all detectable deAMPylation activity. Overexpression of FICD^R371S or knock-in of the mutation at the FICD locus of stressed CHO cells results in inappropriately elevated levels of AMPylated BiP and compromised secretion. These findings, guided by human genetics, highlight the destructive consequences of de-regulated BiP AMPylation and raise the prospect of tuning FICD's antagonistic activities towards therapeutic ends.

Patient stratification for determining optimal second-line and third-line therapy for type 2 diabetes: the TriMaster study

Precision medicine aims to treat an individual based on their clinical characteristics. A differential drug response, critical to using these features for therapy selection, has never been examined directly in type 2 diabetes. In this study, we tested two hypotheses: (1) individuals with body mass index (BMI) > 30 kg/m², compared to BMI ≤ 30 kg/m², have greater glucose lowering with thiazolidinediones than with DPP4 inhibitors, and (2) individuals with estimated glomerular filtration rate (eGFR) 60-90 ml/min/1.73 m², compared to eGFR >90 ml/min/1.73 m², have greater glucose lowering with DPP4 inhibitors than with SGLT2 inhibitors. The primary endpoint for both hypotheses was the achieved HbA1c difference between strata for the two drugs. In total, 525 people with type 2 diabetes participated in this UK-based randomized, double-blind, three-way crossover trial of 16 weeks of treatment with each of sitagliptin 100 mg once daily, canagliflozin 100 mg once daily and pioglitazone 30 mg once daily added to metformin alone or metformin plus sulfonylurea. Overall, the achieved HbA1c was similar for the three drugs: pioglitazone 59.6 mmol/mol, sitagliptin 60.0 mmol/mol and canagliflozin 60.6 mmol/mol (P = 0.2). Participants with BMI > 30 kg/m², compared to BMI ≤ 30 kg/m², had a 2.88 mmol/mol (95% confidence interval (CI): 0.98, 4.79) lower HbA1c on pioglitazone than on sitagliptin (n = 356, P = 0.003). Participants with eGFR 60-90 ml/min/1.73 m², compared to eGFR >90 ml/min/1.73 m², had a 2.90 mmol/mol (95% CI: 1.19, 4.61) lower HbA1c on sitagliptin than on canagliflozin (n = 342, P = 0.001). There were 2,201 adverse events reported, and 447/525 (85%) randomized participants experienced an adverse event on at least one of the study drugs. In this precision medicine trial in type 2 diabetes, our findings support the use of simple, routinely available clinical measures to identify the drug class most likely to deliver the greatest glycemic reduction for a given patient. (ClinicalTrials.gov registration: NCT02653209 ; ISRCTN registration: 12039221 .).

Patient preference for second- and third-line therapies in type 2 diabetes: a prespecified secondary endpoint of the TriMaster study

Patient preference is very important for medication selection in chronic medical conditions, like type 2 diabetes, where there are many different drugs available. Patient preference balances potential efficacy with potential side effects. As both aspects of drug response can vary markedly between individuals, this decision could be informed by the patient personally experiencing the alternative medications, as occurs in a crossover trial. In the TriMaster (NCT02653209, ISRCTN12039221), randomized double-blind, three-way crossover trial patients received three different second- or third-line once-daily type 2 diabetes glucose-lowering drugs (pioglitazone 30 mg, sitagliptin 100 mg and canagliflozin 100 mg). As part of a prespecified secondary endpoint, we examined patients' drug preference after they had tried all three drugs. In total, 448 participants were treated with all three drugs which overall showed similar glycemic control (HbA1c on pioglitazone 59.5 sitagliptin 59.9, canagliflozin 60.5 mmol mol^-1, P = 0.19). In total, 115 patients (25%) preferred pioglitazone, 158 patients (35%) sitagliptin and 175 patients (38%) canagliflozin. The drug preferred by individual patients was associated with a lower HbA1c (mean: 4.6; 95% CI: 3.9, 5.3) mmol mol^-1 lower versus nonpreferred) and fewer side effects (mean: 0.50; 95% CI: 0.35, 0.64) fewer side effects versus nonpreferred). Allocating therapy based on the individually preferred drugs, rather than allocating all patients the overall most preferred drug (canagliflozin), would result in more patients achieving the lowest HbA1c for them (70% versus 30%) and the fewest side effects (67% versus 50%). When precision approaches do not predict a clear optimal therapy for an individual, allowing patients to try potential suitable medications before they choose long-term therapy could be a practical alternative to optimizing treatment for type 2 diabetes.

The relationship between islet autoantibody status and the genetic risk of type 1 diabetes in adult-onset type 1 diabetes

AIMS/HYPOTHESIS

The reason for the observed lower rate of islet autoantibody positivity in clinician-diagnosed adult-onset vs childhood-onset type 1 diabetes is not known. We aimed to explore this by assessing the genetic risk of type 1 diabetes in autoantibody-negative and -positive children and adults.

METHODS

We analysed GAD autoantibodies, insulinoma-2 antigen autoantibodies and zinc transporter-8 autoantibodies (ZnT8A) and measured type 1 diabetes genetic risk by genotyping 30 type 1 diabetes-associated variants at diagnosis in 1814 individuals with clinician-diagnosed type 1 diabetes (1112 adult-onset, 702 childhood-onset). We compared the overall type 1 diabetes genetic risk score (T1DGRS) and non-HLA and HLA (DR3-DQ2, DR4-DQ8 and DR15-DQ6) components with autoantibody status in those with adult-onset and childhood-onset diabetes. We also measured the T1DGRS in 1924 individuals with type 2 diabetes from the Wellcome Trust Case Control Consortium to represent non-autoimmune diabetes control participants.

RESULTS

The T1DGRS was similar in autoantibody-negative and autoantibody-positive clinician-diagnosed childhood-onset type 1 diabetes (mean [SD] 0.274 [0.034] vs 0.277 [0.026], p=0.4). In contrast, the T1DGRS in autoantibody-negative adult-onset type 1 diabetes was lower than that in autoantibody-positive adult-onset type 1 diabetes (mean [SD] 0.243 [0.036] vs 0.271 [0.026], p<0.0001) but higher than that in type 2 diabetes (mean [SD] 0.229 [0.034], p<0.0001). Autoantibody-negative adults were more likely to have the more protective HLA DR15-DQ6 genotype (15% vs 3%, p<0.0001), were less likely to have the high-risk HLA DR3-DQ2/DR4-DQ8 genotype (6% vs 19%, p<0.0001) and had a lower non-HLA T1DGRS (p<0.0001) than autoantibody-positive adults. In contrast to children, autoantibody-negative adults were more likely to be male (75% vs 59%), had a higher BMI (27 vs 24 kg/m²) and were less likely to have other autoimmune conditions (2% vs 10%) than autoantibody-positive adults (all p<0.0001). In both adults and children, type 1 diabetes genetic risk was unaffected by the number of autoantibodies (p>0.3). These findings, along with the identification of seven misclassified adults with monogenic diabetes among autoantibody-negative adults and the results of a sensitivity analysis with and without measurement of ZnT8A, suggest that the intermediate type 1 diabetes genetic risk in autoantibody-negative adults is more likely to be explained by the inclusion of misclassified non-autoimmune diabetes (estimated to represent 67% of all antibody-negative adults, 95% CI 61%, 73%) than by the presence of unmeasured autoantibodies or by a discrete form of diabetes. When these estimated individuals with non-autoimmune diabetes were adjusted for, the prevalence of autoantibody positivity in adult-onset type 1 diabetes was similar to that in children (93% vs 91%, p=0.4).

CONCLUSIONS/INTERPRETATION

The inclusion of non-autoimmune diabetes is the most likely explanation for the observed lower rate of autoantibody positivity in clinician-diagnosed adult-onset type 1 diabetes. Our data support the utility of islet autoantibody measurement in clinician-suspected adult-onset type 1 diabetes in routine clinical practice.

The impact of population-level HbA1c screening on reducing diabetes diagnostic delay in middle-aged adults: a UK Biobank analysis

AIMS/HYPOTHESIS

Screening programmes can detect cases of undiagnosed diabetes earlier than symptomatic or incidental diagnosis. However, the improvement in time to diagnosis achieved by screening programmes compared with routine clinical care is unclear. We aimed to use the UK Biobank population-based study to provide the first population-based estimate of the reduction in time to diabetes diagnosis that could be achieved by HbA_1c-based screening in middle-aged adults.

METHODS

We studied UK Biobank participants aged 40-70 years with HbA_1c measured at enrolment (but not fed back to participants/clinicians) and linked primary and secondary healthcare data (n=179,923) and identified those with a pre-existing diabetes diagnosis (n=13,077, 7.3%). Among the remaining participants (n=166,846) without a diabetes diagnosis, we used an elevated enrolment HbA_1c level (≥48 mmol/mol [≥6.5%]) to identify those with undiagnosed diabetes. For this group, we used Kaplan-Meier analysis to assess the time between enrolment HbA_1c measurement and subsequent clinical diabetes diagnosis up to 10 years, and Cox regression to identify clinical factors associated with delayed diabetes diagnosis.

RESULTS

In total, 1.0% (1703/166,846) of participants without a diabetes diagnosis had undiagnosed diabetes based on calibrated HbA_1c levels at UK Biobank enrolment, with a median HbA_1c level of 51.3 mmol/mol (IQR 49.1-57.2) (6.8% [6.6-7.4]). These participants represented an additional 13.0% of diabetes cases in the study population relative to the 13,077 participants with a diabetes diagnosis. The median time to clinical diagnosis for those with undiagnosed diabetes was 2.2 years, with a median HbA_1c at clinical diagnosis of 58.2 mmol/mol (IQR 51.0-80.0) (7.5% [6.8-9.5]). Female participants with lower HbA_1c and BMI measurements at enrolment experienced the longest delay to clinical diagnosis.

CONCLUSIONS/INTERPRETATION

Our population-based study shows that HbA_1c screening in adults aged 40-70 years can reduce the time to diabetes diagnosis by a median of 2.2 years compared with routine clinical care. The findings support the use of HbA_1c screening to reduce the time for which individuals are living with undiagnosed diabetes.

Identifying type 1 and 2 diabetes in research datasets where classification biomarkers are unavailable: assessing the accuracy of published approaches

OBJECTIVES

We aimed to compare the performance of approaches for classifying insulin-treated diabetes within research datasets without measured classification biomarkers, evaluated against two independent biological definitions of diabetes type.

STUDY DESIGN AND SETTING

We compared accuracy of ten reported approaches for classifying insulin-treated diabetes into type 1 (T1D) and type 2 (T2D) diabetes in two cohorts: UK Biobank (UKBB) n = 26,399 and Diabetes Alliance for Research in England (DARE) n = 1,296. The overall performance for classifying T1D and T2D was assessed using: a T1D genetic risk score and genetic stratification method (UKBB); C-peptide measured at >3 years diabetes duration (DARE).

RESULTS

Approaches' accuracy ranged from 71% to 88% (UKBB) and 68% to 88% (DARE). When classifying all participants, combining early insulin requirement with a T1D probability model (incorporating diagnosis age and body image issue [BMI]), and interview-reported diabetes type (UKBB available in only 15%) consistently achieved high accuracy (UKBB 87% and 87% and DARE 85% and 88%, respectively). For identifying T1D with minimal misclassification, models with high thresholds or young diagnosis age (<20 years) had highest performance. Findings were incorporated into an online tool identifying optimum approaches based on variable availability.

CONCLUSION

Models combining continuous features with early insulin requirement are the most accurate methods for classifying insulin-treated diabetes in research datasets without measured classification biomarkers.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Development of a treatment selection algorithm for SGLT2 and DPP-4 inhibitor therapies in people with type 2 diabetes: a retrospective cohort study

BACKGROUND

Current treatment guidelines do not provide recommendations to support the selection of treatment for most people with type 2 diabetes. We aimed to develop and validate an algorithm to allow selection of optimal treatment based on glycaemic response, weight change, and tolerability outcomes when choosing between SGLT2 inhibitor or DPP-4 inhibitor therapies.

METHODS

In this retrospective cohort study, we identified patients initiating SGLT2 and DPP-4 inhibitor therapies after Jan 1, 2013, from the UK Clinical Practice Research Datalink (CPRD). We excluded those who received SGLT2 or DPP-4 inhibitors as first-line treatment or insulin at the same time, had estimated glomerular filtration rate (eGFR) of less than 45 mL/min per 1·73 m², or did not have a valid baseline glycated haemoglobin (HbA_1c) measure (<53 or ≥120 mmol/mol). The primary efficacy outcome was the HbA_1c value reached 6 months after drug initiation, adjusted for baseline HbA_1c. Clinical features associated with differential HbA_1c outcome on the two therapies were identified in CPRD (n=26 877), and replicated in reanalysis of 14 clinical trials (n=10 414). An algorithm to predict individual-level differential HbA_1c outcome on the two therapies was developed in CPRD (derivation; n=14 069) and validated in head-to-head trials (n=2499) and CPRD (independent validation; n=9376). In CPRD, we further explored heterogeneity in 6-month weight change and treatment discontinuation.

FINDINGS

Among 10 253 patients initiating SGLT2 inhibitors and 16 624 patients initiating DPP-4 inhibitors in CPRD, baseline HbA_1c, age, BMI, eGFR, and alanine aminotransferase were associated with differential HbA_1c outcome with SGLT2 inhibitor and DPP-4 inhibitor therapies. The median age of participants was 62·0 years (IQR 55·0-70·0). 10 016 (37·3%) were women and 16 861 (62·7%) were men. An algorithm based on these five features identified a subgroup, representing around four in ten CPRD patients, with a 5 mmol/mol or greater observed benefit with SGLT2 inhibitors in all validation cohorts (CPRD 8·8 mmol/mol [95% CI 7·8-9·8]; CANTATA-D and CANTATA-D2 trials 5·8 mmol/mol [3·9-7·7]; BI1245.20 trial 6·6 mmol/mol [2·2-11·0]). In CPRD, predicted differential HbA_1c response with SGLT2 inhibitor and DPP-4 inhibitor therapies was not associated with weight change. Overall treatment discontinuation within 6 months was similar in patients predicted to have an HbA_1c benefit with SGLT2 inhibitors over DPP-4 inhibitors (median 15·2% [13·2-20·3] vs 14·4% [12·9-16·7]). A smaller subgroup predicted to have greater HbA_1c reduction with DPP-4 inhibitors were twice as likely to discontinue SGLT2 inhibitors than DPP-4 inhibitors (median 26·8% [23·4-31·0] vs 14·8% [12·9-16·8]).

INTERPRETATION

A validated treatment selection algorithm for SGLT2 inhibitor and DPP-4 inhibitor therapies can support decisions on optimal treatment for people with type 2 diabetes.

FUNDING

BHF-Turing Cardiovascular Data Science Award and the UK Medical Research Council.

Routine Islet Autoantibody Testing in Clinically Diagnosed Adult-Onset Type 1 Diabetes Can Help Identify Misclassification and the Possibility of Successful Insulin Cessation

OBJECTIVE

Recent joint American Diabetes Association and European Association for the Study of Diabetes guidelines recommend routine islet autoantibody testing in all adults newly diagnosed with type 1 diabetes. We aimed to assess the impact of routine islet autoantibody testing in this population.

RESEARCH DESIGN AND METHODS

We prospectively assessed the relationship between islet autoantibody status (GADA, IA-2A, and ZNT8A), clinical and genetic characteristics, and progression (annual change in urine C-peptide-to-creatinine ratio [UCPCR]) in 722 adults (≥18 years old at diagnosis) with clinically diagnosed type 1 diabetes and diabetes duration <12 months. We also evaluated changes in treatment and glycemia over 2 years after informing participants and their clinicians of autoantibody results.

RESULTS

Of 722 participants diagnosed with type 1 diabetes, 24.8% (179) were autoantibody negative. This group had genetic and C-peptide characteristics suggestive of a high prevalence of nonautoimmune diabetes: lower mean type 1 diabetes genetic risk score (islet autoantibody negative vs. positive: 10.85 vs. 13.09 [P < 0.001] [type 2 diabetes 10.12]) and lower annual change in C-peptide (UCPCR), -24% vs. -43% (P < 0.001).After median 24 months of follow-up, treatment change occurred in 36.6% (60 of 164) of autoantibody-negative participants: 22.6% (37 of 164) discontinued insulin, with HbA1c similar to that of participants continuing insulin (57.5 vs. 60.8 mmol/mol [7.4 vs. 7.7%], P = 0.4), and 14.0% (23 of 164) added adjuvant agents to insulin.

CONCLUSIONS

In adult-onset clinically diagnosed type 1 diabetes, negative islet autoantibodies should prompt careful consideration of other diabetes subtypes. When routinely measured, negative antibodies are associated with successful insulin cessation. These findings support recent recommendations for routine islet autoantibody assessment in adult-onset type 1 diabetes.

Reduced penetrance of MODY-associated HNF1A/HNF4A variants but not GCK variants in clinically unselected cohorts

The true prevalence and penetrance of monogenic disease variants are often not known because of clinical-referral ascertainment bias. We comprehensively assess the penetrance and prevalence of pathogenic variants in HNF1A, HNF4A, and GCK that account for >80% of monogenic diabetes. We analyzed clinical and genetic data from 1,742 clinically referred probands, 2,194 family members, clinically unselected individuals from a US health system-based cohort (n = 132,194), and a UK population-based cohort (n = 198,748). We show that one in 1,500 individuals harbor a pathogenic variant in one of these genes. The penetrance of diabetes for HNF1A and HNF4A pathogenic variants was substantially lower in the clinically unselected individuals compared to clinically referred probands and was dependent on the setting (32% in the population, 49% in the health system cohort, 86% in a family member, and 98% in probands for HNF1A). The relative risk of diabetes was similar across the clinically unselected cohorts highlighting the role of environment/other genetic factors. Surprisingly, the penetrance of pathogenic GCK variants was similar across all cohorts (89%-97%). We highlight that pathogenic variants in HNF1A, HNF4A, and GCK are not ultra-rare in the population. For HNF1A and HNF4A, we need to tailor genetic interpretation and counseling based on the setting in which a pathogenic monogenic variant was identified. GCK is an exception with near-complete penetrance in all settings. This along with the clinical implication of diagnosis makes it an excellent candidate for the American College of Medical Genetics secondary gene list.

A biallelic loss-of-function PDIA6 variant in a second patient with polycystic kidney disease, infancy-onset diabetes, and microcephaly

We report a second patient with intrauterine growth retardation, congenital polycystic kidney disease, infancy-onset diabetes, microcephaly, and liver fibrosis caused by a homozygous PDIA6 loss-of-function variant. Our study further defines the genetic and clinical features of this rare syndromic form of infancy-onset diabetes.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Patient-led rapid titration of basal insulin in gestational diabetes is associated with improved glycaemic control and lower birthweight

AIMS

Elevated fasting blood glucose in gestational diabetes (GDM) is a key predictor of high birthweight babies and adverse pregnancy outcomes but is hard to treat. We implemented a simple, patient-led, insulin dose titration algorithm aiming to improve fasting glycaemic control in GDM.

METHODS

In women with GDM, initiating basal insulin, we recommended a daily four-unit dose increase after every fasting glucose value ≥5.0 mmol/mol (90 mg/dl). This approach augmented our pre-existing intensive (weekly) specialist nursing input. Using a before-and-after retrospective observational study design, we examined insulin doses and glucose values at 36 weeks gestation and maternal and neonatal outcomes in 105 women completing pregnancy before and 93 women after the intervention.

RESULTS

The baseline characteristics of women in the before and after groups were the same. Women initiated on insulin after implementation (n = 30 before, n = 43 after) achieved substantially higher doses at 36 weeks (53 vs. 36 units/day; 0.56 vs. 0.37 units/kg/day; p = 0.027). 36-week mean fasting glucose was lower in those on insulin after implementation (4.6 vs. 5.1 mmol/L [83 vs. 92 mg/dl]; p = 0.031). Birthweight was significantly reduced (birthweight Z-scores 0.34 vs. 0.92; p = 0.005). There was no significant difference in macrosomia (after; 2% vs. before; 17% p = 0.078) or caesarean sections (after; 33% vs. before; 47%; p = 0.116). No women experienced severe hypoglycaemia. There were no outcome differences before versus after intervention in women not treated with insulin.

CONCLUSIONS

Patient-led daily insulin titration in gestational diabetes leads to higher insulin dose use lower fasting glucose and is associated with lower birthweight without causing significant hypoglycaemia.

Recurrent 17q12 microduplications contribute to renal disease but not diabetes

BACKGROUND

17q12 microdeletion and microduplication syndromes present as overlapping, multisystem disorders. We assessed the disease phenotypes of individuals with 17q12 CNV in a population-based cohort.

METHODS

We investigated 17q12 CNV using microarray data from 450 993 individuals in the UK Biobank and calculated disease status associations for diabetes, liver and renal function, neurological and psychiatric traits.

RESULTS

We identified 11 17q12 microdeletions and 106 microduplications. Microdeletions were strongly associated with diabetes (p=2×10^-7) but microduplications were not. Estimated glomerular filtration rate (eGFR mL/min/1.73 m²) was consistently lower in individuals with microdeletions (p=3×10^-12) and microduplications (p=6×10^-25). Similarly, eGFR <60, including end-stage renal disease, was associated with microdeletions (p=2×10^-9, p<0.003) and microduplications (p=1×10^-9, p=0.009), respectively, highlighting sometimes substantially reduced renal function in each. Microduplications were associated with decreased fluid intelligence (p=3×10^-4). SNP association analysis in the 17q12 region implicated changes to HNF1B as causing decreased eGFR (NC_000017.11:g.37741642T>G, rs12601991, p=4×10^-21) and diabetes (NC_000017.11:g.37741165C>T, rs7501939, p=6×10^-17). A second locus within the region was also associated with fluid intelligence (NC_000017.11:g.36593168T>C, rs1005552, p=6×10^-9) and decreased eGFR (NC_000017.11:g.36558947T>C, rs12150665, p=4×10^-15).

CONCLUSION

We demonstrate 17q12 microdeletions but not microduplications are associated with diabetes in a population-based cohort, likely caused by HNF1B haploinsufficiency. We show that both 17q12 microdeletions and microduplications are associated with renal disease, and multiple genes within the region likely contribute to renal and neurocognitive phenotypes.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Congenital beta cell defects are not associated with markers of islet autoimmunity, even in the context of high genetic risk for type 1 diabetes

AIMS/HYPOTHESIS

A key unanswered question in type 1 diabetes is whether beta cells initiate their own destruction or are victims of an aberrant immune response (beta cell suicide or homicide?). To investigate this, we assessed islet autoantibodies in individuals with congenital beta cell defects causing neonatal diabetes mellitus (NDM).

METHODS

We measured autoantibodies to GAD (GADA), islet antigen-2 (IA-2A) and zinc transporter 8 (ZnT8A) in 242 individuals with NDM (median age diagnosed 1.8 months [IQR 0.39-2.9 months]; median age collected 4.6 months [IQR 1.8-27.6 months]; median diabetes duration 2 months [IQR 0.6-23 months]), including 75 whose NDM resulted from severe beta cell endoplasmic reticulum (ER) stress. As a control cohort we also tested samples from 69 diabetes-free individuals (median age collected 9.9 months [IQR 9.0-48.6 months]) for autoantibodies.

RESULTS

We found low prevalence of islet autoantibodies in individuals with monogenic NDM; 13/242 (5.4% [95% CI 2.9, 9.0%]) had detectable GADA, IA-2A and/or ZnT8A. This was similar to the proportion in the control participants who did not have diabetes (1/69 positive [1.4%, 95% CI 0.03, 7.8%], p=0.3). Importantly, monogenic individuals with beta cell ER stress had a similar rate of GADA/IA-2A/ZnT8A positivity to non-ER stress aetiologies (2.7% [95% CI 0.3, 9.3%] vs 6.6% [95% CI 3.3, 11.5%] p=0.4). We observed no association between islet autoimmunity and genetic risk, age at testing (including 30 individuals >10 years at testing) or diabetes duration (p>0.4 for all).

CONCLUSIONS/INTERPRETATION

Our data support the hypothesis that beta cell stress/dysfunction alone does not lead to the production of islet autoantibodies, even in the context of high-risk HLA types. This suggests that additional factors are required to trigger an autoimmune response towards beta cells.

Fetal alleles predisposing to metabolically favorable adiposity are associated with higher birth weight

BACKGROUND

Higher birthweight is associated with higher adult body mass index (BMI). Alleles that predispose to greater adult adiposity might act in fetal life to increase fetal growth and birthweight. Whether there are fetal effects of recently identified adult metabolically favorable adiposity alleles on birthweight is unknown.

AIM

We aimed to test the effect on birthweight of fetal genetic predisposition to higher metabolically favorable adult adiposity and compare that with the effect of fetal genetic predisposition to higher adult BMI.

METHODS

We used published genome wide association study data (n = upto 406 063) to estimate fetal effects on birthweight (adjusting for maternal genotype) of alleles known to raise metabolically favorable adult adiposity or BMI. We combined summary data across single nucleotide polymorphisms (SNPs) with random effects meta-analyses. We performed weighted linear regression of SNP-birthweight effects against SNP-adult adiposity effects to test for a dose-dependent association.

RESULTS

Fetal genetic predisposition to higher metabolically favorable adult adiposity and higher adult BMI were both associated with higher birthweight (3 g per effect allele (95% CI: 1-5) averaged over 14 SNPs; P = 0.002; 0.5 g per effect allele (95% CI: 0-1) averaged over 76 SNPs; P = 0.042, respectively). SNPs with greater effects on metabolically favorable adiposity tended to have greater effects on birthweight (R2 = 0.2912, P = 0.027). There was no dose-dependent association for BMI (R2 = -0.0019, P = 0.602).

CONCLUSIONS

Fetal genetic predisposition to both higher adult metabolically favorable adiposity and BMI is associated with birthweight. Fetal effects of metabolically favorable adiposity-raising alleles on birthweight are modestly proportional to their effects on future adiposity, but those of BMI-raising alleles are not.

Islet autoantibody positivity in an adult population with recently diagnosed diabetes in Uganda

Aims

This study aimed to investigate the frequency of islet autoantibody positivity in adult patients with recently diagnosed diabetes in Uganda and its associated characteristics.

Methods

Autoantibodies to glutamic acid decarboxylase-65 (GADA), zinc transporter 8 (ZnT8-A), and tyrosine phosphatase (IA-2A) were measured in 534 adult patients with recently diagnosed diabetes. Islet autoantibody positivity was defined based on diagnostic thresholds derived from a local adult population without diabetes. The socio-demographic, clinical, and metabolic characteristics of islet autoantibody-positive and negative participants were then compared. The differences in these characteristics were analysed using the x² test for categorical data and the Kruskal Wallis test for continuous data. Multivariate analysis was performed to identify predictors of islet autoantibody positivity.

Results

Thirty four (6.4%) participants were positive for ≥1 islet autoantibody. GADA, IA-2A and ZnT8-A positivity was detected in 17 (3.2%), 10 (1.9%), and 7 (1.3%) participants, respectively. Compared with those negative for islet autoantibodies, participants positive for islet autoantibodies were more likely to live in a rural area (n = 18, 52.9% Vs n = 127, 25.5%, p = 0.005), to be initiated on insulin therapy (n = 19, 55.9% Vs n = 134, 26.8%, p<0.001), to have a lower median waist circumference (90 [80–99] cm Vs 96 [87–104.8], p = 0.04), waist circumference: height ratio (0.55 [0.50–0.63] vs 0.59 [0.53–0.65], p = 0.03), and fasting C-peptide concentration (0.9 [0.6–1.8] Vs 1.4 [0.8–2.1] ng/ml, p = 0.01). On multivariate analysis, living in a rural area (odds ratio or OR 3.62, 95%CI 1.68–7.80, p = 0.001) and being initiated on insulin therapy (OR 3.61, 95% CI 1.67–7.83, p = 0.001) were associated with islet autoantibody positivity.

Conclusion

The prevalence of islet autoantibody positivity was relatively low, suggesting that pancreatic autoimmunity is a rare cause of new-onset diabetes in this adult Ugandan population. Living in a rural area and being initiated on insulin therapy were independently associated with islet autoantibody positivity in this study population.

PLIN1 Haploinsufficiency Causes a Favorable Metabolic Profile

CONTEXT

PLIN1 encodes perilipin-1, which coats lipid droplets in adipocytes and is involved in droplet formation, triglyceride storage, and lipolysis. Rare PLIN1 frameshift variants that extend the translated protein have been described to cause lipodystrophy.

OBJECTIVE

This work aimed to test whether PLIN1 protein-truncating variants (PTVs) cause lipodystrophy in a large population-based cohort.

METHODS

We identified individuals with PLIN1 PTVs in individuals with exome data in the UK Biobank. We performed gene-burden testing for individuals with PLIN1 PTVs. We replicated the associations using data from the T2D Knowledge portal. We performed a phenome-wide association study using publicly available association statistics. A total of 362 791 individuals in the UK Biobank, a population-based cohort, and 43 125 individuals in the T2D Knowledge portal, a type 2 diabetes (T2D) case-control study, were included in the analyses. Main outcome measures included 22 diseases and traits relevant to lipodystrophy.

RESULTS

The 735 individuals with PLIN1 PTVs had a favorable metabolic profile. These individuals had increased high-density lipoprotein cholesterol (0.12 mmol/L; 95% CI, 0.09 to 0.14, P = 2 × 10-18), reduced triglycerides (-0.22 mmol/L; 95% CI, -0.29 to -0.14, P = 3 × 10-11), reduced waist-to-hip ratio (-0.02; 95% CI, -0.02 to -0.01, P = 9 × 10-12), and reduced systolic blood pressure (-1.67 mm Hg; 95% CI, -3.25 to -0.09, P = .05). These associations were consistent in the smaller T2D Knowledge portal cohort. In the UK Biobank, PLIN1 PTVs were associated with reduced risk of myocardial infarction (odds ratio [OR] = 0.59; 95% CI, 0.35 to 0.93, P = .02) and hypertension (OR = 0.85; 95% CI, 0.73 to 0.98, P = .03), but not T2D (OR = 0.99; 95% CI, 0.63-1.51, P = .99).

CONCLUSION

Our study suggests that PLIN1 haploinsufficiency causes a favorable metabolic profile and may protect against cardiovascular disease.

Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10-9), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background.