Control of inducible gene expression links cohesin to hematopoietic progenitor self-renewal and differentiation.

Cohesin is important for 3D genome organization. Nevertheless, even the complete removal of cohesin has surprisingly little impact on steady-state gene transcription and enhancer activity. Here we show that cohesin is required for the core transcriptional response of primary macrophages to microbial signals, and for inducible enhancer activity that underpins inflammatory gene expression. Consistent with a role for inflammatory signals in promoting myeloid differentiation of hematopoietic stem and progenitor cells (HPSCs), cohesin mutations in HSPCs led to reduced inflammatory gene expression and increased resistance to differentiation-inducing inflammatory stimuli. These findings uncover an unexpected dependence of inducible gene expression on cohesin, link cohesin with myeloid differentiation, and may help explain the prevalence of cohesin mutations in human acute myeloid leukemia.

Genome sequencing reveals loci under artificial selection that underlie disease phenotypes in the laboratory rat.

Large numbers of inbred laboratory rat strains have been developed for a range of complex disease phenotypes. To gain insights into the evolutionary pressures underlying selection for these phenotypes, we sequenced the genomes of 27 rat strains, including 11 models of hypertension, diabetes, and insulin resistance, along with their respective control strains. Altogether, we identified more than 13 million single-nucleotide variants, indels, and structural variants across these rat strains. Analysis of strain-specific selective sweeps and gene clusters implicated genes and pathways involved in cation transport, angiotensin production, and regulators of oxidative stress in the development of cardiovascular disease phenotypes in rats. Many of the rat loci that we identified overlap with previously mapped loci for related traits in humans, indicating the presence of shared pathways underlying these phenotypes in rats and humans. These data represent a step change in resources available for evolutionary analysis of complex traits in disease models.

PPMS: A framework to Profile Primary MicroRNAs from Single-cell RNA-sequencing datasets.

MOTIVATION: Single-cell/nuclei RNA-sequencing (scRNA-seq) technologies can simultaneously quantify gene expression in thousands of cells across the genome. However, the majority of the noncoding RNAs, such as microRNAs (miRNAs), cannot currently be profiled at the same scale. MiRNAs are a class of small noncoding RNAs and play an important role in gene regulation. MiRNAs originate from the processing of primary transcripts, known as primary-microRNAs (pri-miRNAs). The pri-miRNA transcripts, independent of their cognate miRNAs, can also function as long noncoding RNAs, code for micropeptides or even interact with DNA, acting like enhancers. Therefore, it is apparent that the significance of scRNA-seq pri-miRNA profiling expands beyond using pri-miRNA as proxies of mature miRNAs. However, there are no computational methods that allow profiling and quantification of pri-miRNAs at the single-cell-type resolution. RESULTS: We have developed a simple yet effective computational framework to profile pri-MiRNAs from single-cell RNA-sequencing datasets (PPMS). Based on user input, PPMS can profile pri-miRNAs at cell-type resolution. PPMS can be applied to both newly produced and publicly available datasets obtained via single cell or single-nuclei RNA-seq. It allows users to (i) investigate the distribution of pri-miRNAs across cell types and cell states and (ii) establish a relationship between the number of cells/reads sequenced and the detection of pri-miRNAs. Here, to demonstrate its efficacy, we have applied PPMS to publicly available scRNA-seq data generated from (i) individual chambers (ventricles and atria) of the human heart, (ii) human pluripotent stem cells during their differentiation into cardiomyocytes (the heart beating cells) and (iii) hiPSCs-derived cardiomyocytes infected with severe acute respiratory syndrome coronavirus 2.

Photoreceptor laminin drives differentiation of human pluripotent stem cells to photoreceptor progenitors that partially restore retina function.

Blindness caused by advanced stages of inherited retinal diseases and age-related macular degeneration are characterized by photoreceptor loss. Cell therapy involving replacement with functional photoreceptor-like cells generated from human pluripotent stem cells holds great promise. Here, we generated a human recombinant retina-specific laminin isoform, LN523, and demonstrated the role in promoting the differentiation of human embryonic stem cells into photoreceptor progenitors. This chemically defined and xenogen-free method enables reproducible production of photoreceptor progenitors within 32 days. We observed that the transplantation into rd10 mice were able to protect the host photoreceptor outer nuclear layer (ONL) up to 2 weeks post transplantation as measured by full-field electroretinogram. At 4 weeks post transplantation, the engrafted cells were found to survive, mature, and associate with the host's rod bipolar cells. Visual behavioral assessment using the water maze swimming test demonstrated visual improvement in the cell-transplanted rodents. At 20 weeks post transplantation, the maturing engrafted cells were able to replace the loss of host ONL by extensive association with host bipolar cells and synapses. Post-transplanted rabbit model also provided congruent evidence for synaptic connectivity with the degenerated host retina. The results may pave the way for the development of stem cell-based therapeutics for retina degeneration.

Comparative Analysis of Single-Cell RNA Sequencing Methods with and without Sample Multiplexing.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful technique for investigating biological heterogeneity at the single-cell level in human systems and model organisms. Recent advances in scRNA-seq have enabled the pooling of cells from multiple samples into single libraries, thereby increasing sample throughput while reducing technical batch effects, library preparation time, and the overall cost. However, a comparative analysis of scRNA-seq methods with and without sample multiplexing is lacking. In this study, we benchmarked methods from two representative platforms: Parse Biosciences (Parse; with sample multiplexing) and 10x Genomics (10x; without sample multiplexing). By using peripheral blood mononuclear cells (PBMCs) obtained from two healthy individuals, we demonstrate that demultiplexed scRNA-seq data obtained from Parse showed similar cell type frequencies compared to 10x data where samples were not multiplexed. Despite relatively lower cell capture affecting library preparation, Parse can detect rare cell types (e.g., plasmablasts and dendritic cells) which is likely due to its relatively higher sensitivity in gene detection. Moreover, a comparative analysis of transcript quantification between the two platforms revealed platform-specific distributions of gene length and GC content. These results offer guidance for researchers in designing high-throughput scRNA-seq studies.

IL-11 is a crucial determinant of cardiovascular fibrosis.

Fibrosis is a common pathology in cardiovascular disease. In the heart, fibrosis causes mechanical and electrical dysfunction and in the kidney, it predicts the onset of renal failure. Transforming growth factor ß1 (TGFß1) is the principal pro-fibrotic factor, but its inhibition is associated with side effects due to its pleiotropic roles. We hypothesized that downstream effectors of TGFß1 in fibroblasts could be attractive therapeutic targets and lack upstream toxicity. Here we show, using integrated imaging-genomics analyses of primary human fibroblasts, that upregulation of interleukin-11 (IL-11) is the dominant transcriptional response to TGFß1 exposure and required for its pro-fibrotic effect. IL-11 and its receptor (IL11RA) are expressed specifically in fibroblasts, in which they drive non-canonical, ERK-dependent autocrine signalling that is required for fibrogenic protein synthesis. In mice, fibroblast-specific Il11 transgene expression or Il-11 injection causes heart and kidney fibrosis and organ failure, whereas genetic deletion of Il11ra1 protects against disease. Therefore, inhibition of IL-11 prevents fibroblast activation across organs and species in response to a range of important pro-fibrotic stimuli. These results reveal a central role of IL-11 in fibrosis and we propose that inhibition of IL-11 is a potential therapeutic strategy to treat fibrotic diseases.

Intrinsic mutant HTT-mediated defects in oligodendroglia cause myelination deficits and behavioral abnormalities in Huntington disease.

White matter abnormalities are a nearly universal pathological feature of neurodegenerative disorders including Huntington disease (HD). A long-held assumption is that this white matter pathology is simply a secondary outcome of the progressive neuronal loss that manifests with advancing disease. Using a mouse model of HD, here we show that white matter and myelination abnormalities are an early disease feature appearing before the manifestation of any behavioral abnormalities or neuronal loss. We further show that selective inactivation of mutant huntingtin (mHTT) in the NG2+ oligodendrocyte progenitor cell population prevented myelin abnormalities and certain behavioral deficits in HD mice. Strikingly, the improvements in behavioral outcomes were seen despite the continued expression of mHTT in nonoligodendroglial cells including neurons, astrocytes, and microglia. Using RNA-seq and ChIP-seq analyses, we implicate a pathogenic mechanism that involves enhancement of polycomb repressive complex 2 (PRC2) activity by mHTT in the intrinsic oligodendroglial dysfunction and myelination deficits observed in HD. Our findings challenge the long-held dogma regarding the etiology of white matter pathology in HD and highlight the contribution of epigenetic mechanisms to the observed intrinsic oligodendroglial dysfunction. Our results further suggest that ameliorating white matter pathology and oligodendroglial dysfunction may be beneficial for HD.

Brugada syndrome genetics is associated with phenotype severity.

AIMS: Brugada syndrome (BrS) is associated with an increased risk of sudden cardiac death due to ventricular tachycardia/fibrillation (VT/VF) in young, otherwise healthy individuals. Despite SCN5A being the most commonly known mutated gene to date, the genotype-phenotype relationship is poorly understood and remains uncertain. This study aimed to elucidate the genotype-phenotype correlation in BrS. METHODS AND RESULTS: Brugada syndrome probands deemed at high risk of future arrhythmic events underwent genetic testing and phenotype characterization by the means of epicardial arrhythmogenic substrate (AS) mapping, and were divided into two groups according to the presence or absence of SCN5A mutation. Two-hundred probands (160 males, 80%; mean age 42.6 ± 12.2 years) were included in this study. Patients harbouring SCN5A mutations exhibited a spontaneous type 1 pattern and experienced aborted cardiac arrest or spontaneous VT/VF more frequently than the other subjects. SCN5A-positive patients exhibited a larger epicardial AS area, more prolonged electrograms and more frequently observed non-invasive late potentials. The presence of an SCN5A mutation explained >26% of the variation in the epicardial AS area and was the strongest predictor of a large epicardial area. CONCLUSION: In BrS, the genetic background is the main determinant for the extent of the electrophysiological abnormalities. SCN5A mutation carriers exhibit more pronounced epicardial electrical abnormalities and a more aggressive clinical presentation. These results contribute to the understanding of the genetic determinants of the BrS phenotypic expression and provide possible explanations for the varying degrees of disease expression.

Vascular control of kidney epithelial transporters.

A major pathway in hypertension pathogenesis involves direct activation of ANG II type 1 (AT1) receptors in the kidney, stimulating Na+ reabsorption. AT1 receptors in tubular epithelia control expression and stimulation of Na+ transporters and channels. Recently, we found reduced blood pressure and enhanced natriuresis in mice with cell-specific deletion of AT1 receptors in smooth muscle (SMKO mice). Although impaired vasoconstriction and preserved renal blood flow might contribute to exaggerated urinary Na+ excretion in SMKO mice, we considered whether alterations in Na+ transporter expression might also play a role; therefore, we carried out proteomic analysis of key Na+ transporters and associated proteins. Here, we show that levels of Na+-K+-2Cl- cotransporter isoform 2 (NKCC2) and Na+/H+ exchanger isoform 3 (NHE3) are reduced at baseline in SMKO mice, accompanied by attenuated natriuretic and diuretic responses to furosemide. During ANG II hypertension, we found widespread remodeling of transporter expression in wild-type mice with significant increases in the levels of total NaCl cotransporter, phosphorylated NaCl cotransporter (Ser71), and phosphorylated NKCC2, along with the cleaved, activated forms of the α- and γ-epithelial Na+ channel. However, the increases in α- and γ-epithelial Na+ channel with ANG II were substantially attenuated in SMKO mice. This was accompanied by a reduced natriuretic response to amiloride. Thus, enhanced urinary Na+ excretion observed after cell-specific deletion of AT1 receptors from smooth muscle cells is associated with altered Na+ transporter abundance across epithelia in multiple nephron segments. These findings suggest a system of vascular-epithelial in the kidney, modulating the expression of Na+ transporters and contributing to the regulation of pressure natriuresis.NEW & NOTEWORTHY The use of drugs to block the renin-angiotensin system to reduce blood pressure is common. However, the precise mechanism for how these medications control blood pressure is incompletely understood. Here, we show that mice lacking angiotensin receptors specifically in smooth muscle cells lead to alternation in tubular transporter amount and function. Thus, demonstrating the importance of vascular-tubular cross talk in the control of blood pressure.

Cross-species integration of single-cell RNA-seq resolved alveolar-epithelial transitional states in idiopathic pulmonary fibrosis.

Single-cell transcriptomics analyses of the fibrotic lung uncovered two cell states critical to lung injury recovery in the alveolar epithelium-a reparative transitional cell state in the mouse and a disease-specific cell state (KRT5-/KRT17+) in human idiopathic pulmonary fibrosis (IPF). The murine transitional cell state lies between the differentiation from type 2 (AT2) to type 1 pneumocyte (AT1), and the human KRT5-/KRT17+ cell state may arise from the dysregulation of this differentiation process. We review major findings of single-cell transcriptomics analyses of the fibrotic lung and reanalyzed data from seven single-cell RNA sequencing studies of human and murine models of IPF, focusing on the alveolar epithelium. Our comparative and cross-species single-cell transcriptomics analyses allowed us to further delineate the differentiation trajectories from AT2 to AT1 and AT2 to the KRT5-/KRT17+ cell state. We observed AT1 cells in human IPF retain the transcriptional signature of the murine transitional cell state. Using pseudotime analysis, we recapitulated the differentiation trajectories from AT2 to AT1 and from AT2 to KRT5-/KRT17+ cell state in multiple human IPF studies. We further delineated transcriptional programs underlying cell-state transitions and determined the molecular phenotypes at terminal differentiation. We hypothesize that in addition to the reactivation of developmental programs (SOX4, SOX9), senescence (TP63, SOX4) and the Notch pathway (HES1) are predicted to steer intermediate progenitors to the KRT5-/KRT17+ cell state. Our analyses suggest that activation of SMAD3 later in the differentiation process may explain the fibrotic molecular phenotype typical of KRT5-/KRT17+ cells.

Broad regulation of gene isoform expression by Wnt signaling in cancer.

Differential gene isoform expression is a ubiquitous mechanism to enhance proteome diversity and maintain cell homeostasis. Mechanisms such as splicing that drive gene isoform variability are highly dynamic and responsive to changes in cell signaling pathways. Wnt/ß-catenin signaling has profound effects on cell activity and cell fate and is known to modify several splicing events by altering the expression of individual splicing factors. However, a global assessment of how extensively Wnt signaling regulates splicing and other mechanisms that determine mRNA isoform composition in cancer is lacking. We used deep time-resolved RNA-seq in two independent in vivo Wnt-addicted tumor models during treatment with the potent Wnt inhibitor ETC-159 and examined Wnt regulated splicing events and splicing regulators. We found 1025 genes that underwent Wnt regulated variable exon usage leading to isoform expression changes. This was accompanied by extensive Wnt regulated changes in the expression of splicing regulators. Many of these Wnt regulated events were conserved in multiple human cancers, and many were linked to previously defined cancer-associated splicing quantitative trait loci. This suggests that the Wnt regulated splicing events are components of fundamental oncogenic processes. These findings demonstrate the wide-ranging effects of Wnt signaling on the isoform composition of the cell and provides an extensive resource of expression changes of splicing regulators and gene isoforms regulated by Wnt signaling.

Whole-Genome Sequencing of Finnish Type 1 Diabetic Siblings Discordant for Kidney Disease Reveals DNA Variants associated with Diabetic Nephropathy.

BACKGROUND: Several genetic susceptibility loci associated with diabetic nephropathy have been documented, but no causative variants implying novel pathogenetic mechanisms have been elucidated. METHODS: We carried out whole-genome sequencing of a discovery cohort of Finnish siblings with type 1 diabetes who were discordant for the presence (case) or absence (control) of diabetic nephropathy. Controls had diabetes without complications for 15-37 years. We analyzed and annotated variants at genome, gene, and single-nucleotide variant levels. We then replicated the associated variants, genes, and regions in a replication cohort from the Finnish Diabetic Nephropathy study that included 3531 unrelated Finns with type 1 diabetes. RESULTS: We observed protein-altering variants and an enrichment of variants in regions associated with the presence or absence of diabetic nephropathy. The replication cohort confirmed variants in both regulatory and protein-coding regions. We also observed that diabetic nephropathy-associated variants, when clustered at the gene level, are enriched in a core protein-interaction network representing proteins essential for podocyte function. These genes include protein kinases (protein kinase C isoforms ε and ι) and protein tyrosine kinase 2. CONCLUSIONS: Our comprehensive analysis of a diabetic nephropathy cohort of siblings with type 1 diabetes who were discordant for kidney disease points to variants and genes that are potentially causative or protective for diabetic nephropathy. This includes variants in two isoforms of the protein kinase C family not previously linked to diabetic nephropathy, adding support to previous hypotheses that the protein kinase C family members play a role in diabetic nephropathy and might be attractive therapeutic targets.

Common signalling pathways in macrophage and osteoclast multinucleation.

Macrophage cell fusion and multinucleation are fundamental processes in the formation of multinucleated giant cells (MGCs) in chronic inflammatory disease and osteoclasts in the regulation of bone mass. However, this basic cell phenomenon is poorly understood despite its pathophysiological relevance. Granulomas containing multinucleated giant cells are seen in a wide variety of complex inflammatory disorders, as well as in infectious diseases. Dysregulation of osteoclastic bone resorption underlies the pathogenesis of osteoporosis and malignant osteolytic bone disease. Recent reports have shown that the formation of multinucleated giant cells and osteoclast fusion display a common molecular signature, suggesting shared genetic determinants. In this Review, we describe the background of cell-cell fusion and the similar origin of macrophages and osteoclasts. We specifically focus on the common pathways involved in osteoclast and MGC fusion. We also highlight potential approaches that could help to unravel the core mechanisms underlying bone and granulomatous disorders in humans.

Genome-wide analysis of differential RNA editing in epilepsy.

The recoding of genetic information through RNA editing contributes to proteomic diversity, but the extent and significance of RNA editing in disease is poorly understood. In particular, few studies have investigated the relationship between RNA editing and disease at a genome-wide level. Here, we developed a framework for the genome-wide detection of RNA sites that are differentially edited in disease. Using RNA-sequencing data from 100 hippocampi from mice with epilepsy (pilocarpine-temporal lobe epilepsy model) and 100 healthy control hippocampi, we identified 256 RNA sites (overlapping with 87 genes) that were significantly differentially edited between epileptic cases and controls. The degree of differential RNA editing in epileptic mice correlated with frequency of seizures, and the set of genes differentially RNA-edited between case and control mice were enriched for functional terms highly relevant to epilepsy, including "neuron projection" and "seizures." Genes with differential RNA editing were preferentially enriched for genes with a genetic association to epilepsy. Indeed, we found that they are significantly enriched for genes that harbor nonsynonymous de novo mutations in patients with epileptic encephalopathy and for common susceptibility variants associated with generalized epilepsy. These analyses reveal a functional convergence between genes that are differentially RNA-edited in acquired symptomatic epilepsy and those that contribute risk for genetic epilepsy. Taken together, our results suggest a potential role for RNA editing in the epileptic hippocampus in the occurrence and severity of epileptic seizures.

Genome-wide meta-analysis associates GPSM1 with type 2 diabetes, a plausible gene involved in skeletal muscle function.

Genome-wide association studies (GWASs) have identified many genetic variations associated with type 2 diabetes mellitus (T2DM) in Asians, but understanding the functional genetic variants that influence traits is often a complex process. In this study, fine mapping and other analytical strategies were performed to investigate the effects of G protein signaling modulator 1 (GPSM1) on insulin resistance in skeletal muscle. A total of 128 single-nucleotide polymorphisms (SNPs) within GPSM1 were analysed in 21,897 T2DM cases and 32,710 healthy controls from seven GWASs. The SNP rs28539249 in intron 9 of GPSM1 showed a nominally significant association with T2DM in Asians (OR = 1.07, 95% CI = 1.04-1.10, P < 10-4). The GPSM1 mRNA was increased in skeletal muscle and correlated with T2DM traits across obese mice model. An eQTL for the cis-acting regulation of GPSM1 expression in human skeletal muscle was identified for rs28539249, and the increased GPSM1 expression related with T2DM traits within GEO datasets. Another independent Asian cohort showed that rs28539249 is associated with the skeletal muscle expression of CACFD1, GTF3C5, SARDH, and FAM163B genes, which are functionally enriched for endoplasmic reticulum stress (ERS) and unfolded protein response (UPR) pathways. Moreover, rs28539249 locus was predicted to disrupt regulatory regions in human skeletal muscle with enriched epigenetic marks and binding affinity for CTCF. Supershift EMSA assays followed luciferase assays demonstrated the CTCF specifically binding to rs28539249-C allele leading to decreased transcriptional activity. Thus, the post-GWAS annotation confirmed the Asian-specific association of genetic variant in GPSM1 with T2DM, suggesting a role for the variant in the regulation in skeletal muscle.

Decreased microglial Wnt/ß-catenin signalling drives microglial pro-inflammatory activation in the developing brain.

Microglia of the developing brain have unique functional properties but how their activation states are regulated is poorly understood. Inflammatory activation of microglia in the still-developing brain of preterm-born infants is associated with permanent neurological sequelae in 9 million infants every year. Investigating the regulators of microglial activation in the developing brain across models of neuroinflammation-mediated injury (mouse, zebrafish) and primary human and mouse microglia we found using analysis of genes and proteins that a reduction in Wnt/ß-catenin signalling is necessary and sufficient to drive a microglial phenotype causing hypomyelination. We validated in a cohort of preterm-born infants that genomic variation in the Wnt pathway is associated with the levels of connectivity found in their brains. Using a Wnt agonist delivered by a blood-brain barrier penetrant microglia-specific targeting nanocarrier we prevented in our animal model the pro-inflammatory microglial activation, white matter injury and behavioural deficits. Collectively, these data validate that the Wnt pathway regulates microglial activation, is critical in the evolution of an important form of human brain injury and is a viable therapeutic target.

MWASTools: an R/bioconductor package for metabolome-wide association studies.

Summary: MWASTools is an R package designed to provide an integrated pipeline to analyse metabonomic data in large-scale epidemiological studies. Key functionalities of our package include: quality control analysis; metabolome-wide association analysis using various models (partial correlations, generalized linear models); visualization of statistical outcomes; metabolite assignment using statistical total correlation spectroscopy (STOCSY); and biological interpretation of metabolome-wide association studies results. Availability and implementation: The MWASTools R package is implemented in R (version > =3.4) and is available from Bioconductor: https://bioconductor.org/packages/MWASTools/. Contact: m.dumas@imperial.ac.uk. Supplementary information: Supplementary data are available at Bioinformatics online.

Systems biology of the IMIDIA biobank from organ donors and pancreatectomised patients defines a novel transcriptomic signature of islets from individuals with type 2 diabetes.

AIMS/HYPOTHESIS: Pancreatic islet beta cell failure causes type 2 diabetes in humans. To identify transcriptomic changes in type 2 diabetic islets, the Innovative Medicines Initiative for Diabetes: Improving beta-cell function and identification of diagnostic biomarkers for treatment monitoring in Diabetes (IMIDIA) consortium ( www.imidia.org ) established a comprehensive, unique multicentre biobank of human islets and pancreas tissues from organ donors and metabolically phenotyped pancreatectomised patients (PPP). METHODS: Affymetrix microarrays were used to assess the islet transcriptome of islets isolated either by enzymatic digestion from 103 organ donors (OD), including 84 non-diabetic and 19 type 2 diabetic individuals, or by laser capture microdissection (LCM) from surgical specimens of 103 PPP, including 32 non-diabetic, 36 with type 2 diabetes, 15 with impaired glucose tolerance (IGT) and 20 with recent-onset diabetes (<1 year), conceivably secondary to the pancreatic disorder leading to surgery (type 3c diabetes). Bioinformatics tools were used to (1) compare the islet transcriptome of type 2 diabetic vs non-diabetic OD and PPP as well as vs IGT and type 3c diabetes within the PPP group; and (2) identify transcription factors driving gene co-expression modules correlated with insulin secretion ex vivo and glucose tolerance in vivo. Selected genes of interest were validated for their expression and function in beta cells. RESULTS: Comparative transcriptomic analysis identified 19 genes differentially expressed (false discovery rate ≤0.05, fold change ≥1.5) in type 2 diabetic vs non-diabetic islets from OD and PPP. Nine out of these 19 dysregulated genes were not previously reported to be dysregulated in type 2 diabetic islets. Signature genes included TMEM37, which inhibited Ca2+-influx and insulin secretion in beta cells, and ARG2 and PPP1R1A, which promoted insulin secretion. Systems biology approaches identified HNF1A, PDX1 and REST as drivers of gene co-expression modules correlated with impaired insulin secretion or glucose tolerance, and 14 out of 19 differentially expressed type 2 diabetic islet signature genes were enriched in these modules. None of these signature genes was significantly dysregulated in islets of PPP with impaired glucose tolerance or type 3c diabetes. CONCLUSIONS/INTERPRETATION: These studies enabled the stringent definition of a novel transcriptomic signature of type 2 diabetic islets, regardless of islet source and isolation procedure. Lack of this signature in islets from PPP with IGT or type 3c diabetes indicates differences possibly due to peculiarities of these hyperglycaemic conditions and/or a role for duration and severity of hyperglycaemia. Alternatively, these transcriptomic changes capture, but may not precede, beta cell failure.

Changes in macrophage transcriptome associate with systemic sclerosis and mediate GSDMA contribution to disease risk.

OBJECTIVES: Several common and rare risk variants have been reported for systemic sclerosis (SSc), but the effector cell(s) mediating the function of these genetic variants remains to be elucidated. While innate immune cells have been proposed as the critical targets to interfere with the disease process underlying SSc, no studies have comprehensively established their effector role. Here we investigated the contribution of monocyte-derived macrophages (MDMs) in mediating genetic susceptibility to SSc. METHODS: We carried out RNA sequencing and genome-wide genotyping in MDMs from 57 patients with SSc and 15 controls. Our differential expression and expression quantitative trait locus (eQTL) analysis in SSc was further integrated with epigenetic, expression and eQTL data from skin, monocytes, neutrophils and lymphocytes. RESULTS: We identified 602 genes upregulated and downregulated in SSc macrophages that were significantly enriched for genes previously implicated in SSc susceptibility (P=5×10-4), and 270 cis-regulated genes in MDMs. Among these, GSDMA was reported to carry an SSc risk variant (rs3894194) regulating expression of neighbouring genes in blood. We show that GSDMA is upregulated in SSc MDMs (P=8.4×10-4) but not in the skin, and is a significant eQTL in SSc macrophages and lipopolysaccharide/interferon gamma (IFNγ)-stimulated monocytes. Furthermore, we identify an SSc macrophage transcriptome signature characterised by upregulation of glycolysis, hypoxia and mTOR signalling and a downregulation of IFNγ response pathways. CONCLUSIONS: Our data further establish the link between macrophages and SSc, and suggest that the contribution of the rs3894194 risk variant to SSc susceptibility can be mediated by GSDMA expression in macrophages.