A sensitive assay for dNTPs based on long synthetic oligonucleotides, EvaGreen dye and inhibitor-resistant high-fidelity DNA polymerase

Identification and Validation of Biomarkers to Predict Early Diagnosis of Inflammatory Bowel Disease and Its Progression to Colorectal Cancer

Inflammatory bowel disease (IBD) has become a common global health problem as prevalence continues to rise. It is often associated with increased risk of colorectal cancer (CRC) development. Limitations in current IBD biomarker-based diagnosis hinder the accuracy of early detection of CRC progression. Therefore, in this study, we proposed the use of transcription factor (TF)-based biomarkers that can potentially detect the transition of IBD to CRC. Various bioinformatic analysis and online database validations, and RT-qPCR validations were performed to identify possible diagnostic TFs. RUNX1 was identified as a promising TF that regulates 106 IBD/CRC-related genes. The incorporation of RUNX1 in combination with currently known IBD biomarkers, FEV + NFKB1 + RELA, achieved a comparable sensitivity and specificity scores of 99% and 87%, respectively, while RUNX1 in combination with known CRC markers, CEA + TIMP1 + CA724 + CA199, achieved a sensitivity and specificity score of 97% and 99%, respectively. Furthermore, a small pilot RT-qPCR-based analysis confirmed a demarcated shift in expression profiles in CA724, CEA, RUNX1 and TIMP1 in IBD patients compared to CRC patients' tissue samples. Specifically, CA724 is noticeably elevated in IBD, while the levels of CEA, RUNX1 with TIMP1 are probable genes that may be employed in discerning IBD progression to CRC. Therefore, these preliminary results once validated in large patient cohorts could potentially have a significant impact on CRC disease stratification, resulting in a more precise prediction for treatment and treatment outcomes, especially in South African patients.

Increasing intracellular dNTP levels improves prime editing efficiency

In primary cell types, intracellular deoxynucleotide triphosphate (dNTP) levels are tightly regulated in a cell cycle-dependent manner. We report that prime editing efficiency is increased by mutations that improve the enzymatic properties of Moloney murine leukemia virus reverse transcriptase and treatments that increase intracellular dNTP levels. In combination, these modifications produce substantial increases in precise editing rates.

Single Cell Profiling in the Sox10 Dom/+ Hirschsprung Mouse Implicates Hoxa6 in Enteric Neuron Lineage Allocation

Background & Aims

Enteric nervous system (ENS) development requires migration, proliferation, and appropriate neuronal diversification from progenitors to enable normal gastrointestinal (GI) motility. Sox10 deficit causes aganglionosis, modeling Hirschsprung disease, and disrupts ratios of postnatal enteric neurons in proximal ganglionated bowel. How Sox10 deficiency alters ratios of enteric neuron subtypes is unclear. Sox10's prominent expression in enteric neural crest-derived progenitors (ENCP) and lack of this gene in enteric neurons led us to examine Sox10 Dom effects ENS progenitors and early differentiating enteric neurons.

Methods

ENS progenitors, developing neurons, and enteric glia were isolated from Sox10 +/+ and Sox10 Dom/+ littermates for single-cell RNA sequencing (scRNA-seq). scRNA-seq data was processed to identify cell type-specific markers, differentially expressed genes, cell fate trajectories, and gene regulatory network activity between genotypes. Hybridization chain reaction (HCR) validated expression changes detected in scRNA-seq.

Results

scRNA-seq profiles revealed three neuronal lineages emerging from cycling progenitors via two transition pathways accompanied by elevated activity of Hox gene regulatory networks (GRN) as progenitors transition to neuronal fates. Sox10 Dom/+ scRNA-seq profiles exhibited a novel progenitor cluster, decreased abundance of cells in transitional states, and shifts in cell distributions between two neuronal trajectories. Hoxa6 was differentially expressed in the neuronal lineages impacted in Sox10 Dom/+ mutants and HCR identified altered Hoxa6 expression in early developing neurons of Sox10 Dom/+ ENS.

Conclusions

Sox10 Dom/+ mutation shifts enteric neuron types by altering neuronal trajectories during early ENS lineage segregation. Multiple neurogenic transcription factors are reduced in Sox10 Dom/+ scRNA-seq profiles including multiple Hox genes. This is the first report that implicates Hox genes in lineage diversification of enteric neurons.

Interphase chromosome conformation is specified by distinct folding programs inherited via mitotic chromosomes or through the cytoplasm

Identity-specific interphase chromosome conformation must be re-established each time a cell divides. To understand how interphase folding is inherited, we developed an experimental approach that physically segregates mediators of G1 folding that are intrinsic to mitotic chromosomes from cytoplasmic factors. Proteins essential for nuclear transport, RanGAP1 and Nup93, were degraded in pro-metaphase arrested DLD-1 cells to prevent the establishment of nucleo-cytoplasmic transport during mitotic exit and isolate the decondensing mitotic chromatin of G1 daughter cells from the cytoplasm. Using this approach, we discover a transient folding intermediate entirely driven by chromosome-intrinsic factors. In addition to conventional compartmental segregation, this chromosome-intrinsic folding program leads to prominent genome-scale microcompartmentalization of mitotically bookmarked and cell type-specific cis-regulatory elements. This microcompartment conformation is formed during telophase and subsequently modulated by a second folding program driven by factors inherited through the cytoplasm in G1. This nuclear import-dependent folding program includes cohesin and factors involved in transcription and RNA processing. The combined and inter-dependent action of chromosome-intrinsic and cytoplasmic inherited folding programs determines the interphase chromatin conformation as cells exit mitosis.

Genome-wide discovery and integrative genomic characterization of insulin resistance loci using serum triglycerides to HDL-cholesterol ratio as a proxy

Insulin resistance causes multiple epidemic metabolic diseases, including type 2 diabetes, cardiovascular disease, and fatty liver, but is not routinely measured in epidemiological studies. To discover novel insulin resistance genes in the general population, we conducted genome-wide association studies in 382,129 individuals for triglyceride to HDL-cholesterol ratio (TG/HDL), a surrogate marker of insulin resistance calculable from commonly measured serum lipid profiles. We identified 251 independent loci, of which 62 were more strongly associated with TG/HDL compared to TG or HDL alone, suggesting them as insulin resistance loci. Candidate causal genes at these loci were prioritized by fine mapping with directions-of-effect and tissue specificity annotated through analysis of protein coding and expression quantitative trait variation. Directions-of-effect were corroborated in an independent cohort of individuals with directly measured insulin resistance. We highlight two phospholipase encoding genes, PLA2G12A and PLA2G6, which liberate arachidonic acid and improve insulin sensitivity, and VGLL3, a transcriptional co-factor that increases insulin resistance partially through enhanced adiposity. Finally, we implicate the anti-apoptotic gene TNFAIP8 as a sex-dimorphic insulin resistance factor, which acts by increasing visceral adiposity, specifically in females. In summary, our study identifies several candidate modulators of insulin resistance that have the potential to serve as biomarkers and pharmacological targets.

Variants in tubule epithelial regulatory elements mediate most heritable differences in human kidney function

Kidney failure, the decrease of kidney function below a threshold necessary to support life, is a major cause of morbidity and mortality. We performed a genome-wide association study (GWAS) of 406,504 individuals in the UK Biobank, identifying 430 loci affecting kidney function in middle-aged adults. To investigate the cell types affected by these loci, we integrated the GWAS with human kidney candidate cis-regulatory elements (cCREs) identified using single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq). Overall, 56% of kidney function heritability localized to kidney tubule epithelial cCREs and an additional 7% to kidney podocyte cCREs. Thus, most heritable differences in adult kidney function are a result of altered gene expression in these two cell types. Using enhancer assays, allele-specific scATAC-seq and machine learning, we found that many kidney function variants alter tubule epithelial cCRE chromatin accessibility and function. Using CRISPRi, we determined which genes some of these cCREs regulate, implicating NDRG1, CCNB1 and STC1 in human kidney function.

Osteogenic-Like Microenvironment of Renal Interstitium Induced by Osteomodulin Contributes to Randall's Plaque Formation

Calcium oxalate (CaOx) kidney stones are common and recurrent, lacking pharmacological prevention. Randall's plaques (RPs), calcium deposits in renal papillae, serve as niduses for some CaOx stones. This study explores the role of osteogenic-like cells in RP formation resembling ossification. CaP crystals deposit around renal tubules, interstitium, and blood vessels in RP tissues. Human renal interstitial fibroblasts (hRIFs) exhibit the highest osteogenic-like differentiation potential compared to chloride voltage-gated channel Ka positive tubular epithelial cells, aquaporin 2 positive collecting duct cells, and vascular endothelial cells, echoing the upregulated osteogenic markers primarily in hRIFs within RP tissues. Utilizing RNA-seq, osteomodulin (OMD) is found to be upregulated in hRIFs within RP tissues and hRIFs following osteogenic induction. Furthermore, OMD colocalizes with CaP crystals and calcium vesicles within RP tissues. OMD can enhance osteogenic-like differentiation of hRIFs in vitro and in vivo. Additionally, crystal deposits are attenuated in mice with Omd deletion in renal interstitial fibroblasts following CaOx nephrocalcinosis induction. Mechanically, a positive feedback loop of OMD/BMP2/BMPR1A/RUNX2/OMD drives hRIFs to adopt osteogenic-like fates, by which OMD induces osteogenic-like microenvironment of renal interstitium to participate in RP formation. We identify OMD upregulation as a pathological feature of RP, paving the way for preventing CaOx stones.

BCAA mediated microbiota-liver-heart crosstalk regulates diabetic cardiomyopathy via FGF21

BACKGROUND

Diabetic cardiomyopathy (DCM) is one of leading causes of diabetes-associated mortality. The gut microbiota-derived branched-chain amino acids (BCAA) have been reported to play a central role in the onset and progression of DCM, but the potential mechanisms remain elusive.

RESULTS

We found the type 1 diabetes (T1D) mice had higher circulating BCAA levels due to a reduced BCAA degradation ability of the gut microbiota. Excess BCAA decreased hepatic FGF21 production by inhibiting PPARα signaling pathway and thereby resulted in a higher expression level of cardiac LAT1 via transcription factor Zbtb7c. High cardiac LAT1 increased the levels of BCAA in the heart and then caused mitochondrial damage and myocardial apoptosis through mTOR signaling pathway, leading to cardiac fibrosis and dysfunction in T1D mice. Additionally, transplant of faecal microbiota from healthy mice alleviated cardiac dysfunction in T1D mice, but this effect was abolished by FGF21 knockdown.

CONCLUSIONS

Our study sheds light on BCAA-mediated crosstalk among the gut microbiota, liver and heart to promote DCM and FGF21 serves as a key mediator. Video Abstract.

Integrative multi-omics increase resolution of the sea urchin posterior gut gene regulatory network at single-cell level

Drafting gene regulatory networks (GRNs) requires embryological knowledge pertaining to the cell type families, information on the regulatory genes, causal data from gene knockdown experiments and validations of the identified interactions by cis-regulatory analysis. We use multi-omics involving next-generation sequencing to obtain the necessary information for drafting the Strongylocentrotus purpuratus (Sp) posterior gut GRN. Here, we present an update to the GRN using: (1) a single-cell RNA-sequencing-derived cell atlas highlighting the 2 day-post-fertilization (dpf) sea urchin gastrula cell type families, as well as the genes expressed at the single-cell level; (2) a set of putative cis-regulatory modules and transcription factor-binding sites obtained from chromatin accessibility ATAC-seq data; and (3) interactions directionality obtained from differential bulk RNA sequencing following knockdown of the transcription factor Sp-Pdx1, a key regulator of gut patterning in sea urchins. Combining these datasets, we draft the GRN for the hindgut Sp-Pdx1-positive cells in the 2 dpf gastrula embryo. Overall, our data suggest the complex connectivity of the posterior gut GRN and increase the resolution of gene regulatory cascades operating within it.

Decoding the genomic landscape of chromatin-associated biomolecular condensates

Biomolecular condensates play a significant role in chromatin activities, primarily by concentrating and compartmentalizing proteins and/or nucleic acids. However, their genomic landscapes and compositions remain largely unexplored due to a lack of dedicated computational tools for systematic identification in vivo. To address this, we develop CondSigDetector, a computational framework designed to detect condensate-like chromatin-associated protein co-occupancy signatures (CondSigs), to predict genomic loci and component proteins of distinct chromatin-associated biomolecular condensates. Applying this framework to mouse embryonic stem cells (mESC) and human K562 cells enable us to depict the high-resolution genomic landscape of chromatin-associated biomolecular condensates, and uncover both known and potentially unknown biomolecular condensates. Multi-omics analysis and experimental validation further verify the condensation properties of CondSigs. Additionally, our investigation sheds light on the impact of chromatin-associated biomolecular condensates on chromatin activities. Collectively, CondSigDetector provides an approach to decode the genomic landscape of chromatin-associated condensates, facilitating a deeper understanding of their biological functions and underlying mechanisms in cells.

Single-nucleus chromatin accessibility and transcriptomic map of breast tissues of women of diverse genetic ancestry

Single-nucleus analysis allows robust cell-type classification and helps to establish relationships between chromatin accessibility and cell-type-specific gene expression. Here, using samples from 92 women of several genetic ancestries, we developed a comprehensive chromatin accessibility and gene expression atlas of the breast tissue. Integrated analysis revealed ten distinct cell types, including three major epithelial subtypes (luminal hormone sensing, luminal adaptive secretory precursor (LASP) and basal-myoepithelial), two endothelial and adipocyte subtypes, fibroblasts, T cells, and macrophages. In addition to the known cell identity genes FOXA1 (luminal hormone sensing), EHF and ELF5 (LASP), TP63 and KRT14 (basal-myoepithelial), epithelial subtypes displayed several uncharacterized markers and inferred gene regulatory networks. By integrating breast epithelial cell gene expression signatures with spatial transcriptomics, we identified gene expression and signaling differences between lobular and ductal epithelial cells and age-associated changes in signaling networks. LASP cells and fibroblasts showed genetic ancestry-dependent variability. An estrogen receptor-positive subpopulation of LASP cells with alveolar progenitor cell state was enriched in women of Indigenous American ancestry. Fibroblasts from breast tissues of women of African and European ancestry clustered differently, with accompanying gene expression differences. Collectively, these data provide a vital resource for further exploring genetic ancestry-dependent variability in healthy breast biology.

Sequences within and upstream of the mouse Ets1 gene drive high level expression in B cells, but are not sufficient for consistent expression in T cells

The levels of transcription factor Ets1 are high in resting B and T cells, but are downregulated by signaling through antigen receptors and Toll-like receptors (TLRs). Loss of Ets1 in mice leads to excessive immune cell activation and development of an autoimmune syndrome and reduced Ets1 expression has been observed in human PBMCs in the context of autoimmune diseases. In B cells, Ets1 serves to prevent premature activation and differentiation to antibody-secreting cells. Given these important roles for Ets1 in the immune response, stringent control of Ets1 gene expression levels is required for homeostasis. However, the genetic regulatory elements that control expression of the Ets1 gene remain relatively unknown. Here we identify a topologically-associating domain (TAD) in the chromatin of B cells that includes the mouse Ets1 gene locus and describe an interaction hub that extends over 100 kb upstream and into the gene body. Additionally, we compile epigenetic datasets to find several putative regulatory elements within the interaction hub by identifying regions of high DNA accessibility and enrichment of active enhancer histone marks. Using reporter constructs, we determine that DNA sequences within this interaction hub are sufficient to direct reporter gene expression in lymphoid tissues of transgenic mice. Further analysis indicates that the reporter construct drives faithful expression of the reporter gene in mouse B cells, but variegated expression in T cells, suggesting the existence of T cell regulatory elements outside this region. To investigate how the downregulation of Ets1 transcription is associated with alterations in the epigenetic landscape of stimulated B cells, we performed ATAC-seq in resting and BCR-stimulated primary B cells and identified four regions within and upstream of the Ets1 locus that undergo changes in chromatin accessibility that correlate to Ets1 gene expression. Interestingly, functional analysis of several putative Ets1 regulatory elements using luciferase constructs suggested a high level of functional redundancy. Taken together our studies reveal a complex network of regulatory elements and transcription factors that coordinate the B cell-specific expression of Ets1.

Molecular and cellular dynamics of the developing human neocortex at single-cell resolution

The development of the human neocortex is a highly dynamic process and involves complex cellular trajectories controlled by cell-type-specific gene regulation1. Here, we collected paired single-nucleus chromatin accessibility and transcriptome data from 38 human neocortical samples encompassing both the prefrontal cortex and primary visual cortex. These samples span five main developmental stages, ranging from the first trimester to adolescence. In parallel, we performed spatial transcriptomic analysis on a subset of the samples to illustrate spatial organization and intercellular communication. This atlas enables us to catalog cell type-, age-, and area-specific gene regulatory networks underlying neural differentiation. Moreover, combining single-cell profiling, progenitor purification, and lineage-tracing experiments, we have untangled the complex lineage relationships among progenitor subtypes during the transition from neurogenesis to gliogenesis in the human neocortex. We identified a tripotential intermediate progenitor subtype, termed Tri-IPC, responsible for the local production of GABAergic neurons, oligodendrocyte precursor cells, and astrocytes. Remarkably, most glioblastoma cells resemble Tri-IPCs at the transcriptomic level, suggesting that cancer cells hijack developmental processes to enhance growth and heterogeneity. Furthermore, by integrating our atlas data with large-scale GWAS data, we created a disease-risk map highlighting enriched ASD risk in second-trimester intratelencephalic projection neurons. Our study sheds light on the gene regulatory landscape and cellular dynamics of the developing human neocortex.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

Single-nucleus multi-omics of Parkinson's disease reveals a glutamatergic neuronal subtype susceptible to gene dysregulation via alteration of transcriptional networks

The genetic architecture of Parkinson's disease (PD) is complex and multiple brain cell subtypes are involved in the neuropathological progression of the disease. Here we aimed to advance our understanding of PD genetic complexity at a cell subtype precision level. Using parallel single-nucleus (sn)RNA-seq and snATAC-seq analyses we simultaneously profiled the transcriptomic and chromatin accessibility landscapes in temporal cortex tissues from 12 PD compared to 12 control subjects at a granular single cell resolution. An integrative bioinformatic pipeline was developed and applied for the analyses of these snMulti-omics datasets. The results identified a subpopulation of cortical glutamatergic excitatory neurons with remarkably altered gene expression in PD, including differentially-expressed genes within PD risk loci identified in genome-wide association studies (GWAS). This was the only neuronal subtype showing significant and robust overexpression of SNCA. Further characterization of this neuronal-subpopulation showed upregulation of specific pathways related to axon guidance, neurite outgrowth and post-synaptic structure, and downregulated pathways involved in presynaptic organization and calcium response. Additionally, we characterized the roles of three molecular mechanisms in governing PD-associated cell subtype-specific dysregulation of gene expression: (1) changes in cis-regulatory element accessibility to transcriptional machinery; (2) changes in the abundance of master transcriptional regulators, including YY1, SP3, and KLF16; (3) candidate regulatory variants in high linkage disequilibrium with PD-GWAS genomic variants impacting transcription factor binding affinities. To our knowledge, this study is the first and the most comprehensive interrogation of the multi-omics landscape of PD at a cell-subtype resolution. Our findings provide new insights into a precise glutamatergic neuronal cell subtype, causal genes, and non-coding regulatory variants underlying the neuropathological progression of PD, paving the way for the development of cell- and gene-targeted therapeutics to halt disease progression as well as genetic biomarkers for early preclinical diagnosis.

cBAF generates subnucleosomes that expand OCT4 binding and function beyond DNA motifs at enhancers

The canonical BRG/BRM-associated factor (cBAF) complex is essential for chromatin opening at enhancers in mammalian cells. However, the nature of the open chromatin remains unclear. Here, we show that, in addition to producing histone-free DNA, cBAF generates stable hemisome-like subnucleosomal particles containing the four core histones associated with 50-80 bp of DNA. Our genome-wide analysis indicates that cBAF makes these particles by targeting and splitting fragile nucleosomes. In mouse embryonic stem cells, these subnucleosomes become an in vivo binding substrate for the master transcription factor OCT4 independently of the presence of OCT4 DNA motifs. At enhancers, the OCT4-subnucleosome interaction increases OCT4 occupancy and amplifies the genomic interval bound by OCT4 by up to one order of magnitude compared to the region occupied on histone-free DNA. We propose that cBAF-dependent subnucleosomes orchestrate a molecular mechanism that projects OCT4 function in chromatin opening beyond its DNA motifs.

Single-cell and spatial atlases of spinal cord injury in the Tabulae Paralytica

Here, we introduce the Tabulae Paralytica-a compilation of four atlases of spinal cord injury (SCI) comprising a single-nucleus transcriptome atlas of half a million cells, a multiome atlas pairing transcriptomic and epigenomic measurements within the same nuclei, and two spatial transcriptomic atlases of the injured spinal cord spanning four spatial and temporal dimensions. We integrated these atlases into a common framework to dissect the molecular logic that governs the responses to injury within the spinal cord1. The Tabulae Paralytica uncovered new biological principles that dictate the consequences of SCI, including conserved and divergent neuronal responses to injury; the priming of specific neuronal subpopulations to upregulate circuit-reorganizing programs after injury; an inverse relationship between neuronal stress responses and the activation of circuit reorganization programs; the necessity of re-establishing a tripartite neuroprotective barrier between immune-privileged and extra-neural environments after SCI and a failure to form this barrier in old mice. We leveraged the Tabulae Paralytica to develop a rejuvenative gene therapy that re-established this tripartite barrier, and restored the natural recovery of walking after paralysis in old mice. The Tabulae Paralytica provides a window into the pathobiology of SCI, while establishing a framework for integrating multimodal, genome-scale measurements in four dimensions to study biology and medicine.

Estrogen-related receptor alpha promotes thyroid tumor cell survival via a tumor subtype-specific regulation of target gene networks

Mortalin (encoded by HSPA9) is a mitochondrial chaperone often overexpressed in cancer through as-yet-unknown mechanisms. By searching different RNA-sequencing datasets, we found that ESRRA is a transcription factor highly correlated with HSPA9 in thyroid cancer, especially in follicular, but not C cell-originated, tumors. Consistent with this correlation, ESRRA depletion decreased mortalin expression only in follicular thyroid tumor cells. Further, ESRRA expression and activity were relatively high in thyroid tumors with oncocytic characteristics, wherein ESRRA and mortalin exhibited relatively high functional overlap. Mechanistically, ESRRA directly regulated HSPA9 transcription through a novel ESRRA-responsive element located upstream of the HSPA9 promoter. Physiologically, ESRRA depletion suppressed thyroid tumor cell survival via caspase-dependent apoptosis, which ectopic mortalin expression substantially abrogated. ESRRA depletion also effectively suppressed tumor growth and mortalin expression in the xenografts of oncocytic or ESRRA-overexpressing human thyroid tumor cells in mice. Notably, our Bioinformatics analyses of patient data revealed two ESRRA target gene clusters that contrast oncocytic-like and anaplastic features of follicular thyroid tumors. These findings suggest that ESRRA is a tumor-specific regulator of mortalin expression, the ESRRA-mortalin axis has higher significance in tumors with oncocytic characteristics, and ESRRA target gene networks can refine molecular classification of thyroid cancer.

Variants in tubule epithelial regulatory elements mediate most heritable differences in human kidney function

Kidney disease is highly heritable; however, the causal genetic variants, the cell types in which these variants function, and the molecular mechanisms underlying kidney disease remain largely unknown. To identify genetic loci affecting kidney function, we performed a GWAS using multiple kidney function biomarkers and identified 462 loci. To begin to investigate how these loci affect kidney function, we generated single-cell chromatin accessibility (scATAC-seq) maps of the human kidney and identified candidate cis-regulatory elements (cCREs) for kidney podocytes, tubule epithelial cells, and kidney endothelial, stromal, and immune cells. Kidney tubule epithelial cCREs explained 58% of kidney function SNP-heritability and kidney podocyte cCREs explained an additional 6.5% of SNP-heritability. In contrast, little kidney function heritability was explained by kidney endothelial, stromal, or immune cell-specific cCREs. Through functionally informed fine-mapping, we identified putative causal kidney function variants and their corresponding cCREs. Using kidney scATAC-seq data, we created a deep learning model (which we named ChromKid) to predict kidney cell type-specific chromatin accessibility from sequence. ChromKid and allele specific kidney scATAC-seq revealed that many fine-mapped kidney function variants locally change chromatin accessibility in tubule epithelial cells. Enhancer assays confirmed that fine-mapped kidney function variants alter tubule epithelial regulatory element function. To map the genes which these regulatory elements control, we used CRISPR interference (CRISPRi) to target these regulatory elements in tubule epithelial cells and assessed changes in gene expression. CRISPRi of enhancers harboring kidney function variants regulated NDRG1 and RBPMS expression. Thus, inherited differences in tubule epithelial NDRG1 and RBPMS expression may predispose to kidney disease in humans. We conclude that genetic variants affecting tubule epithelial regulatory element function account for most SNP-heritability of human kidney function. This work provides an experimental approach to identify the variants, regulatory elements, and genes involved in polygenic disease.

Copper(II) Complexes with Isomeric Morpholine-Substituted 2-Formylpyridine Thiosemicarbazone Hybrids as Potential Anticancer Drugs Inhibiting Both Ribonucleotide Reductase and Tubulin Polymerization: The Morpholine Position Matters

The development of copper(II) thiosemicarbazone complexes as potential anticancer agents, possessing dual functionality as inhibitors of R2 ribonucleotide reductase (RNR) and tubulin polymerization by binding at the colchicine site, presents a promising avenue for enhancing therapeutic effectiveness. Herein, we describe the syntheses and physicochemical characterization of four isomeric proligands H2L3-H2L6, with the methylmorpholine substituent at pertinent positions of the pyridine ring, along with their corresponding Cu(II) complexes 3-6. Evidently, the position of the morpholine moiety and the copper(II) complex formation have marked effects on the in vitro antiproliferative activity in human uterine sarcoma MES-SA cells and the multidrug-resistant derivative MES-SA/Dx5 cells. Activity correlated strongly with quenching of the tyrosyl radical (Y•) of mouse R2 RNR protein, inhibition of RNR activity in the cancer cells, and inhibition of tubulin polymerization. Insights into the mechanism of antiproliferative activity, supported by experimental results and molecular modeling calculations, are presented.

Leveraging single cell multiomic analyses to identify factors that drive human chondrocyte cell fate

Cartilage plays a crucial role in skeletal development and function, and abnormal development contributes to genetic and age-related skeletal disease. To better understand how human cartilage develops in vivo , we jointly profiled the transcriptome and open chromatin regions in individual nuclei recovered from distal femurs at 2 fetal timepoints. We used these multiomic data to identify transcription factors expressed in distinct chondrocyte subtypes, link accessible regulatory elements with gene expression, and predict transcription factor-based regulatory networks that are important for growth plate or epiphyseal chondrocyte differentiation. We developed a human pluripotent stem cell platform for interrogating the function of predicted transcription factors during chondrocyte differentiation and used it to test NFATC2 . We expect new regulatory networks we uncovered using multiomic data to be important for promoting cartilage health and treating disease, and our platform to be a useful tool for studying cartilage development in vitro .

Statement of Significance

The identity and integrity of the articular cartilage lining our joints are crucial to pain-free activities of daily living. Here we identified a gene regulatory landscape of human chondrogenesis at single cell resolution, which is expected to open new avenues of research aimed at mitigating cartilage diseases that affect hundreds of millions of individuals world-wide.

Single cell regulatory architecture of human pancreatic islets suggests sex differences in β cell function and the pathogenesis of type 2 diabetes

Biological sex affects the pathogenesis of type 2 and type 1 diabetes (T2D, T1D) including the development of β cell failure observed more often in males. The mechanisms that drive sex differences in β cell failure is unknown. Studying sex differences in islet regulation and function represent a unique avenue to understand the sex-specific heterogeneity in β cell failure in diabetes. Here, we examined sex and race differences in human pancreatic islets from up to 52 donors with and without T2D (including 37 donors from the Human Pancreas Analysis Program [HPAP] dataset) using an orthogonal series of experiments including single cell RNA-seq (scRNA-seq), single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq), dynamic hormone secretion, and bioenergetics. In cultured islets from nondiabetic (ND) donors, in the absence of the in vivo hormonal environment, sex differences in islet cell type gene accessibility and expression predominantly involved sex chromosomes. Of particular interest were sex differences in the X-linked KDM6A and Y-linked KDM5D chromatin remodelers in female and male islet cells respectively. Islets from T2D donors exhibited similar sex differences in differentially expressed genes (DEGs) from sex chromosomes. However, in contrast to islets from ND donors, islets from T2D donors exhibited major sex differences in DEGs from autosomes. Comparing β cells from T2D and ND donors revealed that females had more DEGs from autosomes compared to male β cells. Gene set enrichment analysis of female β cell DEGs showed a suppression of oxidative phosphorylation and electron transport chain pathways, while male β cell had suppressed insulin secretion pathways. Thus, although sex-specific differences in gene accessibility and expression of cultured ND human islets predominantly affect sex chromosome genes, major differences in autosomal gene expression between sexes appear during the transition to T2D and which highlight mitochondrial failure in female β cells.

Spatially resolved multiomics on the neuronal effects induced by spaceflight in mice

Impairment of the central nervous system (CNS) poses a significant health risk for astronauts during long-duration space missions. In this study, we employed an innovative approach by integrating single-cell multiomics (transcriptomics and chromatin accessibility) with spatial transcriptomics to elucidate the impact of spaceflight on the mouse brain in female mice. Our comparative analysis between ground control and spaceflight-exposed animals revealed significant alterations in essential brain processes including neurogenesis, synaptogenesis and synaptic transmission, particularly affecting the cortex, hippocampus, striatum and neuroendocrine structures. Additionally, we observed astrocyte activation and signs of immune dysfunction. At the pathway level, some spaceflight-induced changes in the brain exhibit similarities with neurodegenerative disorders, marked by oxidative stress and protein misfolding. Our integrated spatial multiomics approach serves as a stepping stone towards understanding spaceflight-induced CNS impairments at the level of individual brain regions and cell types, and provides a basis for comparison in future spaceflight studies. For broader scientific impact, all datasets from this study are available through an interactive data portal, as well as the National Aeronautics and Space Administration (NASA) Open Science Data Repository (OSDR).

PPIA dictates NRF2 stability to promote lung cancer progression

Nuclear factor erythroid 2-related factor 2 (NRF2) hyperactivation has been established as an oncogenic driver in a variety of human cancers, including non-small cell lung cancer (NSCLC). However, despite massive efforts, no specific therapy is currently available to target NRF2 hyperactivation. Here, we identify peptidylprolyl isomerase A (PPIA) is required for NRF2 protein stability. Ablation of PPIA promotes NRF2 protein degradation and blocks NRF2-driven growth in NSCLC cells. Mechanistically, PPIA physically binds to NRF2 and blocks the access of ubiquitin/Kelch Like ECH Associated Protein 1 (KEAP1) to NRF2, thus preventing ubiquitin-mediated degradation. Our X-ray co-crystal structure reveals that PPIA directly interacts with a NRF2 interdomain linker via a trans-proline 174-harboring hydrophobic sequence. We further demonstrate that an FDA-approved drug, cyclosporin A (CsA), impairs the interaction of NRF2 with PPIA, inducing NRF2 ubiquitination and degradation. Interestingly, CsA interrupts glutamine metabolism mediated by the NRF2/KLF5/SLC1A5 pathway, consequently suppressing the growth of NRF2-hyperactivated NSCLC cells. CsA and a glutaminase inhibitor combination therapy significantly retard tumor progression in NSCLC patient-derived xenograft (PDX) models with NRF2 hyperactivation. Our study demonstrates that targeting NRF2 protein stability is an actionable therapeutic approach to treat NRF2-hyperactivated NSCLC.

SIMBA: single-cell embedding along with features

Most current single-cell analysis pipelines are limited to cell embeddings and rely heavily on clustering, while lacking the ability to explicitly model interactions between different feature types. Furthermore, these methods are tailored to specific tasks, as distinct single-cell problems are formulated differently. To address these shortcomings, here we present SIMBA, a graph embedding method that jointly embeds single cells and their defining features, such as genes, chromatin-accessible regions and DNA sequences, into a common latent space. By leveraging the co-embedding of cells and features, SIMBA allows for the study of cellular heterogeneity, clustering-free marker discovery, gene regulation inference, batch effect removal and omics data integration. We show that SIMBA provides a single framework that allows diverse single-cell problems to be formulated in a unified way and thus simplifies the development of new analyses and extension to new single-cell modalities. SIMBA is implemented as a comprehensive Python library ( https://simba-bio.readthedocs.io ).

Scalable and unbiased sequence-informed embedding of single-cell ATAC-seq data with CellSpace

Standard scATAC sequencing (scATAC-seq) analysis pipelines represent cells as sparse numeric vectors relative to an atlas of peaks or genomic tiles and consequently ignore genomic sequence information at accessible loci. Here we present CellSpace, an efficient and scalable sequence-informed embedding algorithm for scATAC-seq that learns a mapping of DNA k-mers and cells to the same space, to address this limitation. We show that CellSpace captures meaningful latent structure in scATAC-seq datasets, including cell subpopulations and developmental hierarchies, and can score transcription factor activities in single cells based on proximity to binding motifs embedded in the same space. Importantly, CellSpace implicitly mitigates batch effects arising from multiple samples, donors or assays, even when individual datasets are processed relative to different peak atlases. Thus, CellSpace provides a powerful tool for integrating and interpreting large-scale scATAC-seq compendia.

Conserved regulatory switches for the transition from natal down to juvenile feather in birds

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We report that extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial-mesenchymal interactions for branching morphogenesis. α-SMA (ACTA2) compartmentalizes dermal papilla stem cells for feather renewal cycling. LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We show that this primary feather transition is largely conserved in chicken (precocial) and zebra finch (altricial) and discuss the possibility that this evolutionary adaptation process started in feathered dinosaurs.

Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits

Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.

Hi-Tag: a simple and efficient method for identifying protein-mediated long-range chromatin interactions with low cell numbers

Protein-mediated chromatin interactions can be revealed by coupling proximity-based ligation with chromatin immunoprecipitation. However, these techniques require complex experimental procedures and millions of cells per experiment, which limits their widespread application in life science research. Here, we develop a novel method, Hi-Tag, that identifies high-resolution, long-range chromatin interactions through transposase tagmentation and chromatin proximity ligation (with a phosphorothioate-modified linker). Hi-Tag can be implemented using as few as 100,000 cells, involving simple experimental procedures that can be completed within 1.5 days. Meanwhile, Hi-Tag is capable of using its own data to identify the binding sites of specific proteins, based on which, it can acquire accurate interaction information. Our results suggest that Hi-Tag has great potential for advancing chromatin interaction studies, particularly in the context of limited cell availability.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Temporal coordination of the transcription factor response to H2O2 stress

Oxidative stress from excess H2O2 activates transcription factors that restore redox balance and repair oxidative damage. Although many transcription factors are activated by H2O2, it is unclear whether they are activated at the same H2O2 concentration, or time. Dose-dependent activation is likely as oxidative stress is not a singular state and exhibits dose-dependent outcomes including cell-cycle arrest and cell death. Here, we show that transcription factor activation is both dose-dependent and coordinated over time. Low levels of H2O2 activate p53, NRF2 and JUN. Yet under high H2O2, these transcription factors are repressed, and FOXO1, NF-κB, and NFAT1 are activated. Time-lapse imaging revealed that the order in which these two groups of transcription factors are activated depends on whether H2O2 is administered acutely by bolus addition, or continuously through the glucose oxidase enzyme. Finally, we provide evidence that 2-Cys peroxiredoxins control which group of transcription factors are activated.

Exploring the role of Prx II in mitigating endoplasmic reticulum stress and mitochondrial dysfunction in neurodegeneration

BACKGROUND

Neurodegenerative diseases are increasingly recognized for their association with oxidative stress, which leads to progressive dysfunction and loss of neurons, manifesting in cognitive and motor impairments. This study aimed to elucidate the neuroprotective role of peroxiredoxin II (Prx II) in counteracting oxidative stress-induced mitochondrial damage, a key pathological feature of neurodegeneration.

METHODS

We investigated the impact of Prx II deficiency on endoplasmic reticulum stress and mitochondrial dysfunction using HT22 cell models with knocked down and overexpressed Prx II. We observed alcohol-treated HT22 cells using transmission electron microscopy and monitored changes in the length of mitochondria-associated endoplasmic reticulum membranes and their contact with endoplasmic reticulum mitochondria contact sites (EMCSs). Additionally, RNA sequencing and bioinformatic analysis were conducted to identify the role of Prx II in regulating mitochondrial transport and the formation of EMCSs.

RESULTS

Our results indicated that Prx II preserves mitochondrial integrity by facilitating the formation of EMCSs, which are essential for maintaining mitochondrial Ca2+ homeostasis and preventing mitochondria-dependent apoptosis. Further, we identified a novel regulatory axis involving Prx II, the transcription factor ATF3, and miR-181b-5p, which collectively modulate the expression of Armcx3, a protein implicated in mitochondrial transport. Our findings underscore the significance of Prx II in protecting neuronal cells from alcohol-induced oxidative damage and suggest that modulating the Prx II-ATF3-miR-181b-5p pathway may offer a promising therapeutic strategy against neurodegenerative diseases.

CONCLUSIONS

This study not only expands our understanding of the cytoprotective mechanisms of Prx II but also offers necessary data for developing targeted interventions to bolster mitochondrial resilience in neurodegenerative conditions.

Harnessing alkaline-pH regulatable promoters for efficient methanol-free expression of enzymes of industrial interest in Komagataella Phaffii

BACKGROUND

The yeast Komagataella phaffii has become a very popular host for heterologous protein expression, very often based on the use of the AOX1 promoter, which becomes activated when cells are grown with methanol as a carbon source. However, the use of methanol in industrial settings is not devoid of problems, and therefore, the search for alternative expression methods has become a priority in the last few years.

RESULTS

We recently reported that moderate alkalinization of the medium triggers a fast and wide transcriptional response in K. phaffii. Here, we present the utilization of three alkaline pH-responsive promoters (pTSA1, pHSP12 and pPHO89) to drive the expression of a secreted phytase enzyme by simply shifting the pH of the medium to 8.0. These promoters offer a wide range of strengths, and the production of phytase could be modulated by adjusting the pH to specific values. The TSA1 and PHO89 promoters offered exquisite regulation, with virtually no enzyme production at acidic pH, while limitation of Pi in the medium further potentiated alkaline pH-driven phytase expression from the PHO89 promoter. An evolved strain based on this promoter was able to produce twice as much phytase as the reference pAOX1-based strain. Functional mapping of the TSA1 and HSP12 promoters suggests that both contain at least two alkaline pH-sensitive regulatory regions.

CONCLUSIONS

Our work shows that the use of alkaline pH-regulatable promoters could be a useful alternative to methanol-based expression systems, offering advantages in terms of simplicity, safety and economy.

Non-canonical transcriptional regulation of the poor prognostic factor UGT2B17 in chronic lymphocytic leukemic and normal B cells

BACKGROUND

High expression of the glycosyltransferase UGT2B17 represents an independent adverse prognostic marker in chronic lymphocytic leukemia (CLL). It also constitutes a predictive marker for therapeutic response and a drug resistance mechanism. The key determinants driving expression of the UGT2B17 gene in normal and leukemic B-cells remain undefined. The UGT2B17 transcriptome is complex and is comprised of at least 10 alternative transcripts, identified by previous RNA-sequencing of liver and intestine. We hypothesized that the transcriptional program regulating UGT2B17 in B-lymphocytes is distinct from the canonical expression previously characterized in the liver.

RESULTS

RNA-sequencing and genomics data revealed a specific genomic landscape at the UGT2B17 locus in normal and leukemic B-cells. RNA-sequencing and quantitative PCR data indicated that the UGT2B17 enzyme is solely encoded by alternative transcripts expressed in CLL patient cells and not by the canonical transcript widely expressed in the liver and intestine. Chromatin accessible regions (ATAC-Seq) in CLL cells mapped with alternative promoters and non-coding exons, which may be derived from endogenous retrotransposon elements. By luciferase reporter assays, we identified key cis-regulatory STAT3, RELA and interferon regulatory factor (IRF) binding sequences driving the expression of UGT2B17 in lymphoblastoid and leukemic B-cells. Electrophoretic mobility shift assays and pharmacological inhibition demonstrated key roles for the CLL prosurvival transcription factors STAT3 and NF-κB in the leukemic expression of UGT2B17.

CONCLUSIONS

UGT2B17 expression in B-CLL is driven by key regulators of CLL progression. Our data suggest that a NF-κB/STAT3/IRF/UGT2B17 axis may represent a novel B-cell pathway promoting disease progression and drug resistance.

Cancer-associated fibroblast-derived acetate promotes pancreatic cancer development by altering polyamine metabolism via the ACSS2-SP1-SAT1 axis

The ability of tumour cells to thrive in harsh microenvironments depends on the utilization of nutrients available in the milieu. Here we show that pancreatic cancer-associated fibroblasts (CAFs) regulate tumour cell metabolism through the secretion of acetate, which can be blocked by silencing ATP citrate lyase (ACLY) in CAFs. We further show that acetyl-CoA synthetase short-chain family member 2 (ACSS2) channels the exogenous acetate to regulate the dynamic cancer epigenome and transcriptome, thereby facilitating cancer cell survival in an acidic microenvironment. Comparative H3K27ac ChIP-seq and RNA-seq analyses revealed alterations in polyamine homeostasis through regulation of SAT1 gene expression and enrichment of the SP1-responsive signature. We identified acetate/ACSS2-mediated acetylation of SP1 at the lysine 19 residue that increased SP1 protein stability and transcriptional activity. Genetic or pharmacologic inhibition of the ACSS2-SP1-SAT1 axis diminished the tumour burden in mouse models. These results reveal that the metabolic flexibility imparted by the stroma-derived acetate enabled cancer cell survival under acidosis via the ACSS2-SP1-SAT1 axis.

Distribution-Agnostic Deep Learning Enables Accurate Single-Cell Data Recovery and Transcriptional Regulation Interpretation

Single-cell RNA sequencing (scRNA-seq) is a robust method for studying gene expression at the single-cell level, but accurately quantifying genetic material is often hindered by limited mRNA capture, resulting in many missing expression values. Existing imputation methods rely on strict data assumptions, limiting their broader application, and lack reliable supervision, leading to biased signal recovery. To address these challenges, authors developed Bis, a distribution-agnostic deep learning model for accurately recovering missing sing-cell gene expression from multiple platforms. Bis is an optimal transport-based autoencoder model that can capture the intricate distribution of scRNA-seq data while addressing the characteristic sparsity by regularizing the cellular embedding space. Additionally, they propose a module using bulk RNA-seq data to guide reconstruction and ensure expression consistency. Experimental results show Bis outperforms other models across simulated and real datasets, showcasing superiority in various downstream analyses including batch effect removal, clustering, differential expression analysis, and trajectory inference. Moreover, Bis successfully restores gene expression levels in rare cell subsets in a tumor-matched peripheral blood dataset, revealing developmental characteristics of cytokine-induced natural killer cells within a head and neck squamous cell carcinoma microenvironment.

MYC activity at enhancers drives prognostic transcriptional programs through an epigenetic switch

The transcription factor MYC is overexpressed in most cancers, where it drives multiple hallmarks of cancer progression. MYC is known to promote oncogenic transcription by binding to active promoters. In addition, MYC has also been shown to invade distal enhancers when expressed at oncogenic levels, but this enhancer binding has been proposed to have low gene-regulatory potential. Here, we demonstrate that MYC directly regulates enhancer activity to promote cancer type-specific gene programs predictive of poor patient prognosis. MYC induces transcription of enhancer RNA through recruitment of RNA polymerase II (RNAPII), rather than regulating RNAPII pause-release, as is the case at promoters. This process is mediated by MYC-induced H3K9 demethylation and acetylation by GCN5, leading to enhancer-specific BRD4 recruitment through its bromodomains, which facilitates RNAPII recruitment. We propose that MYC drives prognostic cancer type-specific gene programs through induction of an enhancer-specific epigenetic switch, which can be targeted by BET and GCN5 inhibitors.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Single-cell epigenomic dysregulation of Systemic Sclerosis fibroblasts via CREB1/EGR1 axis in self-assembled human skin equivalents

Systemic sclerosis (SSc) is an autoimmune disease characterized by skin fibrosis, internal organ involvement and vascular dropout. We previously developed and phenotypically characterized an in vitro 3D skin-like tissue model of SSc, and now analyze the transcriptomic (scRNA-seq) and epigenetic (scATAC-seq) characteristics of this model at single-cell resolution. SSc 3D skin-like tissues were fabricated using autologous fibroblasts, macrophages, and plasma from SSc patients or healthy control (HC) donors. SSc tissues displayed increased dermal thickness and contractility, as well as increased α-SMA staining. Single-cell transcriptomic and epigenomic analyses identified keratinocytes, macrophages, and five populations of fibroblasts (labeled FB1 - 5). Notably, FB1 APOE-expressing fibroblasts were 12-fold enriched in SSc tissues and were characterized by high EGR1 motif accessibility. Pseudotime analysis suggests that FB1 fibroblasts differentiate from a TGF-β1-responsive fibroblast population and ligand-receptor analysis indicates that the FB1 fibroblasts are active in macrophage crosstalk via soluble ligands including FGF2 and APP. These findings provide characterization of the 3D skin-like model at single cell resolution and establish that it recapitulates subsets of fibroblasts and macrophage phenotypes observed in skin biopsies.

txci-ATAC-seq: a massive-scale single-cell technique to profile chromatin accessibility

We develop a large-scale single-cell ATAC-seq method by combining Tn5-based pre-indexing with 10× Genomics barcoding, enabling the indexing of up to 200,000 nuclei across multiple samples in a single reaction. We profile 449,953 nuclei across diverse tissues, including the human cortex, mouse brain, human lung, mouse lung, mouse liver, and lung tissue from a club cell secretory protein knockout (CC16-/-) model. Our study of CC16-/- nuclei uncovers previously underappreciated technical artifacts derived from remnant 129 mouse strain genetic material, which cause profound cell-type-specific changes in regulatory elements near many genes, thereby confounding the interpretation of this commonly referenced mouse model.

Aberrant non-canonical NF-κB signalling reprograms the epigenome landscape to drive oncogenic transcriptomes in multiple myeloma

In multiple myeloma, abnormal plasma cells establish oncogenic niches within the bone marrow by engaging the NF-κB pathway to nurture their survival while they accumulate pro-proliferative mutations. Under these conditions, many cases eventually develop genetic abnormalities endowing them with constitutive NF-κB activation. Here, we find that sustained NF-κB/p52 levels resulting from such mutations favours the recruitment of enhancers beyond the normal B-cell repertoire. Furthermore, through targeted disruption of p52, we characterise how such enhancers are complicit in the formation of super-enhancers and the establishment of cis-regulatory interactions with myeloma dependencies during constitutive activation of p52. Finally, we functionally validate the pathological impact of these cis-regulatory modules on cell and tumour phenotypes using in vitro and in vivo models, confirming RGS1 as a p52-dependent myeloma driver. We conclude that the divergent epigenomic reprogramming enforced by aberrant non-canonical NF-κB signalling potentiates transcriptional programs beneficial for multiple myeloma progression.

Id2 epigenetically controls CD8+ T-cell exhaustion by disrupting the assembly of the Tcf3-LSD1 complex

CD8+ T-cell exhaustion is a state of dysfunction that promotes tumor progression and is marked by the generation of Slamf6+ progenitor exhausted (Texprog) and Tim-3+ terminally exhausted (Texterm) subpopulations. Inhibitor of DNA binding protein 2 (Id2) has been shown to play important roles in T-cell development and CD8+ T-cell immunity. However, the role of Id2 in CD8+ T-cell exhaustion is unclear. Here, we found that Id2 transcriptionally and epigenetically regulates the generation of Texprog cells and their conversion to Texterm cells. Genetic deletion of Id2 dampens CD8+ T-cell-mediated immune responses and the maintenance of stem-like CD8+ T-cell subpopulations, suppresses PD-1 blockade and increases tumor susceptibility. Mechanistically, through its HLH domain, Id2 binds and disrupts the assembly of the Tcf3-Tal1 transcriptional regulatory complex, and thus modulates chromatin accessibility at the Slamf6 promoter by preventing the interaction of Tcf3 with the histone lysine demethylase LSD1. Therefore, Id2 increases the abundance of the permissive H3K4me2 mark on the Tcf3-occupied E-boxes in the Slamf6 promoter, modulates chromatin accessibility at the Slamf6 promoter and epigenetically regulates the generation of Slamf6+ Texprog cells. An LSD1 inhibitor GSK2879552 can rescue the Id2 knockout phenotype in tumor-bearing mice. Inhibition of LSD1 increases the abundance of Slamf6+Tim-3- Texprog cells in tumors and the expression level of Tcf1 in Id2-deleted CD8+ T cells. This study demonstrates that Id2-mediated transcriptional and epigenetic modification drives hierarchical CD8+ T-cell exhaustion, and the mechanistic insights gained may have implications for therapeutic intervention with tumor immune evasion.

Hypoxia-induced immortalization of primary cells depends on Tfcp2L1 expression

Cellular senescence is a stress response mechanism that induces proliferative arrest. Hypoxia can bypass senescence and extend the lifespan of primary cells, mainly by decreasing oxidative damage. However, how hypoxia promotes these effects prior to malignant transformation is unknown. Here we observed that the lifespan of mouse embryonic fibroblasts (MEFs) is increased when they are cultured in hypoxia by reducing the expression of p16INK4a, p15INK4b and p21Cip1. We found that proliferating MEFs in hypoxia overexpress Tfcp2l1, which is a main regulator of pluripotency and self-renewal in embryonic stem cells, as well as stemness genes including Oct3/4, Sox2 and Nanog. Tfcp2l1 expression is lost during culture in normoxia, and its expression in hypoxia is regulated by Hif1α. Consistently, its overexpression in hypoxic levels increases the lifespan of MEFs and promotes the overexpression of stemness genes. ATAC-seq and Chip-seq experiments showed that Tfcp2l1 regulates genes that control proliferation and stemness such as Sox2, Sox9, Jarid2 and Ezh2. Additionally, Tfcp2l1 can replicate the hypoxic effect of increasing cellular reprogramming. Altogether, our data suggest that the activation of Tfcp2l1 by hypoxia contributes to immortalization prior to malignant transformation, facilitating tumorigenesis and dedifferentiation by regulating Sox2, Sox9, and Jarid2.

CHOP-mediated Gasdermin E expression promotes pyroptosis, inflammation, and mitochondrial damage in renal ischemia-reperfusion injury

In clinical practice, renal ischemia-reperfusion injury (IRI) is a common cause of acute kidney injury (AKI), often leading to acute renal failure or end-stage renal disease (ESRD). The current understanding of renal IRI mechanisms remains unclear, and effective therapeutic strategies and clear targets are lacking. Therefore, the need to find explicit and effective ways to reduce renal IRI remains a scientific challenge. The current study explored pyroptosis, a type of inflammation-regulated programmed cell death, and the role of Gasdermins E (GSDME)-mediated pyroptosis, mitochondrial damage, and inflammation in renal IRI. The analysis of human samples showed that the expression levels of GSDME in normal human renal tissues were higher than those of GSDMD. Moreover, our study demonstrated that GSDME played an important role in mediating pyroptosis, inflammation, and mitochondrial damage in renal IRI. Subsequently, GSDME-N accumulated in the mitochondrial membrane, leading to mitochondrial damage and activation of caspase3, which generated a feed-forward loop of self-amplification injury. However, GSDME knockout resulted in the amelioration of renal IRI. Moreover, the current study found that the transcription factor CHOP was activated much earlier in renal IRI. Inhibition of BCL-2 by CHOP leaded to casapse3 activation, resulting in mitochondrial damage and apoptosis; not only that, but CHOP positively regulated GSDME thereby causing pyroptosis. Therefore, this study explored the transcriptional mechanisms of GSDME during IRI development and the important role of CHOP/Caspase3/GSDME mechanistic axis in regulating pyroptosis in renal IRI. This axis might serve as a potential therapeutic target.

Cholesterol lowering depletes atherosclerotic lesions of smooth muscle cell-derived fibromyocytes and chondromyocytes

Drugs that lower plasma apolipoprotein B (ApoB)-containing lipoproteins are central to treating advanced atherosclerosis and provide partial protection against clinical events. Previous research showed that lowering ApoB-containing lipoproteins stops plaque inflammation, but how these drugs affect the heterogeneous population of plaque cells derived from smooth muscle cells (SMCs) is unknown. SMC-derived cells are the main cellular component of atherosclerotic lesions and the source of structural components that determine the size of plaques and their propensity to rupture and trigger thrombosis, the proximate cause of heart attack and stroke. Using lineage tracing and single-cell techniques to investigate the full SMC-derived cellular compartment in progressing and regressing plaques in mice, here we show that lowering ApoB-containing lipoproteins reduces nuclear factor kappa-light-chain-enhancer of activated B cells signaling in SMC-derived fibromyocytes and chondromyocytes and leads to depletion of these abundant cell types from plaques. These results uncover an important mechanism through which cholesterol-lowering drugs can achieve plaque regression.

A type 1 immunity-restricted promoter of the IL-33 receptor gene directs antiviral T-cell responses

The pleiotropic alarmin interleukin-33 (IL-33) drives type 1, type 2 and regulatory T-cell responses via its receptor ST2. Subset-specific differences in ST2 expression intensity and dynamics suggest that transcriptional regulation is key in orchestrating the context-dependent activity of IL-33-ST2 signaling in T-cell immunity. Here, we identify a previously unrecognized alternative promoter in mice and humans that is located far upstream of the curated ST2-coding gene and drives ST2 expression in type 1 immunity. Mice lacking this promoter exhibit a selective loss of ST2 expression in type 1- but not type 2-biased T cells, resulting in impaired expansion of cytotoxic T cells (CTLs) and T-helper 1 cells upon viral infection. T-cell-intrinsic IL-33 signaling via type 1 promoter-driven ST2 is critical to generate a clonally diverse population of antiviral short-lived effector CTLs. Thus, lineage-specific alternative promoter usage directs alarmin responsiveness in T-cell subsets and offers opportunities for immune cell-specific targeting of the IL-33-ST2 axis in infections and inflammatory diseases.

Nucleus pulposus cells regulate macrophages in degenerated intervertebral discs via the integrated stress response-mediated CCL2/7-CCR2 signaling pathway

Lower back pain (LBP), which is a primary cause of disability, is largely attributed to intervertebral disc degeneration (IDD). Macrophages (MΦs) in degenerated intervertebral discs (IVDs) form a chronic inflammatory microenvironment, but how MΦs are recruited to degenerative segments and transform into a proinflammatory phenotype remains unclear. We evaluated chemokine expression in degenerated nucleus pulposus cells (NPCs) to clarify the role of NPCs in the establishment of an inflammatory microenvironment in IDD and explored the mechanisms. We found that the production of C-C motif chemokine ligand 2 (CCL2) and C-C motif chemokine ligand 7 (CCL7) was significantly increased in NPCs under inflammatory conditions, and blocking CCL2/7 and their receptor, C-C chemokine receptor type 2(CCR2), inhibited the inductive effects of NPCs on MΦ infiltration and proinflammatory polarization. Moreover, activation of the integrated stress response (ISR) was obvious in IDD, and ISR inhibition reduced the production of CCL2/7 in NPCs. Further investigation revealed that activating Transcription Factor 3 (ATF3) responded to ISR activation, and ChIP-qPCR verified the DNA-binding activity of ATF3 on CCL2/7 promoters. In addition, we found that Toll-like receptor 4 (TLR4) inhibition modulated ISR activation, and TLR4 regulated the accumulation of mitochondrial reactive oxygen species (mtROS) and double-stranded RNA (dsRNA). Downregulating the level of mtROS reduced the amount of dsRNA and ISR activation. Deactivating the ISR or blocking CCL2/7 release alleviated inflammation and the progression of IDD in vivo. Moreover, MΦ infiltration and IDD were inhibited in CCR2-knockout mice. In conclusion, this study highlights the critical role of TLR4/mtROS/dsRNA axis-mediated ISR activation in the production of CCL2/7 and the progression of IDD, which provides promising therapeutic strategies for discogenic LBP.

The nuclear factor ID3 endows macrophages with a potent anti-tumour activity

Macrophage activation is controlled by a balance between activating and inhibitory receptors1-7, which protect normal tissues from excessive damage during infection8,9 but promote tumour growth and metastasis in cancer7,10. Here we report that the Kupffer cell lineage-determining factor ID3 controls this balance and selectively endows Kupffer cells with the ability to phagocytose live tumour cells and orchestrate the recruitment, proliferation and activation of natural killer and CD8 T lymphoid effector cells in the liver to restrict the growth of a variety of tumours. ID3 shifts the macrophage inhibitory/activating receptor balance to promote the phagocytic and lymphoid response, at least in part by buffering the binding of the transcription factors ELK1 and E2A at the SIRPA locus. Furthermore, loss- and gain-of-function experiments demonstrate that ID3 is sufficient to confer this potent anti-tumour activity to mouse bone-marrow-derived macrophages and human induced pluripotent stem-cell-derived macrophages. Expression of ID3 is therefore necessary and sufficient to endow macrophages with the ability to form an efficient anti-tumour niche, which could be harnessed for cell therapy in cancer.

Single-nucleus multi-omic profiling of human placental syncytiotrophoblasts identifies cellular trajectories during pregnancy

The human placenta has a vital role in ensuring a successful pregnancy. Despite the growing body of knowledge about its cellular compositions and functions, there has been limited research on the heterogeneity of the billions of nuclei within the syncytiotrophoblast (STB), a multinucleated entity primarily responsible for placental function. Here we conducted integrated single-nucleus RNA sequencing and single-nucleus ATAC sequencing analyses of human placentas from early and late pregnancy. Our findings demonstrate the dynamic heterogeneity and developmental trajectories of STB nuclei and their correspondence with human trophoblast stem cell (hTSC)-derived STB. Furthermore, we identified transcription factors associated with diverse STB nuclear lineages through their gene regulatory networks and experimentally confirmed their function in hTSC and trophoblast organoid-derived STBs. Together, our data provide insights into the heterogeneity of human STB and represent a valuable resource for interpreting associated pregnancy complications.

Single-cell multiomics decodes regulatory programs for mouse secondary palate development

Perturbations in gene regulation during palatogenesis can lead to cleft palate, which is among the most common congenital birth defects. Here, we perform single-cell multiome sequencing and profile chromatin accessibility and gene expression simultaneously within the same cells (n = 36,154) isolated from mouse secondary palate across embryonic days (E) 12.5, E13.5, E14.0, and E14.5. We construct five trajectories representing continuous differentiation of cranial neural crest-derived multipotent cells into distinct lineages. By linking open chromatin signals to gene expression changes, we characterize the underlying lineage-determining transcription factors. In silico perturbation analysis identifies transcription factors SHOX2 and MEOX2 as important regulators of the development of the anterior and posterior palate, respectively. In conclusion, our study charts epigenetic and transcriptional dynamics in palatogenesis, serving as a valuable resource for further cleft palate research.