Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position

Novel Approach to Overcome Osimertinib Resistance Using Bromodomain and Extra-Terminal Domain Inhibitors

Osimertinib, a third-generation EGFR-tyrosine kinase inhibitor, is the first-line therapy for lung cancer harboring EGFR mutations. The mechanisms underlying osimertinib resistance are diverse, with approximately half remaining unknown. Epigenetic dysregulation is implicated in drug resistance; however, the mechanisms remain unclear. Therefore, we investigated epigenetic involvement in osimertinib resistance and its therapeutic potential. We established osimertinib-resistant cells and used an assay for transposase-accessible chromatin using sequencing to evaluate chromatin accessibility, finding significant changes post-resistance. Combining the assay for transposase-accessible chromatin and RNA sequencing data, we identified FGF1 as a resistance-related gene regulated by histone modifications. FGF1 induced osimertinib resistance, and its suppression attenuated resistance. Bromodomain and extra-terminal domain inhibitors combined with osimertinib overcame osimertinib resistance by reducing FGF1 expression. Increased FGF1 expression was observed in osimertinib-resistant clinical samples. This combination therapy was effective in cell lines and mouse xenograft models. These results suggest targeting histone modifications using bromodomain and extra-terminal domain inhibitors as a novel approach to overcoming osimertinib resistance.

Matrix-producing neutrophils populate and shield the skin

Defence from environmental threats is provided by physical barriers that confer mechanical protection and prevent the entry of microorganisms1. If microorganisms overcome those barriers, however, innate immune cells use toxic chemicals to kill the invading cells2,3. Here we examine immune diversity across tissues and identify a population of neutrophils in the skin that expresses a broad repertoire of proteins and enzymes needed to build the extracellular matrix. In the naive skin, these matrix-producing neutrophils contribute to the composition and structure of the extracellular matrix, reinforce its mechanical properties and promote barrier function. After injury, these neutrophils build 'rings' of matrix around wounds, which shield against foreign molecules and bacteria. This structural program relies on TGFβ signalling; disabling the TGFβ receptor in neutrophils impaired ring formation around wounds and facilitated bacterial invasion. We infer that the innate immune system has evolved diverse strategies for defence, including one that physically shields the host from the outside world.

Chromatin Accessibility and Translational Landscapes of Bermudagrass (Cynodon dactylon) Under Chilling Stress

Cold tolerance in plants is a complex trait regulated by a network of transcription factors (TFs) and their downstream genes. While C-repeat binding factors (CBFs) are well-known for their role in cold tolerance, other regulatory networks remain largely unexplored. This study utilizes a combined approach of Assay for Transposase-Accessible Chromatin with Sequencing (ATAC-seq) and RNA-seq to identify cold-responsive TFs in Cynodon dactylon (Bermudagrass), a species widely grown in southern China but limited by cold stress. Early-stage cold stress was found to induce dynamic changes in 331 differentially accessible regions (DARs), with over 45% located in gene promoter regions. Key TFs, including CAMTA1, CAMTA2, WRKY43, WRKY48, WRKY21, and DREB1G, were associated with gained DARs, highlighting their roles in cold response regulation. In contrast, several Heat Shock Factor (HSF) family members, such as HSFA6B, HSFB2A, HSFC1, and HSFB2B, were linked to lost DARs, underscoring their involvement in cold stress regulation. The correlation between chromatin accessibility and gene expression emphasizes the critical role of TFs in plant cold stress adaptation. This study highlights the potential of ATAC-seq and RNA-seq as powerful tools for uncovering novel TFs essential for cold tolerance in plants, offering insights for future research and breeding efforts.

Tri-omic single-cell mapping of the 3D epigenome and transcriptome in whole mouse brains throughout the lifespan

Exploring the genomic basis of transcriptional programs has been a long-standing research focus. Here we report a single-cell method, ChAIR, to map chromatin accessibility, chromatin interactions and RNA expression simultaneously. After validating in cultured cells, we applied ChAIR to whole mouse brains and delineated the concerted dynamics of epigenome, three-dimensional (3D) genome and transcriptome during maturation and aging. In particular, gene-centric chromatin interactions and open chromatin states provided 3D epigenomic mechanism underlying cell-type-specific transcription and revealed spatially resolved specificity. Importantly, the composition of short-range and ultralong chromatin contacts in individual cells is remarkably correlated with transcriptional activity, open chromatin state and genome folding density. This genomic property, along with associated cellular properties, differs in neurons and non-neuronal cells across different anatomic regions throughout the lifespan, implying divergent nuclear mechano-genomic mechanisms at play in brain cells. Our results demonstrate ChAIR's robustness in revealing single-cell 3D epigenomic states of cell-type-specific transcription in complex tissues.

Advances in Spatial Omics Technologies

Rapidly developing spatial omics technologies provide us with new approaches to deeply understanding the diversity and functions of cell types within organisms. Unlike traditional approaches, spatial omics technologies enable researchers to dissect the complex relationships between tissue structure and function at the cellular or even subcellular level. The application of spatial omics technologies provides new perspectives on key biological processes such as nervous system development, organ development, and tumor microenvironment. This review focuses on the advancements and strategies of spatial omics technologies, summarizes their applications in biomedical research, and highlights the power of spatial omics technologies in advancing the understanding of life sciences related to development and disease.

Glucose modulates IRF6 transcription factor dimerization to enable epidermal differentiation

Non-energetic roles for glucose are largely unclear, as is the interplay between transcription factors (TFs) and ubiquitous biomolecules. Metabolomic analyses uncovered elevation of intracellular glucose during differentiation of diverse cell types. Human and mouse tissue engineered with glucose sensors detected a glucose gradient that peaked in the outermost differentiated layers of the epidermis. Free glucose accumulation was essential for epidermal differentiation and required the SGLT1 glucose transporter. Glucose affinity chromatography uncovered glucose binding to diverse regulatory proteins, including the IRF6 TF. Direct glucose binding enabled IRF6 dimerization, DNA binding, genomic localization, and induction of IRF6 target genes, including essential pro-differentiation TFs GRHL1, GRHL3, HOPX, and PRDM1. These data identify a role for glucose as a gradient morphogen that modulates protein multimerization in cellular differentiation.

An ISWI-related chromatin remodeller regulates stage-specific gene expression in Toxoplasma gondii

ATP-dependent chromatin remodellers are specialized multiprotein machines that organize the genome in eukaryotic cells and regulate its accessibility by repositioning, ejecting or modifying nucleosomes. However, their role in Toxoplasma gondii is poorly understood. Here we show that T. gondii has evolved two divergent proteins within the imitation switch (ISWI) family: TgSNF2h and TgSNF2L. TgSNF2h specifically forms a core complex with the transcription factor AP2VIII-2 and the scaffold protein TgRFTS. Depletion of TgRFTS phenocopies the knockdown of TgSNF2h, restricting access to chromatin and altering local gene expression. At the genomic level, TgSNF2h insulates highly transcribed genes from silenced neighbours, ensuring stage-specific gene regulation. By modulating chromatin accessibility to transcription factors, TgSNF2h exerts epistatic control over MORC, a key regulator of sexual commitment. Our findings show that a specific ISWI complex orchestrates the partitioning of developmental genes and ensures transcriptional fidelity throughout the parasite life cycle.

Heterogeneity of the tumor immune cell microenvironment revealed by single-cell sequencing in head and neck cancer

Head and neck cancer (HNC) is the sixth most common disease in the world. The recurrence rate of patients is relatively high, and the heterogeneity of tumor immune microenvironment (TIME) cells may be an important reason for this. Single-cell sequencing (SCS) is currently the most promising and mature application in cancer research. It can identify unique genes expressed in cells and study tumor heterogeneity. According to current research, the heterogeneity of immune cells has become an important factor affecting the occurrence and development of HNC. SCSs can provide effective therapeutic targets and prognostic factors for HNC patients through analyses of gene expression levels and cell heterogeneity. Therefore, this study analyzes the basic theory of HNC and the development of SCS technology, elaborating on the application of SCS technology in HNC and its potential value in identifying HNC therapeutic targets and biomarkers.

Single-cell parallel analysis of DNA damage and transcriptome reveals selective genome vulnerability

Maintenance of genome integrity is paramount to molecular programs in multicellular organisms. Throughout the lifespan, various endogenous and environmental factors pose persistent threats to the genome, which can result in DNA damage. Understanding the functional consequences of DNA damage requires investigating their preferred genomic distributions and influences on gene regulatory programs. However, such analysis is hindered by both the complex cell-type compositions within organs and the high background levels due to the stochasticity of damage formation. To address these challenges, we developed Paired-Damage-seq for joint analysis of oxidative and single-stranded DNA damage with gene expression in single cells. We applied this approach to cultured HeLa cells and the mouse brain as a proof of concept. Our results indicated the associations between damage formation and epigenetic changes. The distribution of oxidative DNA damage hotspots exhibits cell-type-specific patterns; this selective genome vulnerability, in turn, can predict cell types and dysregulated molecular programs that contribute to disease risks.

Multipronged SMAD pathway targeting by lipophilic poly(β-amino ester) miR-590-3p nanomiRs inhibits mesenchymal glioblastoma growth and prolongs survival

ASBSTRACT

Despite aggressive therapy, glioblastoma (GBM) recurs in almost all patients and treatment options are very limited. Despite our growing understanding of how cellular transitions associate with relapse in GBM, critical gaps remain in our ability to block these molecular changes and treat recurrent disease. In this study we combine computational biology, forward-thinking understanding of miRNA biology and cutting-edge nucleic acid delivery vehicles to advance targeted therapeutics for GBM. Computational analysis of RNA sequencing from clinical GBM specimens identified TGFβ type II receptor (TGFBR2) as a key player in the mesenchymal transition associated with worse outcome in GBM. Mechanistically, we show that elevated levels of TGFBR2 is conducive to reduced temozolomide (TMZ) sensitivity. This effect is, at least partially, induced by stem-cell driving events coordinated by the reprogramming transcription factors Oct4 and Sox2 that lead to open chromatin states. We show that blocking TGFBR2 via molecular and pharmacological approaches decreases stem cell capacity and sensitivity of clinical recurrent GBM (rGBM) isolates to TMZ in vitro. Network analysis uncovered miR-590-3p as a tumor suppressor that simultaneously inhibits multiple oncogenic nodes downstream of TGFBR2. We also developed novel biodegradable lipophilic poly(β-amino ester) nanoparticles (LiPBAEs) for in vivo microRNA (miRNAs) delivery. Following direct intra-tumoral infusion, these nanomiRs efficiently distribute through the tumors. Importantly, miR-590-3p nanomiRs inhibited the growth and extended survival of mice bearing orthotopic human rGBM xenografts, with an apparent 30% cure rate. These results show that miRNA-based targeted therapeutics provide new opportunities to treat rGBM and bypass the resistance to standard of care therapy.

Rewiring of SINE-MIR enhancer topology and Esrrb modulation in expanded and naive pluripotency

BACKGROUND

The interplay between 3D genomic structure and transposable elements (TE) in regulating cell state-specific gene expression program is largely unknown. Here, we explore the utilization of TE-derived enhancers in naïve and expanded pluripotent states by integrative analysis of genome-wide Hi-C-defined enhancer interactions, H3K27ac HiChIP profiling and CRISPR-guided TE proteomics landscape.

RESULTS

We find that short interspersed nuclear elements (SINEs) are the more involved TEs in the active chromatin and 3D genome architecture. In particular, mammalian-wide interspersed repeat (MIR), a SINE family member, is highly associated with naïve-specific genomic interactions compared to the expanded state. Primarily, in the naïve pluripotent state, MIR enhancer is co-opted by ESRRB for naïve-specific gene expression program. This ESRRB and MIR enhancer interaction is crucial for the formation of loops that build a network of enhancers and super-enhancers regulating pluripotency genes. We demonstrate that loss of a ESRRB-bound MIR enhancer impairs self-renewal. We also find that MIR is co-bound by structural protein complex, ESRRB-YY1, in the naïve pluripotent state.

CONCLUSIONS

Altogether, our study highlights the topological regulation of ESRRB on MIR in the naïve potency state.

Chromatin profiling to identify biomarkers in inflammatory bowel diseases

PURPOSE OF REVIEW

Chromatin plays a critical role in gene regulation and disease pathogenesis. In inflammatory bowel disease (IBD), alterations in chromatin structure contribute to disease heterogeneity and impact treatment responses. This review explores chromatin accessibility and chromatin-associated proteins as biomarkers for IBD and highlights recent technological advancements enabling targeted biomarker discovery and novel therapies.

RECENT FINDINGS

Advancements in high-throughput sequencing have enabled genome-wide profiling of chromatin interactions in IBD. Studies have identified distinct chromatin landscapes in Crohn's disease (CD) and ulcerative colitis (UC), revealing stable regulatory shifts independent of inflammation.

SUMMARY

Chromatin profiling offers a novel approach for identifying biomarkers and therapeutic targets in IBD. Integrating chromatin accessibility data with transcriptomic and epigenomic analyses can refine disease classification and guide personalized treatment strategies. Emerging techniques compatible with formalin-fixed paraffin-embedded (FFPE) samples enhance clinical applicability, bridging the gap between molecular research and precision gastroenterology.

Altering metabolism programs cell identity via NAD+-dependent deacetylation

Cells change their metabolic profiles in response to underlying gene regulatory networks, but how can alterations in metabolism encode specific transcriptional instructions? Here, we show that forcing a metabolic change in embryonic stem cells (ESCs) promotes a developmental identity that better approximates the inner cell mass (ICM) of the early mammalian blastocyst in cultures. This shift in cellular identity depends on the inhibition of glycolysis and stimulation of oxidative phosphorylation (OXPHOS) triggered by the replacement of D-glucose by D-galactose in ESC media. Enhanced OXPHOS in turn activates NAD + -dependent deacetylases of the Sirtuin family, resulting in the deacetylation of histones and key transcription factors to focus enhancer activity while reducing transcriptional noise, which results in a robustly enhanced ESC phenotype. This exploitation of a NAD + /NADH coenzyme coupled to OXPHOS as a means of programming lineage-specific transcription suggests new paradigms for how cells respond to alterations in their environment, and implies cellular rejuvenation exploits enzymatic activities for simultaneous activation of a discrete enhancer set alongside silencing genome-wide transcriptional noise.

IKAROS levels are associated with antigen escape in CD19- and CD22-targeted therapies for B-cell malignancies

Antigen escape relapse is a major challenge in targeted immunotherapies, including CD19- and CD22-directed chimeric antigen receptor (CAR) T-cell for B-cell acute lymphoblastic leukemia (B-ALL). To identify tumor-intrinsic factors driving antigen loss, we perform single-cell analyses on 61 B-ALL patient samples treated with CAR T cells. Here we show that low levels of IKAROS in pro-B-like B-ALL cells before CAR T treatment correlate with antigen escape. IKAROSlow B-ALL cells undergo epigenetic and transcriptional changes that diminish B-cell identity, making them resemble progenitor cells. This shift leads to reduced CD19 and CD22 surface expression. We demonstrate that CD19 and CD22 expression is IKAROS dose-dependent and reversible. Furthermore, IKAROSlow cells exhibit higher resistance to CD19- and CD22-targeted therapies. These findings establish a role for IKAROS as a regulator of antigens targeted by widely used immunotherapies and in the risk of antigen escape relapse, identifying it as a potential prognostic target.

Polymer model integrates imaging and sequencing to reveal how nanoscale heterochromatin domains influence gene expression

Chromatin organization regulates gene expression, with nanoscale heterochromatin domains playing a fundamental role. Their size varies with microenvironmental stiffness and epigenetic interventions, but how these factors regulate their formation and influence transcription remains unclear. To address this, we developed a sequencing-informed copolymer model that simulates chromatin evolution through diffusion and active epigenetic reactions. Our model predicts the formation of nanoscale heterochromatin domains and quantifies how domain size scales with epigenetic reaction rates, showing that epigenetic and compaction changes primarily occur at domain boundaries. We validated these predictions via Hi-C and super-resolution imaging of hyperacetylated melanoma cells and identified differential expression of metastasis-related genes through RNA-seq. We validated our findings in hMSCs, where epigenetic reaction rates respond to microenvironmental stiffness. Conclusively, our simulations reveal that heterochromatin domain boundaries regulate gene expression and epigenetic memory. These findings demonstrate how external cues drive chromatin organization and transcriptional memory in development and disease.

Integrative multi-omics data provide insights into the biosynthesis of furanocoumarins and mechanisms regulating their accumulation in Angelica dahurica

Furocoumarins (FCs), important natural compounds with biodefense roles and pharmacological activities, are notably abundant in medicinal plant Angelica dahurica. However, its accumulation patterns over development stages in FC-enriched tissue, biosynthetic pathways, and regulatory mechanisms in A. dahurica remain elusive. Here, we quantified the concentration dynamics of 17 coumarins across six developmental stages of root and found a gradual decrease in FC concentration as the roots develop. Using a de-novo assembled chromosome-level genome for A. dahurica, we conducted integrative multi-omics analyses to screen out candidate genes to fill in the sole missing step in the biosynthesis of imperatorin and isoimperatorin. This revealed that CYP71AZ18 catalyzes hydroxylation at the C-5 position of psoralen to generate bergaptol, while CYP71AZ19 and CYP83F95 catalyze hydroxylation at the C-8 position to produce xanthotoxol, notably indicating that a single step is catalyzed by two genes from distinct CYP450 subfamilies in this species. CYP71AZ19 originated from a proximal duplication event of CYP71AZ18, specific to A. dahurica, and subsequently underwent neofunctionalization. Accessible chromatin regions (ACRs), especially proximal ACRs, correlated with high gene expression levels, and the three validated genes exhibited strong signals of ACRs, showing the importance of chromosomal accessibility in regulating metabolite biosynthesis.

Orchestration of pluripotent stem cell genome reactivation during mitotic exit

Cell identity maintenance faces many challenges during mitosis, as most DNA-binding proteins are evicted from DNA and transcription is virtually abolished. How cells maintain their identity through division and faithfully re-initiate gene expression during mitotic exit is unclear. Here, we develop a novel reporter system enabling cell cycle synchronization-free separation of pluripotent stem cells in temporal bins of <30 min during mitotic exit. This allows us to quantify genome-wide reactivation of transcription, sequential changes in chromatin accessibility and transcription factor footprints, and re-binding of the pluripotency transcription factors OCT4, SOX2, and NANOG (OSN). We find that transcriptional activity progressively ramps up after mitosis and that OSN rapidly reoccupy the genome during the anaphase-telophase transition. We also demonstrate transcription factor-specific, dynamic relocation patterns and a hierarchical reorganization of the OSN binding landscape governed by OCT4 and SOX2. Our study sheds light on the dynamic orchestration of transcriptional reactivation after mitosis.

Maternal transcription factor OTX2 directly induces SETD1A and promotes embryonic genome activation in human pre-implantation embryos

Early embryonic development is controlled by maternal factors originating from mature oocytes. The zygotic genome is activated from a transcriptionally quiescent state through a process called embryonic genome activation (EGA), which involves the depletion and clearance of maternal factors. However, the mechanism by which maternal factors regulate EGA and embryonic development, particularly in humans, remains elusive. In this study, using tri-pronuclear (3PN) embryos and human embryonic stem cells (hESCs), we demonstrated that the maternal transcription factor Orthodenticle Homeobox 2 (OTX2), a paired-like homeobox gene, promotes EGA in human pre-implantation embryos. Knockdown of OTX2 through Trim-Away technology blocked embryonic development and minor EGA gene expression. Overexpression of OTX2 (OTX2OE) in hESCs increased transcript products, primarily at the 2-cell embryo stage genes, including genes encoding methyltransferase of histone H3K4. OTX2OE increased the level of H3K4me3 and increased the open chromatin region that co-occurs with the H3K4me3 region at the 4-cell stage in hESCs. Based on these findings in hESCs, we further verified that OTX2 directly induced the expression of SETD1A by binding to its promoter, leading to increased H3K4me3 levels in both hESCs and 3PN embryos. These findings suggest that the maternal transcription factor OTX2 regulates EGA and early embryogenesis via epigenetic mechanisms.

DUX-mediated configuration of p300/CBP drives minor zygotic genome activation independent of its catalytic activity

Maternal-deposited factors initiate zygotic genome activation (ZGA), driving the maternal-to-zygotic transition; however, the coordination between maternal coactivators and transcription factors (TFs) in this process remains unclear. In this study, by profiling the dynamic landscape of p300 during mouse ZGA, we reveal its role in promoting RNA polymerase II (Pol II) pre-configuration at ZGA gene regions and sequentially establishing enhancer activity and regulatory networks. Moreover, p300/CBP-catalyzed acetylation drives Pol II elongation and minor ZGA gene expression by inducing pivotal TFs such as Dux. Remarkably, the supplementation of exogenous Dux rescues ZGA failure and developmental defects caused by the loss of p300/CBP acetylation. DUX functions as a pioneer factor, guiding p300 and Pol II to minor ZGA gene regions and activating them in a manner dependent on the non-catalytic functions of p300/CBP. Together, our findings reveal a mutual dependency between p300/CBP and DUX, highlighting their coordinated role in regulating minor ZGA activation.

Microglia endotoxin tolerance is retained after enforced repopulation

Microglia are crucial for CNS homeostasis and are involved in a wide range of neurodegenerative and neuroinflammatory diseases. Systemic inflammation and infections can contribute to neurodegeneration later in life by affecting microglia. Like other innate immune cells, microglia can develop innate immune memory (IIM) in response to an inflammatory challenge, altering their response to subsequent stimuli. IIM can ameliorate or worsen CNS pathology, but it is unclear if IIM can be reversed to restore microglia functions. Here, we investigated whether microglia depletion-repopulation by inhibition of the colony-stimulating factor 1 receptor with BLZ945 reversed LPS-induced microglia endotoxin tolerance in mice. Repopulated microglia displayed a reduced expression of homeostatic genes and genes related to mitochondrial respiration and TCA cycle metabolism and an increased expression of immune effector and activation genes. Nonetheless, the blunted inflammatory gene response after LPS-preconditioning was retained after a depletion-repopulation cycle. Our study highlights the persistence of endotoxin tolerance in microglia after a depletion-repopulation cycle, which might impact the potential effectiveness of strategies targeted at microglia depletion for clinical applications.

Unified molecular approach for spatial epigenome, transcriptome, and cell lineages

Spatial epigenomics and multiomics can provide fine-grained insights into cellular states but their widespread adoption is limited by the requirement for bespoke slides and capture chemistries for each data modality. Here, we present SPatial assay for Accessible chromatin, Cell lineages, and gene Expression with sequencing (SPACE-seq), a method that utilizes polyadenine-tailed epigenomic libraries to enable facile spatial multiomics using standard whole transcriptome reagents. Applying SPACE-seq to a human glioblastoma specimen, we reveal the state of the tumor microenvironment, extrachromosomal DNA copy numbers, and identify putative mitochondrial DNA variants.

Cathelicidin antimicrobial peptides mediate immune protection in marsupial neonates

Marsupial neonates are born with immature immune systems, making them vulnerable to pathogens. While neonates receive maternal protection, they can also independently combat pathogens, although the mechanisms remain unknown. Using the sugar glider (Petaurus breviceps) as a model, we investigated immunological defense strategies of marsupial neonates. Cathelicidins-a family of antimicrobial peptides expanded in the genomes of marsupials-are highly expressed in developing neutrophils. Sugar glider cathelicidins reside in two genomic clusters, and their coordinated expression is achieved by enhancer sharing within clusters and long-range physical interactions between clusters. Functionally, cathelicidins modulate immune responses and have potent antibacterial effects, sufficient to provide protection in a mouse model of sepsis. Evolutionarily, cathelicidins have a complex history, with marsupials and monotremes uniquely retaining both clusters among tetrapods. Thus, cathelicidins are critical mediators of marsupial immunity, and their evolution may reflect the life history-specific immunological needs of these animals.

A Sox2 enhancer cluster regulates region-specific neural fates from mouse embryonic stem cells

Sex-determining region Y box 2 (Sox2) is a critical transcription factor for embryogenesis and neural stem and progenitor cell (NSPC) maintenance. While distal enhancers control Sox2 in embryonic stem cells (ESCs), enhancers closer to the gene are implicated in Sox2 transcriptional regulation in neural development. We hypothesize that a downstream enhancer cluster, termed Sox2 regulatory regions 2-18 (SRR2-18), regulates Sox2 transcription in neural stem cells and we investigate this in NSPCs derived from mouse ESCs. Using functional genomics and CRISPR-Cas9-mediated deletion analyses, we investigate the role of SRR2-18 in Sox2 regulation during neural differentiation. Transcriptome analyses demonstrate that the loss of even 1 copy of SRR2-18 disrupts the region-specific identity of NSPCs, reducing the expression of genes associated with more anterior regions of the embryonic nervous system. Homozygous deletion of this Sox2 neural enhancer cluster causes reduced SOX2 protein, less frequent interaction with transcriptional machinery, and leads to perturbed chromatin accessibility genome-wide further affecting the expression of neurodevelopmental and anterior-posterior regionalization genes. Furthermore, homozygous NSPC deletants exhibit self-renewal defects and impaired differentiation into cell types found in the brain. Altogether, our data define a cis-regulatory enhancer cluster controlling Sox2 transcription in NSPCs and highlight the sensitivity of neural differentiation processes to decreased Sox2 transcription, which causes differentiation into posterior neural fates, specifically the caudal neural tube. This study highlights the importance of precise Sox2 regulation by SRR2-18 in neural differentiation.

Hmgb2 improves astrocyte to neuron conversion by increasing the chromatin accessibility of genes associated with neuronal maturation in a proneuronal factor-dependent manner

BACKGROUND

Direct conversion of reactive glial cells to neurons is a promising avenue for neuronal replacement therapies after brain injury or neurodegeneration. The overexpression of neurogenic fate determinants in glial cells results in conversion to neurons. For repair purposes, the conversion should ideally be induced in the pathology-induced neuroinflammatory environment. However, very little is known regarding the influence of the injury-induced neuroinflammatory environment and released growth factors on the direct conversion process.

RESULTS

We establish a new in vitro culture system of postnatal astrocytes without epidermal growth factor that reflects the direct conversion rate in the injured, neuroinflammatory environment in vivo. We demonstrate that the growth factor combination corresponding to the injured environment defines the ability of glia to be directly converted to neurons. Using this culture system, we show that chromatin structural protein high mobility group box 2 (HMGB2) regulates the direct conversion rate downstream of the growth factor combination. We further demonstrate that Hmgb2 cooperates with neurogenic fate determinants, such as Neurog2, in opening chromatin at the loci of genes regulating neuronal maturation and synapse formation. Consequently, early chromatin rearrangements occur during direct fate conversion and are necessary for full fate conversion.

CONCLUSIONS

Our data demonstrate novel growth factor-controlled regulation of gene expression during direct fate conversion. This regulation is crucial for proper maturation of induced neurons and could be targeted to improve the repair process.

A nucleosome switch primes hepatitis B virus infection

Chronic hepatitis B virus (HBV) infection is an incurable pathogen responsible for causing liver disease and hepatocellular carcinoma. During the genesis of infection, HBV establishes an independent minichromosome consisting of the viral covalently closed circular DNA (cccDNA) genome and host histones. The viral X gene must be expressed immediately upon infection to induce degradation of the host silencing factor, the Smc5/6 complex. However, the relationship between cccDNA chromatinization and X gene transcription remains poorly understood. By establishing a reconstituted viral minichromosome platform, we found that nucleosome occupancy in cccDNA regulates X transcription. We corroborated these findings in situ and further showed that the chromatin-destabilizing molecule CBL137 inhibits full-length X transcription and HBV infection in primary human hepatocytes. Our results shed light on a long-standing paradox and represent a potential therapeutic approach for the treatment of chronic HBV infection.

PRDM16-dependent antigen-presenting cells induce tolerance to gut antigens

The gastrointestinal tract is continuously exposed to foreign antigens in food and commensal microorganisms with potential to induce adaptive immune responses. Peripherally induced T regulatory (pTreg) cells are essential for mitigating inflammatory responses to these agents1-4. Although RORγt+ antigen-presenting cells (APCs) have been shown to programme gut microbiota-specific pTreg cells5-7, their definition remains incomplete, and the APC responsible for food tolerance has remained unknown. Here we identify an APC subset that is required for differentiation of both food- and microbiota-specific pTreg cells and for establishment of oral tolerance. Development and function of these APCs require expression of the transcription factors PRDM16 and RORγt, as well as a unique Rorc(t) cis-regulatory element. Gene expression, chromatin accessibility, and surface marker analysis establish the pTreg-inducing APCs as myeloid in origin, distinct from type 3 innate lymphoid cells, and sharing epigenetic profiles with classical dendritic cells, and designate them PRDM16+RORγt+ tolerizing dendritic cells (tolDCs). Upon genetic perturbation of tolDCs, we observe a substantial increase in food antigen-specific T helper 2 cells in lieu of pTreg cells, leading to compromised tolerance in mouse models of asthma and food allergy. Single-cell analyses of freshly resected mesenteric lymph nodes from a human organ donor, as well as multiple specimens of human intestine and tonsil, reveal candidate tolDCs with co-expression of PRDM16 and RORC and an extensive transcriptome shared with tolDCs from mice, highlighting an evolutionarily conserved role across species. Our findings suggest that a better understanding of how tolDCs develop and how they regulate T cell responses to food and microbial antigens could offer new insights into developing therapeutic strategies for autoimmune and allergic diseases as well as organ transplant tolerance.

Epigenomic signatures of cis-regulatory elements in the developing mouse and pig forelimb

Cis-regulatory elements (CREs) orchestrate the spatiotemporal regulation of key transcriptional programs. These genomic regions-including promoters, enhancers, and insulators-play a crucial role during embryonic development, and their functional diversification is thought to contribute significantly to the morphological evolution of animal body plans. We performed chromatin immunoprecipitation for various histone modifications (H3K4me3, H3K27ac, and H3K4me1) from equivalent developmental stages of mouse (E11.5) and pig (day 23; D23) to identify active regulatory regions during forelimb development in both species. The overlap of these epigenomic signatures with the pattern of open chromatin allowed us to classify these putative regulatory regions into different chromatin states in mouse and pig limb primordia. Our profiling of the regulatory genome in mouse and pig limb buds offers a valuable resource in the evo-devo field for exploring mechanisms underlying the morphological evolution of the tetrapod limb.

Inferring differential protein binding from time-series chromatin accessibility data

Motivation

Due to internal and external factors, the epigenomic landscape is constantly changing in ways that are linked to changes in gene expression. Chromatin accessibility data, such as MNase-seq, provide valuable insights into this landscape and have been used to compute chromatin occupancy profiles. Multiple datasets generated over time or under different conditions can thus be used to study dynamic changes in chromatin occupancy across the genome.

Results

Our existing model, RoboCOP, computes a genome-wide chromatin occupancy profile for nucleosomes and hundreds of transcription factors. Here, we present a new method called DynaCOP that takes multiple chromatin occupancy profiles and uses them to generate a series of nucleosome-guided difference profiles. These profiles identify differentially binding transcription factors and reveal changes in nucleosome occupancy and positioning. We apply DynaCOP to chromatin occupancy profiles derived from deeply sequenced time-series MNase-seq data to study differential chromatin occupancy in the yeast genome under cadmium stress. We find strong correlations between the observed chromatin changes and changes in transcription.

Availability and implementation

https://github.com/HarteminkLab/RoboCOP.

Characterization of non-coding variants associated with transcription-factor binding through ATAC-seq-defined footprint QTLs in liver

Non-coding variants discovered by genome-wide association studies (GWASs) are enriched in regulatory elements harboring transcription factor (TF) binding motifs, strongly suggesting a connection between disease association and the disruption of cis-regulatory sequences. Occupancy of a TF inside a region of open chromatin can be detected in ATAC-seq where bound TFs block the transposase Tn5, leaving a pattern of relatively depleted Tn5 insertions known as a "footprint." Here, we sought to identify variants associated with TF binding, or "footprint quantitative trait loci" (fpQTLs), in ATAC-seq data generated from 170 human liver samples. We used computational tools to scan the ATAC-seq reads to quantify TF binding likelihood as "footprint scores" at variants derived from whole-genome sequencing generated in the same samples. We tested for association between genotype and footprint score and observed 809 fpQTLs associated with footprint-inferred TF binding (FDR < 5%). Given that Tn5 insertion sites are measured with base-pair resolution, we show that fpQTLs can aid GWAS and QTL fine-mapping by precisely pinpointing TF activity within broad trait-associated loci where the underlying causal variant is unknown. Liver fpQTLs were strongly enriched across ChIP-seq peaks, liver expression QTLs (eQTLs), and liver-related GWAS loci, and their inferred effect on TF binding was concordant with their effect on underlying sequence motifs in 78% of cases. We conclude that fpQTLs can reveal causal GWAS variants, define the role of TF binding-site disruption in complex traits, and provide functional insights into non-coding variants, ultimately informing novel treatments for common diseases.

Interpreting regulatory mechanisms of Hippo signaling through a deep learning sequence model

Signaling pathway components are well studied, but how they mediate cell-type-specific transcription responses is an unresolved problem. Using the Hippo pathway in mouse trophoblast stem cells as a model, we show that the DNA binding of signaling effectors is driven by cell-type-specific sequence rules that can be learned genome wide by deep learning models. Through model interpretation and experimental validation, we show that motifs for the cell-type-specific transcription factor TFAP2C enhance TEAD4/YAP1 binding in a nucleosome-range and distance-dependent manner, driving synergistic enhancer activation. We also discovered that Tead double motifs are widespread, highly active canonical response elements. Molecular dynamics simulations suggest that TEAD4 binds them cooperatively through surprisingly labile protein-protein interactions that depend on the DNA template. These results show that the response to signaling pathways is encoded in the cis-regulatory sequences and that interpreting the rules reveals insights into the mechanisms by which signaling effectors influence cell-type-specific enhancer activity.

An in vivo systemic massively parallel platform for deciphering animal tissue-specific regulatory function

Introduction: Transcriptional regulation is an important process wherein non-protein coding enhancer sequences play a key role in determining cell type identity and phenotypic diversity. In neural tissue, these gene regulatory processes are crucial for coordinating a plethora of interconnected and regionally specialized cell types, ensuring their synchronized activity in generating behavior. Recognizing the intricate interplay of gene regulatory processes in the brain is imperative, as mounting evidence links neurodevelopment and neurological disorders to non-coding genome regions. While genome-wide association studies are swiftly identifying non-coding human disease-associated loci, decoding regulatory mechanisms is challenging due to causal variant ambiguity and their specific tissue impacts. Methods: Massively parallel reporter assays (MPRAs) are widely used in cell culture to study the non-coding enhancer regions, linking genome sequence differences to tissue-specific regulatory function. However, widespread use in animals encounters significant challenges, including insufficient viral library delivery and library quantification, irregular viral transduction rates, and injection site inflammation disrupting gene expression. Here, we introduce a systemic MPRA (sysMPRA) to address these challenges through systemic intravenous AAV viral delivery. Results: We demonstrate successful transduction of the MPRA library into diverse mouse tissues, efficiently identifying tissue specificity in candidate enhancers and aligning well with predictions from machine learning models. We highlight that sysMPRA effectively uncovers regulatory effects stemming from the disruption of MEF2C transcription factor binding sites, single-nucleotide polymorphisms, and the consequences of genetic variations associated with late-onset Alzheimer's disease. Conclusion: SysMPRA is an effective library delivering method that simultaneously determines the transcriptional functions of hundreds of enhancers in vivo across multiple tissues.

The Chromosome 19 miRNA Cluster Guards Trophoblasts Against Overacting Innate Immunity

To maintain pregnancy health, the human placenta delicately balances protection of the developing fetus from invading pathogens with suppression of excessive inflammation that could lead to fetal and neonatal autoimmune disorders. Previous research, including our own, has shown that small RNA products of the Chromosome 19 MicroRNA Cluster (C19MC) promote viral resistance in non-trophoblastic cells. However, the role of C19MC products in placental trophoblasts remained unclear. Here, we analyzed chromatin accessibility in the C19MC enhancer and identified a previously unknown regulatory domain. Deletion of this domain silenced the expression of C19MC microRNA and Alu elements in trophoblasts. This silencing unexpectedly led to marked activation of cellular innate immune response and strikingly increased Toll-like receptor 3 (TLR3)-mediated sensitivity to poly(I:C), a viral RNA mimic. Our data suggest that C19MC non-coding RNAs interfere with endosomal TLR3 activation in trophoblasts, highlighting a previously unrecognized mechanism for hindrance of excessive innate immune activation.

NAD + Sensing by PARP7 Regulates the C/EBPβ-Dependent Transcription Program in Adipose Tissue In Vivo

We have identified PARP7, an NAD + -dependent mono(ADP-ribosyl) transferase, as a key regulator of the C/EBPβ-dependent proadipogenic transcription program. Moreover, PARP7 is required for efficient adipogenesis and downstream biological functions, including involution of the lactating mammary gland. PARP7 serves as a coregulator of C/EBPβ, and depletion of PARP7 causes a dramatic reduction in C/EBPβ binding across the genome. PARP7 functions as a sensor of nuclear NAD + levels to control gene expression. At the relatively high nuclear NAD + concentrations in undifferentiated preadipocytes, PARP7 is catalytically active for auto- mono(ADP-ribosyl)ation (autoMARylation). As nuclear NAD + concentrations decline post- differentiation, autoMARylation decreases dramatically. AutoMARylation promotes instability of PARP7 through an E3 ligase-ubiquitin-proteasome pathway mediated by the ADP-ribose (ADPR)-binding ubiquitin E3 ligases DTX2 and RNF114. Genetic depletion of PARP7 in mice promotes a dramatic reduction in a wide array of lipids in the mammary gland fat pads and milk from lactating females, as well as a significant decrease in nicotinamide mononucleotide (NMN), a key nutrient in mother's milk. The latter is due to reduced expression of Nampt , the gene encoding NAMPT, the enzyme that produces NMN, which is a direct transcriptional target of PARP7 and C/EBPβ. Collectively, our results extend the biology of PARP7 to adipogenesis and perinatal health. Moreover, our results describe the molecular events that regulate these downstream biological functions.

Active regulatory elements recruit cohesin to establish cell specific chromatin domains

As the 3D structure of the genome is analysed at ever increasing resolution it is clear that there is considerable variation in the 3D chromatin architecture across different cell types. It has been proposed that this may, in part, be due to increased recruitment of cohesin to activated cis-elements (enhancers and promoters) leading to cell-type specific loop extrusion underlying the formation of new sub-TADs. Here we show that cohesin correlates well with the presence of active enhancers and that this varies in an allele-specific manner with the presence or absence of polymorphic enhancers which vary from one individual to another. Using the alpha globin cluster as a model, we show that when all enhancers are removed, peaks of cohesin disappear from these regions and the erythroid specific sub-TAD is no longer formed. Re-insertion of the major alpha globin enhancer (R2) is associated with re-establishment of recruitment and increased interactions. In complementary experiments insertion of the R2 enhancer element into a "neutral" region of the genome recruits cohesin, induces transcription and creates a new large (75 kb) erythroid-specific domain. Together these findings support the proposal that active enhancers recruit cohesin, stimulate loop extrusion and promote the formation of cell specific sub-TADs.

Chromatin unfolding via loops can drive clustered transposon insertion

Transposons, DNA sequences capable of relocating within the genome, make up a significant portion of eukaryotic genomes and are often found in clusters. Within the cell nucleus, the genome is organized into chromatin, a structure with varying degrees of compaction due to three-dimensional folding. Transposon insertion or activation can lead to chromatin decompaction, increasing accessibility and potentially facilitating further nearby insertions. This positive feedback between chromatin unfolding and transposon insertion may result in transposon clustering. Here, we combine bioinformatics with polymer modeling to explore possible mechanisms and conditions that promote clustered transposon insertions. Our analysis of human cell line genomic repeat data reveals extensive clustering of heterochromatic LINE-1 elements and euchromatic Alu elements. For Alu elements, this clustering correlates with increased chromatin accessibility. Both Alu and LINE-1 deviate in their sequence-inherent flexibility from the overall genome, with above-average flexibility for Alu and below-average flexibility for most LINE-1 sequences. Flexibility was highest in young transposons, so that young Alu and LINE-1 exceed overall genome flexibility. We developed an according polymer model of transposon insertion, consisting of a self-attracting chromatin domain. Transposon insertions locally disrupt self-attraction, leading to unfolding of the domain as more transposons are inserted. In simulations where transposons are inserted adjacent to existing ones, we observed gradual unfolding through loop extensions from a folded core. Including transposases as explicit particles, our model shows that adjacent transposon insertion occurs when densely packed chromatin excludes transposases or when insertion rates exceed the thermal equilibration rate of polymer configurations. We conclude that 1) dense chromatin packing that hinders transposase access as well as 2) a local loss of compaction upon transposon insertion favor clustered transposon insertion via loop formation. This biophysical mechanism of clustered insertion site preference would act in combination with selective pressures shaping transposon distribution over evolutionary timescales.

Lactate activates trained immunity by fueling the tricarboxylic acid cycle and regulating histone lactylation

Trained immunity refers to the long-term memory of the innate immune cells. However, little is known about how environmental nutrient availability influences trained immunity. This study finds that physiologic carbon sources impact glucose contribution to the tricarboxylic acid (TCA) cycle and enhance cytokine production of trained monocytes. Our experiments demonstrate that trained monocytes preferentially employe lactate over glucose as a TCA cycle substrate, and lactate metabolism is required for trained immune cell responses to bacterial and fungal infection. Except for the contribution to the TCA cycle, endogenous lactate or exogenous lactate also supports trained immunity by regulating histone lactylation. Further transcriptome analysis, ATAC-seq, and CUT&Tag-seq demonstrate that lactate enhance chromatin accessibility in a manner dependent histone lactylation. Inhibiting lactate-dependent metabolism by silencing lactate dehydrogenase A (LDHA) impairs both lactate fueled the TCA cycle and histone lactylation. These findings suggest that lactate is the hub of immunometabolic and epigenetic programs in trained immunity.

Three-dimensional regulatory hubs support oncogenic programs in glioblastoma

Dysregulation of enhancer-promoter communication in the three-dimensional (3D) nucleus is increasingly recognized as a potential driver of oncogenic programs. Here, we profiled the 3D enhancer-promoter networks of patient-derived glioblastoma stem cells to identify central regulatory nodes. We focused on hyperconnected 3D hubs and demonstrated that hub-interacting genes exhibit high and coordinated expression at the single-cell level and are associated with oncogenic programs that distinguish glioblastoma from low-grade glioma. Epigenetic silencing of a recurrent hub-with an uncharacterized role in glioblastoma-was sufficient to cause downregulation of hub-connected genes, shifts in transcriptional states, and reduced clonogenicity. Integration of datasets across 16 cancers identified "universal" and cancer-type-specific 3D hubs that enrich for oncogenic programs and factors associated with worse prognosis. Genetic alterations could explain only a small fraction of hub hyperconnectivity and increased activity. Overall, our study provides strong support for the potential central role of 3D regulatory hubs in controlling oncogenic programs and properties.

A non-canonical role for the tyrosyl tRNA synthetase: YARS regulates senescence induction and escape and controls the transcription of LIN9

Senescence is a tumor suppressor mechanism triggered by oncogene expression and chemotherapy treatment. It orchestrates a definitive cessation of cell proliferation through the activation of the p53-p21 and p16-Rb pathways, coupled with the compaction of proliferative genes within heterochromatin regions. Some cancer cells have the ability to elude this proliferative arrest but the signaling pathways involved in circumventing senescence remain to be characterized. We have recently described that malignant cells capable of evading senescence have an increased expression of specific tRNAs, such as tRNA-Leu-CAA and tRNA-Tyr-GTA, alongside the activation of their corresponding tRNA ligases, namely LARS and YARS. We have previously shown that YARS promotes senescence escape by activating proliferation and cell cycle genes but its functions during this proliferative arrest remain largely unknown. In this study, we have continued to characterize the functions of YARS, describing non-canonical transcriptional functions of the ligase. Our results show that YARS is present in the nucleus of proliferating and senescent cells and interacts with the Trim28 transcriptional regulator. Importantly, YARS binds to the LIN9 promoter, a critical member of the Dream complex responsible for regulating cell cycle gene transcription. The ligase facilitates the binding and the phosphorylation of the type II RNA polymerase and promotes the deposition of activating epigenetic marks on the LIN9 promoter. Consequently, during senescence escape, YARS activates LIN9 expression and both proteins are necessary to induce the proliferation of emergent cells. These results underscore unconventional transcriptional functions of YARS in activating LIN9 expression in proliferating cells and during senescence escape.

Assessing the virological response to direct-acting antiviral therapies in the HBV cure programme

Hepatitis B virus (HBV) is a global health problem with over 250 million people affected worldwide. Nucleos(t)ide analogues remain the standard of care and suppress production of progeny virions; however, they have limited effect on the viral transcriptome and long-term treatment is associated with off-target toxicities. Promising results are emerging from clinical trials and several drug classes have been evaluated, including capsid assembly modulators and RNA interfering agents. Whilst peripheral biomarkers are used to monitor responses and define treatment endpoints, they fail to reflect the full reservoir of infected hepatocytes. Given these limitations, consideration should be given to the merits of sampling liver tissue, especially in the context of clinical trials. In this review article, we will discuss methods for profiling HBV in liver tissue and their value to the HBV cure programme.

Towards multimodal foundation models in molecular cell biology

The rapid advent of high-throughput omics technologies has created an exponential growth in biological data, often outpacing our ability to derive molecular insights. Large-language models have shown a way out of this data deluge in natural language processing by integrating massive datasets into a joint model with manifold downstream use cases. Here we envision developing multimodal foundation models, pretrained on diverse omics datasets, including genomics, transcriptomics, epigenomics, proteomics, metabolomics and spatial profiling. These models are expected to exhibit unprecedented potential for characterizing the molecular states of cells across a broad continuum, thereby facilitating the creation of holistic maps of cells, genes and tissues. Context-specific transfer learning of the foundation models can empower diverse applications from novel cell-type recognition, biomarker discovery and gene regulation inference, to in silico perturbations. This new paradigm could launch an era of artificial intelligence-empowered analyses, one that promises to unravel the intricate complexities of molecular cell biology, to support experimental design and, more broadly, to profoundly extend our understanding of life sciences.

Reframing Formalin: A Molecular Opportunity Enabling Historical Epigenomics and Retrospective Gene Expression Studies

Formalin preservation of museum specimens has long been considered a barrier to molecular research due to extensive crosslinking and chemical modification. However, recent optimisation of hot alkaline lysis and proteinase K digestion DNA extraction methods have enabled a growing number of studies to overcome these challenges and conduct genome-wide re-sequencing and targeted locus-specific sequencing. The newest, and perhaps most unexpected utility of formalin preservation in archival samples is its ability to preserve in situ DNA-protein interactions at a molecular level. Retrieving this signal provides information about the relative compaction or accessibility of the genome to the transcriptional machinery required for gene expression. Thus, exposure to formalin essentially corresponds to taking a snapshot of organism-wide gene expression at the time of death. While DNA methylation and RNA-Seq analyses of dried tissues have provided glimpses into historical gene regulation, these techniques were previously limited to skeletal or desiccated remains, offering only partial insights. By examining fluid-preserved specimens, molecular tools can now be applied to a broader range of tissues, enabling more detailed tissue-specific gene regulation profiling across vertebrates. In this review, we chronicle the historical use of formaldehyde in collections and discuss how targeted chromatin profiling with assays like MNase-seq and FAIRE-seq are surmounting fixation challenges and unlocking invaluable insights into historical genomes and gene expression profiles. The deeper integration of molecular genetics with museum collections bridges the gap between past and present and provides a vital tool that could help us predict and mitigate some of the impacts of future environmental change, novel pathogens, or invasive species.

Crypt density and recruited enhancers underlie intestinal tumour initiation

Oncogenic mutations that drive colorectal cancer can be present in healthy intestines for long periods without overt consequence1,2. Mutation of Apc, the most common initiating event in conventional adenomas3, activates Wnt signalling, thus conferring fitness on mutant intestinal stem cells (ISCs)4,5. Apc mutations may occur in ISCs that arise by routine self-renewal or by dedifferentiation of their progeny. Although ISCs of these different origins are fundamentally similar6,7, it is unclear whether both generate tumours equally well in uninjured intestines. It is also unknown whether cis-regulatory elements are substantively modulated upon Wnt hyperactivation or as a feature of subsequent tumours. Here we show in two mouse models that adenomas are not an obligatory outcome of Apc deletion in either ISC source, but require proximity of mutant intestinal crypts. Reduced crypt density abrogates, and aggregation of mutant colonic crypts augments, adenoma formation. Moreover, adenoma-resident ISCs open chromatin at thousands of enhancers that are inaccessible in Apc-null ISCs that are not associated with adenomas. These cis elements explain adenoma-selective gene activity and persist, with little further expansion of the repertoire, as other oncogenic mutations accumulate. Thus, cooperativity between neighbouring mutant crypts and new accessibility at specific enhancers are key steps early in intestinal tumorigenesis.

Scalable high-performance single cell data analysis with BPCells

The growth of single-cell datasets to multi-million cell atlases has uncovered major scalability problems for single-cell analysis software. Here, we present BPCells, a package for high-performance single-cell analysis of RNA-seq and ATAC-seq datasets. BPCells uses disk-backed streaming compute algorithms to reduce memory requirements by nearly 70-fold compared to in-memory workflows with little to no loss of execution speed. BPCells also introduces high-performance compressed formats based on bitpacking compression for ATAC-seq fragment files and single-cell sparse matrices. These novel compression algorithms help to accelerate disk-backed analysis by reducing data transfer from disk, while providing the lowest computational overhead of all compression algorithms tested. Using BPCells, we perform normalization and PCA of a 44 million cell dataset on a laptop, demonstrating that BPCells makes working with the largest contemporary single-cell datasets feasible on modest hardware, while leaving headroom on servers for future datasets an order of magnitude larger.

Time-course transcriptome and chromatin accessibility analyses reveal the dynamic transcriptional regulation shaping spikelet hull size

The grains of rice (Oryza sativa) are enclosed by a spikelet hull comprising the lemma and palea. Development of the spikelet hull determines the storage capacity of the grain, thus affecting grain yield and quality. Although multiple signaling pathways controlling grain size have been identified, the transcriptional regulatory mechanisms underlying grain development remain limited. Here, we used RNA-seq and ATAC-seq to characterize the transcription and chromatin accessibility dynamics during the development of spikelet hulls. A time-course analysis showed that more than half of the genes were sequentially expressed during hull development and that the accessibility of most open chromatin regions (OCRs) changed moderately, although some regions positively or negatively affected the expression of their closest genes. We revealed a crucial role of GROWTH-REGULATING FACTORs in shaping grain size by influencing multiple metabolic and signaling pathways, and a coordinated transcriptional regulation in response to auxin and cytokinin signaling. We also demonstrated the function of SCL6-IIb, a member of the GRAS family transcription factors, in regulating grain size, with SCL6-IIb expression being activated by SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 18 (OsSPL18). When we edited the DNA sequences within OCRs upstream of the start codon of BRASSINAZOLE-RESISTANT 1 (BZR1) and SCL6-IIb, we generated multiple mutant lines with longer grains. These findings offer a comprehensive overview of the cis-regulatory landscape involved in forming grain capacity and a valuable resource for exploring the regulatory network behind grain development.

The epigenetic landscape of fate decisions in T cells

Specialized T cell subsets mediate adaptive immunity in response to cytokine signaling and specific transcription factor activity. The epigenetic landscape of T cells has an important role in facilitating and establishing T cell fate decisions. Here, we review the interplay between transcription factors, histone modifications, DNA methylation and three-dimensional chromatin organization to define key elements of the epigenetic landscape in T cells. We introduce key technologies in the areas of sequencing, microscopy and proteomics that have enabled the multi-scale profiling of the epigenetic landscape. We highlight the dramatic changes of the epigenetic landscape as multipotent progenitor cells commit to the T cell lineage during development and discuss the epigenetic changes that favor the emergence of CD4+ and CD8+ T cells. Finally, we discuss the inheritance of epigenetic marks and its potential effects on immune responses as well as therapeutic strategies with potential for epigenetic regulation.

An organoid library of human esophageal squamous cell carcinomas (ESCCs) uncovers the chemotherapy-resistant ESCC features

Esophageal squamous cell carcinoma (ESCC) is a deadly cancer with a poor prognosis and a high recurrence rate after chemotherapy, posing a significant clinical challenge. To elucidate the molecular basis of chemotherapy (chemo)-resistance and to develop methods to effectively eliminate chemo-resistant tumor clones, we established an ESCC organoid (ESCCO) library from 24 ESCC patients of various stages, ages, and treatments. These ESCCOs faithfully recapitulate the oncogenic mutations observed in the original ESCC tissues and manifest tumorigenic properties when xenografted. The ESCCOs respond differently to cisplatin and 5-fluorouracil, chemotherapeutic agents commonly used to treat ESCC patients, with 7 ESCCOs exhibiting potent chemo-resistance. Notably, the chemo-resistant ESCCOs show higher genes involved in antioxidant stress response pathways and more accessible chromatin at their loci than the sensitive ESCCOs. These genes can serve as valuable biomarkers to stratify chemo-resistant ESCCs in histopathological specimens. Through drug screening using the ESCCO library, we reveal that fedratinib effectively induces cell death in chemo-resistant ESCCOs. Collectively, our human ESCCO model offers novel insights into the mechanism of chemo-resistance in ESCCs, which is critical for developing effective therapeutic approaches to eradicate the recurrence of ESCCs.

Changes in the Immune Profile and Chromatin Accessibility of Peripheral Regulatory T Cells in Psoriasis Patients Before and After Treatment With Biologics

Psoriasis is a chronic inflammatory skin disease. The excessive activation of proinflammatory cytokines interleukin-17 (IL-17), IL-23 and T helper cell 17(Th17) is the main pathogenic factor. In addition, the dysfunction of suppressor cells such as regulatory T cells (Tregs) and the imbalance of the Th17/Treg ratio also play important roles in the pathogenesis of psoriasis. By testing the immune function of peripheral Tregs in psoriasis, psoriasis treated with anti-IL-17 biologics, and healthy controls, we found that the number and function of psoriatic peripheral Tregs were abnormal, and Tregs differentiated from 'inhibitory' to 'inflammatory' cells in the inflammatory environment, which may be the cause of Tregs dysfunction in psoriasis. We also found through the assay for targeting accessible chromatin with high-throughput sequencing (ATAC-seq) analysis that the chromatin accessibility of psoriatic peripheral Tregs was significantly higher than that of healthy controls and decreased after treatment, which may be related to INO80, a gene that controls changes in chromatin tightness or relaxation status. In addition, the differentially expressed genes (DEGs) of three groups, such as NCAM2, CDH18, ZEB1 and CCDC22, were mainly concentrated in the signalling pathways related to effector T(Teff) cell aggregation and Tregs dysfunction. This study provides an important basis for the study of peripheral Tregs dysfunction in psoriasis.

Deciphering T-cell exhaustion in the tumor microenvironment: paving the way for innovative solid tumor therapies

In solid tumors, the tumor microenvironment (TME) is a complex mix of tumor, immune, stromal cells, fibroblasts, and the extracellular matrix. Cytotoxic T lymphocytes (CTLs) constitute a fraction of immune cells that may infiltrate into the TME. The primary function of these T-cells is to detect and eliminate tumor cells. However, due to the immunosuppressive factors present in the TME primarily mediated by Myeloid-Derived Suppressor Cells (MDSCs), Tumor associated macrophages (TAMs), Cancer Associated Fibroblasts (CAFs) as well as the tumor cells themselves, T-cells fail to differentiate into effector cells or become dysfunctional and are unable to eliminate the tumor. In addition, chronic antigen stimulation within the TME also leads to a phenomenon, first identified in chronic lymphocytic choriomeningitis virus (LCMV) infection in mice, where the T-cells become exhausted and lose their effector functions. Exhausted T-cells (Tex) are characterized by the presence of remarkably conserved inhibitory receptors, transcription and signaling factors and the downregulation of key effector molecules. Tex cells have been identified in various malignancies, including melanoma, colorectal and hepatocellular cancers. Recent studies have indicated novel strategies to reverse T-cell exhaustion. These include checkpoint inhibitor blockade targeting programmed cell death protein 1 (PD-1), T-cell immunoglobulin and mucin-domain containing-3 (Tim-3), cytotoxic T-lymphocyte associated protein 4 (CTLA-4), or combinations of different immune checkpoint therapies (ICTs) or combination of ICTs with cytokine co-stimulation. In this review, we discuss aspects of T-cell dysfunction within the TME with a focus on T-cell exhaustion. We believe that gaining insight into the mechanisms of T-cell exhaustion within the TME of human solid tumors will pave the way for developing therapeutic strategies to target and potentially re-invigorate exhausted T-cells in cancer.

Intestinal secretory differentiation reflects niche-driven phenotypic and epigenetic plasticity of a common signal-responsive terminal cell

Enterocytes and four classic secretory cell types derive from intestinal epithelial stem cells. Based on morphology, location, and canonical markers, goblet and Paneth cells are considered distinct secretory types. Here, we report high overlap in their transcripts and sites of accessible chromatin, in marked contrast to those of their enteroendocrine or tuft cell siblings. Mouse and human goblet and Paneth cells express extraordinary fractions of few antimicrobial genes, which reflect specific responses to local niches. Wnt signaling retains some ATOH1+ secretory cells in crypt bottoms, where the absence of BMP signaling potently induces Paneth features. Cells that migrate away from crypt bottoms encounter BMPs and thereby acquire goblet properties. These phenotypes and underlying accessible cis-elements interconvert in post-mitotic cells. Thus, goblet and Paneth properties represent alternative phenotypic manifestations of a common signal-responsive terminal cell type. These findings reveal exquisite niche-dependent cell plasticity and cis-regulatory dynamics in likely response to antimicrobial needs.

Dynamic single cell transcriptomics defines kidney FGF23/KL bioactivity and novel segment-specific inflammatory targets

Fibroblast growth factor 23 (FGF23) via its coreceptor αKlotho (KL) provides critical control of phosphate metabolism, which is altered in both rare and very common syndromes. However, the spatial-temporal mechanisms dictating kidney FGF23 functions remain poorly understood. Thus, developing approaches to modify specific FGF23-dictated pathways has proven problematic. Herein, wild type mice were injected with rFGF23 for one, four and 12h and kidney FGF23 bioactivity was determined at single cell resolution. Computational analysis identified distinct epithelial, endothelial, stromal, and immune cell clusters, with differential expressional analysis uniquely tracking FGF23 bioactivity at each time point. FGF23 actions were sex independent but critically relied upon constitutive KL expression mapped within proximal tubule (segments S1-S3) and distal convoluted tub/connecting tubule cell sub-populations. Temporal KL-dependent FGF23 responses drove unique and transient cellular identities, including genes in key MAPK-signaling and vitamin D-metabolic pathways via early- (transcription factor AP-1-related) and late-phase (initiation factor EIF2 signaling) transcriptional regulons. Combining ATACseq/RNAseq data from a cell line stably expressing KL with the in vivo scRNAseq pinpointed genomic accessibility changes in MAPK-dependent genes, including the identification of FGF23-dependent early growth factor-1 distal enhancers. Finally, we identified unexpected crosstalk between FGF23-mediated MAPK signaling and pro inflammatory TNF receptor activation via transcription factor NF-κB, which blocked FGF23 bioactivity in vitro and in vivo. Collectively, our findings have uncovered novel pathways at the single cell level that likely influence FGF23-dependent disease mechanisms.