Single-cell Transcriptome Profiling Research Articles

The stress response regulator HSF1 modulates natural killer cell anti-tumour immunity.

Diverse cellular insults converge on activation of the heat shock factor 1 (HSF1), which regulates the proteotoxic stress response to maintain protein homoeostasis. HSF1 regulates numerous gene programmes beyond the proteotoxic stress response in a cell-type- and context-specific manner to promote malignancy. However, the role(s) of HSF1 in immune populations of the tumour microenvironment remain elusive. Here, we leverage an in vivo model of HSF1 activation and single-cell transcriptomic tumour profiling to show that augmented HSF1 activity in natural killer (NK) cells impairs cytotoxicity, cytokine production and subsequent anti-tumour immunity. Mechanistically, HSF1 directly binds and regulates the expression of key mediators of NK cell effector function. This work demonstrates that HSF1 regulates the immune response under the stress conditions of the tumour microenvironment. These findings have important implications for enhancing the efficacy of adoptive NK cell therapies and for designing combinatorial strategies including modulators of NK cell-mediated tumour killing.

Decoding the pathogenesis of spermatogenic failure in cryptorchidism through single-cell transcriptomic profiling

Cryptorchidism, commonly known as undescended testis, affects 1%-9% of male newborns, posing infertility and testis tumor risks. Despite its prevalence, the detailed pathophysiology underlying male infertility within cryptorchidism remains unclear. Here, we profile and analyze 46,644 single-cell transcriptomes from individual testicular cells obtained from adult males diagnosed with cryptorchidism and healthy controls. Spermatogenesis compromise in cryptorchidism links primarily to spermatogonium self-renewal and differentiation dysfunctions. We illuminate the involvement of testicular somatic cells, including immune cells, thereby unveiling the activation and degranulation of mast cells in cryptorchidism. Mast cells are identified as contributors to interstitial fibrosis via transforming growth factor β1 (TGF-β1) and cathepsin G secretion. Furthermore, significantly increased levels of secretory proteins indicate mast cell activation and testicular fibrosis in the seminal plasma of individuals with cryptorchidism compared to controls. These insights serve as valuable translational references, enriching our comprehension of testicular pathogenesis and informing more precise diagnosis and targeted therapeutic strategies for cryptorchidism.

CD9 and Aryl Hydrocarbon Receptor Are Markers of Human CD19+CD14+ Atypical B Cells and Are Dysregulated in Systemic Lupus Erythematous Disease.

Aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor whose expression regulates immune cell differentiation. Single-cell transcriptomic profiling was used to ascertain the heterogeneity of AHR expression in human B cell subpopulations. We identified a unique population of B cells marked by expression of AHR, CD9, and myeloid genes such as CD14 and CXCL8. Results were confirmed directly in human PBMCs and purified B cells at the protein level. TLR9 signaling induced CD14, CD9, and IL-8 protein expression in CD19+ B cells. CD14-expressing CD9+ B cells also highly expressed AHR and atypical B cell markers such as CD11c and TBET. In patients with active lupus disease, CD14+ and CD9+ B cells are dysregulated, with loss of CD9+ B cells strongly predicting disease severity and demonstrating the relevance of CD9+ B cells in systemic lupus erythematosus and autoimmune disease.

Innate immune responses against mRNA vaccine promote cellular immunity through IFN-β at the injection site

mRNA vaccines against SARS-CoV-2 have revolutionized vaccine development, but their immunological mechanisms are not fully understood. Here, we investigate injection site responses of mRNA vaccines by generating a comprehensive single-cell transcriptome profile upon lipid nanoparticle (LNP) or LNP-mRNA challenge in female BALB/c mice. We show that LNP-induced stromal pro-inflammatory responses and mRNA-elicited type I interferon responses dominate the initial injection site responses. By tracking the fate of delivered mRNA, we discover that injection site fibroblasts are highly enriched with the delivered mRNA and that they express IFN-β specifically in response to the mRNA component, not to the LNP component of mRNA vaccines. Moreover, the mRNA-LNP, but not LNP alone, induces migratory dendritic cells highly expressing IFN-stimulated genes (mDC_ISGs) at the injection site and draining lymph nodes. When co-injected with LNP-subunit vaccine, IFN-β induces mDC_ISGs at the injection site, and importantly, it substantially enhances antigen-specific cellular immune responses. Furthermore, blocking IFN-β signaling at the injection site significantly decreases mRNA vaccine-induced cellular immune responses. Collectively, these data highlight the importance of injection site fibroblasts and IFN-β signaling during early immune responses against the mRNA vaccine and provide detailed information on the initial chain of immune reactions elicited by mRNA vaccine injection.

Single-Cell Time-Resolved Metabolomics and Lipidomics Reveal Apoptotic and Ferroptotic Heterogeneity during Foam Cell Formation.

Macrophage-derived foam cells play a crucial role in plaque formation and rupture during the progression of atherosclerosis. Traditional studies have often overlooked the heterogeneity of foam cells, focusing instead on populations of cells. To address this, we have developed time-resolved, single-cell metabolomics and lipidomics approaches to explore the heterogeneity of macrophages during foam cell formation. Our dynamic metabolomic and lipidomic analyses revealed a dual regulatory axis involving inflammation and ferroptosis. Further, single-cell metabolomics and lipidomics have delineated a continuum of macrophage states, with varied susceptibilities to apoptosis and ferroptosis. Single-cell transcriptomic profiling confirmed these divergent fates, both in established cell lines and in macrophages derived from peripheral blood monocytes. This research has uncovered the complex molecular interactions that dictate these divergent cell fates, providing crucial insights into the pathogenesis of atherosclerosis.

Patient iPSC models reveal glia-intrinsic phenotypes in multiple sclerosis

Multiple sclerosis (MS) is an inflammatory and neurodegenerative disease of the central nervous system (CNS), resulting in neurological disability that worsens over time. While progress has been made in defining the immune system's role in MS pathophysiology, the contribution of intrinsic CNS cell dysfunction remains unclear. Here, we generated a collection of induced pluripotent stem cell (iPSC) lines from people with MS spanning diverse clinical subtypes and differentiated them into glia-enriched cultures. Using single-cell transcriptomic profiling and orthogonal analyses, we observed several distinguishing characteristics of MS cultures pointing to glia-intrinsic disease mechanisms. We found that primary progressive MS-derived cultures contained fewer oligodendrocytes. Moreover, MS-derived oligodendrocyte lineage cells and astrocytes showed increased expression of immune and inflammatory genes, matching those of glia from MS postmortem brains. Thus, iPSC-derived MS models provide a unique platform for dissecting glial contributions to disease phenotypes independent of the peripheral immune system and identify potential glia-specific targets for therapeutic intervention.

Single-Cell Transcriptomic Profiling of Cholangiocyte Organoids Derived from Bile Ducts of Primary Sclerosing Cholangitis Patients

Background and AimsPrimary sclerosing cholangitis (PSC) is a chronic inflammatory liver disorder without effective medical treatment which is characterized by inflammation and fibrotic structures around the bile ducts. Biliary epithelial cells (cholangiocytes) are the target and potential disease drivers in PSC, yet little is known if cholangiocytes from PSC patients differ from non-PSC controls. To characterize cholangiocytes at early rather than end-stage disease, cholangiocyte organoids (COs) were derived from diseased bile ducts of PSC patients and compared to organoids generated from disease controls.MethodsCholangiocytes were obtained during endoscopic retrograde cholangiopancreatography (ERCP) brushing of diseased bile duct areas and expanded as organoids using previously established culture methods. Stable CO lines were analyzed for cell type identity, basic cholangiocyte function, and transcriptomic signature.ResultsWe demonstrate that cholangiocytes, derived from the damaged area within the bile ducts of PSC patients, can be expanded in culture without displaying functional or genetic disease-related features. We further show that COs from patients who later were diagnosed with dysplasia exhibit higher expression of the cancer-associated genes PGC, FXYD2, MIR4435-2HG, and HES1.ConclusionsOur results demonstrate that PSC organoids are largely similar to control organoids after culture and highlight the significance of COs as a tool for regenerative medicine approaches as well as their potential for discovering new potential biomarkers for diagnosing cholangiocarcinoma.

Cross-species imputation and comparison of single-cell transcriptomic profiles.

Cross-species comparison and prediction of gene expression profiles are important to understand regulatory changes during evolution and to transfer knowledge learned from model organisms to humans. Single-cell RNA-seq (scRNA-seq) profiles enable us to capture gene expression profiles with respect to variations among individual cells; however, cross-species comparison of scRNA-seq profiles is challenging because of data sparsity, batch effects, and the lack of one-to-one cell matching across species. Moreover, single-cell profiles are challenging to obtain in certain biological contexts, limiting the scope of hypothesis generation. Here we developed Icebear, a neural network framework that decomposes single-cell measurements into factors representing cell identity, species, and batch factors. Icebear enables accurate prediction of single-cell gene expression profiles across species, thereby providing high-resolution cell type and disease profiles in under-characterized contexts. Icebear also facilitates direct cross-species comparison of single-cell expression profiles for conserved genes that are located on the X chromosome in eutherian mammals but on autosomes in chicken. This comparison, for the first time, revealed evolutionary and diverse adaptations of X-chromosome upregulation in mammals.

Single-cell transcriptome profiles the heterogeneity of tumor cells and microenvironments for different pathological endometrial cancer and identifies specific sensitive drugs

Endometrial cancer (EC) is a highly heterogeneous malignancy characterized by varied pathology and prognoses, and the heterogeneity of its cancer cells and the tumor microenvironment (TME) remains poorly understood. We conducted single-cell RNA sequencing (scRNA-seq) on 18 EC samples, encompassing various pathological types to delineate their specific unique transcriptional landscapes. Cancer cells from diverse pathological sources displayed distinct hallmarks labeled as immune-modulating, proliferation-modulating, and metabolism-modulating cancer cells in uterine clear cell carcinomas (UCCC), well-differentiated endometrioid endometrial carcinomas (EEC-I), and uterine serous carcinomas (USC), respectively. Cancer cells from the UCCC exhibited the greatest heterogeneity. We also identified potential effective drugs and confirmed their effectiveness using patient-derived EC organoids for each pathological group. Regarding the TME, we observed that prognostically favorable CD8+ Tcyto and NK cells were prominent in normal endometrium, whereas CD4+ Treg, CD4+ Tex, and CD8+ Tex cells dominated the tumors. CXCL3+ macrophages associated with M2 signature and angiogenesis were exclusively found in tumors. Prognostically relevant epithelium-specific cancer-associated fibroblasts (eCAFs) and SOD2+ inflammatory CAFs (iCAFs) predominated in EEC-I and UCCC groups, respectively. We also validated the oncogenic effects of SOD2+ iCAFs in vitro. Our comprehensive study has yielded deeper insights into the pathogenesis of EC, potentially facilitating personalized treatments for its varied pathological types.

BCAA-producing Clostridium symbiosum promotes colorectal tumorigenesis through the modulation of host cholesterol metabolism

Identification of potential bacterial players in colorectal tumorigenesis has been a focus of intense research. Herein, we find that Clostridium symbiosum (C.symbiosum) is selectively enriched in tumor tissues of patients with colorectal cancer (CRC) and associated with higher colorectal adenoma recurrence after endoscopic polypectomy. The tumorigenic effect of C.symbiosum is observed in multiple murine models. Single-cell transcriptome profiling along with functional assays demonstrates that C.symbiosum promotes the proliferation of colonic stem cells and enhances cancer stemness. Mechanistically, C.symbiosum intensifies cellular cholesterol synthesis by producing branched-chain amino acids (BCAAs), which sequentially activates Sonic hedgehog signaling. Low dietary BCAA intake or blockade of cholesterol synthesis by statins could partially abrogate the C.symbiosum-induced cell proliferation invivo and invitro. Collectively, we reveal C.symbiosum as a bacterial driver of colorectal tumorigenesis, thus identifying a potential target in CRC prediction, prevention, and treatment.

Atlas of metastatic gastric cancer links ferroptosis to disease progression and immunotherapy response

BACKGROUND & AIMSMetastases from gastric adenocarcinoma (GAC) lead to high morbidity and mortality. Developing innovative and effective therapies requires a comprehensive understanding of the tumor and immune biology of advanced GAC. Yet, collecting matched specimens from advanced, treatment-naïve GAC patients poses a significant challenge, limiting the scope of current research, which has predominantly focused on localized tumors. This gap hinders a deeper insight into the metastatic dynamics of GAC. METHODSWe performed in-depth single-cell transcriptome and immune profiling on 68 paired, treatment-naïve, primary-metastatic tumors to delineate alterations in cancer cells and their tumor microenvironment (TME) during metastatic progression. To validate our observations, we conducted comprehensive functional studies both in vitro and in vivo, employing cell lines, multiple PDX and novel mouse models of GAC. RESULTSLiver and peritoneal metastases exhibited distinct properties in cancer cells and dynamics of TME phenotypes, supporting the notion that cancer cells and their local TMEs co-evolve at metastatic sites. Our study also revealed differential activation of cancer meta-programs across metastases. We observed evasion of cancer cell ferroptosis via GPX4 upregulation during GAC progression. Conditional depletion of Gpx4 or pharmacological inhibition of ferroptosis resistance significantly attenuated tumor growth and metastatic progression. Additionally, ferroptosis-resensitizing treatments augmented the efficacy of CAR T-cell therapy. CONCLUSIONSThis study represents the largest single-cell dataset of metastatic GACs to date. High-resolution mapping of the molecular and cellular dynamics of GAC metastasis has revealed a rationale for targeting ferroptosis defense in combination with CAR T-cell therapy as a novel therapeutic strategy with potential immense clinical implications.

Maternal Innate Immune Reprogramming After Complicated Pregnancy.

Preeclampsia (PE) and fetal growth restriction (FGR) are often associated with maternal inflammation and an increased risk of cardiovascular and metabolic disease in the affected mothers. The mechanism responsible for this increased risk of subsequent disease may involve reprogramming of innate immune cells, characterized by epigenetic modifications. Circulating monocytes from women with PE, FGR, or uncomplicated pregnancies (control) were isolated before labor. Cytokine release from monocytes following exposure to lipopolysaccharide (LPS) and the presence of lysine 4-trimethylated histone 3 (H3K4me3) within TNF promoter sequences were evaluated. Single-cell transcriptomic profiles of circulating monocytes from women with PE or uncomplicated pregnancies were assessed. Monocytes from women with PE or FGR exhibited increased IL-10 secretion and decreased IL-1β and GM-CSF secretion in response to LPS. While TNFα secretion was not significantly different in cultures of control monocytes versus those from complicated pregnancies with or without LPS exposure, monocytes from complicated pregnancies had significantly decreased levels of H3K4me3 associated with TNF promoter sequences. Cluster quantification and pathway analysis of differentially expressed genes revealed an increased proportion of anti-inflammatory myeloid cells and a lower proportion of inflammatory non-classical monocytes among the circulating monocyte population in women with PE. Monocytes from women with PE and FGR exhibit an immune tolerance phenotype before initiation of labor. Further investigation is required to determine whether this tolerogenic phenotype persists after the affected pregnancy and contributes to increased risk of subsequent disease.

Decoding the genetic landscape of juvenile dermatomyositis: insights from phosphorylation-associated single nucleotide polymorphisms.

Genome-wide association studies (GWASs) have identified genetic susceptibility loci associated with juvenile dermatomyositis (JDM). Single nucleotide polymorphisms related to phosphorylation (phosSNPs) are critical nonsynonymous mutations exerting substantial influence on gene expression regulation. The aim of this study was to identify JDM susceptibility genes in the GWAS loci by the use of phosSNPs. We explored quantitative trait loci (QTLs) among the phosSNPs associated with JDM using data from eQTL (bulk tissues and single-cell) and pQTL studies. For gene expression and protein levels significantly influenced by JDM-associated phosSNPs, we assessed their associations with JDM through MR analyses. Additionally, we conducted differential expression gene analyses, incorporating single-cell transcriptomic profiling of 6 JDM cases and 11 juvenile controls (99,396 cells). We identified 31 phosSNPs situated in the 6p21 locus that were associated with JDM. Half of these phosSNPs showed effects on gene expression in various cells and circulating protein levels. In MR analyses, we established associations between the expression levels of pivotal JDM-associated genes, including MICB, C4A, HLA-DRB1, HLA-DRB5, and PSMB9, in skin, muscle, or blood cells and circulating levels of C4A, with JDM. Utilizing single-cell eQTL data, we identified a total of 276 association signals across 14 distinct immune cell types for 28 phosSNPs. Further insights were gained through single-cell differential expression analysis, revealing differential expression of PSMB9, HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQA1, HLA-DQB1, and HLA-DRB1 in immune cells. The present study pinpointed phosSNPs within susceptibility genes for JDM and unraveled the intricate relationships among these SNPs, gene expression levels, and JDM.

Unraveling the Prefrontal Cortex-Basolateral Amygdala Pathway's Role on Schizophrenia's Cognitive Impairments: A Multimodal Study in Patients and Mouse Models.

This study investigated the role of the medial prefrontal cortex (mPFC)-basolateral amygdala (BLA) pathway in schizophrenia (SCZ)-related cognitive impairments using various techniques. This study utilized clinical scales, magnetic resonance imaging, single-cell RNA sequencing, and optogenetics to investigate the mPFC-BLA pathway in SCZ patients. In the mouse model, 6-week-old methylazoxymethanol acetate-induced mice demonstrated significant cognitive deficits, which were addressed through stereotaxic injections of an adeno-associated viral vector to unveil the neural connection between the mPFC and BLA. Significant disparities in brain volume and neural activity, particularly in the dorsolateral prefrontal cortex (DLPFC) and BLA regions, were found between SCZ patients and healthy controls. Additionally, we observed correlations indicating that reduced volumes of the DLPFC and BLA were associated with lower cognitive function scores. Activation of the mPFC-BLA pathway notably improved cognitive performance in the SCZ model mice, with the targeting of excitatory or inhibitory neurons alone failing to replicate this effect. Single-cell transcriptomic profiling revealed gene expression differences in excitatory and inhibitory neurons in the BLA of SCZ model mice. Notably, genes differentially expressed in the BLA of these model mice were also found in the blood exosomes of SCZ patients. Our research provides a comprehensive understanding of the role of the PFC-BLA pathway in SCZ, underscoring its significance in cognitive impairment and offering novel diagnostic and therapeutic avenues. Additionally, our research highlights the potential of blood exosomal mRNAs as noninvasive biomarkers for SCZ diagnosis, underscoring the clinical feasibility and utility of this method.

Single cell transcriptomic profiling identifies tumor-acquired and therapy-resistant cell states in pediatric rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric tumor that resembles undifferentiated muscle cells; yet the extent to which cell state heterogeneity is shared with human development has not been described. Using single-cell/nucleus RNA sequencing from patient tumors, patient-derived xenografts, primary in vitro cultures, and cell lines, we identify four dominant muscle-lineage cell states: progenitor, proliferative, differentiated, and ground cells. We stratify these RMS cells/nuclei along the continuum of human muscle development and show that they share expression patterns with fetal/embryonal myogenic precursors rather than postnatal satellite cells. Fusion-negative RMS (FN-RMS) have a discrete stem cell hierarchy that recapitulates fetal muscle development and contain therapy-resistant FN-RMS progenitors that share transcriptomic similarity with bipotent skeletal mesenchymal cells. Fusion-positive RMS have tumor-acquired cells states, including a neuronal cell state, that are not found in myogenic development. This work identifies previously underappreciated cell state heterogeneity including unique treatment-resistant and tumor-acquired cell states that differ across RMS subtypes.

The 9p21.3 coronary artery disease risk locus drives vascular smooth muscle cells to an osteochondrogenic state.

Genome-wide association studies have identified common genetic variants at ~400 human genomic loci linked to coronary artery disease (CAD) susceptibility. Among these genomic regions, the most impactful is the 9p21.3 CAD risk locus, which spans a 60 kb gene desert and encompasses ~80 SNPs in high linkage disequilibrium. Despite nearly two decades since its discovery, the functional mechanism of this genomic region remains incompletely resolved. To investigate the transcriptional gene programs mediated by 9p21.3 risk locus, we applied a model of induced pluripotent stem cells (iPSCs) from risk and non-risk donors at 9p21.3, as well as isogenic lines with a full haplotype deletion. Upon differentiation to vascular smooth muscle cells (VSMC), single-cell transcriptomic profiling demonstrated iPSC-VSMC phenotypes resembling those from native human coronary arteries, confirming the robustness of this model. Remarkably, our analyses revealed that VSMCs harboring the 9p21.3 risk haplotype preferentially adopt an osteochondrogenic state. Importantly, we identified LIMCH1 and CRABP1 as signature genes critical for defining this transcriptional program. Our study provides new insights into the mechanism at the 9p21.3 risk locus and defines its role in driving a disease-prone transcriptional state in VSMCs.

Redefinition of Synovial Fibroblasts in Rheumatoid Arthritis.

The breakdown of immune tolerance and the rise in autoimmunity contribute to the onset of rheumatoid arthritis (RA), driven by significant changes in immune components. Recent advances in single-cell and spatial transcriptome profiling have revealed shifts in cell distribution and composition, expanding our understanding beyond molecular-level changes in inflammatory cytokines, autoantibodies, and autoantigens in RA. Surprisingly, synovial fibroblasts (SFs) play an active immunopathogenic role rather than remaining passive bystanders in RA, with notable alterations in their subpopulation distribution and composition. This study examines these changes in SF heterogeneity, assesses their impact on RA progression, and elucidates the immune characteristics and functions of SF subsets in the RA autoimmunity, encompassing both intrinsic and adaptive immunity. Additionally, this review discusses therapeutic strategies targeting immune SF subsets, highlighting the potential of future interventions in SF phenotypic reprogramming. Overall, this review redefines the role of SFs in RA and suggests targeting SF phenotypic reprogramming and its upstream molecules as a promising therapeutic approach to restore immune balance and modulate immune tolerance in RA.

TWIST1+FAP+ fibroblasts in the pathogenesis of intestinal fibrosis in Crohn's disease.

Intestinal fibrosis, a severe complication of Crohn's disease (CD), is characterized by excessive extracellular matrix (ECM) deposition and induces intestinal strictures, but there are no effective antifibrosis drugs available for clinical application. We performed single-cell RNA sequencing (scRNA-Seq) of fibrotic and nonfibrotic ileal tissues from patients with CD with intestinal obstruction. Analysis revealed mesenchymal stromal cells (MSCs) as the major producers of ECM and the increased infiltration of its subset FAP+ fibroblasts in fibrotic sites, which was confirmed by immunofluorescence and flow cytometry. Single-cell transcriptomic profiling of chronic dextran sulfate sodium salt murine colitis model revealed that CD81+Pi16- fibroblasts exhibited transcriptomic and functional similarities to human FAP+ fibroblasts. Consistently, FAP+ fibroblasts were identified as the key subtype with the highest level of ECM production in fibrotic intestines. Furthermore, specific knockout or pharmacological inhibition of TWIST1, which was highly expressed by FAP+ fibroblasts, could significantly ameliorate fibrosis in mice. In addition, TWIST1 expression was induced by CXCL9+ macrophages enriched in fibrotic tissues via IL-1β and TGF-β signal. These findings suggest the inhibition of TWIST1 as a promising strategy for CD fibrosis treatment.

PHLDA1-PRDM1 mediates the effect of lentiviral vectors on fate-determination of human retinal progenitor cells

Lentiviral vectors have markedly enhanced gene therapy efficiency in treating congenital diseases, but their long-term safety remains controversial. Most gene therapies for congenital eye diseases need to be carried out at early ages, yet the assessment of related risks to ocular development posed by lentiviral vectors is challenging. Utilizing single-cell transcriptomic profiling on human retinal organoids, this study explored the impact of lentiviral vectors on the retinal development and found that lentiviral vectors can cause retinal precursor cells to shift toward photoreceptor fate through the up-regulation of key fate-determining genes such as PRDM1. Further investigation demonstrated that the intron and intergenic region of PRDM1 was bound by PHLDA1, which was also up-regulated by lentiviral vectors exposure. Importantly, knockdown of PHLDA1 successfully suppressed the lentivirus-induced differentiation bias of photoreceptor cells. The findings also suggest that while lentiviral vectors may disrupt the fate determination of retinal precursor cells, posing risks in early-stage retinal gene therapy, these risks could potentially be reduced by inhibiting the PHLDA1-PRDM1 axis.

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

Aims/hypothesisDisruption 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.MethodsWe 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.ResultsIn 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/interpretationThe 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 availabilityThe 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.Graphical