Single-cell transcriptomic profiling of human pancreatic islets reveals genes responsive to glucose exposure over 24 h.

AIMS/HYPOTHESIS: Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycaemia, beta cell glucotoxicity and subsequently type 2 diabetes. In this study, we explored the effects of in vitro hyperglycaemic conditions on human pancreatic islet gene expression across 24 h in six pancreatic cell types: alpha; beta; gamma; delta; ductal; and acinar. We hypothesised that genes associated with hyperglycaemic conditions may be relevant to the onset and progression of diabetes. METHODS: We exposed human pancreatic islets from two donors to low (2.8 mmol/l) and high (15.0 mmol/l) glucose concentrations over 24 h in vitro. To assess the transcriptome, we performed single-cell RNA-seq (scRNA-seq) at seven time points. We modelled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Additionally, we integrated genomic features and genetic summary statistics to nominate candidate effector genes. For three of these genes, we functionally characterised the effect on insulin production and secretion using CRISPR interference to knock down gene expression in EndoC-ßH1 cells, followed by a glucose-stimulated insulin secretion assay. RESULTS: In the discrete time models, we identified 1344 genes associated with time and 668 genes associated with glucose exposure across all cell types and time points. In the continuous time models, we identified 1311 genes associated with time, 345 genes associated with glucose exposure and 418 genes associated with interaction effects between time and glucose across all cell types. By integrating these expression profiles with summary statistics from genetic association studies, we identified 2449 candidate effector genes for type 2 diabetes, HbA1c, random blood glucose and fasting blood glucose. Of these candidate effector genes, we showed that three (ERO1B, HNRNPA2B1 and RHOBTB3) exhibited an effect on glucose-stimulated insulin production and secretion in EndoC-ßH1 cells. CONCLUSIONS/INTERPRETATION: The findings of our study provide an in-depth characterisation of the 24 h transcriptomic response of human pancreatic islets to glucose exposure at a single-cell resolution. By integrating differentially expressed genes with genetic signals for type 2 diabetes and glucose-related traits, we provide insights into the molecular mechanisms underlying glucose homeostasis. Finally, we provide functional evidence to support the role of three candidate effector genes in insulin secretion and production. DATA AVAILABILITY: The scRNA-seq data from the 24 h glucose exposure experiment performed in this study are available in the database of Genotypes and Phenotypes (dbGap; https://www.ncbi.nlm.nih.gov/gap/ ) with accession no. phs001188.v3.p1. Study metadata and summary statistics for the differential expression, gene set enrichment and candidate effector gene prediction analyses are available in the Zenodo data repository ( https://zenodo.org/ ) under accession number 11123248. The code used in this study is publicly available at https://github.com/CollinsLabBioComp/publication-islet_glucose_timecourse .

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic human embryonic stem cell-derived ß-like cells.

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional roles of many loci remain unexplored. Here, we engineered isogenic knockout human embryonic stem cell lines for 20 genes associated with T2D risk. We examined the impacts of each knockout on ß cell differentiation, functions, and survival. We generated gene expression and chromatin accessibility profiles on ß cells derived from each knockout line. Analyses of T2D-association signals overlapping HNF4A-dependent ATAC peaks identified a likely causal variant at the FAIM2 T2D-association signal. Additionally, the integrative association analyses identified four genes (CP, RNASE1, PCSK1N, and GSTA2) associated with insulin production, and two genes (TAGLN3 and DHRS2) associated with ß cell sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental line and identified a single likely functional variant at each of 23 T2D-association signals.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic hESC-derived ß-like cells.

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional role of many loci has remained unexplored. In this study, we engineered isogenic knockout human embryonic stem cell (hESC) lines for 20 genes associated with T2D risk. We systematically examined ß-cell differentiation, insulin production and secretion, and survival. We performed RNA-seq and ATAC-seq on hESC-ß cells from each knockout line. Analyses of T2D GWAS signals overlapping with HNF4A-dependent ATAC peaks identified a specific SNP as a likely causal variant. In addition, we performed integrative association analyses and identified four genes ( CP, RNASE1, PCSK1N and GSTA2 ) associated with insulin production, and two genes ( TAGLN3 and DHRS2 ) associated with sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental hESC line, to identify a single likely functional variant at each of 23 T2D GWAS signals.

A multi-organoid platform identifies CIART as a key factor for SARS-CoV-2 infection.

COVID-19 is a systemic disease involving multiple organs. We previously established a platform to derive organoids and cells from human pluripotent stem cells to model SARS-CoV-2 infection and perform drug screens1,2. This provided insight into cellular tropism and the host response, yet the molecular mechanisms regulating SARS-CoV-2 infection remain poorly defined. Here we systematically examined changes in transcript profiles caused by SARS-CoV-2 infection at different multiplicities of infection for lung airway organoids, lung alveolar organoids and cardiomyocytes, and identified several genes that are generally implicated in controlling SARS-CoV-2 infection, including CIART, the circadian-associated repressor of transcription. Lung airway organoids, lung alveolar organoids and cardiomyocytes derived from isogenic CIART-/- human pluripotent stem cells were significantly resistant to SARS-CoV-2 infection, independently of viral entry. Single-cell RNA-sequencing analysis further validated the decreased levels of SARS-CoV-2 infection in ciliated-like cells of lung airway organoids. CUT&RUN, ATAC-seq and RNA-sequencing analyses showed that CIART controls SARS-CoV-2 infection at least in part through the regulation of NR4A1, a gene also identified from the multi-organoid analysis. Finally, transcriptional profiling and pharmacological inhibition led to the discovery that the Retinoid X Receptor pathway regulates SARS-CoV-2 infection downstream of CIART and NR4A1. The multi-organoid platform identified the role of circadian-clock regulation in SARS-CoV-2 infection, which provides potential therapeutic targets for protection against COVID-19 across organ systems.

Ugi Reaction Mediated Detergent Assembly for Membrane Protein Studies.

Despite the continuous efforts, the current repertoire of detergents is still far from sufficient for the biophysics studies of membrane proteins (MPs). Toward the rapid expansion of detergent diversity, we herein report a new strategy based on Ugi reaction mediated modular assembly. Structural varieties, including hydrophobic tails and hydrophilic heads, could be conveniently introduced from the multiple reaction components. New detergents then were comprehensively evaluated in the physical properties and preliminarily screened by the thermal stabilization for a transporter MsbA and a spectrum of G protein-coupled receptors (GPCRs). For the glucagon-like peptide-1 receptor (GLP-1R), a class B GPCR, detergent M-23-M finally stood out in a second evaluation for the maintenance of homogeneity and was further illustrated its application in the improvement of NMR study. Besides the promising utility in the MP study, the current results exhibit intriguing structural-physical relationship that would allow the guidance in the tuning of detergent properties in the future.

Two-Dimensional Detergent Expansion Strategy for Membrane Protein Studies.

Detergents are the most frequently applied reagents in membrane protein (MP) studies. The limited diversity of one-head-one-tailed traditional detergents, however, is far from sufficient for structurally distinct MPs. Expansion of detergent repertoire has a continuous momentum. In line with the speculation that detergent pre-assembly exerts superiority, herein we report for the first time cross-conjugation of two series of monomeric detergents for constructing a two-dimensional library of dimeric detergents. Optimum detergents stood out with unique preferences in the systematic evaluation of individual MPs. Furthermore, unprecedented hybrid detergents 14M8G and 14M9G enabled high-quality EM study of transporter MsbA and NMR study of G protein-coupled receptor A2A AR, respectively. Given the abundance of cross-coupling chemistries, comprehensive diversity could be readily covered that would facilitate the finding of new detergents for the manipulation of thorny MPs and innovation of the functional and structural study in future.

Structure-Activity Relationship Studies of Hydantoin-Cored Ligands for Smoothened Receptor.

An underside binding site was recently identified in the transmembrane domain of smoothened receptor (SMO). Herein, we report efforts in the exploration of new insights into the interactions between the ligand and SMO. The hydantoin core in the middle of the parent compound was found to be highly conservative in chirality, ring size, and substituents. On each benzene at two ends, a plethora of variations, particularly halogen substitutions, were introduced and investigated. Analysis of the structure-activity relationship revealed miscellaneous halogen effects. The ligands with double halogen substituents exhibit remarkably enhanced potency, providing promising candidates that potentially overcome the common drug resistance and useful heavy-atom labeled chemical tools for co-crystallization studies of SMO.

Elucidation of Distinct Modular Assemblies of Smoothened Receptor by Bitopic Ligand Measurement.

Class F G protein-coupled receptors are characterized by a large extracellular domain (ECD) in addition to the common transmembrane domain (TMD) with seven α-helixes. For smoothened receptor (SMO), structural studies revealed dissected ECD and TMD, and their integrated assemblies. However, distinct assemblies were reported under different circumstances. Using an unbiased approach based on four series of cross-conjugated bitopic ligands, we explore the relationship between the active status and receptor assembly. Different activity dependency on the linker length for these bitopic ligands corroborates the various occurrences of SMO assembly. These results reveal a rigid "near" assembly for active SMO, which is in contrast to previous results. Conversely, inactive SMO adopts a free ECD, which would be remotely captured at "far" assembly by cholesterol. Altogether, we propose a mechanism of cholesterol flow-caused SMO activation involving an erection of ECD from far to near assembly.

Rational Remodeling of Atypical Scaffolds for the Design of Photoswitchable Cannabinoid Receptor Tools.

Azobenzene-embedded photoswitchable ligands are the widely used chemical tools in photopharmacological studies. Current approaches to azobenzene introduction rely mainly on the isosteric replacement of typical azologable groups. However, atypical scaffolds may offer more opportunities for photoswitch remodeling, which are chemically in an overwhelming majority. Herein, we investigate the rational remodeling of atypical scaffolds for azobenzene introduction, as exemplified in the development of photoswitchable ligands for the cannabinoid receptor 2 (CB2). Based on the analysis of residue-type clusters surrounding the binding pocket, we conclude that among the three representative atypical arms of the CB2 antagonist, AM10257, the adamantyl arm is the most appropriate for azobenzene remodeling. The optimizing spacer length and attachment position revealed AzoLig 9 with excellent thermal bistability, decent photopharmacological switchability between its two configurations, and high subtype selectivity. This structure-guided approach gave new impetus in the extension of new chemical spaces for tool customization for increasingly diversified photo-pharmacological studies and beyond.

Catalytically Cleavable Detergent for Membrane Protein Studies.

Throughout the in vitro studies of membrane proteins (MPs), proper detergents are essential for the preparation of stable aqueous samples. To date, universally applicable detergents have not yet been reported to accommodate the distinct requirements for the highly diversified MPs and at the different stages of MP manipulation. Detergent exchange often has to be performed. We report herein the catalytically cleavable detergents (CatCDs) that can be efficiently removed to facilitate a complete exchange. To this end, functional groups, like propargyl and allyl, are introduced as branched chains or built in the hydrophobic chain close to the hydrophilic head. The representative CatCDs can be used as usual detergents in the extraction and purification of MPs and later be removed upon the addition of catalytic palladium. Mediated by CatCD-1, reconstitution of a transporter protein MsbA into a series of detergents was achieved. The extension of these designs could facilitate the future optimization of other biophysics studies.

SARS-CoV-2 infection induces beta cell transdifferentiation.

Recent clinical data have suggested a correlation between coronavirus disease 2019 (COVID-19) and diabetes. Here, we describe the detection of SARS-CoV-2 viral antigen in pancreatic beta cells in autopsy samples from individuals with COVID-19. Single-cell RNA sequencing and immunostaining from ex vivo infections confirmed that multiple types of pancreatic islet cells were susceptible to SARS-CoV-2, eliciting a cellular stress response and the induction of chemokines. Upon SARS-CoV-2 infection, beta cells showed a lower expression of insulin and a higher expression of alpha and acinar cell markers, including glucagon and trypsin1, respectively, suggesting cellular transdifferentiation. Trajectory analysis indicated that SARS-CoV-2 induced eIF2-pathway-mediated beta cell transdifferentiation, a phenotype that could be reversed with trans-integrated stress response inhibitor (trans-ISRIB). Altogether, this study demonstrates an example of SARS-CoV-2 infection causing cell fate change, which provides further insight into the pathomechanisms of COVID-19.

Disulfide-Containing Detergents (DCDs) for the Structural Biology of Membrane Proteins.

Disulfide-containing detergents (DCDs) are introduced, which contain a disulfide bond in the hydrophobic tail. DCDs form smaller micelles than corresponding detergents with linear hydrocarbon chains, while providing good solubilization and reconstitution of membrane proteins. The use of this new class of detergents in structural biology is illustrated with solution NMR spectra of the human G protein-coupled receptor A2A AR, which is an α-helical protein, and the ß-barrel protein OmpX from E. coli.

The structure-based traceless specific fluorescence labeling of the smoothened receptor.

The smoothened receptor (SMO) mediates the hedgehog (Hh) signaling pathway and plays a vital role in embryonic development and tumorigenesis. The visualization of SMO has the potential to provide new insights into its enigmatic mechanisms and associated disease pathogenesis. Based on recent progress in structural studies of SMO, we have designed and characterized a group of affinity probes to facilitate the turn-on fluorescence labeling of SMO at the ε-amine of K395. These chemical probes were derived from a potent SMO antagonist skeleton by the conjugation of a small non-fluorescent unit, O-nitrobenzoxadiazole (O-NBD). In this context, optimal probes were developed to be capable of efficiently and selectively lighting up SMO regardless of whether it is in micelles or in native membranes. More importantly, the resulting labeled SMO only bears a very small fluorophore and allows for the recovery of the unoccupied pocket by dissociation of the residual ligand module. These advantages should allow the probe to serve as a potential tool for monitoring SMO trafficking, understanding Hh activation mechanisms, and even the diagnosis of tumorigenesis in the future.

A Chemical Strategy for Amphiphile Replacement in Membrane Protein Research.

Membrane mimics are indispensable tools in the structural and functional understanding of membrane proteins (MPs). Given stringent requirements of integral MP manipulations, amphiphile replacement is often required in sample preparation for various biophysical purposes. Current protocols generally rely on physical methodologies and rarely reach complete replacement. In comparison, we report herein a chemical alternative that facilitates the exhaustive exchange of membrane-mimicking systems for MP reconstitution. This method, named sacrifice-replacement strategy, was enabled by a class of chemically cleavable detergents (CCDs), derived from the disulfide incorporation in the traditional detergent n-dodecyl-ß-d-maltopyranoside. The representative CCD behaved well in both solubilizing the diverse α-helical human G protein-coupled receptors and refolding of the ß-barrel bacterial outer membrane protein X, and more importantly, it could also be readily degraded under mild conditions. By this means, the A2A adenosine receptor was successfully reconstituted into a series of commercial detergents for stabilization screening and nanodiscs for electron microscopy analysis. Featured by the simplicity and compatibility, this CCD-mediated strategy would later find more applications when being integrated in other biophysics studies.

A structurally guided dissection-then-evolution strategy for ligand optimization of smoothened receptor.

We present herein a novel dissection-then-evolution strategy for ligand optimization. Using the co-crystal structure of the smoothened receptor (SMO) as a guide, we studied the modular contribution of LY2940680 by systematically "silencing" the specific interaction between the individual residue(s) and the fragment in the ligand. Following evolution by focusing on the benzoyl part finally yielded an improved ligand 21.