Single-Cell Analysis of Human Pancreas Reveals Transcriptional Signatures of Aging and Somatic Mutation Patterns

β cell acetate production and release are negligible

BACKGROUND

Studies suggest that short chain fatty acids (SCFAs), which are primarily produced from fermentation of fiber, regulate insulin secretion through free fatty acid receptors 2 and 3 (FFA2 and FFA3). As these are G-protein coupled receptors (GPCRs), they have potential therapeutic value as targets for treating type 2 diabetes (T2D). The exact mechanism by which these receptors regulate insulin secretion and other aspects of pancreatic β cell function is unclear. It has been reported that glucose-dependent release of acetate from pancreatic β cells negatively regulates glucose stimulated insulin secretion. While these data raise the possibility of acetate's potential autocrine action on these receptors, these findings have not been independently confirmed, and multiple concerns exist with this observation, particularly the lack of specificity and precision of the acetate detection methodology used.

METHODS

Using Min6 cells and mouse islets, we assessed acetate and pyruvate production and secretion in response to different glucose concentrations, via liquid chromatography mass spectrometry.

RESULTS

Using Min6 cells and mouse islets, we showed that both intracellular pyruvate and acetate increased with high glucose conditions; however, intracellular acetate level increased only slightly and exclusively in Min6 cells but not in the islets. Further, extracellular acetate levels were not affected by the concentration of glucose in the incubation medium of either Min6 cells or islets.

CONCLUSIONS

Our findings do not substantiate the glucose-dependent release of acetate from pancreatic β cells, and therefore, invalidate the possibility of an autocrine inhibitory effect on glucose stimulated insulin secretion.

Regulation of endocrine cell alternative splicing revealed by single-cell RNA sequencing in type 2 diabetes pathogenesis

The prevalent RNA alternative splicing (AS) contributes to molecular diversity, which has been demonstrated in cellular function regulation and disease pathogenesis. However, the contribution of AS in pancreatic islets during diabetes progression remains unclear. Here, we reanalyze the full-length single-cell RNA sequencing data from the deposited database to investigate AS regulation across human pancreatic endocrine cell types in non-diabetic (ND) and type 2 diabetic (T2D) individuals. Our analysis demonstrates the significant association between transcriptomic AS profiles and cell-type-specificity, which could be applied to distinguish the clustering of major endocrine cell types. Moreover, AS profiles are enabled to clearly define the mature subset of β-cells in healthy controls, which is completely lost in T2D. Further analysis reveals that RNA-binding proteins (RBPs), heterogeneous nuclear ribonucleoproteins (hnRNPs) and FXR1 family proteins are predicted to induce the functional impairment of β-cells through regulating AS profiles. Finally, trajectory analysis of endocrine cells suggests the β-cell identity shift through dedifferentiation and transdifferentiation of β-cells during the progression of T2D. Together, our study provides a mechanism for regulating β-cell functions and suggests the significant contribution of AS program during diabetes pathogenesis.

HumanIslets: An integrated platform for human islet data access and analysis

Comprehensive molecular and cellular phenotyping of human islets can enable deep mechanistic insights for diabetes research. We established the Human Islet Data Analysis and Sharing (HI-DAS) consortium to advance goals in accessibility, usability, and integration of data from human islets isolated from donors with and without diabetes at the Alberta Diabetes Institute (ADI) IsletCore. Here we introduce HumanIslets.com , an open resource for the research community. This platform, which presently includes data on 547 human islet donors, allows users to access linked datasets describing molecular profiles, islet function and donor phenotypes, and to perform various statistical and functional analyses at the donor, islet and single-cell levels. As an example of the analytic capacity of this resource we show a dissociation between cell culture effects on transcript and protein expression, and an approach to correct for exocrine contamination found in hand-picked islets. Finally, we provide an example workflow and visualization that highlights links between type 2 diabetes status, SERCA3b Ca 2+ -ATPase levels at the transcript and protein level, insulin secretion and islet cell phenotypes. HumanIslets.com provides a growing and adaptable set of resources and tools to support the metabolism and diabetes research community.

Gene regulatory networks in disease and ageing

The precise control of gene expression is required for the maintenance of cellular homeostasis and proper cellular function, and the declining control of gene expression with age is considered a major contributor to age-associated changes in cellular physiology and disease. The coordination of gene expression can be represented through models of the molecular interactions that govern gene expression levels, so-called gene regulatory networks. Gene regulatory networks can represent interactions that occur through signal transduction, those that involve regulatory transcription factors, or statistical models of gene-gene relationships based on the premise that certain sets of genes tend to be coexpressed across a range of conditions and cell types. Advances in experimental and computational technologies have enabled the inference of these networks on an unprecedented scale and at unprecedented precision. Here, we delineate different types of gene regulatory networks and their cell-biological interpretation. We describe methods for inferring such networks from large-scale, multi-omics datasets and present applications that have aided our understanding of cellular ageing and disease mechanisms.

Remodeling ceramide homeostasis promotes functional maturation of human pluripotent stem cell-derived β cells

Human pluripotent stem cell-derived β cells (hPSC-β cells) show the potential to restore euglycemia. However, the immature functionality of hPSC-β cells has limited their efficacy in application. Here, by deciphering the continuous maturation process of hPSC-β cells post transplantation via single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), we show that functional maturation of hPSC-β cells is an orderly multistep process during which cells sequentially undergo metabolic adaption, removal of negative regulators of cell function, and establishment of a more specialized transcriptome and epigenome. Importantly, remodeling lipid metabolism, especially downregulating the metabolic activity of ceramides, the central hub of sphingolipid metabolism, is critical for β cell maturation. Limiting intracellular accumulation of ceramides in hPSC-β cells remarkably enhanced their function, as indicated by improvements in insulin processing and glucose-stimulated insulin secretion. In summary, our findings provide insights into the maturation of human pancreatic β cells and highlight the importance of ceramide homeostasis in function acquisition.

Application potential of senolytics in clinical treatment

Of the factors studied in individual ageing, the accumulation of senescent cells has been considered as an essential cause of organ degeneration to eventually initiate age-related diseases. Cellular senescence is attributed to the accumulation of damage for an inducement in the activation of cell cycle inhibitory pathways, resulting the cell permanently withdraw from the cell proliferation cycle. Further, senescent cells will activate the inflammatory factor secretion pathway to promote the development of various age-related diseases. Senolytics, a small molecule compound, can delay disease development and extend mammalian lifespan. The evidence from multiple trials shows that the targeted killing of senescent cells has a significant clinical application for the treatment of age-related diseases. In addition, senolytics are also significant for the development of ageing research in solid organ transplantation, which can fully develop the potential of elderly organs and reduce the age gap between demand and supply. We conclude that the main characteristics of cellular senescence, the anti-ageing drug senolytics in the treatment of chronic diseases and organ transplantation, and the latest clinical progress of related researches in order to provide a theoretical basis for the prevention and treatment of ageing and related diseases.

The age-dependent regulation of pancreatic islet landscape is fueled by a HNF1a-immune signaling loop

Animal longevity is a function of global vital organ functionality and, consequently, a complex polygenic trait. Yet, monogenic regulators controlling overall or organ-specific ageing exist, owing their conservation to their function in growth and development. Here, by using pathway analysis combined with wet-biology methods on several dynamic timelines, we identified Hnf1a as a novel master regulator of the maturation and ageing in the adult pancreatic islet during the first year of life. Conditional transgenic mice bearing suboptimal levels of this transcription factor in the pancreatic islets displayed age-dependent changes, with a profile echoing precocious maturation. Additionally, the comparative pathway analysis revealed a link between Hnf1a age-dependent regulation and immune signaling, which was confirmed in the ageing timeline of an overly immunodeficient mouse model. Last, the global proteome analysis of human islets spanning three decades of life largely backed the age-specific regulation observed in mice. Collectively, our results suggest a novel role of Hnf1a as a monogenic regulator of the maturation and ageing process in the pancreatic islet via a direct or indirect regulatory loop with immune signaling.

Spatiotemporal multi-omics: exploring molecular landscapes in aging and regenerative medicine

Aging and regeneration represent complex biological phenomena that have long captivated the scientific community. To fully comprehend these processes, it is essential to investigate molecular dynamics through a lens that encompasses both spatial and temporal dimensions. Conventional omics methodologies, such as genomics and transcriptomics, have been instrumental in identifying critical molecular facets of aging and regeneration. However, these methods are somewhat limited, constrained by their spatial resolution and their lack of capacity to dynamically represent tissue alterations. The advent of emerging spatiotemporal multi-omics approaches, encompassing transcriptomics, proteomics, metabolomics, and epigenomics, furnishes comprehensive insights into these intricate molecular dynamics. These sophisticated techniques facilitate accurate delineation of molecular patterns across an array of cells, tissues, and organs, thereby offering an in-depth understanding of the fundamental mechanisms at play. This review meticulously examines the significance of spatiotemporal multi-omics in the realms of aging and regeneration research. It underscores how these methodologies augment our comprehension of molecular dynamics, cellular interactions, and signaling pathways. Initially, the review delineates the foundational principles underpinning these methods, followed by an evaluation of their recent applications within the field. The review ultimately concludes by addressing the prevailing challenges and projecting future advancements in the field. Indubitably, spatiotemporal multi-omics are instrumental in deciphering the complexities inherent in aging and regeneration, thus charting a course toward potential therapeutic innovations.

scPLAN: a hierarchical computational framework for single transcriptomics data annotation, integration and cell-type label refinement

MOTIVATION

In the past decade, single-cell RNA sequencing (scRNA-seq) has emerged as a pivotal method for transcriptomic profiling in biomedical research. Precise cell-type identification is crucial for subsequent analysis of single-cell data. And the integration and refinement of annotated data are essential for building comprehensive databases. However, prevailing annotation techniques often overlook the hierarchical organization of cell types, resulting in inconsistent annotations. Meanwhile, most existing integration approaches fail to integrate datasets with different annotation depths and none of them can enhance the labels of outdated data with lower annotation resolutions using more intricately annotated datasets or novel biological findings.

RESULTS

Here, we introduce scPLAN, a hierarchical computational framework designed for scRNA-seq data analysis. scPLAN excels in annotating unlabeled scRNA-seq data using a reference dataset structured along a hierarchical cell-type tree. It identifies potential novel cell types in a systematic, layer-by-layer manner. Additionally, scPLAN effectively integrates annotated scRNA-seq datasets with varying levels of annotation depth, ensuring consistent refinement of cell-type labels across datasets with lower resolutions. Through extensive annotation and novel cell detection experiments, scPLAN has demonstrated its efficacy. Two case studies have been conducted to showcase how scPLAN integrates datasets with diverse cell-type label resolutions and refine their cell-type labels.

AVAILABILITY

https://github.com/michaelGuo1204/scPLAN.

Netrin G1 Ligand is a new stromal immunomodulator that promotes pancreatic cancer

Understanding pancreatic cancer biology is fundamental for identifying new targets and for developing more effective therapies. In particular, the contribution of the stromal microenvironment to pancreatic cancer tumorigenesis requires further exploration. Here, we report the stromal roles of the synaptic protein Netrin G1 Ligand (NGL-1) in pancreatic cancer, uncovering its pro-tumor functions in cancer-associated fibroblasts and in immune cells. We observed that the stromal expression of NGL-1 inversely correlated with patients' overall survival. Moreover, germline knockout (KO) mice for NGL-1 presented decreased tumor burden, with a microenvironment that is less supportive of tumor growth. Of note, tumors from NGL-1 KO mice produced less immunosuppressive cytokines and displayed an increased percentage of CD8 + T cells than those from control mice, while preserving the physical structure of the tumor microenvironment. These effects were shown to be mediated by NGL-1 in both immune cells and in the local stroma, in a TGF-β-dependent manner. While myeloid cells lacking NGL-1 decreased the production of immunosuppressive cytokines, NGL-1 KO T cells showed increased proliferation rates and overall polyfunctionality compared to control T cells. CAFs lacking NGL-1 were less immunosuppressive than controls, with overall decreased production of pro-tumor cytokines and compromised ability to inhibit CD8 + T cells activation. Mechanistically, these CAFs downregulated components of the TGF-β pathway, AP-1 and NFAT transcription factor families, resulting in a less tumor-supportive phenotype. Finally, targeting NGL-1 genetically or using a functionally antagonistic small peptide phenocopied the effects of chemotherapy, while modulating the immunosuppressive tumor microenvironment (TME), rather than eliminating it. We propose NGL-1 as a new local stroma and immunomodulatory molecule, with pro-tumor roles in pancreatic cancer.

Statement of Significance

Here we uncovered the pro-tumor roles of the synaptic protein NGL-1 in the tumor microenvironment of pancreatic cancer, defining a new target that simultaneously modulates tumor cell, fibroblast, and immune cell functions. This study reports a new pathway where NGL-1 controls TGF-β, AP-1 transcription factor members and NFAT1, modulating the immunosuppressive microenvironment in pancreatic cancer. Our findings highlight NGL-1 as a new stromal immunomodulator in pancreatic cancer.

Human skeletal muscle aging atlas

Skeletal muscle aging is a key contributor to age-related frailty and sarcopenia with substantial implications for global health. Here we profiled 90,902 single cells and 92,259 single nuclei from 17 donors to map the aging process in the adult human intercostal muscle, identifying cellular changes in each muscle compartment. We found that distinct subsets of muscle stem cells exhibit decreased ribosome biogenesis genes and increased CCL2 expression, causing different aging phenotypes. Our atlas also highlights an expansion of nuclei associated with the neuromuscular junction, which may reflect re-innervation, and outlines how the loss of fast-twitch myofibers is mitigated through regeneration and upregulation of fast-type markers in slow-twitch myofibers with age. Furthermore, we document the function of aging muscle microenvironment in immune cell attraction. Overall, we present a comprehensive human skeletal muscle aging resource ( https://www.muscleageingcellatlas.org/ ) together with an in-house mouse muscle atlas to study common features of muscle aging across species.

NKX2-2 based nuclei sorting on frozen human archival pancreas enables the enrichment of islet endocrine populations for single-nucleus RNA sequencing

BACKGROUND

Current approaches to profile the single-cell transcriptomics of human pancreatic endocrine cells almost exclusively rely on freshly isolated islets. However, human islets are limited in availability. Furthermore, the extensive processing steps during islet isolation and subsequent single cell dissolution might alter gene expressions. In this work, we report the development of a single-nucleus RNA sequencing (snRNA-seq) approach with targeted islet cell enrichment for endocrine-population focused transcriptomic profiling using frozen archival pancreatic tissues without islet isolation.

RESULTS

We cross-compared five nuclei isolation protocols and selected the citric acid method as the best strategy to isolate nuclei with high RNA integrity and low cytoplasmic contamination from frozen archival human pancreata. We innovated fluorescence-activated nuclei sorting based on the positive signal of NKX2-2 antibody to enrich nuclei of the endocrine population from the entire nuclei pool of the pancreas. Our sample preparation procedure generated high-quality single-nucleus gene-expression libraries while preserving the endocrine population diversity. In comparison with single-cell RNA sequencing (scRNA-seq) library generated with live cells from freshly isolated human islets, the snRNA-seq library displayed comparable endocrine cellular composition and cell type signature gene expression. However, between these two types of libraries, differential enrichments of transcripts belonging to different functional classes could be observed.

CONCLUSIONS

Our work fills a technological gap and helps to unleash frozen archival pancreatic tissues for molecular profiling targeting the endocrine population. This study opens doors to retrospective mappings of endocrine cell dynamics in pancreatic tissues of complex histopathology. We expect that our protocol is applicable to enrich nuclei for transcriptomics studies from various populations in different types of frozen archival tissues.

LDHB contributes to the regulation of lactate levels and basal insulin secretion in human pancreatic β cells

Using 13C6 glucose labeling coupled to gas chromatography-mass spectrometry and 2D 1H-13C heteronuclear single quantum coherence NMR spectroscopy, we have obtained a comparative high-resolution map of glucose fate underpinning β cell function. In both mouse and human islets, the contribution of glucose to the tricarboxylic acid (TCA) cycle is similar. Pyruvate fueling of the TCA cycle is primarily mediated by the activity of pyruvate dehydrogenase, with lower flux through pyruvate carboxylase. While the conversion of pyruvate to lactate by lactate dehydrogenase (LDH) can be detected in islets of both species, lactate accumulation is 6-fold higher in human islets. Human islets express LDH, with low-moderate LDHA expression and β cell-specific LDHB expression. LDHB inhibition amplifies LDHA-dependent lactate generation in mouse and human β cells and increases basal insulin release. Lastly, cis-instrument Mendelian randomization shows that low LDHB expression levels correlate with elevated fasting insulin in humans. Thus, LDHB limits lactate generation in β cells to maintain appropriate insulin release.

Contributions of transcriptional noise to leukaemia evolution: KAT2A as a case-study

Transcriptional noise is proposed to participate in cell fate changes, but contributions to mammalian cell differentiation systems, including cancer, remain associative. Cancer evolution is driven by genetic variability, with modulatory or contributory participation of epigenetic variants. Accumulation of epigenetic variants enhances transcriptional noise, which can facilitate cancer cell fate transitions. Acute myeloid leukaemia (AML) is an aggressive cancer with strong epigenetic dependencies, characterized by blocked differentiation. It constitutes an attractive model to probe links between transcriptional noise and malignant cell fate regulation. Gcn5/KAT2A is a classical epigenetic transcriptional noise regulator. Its loss increases transcriptional noise and modifies cell fates in stem and AML cells. By reviewing the analysis of KAT2A-depleted pre-leukaemia and leukaemia models, I discuss that the net result of transcriptional noise is diversification of cell fates secondary to alternative transcriptional programmes. Cellular diversification can enable or hinder AML progression, respectively, by differentiation of cell types responsive to mutations, or by maladaptation of leukaemia stem cells. KAT2A-dependent noise-responsive genes participate in ribosome biogenesis and KAT2A loss destabilizes translational activity. I discuss putative contributions of perturbed translation to AML biology, and propose KAT2A loss as a model for mechanistic integration of transcriptional and translational control of noise and fate decisions. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

A single-cell atlas of lung homeostasis reveals dynamic changes during development and aging

Aging is a global challenge, marked in the lungs by function decline and structural disorders, which affects the health of the elderly population. To explore anti-aging strategies, we develop a dynamic atlas covering 45 cell types in human lungs, spanning from embryonic development to aging. We aim to apply the discoveries of lung's development to address aging-related issues. We observe that both epithelial and immune cells undergo a process of acquisition and loss of essential function as they transition from development to aging. During aging, we identify cellular phenotypic alternations that result in reduced pulmonary compliance and compromised immune homeostasis. Furthermore, we find a distinctive expression pattern of the ferritin light chain (FTL) gene, which increases during development but decreases in various types of lung cells during the aging process.

Polyploidisation pleiotropically buffers ageing in hepatocytes

BACKGROUND & AIMS

Polyploidy in hepatocytes has been proposed as a genetic mechanism to buffer against transcriptional dysregulation. Here, we aim to demonstrate the role of polyploidy in modulating gene regulatory networks in hepatocytes during ageing.

METHODS

We performed single-nucleus RNA sequencing in hepatocyte nuclei of different ploidy levels isolated from young and old wild-type mice. Changes in the gene expression and regulatory network were compared to three independent strains that were haploinsufficient for HNF4A, CEBPA or CTCF, representing non-deleterious perturbations. Phenotypic characteristics of the liver section were additionally evaluated histologically, whereas the genomic allele composition of hepatocytes was analysed by BaseScope.

RESULTS

We observed that ageing in wild-type mice results in nuclei polyploidy and a marked increase in steatosis. Haploinsufficiency of liver-specific master regulators (HFN4A or CEBPA) results in the enrichment of hepatocytes with tetraploid nuclei at a young age, affecting the genomic regulatory network, and dramatically suppressing ageing-related steatosis tissue wide. Notably, these phenotypes are not the result of subtle disruption to liver-specific transcriptional networks, since haploinsufficiency in the CTCF insulator protein resulted in the same phenotype. Further quantification of genotypes of tetraploid hepatocytes in young and old HFN4A-haploinsufficient mice revealed that during ageing, tetraploid hepatocytes lead to the selection of wild-type alleles, restoring non-deleterious genetic perturbations.

CONCLUSIONS

Our results suggest a model whereby polyploidisation leads to fundamentally different cell states. Polyploid conversion enables pleiotropic buffering against age-related decline via non-random allelic segregation to restore a wild-type genome.

IMPACT AND IMPLICATIONS

The functional role of hepatocyte polyploidisation during ageing is poorly understood. Using single-nucleus RNA sequencing and BaseScope approaches, we have studied ploidy dynamics during ageing in murine livers with non-deleterious genetic perturbations. We have identified that hepatocytes present different cellular states and the ability to buffer ageing-associated dysfunctions. Tetraploid nuclei exhibit robust transcriptional networks and are better adapted to genomically overcome perturbations. Novel therapeutic interventions aimed at attenuating age-related changes in tissue function could be exploited by manipulation of ploidy dynamics during chronic liver conditions.

Methylation entropy landscape of Chinese long-lived individuals reveals lower epigenetic noise related to human healthy aging

The transition from ordered to noisy is a significant epigenetic signature of aging and age-related disease. As a paradigm of healthy human aging and longevity, long-lived individuals (LLI, >90 years old) may possess characteristic strategies in coping with the disordered epigenetic regulation. In this study, we constructed high-resolution blood epigenetic noise landscapes for this cohort by a methylation entropy (ME) method using whole genome bisulfite sequencing (WGBS). Although a universal increase in global ME occurred with chronological age in general control samples, this trend was suppressed in LLIs. Importantly, we identified 38,923 genomic regions with LLI-specific lower ME (LLI-specific lower entropy regions, for short, LLI-specific LERs). These regions were overrepresented in promoters, which likely function in transcriptional noise suppression. Genes associated with LLI-specific LERs have a considerable impact on SNP-based heritability of some aging-related disorders (e.g., asthma and stroke). Furthermore, neutrophil was identified as the primary cell type sustaining LLI-specific LERs. Our results highlight the stability of epigenetic order in promoters of genes involved with aging and age-related disorders within LLI epigenomes. This unique epigenetic feature reveals a previously unknown role of epigenetic order maintenance in specific genomic regions of LLIs, which helps open a new avenue on the epigenetic regulation mechanism in human healthy aging and longevity.

Somatic mutations in human ageing: New insights from DNA sequencing and inherited mutations

The accumulation of somatic mutations is a driver of cancer and has long been associated with ageing. Due to limitations in quantifying mutation burden with age in non-cancerous tissues, the impact of somatic mutations in other ageing phenotypes is unclear. Recent advances in DNA sequencing technologies have allowed the large-scale quantification of somatic mutations in ageing tissues. These studies have revealed a gradual accumulation of mutations in normal tissues with age as well as a substantial clonal expansion driven mostly by cancer-related mutations. Nevertheless, it is difficult to envision how the burden and stochastic nature of age-related somatic mutations identified so far can explain most ageing phenotypes that develop gradually. Studies across species have also found that longer-lived species have lower somatic mutation rates, though these could be due to selective pressures acting on other phenotypes such as perhaps cancer. Recent studies in patients with higher somatic mutation burden and no signs of accelerated ageing further question the role of somatic mutations in ageing. Overall, with a few exceptions like cancer, recent DNA sequencing studies and inherited mutations do not support the idea that somatic mutations accumulating with age drive ageing phenotypes, and the phenotypic role, if any, of somatic mutations in ageing remains unclear.

Effects of a bitter substance, denatonium benzoate, on pancreatic hormone secretion

There is increasing evidence linking bitter taste receptor (BTR) signaling to gut hormone secretion and glucose homeostasis. However, its effect on islet hormone secretion has been poorly characterized. This study investigated the effect of the bitter substance, denatonium benzoate (DB), on hormone secretion from mouse pancreatic islets and INS-1 832/13 cells. DB (0.5-1 mM) augmented insulin secretion at both 2.8 mM and 16.7 mM glucose. This effect was no longer present at 5 mM DB likely due to the greater levels of cellular apoptosis. DB-stimulated insulin secretion involved closure of the KATP channel, activation of T2R signaling in beta-cells, and intraislet glucagon-like peptide-1 (GLP-1) release. DB also enhanced glucagon and somatostatin secretion, but the underlying mechanism was less clear. Together, this study demonstrates that the bitter substance, DB, is a strong potentiator of islet hormone secretion independent of glucose. This observation highlights the potential for widespread off-target effects associated with the clinical use of bitter-tasting substances.NEW & NOTEWORTHY We show that the bitter substance, denatonium benzoate (DB), stimulates insulin, glucagon, somatostatin, and GLP-1 secretion from pancreatic islets, independent of glucose, and that DB augments insulin release via the KATP channel, bitter taste receptor signaling, and intraislet GLP-1 secretion. Exposure to a high dose of DB (5 mM) induces cellular apoptosis in pancreatic islets. Therefore, clinical use of bitter substances to improve glucose homeostasis may have unintended negative impacts beyond the gut.

Somatic mutations in aging and disease

Time always leaves its mark, and our genome is no exception. Mutations in the genome of somatic cells were first hypothesized to be the cause of aging in the 1950s, shortly after the molecular structure of DNA had been described. Somatic mutation theories of aging are based on the fact that mutations in DNA as the ultimate template for all cellular functions are irreversible. However, it took until the 1990s to develop the methods to test if DNA mutations accumulate with age in different organs and tissues and estimate the severity of the problem. By now, numerous studies have documented the accumulation of somatic mutations with age in normal cells and tissues of mice, humans, and other animals, showing clock-like mutational signatures that provide information on the underlying causes of the mutations. In this review, we will first briefly discuss the recent advances in next-generation sequencing that now allow quantitative analysis of somatic mutations. Second, we will provide evidence that the mutation rate differs between cell types, with a focus on differences between germline and somatic mutation rate. Third, we will discuss somatic mutational signatures as measures of aging, environmental exposure, and activities of DNA repair processes. Fourth, we will explain the concept of clonally amplified somatic mutations, with a focus on clonal hematopoiesis. Fifth, we will briefly discuss somatic mutations in the transcriptome and in our other genome, i.e., the genome of mitochondria. We will end with a brief discussion of a possible causal contribution of somatic mutations to the aging process.

scFSNN: a feature selection method based on neural network for single-cell RNA-seq data

While single-cell RNA sequencing (scRNA-seq) allows researchers to analyze gene expression in individual cells, its unique characteristics like over-dispersion, zero-inflation, high gene-gene correlation, and large data volume with many features pose challenges for most existing feature selection methods. In this paper, we present a feature selection method based on neural network (scFSNN) to solve classification problem for the scRNA-seq data. scFSNN is an embedded method that can automatically select features (genes) during model training, control the false discovery rate of selected features and adaptively determine the number of features to be eliminated. Extensive simulation and real data studies demonstrate its excellent feature selection ability and predictive performance.

Glucagon-producing α-cell transcriptional identity and reprogramming towards insulin production

β-Cell replacement by in situ reprogramming of non-β-cells is a promising diabetes therapy. Following the observation that near-total β-cell ablation in adult mice triggers the reprogramming of pancreatic α-, δ-, and γ-cells into insulin (INS)-producing cells, recent studies are delving deep into the mechanisms controlling adult α-cell identity. Systematic analyses of the α-cell transcriptome and epigenome have started to pinpoint features that could be crucial for maintaining α-cell identity. Using different transgenic and chemical approaches, significant advances have been made in reprogramming α-cells in vivo into INS-secreting cells in mice. The recent reprogramming of human α-cells in vitro is an important step forward that must now be complemented with a comprehensive molecular dissection of the mechanisms controlling α-cell identity.

Aging induces region-specific dysregulation of hormone synthesis in the primate adrenal gland

Adrenal glands, vital for steroid secretion and the regulation of metabolism, stress responses and immune activation, experience age-related decline, impacting systemic health. However, the regulatory mechanisms underlying adrenal aging remain largely uninvestigated. Here we established a single-nucleus transcriptomic atlas of both young and aged primate suprarenal glands, identifying lipid metabolism and steroidogenic pathways as core processes impacted by aging. We found dysregulation in centripetal adrenocortical differentiation in aged adrenal tissues and cells in the zona reticularis region, responsible for producing dehydroepiandrosterone sulfate (DHEA-S), were highly susceptible to aging, reflected by senescence, exhaustion and disturbed hormone production. Remarkably, LDLR was downregulated in all cell types of the outer cortex, and its targeted inactivation in human adrenal cells compromised cholesterol uptake and secretion of dehydroepiandrosterone sulfate, as observed in aged primate adrenal glands. Our study provides crucial insights into endocrine physiology, holding therapeutic promise for addressing aging-related adrenal insufficiency and delaying systemic aging.

Combining Global-Constrained Concept Factorization and a Regularized Gaussian Graphical Model for Clustering Single-Cell RNA-seq Data

Single-cell RNA sequencing technology is one of the most cost-effective ways to uncover transcriptomic heterogeneity. With the rapid rise of this technology, enormous amounts of scRNA-seq data have been produced. Due to the high dimensionality, noise, sparsity and missing features of the available scRNA-seq data, accurately clustering the scRNA-seq data for downstream analysis is a significant challenge. Many computational methods have been designed to address this issue; nevertheless, the efficacy of the available methods is still inadequate. In addition, most similarity-based methods require a number of clusters as input, which is difficult to achieve in real applications. In this study, we developed a novel computational method for clustering scRNA-seq data by considering both global and local information, named GCFG. This method characterizes the global properties of data by applying concept factorization, and the regularized Gaussian graphical model is utilized to evaluate the local embedding relationship of data. To learn the cell-cell similarity matrix, we integrated the two components, and an iterative optimization algorithm was developed. The categorization of single cells is obtained by applying Louvain, a modularity-based community discovery algorithm, to the similarity matrix. The behavior of the GCFG approach is assessed on 14 real scRNA-seq datasets in terms of ACC and ARI, and comparison results with 17 other competitive methods suggest that GCFG is effective and robust.

New generative methods for single-cell transcriptome data in bulk RNA sequence deconvolution

Numerous methods for bulk RNA sequence deconvolution have been developed to identify cellular targets of diseases by understanding the composition of cell types in disease-related tissues. However, issues of heterogeneity in gene expression between subjects and the shortage of reference single-cell RNA sequence data remain to achieve accurate bulk deconvolution. In our study, we investigated whether a new data generative method named sc-CMGAN and benchmarking generative methods (Copula, CTGAN and TVAE) could solve these issues and improve the bulk deconvolutions. We also evaluated the robustness of sc-CMGAN using three deconvolution methods and four public datasets. In almost all conditions, the generative methods contributed to improved deconvolution. Notably, sc-CMGAN outperformed the benchmarking methods and demonstrated higher robustness. This study is the first to examine the impact of data augmentation on bulk deconvolution. The new generative method, sc-CMGAN, is expected to become one of the powerful tools for the preprocessing of bulk deconvolution.

scEM: A New Ensemble Framework for Predicting Cell Type Composition Based on scRNA-Seq Data

With the advent of single-cell RNA sequencing (scRNA-seq) technology, many scRNA-seq data have become available, providing an unprecedented opportunity to explore cellular composition and heterogeneity. Recently, many computational algorithms for predicting cell type composition have been developed, and these methods are typically evaluated on different datasets and performance metrics using diverse techniques. Consequently, the lack of comprehensive and standardized comparative analysis makes it difficult to gain a clear understanding of the strengths and weaknesses of these methods. To address this gap, we reviewed 20 cutting-edge unsupervised cell type identification methods and evaluated these methods comprehensively using 24 real scRNA-seq datasets of varying scales. In addition, we proposed a new ensemble cell-type identification method, named scEM, which learns the consensus similarity matrix by applying the entropy weight method to the four representative methods are selected. The Louvain algorithm is adopted to obtain the final classification of individual cells based on the consensus matrix. Extensive evaluation and comparison with 11 other similarity-based methods under real scRNA-seq datasets demonstrate that the newly developed ensemble algorithm scEM is effective in predicting cellular type composition.

Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

The cycling and aging mouse female reproductive tract at single-cell resolution

The female reproductive tract (FRT) undergoes extensive remodeling during reproductive cycling. This recurrent remodeling and how it shapes organ-specific aging remains poorly explored. Using single-cell and spatial transcriptomics, we systematically characterized morphological and gene expression changes occurring in ovary, oviduct, uterus, cervix, and vagina at each phase of the mouse estrous cycle, during decidualization, and into aging. These analyses reveal that fibroblasts play central-and highly organ-specific-roles in FRT remodeling by orchestrating extracellular matrix (ECM) reorganization and inflammation. Our results suggest a model wherein recurrent FRT remodeling over reproductive lifespan drives the gradual, age-related development of fibrosis and chronic inflammation. This hypothesis was directly tested using chemical ablation of cycling, which reduced fibrotic accumulation during aging. Our atlas provides extensive detail into how estrus, pregnancy, and aging shape the organs of the female reproductive tract and reveals the unexpected cost of the recurrent remodeling required for reproduction.

Establishing evidence for immune surveillance of β-cell senescence

Cellular senescence is a programmed state of cell cycle arrest that involves a complex immunogenic secretome, eliciting immune surveillance and senescent cell clearance. Recent work has shown that a subpopulation of pancreatic β-cells becomes senescent in the context of diabetes; however, it is not known whether these cells are normally subject to immune surveillance. In this opinion article, we advance the hypothesis that immune surveillance of β-cells undergoing a senescence stress response normally limits their accumulation during aging and that the breakdown of these mechanisms is a driver of senescent β-cell accumulation in diabetes. Elucidation and therapeutic activation of immune surveillance mechanisms in the pancreas holds promise for the improvement of approaches to target stressed senescent β-cells in the treatment of diabetes.

Deciphering Immune Landscape Remodeling Unravels the Underlying Mechanism for Synchronized Muscle and Bone Aging

Evidence from numerous studies has revealed the synchronous progression of aging in bone and muscle; however, little is known about the underlying mechanisms. To this end, human muscles and bones are harvested and the aging-associated transcriptional dynamics of two tissues in parallel using single-cell RNA sequencing are surveyed. A subset of lipid-associated macrophages (triggering receptor expressed on myeloid cells 2, TREM2+ Macs) is identified in both aged muscle and bone. Genes responsible for muscle dystrophy and bone loss, such as secreted phosphoprotein 1 (SPP1), are also highly expressed in TREM2+ Macs, suggesting its conserved role in aging-related features. A common transition toward pro-inflammatory phenotypes in aged CD4+ T cells across tissues is also observed, activated by the nuclear factor kappa B subunit 1 (NFKB1). CD4+ T cells in aged muscle experience Th1-like differentiation, whereas, in bone, a skewing toward Th17 cells is observed. Furthermore, these results highlight that degenerated myocytes produce BAG6-containing exosomes that can communicate with Th17 cells in the bone through its receptor natural cytotoxicity triggering receptor 3 (NCR3). This communication upregulates CD6 expression in Th17 cells, which then interact with TREM2+ Macs through CD6-ALCAM signaling, ultimately stimulating the transcription of SPP1 in TREM2+ Macs. The negative correlation between serum exosomal BCL2-associated athanogene 6 (BAG6) levels and bone mineral density further supports its role in mediating muscle and bone synchronization with aging.

Deciphering Drug Targets and Actions with Single-Cell and Spatial Resolution

Recent advances in chemical, molecular, and genetic approaches have provided us with an unprecedented capacity to identify drug-target interactions across the whole proteome and genome. Meanwhile, rapid developments of single-cell and spatial omics technologies are revolutionizing our understanding of the molecular architecture of biological systems. However, a significant gap remains in how we align our understanding of drug actions, traditionally based on molecular affinities, with the in vivo cellular and spatial tissue heterogeneity revealed by these newer techniques. Here, we review state-of-the-art methods for profiling drug-target interactions and emerging multiomics tools to delineate the tissue heterogeneity at single-cell resolution. Highlighting the recent technical advances enabling high-resolution, multiplexable in situ small-molecule drug imaging (clearing-assisted tissue click chemistry, or CATCH), we foresee the integration of single-cell and spatial omics platforms, data, and concepts into the future framework of defining and understanding in vivo drug-target interactions and mechanisms of actions.

Continually adapting pre-trained language model to universal annotation of single-cell RNA-seq data

MOTIVATION

Cell-type annotation of single-cell RNA-sequencing (scRNA-seq) data is a hallmark of biomedical research and clinical application. Current annotation tools usually assume the simultaneous acquisition of well-annotated data, but without the ability to expand knowledge from new data. Yet, such tools are inconsistent with the continuous emergence of scRNA-seq data, calling for a continuous cell-type annotation model. In addition, by their powerful ability of information integration and model interpretability, transformer-based pre-trained language models have led to breakthroughs in single-cell biology research. Therefore, the systematic combining of continual learning and pre-trained language models for cell-type annotation tasks is inevitable.

RESULTS

We herein propose a universal cell-type annotation tool, called CANAL, that continuously fine-tunes a pre-trained language model trained on a large amount of unlabeled scRNA-seq data, as new well-labeled data emerges. CANAL essentially alleviates the dilemma of catastrophic forgetting, both in terms of model inputs and outputs. For model inputs, we introduce an experience replay schema that repeatedly reviews previous vital examples in current training stages. This is achieved through a dynamic example bank with a fixed buffer size. The example bank is class-balanced and proficient in retaining cell-type-specific information, particularly facilitating the consolidation of patterns associated with rare cell types. For model outputs, we utilize representation knowledge distillation to regularize the divergence between previous and current models, resulting in the preservation of knowledge learned from past training stages. Moreover, our universal annotation framework considers the inclusion of new cell types throughout the fine-tuning and testing stages. We can continuously expand the cell-type annotation library by absorbing new cell types from newly arrived, well-annotated training datasets, as well as automatically identify novel cells in unlabeled datasets. Comprehensive experiments with data streams under various biological scenarios demonstrate the versatility and high model interpretability of CANAL.

AVAILABILITY

An implementation of CANAL is available from https://github.com/aster-ww/CANAL-torch.

CONTACT

dengmh@pku.edu.cn.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Journal Name online.

PancrESS - a meta-analysis resource for understanding cell-type specific expression in the human pancreas

BACKGROUND

The human pancreas is composed of specialized cell types producing hormones and enzymes critical to human health. These specialized functions are the result of cell type-specific transcriptional programs which manifest in cell-specific gene expression. Understanding these programs is essential to developing therapies for pancreatic disorders. Transcription in the human pancreas has been widely studied by single-cell RNA technologies, however the diversity of protocols and analysis methods hinders their interpretability in the aggregate.

RESULTS

In this work, we perform a meta-analysis of pancreatic single-cell RNA sequencing data. We present a database for reference transcriptome abundances and cell-type specificity metrics. This database facilitates the identification and definition of marker genes within the pancreas. Additionally, we introduce a versatile tool which is freely available as an R package, and should permit integration into existing workflows. Our tool accepts count data files generated by widely-used single-cell gene expression platforms in their original format, eliminating an additional pre-formatting step. Although we designed it to calculate expression specificity of pancreas cell types, our tool is agnostic to the biological source of count data, extending its applicability to other biological systems.

CONCLUSIONS

Our findings enhance the current understanding of expression specificity within the pancreas, surpassing previous work in terms of scope and detail. Furthermore, our database and tool enable researchers to perform similar calculations in diverse biological systems, expanding the applicability of marker gene identification and facilitating comparative analyses.

Gut epithelial Interleukin-17 receptor A signaling can modulate distant tumors growth through microbial regulation

Microbes influence cancer initiation, progression and therapy responsiveness. IL-17 signaling contributes to gut barrier immunity by regulating microbes but also drives tumor growth. A knowledge gap remains regarding the influence of enteric IL-17-IL-17RA signaling and their microbial regulation on the behavior of distant tumors. We demonstrate that gut dysbiosis induced by systemic or gut epithelial deletion of IL-17RA induces growth of pancreatic and brain tumors due to excessive development of Th17, primary source of IL-17 in human and mouse pancreatic ductal adenocarcinoma, as well as B cells that circulate to distant tumors. Microbial dependent IL-17 signaling increases DUOX2 signaling in tumor cells. Inefficacy of pharmacological inhibition of IL-17RA is overcome with targeted microbial ablation that blocks the compensatory loop. These findings demonstrate the complexities of IL-17-IL-17RA signaling in different compartments and the relevance for accounting for its homeostatic host defense function during cancer therapy.

The entrancing role of dietary polyphenols against the most frequent aging-associated diseases

Aging, a fundamental physiological process influenced by innumerable biological and genetic pathways, is an important driving factor for several aging-associated disorders like diabetes mellitus, osteoporosis, cancer, and neurodegenerative diseases including Alzheimer's and Parkinson's diseases. In the modern era, the several mechanisms associated with aging have been deeply studied. Treatment and therapeutics for age-related diseases have also made considerable advances; however, for the effective and long-lasting treatment, nutritional therapy particularly including dietary polyphenols from the natural origin are endorsed. These dietary polyphenols (e.g., apigenin, baicalin, curcumin, epigallocatechin gallate, kaempferol, quercetin, resveratrol, and theaflavin), and many other phytochemicals target certain molecular, genetic mechanisms. The most common pathways of age-associated diseases are mitogen-activated protein kinase, reactive oxygen species production, nuclear factor kappa light chain enhancer of activated B cells signaling pathways, metal chelation, c-Jun N-terminal kinase, and inflammation. Polyphenols slow down the course of aging and help in combatting age-linked disorders. This exemplified in the form of clinical trials on specific dietary polyphenols in various aging-associated diseases. With this context in mind, this review reveals the new insights to slow down the aging process, and consequently reduce some classic diseases associated with age such as aforementioned, and targeting age-associated diseases by the activities of dietary polyphenols of natural origin.

Apolipoprotein L genes are novel mediators of inflammation in beta cells

AIMS/HYPOTHESIS

Inflammation induces beta cell dysfunction and demise but underlying molecular mechanisms remain unclear. The apolipoprotein L (APOL) family of genes has been associated with innate immunity and apoptosis in non-pancreatic cell types, but also with metabolic syndrome and type 2 diabetes mellitus. Here, we hypothesised that APOL genes play a role in inflammation-induced beta cell damage.

METHODS

We used single-cell transcriptomics datasets of primary human pancreatic islet cells to study the expression of APOL genes upon specific stress conditions. Validation of the findings was carried out in EndoC-βH1 cells and primary human islets. Finally, we performed loss- and gain-of-function experiments to investigate the role of APOL genes in beta cells.

RESULTS

APOL genes are expressed in primary human beta cells and APOL1, 2 and 6 are strongly upregulated upon inflammation via the Janus kinase (JAK)-signal transducer and activator of transcription (STAT) pathway. APOL1 overexpression increases endoplasmic reticulum stress while APOL1 knockdown prevents cytokine-induced beta cell death and interferon-associated response. Furthermore, we found that APOL genes are upregulated in beta cells from donors with type 2 diabetes compared with donors without diabetes mellitus.

CONCLUSIONS/INTERPRETATION

APOLs are novel regulators of islet inflammation and may contribute to beta cell damage during the development of diabetes.

DATA AVAILABILITY

scRNAseq data generated by our laboratory and used in this study are available in the Gene Expression Omnibus (GEO; www.ncbi.nlm.nih.gov/geo/ ), accession number GSE218316.

De novo identification of expressed cancer somatic mutations from single-cell RNA sequencing data

Identifying expressed somatic mutations from single-cell RNA sequencing data de novo is challenging but highly valuable. We propose RESA - Recurrently Expressed SNV Analysis, a computational framework to identify expressed somatic mutations from scRNA-seq data. RESA achieves an average precision of 0.77 on three in silico spike-in datasets. In extensive benchmarking against existing methods using 19 datasets, RESA consistently outperforms them. Furthermore, we applied RESA to analyze intratumor mutational heterogeneity in a melanoma drug resistance dataset. By enabling high precision detection of expressed somatic mutations, RESA substantially enhances the reliability of mutational analysis in scRNA-seq. RESA is available at https://github.com/ShenLab-Genomics/RESA .

Cellular age explains variation in age-related cell-to-cell transcriptome variability

Organs and tissues age at different rates within a single individual. Such asynchrony in aging has been widely observed at multiple levels, from functional hallmarks, such as anatomical structures and physiological processes, to molecular endophenotypes, such as the transcriptome and metabolome. However, we lack a conceptual framework to understand why some components age faster than others. Just as demographic models explain why aging evolves, here we test the hypothesis that demographic differences among cell types, determined by cell-specific differences in turnover rate, can explain why the transcriptome shows signs of aging in some cell types but not others. Through analysis of mouse single-cell transcriptome data across diverse tissues and ages, we find that cellular age explains a large proportion of the variation in the age-related increase in transcriptome variance. We further show that long-lived cells are characterized by relatively high expression of genes associated with proteostasis and that the transcriptome of long-lived cells shows greater evolutionary constraint than short-lived cells. In contrast, in short-lived cell types, the transcriptome is enriched for genes associated with DNA repair. Based on these observations, we develop a novel heuristic model that explains how and why aging rates differ among cell types.

Dynamic scRNA-seq of live human pancreatic slices reveals functional endocrine cell neogenesis through an intermediate ducto-acinar stage

Human pancreatic plasticity is implied from multiple single-cell RNA sequencing (scRNA-seq) studies. However, these have been invariably based on static datasets from which fate trajectories can only be inferred using pseudotemporal estimations. Furthermore, the analysis of isolated islets has resulted in a drastic underrepresentation of other cell types, hindering our ability to interrogate exocrine-endocrine interactions. The long-term culture of human pancreatic slices (HPSs) has presented the field with an opportunity to dynamically track tissue plasticity at the single-cell level. Combining datasets from same-donor HPSs at different time points, with or without a known regenerative stimulus (BMP signaling), led to integrated single-cell datasets storing true temporal or treatment-dependent information. This integration revealed population shifts consistent with ductal progenitor activation, blurring of ductal/acinar boundaries, formation of ducto-acinar-endocrine differentiation axes, and detection of transitional insulin-producing cells. This study provides the first longitudinal scRNA-seq analysis of whole human pancreatic tissue, confirming its plasticity in a dynamic fashion.

Smurfness-based two-phase model of ageing helps deconvolve the ageing transcriptional signature

Ageing is characterised at the molecular level by six transcriptional 'hallmarks of ageing', that are commonly described as progressively affected as time passes. By contrast, the 'Smurf' assay separates high-and-constant-mortality risk individuals from healthy, zero-mortality risk individuals, based on increased intestinal permeability. Performing whole body total RNA sequencing, we found that Smurfness distinguishes transcriptional changes associated with chronological age from those associated with biological age. We show that transcriptional heterogeneity increases with chronological age in non-Smurf individuals preceding the other five hallmarks of ageing that are specifically associated with the Smurf state. Using this approach, we also devise targeted pro-longevity genetic interventions delaying entry in the Smurf state. We anticipate that increased attention to the evolutionary conserved Smurf phenotype will bring about significant advances in our understanding of the mechanisms of ageing.

Failure to repair endogenous DNA damage in β-cells causes adult-onset diabetes in mice

Age is the greatest risk factor for the development of type 2 diabetes mellitus (T2DM). Age-related decline in organ function is attributed to the accumulation of stochastic damage, including damage to the nuclear genome. Islets of T2DM patients display increased levels of DNA damage. However, whether this is a cause or consequence of the disease has not been elucidated. Here, we asked if spontaneous, endogenous DNA damage in β-cells can drive β-cell dysfunction and diabetes, via deletion of Ercc1, a key DNA repair gene, in β-cells. Mice harboring Ercc1-deficient β-cells developed adult-onset diabetes as demonstrated by increased random and fasted blood glucose levels, impaired glucose tolerance, and reduced insulin secretion. The inability to repair endogenous DNA damage led to an increase in oxidative DNA damage and apoptosis in β-cells and a significant loss of β-cell mass. Using electron microscopy, we identified β-cells in clear distress that showed an increased cell size, enlarged nuclear size, reduced number of mature insulin granules, and decreased number of mitochondria. Some β-cells were more affected than others consistent with the stochastic nature of spontaneous DNA damage. Ercc1-deficiency in β-cells also resulted in loss of β-cell function as glucose-stimulated insulin secretion and mitochondrial function were impaired in islets isolated from mice harboring Ercc1-deficient β-cells. These data reveal that unrepaired endogenous DNA damage is sufficient to drive β-cell dysfunction and provide a mechanism by which age increases the risk of T2DM.

Measuring cell-to-cell expression variability in single-cell RNA-sequencing data: a comparative analysis and applications to B cell aging

BACKGROUND

Single-cell RNA-sequencing (scRNA-seq) technologies enable the capture of gene expression heterogeneity and consequently facilitate the study of cell-to-cell variability at the cell type level. Although different methods have been proposed to quantify cell-to-cell variability, it is unclear what the optimal statistical approach is, especially in light of challenging data structures that are unique to scRNA-seq data like zero inflation.

RESULTS

We systematically evaluate the performance of 14 different variability metrics that are commonly applied to transcriptomic data for measuring cell-to-cell variability. Leveraging simulations and real datasets, we benchmark the metric performance based on data-specific features, sparsity and sequencing platform, biological properties, and the ability to recapitulate true levels of biological variability based on known gene sets. Next, we use scran, the metric with the strongest all-round performance, to investigate changes in cell-to-cell variability that occur during B cell differentiation and the aging processes. The analysis of primary cell types from hematopoietic stem cells (HSCs) and B lymphopoiesis reveals unique gene signatures with consistent patterns of variable and stable expression profiles during B cell differentiation which highlights the significance of these methods. Identifying differentially variable genes between young and old cells elucidates the regulatory changes that may be overlooked by solely focusing on mean expression changes and we investigate this in the context of regulatory networks.

CONCLUSIONS

We highlight the importance of capturing cell-to-cell gene expression variability in a complex biological process like differentiation and aging and emphasize the value of these findings at the level of individual cell types.

Single-cell metabolic profiling reveals subgroups of primary human hepatocytes with heterogeneous responses to drug challenge

BACKGROUND

Xenobiotics are primarily metabolized by hepatocytes in the liver, and primary human hepatocytes are the gold standard model for the assessment of drug efficacy, safety, and toxicity in the early phases of drug development. Recent advances in single-cell genomics demonstrate liver zonation and ploidy as main drivers of cellular heterogeneity. However, little is known about the impact of hepatocyte specialization on liver function upon metabolic challenge, including hepatic metabolism, detoxification, and protein synthesis.

RESULTS

Here, we investigate the metabolic capacity of individual human hepatocytes in vitro. We assess how chronic accumulation of lipids enhances cellular heterogeneity and impairs the metabolisms of drugs. Using a phenotyping five-probe cocktail, we identify four functional subgroups of hepatocytes responding differently to drug challenge and fatty acid accumulation. These four subgroups display differential gene expression profiles upon cocktail treatment and xenobiotic metabolism-related specialization. Notably, intracellular fat accumulation leads to increased transcriptional variability and diminishes the drug-related metabolic capacity of hepatocytes.

CONCLUSIONS

Our results demonstrate that, upon a metabolic challenge such as exposure to drugs or intracellular fat accumulation, hepatocyte subgroups display different and heterogeneous transcriptional responses.

Reprogramming of pancreatic islet cells for regeneration and rejuvenation

The pancreatic β cell, which produces insulin, is a terminally differentiated cell type that divides rarely. Consequently, the regenerative ability of β cells is limited and irreversible diabetes occurs after severe loss of β-cell function. In view of such poor regenerative capacity, considerable research efforts have been made to promote the expansion of functional insulin-producing cells as a regenerative therapy for diabetes. Here, we discuss recent findings regarding the robust expansion of functional mature islet cells both in vivo and ex vivo through MYCL-mediated reprogramming. We also describe the potential prospects for the application of reprogramming technologies to regenerative therapy and rejuvenation of islet cells.

Decoding Human Biology and Disease Using Single-cell Omics Technologies

Over the past decade, advances in single-cell omics (SCO) technologies have enabled the investigation of cellular heterogeneity at an unprecedented resolution and scale, opening a new avenue for understanding human biology and disease. In this review, we summarize the developments of sequencing-based SCO technologies and computational methods, and focus on considerable insights acquired from SCO sequencing studies to understand normal and diseased properties, with a particular emphasis on cancer research. We also discuss the technological improvements of SCO and its possible contribution to fundamental research of the human, as well as its great potential in clinical diagnoses and personalized therapies of human disease.

IFNγ-Stat1 axis drives aging-associated loss of intestinal tissue homeostasis and regeneration

The influence of aging on intestinal stem cells and their niche can explain underlying causes for perturbation in their function observed during aging. Molecular mechanisms for such a decrease in the functionality of intestinal stem cells during aging remain largely undetermined. Using transcriptome-wide approaches, our study demonstrates that aging intestinal stem cells strongly upregulate antigen presenting pathway genes and over-express secretory lineage marker genes resulting in lineage skewed differentiation into the secretory lineage and strong upregulation of MHC class II antigens in the aged intestinal epithelium. Mechanistically, we identified an increase in proinflammatory cells in the lamina propria as the main source of elevated interferon gamma (IFNγ) in the aged intestine, that leads to the induction of Stat1 activity in intestinal stem cells thus priming the aberrant differentiation and elevated antigen presentation in epithelial cells. Of note, systemic inhibition of IFNγ-signaling completely reverses these aging phenotypes and reinstalls regenerative capacity of the aged intestinal epithelium.

Bioinformatics Analysis and Experimental Verification of Exercise for Aging Mice in Different Brain Regions Based on Transcriptome Sequencing

PURPOSE

Physical exercise mitigates the effects of aging and cognitive decline. However, the precise neurobiological mechanisms underlying this phenomenon remain unclear. The primary aim of this study was to investigate the protective effect of exercise on age-related memory deficits in the prefrontal cortex (PFC) and hippocampus using bioinformatic analysis and biochemical verification.

METHODS

Young and aging mice were subjected to natural feeding or treadmill exercise (12 m/min, 8 weeks). Cognitive function was accessed using the Barnes maze and novel object recognition. Bioinformatic analysis was performed to identify co-expressed genes in different groups and brain regions. The selected genes and pathways were validated using RT-qPCR.

RESULTS

Regular exercise significantly ameliorated age-related cognitive deficits. Four up-regulated targets (Ifi27l2a, Irf7, Oas1b, Ifit1) and one down-regulation (Septin2) were reversed by exercise, demonstrating the underlying mechanisms of cognitive functions induced by aging with exercise in the hippocampus and PFC. The Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses indicated that the NOD-like receptor signaling pathway was inhibited in the neuroinflammation effects of exercise in aging mice in both brain regions.

CONCLUSION

Exercise enhances age-related learning and memory deficits. This beneficial effect may be attributed to the changes in five up/down-regulated genes and the NOD-like receptor signaling pathway in both the hippocampus and PFC. These findings establish the modulation of neuroinflammation as a pivotal molecular mechanism supporting exercise intervention in the brain aging process.

GOWDL: gene ontology-driven wide and deep learning model for cell typing of scRNA-seq data

Single-cell RNA-sequencing (scRNA-seq) allows for obtaining genomic and transcriptomic profiles of individual cells. That data make it possible to characterize tissues at the cell level. In this context, one of the main analyses exploiting scRNA-seq data is identifying the cell types within tissue to estimate the quantitative composition of cell populations. Due to the massive amount of available scRNA-seq data, automatic classification approaches for cell typing, based on the most recent deep learning technology, are needed. Here, we present the gene ontology-driven wide and deep learning (GOWDL) model for classifying cell types in several tissues. GOWDL implements a hybrid architecture that considers the functional annotations found in Gene Ontology and the marker genes typical of specific cell types. We performed cross-validation and independent external testing, comparing our algorithm with 12 other state-of-the-art predictors. Classification scores demonstrated that GOWDL reached the best results over five different tissues, except for recall, where we got about 92% versus 97% of the best tool. Finally, we presented a case study on classifying immune cell populations in breast cancer using a hierarchical approach based on GOWDL.

Caloric restriction promotes beta cell longevity and delays aging and senescence by enhancing cell identity and homeostasis mechanisms

Caloric restriction (CR) extends organismal lifespan and health span by improving glucose homeostasis mechanisms. How CR affects organellar structure and function of pancreatic beta cells over the lifetime of the animal remains unknown. Here, we used single nucleus transcriptomics to show that CR increases the expression of genes for beta cell identity, protein processing, and organelle homeostasis. Gene regulatory network analysis link this transcriptional phenotype to transcription factors involved in beta cell identity (Mafa) and homeostasis (Atf6). Imaging metabolomics further demonstrates that CR beta cells are more energetically competent. In fact, high-resolution light and electron microscopy indicates that CR reduces beta cell mitophagy and increases mitochondria mass, increasing mitochondrial ATP generation. Finally, we show that long-term CR delays the onset of beta cell aging and senescence to promote longevity by reducing beta cell turnover. Therefore, CR could be a feasible approach to preserve compromised beta cells during aging and diabetes.

Delineating mouse β-cell identity during lifetime and in diabetes with a single cell atlas

Although multiple pancreatic islet single-cell RNA-sequencing (scRNA-seq) datasets have been generated, a consensus on pancreatic cell states in development, homeostasis and diabetes as well as the value of preclinical animal models is missing. Here, we present an scRNA-seq cross-condition mouse islet atlas (MIA), a curated resource for interactive exploration and computational querying. We integrate over 300,000 cells from nine scRNA-seq datasets consisting of 56 samples, varying in age, sex and diabetes models, including an autoimmune type 1 diabetes model (NOD), a glucotoxicity/lipotoxicity type 2 diabetes model (db/db) and a chemical streptozotocin β-cell ablation model. The β-cell landscape of MIA reveals new cell states during disease progression and cross-publication differences between previously suggested marker genes. We show that β-cells in the streptozotocin model transcriptionally correlate with those in human type 2 diabetes and mouse db/db models, but are less similar to human type 1 diabetes and mouse NOD β-cells. We also report pathways that are shared between β-cells in immature, aged and diabetes models. MIA enables a comprehensive analysis of β-cell responses to different stressors, providing a roadmap for the understanding of β-cell plasticity, compensation and demise.