Genome-wide aggregated trans-effects on risk of type 1 diabetes: A test of the "omnigenic" sparse effector hypothesis of complex trait genetics.

The "omnigenic" hypothesis postulates that the polygenic effects of common SNPs on a typical complex trait are mediated through trans-effects on expression of a relatively sparse set of effector ("core") genes. We tested this hypothesis in a study of 4,964 cases of type 1 diabetes (T1D) and 7,497 controls by using summary statistics to calculate aggregated (excluding the HLA region) trans-scores for gene expression in blood. From associations of T1D with aggregated trans-scores, nine putative core genes were identified, of which three-STAT1, CTLA4 and FOXP3-are genes in which variants cause monogenic forms of autoimmune diabetes. Seven of these genes affect the activity of regulatory T cells, and two are involved in immune responses to microbial lipids. Four T1D-associated genomic regions could be identified as master regulators via trans-effects on gene expression. These results support the sparse effector hypothesis and reshape our understanding of the genetic architecture of T1D.

Risk factors and prediction of hypoglycaemia using the Hypo-RESOLVE cohort: a secondary analysis of pooled data from insulin clinical trials.

AIMS/HYPOTHESIS: The objective of the Hypoglycaemia REdefining SOLutions for better liVES (Hypo-RESOLVE) project is to use a dataset of pooled clinical trials across pharmaceutical and device companies in people with type 1 or type 2 diabetes to examine factors associated with incident hypoglycaemia events and to quantify the prediction of these events. METHODS: Data from 90 trials with 46,254 participants were pooled. Analyses were done for type 1 and type 2 diabetes separately. Poisson mixed models, adjusted for age, sex, diabetes duration and trial identifier were fitted to assess the association of clinical variables with hypoglycaemia event counts. Tree-based gradient-boosting algorithms (XGBoost) were fitted using training data and their predictive performance in terms of area under the receiver operating characteristic curve (AUC) evaluated on test data. Baseline models including age, sex and diabetes duration were compared with models that further included a score of hypoglycaemia in the first 6 weeks from study entry, and full models that included further clinical variables. The relative predictive importance of each covariate was assessed using XGBoost's importance procedure. Prediction across the entire trial duration for each trial (mean of 34.8 weeks for type 1 diabetes and 25.3 weeks for type 2 diabetes) was assessed. RESULTS: For both type 1 and type 2 diabetes, variables associated with more frequent hypoglycaemia included female sex, white ethnicity, longer diabetes duration, treatment with human as opposed to analogue-only insulin, higher glucose variability, higher score for hypoglycaemia across the 6 week baseline period, lower BP, lower lipid levels and treatment with psychoactive drugs. Prediction of any hypoglycaemia event of any severity was greater than prediction of hypoglycaemia requiring assistance (level 3 hypoglycaemia), for which events were sparser. For prediction of level 1 or worse hypoglycaemia during the whole follow-up period, the AUC was 0.835 (95% CI 0.826, 0.844) in type 1 diabetes and 0.840 (95% CI 0.831, 0.848) in type 2 diabetes. For level 3 hypoglycaemia, the AUC was lower at 0.689 (95% CI 0.667, 0.712) for type 1 diabetes and 0.705 (95% CI 0.662, 0.748) for type 2 diabetes. Compared with the baseline models, almost all the improvement in prediction could be captured by the individual's hypoglycaemia history, glucose variability and blood glucose over a 6 week baseline period. CONCLUSIONS/INTERPRETATION: Although hypoglycaemia rates show large variation according to sociodemographic and clinical characteristics and treatment history, looking at a 6 week period of hypoglycaemia events and glucose measurements predicts future hypoglycaemia risk.

Expression quantitative trait loci analysis in rheumatoid arthritis identifies tissue specific variants associated with severity and outcome.

OBJECTIVE: Genome-wide association studies have successfully identified more than 100 loci associated with susceptibility to rheumatoid arthritis (RA). However, our understanding of the functional effects of genetic variants in causing RA and their effects on disease severity and response to treatment remains limited. METHODS: In this study, we conducted expression quantitative trait locus (eQTL) analysis to dissect the link between genetic variants and gene expression comparing the disease tissue against blood using RNA-Sequencing of synovial biopsies (n=85) and blood samples (n=51) from treatment-naïve patients with RA from the Pathobiology of Early Arthritis Cohort. RESULTS: This identified 898 eQTL genes in synovium and genes loci in blood, with 232 genes in common to both synovium and blood, although notably many eQTL were tissue specific. Examining the HLA region, we uncovered a specific eQTL at HLA-DPB2 with the critical triad of single-nucleotide polymorphisms (SNPs) rs3128921 driving synovial HLA-DPB2 expression, and both rs3128921 and HLA-DPB2 gene expression correlating with clinical severity and increasing probability of the lympho-myeloid pathotype. CONCLUSIONS: This analysis highlights the need to explore functional consequences of genetic associations in disease tissue. HLA-DPB2 SNP rs3128921 could potentially be used to stratify patients to more aggressive treatment immediately at diagnosis.

Flash monitor initiation is associated with improvements in HbA1c levels and DKA rates among people with type 1 diabetes in Scotland: a retrospective nationwide observational study.

AIMS/HYPOTHESIS: We assessed the real-world effect of flash monitor (FM) usage on HbA1c levels and diabetic ketoacidosis (DKA) and severe hospitalised hypoglycaemia (SHH) rates among people with type 1 diabetes in Scotland and across sociodemographic strata within this population. METHODS: This study was retrospective, observational and registry based. Using the national diabetes registry, 14,682 individuals using an FM at any point between 2014 and mid-2020 were identified. Within-person change from baseline in HbA1c following FM initiation was modelled using linear mixed models accounting for within-person pre-exposure trajectory. DKA and SHH events were captured through linkage to hospital admission and mortality data. The difference in DKA and SHH rates between FM-exposed and -unexposed person-time was assessed among users, using generalised linear mixed models with a Poisson likelihood. In a sensitivity analysis, we tested whether changes in these outcomes were seen in an age-, sex- and baseline HbA1c-matched sample of non-users over the same time period. RESULTS: Prevalence of ever-FM use was 45.9% by mid-2020, with large variations by age and socioeconomic status: 64.3% among children aged <13 years vs 32.7% among those aged ≥65 years; and 54.4% vs 36.2% in the least-deprived vs most-deprived quintile. Overall, the median (IQR) within-person change in HbA1c in the year following FM initiation was -2.5 (-9.0, 2.5) mmol/mol (-0.2 [-0.8, 0.2]%). The change varied widely by pre-usage HbA1c: -15.5 (-31.0, -4.0) mmol/mol (-1.4 [-2.8, -0.4]%) in those with HbA1c > 84 mmol/mol [9.8%] and 1.0 (-2.0, 5.5) mmol/mol (0.1 [-0.2, 0.5]%) in those with HbA1c < 54 mmol/mol (7.1%); the corresponding estimated fold change (95% CI) was 0.77 (0.76, 0.78) and 1.08 (1.07, 1.09). Significant reductions in HbA1c were found in all age bands, sexes and socioeconomic strata, and regardless of prior/current pump use, completion of a diabetes education programme or early FM adoption. Variation between the strata of these factors beyond that driven by differing HbA1c at baseline was slight. No change in HbA1c in matched non-users was observed in the same time period (median [IQR] within-person change = 0.5 [-5.0, 5.5] mmol/mol [0.0 (-0.5, 0.5)%]). DKA rates decreased after FM initiation overall and in all strata apart from the adolescents. Estimated overall reduction in DKA event rates (rate ratio) was 0.59 [95% credible interval (CrI) 0.53, 0.64]) after FM vs before FM initiation, accounting for pre-exposure trend. Finally, among those at higher risk for SHH, estimated reduction in event rates was rate ratio 0.25 (95%CrI 0.20, 0.32) after FM vs before FM initiation. CONCLUSIONS/INTERPRETATION: FM initiation is associated with clinically important reductions in HbA1c and striking reduction in DKA rate. Increasing uptake among the socioeconomically disadvantaged offers considerable potential for tightening the current socioeconomic disparities in glycaemia-related outcomes.

Development and validation of a cardiovascular risk prediction model in type 1 diabetes.

AIMS/HYPOTHESIS: We aimed to report current rates of CVD in type 1 diabetes and to develop a CVD risk prediction tool for type 1 diabetes. METHODS: A cohort of 27,527 people with type 1 diabetes without prior CVD was derived from the national register in Scotland. Incident CVD events during 199,552 person-years of follow-up were ascertained using hospital admissions and death registers. A Poisson regression model of CVD was developed and then validated in the Swedish National Diabetes Register (n = 33,183). We compared the percentage with a high 10 year CVD risk (i.e., ≥10%) using the model with the percentage eligible for statins using current guidelines by age. RESULTS: The age-standardised rate of CVD per 100,000 person-years was 4070 and 3429 in men and women, respectively, with type 1 diabetes in Scotland, and 4014 and 3956 in men and women in Sweden. The final model was well calibrated (Hosmer-Lemeshow test p > 0.05) and included a further 22 terms over a base model of age, sex and diabetes duration (C statistic 0.82; 95% CI 0.81, 0.83). The model increased the base model C statistic from 0.66 to 0.80, from 0.60 to 0.75 and from 0.62 to 0.68 in those aged <40, 40-59 and ≥ 60 years, respectively (all p values <0.005). The model required minimal calibration in Sweden and had a C statistic of 0.85. Under current guidelines, >90% of those aged 20-39 years and 100% of those ≥40 years with type 1 diabetes were eligible for statins, but it was not until age 65 upwards that 100% had a modelled risk of CVD ≥10% in 10 years. CONCLUSIONS/INTERPRETATION: A prediction tool such as that developed here can provide individualised risk predictions. This 10 year CVD risk prediction tool could facilitate patient discussions regarding appropriate statin prescribing. Apart from 10 year risk, such discussions may also consider longer-term CVD risk, the potential for greater benefits from early vs later statin intervention, the potential impact on quality of life of an early CVD event and evidence on safety, all of which could influence treatment decisions, particularly in younger people with type 1 diabetes.

The association of polypharmacy and high-risk drug classes with adverse health outcomes in the Scottish population with type 1 diabetes.

AIMS/HYPOTHESIS: The aim of this work was to map the number of prescribed drugs over age, sex and area-based socioeconomic deprivation, and to examine the association between the number of drugs and particular high-risk drug classes with adverse health outcomes among a national cohort of individuals with type 1 diabetes. METHODS: Utilising linked healthcare records from the population-based diabetes register of Scotland, we identified 28,245 individuals with a diagnosis of type 1 diabetes on 1 January 2017. For this population, we obtained information on health status, predominantly reflecting diabetes-related complications, and information on the total number of drugs and particular high-risk drug classes prescribed. We then studied the association of these baseline-level features with hospital admissions for falls, diabetic ketoacidosis (DKA), and hypoglycaemia or death within the subsequent year using multivariate Cox proportional hazards models. RESULTS: Not considering insulin and treatment for hypoglycaemia, the mean number of prescribed drugs was 4.00 (SD 4.35). The proportion of individuals being prescribed five or more drugs at baseline consistently increased with age (proportion [95% CI]: 0-19 years 2.04% [1.60, 2.49]; 40-49 years 28.50% [27.08, 29.93]; 80+ years 76.04% [67.73, 84.84]). Controlling for age, sex, area-based socioeconomic deprivation and health status, each additional drug at baseline was associated with an increase in the hazard for hospitalisation for falls, hypoglycaemia and death but not for DKA admissions (HR [95% CI]: falls 1.03 [1.01, 1.06]; DKA 1.01 [1.00, 1.03]; hypoglycaemia 1.05 [1.02, 1.07]; death 1.04 [1.02, 1.06]). We found a number of drug classes to be associated with an increased hazard of one or more of these adverse health outcomes, including antithrombotic/anticoagulant agents, corticosteroids, opioids, antiepileptics, antipsychotics, hypnotics and sedatives, and antidepressants. CONCLUSIONS: Polypharmacy is common among the Scottish population with type 1 diabetes and is strongly patterned by sociodemographic factors. The number of prescribed drugs and the prescription of particular high-risk drug classes are strong markers of an increased risk of adverse health outcomes, including acute complications of diabetes.

Marked improvements in glycaemic outcomes following insulin pump therapy initiation in people with type 1 diabetes: a nationwide observational study in Scotland.

AIMS/HYPOTHESIS: Our aim was to assess the use of continuous subcutaneous insulin infusion (CSII) in people with type 1 diabetes in Scotland and its association with glycaemic control, as measured by HbA1c levels, frequency of diabetic ketoacidosis (DKA) and severe hospitalised hypoglycaemia (SHH), overall and stratified by baseline HbA1c. METHODS: We included 4684 individuals with type 1 diabetes from the national Scottish register, who commenced CSII between 2004 and 2019. We presented crude within-person differences from baseline HbA1c over time since initiation, crude DKA and SHH event-rates pre-/post-CSII exposure. We then used mixed models to assess the significance of CSII exposure, taking into account: (1) the diffuse nature of the intervention (i.e. structured education often precedes initiation); (2) repeated within-person measurements; and (3) background time-trends occurring pre-intervention. RESULTS: HbA1c decreased after CSII initiation, with a median within-person change of -5.5 mmol/mol (IQR -12.0, 0.0) (-0.5% [IQR -1.1, 0.0]). Within-person changes were most substantial in those with the highest baseline HbA1c, with median -21.0 mmol/mol (-30.0, -11.0) (-1.9% [-2.7, -1.0]) change in those with a baseline >84 mmol/mol (9.8%) within a year of exposure, that was sustained: -19.0 mmol/mol (-27.6, -6.5) (-1.7% [-2.5, -0.6]) at ≥5 years. Statistical significance and magnitude of change were supported by the mixed models results. The crude DKA event-rate was significantly lower in post-CSII person-time compared with pre-CSII person-time: 49.6 events (95% CI 46.3, 53.1) per 1000 person-years vs 67.9 (64.1, 71.9); rate ratio from Bayesian mixed models adjusting for pre-exposure trend: 0.61 (95% credible interval [CrI] 0.47, 0.77; posterior probability of reduction pp = 1.00). The crude overall SHH event-rate in post-CSII vs pre-CSII person-time was also lower: 17.8 events (95% CI 15.8, 19.9) per 1000 person-years post-exposure vs 25.8 (23.5, 28.3) pre-exposure; rate ratio from Bayesian mixed models adjusting for pre-exposure trend: 0.67 (95% CrI 0.45, 1.01; pp = 0.97). CONCLUSIONS/INTERPRETATION: CSII therapy was associated with marked falls in HbA1c especially in those with high baseline HbA1c. CSII was independently associated with reduced DKA and SHH rates. CSII appears to be an effective option for intensive insulin therapy in people with diabetes for improving suboptimal glycaemic control.

Quantitative levels of serum N-glycans in type 1 diabetes and their association with kidney disease.

We investigated associations of quantitative levels of N-glycans with hemoglobin A1c (HbA1c), renal function and renal function decline in type 1 diabetes. We measured 46 total N-glycan peaks (GPs) on 1565 serum samples from the Scottish Diabetes Research Network Type 1 Bioresource Study (SDRNT1BIO) and a pool of healthy donors. Quantitation of absolute abundance of each GP used 2AB-labeled mannose-3 as a standard. We studied cross-sectional associations of GPs and derived measures with HbA1c, albumin/creatinine ratio (ACR) and estimated glomerular filtration rate (eGFR), and prospective associations with incident albuminuria and final eGFR. All GPs were 1.4 to 3.2 times more abundant in SDRTN1BIO than in the healthy samples. Absolute levels of all GPs were slightly higher with higher HbA1c, with strongest associations for triantennary trigalactosylated disialylated, triantennary trigalactosylated trisialylated structures with core or outer arm fucose, and tetraantennary tetragalactosylated trisialylated glycans. Most GPs showed increased abundance with worsening ACR. Lower eGFR was associated with higher absolute GP levels, most significantly with biantennary digalactosylated disialylated glycans with and without bisect, triantennary trigalactosylated trisialylated glycans with and without outer arm fucose, and core fucosylated biantennary monogalactosylated monosialylated glycans. Although several GPs were inversely associated prospectively with final eGFR, cross-validated multivariable models did not improve prediction beyond clinical covariates. Elevated HbA1c is associated with an altered N-glycan profile in type 1 diabetes. Although we could not establish GPs to be prognostic of future renal function decline independently of HbA1c, further studies to evaluate their impact in the pathogenesis of diabetic kidney disease are warranted.

Relation of severe COVID-19 in Scotland to transmission-related factors and risk conditions eligible for shielding support: REACT-SCOT case-control study.

BACKGROUND: Clinically vulnerable individuals have been advised to shield themselves during the COVID-19 epidemic. The objectives of this study were to investigate (1) the rate ratio of severe COVID-19 associated with eligibility for the shielding programme in Scotland across the first and second waves of the epidemic and (2) the relation of severe COVID-19 to transmission-related factors in those in shielding and the general population. METHODS: In a matched case-control design, all 178,578 diagnosed cases of COVID-19 in Scotland from 1 March 2020 to 18 February 2021 were matched for age, sex and primary care practice to 1,744,283 controls from the general population. This dataset (REACT-SCOT) was linked to the list of 212,702 individuals identified as eligible for shielding. Severe COVID-19 was defined as cases that entered critical care or were fatal. Rate ratios were estimated by conditional logistic regression. RESULTS: With those without risk conditions as reference category, the univariate rate ratio for severe COVID-19 was 3.21 (95% CI 3.01 to 3.41) in those with moderate risk conditions and 6.3 (95% CI 5.8 to 6.8) in those eligible for shielding. The highest rate was in solid organ transplant recipients: rate ratio 13.4 (95% CI 9.6 to 18.8). Risk of severe COVID-19 increased with the number of adults but decreased with the number of school-age children in the household. Severe COVID-19 was strongly associated with recent exposure to hospital (defined as 5 to 14 days before presentation date): rate ratio 12.3 (95% CI 11.5 to 13.2) overall. The population attributable risk fraction for recent exposure to hospital peaked at 50% in May 2020 and again at 65% in December 2020. CONCLUSIONS: The effectiveness of shielding vulnerable individuals was limited by the inability to control transmission in hospital and from other adults in the household. Mitigating the impact of the epidemic requires control of nosocomial transmission.

Relation of severe COVID-19 to polypharmacy and prescribing of psychotropic drugs: the REACT-SCOT case-control study.

BACKGROUND: The objective of this study was to investigate the relation of severe COVID-19 to prior drug prescribing. METHODS: Severe cases were defined by entry to critical care or fatal outcome. For this matched case-control study (REACT-SCOT), all 4251 cases of severe COVID-19 in Scotland since the start of the epidemic were matched for age, sex and primary care practice to 36,738 controls from the population register. Records were linked to hospital discharges since June 2015 and dispensed prescriptions issued in primary care during the last 240 days. RESULTS: Severe COVID-19 was strongly associated with the number of non-cardiovascular drug classes dispensed. This association was strongest in those not resident in a care home, in whom the rate ratio (95% CI) associated with dispensing of 12 or more drug classes versus none was 10.8 (8.8, 13.3), and in those without any of the conditions designated as conferring increased risk of COVID-19. Of 17 drug classes postulated at the start of the epidemic to be "medications compromising COVID", all were associated with increased risk of severe COVID-19 and these associations were present in those without any of the designated risk conditions. The fraction of cases in the population attributable to exposure to these drug classes was 38%. The largest effect was for antipsychotic agents: rate ratio 4.18 (3.42, 5.11). Other drug classes with large effects included proton pump inhibitors (rate ratio 2.20 (1.72, 2.83) for = 2 defined daily doses/day), opioids (3.66 (2.68, 5.01) for = 50 mg morphine equivalent/day) and gabapentinoids. These associations persisted after adjusting for covariates and were stronger with recent than with non-recent exposure. CONCLUSIONS: Severe COVID-19 is associated with polypharmacy and with drugs that cause sedation, respiratory depression, or dyskinesia; have anticholinergic effects; or affect the gastrointestinal system. These associations are not easily explained by co-morbidity. Measures to reduce the burden of mortality and morbidity from COVID-19 should include reinforcing existing guidance on reducing overprescribing of these drug classes and limiting inappropriate polypharmacy. REGISTRATION: ENCEPP number EUPAS35558.

Association of cerebral venous thrombosis with recent COVID-19 vaccination: case-crossover study using ascertainment through neuroimaging in Scotland.

BACKGROUND: To investigate the association of primary acute cerebral venous thrombosis (CVT) with COVID-19 vaccination through complete ascertainment of all diagnosed CVT in the population of Scotland. METHODS: Case-crossover study comparing cases of CVT recently exposed to vaccination (1-14 days after vaccination) with cases less recently exposed. Cases in Scotland from 1 December 2020 were ascertained through neuroimaging studies up to 17 May 2021 and diagnostic coding of hospital discharges up to 28 April 2021, linked to national vaccination records. The main outcome measure was primary acute CVT. RESULTS: Of 50 primary acute CVT cases, 29 were ascertained only from neuroimaging studies, 2 were ascertained only from hospital discharges, and 19 were ascertained from both sources. Of these 50 cases, 14 had received the Astra-Zeneca ChAdOx1 vaccine and 3 the Pfizer BNT162b2 vaccine. The incidence of CVT per million doses in the first 14 days after vaccination was 2.2 (95% credible interval 0.9 to 4.1) for ChAdOx1 and 1 (95% credible interval 0.1 to 2.9) for BNT162b2. The rate ratio for CVT associated with exposure to ChAdOx1 in the first 14 days compared with exposure 15-84 days after vaccination was 3.2 (95% credible interval 1.1 to 9.5). CONCLUSIONS: These findings support a causal association between CVT and the AstraZeneca vaccine. The absolute risk of post-vaccination CVT in this population-wide study in Scotland was lower than has been reported for populations in Scandinavia and Germany; the explanation for this is not clear.

Predicting renal disease progression in a large contemporary cohort with type 1 diabetes mellitus.

AIMS/HYPOTHESIS: The aim of this study was to provide data from a contemporary population-representative cohort on rates and predictors of renal decline in type 1 diabetes. METHODS: We used data from a cohort of 5777 people with type 1 diabetes aged 16 and older, diagnosed before the age of 50, and representative of the adult population with type 1 diabetes in Scotland (Scottish Diabetes Research Network Type 1 Bioresource; SDRNT1BIO). We measured serum creatinine and urinary albumin/creatinine ratio (ACR) at recruitment and linked the data to the national electronic healthcare records. RESULTS: Median age was 44.1 years and diabetes duration 20.9 years. The prevalence of CKD stages G1, G2, G3 and G4 and end-stage renal disease (ESRD) was 64.0%, 29.3%, 5.4%, 0.6%, 0.7%, respectively. Micro/macroalbuminuria prevalence was 8.6% and 3.0%, respectively. The incidence rate of ESRD was 2.5 (95% CI 1.9, 3.2) per 1000 person-years. The majority (59%) of those with chronic kidney disease stages G3-G5 did not have albuminuria on the day of recruitment or previously. Over 11.6 years of observation, the median annual decline in eGFR was modest at -1.3 ml min-1 [1.73 m]-2 year-1 (interquartile range [IQR]: -2.2, -0.4). However, 14% experienced a more significant loss of at least 3 ml min-1 [1.73 m]-2. These decliners had more cardiovascular disease (OR 1.9, p = 5 × 10-5) and retinopathy (OR 1.3 p = 0.02). Adding HbA1c, prior cardiovascular disease, recent mean eGFR and prior trajectory of eGFR to a model with age, sex, diabetes duration, current eGFR and ACR maximised the prediction of final eGFR (r2 increment from 0.698 to 0.745, p < 10-16). Attempting to model nonlinearity in eGFR decline or to detect latent classes of decliners did not improve prediction. CONCLUSIONS: These data show much lower levels of kidney disease than historical estimates. However, early identification of those destined to experience significant decline in eGFR remains challenging.

Time trends in deaths before age 50 years in people with type 1 diabetes: a nationwide analysis from Scotland 2004-2017.

AIMS/HYPOTHESIS: We aimed to examine whether crude mortality and mortality relative to the general population below 50 years of age have improved in recent years in those with type 1 diabetes. METHODS: Individuals with type 1 diabetes aged below 50 and at least 1 year old at any time between 2004 and 2017 in Scotland were identified using the national register. Death data were obtained by linkage to Scottish national death registrations. Indirect age standardisation was used to calculate sex-specific standardised mortality ratios (SMRs). Poisson regression was used to test for calendar-time effects as incidence rate ratios (IRRs). RESULTS: There were 1138 deaths in 251,143 person-years among 27,935 people with type 1 diabetes. There was a significant decline in mortality rate over time (IRR for calendar year 0.983 [95% CI 0.967, 0.998], p = 0.03), but the SMR remained approximately stable at 3.1 and 3.6 in men and 4.09 and 4.16 in women for 2004 and 2017, respectively. Diabetic ketoacidosis or coma (DKAoC) accounted for 22% of deaths and the rate did not decline significantly (IRR 0.975 [95% CI 0.94, 1.011], p = 0.168); 79.3% of DKAoC deaths occurred out of hospital. Circulatory diseases accounted for 27% of deaths and did decline significantly (IRR 0.946 [95% CI 0.914, 0.979], p = 0.002). CONCLUSIONS/INTERPRETATION: Absolute mortality has fallen, but the relative impact of type 1 diabetes on mortality below 50 years has not improved. There is scope to improve prevention of premature circulatory diseases and DKAoC and to develop more effective strategies for enabling people with type 1 diabetes to avoid clinically significant hyper- or hypoglycaemia. Graphical abstract.

Comparison of serum and urinary biomarker panels with albumin/creatinine ratio in the prediction of renal function decline in type 1 diabetes.

AIMS/HYPOTHESIS: We examined whether candidate biomarkers in serum or urine can improve the prediction of renal disease progression in type 1 diabetes beyond prior eGFR, comparing their performance with urinary albumin/creatinine ratio (ACR). METHODS: From the population-representative Scottish Diabetes Research Network Type 1 Bioresource (SDRNT1BIO) we sampled 50% and 25% of those with starting eGFR below and above 75 ml min-1 [1.73 m]-2, respectively (N = 1629), and with median 5.1 years of follow-up. Multiplexed ELISAs and single molecule array technology were used to measure nine serum biomarkers and 13 urine biomarkers based on our and others' prior work using large discovery and candidate studies. Associations with final eGFR and with progression to <30 ml min-1 [1.73] m-2, both adjusted for baseline eGFR, were tested using linear and logistic regression models. Parsimonious biomarker panels were identified using a penalised Bayesian approach, and their performance was evaluated through tenfold cross-validation and compared with using urinary ACR and other clinical record data. RESULTS: Seven serum and seven urine biomarkers were strongly associated with either final eGFR or progression to <30 ml min-1 [1.73 m]-2, adjusting for baseline eGFR and other covariates (all at p<2.3 × 10-3). Of these, associations of four serum biomarkers were independent of ACR for both outcomes. The strongest associations with both final eGFR and progression to <30 ml min-1 [1.73 m]-2 were for serum TNF receptor 1, kidney injury molecule 1, CD27 antigen, α-1-microglobulin and syndecan-1. These serum associations were also significant in normoalbuminuric participants for both outcomes. On top of baseline covariates, the r2 for prediction of final eGFR increased from 0.702 to 0.743 for serum biomarkers, and from 0.702 to 0.721 for ACR alone. The area under the receiver operating characteristic curve for progression to <30 ml min-1 [1.73 m]-2 increased from 0.876 to 0.953 for serum biomarkers, and to 0.911 for ACR alone. Other urinary biomarkers did not outperform ACR. CONCLUSIONS/INTERPRETATION: A parsimonious panel of serum biomarkers easily measurable along with serum creatinine may outperform ACR for predicting renal disease progression in type 1 diabetes, potentially obviating the need for urine testing.

Prediction of colorectal cancer risk based on profiling with common genetic variants.

Increasing numbers of common genetic variants associated with colorectal cancer (CRC) have been identified. Our study aimed to determine whether risk prediction based on common genetic variants might enable stratification for CRC risk. Meta-analysis of 11 genome-wide association studies comprising 16 871 cases and 26 328 controls was performed to capture CRC susceptibility variants. Genetic prediction models with several candidate polygenic risk scores (PRSs) were generated from Scottish CRC case-control studies (6478 cases and 11 043 controls) and the score with the best performance was then tested in UK Biobank (UKBB) (4800 cases and 20 287 controls). A weighted PRS of 116 CRC single nucleotide polymorphisms (wPRS116 ) was found with the best predictive performance, reporting a c-statistics of 0.60 and an odds ratio (OR) of 1.46 (95% confidence interval [CI] = 1.41-1.50, per SD increase) in Scottish data set. The predictive performance of this wPRS116 was consistently validated in UKBB data set with c-statistics of 0.61 and an OR of 1.49 (95% CI = 1.44-1.54, per SD increase). Modeling the levels of PRS with age and sex in the general UK population shows that employing genetic risk profiling can achieve a moderate degree of risk discrimination that could be helpful to identify a subpopulation with higher CRC risk due to genetic susceptibility.

Rapid Epidemiological Analysis of Comorbidities and Treatments as risk factors for COVID-19 in Scotland (REACT-SCOT): A population-based case-control study.

BACKGROUND: The objectives of this study were to identify risk factors for severe coronavirus disease 2019 (COVID-19) and to lay the basis for risk stratification based on demographic data and health records. METHODS AND FINDINGS: The design was a matched case-control study. Severe COVID-19 was defined as either a positive nucleic acid test for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the national database followed by entry to a critical care unit or death within 28 days or a death certificate with COVID-19 as underlying cause. Up to 10 controls per case matched for sex, age, and primary care practice were selected from the national population register. For this analysis-based on ascertainment of positive test results up to 6 June 2020, entry to critical care up to 14 June 2020, and deaths registered up to 14 June 2020-there were 36,948 controls and 4,272 cases, of which 1,894 (44%) were care home residents. All diagnostic codes from the past 5 years of hospitalisation records and all drug codes from prescriptions dispensed during the past 240 days were extracted. Rate ratios for severe COVID-19 were estimated by conditional logistic regression. In a logistic regression using the age-sex distribution of the national population, the odds ratios for severe disease were 2.87 for a 10-year increase in age and 1.63 for male sex. In the case-control analysis, the strongest risk factor was residence in a care home, with rate ratio 21.4 (95% CI 19.1-23.9, p = 8 × 10-644). Univariate rate ratios for conditions listed by public health agencies as conferring high risk were 2.75 (95% CI 1.96-3.88, p = 6 × 10-9) for type 1 diabetes, 1.60 (95% CI 1.48-1.74, p = 8 × 10-30) for type 2 diabetes, 1.49 (95% CI 1.37-1.61, p = 3 × 10-21) for ischemic heart disease, 2.23 (95% CI 2.08-2.39, p = 4 × 10-109) for other heart disease, 1.96 (95% CI 1.83-2.10, p = 2 × 10-78) for chronic lower respiratory tract disease, 4.06 (95% CI 3.15-5.23, p = 3 × 10-27) for chronic kidney disease, 5.4 (95% CI 4.9-5.8, p = 1 × 10-354) for neurological disease, 3.61 (95% CI 2.60-5.00, p = 2 × 10-14) for chronic liver disease, and 2.66 (95% CI 1.86-3.79, p = 7 × 10-8) for immune deficiency or suppression. Seventy-eight percent of cases and 52% of controls had at least one listed condition (51% of cases and 11% of controls under age 40). Severe disease was associated with encashment of at least one prescription in the past 9 months and with at least one hospital admission in the past 5 years (rate ratios 3.10 [95% CI 2.59-3.71] and 2.75 [95% CI 2.53-2.99], respectively) even after adjusting for the listed conditions. In those without listed conditions, significant associations with severe disease were seen across many hospital diagnoses and drug categories. Age and sex provided 2.58 bits of information for discrimination. A model based on demographic variables, listed conditions, hospital diagnoses, and prescriptions provided an additional 1.07 bits (C-statistic 0.804). A limitation of this study is that records from primary care were not available. CONCLUSIONS: We have shown that, along with older age and male sex, severe COVID-19 is strongly associated with past medical history across all age groups. Many comorbidities beyond the risk conditions designated by public health agencies contribute to this. A risk classifier that uses all the information available in health records, rather than only a limited set of conditions, will more accurately discriminate between low-risk and high-risk individuals who may require shielding until the epidemic is over.

Investigation of genetically regulated gene expression and response to treatment in rheumatoid arthritis highlights an association between IL18RAP expression and treatment response.

OBJECTIVES: In this study, we sought to investigate whether there was any association between genetically regulated gene expression (as predicted using various reference panels) and anti-tumour necrosis factor (anti-TNF) treatment response (change in erythrocyte sedimentation rate (ESR)) using 3158 European ancestry patients with rheumatoid arthritis. METHODS: The genetically regulated portion of gene expression was estimated in the full cohort of 3158 subjects (as well as within a subcohort consisting of 1575 UK patients) using the PrediXcan software package with three different reference panels. Estimated expression was tested for association with anti-TNF treatment response. As a replication/validation experiment, we also investigated the correlation between change in ESR with measured gene expression at the Interleukin 18 Receptor Accessory Protein (IL18RAP) gene in whole blood and synovial tissue, using an independent replication data set of patients receiving conventional synthetic disease modifying anti-rheumatic drugs, with directly measured (via RNA sequencing) gene expression. RESULTS: We found that predicted expression of IL18RAP showed a consistent signal of association with treatment response across the reference panels. In our independent replication data set, IL18RAP expression in whole blood showed correlation with the change in ESR between baseline and follow-up (r=-0.35, p=0.0091). Change in ESR was also correlated with the expression of IL18RAP in synovial tissue (r=-0.28, p=0.02). CONCLUSION: Our results suggest that IL18RAP expression is worthy of further investigation as a potential predictor of treatment response in rheumatoid arthritis that is not specific to a particular drug type.

Biomarkers associated with early stages of kidney disease in adolescents with type 1 diabetes.

OBJECTIVES: To identify biomarkers of renal disease in adolescents with type 1 diabetes (T1D) and to compare findings in adults with T1D. METHODS: Twenty-five serum biomarkers were measured, using a Luminex platform, in 553 adolescents (median [interquartile range] age: 13.9 [12.6, 15.2] years), recruited to the Adolescent Type 1 Diabetes Cardio-Renal Intervention Trial. Associations with baseline and final estimated glomerular filtration rate (eGFR), rapid decliner and rapid increaser phenotypes (eGFR slopes <-3 and > 3 mL/min/1.73m2 /year, respectively), and albumin-creatinine ratio (ACR) were assessed. Results were also compared with those obtained in 859 adults (age: 55.5 [46.1, 64.4) years) from the Scottish Diabetes Research Network Type 1 Bioresource. RESULTS: In the adolescent cohort, baseline eGFR was negatively associated with trefoil factor-3, cystatin C, and beta-2 microglobulin (B2M) (B coefficient[95%CI]: -0.19 [-0.27, -0.12], P = 7.0 × 10-7 ; -0.18 [-0.26, -0.11], P = 5.1 × 10-6 ; -0.12 [-0.20, -0.05], P = 1.6 × 10-3 ), in addition to clinical covariates. Final eGFR was negatively associated with osteopontin (-0.21 [-0.28, -0.14], P = 2.3 × 10-8 ) and cystatin C (-0.16 [-0.22, -0.09], P = 1.6 × 10-6 ). Rapid decliner phenotype was associated with osteopontin (OR: 1.83 [1.42, 2.41], P = 7.3 × 10-6 ), whereas rapid increaser phenotype was associated with fibroblast growth factor-23 (FGF-23) (1.59 [1.23, 2.04], P = 2.6 × 10-4 ). ACR was not associated with any of the biomarkers. In the adult cohort similar associations with eGFR were found; however, several additional biomarkers were associated with eGFR and ACR. CONCLUSIONS: In this young population with T1D and high rates of hyperfiltration, osteopontin was the most consistent biomarker associated with prospective changes in eGFR. FGF-23 was associated with eGFR increases, whereas trefoil factor-3, cystatin C, and B2M were associated with baseline eGFR.

Genome-Wide Association Study of Diabetic Kidney Disease Highlights Biology Involved in Glomerular Basement Membrane Collagen.

BACKGROUND: Although diabetic kidney disease demonstrates both familial clustering and single nucleotide polymorphism heritability, the specific genetic factors influencing risk remain largely unknown. METHODS: To identify genetic variants predisposing to diabetic kidney disease, we performed genome-wide association study (GWAS) analyses. Through collaboration with the Diabetes Nephropathy Collaborative Research Initiative, we assembled a large collection of type 1 diabetes cohorts with harmonized diabetic kidney disease phenotypes. We used a spectrum of ten diabetic kidney disease definitions based on albuminuria and renal function. RESULTS: Our GWAS meta-analysis included association results for up to 19,406 individuals of European descent with type 1 diabetes. We identified 16 genome-wide significant risk loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain (COL4A3) gene, which encodes a major structural component of the glomerular basement membrane (GBM). Mutations in COL4A3 are implicated in heritable nephropathies, including the progressive inherited nephropathy Alport syndrome. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of diabetic kidney disease, including albuminuria and ESKD, and demonstrated a significant association with GBM width; protective allele carriers had thinner GBM before any signs of kidney disease, and its effect was dependent on glycemia. Three other loci are in or near genes with known or suggestive involvement in this condition (BMP7) or renal biology (COLEC11 and DDR1). CONCLUSIONS: The 16 diabetic kidney disease-associated loci may provide novel insights into the pathogenesis of this condition and help identify potential biologic targets for prevention and treatment.

Biomarker panels associated with progression of renal disease in type 1 diabetes.

AIMS/HYPOTHESIS: We aimed to identify a sparse panel of biomarkers for improving the prediction of renal disease progression in type 1 diabetes. METHODS: We considered 859 individuals recruited from the Scottish Diabetes Research Network Type 1 Bioresource (SDRNT1BIO) and 315 individuals from the Finnish Diabetic Nephropathy (FinnDiane) study. All had an entry eGFR between 30 and 75 ml min-1[1.73 m]-2, with those from FinnDiane being oversampled for albuminuria. A total of 297 circulating biomarkers (30 proteins, 121 metabolites, 146 tryptic peptides) were measured in non-fasting serum samples using the Luminex platform and LC electrospray tandem MS (LC-MS/MS). We investigated associations with final eGFR adjusted for baseline eGFR and with rapid progression (a loss of more than 3 ml min-1[1.73 m]-2 year-1) using linear and logistic regression models. Panels of biomarkers were identified using a penalised Bayesian approach, and their performance was evaluated through 10-fold cross-validation and compared with using clinical record data alone. RESULTS: For final eGFR, 16 proteins and 30 metabolites or tryptic peptides showed significant association in SDRNT1BIO, and nine proteins and five metabolites or tryptic peptides in FinnDiane, beyond age, sex, diabetes duration, study day eGFR and length of follow-up (all at p < 10-4). The strongest associations were with CD27 antigen (CD27), kidney injury molecule 1 (KIM-1) and α1-microglobulin. Including the Luminex biomarkers on top of baseline covariates increased the r2 for prediction of final eGFR from 0.47 to 0.58 in SDRNT1BIO and from 0.33 to 0.48 in FinnDiane. At least 75% of the increment in r2 was attributable to CD27 and KIM-1. However, using the weighted average of historical eGFR gave similar performance to biomarkers. The LC-MS/MS platform performed less well. CONCLUSIONS/INTERPRETATION: Among a large set of associated biomarkers, a sparse panel of just CD27 and KIM-1 contains most of the predictive information for eGFR progression. The increment in prediction beyond clinical data was modest but potentially useful for oversampling individuals with rapid disease progression into clinical trials, especially where there is little information on prior eGFR trajectories.