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.

Upregulation of Anti-Angiogenic miR-106b-3p Correlates Negatively with IGF-1 and Vascular Health Parameters in a Model of Subclinical Cardiovascular Disease: Study with Metformin Therapy

Well-controlled type 1 diabetes mellitus (T1DM) is regarded as a model of subclinical cardiovascular disease (CVD), characterized by inflammation and adverse vascular health. However, the underlying mechanisms are not fully understood. We investigated insulin-like growth factor-1 (IGF-1) and IGF-binding protein-3 (IGFBP-3) levels, their correlation to miR-106b-3p expression in a subclinical CVD model, and the cardioprotective effect of metformin. A total of 20 controls and 29 well-controlled T1DM subjects were studied. Plasma IGF-1, IGFBP-3 levels, and miR-106b-3p expression in colony-forming unit-Hills were analyzed and compared with vascular markers. miR-106b-3p was upregulated in T1DM (p < 0.05) and negatively correlated with pro-angiogenic markers CD34+/100-lymphocytes (p < 0.05) and IGF-1 (p < 0.05). IGF-1 was downregulated in T1DM (p < 0.01), which was associated with increased inflammatory markers TNF-α, CRP, and IL-10 and reduced CD34+/100-lymphocytes. IGFBP-3 had no significant results. Metformin had no effect on IGF-1 but significantly reduced miR-106b-3p (p < 0.0001). An Ingenuity Pathway analysis predicted miR-106b-3p to inhibit PDGFA, PIK3CG, GDNF, and ADAMTS13, which activated CVD. Metformin was predicted to be cardioprotective by inhibiting miR-106b-3p. In conclusion: Subclinical CVD is characterized by a cardio-adverse profile of low IGF-1 and upregulated miR-106b-3p. We demonstrated that the cardioprotective effect of metformin may be via downregulation of upregulated miR-106b-3p and its effect on downstream targets.

A brief comparison of polygenic risk scores and Mendelian randomisation

Mendelian randomisation and polygenic risk score analysis have become increasingly popular in the last decade due to the advent of large-scale genome-wide association studies. Each approach has valuable applications, some of which are overlapping, yet there are important differences which we describe here.

Chronic inflammation does not mediate the effect of adiposity on grip strength: results from a multivariable Mendelian randomization study

The relationship between adiposity and grip strength (GS) is complex. We investigated whether one pathway through which adiposity affects GS was via chronic inflammation. 367,583 UK Biobank participants had body mass index (BMI), waist-hip-ratio (WHR), C-reactive protein (CRP) and GS data. Univariable Mendelian randomization (MR) and multivariable Mendelian randomization (MVMR) analyses (using inverse variance weighted (IVW) weighted median estimates (WME) and MR-Egger models) estimated total, direct and indirect effects of adiposity traits on GS using genetic instruments for BMI and WHR (exposures) and CRP (mediator). Observational findings suggested higher BMI was associated with stronger grip, e.g., in males, per standard deviation (SD) higher BMI, GS was higher by 0.48 kg (95% confidence interval(CI):0.44,0.51), independent of CRP. For males MR estimates were directionally consistent; for females, estimates were consistent with the null. Observational findings for WHR suggested that higher WHR was associated with weaker grip. In multivariable MR-IVW analyses, effects in males were consistent with the null. In females, there were consistent effects such that higher WHR was associated with stronger grip, e.g., 1-SD higher WHR was associated with 1.25 kg (MVMR-Egger; 95% CI:0.72,1.78) stronger grip, independent of CRP. Across sexes and adiposity indicators, CRP's mediating role was minor. Greater adiposity may increase GS in early old age, but effects vary by sex and adiposity location. There was no evidence that inflammation mediated these effects.

Decoding of miR-7-5p in Colony Forming Unit-Hill Colonies as a Biomarker of Subclinical Cardiovascular Disease-A MERIT Study

Colony forming unit-Hill (CFU-Hill) colonies were established to serve as a sensitive biomarker for vascular health. In animals, the overexpression of miR-7-5p was shown to be pro-atherogenic and associated with increased cardiovascular disease (CVD) risk. In a MERIT study, we aimed to explore the role of miR-7-5p expression in CFU-Hill colonies in type 1 diabetes mellitus (T1DM) and the effect of metformin in subclinical CVD. The expression of miR-7-5p in CFU-Hill colonies in 29 T1DM subjects without CVD and 20 healthy controls (HC) was measured. Metformin was administered to T1DM subjects for eight weeks. MiR-7-5p was upregulated in T1DM whereas metformin reduced it to HC levels. MiR-7-5p was positively correlated with c-reactive protein, and C-X-C motif chemokine ligand 10. The receiver operating characteristic curve revealed miR-7-5p as a biomarker of CVD, and upregulated miR-7-5p, defining subclinical CVD at a HbA1c level of 44.3 mmol/mol. Ingenuity pathway analysis predicted miR-7-5p to inhibit the mRNA expression of Krüppel-like factor 4, epidermal growth factor receptor, insulin-like growth factor 1 receptor, v-raf-1 murine leukemia viral oncogene homolog 1 and insulin receptor substrate ½, and insulin receptor, while metformin activated these miRNAs via transforming growth factor-β1 and Smad2/3. We proved the pro-atherogenic effect of miR-7-5p that maybe used as a prognostic biomarker.

Altered microvasculature in pancreatic islets from subjects with type 1 diabetes

AIMS

The transcriptome of different dissociated pancreatic islet cells has been described in enzymatically isolated islets in both health and disease. However, the isolation, culturing, and dissociation procedures likely affect the transcriptome profiles, distorting the biological conclusions. The aim of the current study was to characterize the cells of the islets of Langerhans from subjects with and without type 1 diabetes in a way that reflects the in vivo situation to the highest possible extent.

METHODS

Islets were excised using laser capture microdissection directly from frozen pancreatic tissue sections obtained from organ donors with (n = 7) and without (n = 8) type 1 diabetes. Transcriptome analysis of excised samples was performed using AmpliSeq. Consecutive pancreatic sections were used to estimate the proportion of beta-, alpha-, and delta cells using immunofluorescence and to examine the presence of CD31 positive endothelial regions using immunohistochemistry.

RESULTS

The proportion of beta cells in islets from subjects with type 1 diabetes was reduced to 0% according to both the histological and transcriptome data, and several alterations in the transcriptome were derived from the loss of beta cells. In total, 473 differentially expressed genes were found in the islets from subjects with type 1 diabetes. Functional enrichment analysis showed that several of the most upregulated gene sets were related to vasculature and angiogenesis, and histologically, vascular density was increased in subjects with type 1 diabetes. Downregulated in type 1 diabetes islets was the gene set epithelial mesenchymal transition.

CONCLUSION

A number of transcriptional alterations are present in islets from subjects with type 1 diabetes. In particular, several gene sets related to vasculature and angiogenesis are upregulated and there is an increased vascular density, suggesting an altered microvasculature in islets from subjects with type 1 diabetes. By studying pancreatic islets extracted directly from snap-frozen pancreatic tissue, this study reflects the in vivo situation to a high degree and gives important insights into islet pathophysiology in type 1 diabetes.

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.

Integrating polygenic risk scores in the prediction of type 2 diabetes risk and subtypes in British Pakistanis and Bangladeshis: A population-based cohort study

BACKGROUND

Type 2 diabetes (T2D) is highly prevalent in British South Asians, yet they are underrepresented in research. Genes & Health (G&H) is a large, population study of British Pakistanis and Bangladeshis (BPB) comprising genomic and routine health data. We assessed the extent to which genetic risk for T2D is shared between BPB and European populations (EUR). We then investigated whether the integration of a polygenic risk score (PRS) for T2D with an existing risk tool (QDiabetes) could improve prediction of incident disease and the characterisation of disease subtypes.

METHODS AND FINDINGS

In this observational cohort study, we assessed whether common genetic loci associated with T2D in EUR individuals were replicated in 22,490 BPB individuals in G&H. We replicated fewer loci in G&H (n = 76/338, 22%) than would be expected given power if all EUR-ascertained loci were transferable (n = 101, 30%; p = 0.001). Of the 27 transferable loci that were powered to interrogate this, only 9 showed evidence of shared causal variants. We constructed a T2D PRS and combined it with a clinical risk instrument (QDiabetes) in a novel, integrated risk tool (IRT) to assess risk of incident diabetes. To assess model performance, we compared categorical net reclassification index (NRI) versus QDiabetes alone. In 13,648 patients free from T2D followed up for 10 years, NRI was 3.2% for IRT versus QDiabetes (95% confidence interval (CI): 2.0% to 4.4%). IRT performed best in reclassification of individuals aged less than 40 years deemed low risk by QDiabetes alone (NRI 5.6%, 95% CI 3.6% to 7.6%), who tended to be free from comorbidities and slim. After adjustment for QDiabetes score, PRS was independently associated with progression to T2D after gestational diabetes (hazard ratio (HR) per SD of PRS 1.23, 95% CI 1.05 to 1.42, p = 0.028). Using cluster analysis of clinical features at diabetes diagnosis, we replicated previously reported disease subgroups, including Mild Age-Related, Mild Obesity-related, and Insulin-Resistant Diabetes, and showed that PRS distribution differs between subgroups (p = 0.002). Integrating PRS in this cluster analysis revealed a Probable Severe Insulin Deficient Diabetes (pSIDD) subgroup, despite the absence of clinical measures of insulin secretion or resistance. We also observed differences in rates of progression to micro- and macrovascular complications between subgroups after adjustment for confounders. Study limitations include the absence of an external replication cohort and the potential biases arising from missing or incorrect routine health data.

CONCLUSIONS

Our analysis of the transferability of T2D loci between EUR and BPB indicates the need for larger, multiancestry studies to better characterise the genetic contribution to disease and its varied aetiology. We show that a T2D PRS optimised for this high-risk BPB population has potential clinical application in BPB, improving the identification of T2D risk (especially in the young) on top of an established clinical risk algorithm and aiding identification of subgroups at diagnosis, which may help future efforts to stratify care and treatment of the disease.

Arginine-vasopressin mediates counter-regulatory glucagon release and is diminished in type 1 diabetes

Insulin-induced hypoglycemia is a major treatment barrier in type-1 diabetes (T1D). Accordingly, it is important that we understand the mechanisms regulating the circulating levels of glucagon. Varying glucose over the range of concentrations that occur physiologically between the fed and fuel-deprived states (8 to 4 mM) has no significant effect on glucagon secretion in the perfused mouse pancreas or in isolated mouse islets (in vitro), and yet associates with dramatic increases in plasma glucagon. The identity of the systemic factor(s) that elevates circulating glucagon remains unknown. Here, we show that arginine-vasopressin (AVP), secreted from the posterior pituitary, stimulates glucagon secretion. Alpha-cells express high levels of the vasopressin 1b receptor (V1bR) gene (Avpr1b). Activation of AVP neurons in vivo increased circulating copeptin (the C-terminal segment of the AVP precursor peptide) and increased blood glucose; effects blocked by pharmacological antagonism of either the glucagon receptor or V1bR. AVP also mediates the stimulatory effects of hypoglycemia produced by exogenous insulin and 2-deoxy-D-glucose on glucagon secretion. We show that the A1/C1 neurons of the medulla oblongata drive AVP neuron activation in response to insulin-induced hypoglycemia. AVP injection increased cytoplasmic Ca2+ in alpha-cells (implanted into the anterior chamber of the eye) and glucagon release. Hypoglycemia also increases circulating levels of AVP/copeptin in humans and this hormone stimulates glucagon secretion from human islets. In patients with T1D, hypoglycemia failed to increase both copeptin and glucagon. These findings suggest that AVP is a physiological systemic regulator of glucagon secretion and that this mechanism becomes impaired in T1D.

Exploiting collider bias to apply two-sample summary data Mendelian randomization methods to one-sample individual level data

Over the last decade the availability of SNP-trait associations from genome-wide association studies has led to an array of methods for performing Mendelian randomization studies using only summary statistics. A common feature of these methods, besides their intuitive simplicity, is the ability to combine data from several sources, incorporate multiple variants and account for biases due to weak instruments and pleiotropy. With the advent of large and accessible fully-genotyped cohorts such as UK Biobank, there is now increasing interest in understanding how best to apply these well developed summary data methods to individual level data, and to explore the use of more sophisticated causal methods allowing for non-linearity and effect modification. In this paper we describe a general procedure for optimally applying any two sample summary data method using one sample data. Our procedure first performs a meta-analysis of summary data estimates that are intentionally contaminated by collider bias between the genetic instruments and unmeasured confounders, due to conditioning on the observed exposure. These estimates are then used to correct the standard observational association between an exposure and outcome. Simulations are conducted to demonstrate the method's performance against naive applications of two sample summary data MR. We apply the approach to the UK Biobank cohort to investigate the causal role of sleep disturbance on HbA1c levels, an important determinant of diabetes. Our approach can be viewed as a generalization of Dudbridge et al. (Nat. Comm. 10: 1561), who developed a technique to adjust for index event bias when uncovering genetic predictors of disease progression based on case-only data. Our work serves to clarify that in any one sample MR analysis, it can be advantageous to estimate causal relationships by artificially inducing and then correcting for collider bias.

Zinc-alpha2-glycoprotein, dysglycaemia and insulin resistance: a systematic review and meta-analysis

To systematically review the current literature investigating associations between zinc-alpha2-glycoprotein (ZAG) and dysglycaemia (including type 2 diabetes (T2DM), poly-cystic-ovary syndrome (PCOS), pre-diabetes or insulin resistance). This included relationships between ZAG and continuous measures of insulin and glucose. Additionally, we performed a meta-analysis to estimate the extent that ZAG differs between individuals with or without dysglycaemia; whilst examining the potential influence of adiposity. A systematic search was performed on four databases for studies on circulating ZAG concentrations in adult human populations, comparing healthy controls to individuals with dysglycaemia. Key characteristics, including the mean ZAG concentrations (mg∙L^-1), and any correlational statistics between ZAG and continuous measures of glucose, glycated haemoglobin (HbA1c) or insulin were extracted. Meta-analyses were performed to compare metabolically healthy controls to cases, and on studies that compared controls and cases considered overweight or obese (body mass index (BMI) ≥25 kg^.m²). 1575 papers were identified and 14 studies (16 cohorts) were considered eligible for inclusion. Circulating ZAG was lower in individuals with dysglycaemia compared to metabolically healthy controls (-4.14 [-8.17, -0.11] mg^.L^-1; I² = 98.5%; p < 0.001). When using data from only studies with overweight or obese groups with or without dysglycaemia (three studies (four cohorts); pooled n = 332), the difference in circulating ZAG was no longer significant (-0.30 [-3.67, 3.07] mg^. L^-1; I² = 28.0%; p = 0.225). These data suggest that ZAG may be implicated in dysglycaemia, although there was significant heterogeneity across different studies and the mediating effect of adiposity cannot be excluded. Therefore, more research is needed before robust conclusions can be drawn.

GABAA Receptor-Mediated Currents and Hormone mRNAs in Cells Expressing More Than One Hormone Transcript in Intact Human Pancreatic Islets

In pancreatic islets, the major cell-types are α, β and δ cells. The γ-aminobutyric acid (GABA) signalling system is expressed in human pancreatic islets. In single hormone transcript-expressing cells, we have previously characterized the functional properties of islet GABAA receptors (iGABAARs). Here, we extended these studies to islet cells expressing mRNAs for more than one hormone and sought for correlation between iGABAAR activity level and relative mRNA expression ratio. The single-cell RT-PCR in combination with the patch-clamp current recordings was used to examine functional properties of iGABAARs in the multiple hormone mRNA-expressing cells. We detected cells expressing double (α/β, α/δ, β/δ cell-types) and triple (α/β/δ cell-type) hormone transcripts. The most common mixed-identity cell-type was the α/β group where the cells could be grouped into β- and α-like subgroups. The β-like cells had low GCG/INS expression ratio (<0.6) and significantly higher frequency of iGABAAR single-channel openings than the α-like cells where the GCG/INS expression ratio was high (>1.2). The hormone expression levels and iGABAAR single-channel characteristics varied in the α/β/δ cell-type. Clearly, multiple hormone transcripts can be expressed in islet cells whereas iGABAAR single-channel functional properties appear to be α or β cell specific.

Genome-Wide Association Analyses in 128,266 Individuals Identifies New Morningness and Sleep Duration Loci

Disrupted circadian rhythms and reduced sleep duration are associated with several human diseases, particularly obesity and type 2 diabetes, but until recently, little was known about the genetic factors influencing these heritable traits. We performed genome-wide association studies of self-reported chronotype (morning/evening person) and self-reported sleep duration in 128,266 white British individuals from the UK Biobank study. Sixteen variants were associated with chronotype (P<5x10^-8), including variants near the known circadian rhythm genes RGS16 (1.21 odds of morningness, 95% CI [1.15, 1.27], P = 3x10^-12) and PER2 (1.09 odds of morningness, 95% CI [1.06, 1.12], P = 4x10^-10). The PER2 signal has previously been associated with iris function. We sought replication using self-reported data from 89,283 23andMe participants; thirteen of the chronotype signals remained associated at P<5x10^-8 on meta-analysis and eleven of these reached P<0.05 in the same direction in the 23andMe study. We also replicated 9 additional variants identified when the 23andMe study was used as a discovery GWAS of chronotype (all P<0.05 and meta-analysis P<5x10^-8). For sleep duration, we replicated one known signal in PAX8 (2.6 minutes per allele, 95% CI [1.9, 3.2], P = 5.7x10^-16) and identified and replicated two novel associations at VRK2 (2.0 minutes per allele, 95% CI [1.3, 2.7], P = 1.2x10^-9; and 1.6 minutes per allele, 95% CI [1.1, 2.2], P = 7.6x10^-9). Although we found genetic correlation between chronotype and BMI (rG = 0.056, P = 0.05); undersleeping and BMI (rG = 0.147, P = 1x10^-5) and oversleeping and BMI (rG = 0.097, P = 0.04), Mendelian Randomisation analyses, with limited power, provided no consistent evidence of causal associations between BMI or type 2 diabetes and chronotype or sleep duration. Our study brings the total number of loci associated with chronotype to 22 and with sleep duration to three, and provides new insights into the biology of sleep and circadian rhythms in humans.

Numerous studies have identified links between too little or too much sleep and circadian misalignment with metabolic disorders such as obesity and type 2 diabetes. However, cause-and-effect is not easily determined, because of multiple confounding factors affecting both sleep patterns and disease risk. Using the first release of the UK Biobank study, which combines detailed measurements and questionnaire data with genetic data, we investigate the genetics of two self-report sleep measures, chronotype and average sleep duration, in 128,266 white British individuals. We replicate previous genetic associations and identify seven and two novel genetic variants influencing chronotype and sleep duration, respectively. Associated variants are located near genes implicated in circadian rhythm regulation (RGS16, PER2), near a serotonin receptor gene (HTR6) and another gene (INADL) encoding a protein thought to be important in photosensitive retinal cells, cells known to communicate with the brain’s primary circadian pacemaker. Using the genetic risk factors, we estimate the unconfounded causal associations of BMI and type 2 diabetes on sleep patterns (and vice versa) through Mendelian Randomisation. However, we find no evidence for causal associations in either direction. The full UK Biobank release of 500,000 individuals will boost our power to detect causal associations.