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.

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.

Safety and effectiveness of SGLT2 inhibitors in a UK population with type 2 diabetes and aged over 70 years: an instrumental variable approach.

AIMS/HYPOTHESIS: Older adults are under-represented in trials, meaning the benefits and risks of glucose-lowering agents in this age group are unclear. The aim of this study was to assess the safety and effectiveness of sodium-glucose cotransporter 2 inhibitors (SGLT2i) in people with type 2 diabetes aged over 70 years using causal analysis. METHODS: Hospital-linked UK primary care data (Clinical Practice Research Datalink, 2013-2020) were used to compare adverse events and effectiveness in individuals initiating SGLT2i compared with dipeptidyl peptidase-4 inhibitors (DPP4i). Analysis was age-stratified: <70 years (SGLT2i n=66,810, DPP4i n=76,172), ≥70 years (SGLT2i n=10,419, DPP4i n=33,434). Outcomes were assessed using the instrumental variable causal inference method and prescriber preference as the instrument. RESULTS: Risk of diabetic ketoacidosis was increased with SGLT2i in those aged ≥70 (incidence rate ratio compared with DPP4i: 3.82 [95% CI 1.12, 13.03]), but not in those aged <70 (1.12 [0.41, 3.04]). However, incidence rates with SGLT2i in those ≥70 was low (29.6 [29.5, 29.7]) per 10,000 person-years. SGLT2i were associated with similarly increased risk of genital infection in both age groups (incidence rate ratio in those <70: 2.27 [2.03, 2.53]; ≥70: 2.16 [1.77, 2.63]). There was no evidence of an increased risk of volume depletion, poor micturition control, urinary frequency, falls or amputation with SGLT2i in either age group. In those ≥70, HbA1c reduction was similar between SGLT2i and DPP4i (-0.3 mmol/mol [-1.6, 1.1], -0.02% [0.1, 0.1]), but in those <70, SGLT2i were more effective (-4 mmol/mol [4.8, -3.1], -0.4% [-0.4, -0.3]). CONCLUSIONS/INTERPRETATION: Causal analysis suggests SGLT2i are effective in adults aged ≥70 years, but increase risk for genital infections and diabetic ketoacidosis. Our study extends RCT evidence to older adults with type 2 diabetes.

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.

Developmentally dynamic changes in DNA methylation in the human pancreas.

Development of the human pancreas requires the precise temporal control of gene expression via epigenetic mechanisms and the binding of key transcription factors. We quantified genome-wide patterns of DNA methylation in human fetal pancreatic samples from donors aged 6 to 21 post-conception weeks. We found dramatic changes in DNA methylation across pancreas development, with > 21% of sites characterized as developmental differentially methylated positions (dDMPs) including many annotated to genes associated with monogenic diabetes. An analysis of DNA methylation in postnatal pancreas tissue showed that the dramatic temporal changes in DNA methylation occurring in the developing pancreas are largely limited to the prenatal period. Significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small proportion of sites showing sex-specific DNA methylation trajectories across pancreas development. Pancreas dDMPs were not distributed equally across the genome and were depleted in regulatory domains characterized by open chromatin and the binding of known pancreatic development transcription factors. Finally, we compared our pancreas dDMPs to previous findings from the human brain, identifying evidence for tissue-specific developmental changes in DNA methylation. This study represents the first systematic exploration of DNA methylation patterns during human fetal pancreas development and confirms the prenatal period as a time of major epigenomic plasticity.

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.

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.

Comparison of Bayesian approaches for developing prediction models in rare disease: application to the identification of patients with Maturity-Onset Diabetes of the Young.

BACKGROUND: Clinical prediction models can help identify high-risk patients and facilitate timely interventions. However, developing such models for rare diseases presents challenges due to the scarcity of affected patients for developing and calibrating models. Methods that pool information from multiple sources can help with these challenges. METHODS: We compared three approaches for developing clinical prediction models for population screening based on an example of discriminating a rare form of diabetes (Maturity-Onset Diabetes of the Young - MODY) in insulin-treated patients from the more common Type 1 diabetes (T1D). Two datasets were used: a case-control dataset (278 T1D, 177 MODY) and a population-representative dataset (1418 patients, 96 MODY tested with biomarker testing, 7 MODY positive). To build a population-level prediction model, we compared three methods for recalibrating models developed in case-control data. These were prevalence adjustment ("offset"), shrinkage recalibration in the population-level dataset ("recalibration"), and a refitting of the model to the population-level dataset ("re-estimation"). We then developed a Bayesian hierarchical mixture model combining shrinkage recalibration with additional informative biomarker information only available in the population-representative dataset. We developed a method for dealing with missing biomarker and outcome information using prior information from the literature and other data sources to ensure the clinical validity of predictions for certain biomarker combinations. RESULTS: The offset, re-estimation, and recalibration methods showed good calibration in the population-representative dataset. The offset and recalibration methods displayed the lowest predictive uncertainty due to borrowing information from the fitted case-control model. We demonstrate the potential of a mixture model for incorporating informative biomarkers, which significantly enhanced the model's predictive accuracy, reduced uncertainty, and showed higher stability in all ranges of predictive outcome probabilities. CONCLUSION: We have compared several approaches that could be used to develop prediction models for rare diseases. Our findings highlight the recalibration mixture model as the optimal strategy if a population-level dataset is available. This approach offers the flexibility to incorporate additional predictors and informed prior probabilities, contributing to enhanced prediction accuracy for rare diseases. It also allows predictions without these additional tests, providing additional information on whether a patient should undergo further biomarker testing before genetic testing.

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 m2) 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.

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.