ATAC Sequencing Research Articles

IKZF2 Specifically Participates in AML1-ETO Mediated Blockage of Myeloid Differentiation

Objectives:The chromosomal translocation (8;21) is one of the most recurrent cytogenetic aberrations in acute myeloid leukemia (AML). Although most patients achieve complete remission, one third still relapse due to the persistent leukemia stem cells (LSCs). The lack of significant improvements in outcomes, especially in elderly patients ineligible for standard chemotherapy, has spurred further investigation in t(8;21) leukemia. In this study, we investigate the mechanism of AML1-ETO (AE) fusion gene induced differentiation blockage from earlier stages of HSPCs, aiming to identify potential therapeutic targets that can effectively overcome treatment resistance and minimize long-term side effects. Methods:The AML1-ETO conditional knock-in mouse (Aml1 Eto/+; Mx1-Cre) was previously established, hereafter called AE KI. The Cas9-expressing AE KI model, generated by crossing AE KI and Cas9 mice strains, was employed to define the regulatory program of AE fusion gene. For knockout of AE and IKZF2, Kasumi-1 cells were transduced with the single-vector lentiviral system containing sgRNAs for the respective genes, or a nontargeting control. Colony formation and serial bone marrow transplantation assay were performed to evaluate the function of HSPCs. ATAC sequencing, transcriptional profiling (RNA-seq) and Cut&Tag for H3K4me3, H3K27me3, H3K27ac were performed and analyzed comprehensively to characterize the regulatory profile of AE fusion gene. Results:Our previous research demonstrated that induction expression of AE fusion gene in vivo resulted in excessive accumulation of Lin -Sca-1 +c-Kit + (LSK) cells, leading to differentiation blockage from HSCs to HPCs. IKZF2 showed significantly higher expression in LSK cells following PIPC induction of AE expression compared to that of the control group. IKZF2 was highly increased in AML patients harboring AE fusion gene. And the level of IKZF2 was significantly higher in human t(8;21) cell lines Kasumi-1 and SKNO-1 when compared to that of non-AE controls. Depletion of the AE fusion gene led to a dramatic decrease in IKZF2 expression, while constitutive expression of the AE fusion protein in U-937 cells resulted in elevated levels of IKZF2. Loss of IKZF2 resulted in apoptosis in AE-expressing leukemia cells, and led to an increase in the myeloid differentiation markers CD11b and CD14. These findings suggested a functional relevance of IKZF2 in AE-driven AML. In our Cas9-expressing AE KI model, a significant improvement in the peripheral chimerism rate was observed following IKZF2 knockout at various time points. The population of HSCs, mainly the ST-HSCs, was reduced, while the more differentiated MPP and Lin -Sca-1 -c-Kit + (LK) pools, especially the GMP and CMP groups, were expanded. These findings indicated that the differentiation block triggered by AE was partially relieved by IKZF2 knockout. Furthermore, the deletion of IKZF2 had minimal impact on normal myeloid differentiation, highlighting its specific role in AE settings. We further investigated the underlying mechanisms through multi-omics analysis. The induction of AE led to increased accessibility of the myeloid differentiation inhibitory elements, including IKZF2. CUT&Tag data confirmed a significant decrease in H3K27me3 signals at the promoter regions of IKZF2 compared to that of the control group. Regarding the role of IKZF2 in chromatin remodeling, it is likely to interact with AE to form an epigenetic-transcriptional control module. This may expand the signal for differentiation blockage and drive critical leukemogenic gene expression profiles. Conclusions: Our results demonstrate that IKZF2 is a target gene of AML1-ETO crucial for AML1-ETO-induced differentiation blockage from HSCs to HPCs. The epigenetic regulation mechanism may contribute to the AML1-ETO-IKZF2 regulatory axis. And targeting IKZF2 may become a potential therapeutic strategy in human t(8;21) AML especially for the eradication of AML1-ETO-associated LSCs.

A Bayesian method to infer copy number clones from single-cell RNA and ATAC sequencing.

Single-cell RNA and ATAC sequencing technologies enable the examination of gene expression and chromatin accessibility in individual cells, providing insights into cellular phenotypes. In cancer research, it is important to consistently analyze these states within an evolutionary context on genetic clones. Here we present CONGAS+, a Bayesian model to map single-cell RNA and ATAC profiles onto the latent space of copy number clones. CONGAS+ clusters cells into tumour subclones with similar ploidy, rendering straightforward to compare their expression and chromatin profiles. The framework, implemented on GPU and tested on real and simulated data, scales to analyse seamlessly thousands of cells, demonstrating better performance than single-molecule models, and supporting new multi-omics assays. In prostate cancer, lymphoma and basal cell carcinoma, CONGAS+ successfully identifies complex subclonal architectures while providing a coherent mapping between ATAC and RNA, facilitating the study of genotype-phenotype maps and their connection to genomic instability.

Epigenetic Classification of Acute Myeloid Leukemia Revealed By Genome-Wide Chromatin Profiling

Background Acute myeloid leukemia (AML) is a clinically and molecularly heterogeneous disease and is predominantly defined by genetic abnormalities and morphology. However, many lines of evidence suggest the role of epigenetic alterations in the pathogenesis and the classification of AML, whereas no large-scale comprehensive analyses of both genetic and epigenetic alterations in AML have been reported. Methods We enrolled a total of 1,564 primary AML cases (524 from Japanese and 1,040 from Swedish cohorts), whose diagnostic samples were subjected to ATAC sequencing together with targeted capture sequencing using an in-house gene panel including 331 known AML driver genes and an additional 1,158-1,317 probes for copy number detection. We also performed RNA (N=1,398) and H3K27ac-targeted ChIP sequencing (N=80). Normal bone marrow samples from 25 donors as well as 15 remission samples were analyzed by ATAC-seq as a control. Results ATAC-seq revealed approximately 300,000 reproducible peaks in primary AML. Most peaks were found in gene-distal regions, and explained a large variance across samples. Among all peaks, 72% were not detected in control samples and thus, specific to AML. Analysis of the correlations between ATAC-seq accessibility and gene expression revealed a subset of regulatory elements that regulate nearby gene expression; each peak was linked to a median of 1 gene (0-31), while each gene was associated with a median of 11 peaks (0-110). We next estimated the cellular composition of the bulk AML samples by deconvolution analysis using ATAC-seq data. The predominance of hematopoietic stem cell (HSC), monocyte, and erythroid signatures was associated with French-American-British (FAB) subtypes of minimal differentiation/without maturation subtypes (M0/M1), myelomonocytic/monoblastic/monocytic subtypes (M4/M5), and erythroid subtype (M6), respectively, suggesting bulk ATAC-seq of AML samples reflected the differentiation potential of leukemic cells. The ATAC-based clustering identified 11 distinct epigenetic subgroups, including three well-known genetic classes defined by gene fusions, such as PML::RARA, RUNX1::RUNX1T1, and CBFB::MYH11. By contrast, some genetic classes defined by a single mutation, such as NPM1 or CEBPA mutation, were further classified into multiple distinct epigenetic subgroups. For example, the NPM1-mutated subtype was separated into four subgroups based on the varying degrees of ATAC signatures associated with HSC and monocytes. Each subgroup of NPM1-mutated AML showed different clinical characteristics such as white blood cell counts, age, and prognosis, and co-mutation pattern in FLT3, DNMT3A, and IDH1/ IDH2 mutations. Gene ontology analysis revealed that the subgroup of NPM1-mutated AML with high HSC signature showed upregulated TGF-β signaling, while one with monocytic signature showed upregulated TNF-α signaling and interferon-γ response. Next, we evaluated transcription factor (TF) binding sites in a genome-wide manner, using TF motifs and foot printing in upregulated ATAC-peaks in each subgroup and revealed distinct profiles of activated TFs in each epigenetic subgroup of NPM1-mutated AML. In the subgroup characterized by high HSC signature, HSC-related TFs that belong to MECOM, HOX, and RUNX families, were more active. By contrast, other TFs, such as IRF, CEBP, and JUN/FOS families, were more enriched in the subgroup showing high monocyte signature. By developing a prediction model for these subgroups of NPM1-mutated AML based on the gene expression profile, we reproduced these subgroups and validated clinical and molecular features in an external Beat AML cohort. Utilizing ex vivo drug sensitivity dataset in the Beat AML, we identified candidate drugs showing subgroup-specific effectiveness such as panobinostat for NPM1-mutated AML with monocytic differentiation. Finally, ChIP-seq analysis identified a number of super-enhancers in the vicinity of leukemia-related TFs and oncogenes. Conclusion Through comprehensive profiling of chromatin accessibility and gene mutations/expression in a large cohort of AML patients, we demonstrate the epigenetic heterogeneity and identify novel AML subgroups that had unique epigenetic and clinical features. Our findings highlight the role of epigenetic alterations in AML pathogenesis and impact on molecular classification of AML.

Low CD117 Expression Identifies Hematopoietic Stem Cell Subsets with Enhanced Lymphoid Potential in Mice and Humans

Background: Hematopoietic stem cells (HSCs) with multilineage potential are critical for effective de novo T-cell generation and thymic recovery, restoring the adaptive immune system after Hematopoietic Cell Transplantation (HCT). Recent studies have brought new insights into functional heterogeneity within HSCs, revealing an organized and predictable framework governing the adoption of lineage-restricted fates. However, the gene-regulatory networks underlying lymphoid determination and their conservation with aging in HSCs need to be elucidated. This study aims to comprehensively investigate the molecular profile of young and old HSCs and identify HSC subsets with multilineage transcriptional programs. Results: We conducted single-cell multiomic RNA and ATAC sequencing of young and old HSCs defined by Lineage-CD34-CD48-CD150+Sca-1+Kit+. We identified two CD117 (Kit) HSC subsets, Kit lo and Kit hi, with unique transcriptional profiles. Kit hi HSCs were characterized by quiescence, platelet bias, and low-output gene signatures, while Kit lo HSCs enriched for multi-lineage and high-output gene signatures. We observed a decreased frequency of Kit lo HSCs in aged mice. Although Kit lo HSCs have been previously described to exhibit increased self-renewal (Shin et al. JEM, 2014), their T-cell potential and the molecular underpinnings governing lymphoid differentiation programs are unknown. We performed in vitro studies, S17 stromal assay, and artificial thymic organoids (ATOs), which revealed an enhanced output of lymphoid progenitors and T-cells from Kit lo HSCs, consistent across all age cohorts. We next compared Kit lo vs. Kit hi HSCs in a competitive allogeneic HCT model. Kit lo HSCs demonstrated improved thymic recovery and post-HCT T-cell reconstitution, independent of age. To orthogonally evaluate post-HCT thymic function, we assessed Recent Thymic Emigrants (RTEs) output, which further substantiated superior recovery in mice receiving Kit lo HSCs. We further observed that Kit lo HSC-derived T-cells exhibited better proliferation in response to an acute Listeria monocytogenes infection. Next, we transplanted middle-aged mice to evaluate the rejuvenation potential of different Kit HSC subsets. We observed that Kit lo HSCs partly mitigated age-related changes in the thymic microenvironment, specifically by enhancing the regeneration of thymic epithelial cells (TECs) in middle-aged mice, thereby reversing the decline in T-cell production. We found differential expression and activity of key transcription factors (TFs), including Runx3, Zbtb7a, and Ezh1, associated with lymphoid differentiation through integrative transcriptional and chromatin accessibility analyses of Kit lo HSCs. These TFs notably showed distinct patterns in Kit lo HSCs independent of age, indicating their potential role in driving their enhanced T-cell potential. To establish the existence of a comparable human subset, we interrogated a human BM CITE-Seq dataset (Sommarin et al. Biorxiv 2021) for the mouse Kit lo gene signature. In concordance with our findings in mice, we observed an enrichment of the Kit lo gene program in young HSCs and decreased frequency of Kit lo HSCs in old BM. Next, we used ATOs to validate differential lymphoid potential. We found increased T-cell output originating from Kit lo HSCs, thus underscoring their enhanced lymphoid capacity and potential clinical relevance. Conclusion: Collectively, we demonstrate a distinct HSC subset, Kit lo HSC, with an age-conserved lymphoid gene program that enhances T-cell production and facilitates thymic regeneration. Importantly, we demonstrate the presence of an analogous HSC subset in humans, revealing the relevance of these insights to human health and disease. Thus, Kit lo HSCs have therapeutic potential to counteract age-associated immune senescence and treatment-related immune suppression.

Multiplexed Single-Nucleus RNA and ATAC Sequencing in the Renal Biopsy Specimen

Multiplexed Single-Nucleus RNA and ATAC Sequencing in the Renal Biopsy Specimen

Genome-wide profiling of transcription factor activity in primary liver cancer using single-cell ATAC sequencing

Primary liver cancer (PLC) consists of two main histological subtypes; hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA). The role of transcription factors (TFs) in malignant hepatobiliary lineage commitment between HCC and iCCA remains underexplored. Here, we present genome-wide profiling of transcription regulatory elements of 16 PLC patients using single-cell assay for transposase accessible chromatin sequencing. Single-cell open chromatin profiles reflect the compositional diversity of liver cancer, identifying both malignant and microenvironment component cells. TF motif enrichment levels of 31 TFs strongly discriminate HCC from iCCA tumors. These TFs are members of the nuclear/retinoid receptor, POU, or ETS motif families. POU factors are associated with prognostic features in iCCA. Overall, nuclear receptors, ETS and POU TF motif families delineate transcription regulation between HCC and iCCA tumors, which may be relevant to development and selection of PLC subtype-specific therapeutics.

Exploring intrinsic variability between cultured nasal and bronchial epithelia in cystic fibrosis

The nasal and bronchial epithelium are unified parts of the respiratory tract that are affected in the monogenic disorder cystic fibrosis (CF). Recent studies have uncovered that nasal and bronchial tissues exhibit intrinsic variability, including differences in mucociliary cell composition and expression of unique transcriptional regulatory proteins which relate to germ layer origin. In the present study, we explored whether intrinsic differences between nasal and bronchial epithelial cells persist in cell cultures and affect epithelial cell functioning in CF. Comparison of air–liquid interface (ALI) differentiated epithelial cells from subjects with CF revealed distinct mucociliary differentiation states of nasal and bronchial cultures. Moreover, using RNA sequencing we identified cell type-specific signature transcription factors in differentiated nasal and bronchial epithelial cells, some of which were already poised for expression in basal progenitor cells as evidenced by ATAC sequencing. Analysis of differentiated nasal and bronchial epithelial 3D organoids revealed distinct capacities for fluid secretion, which was linked to differences in ciliated cell differentiation. In conclusion, we show that unique phenotypical and functional features of nasal and bronchial epithelial cells persist in cell culture models, which can be further used to investigate the effects of tissue-specific features on upper and lower respiratory disease development in CF.

Tumor heterogeneity and tumor-microglia interactions in primary and recurrent IDH1-mutant gliomas

The isocitrate dehydrogenase (IDH) gene is recurrently mutated in adult diffuse gliomas. IDH-mutant gliomas are categorized into oligodendrogliomas and astrocytomas, each with unique pathological features. Here, we use single-nucleus RNA and ATAC sequencing to compare the molecular heterogeneity of these glioma subtypes. In addition to astrocyte-like, oligodendrocyte progenitor-like, and cycling tumor subpopulations, a tumor population enriched for ribosomal genes and translation elongation factors is primarily present in oligodendrogliomas. Longitudinal analysis of astrocytomas indicates that the proportion of tumor subpopulations remains stable in recurrent tumors. Analysis of tumor-associated microglia/macrophages (TAMs) reveals significant differences between oligodendrogliomas, with astrocytomas harboring inflammatory TAMs expressing phosphorylated STAT1, as confirmed by immunohistochemistry. Furthermore, inferred receptor-ligand interactions between tumor subpopulations and TAMs may contribute to TAM state diversity. Overall, our study sheds light on distinct tumor populations, TAM heterogeneity, TAM-tumor interactions in IDH-mutant glioma subtypes, and the relative stability of tumor subpopulations in recurrent astrocytomas.

The BAF chromatin remodeling complex licenses planarian stem cells access to ectodermal and mesodermal cell fates

BackgroundThe flatworm planarian, Schmidtea mediterranea, has a large population of adult stem cells (ASCs) that replace any cell type during tissue turnover or regeneration. How planarian ASCs (called neoblasts) manage self-renewal with the ability to produce daughter cells of different cell lineages (multipotency) is not well understood. Chromatin remodeling complexes ultimately control access to DNA regions of chromosomes and together with specific transcription factors determine whether a gene is transcribed in a given cell type. Previous work in planarians determined that RNAi of core components of the BAF chromatin remodeling complex, brg1 and smarcc2, caused increased ASCs and failed regeneration, but how these cellular defects arise at the level of gene regulation in neoblasts is unknown.ResultsHere, we perform ATAC and RNA sequencing on purified neoblasts, deficient for the BAF complex subunits brg-1 and smarcc2. The data demonstrate that the BAF complex promotes chromatin accessibility and facilitates transcription at target loci, as in other systems. Interestingly, we find that the BAF complex enables access to genes known to be required for the generation of mesoderm- and ectoderm-derived lineages, including muscle, parenchymal cathepsin, neural, and epithelial lineages. BAF complex knockdowns result in disrupted differentiation into these cell lineages and functional consequences on planarian regeneration and tissue turnover. Notably, we did not detect a role for the BAF complex in neoblasts making endodermal lineages.ConclusionsOur study provides functional insights into how the BAF complex contributes to cell fate decisions in planarian ASCs in vivo.

Genome-Wide Characterization of TAZ Binding Sites in Mammary Epithelial Cells.

The transcriptional co-activator with PDZ binding motif (TAZ) is a key effector of the Hippo signaling pathway. We and others previously reported that high expression levels of TAZ are positively associated with decreased survival rates and shorter times to relapse in basal-like breast cancer (BLBC) patients. The oncogenic activity of TAZ involves the regulation of diverse signal transduction pathways that direct processes such as cell proliferation, migration, and resistance to apoptosis, albeit through poorly characterized gene expression programs. Here, using a tet-inducible system in mammary epithelial MCF10A cells, we have characterized the TAZ-regulated transcription program using RNA sequencing in a temporal and spatial manner. We further identified global TAZ binding sites at different TAZ activation time points by chromatin immunoprecipitation (ChIP) sequencing analysis. We found that the vast majority of TAZ was rapidly localized in enhancer regions at the early TAZ activation time point and then gradually spread to promoter regions. TAZ bound to enhancer regions following a switch in potential TEAD and FOSL2 transcription factor motifs. Furthermore, the ATAC sequencing analysis indicated that TAZ activation led to chromatin structural alterations. Together, our results have revealed the landscape of genome-wide TAZ binding sites and may lead to improvements in the current understanding of how TAZ regulates the gene expression program that contributes to the development of breast cancer.

Single-cell ATAC sequencing illuminates the cis-regulatory differentiation of taxol biosynthesis between leaf mesophyll and leaf epidermal cells in Taxus mairei

Taxol is a widely used anticancer drug; however, its accumulation pattern can vary largely among leaf tissues. Taxol biosynthesis is driven by a series of transcription factors (TFs), but the extent of chromatin changes that functionally regulate gene expression associated with taxol biosynthesis remain poorly understood. To assess the effects of tissue heterogeneity on taxol biosynthesis, we applied single-cell transposase-accessible chromatin sequencing to Taxus mairei leaves and identified differentially accessible regions. We captured 9488 independent cells with 42,192 fragment reads per cell. In total, 66.64% were confidently mapped onto the T. mairei genome. Projection analysis divided all the cells into nine clusters, which suggests a high degree of cell heterogeneity in T. mairei leaves. Based on the cell type markers, Clusters 2 and 7 were annotated as leaf mesophyll cells, and Clusters 5 and 6 were annotated as leaf epidermal cells. A differential accessibility analysis identified 9600 and 8538 leaf mesophyll and epidermal, respectively, cell-specific open chromatin regions. We predicted the genomic features of the genes involved in different stages of taxol biosynthesis, and there was a close correlation between chromatin open peaks and the expression levels of several taxol biosynthesis-related genes. We further identified a number of cell type-specific TF motifs and reference TFs that appeared to be involved in taxol biosynthesis. The study of T. mairei leaves at a single-cell resolution provides a valuable resource for determining the basic principles of the cell type-specific regulation of taxol biosynthesis.

A scATAC-seq atlas of chromatin accessibility in axolotl brain regions

Axolotl (Ambystoma mexicanum) is an excellent model for investigating regeneration, the interaction between regenerative and developmental processes, comparative genomics, and evolution. The brain, which serves as the material basis of consciousness, learning, memory, and behavior, is the most complex and advanced organ in axolotl. The modulation of transcription factors is a crucial aspect in determining the function of diverse regions within the brain. There is, however, no comprehensive understanding of the gene regulatory network of axolotl brain regions. Here, we utilized single-cell ATAC sequencing to generate the chromatin accessibility landscapes of 81,199 cells from the olfactory bulb, telencephalon, diencephalon and mesencephalon, hypothalamus and pituitary, and the rhombencephalon. Based on these data, we identified key transcription factors specific to distinct cell types and compared cell type functions across brain regions. Our results provide a foundation for comprehensive analysis of gene regulatory programs, which are valuable for future studies of axolotl brain development, regeneration, and evolution, as well as on the mechanisms underlying cell-type diversity in vertebrate brains.

The AD odyssey 2023: Tales of single cell

Deciphering cellular changes in Alzheimer's disease (AD) using large cohorts with defined clinical stages isessential for understanding the diverse trajectories of AD progression. In this issue of Cell, five studies harnessed the power of single-nuclei RNA sequencing (snRNA-seq) and single-nuclei ATAC sequencing (snATAC-seq) at unprecedented scale and revealed exciting insights into cell-type-specific mechanisms underlying the progression of AD pathogenesis.

Coordinated Biophysical Stimulation of MSCs via Electromagnetized Au‐Nanofiber Matrix Regulates Cytoskeletal‐to‐Nuclear Mechanoresponses and Lineage Specification

AbstractBiophysical stimulation regulates stem cell functions, including proliferation and differentiation. Matrix nanotopography and external forces, such as electromagnetic fields (EMF), can enhance this stimulation. Here, it is demonstrated that biophysical multiple cues coordinated from electromagnetized Au‐nanoparticles‐decorated polymer nanofiber under EMF significantly regulate the adhesion, alignment, proliferation, and lineage commitment of hMSCs. Without EMF, matrix cues of electrical conductivity and nanodotted fibrous topography accelerate the anchorage and spreading of hMSCs. Of note, EMF synergizes with the matrix cues to enhance cellular behaviors, resulting in elongated and aligned cells along the field direction. Microtubules are highly polymerized, acetylated, and aligned, playing an active role in these events. Actin filaments also develop in parallel with the microtubules, facilitating actin‐microtubule crosstalks. These phenomena lead to changes in the nuclear mechanics of hMSCs, including elongated nuclear shape and decondensed chromatins with histone acetylation. The EMF+matrix‐stimulated hMSCs express genes related to microtubule organization and euchromatin, as revealed by RNA sequencing, and show chromatin accessibility with enrichment of genes related to mechanotransduction and lineage specification, as analyzed by ATAC sequencing. The EMF+matrix biophysical stimulation further increases the capacity for lineage specification (predominantly towards osteogenic, myogenic, and tenogenic), offering a promising bioengineering platform for stem cell engineering and therapies.

Abstract P1077: Transcription Factor Nfya Regulates Mitochondrial Metabolism And Cardiomyocyte Proliferation In The Fetal Heart

We previously identified a neonatal cardiomyocyte population that enters the cell-cycle following injury and disappears as the heart loses the ability to regenerate. These regenerative neonatal cardiomyocytes display a unique transcriptional program regulated by NFYA, a subunit of the trimeric complex Nuclear Transcription Factor Y (NFY) that functions as a pioneer factor to promote cell proliferation. To study whether NFYA is required for cardiomyocyte proliferation, we generated mice with cardiomyocyte-specific deletion of NFYA using alpha-MHC-Cre. NFYA KO mice died before birth; their hearts were smaller in ventricular volumes and developed a distinct cardiac noncompaction phenotype, with hyper trabeculation and reduced cardiomyocyte proliferation. Integration of single-nuclear RNA and ATAC sequencing with spatial transcriptomics identified different embryonic cardiomyocyte subtypes in wildtype and NFYA KO hearts, including those from atria, ventricles, trabeculae, and the conduction system. We found that NFYA deletion reduced the reservoir of immature regenerative cardiomyocytes in the fetal heart and increased the number of trabecular cardiomyocytes. Mechanistically, NFYA deletion reduced the expression of genes regulating oxidative phosphorylation metabolism. The defect in mitochondrial metabolism was further supported by decreased activity of NADH oxidation at mitochondrial membranes and reduced levels of TCA cycle metabolites in NFYA KO hearts. Using chromatin immunoprecipitation and reporter assays, we show that NFYA directly regulates the transcription of those mitochondrial metabolic genes as well as cell-cycle genes in cardiomyocytes, thus providing a transcriptional coupling between mitochondrial metabolism and cardiomyocyte proliferation. Our study identified a previously unrecognized transcriptional control mechanism of metabolic genes controlled by NFYA and highlights the importance of oxidative metabolism during fetal heart development and regeneration.

Prmt5 deficiency inhibits CD4+ T-cell Klf2/S1pr1 expression and ameliorates EAE disease

BackgroundProtein arginine methyltransferase 5 (Prmt5) is the main type II methyltransferase, catalyzes protein arginine residue symmetric dimethylation, and modulates normal cellular physiology and disease progression. Prmt5 inhibition or deletion in CD4+ T cells has been reported to ameliorate experimental autoimmune encephalomyelitis (EAE), but the detailed molecular mechanisms have not yet been elucidated.MethodsEAE was induced by administration of myelin oligodendrocyte glycoprotein (MOG35–55) in T cells Prmt5 conditional knockout (CD4-cre-Prmt5fl/fl, Prmt5cko) and Prmt5fl/fl (WT) mice. Flow cytometry, single-cell RNA sequencing, ATAC sequencing and chromatin immunoprecipitation assay (ChIP) approaches were used to explore the detail mechanisms.ResultsWe find that Prmt5cko mice are resistant to EAE; infiltrating inflammatory CD4+ T cells in the central nervous system (CNS) are greatly reduced. However, in Prmt5cko mice, T cells in the spleen show much more proliferation and activation properties, the total number of CD4+ T cells in the spleen is not reduced, and the percentage of Rora+ CD4+ T cells is elevated. Also, CD4+ T cells express lower levels of S1pr1 and Klf2 than WT mice, which may influence pathogenic CD4+ T-cell egress from the spleen and migration to the CNS. Moreover, the single-cell ATAC sequence and ChIP assay reveal that the transcription factor Klf2 is enriched at the S1pr1 promoter and that Klf2 motif activity is reduced in Prmt5cko mice.ConclusionsOur study delineates the undiscovered role of Prmt5 in T-cell biology in which Prmt5 may inhibit Klf2-S1pr1 pathway to ameliorate EAE disease. Controlling T-cell Prmt5 expression may be helpful for the treatment of autoimmune diseases.

Unraveling the heterogeneity of intercellular crosstalk in intrahepatic and extrahepatic cholangiocarcinoma using single-cell and spatial transcriptome profiling.

50 Background: Extrahepatic cholangiocarcinoma (ECC) and intrahepatic cholangiocarcinoma (ICC) are two types of malignant tumors that arise from the bile duct epithelium at different locations. Currently, effective treatment options for these tumors are limited, and their microenvironmental differences remain poorly understood. Methods: We collected retrospective clinicopathological data of ECC (N=64) and ICC (N=153) patients from our hospital and integrated public database information to analyze the clinical characteristics, survival, and prognosis of ECC and ICC. We performed single-cell sequencing (scRNA-seq) and spatial transcriptomic techniques on three ECC and three ICC tumor tissues. Verification of scRNA-seq results was done using immunohistochemistry, flow cytometry, and polychromatic immunofluorescence. Additionally, we analyzed RNA and ATAC sequencing data for exhausted CD8+ T cells TEX obtained from ICC and ECC tissue samples. We used InferCNV to infer copy number variations (CNVs) of epithelial cells from ICC and ECC. Results: Our results showed that ECC and ICC differ in their clinical characteristics such as age, gender, bilirubin, CEA, and outcome. ECC had higher copy number amplification variation and exhibited inflammatory signaling and ongoing immunoediting, reflected in higher immune checkpoint expression such as PDCD1, CTLA4, and LAG3, and tumor infiltration with highly differentiated dysfunctional PDCD1+ LAG3+ CXCL13+ exhausted CD8+ T cells (TEX) and precursors of exhausted T (TPEX). RNA velocity analysis revealed different developmental trajectories of T cells in ECC and ICC, predominantly in naive/stem-like and memory T cells. We also identified a trefoil factor 3 (TFF3)-positive epithelial subpopulation in ECC. TFF3+ epithelial cells had closer interaction with CTHRC1+ fibroblasts, which induced more CD8+ T cells to TEX and TPEX. TFF3+ cells exhibited higher CNV scores and higher PD-L1 and PD-L2 expression, and co-localized with TEX and TPEX, providing one explanation for the immunosuppressive environment in ECC and suggesting a potential response to immune checkpoint therapy (ICT). Conclusions: Our study provides a comprehensive characterization of the cellular heterogeneity landscape of ECC and ICC, which can guide the optimization of personalized cancer immunotherapy. We identified potential biomarkers and targets for immunotherapies, and our findings shed light on the underlying mechanisms of tumor immune evasion and provide a foundation for further research in this area.

Multimodal single-cell analysis of non-random heteroplasmy distribution in human retinal mitochondrial disease.

Variants within the high copy number mitochondrial genome (mtDNA) can disrupt organelle function and lead to severe multi-system disease. The wide range of manifestations observed in mitochondrial disease patients results from varying fractions of abnormal mtDNA molecules in different cells and tissues, a phenomenon termed heteroplasmy. However, the landscape of heteroplasmy across cell types within tissues and its influence on phenotype expression in affected patients remains largely unexplored. Here, we identify non-random distribution of a pathogenic mtDNA variant across a complex tissue using single-cell RNA sequencing, mitochondrial single-cell ATAC sequencing, and multimodal single-cell sequencing. We profile the transcriptome, chromatin accessibility state, and heteroplasmy in cells from the eyes of a patient with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and healthy control donors. Utilizing the retina as a model for complex multi-lineage tissues, we found that the proportion of the pathogenic m.3243A>G allele was neither evenly nor randomly distributed across diverse cell types. All neuroectoderm-derived neural cells exhibited a high percentage of the mutant variant. However, a subset of mesoderm-derived lineage, namely the vasculature of the choroid, was near homoplasmic for the wildtype allele. Gene expression and chromatin accessibility profiles of cell types with high and low proportions of m.3243A>G implicate mTOR signaling in the cellular response to heteroplasmy. We further found by multimodal single-cell sequencing of retinal pigment epithelial cells that a high proportion of the pathogenic mtDNA variant was associated with transcriptionally and morphologically abnormal cells. Together, these findings show the non-random nature of mitochondrial variant partitioning in human mitochondrial disease and underscore its implications for mitochondrial disease pathogenesis and treatment.

P53 governs an AT1 differentiation programme in lung cancer suppression.

Lung cancer is the leading cause of cancer deaths worldwide1. Mutations in the tumour suppressor gene TP53 occur in 50% of lung adenocarcinomas (LUADs) and are linked to poor prognosis1-4, but how p53 suppresses LUAD development remains enigmatic. We show here that p53 suppresses LUAD by governing cell state, specifically by promoting alveolar type 1 (AT1) differentiation. Using mice that express oncogenic Kras and null, wild-type or hypermorphic Trp53 alleles in alveolar type 2 (AT2) cells, we observed graded effects of p53 on LUAD initiation and progression. RNA sequencing and ATAC sequencing of LUAD cells uncovered a p53-induced AT1 differentiation programme during tumour suppression in vivo through direct DNA binding, chromatin remodelling and inductionof genes characteristic of AT1 cells. Single-cell transcriptomics analyses revealed that during LUAD evolution, p53 promotes AT1 differentiation throughaction in a transitional cell state analogous to a transient intermediary seen during AT2-to-AT1 cell differentiation in alveolar injury repair. Notably, p53 inactivation results in the inappropriate persistence of these transitional cancer cells accompanied by upregulated growth signalling and divergence from lung lineage identity, characteristics associated with LUAD progression. Analysis of Trp53 wild-type and Trp53-null mice showed that p53 also directs alveolar regeneration after injury by regulating AT2 cell self-renewal and promoting transitional cell differentiation into AT1 cells. Collectively, these findings illuminate mechanisms of p53-mediated LUAD suppression, in which p53 governs alveolar differentiation, and suggest that tumour suppression reflects a fundamental role of p53 in orchestrating tissue repair after injury.

Single Cell Multiomics Identifies Cells and Genetic Networks Underlying Alveolar Capillary Dysplasia.

Rationale: Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a lethal developmental disorder of lung morphogenesis caused by insufficiency of FOXF1 (forkhead box F1) transcription factor function. The cellular and transcriptional mechanisms by which FOXF1 deficiency disrupts human lung formation are unknown. Objectives: To identify cell types, gene networks, and cell-cell interactions underlying the pathogenesis of ACDMPV. Methods: We used single-nucleus RNA and assay for transposase-accessible chromatin sequencing, immunofluorescence confocal microscopy, and RNA in situ hybridization to identify cell types and molecular networks influenced by FOXF1 in ACDMPV lungs. Measurements and Main Results: Pathogenic single-nucleotide variants and copy-number variant deletions involving the FOXF1 gene locus in all subjects with ACDMPV (n = 6) were accompanied by marked changes in lung structure, including deficient alveolar development and a paucity of pulmonary microvasculature. Single-nucleus RNA and assay for transposase-accessible chromatin sequencing identified alterations in cell number and gene expression in endothelial cells (ECs), pericytes, fibroblasts, and epithelial cells in ACDMPV lungs. Distinct cell-autonomous roles for FOXF1 in capillary ECs and pericytes were identified. Pathogenic variants involving the FOXF1 gene locus disrupt gene expression in EC progenitors, inhibiting the differentiation or survival of capillary 2 ECs and cell-cell interactions necessary for both pulmonary vasculogenesis and alveolar type 1 cell differentiation. Loss of the pulmonary microvasculature was associated with increased VEGFA (vascular endothelial growth factor A) signaling and marked expansion of systemic bronchial ECs expressing COL15A1 (collagen type XV α 1 chain). Conclusions: Distinct FOXF1 gene regulatory networks were identified in subsets of pulmonary endothelial and fibroblast progenitors, providing both cellular and molecular targets for the development of therapies for ACDMPV and other diffuse lung diseases of infancy.