Single-cell and spatial transcriptomics analysis of human adrenal aging.

OBJECTIVE: The human adrenal cortex comprises three functionally and structurally distinct layers that produce layer-specific steroid hormones. With aging, the human adrenal cortex undergoes functional and structural alteration or "adrenal aging", leading to the unbalanced production of steroid hormones. Given the marked species differences in adrenal biology, the underlying mechanisms of human adrenal aging have not been sufficiently studied. This study was designed to elucidate the mechanisms linking the functional and structural alterations of the human adrenal cortex. METHODS: We conducted single-cell RNA sequencing and spatial transcriptomics analysis of the aged human adrenal cortex. RESULTS: The data of this study suggest that the layer-specific alterations of multiple signaling pathways underlie the abnormal layered structure and layer-specific changes in steroidogenic cells. We also highlighted that macrophages mediate age-related adrenocortical cell inflammation and senescence. CONCLUSIONS: This study is the first detailed analysis of the aged human adrenal cortex at single-cell resolution and helps to elucidate the mechanism of human adrenal aging, thereby leading to a better understanding of the pathophysiology of age-related disorders associated with adrenal aging.

Steroids-producing nodules: a two-layered adrenocortical nodular structure as a precursor lesion of cortisol-producing adenoma.

BACKGROUND: The human adrenal cortex consists of three functionally and structurally distinct layers; zona glomerulosa, zona fasciculata (zF), and zona reticularis (zR), and produces adrenal steroid hormones in a layer-specific manner; aldosterone, cortisol, and adrenal androgens, respectively. Cortisol-producing adenomas (CPAs) occur mostly as a result of somatic mutations associated with the protein kinase A pathway. However, how CPAs develop after adrenocortical cells acquire genetic mutations, remains poorly understood. METHODS: We conducted integrated approaches combining the detailed histopathologic studies with genetic, RNA-sequencing, and spatially resolved transcriptome (SRT) analyses for the adrenal cortices adjacent to human adrenocortical tumours. FINDINGS: Histopathological analysis revealed an adrenocortical nodular structure that exhibits the two-layered zF- and zR-like structure. The nodular structures harbour GNAS somatic mutations, known as a driver mutation of CPAs, and confer cell proliferative and autonomous steroidogenic capacities, which we termed steroids-producing nodules (SPNs). RNA-sequencing coupled with SRT analysis suggests that the expansion of the zF-like structure contributes to the formation of CPAs, whereas the zR-like structure is characterised by a macrophage-mediated immune response. INTERPRETATION: We postulate that CPAs arise from a precursor lesion, SPNs, where two distinct cell populations might contribute differently to adrenocortical tumorigenesis. Our data also provide clues to the molecular mechanisms underlying the layered structures of human adrenocortical tissues. FUNDING: KAKENHI, The Uehara Memorial Foundation, Daiwa Securities Health Foundation, Kaibara Morikazu Medical Science Promotion Foundation, Secom Science and Technology Foundation, ONO Medical Research Foundation, and Japan Foundation for Applied Enzymology.

Downregulation of LAMB3 Altered the Carcinogenic Properties of Human Papillomavirus 16-Positive Cervical Cancer Cells.

Nearly all cervical cancer cases are caused by infection with high-risk human papillomavirus (HR-HPV) types. The mechanism of cervical cell transformation is related to the powerful action of viral oncoproteins and cellular gene alterations. Transcriptomic data from cervical cancer and normal cervical cells were utilized to identify upregulated genes and their associated pathways. The laminin subunit beta-3 (LAMB3) mRNAwas overexpressed in cervical cancer and was chosen for functional analysis. The LAMB3 was predominantly expressed in the extracellular region and the plasma membrane, which play a role in protein binding and cell adhesion molecule binding, leading to cell migration and tissue development. LAMB3 was found to be implicated in the pathway in cancer and the PI3K-AKT signaling pathway. LAMB3 knockdown decreased cell migration, invasion, anchorage-dependent and anchorage-independent cell growth and increased the number of apoptotic cells. These effects were linked to a decrease in protein levels involved in the PI3K-AKT signaling pathway and an increase in p53 protein. This study demonstrated that LAMB3 could promote cervical cancer cell migration, invasion and survival.

Mechanisms of action and resistance in histone methylation-targeted therapy.

Epigenomes enable the rectification of disordered cancer gene expression, thereby providing new targets for pharmacological interventions. The clinical utility of targeting histone H3 lysine trimethylation (H3K27me3) as an epigenetic hallmark has been demonstrated1-7. However, in actual therapeutic settings, the mechanism by which H3K27me3-targeting therapies exert their effects and the response of tumour cells remain unclear. Here we show the potency and mechanisms of action and resistance of the EZH1-EZH2 dual inhibitor valemetostat in clinical trials of patients with adult T cell leukaemia/lymphoma. Administration of valemetostat reduced tumour size and demonstrated durable clinical response in aggressive lymphomas with multiple genetic mutations. Integrative single-cell analyses showed that valemetostat abolishes the highly condensed chromatin structure formed by the plastic H3K27me3 and neutralizes multiple gene loci, including tumour suppressor genes. Nevertheless, subsequent long-term treatment encounters the emergence of resistant clones with reconstructed aggregate chromatin that closely resemble the pre-dose state. Acquired mutations at the PRC2-compound interface result in the propagation of clones with increased H3K27me3 expression. In patients free of PRC2 mutations, TET2 mutation or elevated DNMT3A expression causes similar chromatin recondensation through de novo DNA methylation in the H3K27me3-associated regions. We identified subpopulations with distinct metabolic and gene translation characteristics implicated in primary susceptibility until the acquisition of the heritable (epi)mutations. Targeting epigenetic drivers and chromatin homeostasis may provide opportunities for further sustained epigenetic cancer therapies.

An improved method for the highly specific detection of transcription start sites.

Precise detection of the transcriptional start site (TSS) is a key for characterizing transcriptional regulation of genes and for annotation of newly sequenced genomes. Here, we describe the development of an improved method, designated 'TSS-seq2.' This method is an iterative improvement of TSS-seq, a previously published enzymatic cap-structure conversion method to detect TSSs in base sequences. By modifying the original procedure, including by introducing split ligation at the key cap-selection step, the yield and the accuracy of the reaction has been substantially improved. For example, TSS-seq2 can be conducted using as little as 5 ng of total RNA with an overall accuracy of 96%; this yield a less-biased and more precise detection of TSS. We then applied TSS-seq2 for TSS analysis of four plant species that had not yet been analyzed by any previous TSS method.

Mint3-depletion-induced energy stress sensitizes triple-negative breast cancer to chemotherapy via HSF1 inactivation.

Given the lack of therapeutic targets, the conventional approach for managing triple-negative breast cancer (TNBC) involves the utilization of cytotoxic chemotherapeutic agents. However, most TNBCs acquire resistance to chemotherapy, thereby lowering the therapeutic outcome. In addition to oncogenic mutations in TNBC, microenvironment-induced mechanisms render chemoresistance more complex and robust in vivo. Here, we aimed to analyze whether depletion of Munc18-1 interacting protein 3 (Mint3), which activates hypoxia-inducible factor 1 (HIF-1) during normoxia, sensitizes TNBC to chemotherapy. We found that Mint3 promotes the chemoresistance of TNBC in vivo. Mint3 depletion did not affect the sensitivity of human TNBC cell lines to doxorubicin and paclitaxel in vitro but sensitized tumors of these cells to chemotherapy in vivo. Transcriptome analyses revealed that the Mint3-HIF-1 axis enhanced heat shock protein 70 (HSP70) expression in tumors of TNBC cells. Administering an HSP70 inhibitor enhanced the antitumor activity of doxorubicin in TNBC tumors, similar to Mint3 depletion. Mint3 expression was also correlated with HSP70 expression in human TNBC specimens. Mechanistically, Mint3 depletion induces glycolytic maladaptation to the tumor microenvironment in TNBC tumors, resulting in energy stress. This energy stress by Mint3 depletion inactivated heat shock factor 1 (HSF-1), the master regulator of HSP expression, via the AMP-activated protein kinase/mechanistic target of the rapamycin pathway following attenuated HSP70 expression. In conclusion, Mint3 is a unique regulator of TNBC chemoresistance in vivo via metabolic adaptation to the tumor microenvironment, and a combination of Mint3 inhibition and chemotherapy may be a good strategy for TNBC treatment.

Whole-genome sequencing reveals the molecular implications of the stepwise progression of lung adenocarcinoma.

The mechanism underlying the development of tumors, particularly at early stages, still remains mostly elusive. Here, we report whole-genome long and short read sequencing analysis of 76 lung cancers, focusing on very early-stage lung adenocarcinomas such as adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma. The obtained data is further integrated with bulk and spatial transcriptomic data and epigenomic data. These analyses reveal key events in lung carcinogenesis. Minimal somatic mutations in pivotal driver mutations and essential proliferative factors are the only detectable somatic mutations in the very early-stage of AIS. These initial events are followed by copy number changes and global DNA hypomethylation. Particularly, drastic changes are initiated at the later AIS stage, i.e., in Noguchi type B tumors, wherein cancer cells are exposed to the surrounding microenvironment. This study sheds light on the pathogenesis of lung adenocarcinoma from integrated pathological and molecular viewpoints.

The dataset of de novo assembly and inferred functional annotation of the transcriptome of Heterosigma akashiwo, a bloom-forming, cosmopolitan raphidophyte.

Heterosigma akashiwo is a eukaryotic, cosmopolitan, and unicellular alga (class: Raphidophyceae), and produces fish-killing blooms. There is a substantial scientific and practical interest in its ecophysiological characteristics that determine bloom dynamics and its adaptation to broad climate zones. A well-annotated genomic/genetic sequence information enables researchers to characterize organisms using modern molecular technology. In the present study, we conducted H. akashiwo RNA sequencing, a de novo transcriptome assembly of 84,693,530 high-quality deduplicated short-read sequences. Obtained RNA reads were assembled by Trinity assembler and 144,777 contigs were identified with N50 values of 1085. Total 60,877 open reading frames with the length of 150 bp or greater were predicted. For further analyses, top Gene Ontology terms, pfam hits, and blast hits were annotated for all the predicted genes. The raw data were deposited in the NCBI SRA database (BioProject PRJDB6241 and PRJDB15108), and the assemblies are available in NCBI TSA database (ICRV01). The annotation information can be obtained in Dryad and can be accessed via doi: 10.5061/dryad.m0cfxpp56.

Stem cell-based modeling and single-cell multiomics reveal gene-regulatory mechanisms underlying human skeletal development.

Although the skeleton is essential for locomotion, endocrine functions, and hematopoiesis, the molecular mechanisms of human skeletal development remain to be elucidated. Here, we introduce an integrative method to model human skeletal development by combining in vitro sclerotome induction from human pluripotent stem cells and in vivo endochondral bone formation by implanting the sclerotome beneath the renal capsules of immunodeficient mice. Histological and scRNA-seq analyses reveal that the induced bones recapitulate endochondral ossification and are composed of human skeletal cells and mouse circulatory cells. The skeletal cell types and their trajectories are similar to those of human embryos. Single-cell multiome analysis reveals dynamic changes in chromatin accessibility associated with multiple transcription factors constituting cell-type-specific gene-regulatory networks (GRNs). We further identify ZEB2, which may regulate the GRNs in human osteogenesis. Collectively, these results identify components of GRNs in human skeletal development and provide a valuable model for its investigation.

Direct Reprogramming Improves Cardiac Function and Reverses Fibrosis in Chronic Myocardial Infarction.

BACKGROUND: Because adult cardiomyocytes have little regenerative capacity, resident cardiac fibroblasts (CFs) synthesize extracellular matrix after myocardial infarction (MI) to form fibrosis, leading to cardiac dysfunction and heart failure. Therapies that can regenerate the myocardium and reverse fibrosis in chronic MI are lacking. The overexpression of cardiac transcription factors, including Mef2c/Gata4/Tbx5/Hand2 (MGTH), can directly reprogram CFs into induced cardiomyocytes (iCMs) and improve cardiac function under acute MI. However, the ability of in vivo cardiac reprogramming to repair chronic MI with established scars is undetermined. METHODS: We generated a novel Tcf21iCre/reporter/MGTH2A transgenic mouse system in which tamoxifen treatment could induce both MGTH and reporter expression in the resident CFs for cardiac reprogramming and fibroblast lineage tracing. We first tested the efficacy of this transgenic system in vitro and in vivo for acute MI. Next, we analyzed in vivo cardiac reprogramming and fusion events under chronic MI using Tcf21iCre/Tomato/MGTH2A and Tcf21iCre/mTmG/MGTH2A mice, respectively. Microarray and single-cell RNA sequencing were performed to determine the mechanism of cardiac repair by in vivo reprogramming. RESULTS: We confirmed the efficacy of transgenic in vitro and in vivo cardiac reprogramming for acute MI. In chronic MI, in vivo cardiac reprogramming converted ≈2% of resident CFs into iCMs, in which a majority of iCMs were generated by means of bona fide cardiac reprogramming rather than by fusion with cardiomyocytes. Cardiac reprogramming significantly improved myocardial contraction and reduced fibrosis in chronic MI. Microarray analyses revealed that the overexpression of MGTH activated cardiac program and concomitantly suppressed fibroblast and inflammatory signatures in chronic MI. Single-cell RNA sequencing demonstrated that resident CFs consisted of 7 subclusters, in which the profibrotic CF population increased under chronic MI. Cardiac reprogramming suppressed fibroblastic gene expression in chronic MI by means of conversion of profibrotic CFs to a quiescent antifibrotic state. MGTH overexpression induced antifibrotic effects partly by suppression of Meox1, a central regulator of fibroblast activation. CONCLUSIONS: These results demonstrate that cardiac reprogramming could repair chronic MI by means of myocardial regeneration and reduction of fibrosis. These findings present opportunities for the development of new therapies for chronic MI and heart failure.

Establishment of experimental salivary gland cancer models using organoid culture and patient-derived xenografting.

PURPOSE: Depending on its histological subtype, salivary gland carcinoma (SGC) may have a poor prognosis. Due to the scarcity of preclinical experimental models, its molecular biology has so far remained largely unknown, hampering the development of new treatment modalities for patients with these malignancies. The aim of this study was to generate experimental human SGC models of multiple histological subtypes using patient-derived xenograft (PDX) and organoid culture techniques. METHODS: Tumor specimens from surgically resected SGCs were processed for the preparation of PDXs and patient-derived organoids (PDOs). Specimens from SGC PDXs were also processed for PDX-derived organoid (PDXO) generation. In vivo tumorigenicity was assessed using orthotopic transplantation of SGC organoids. The pathological characteristics of each model were compared to those of the original tumors using immunohistochemistry. RNA-seq was used to analyze the genetic traits of our models. RESULTS: Three series of PDOs, PDXs and PDXOs of salivary duct carcinomas, one series of PDOs, PDXs and PDXOs of mucoepidermoid carcinomas and PDXs of myoepithelial carcinomas were successfully generated. We found that PDXs and orthotopic transplants from PDOs/PDXOs showed similar histological features as the original tumors. Our models also retained their genetic traits, i.e., transcription profiles, genomic variants and fusion genes of the corresponding histological subtypes. CONCLUSION: We report the generation of SGC PDOs, PDXs and PDXOs of multiple histological subtypes, recapitulating the histological and genetical characteristics of the original tumors. These experimental SGC models may serve as a useful resource for the development of novel therapeutic strategies and for investigating the molecular mechanisms underlying the development of these malignancies.

UTX deficiency in neural stem/progenitor cells results in impaired neural development, fetal ventriculomegaly, and postnatal death.

Recent studies have demonstrated that epigenetic modifications are deeply involved in neurogenesis; however, the precise mechanisms remain largely unknown. To determine the role of UTX (also known as KDM6A), a demethylase of histone H3K27, in neural development, we generated Utx-deficient mice in neural stem/progenitor cells (NSPCs). Since Utx is an X chromosome-specific gene, the genotypes are sex-dependent; female mice lose both Utx alleles (UtxΔ/Δ ), and male mice lose one Utx allele yet retain one Uty allele, the counterpart of Utx on the Y chromosome (UtxΔ/Uty ). We found that UtxΔ/Δ mice exhibited fetal ventriculomegaly and died soon after birth. Immunofluorescence staining and EdU labeling revealed a significant increase in NSPCs and a significant decrease in intermediate-progenitor and differentiated neural cells. Molecular analyses revealed the downregulation of pathways related to DNA replication and increased H3K27me3 levels around the transcription start sites in UtxΔ/Δ NSPCs. These results indicate that UTX globally regulates the expression of genes required for proper neural development in NSPCs, and UTX deficiency leads to impaired cell cycle exit, reduced differentiation, and neonatal death. Interestingly, although UtxΔ/Uty mice survived the postnatal period, most died of hydrocephalus, a clinical feature of Kabuki syndrome, a congenital anomaly involving UTX mutations. Our findings provide novel insights into the role of histone modifiers in neural development and suggest that UtxΔ/Uty mice are a potential disease model for Kabuki syndrome.

Runx2 regulates chromatin accessibility to direct the osteoblast program at neonatal stages.

The transcriptional regulator Runx2 (runt-related transcription factor 2) has essential but distinct roles in osteoblasts and chondrocytes in skeletal development. However, Runx2-mediated regulatory mechanisms underlying the distinctive programming of osteoblasts and chondrocytes are not well understood. Here, we perform an integrative analysis to investigate Runx2-DNA binding and chromatin accessibility ex vivo using neonatal osteoblasts and chondrocytes. We find that Runx2 engages with cell-type-distinct chromatin-accessible regions, potentially interacting with different combinations of transcriptional regulators, forming cell-type-specific hotspots, and potentiating chromatin accessibility. Genetic analysis and direct cellular reprogramming studies suggest that Runx2 is essential for establishment of chromatin accessibility in osteoblasts. Functional enhancer studies identify an Sp7 distal enhancer driven by Runx2-dependent binding and osteoblast-specific chromatin accessibility, contributing to normal osteoblast differentiation. Our findings provide a framework for understanding the regulatory landscape encompassing Runx2-mediated and cell-type-distinct enhancer networks that underlie the specification of osteoblasts.

Unraveling unique features of plasma cell clones in POEMS syndrome with single-cell analysis.

POEMS syndrome is a rare monoclonal plasma cell disorder, with unique symptoms distinct from those of other plasma cell neoplasms, including high serum VEGF levels. Because the prospective isolation of POEMS clones has not yet been successful, their real nature remains unclear. Herein, we performed single-cell RNA-Seq of BM plasma cells from patients with POEMS syndrome and identified POEMS clones that had Ig λ light chain (IGL) sequences (IGLV1-36, -40, -44, and -47) with amino acid changes specific to POEMS syndrome. The proportions of POEMS clones in plasma cells were markedly smaller than in patients with multiple myeloma (MM) and patients with monoclonal gammopathy of undetermined significance (MGUS). Single-cell transcriptomes revealed that POEMS clones were CD19+, CD138+, and MHC class IIlo, which allowed for their prospective isolation. POEMS clones expressed significantly lower levels of c-MYC and CCND1 than MM clones, accounting for their small size. VEGF mRNA was not upregulated in POEMS clones, directly indicating that VEGF is not produced by POEMS clones. These results reveal unique features of POEMS clones and enhance our understanding of the pathogenesis of POEMS syndrome.

Semibulk RNA-seq analysis as a convenient method for measuring gene expression statuses in a local cellular environment.

When biologically interpretation of the data obtained from the single-cell RNA sequencing (scRNA-seq) analysis is attempted, additional information on the location of the single cells, behavior of the surrounding cells, and the microenvironment they generate, would be very important. We developed an inexpensive, high throughput application while preserving spatial organization, named "semibulk RNA-seq" (sbRNA-seq). We utilized a microfluidic device specifically designed for the experiments to encapsulate both a barcoded bead and a cell aggregate (a semibulk) into a single droplet. Using sbRNA-seq, we firstly analyzed mouse kidney specimens. In the mouse model, we could associate the pathological information with the gene expression information. We validated the results using spatial transcriptome analysis and found them highly consistent. When we applied the sbRNA-seq analysis to the human breast cancer specimens, we identified spatial interactions between a particular population of immune cells and that of cancer-associated fibroblast cells, which were not precisely represented solely by the single-cell analysis. Semibulk analysis may provide a convenient and versatile method, compared to a standard spatial transcriptome sequencing platform, to associate spatial information with transcriptome information.

A bifurcation concept for B-lymphoid/plasmacytoid dendritic cells with largely fluctuating transcriptome dynamics.

Hematopoiesis was considered a hierarchical stepwise process but was revised to a continuous process following single-cell RNA sequencing. However, the uncertainty or fluctuation of single-cell transcriptome dynamics during differentiation was not considered, and the dendritic cell (DC) pathway in the lymphoid context remains unclear. Here, we identify human B-plasmacytoid DC (pDC) bifurcation as large fluctuating transcriptome dynamics in the putative B/NK progenitor region by dry and wet methods. By converting splicing kinetics into diffusion dynamics in a deep generative model, our original computational methodology reveals strong fluctuation at B/pDC bifurcation in IL-7Rα+ regions, and LFA-1 fluctuates positively in the pDC direction at the bifurcation. These expectancies are validated by the presence of B/pDC progenitors in the IL-7Rα+ fraction and preferential expression of LFA-1 in pDC-biased progenitors with a niche-like culture system. We provide a model of fluctuation-based differentiation, which reconciles continuous and discrete models and is applicable to other developmental systems.

RSPO2 defines a distinct undifferentiated progenitor in the tendon/ligament and suppresses ectopic ossification.

Ectopic endochondral ossification in the tendon/ligament is caused by repetitive mechanical overload or inflammation. Tendon stem/progenitor cells (TSPCs) contribute to tissue repair, and some express lubricin [proteoglycan 4 (PRG4)]. However, the mechanisms of ectopic ossification and association of TSPCs are not yet known. Here, we investigated the characteristics of Prg4-positive (+) cells and identified that R-spondin 2 (RSPO2), a WNT activator, is specifically expressed in a distinct Prg4+ TSPC cluster. The Rspo2+ cluster was characterized as mostly undifferentiated, and RSPO2 overexpression suppressed ectopic ossification in a mouse Achilles tendon puncture model via chondrogenic differentiation suppression. RSPO2 expression levels in patients with ossification of the posterior longitudinal ligament were lower than those in spondylosis patients, and RSPO2 protein suppressed chondrogenic differentiation of human ligament cells. RSPO2 was induced by inflammatory stimulation and mechanical loading via nuclear factor κB. Rspo2+ cells may contribute to tendon/ligament homeostasis under pathogenic conditions.

Single mutation makes Escherichia coli an insect mutualist.

Microorganisms often live in symbiosis with their hosts, and some are considered mutualists, where all species involved benefit from the interaction. How free-living microorganisms have evolved to become mutualists is unclear. Here we report an experimental system in which non-symbiotic Escherichia coli evolves into an insect mutualist. The stinkbug Plautia stali is typically associated with its essential gut symbiont, Pantoea sp., which colonizes a specialized symbiotic organ. When sterilized newborn nymphs were infected with E. coli rather than Pantoea sp., only a few insects survived, in which E. coli exhibited specific localization to the symbiotic organ and vertical transmission to the offspring. Through transgenerational maintenance with P. stali, several hypermutating E. coli lines independently evolved to support the host's high adult emergence and improved body colour; these were called 'mutualistic' E. coli. These mutants exhibited slower bacterial growth, smaller size, loss of flagellar motility and lack of an extracellular matrix. Transcriptomic and genomic analyses of 'mutualistic' E. coli lines revealed independent mutations that disrupted the carbon catabolite repression global transcriptional regulator system. Each mutation reproduced the mutualistic phenotypes when introduced into wild-type E. coli, confirming that single carbon catabolite repression mutations can make E. coli an insect mutualist. These findings provide an experimental system for future work on host-microbe symbioses and may explain why microbial mutualisms are omnipresent in nature.

Characterization of Aldosterone-producing Cell Cluster (APCC) at Single-cell Resolution.

CONTEXT: The adrenal cortex consists of zona glomerulosa (ZG), fasciculata (ZF), and reticularis. Aldosterone-producing cell clusters (APCCs) that strongly express aldosterone synthase (CYP11B2) are frequently found in adult adrenals and harbor somatic mutations that are also detected in aldosterone-producing adenomas (APAs). Primary aldosteronism is mainly caused by APAs or idiopathic hyperaldosteronism (IHA). We presume that APCCs are causing IHA and are precursors of APAs. However, the gene expression characteristics and especially the development of APCCs are not well understood. OBJECTIVE: This study aimed to analyze the transcriptome of APCCs at single-cell resolution and infer the developmental trajectory. METHODS: Single-cell RNA sequencing (scRNA-seq) of 2 adult adrenals was performed. RESULTS: Immunohistochemical analyses confirmed the 2 adrenals had APCCs. scRNA-seq data of 2928 adrenal cells were obtained and 1765 adrenocortical cells were identified based on unsupervised clustering and the marker gene expression. The adrenocortical cells were divided into 6 clusters, of which 3 clusters (923 cells) were composed of APCC/ZG cells. By further subclustering, the APCC/ZG cells were divided into 3 clusters (clusters C1, C2, and C3), we finally identified APCC cluster (C3) and ZG cluster (C1). Cluster C2 seemed to be ZG-to-ZF transitional cells. RNA velocity analysis inferred the developmental direction from cluster ZG-cluster-C1 to APCC-cluster-C3. The scRNA-seq additionally revealed that many CYP11B2-positive cells were positive for CYP11B1 and/or CYP17A1, which were essential for cortisol but not for aldosterone production. CONCLUSIONS: Our results revealed the gene expression characteristics of APCC at single-cell resolution and show that some ZG cells remodel to APCC.