Search

Abstract Background Domains are the building blocks of proteins. During evolution, they have been duplicated, fused and recombined, to produce proteins with novel structures and functions. Structural and genome-scale studies have shown that pairs or groups of domains observed together in a protein are almost always found in only one N to C terminal order and are the result of a single recombination event that has been propagated by duplication of the multi-domain unit. Previous studies of domain organisation have used graph theory to represent the co-occurrence of domains within proteins. We build on this approach by adding directionality to the graphs and connecting nodes based on their relative order in the protein. Most of the time, the linear order of domains is conserved. However, using the directed graph representation we have identified non-linear features of domain organization that are over-represented in genomes. Recognising these patterns and unravelling how they have arisen may allow us to understand the functional relationships between domains and understand how the protein repertoire has evolved. Results We identify groups of domains that are not linearly conserved, but instead have been shuffled during evolution so that they occur in multiple different orders. We consider 192 genomes across all three kingdoms of life and use domain and protein annotation to understand their functional significance. To identify these features and assess their statistical significance, we represent the linear order of domains in proteins as a directed graph and apply graph theoretical methods. We describe two higher-order patterns of domain organisation: clusters and bi-directionally associated domain pairs and explore their functional importance and phylogenetic conservation. Conclusion Taking into account the order of domains, we have derived a novel picture of global protein organization. We found that all genomes have a higher than expected degree of clustering and more domain pairs in forward and reverse orientation in different proteins relative to random graphs with identical degree distributions. While these features were statistically over-represented, they are still fairly rare. Looking in detail at the proteins involved, we found strong functional relationships within each cluster. In addition, the domains tended to be involved in protein-protein interaction and are able to function as independent structural units. A particularly striking example was the human Jak-STAT signalling pathway which makes use of a set of domains in a range of orders and orientations to provide nuanced signaling functionality. This illustrated the importance of functional and structural constraints (or lack thereof) on domain organisation.

Abstract Background Different regions in a genome evolve at different rates depending on structural and functional constraints. Some genomic regions are highly conserved during metazoan evolution, while other regions may evolve rapidly, either in all species or in a lineage-specific manner. A strong or even moderate change in constraints in functional regions, for example in coding regions, can have significant evolutionary consequences. Results Here we discuss a novel framework, 'BaseDiver', to classify groups of genes in humans based on the patterns of evolutionary constraints on polymorphic positions in their coding regions. Comparing the nucleotide-level divergence among mammals with the extent of deviation from the ancestral base in the human lineage, we identify patterns of evolutionary pressure on nonsynonymous base-positions in groups of genes belonging to the same functional category. Focussing on groups of genes in functional categories, we find that transcription factors contain a significant excess of nonsynonymous base-positions that are conserved in other mammals but changed in human, while immunity related genes harbour mutations at base-positions that evolve rapidly in all mammals including humans due to strong preference for advantageous alleles. Genes involved in olfaction also evolve rapidly in all mammals, and in humans this appears to be due to weak negative selection. Conclusion While recent studies have identified genes under positive selection in humans, our approach identifies evolutionary constraints on Gene Ontology groups identifying changes in humans relative to some of the other mammals.

The Human Cell Atlas (HCA) is a global partnership "to create comprehensive reference maps of all human cells-the fundamental units of life - as a basis for both understanding human health and diagnosing, monitoring, and treating disease." ( https://www.humancellatlas.org/ ) The atlas shall characterize cells from diverse individuals across the globe to better understand human biology. HCA proactively considers the priorities of, and benefits accrued to, contributing communities. Here, we lay out principles and action items that have been adopted to affirm HCA's commitment to equity so that the atlas is beneficial to all of humanity., Competing Interests: Competing interests: A.R. is an employee of Genentech, a member of the Roche group, and has equity in Roche. A.R. was a co-founder and equity holder of Celsius Therapeutics, is an equity holder in Immunitas, and until July 31, 2020 was an S.A.B. member of Thermo Fisher Scientific, Syros Pharmaceuticals, Neogene Therapeutics and Asimov. A.R. is an inventor on multiple patents related to single cell and spatial genomics. All other authors declare no competing interests., (© 2024. The Author(s).)

Positional coding along the anterior-posterior axis is regulated by HOX genes, whose 3' to 5' expression correlates with location along this axis. The precise utilisation of HOX genes in different human cell types is not fully understood. Here, we use single-cell and spatial-transcriptomics, along with in-situ sequencing, to create a developmental atlas of the human fetal spine. We analyse HOX gene expression across cell types during development, finding that neural-crest derivatives unexpectedly retain the anatomical HOX code of their origin while also adopting the code of their destination. This trend is confirmed across multiple organs. In the axial plane of the spinal cord, we find distinct patterns in the ventral and dorsal domains, providing insights into motor pool organisation and loss of collinearity in HOXB genes. Our findings shed new light on HOX gene expression in the developing spine, highlighting a HOX gene 'source code' in neural-crest cell derivatives., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

With the convergence in exciting advances in molecular and spatial profiling methods and new computational approaches leveraging artificial intelligence and machine learning (AI/ML), the construction of cell atlases is progressing from data collection to atlas integration and beyond. Here, we explore five ways in which cell atlases, including the Human Cell Atlas, are already revealing valuable biological insights, and how they are poised to provide even greater benefits in the coming years. In particular, we discuss cell atlases as censuses of cells; as 3D maps of cells in the body, across modalities and scales; as maps connecting genotype causes to phenotype effects; as 4D maps of development; and, ultimately, as foundation models of biology unifying all these aspects and helping to transform medicine., (© 2024. Springer Nature Limited.)

Recent technological advancements in single-cell genomics have enabled joint profiling of gene expression and alternative modalities at unprecedented scale. Consequently, the complexity of multi-omics data sets is increasing massively. Existing models for multi-modal data are typically limited in functionality or scalability, making data integration and downstream analysis cumbersome. We present multiDGD, a scalable deep generative model providing a probabilistic framework to learn shared representations of transcriptome and chromatin accessibility. It shows outstanding performance on data reconstruction without feature selection. We demonstrate on several data sets from human and mouse that multiDGD learns well-clustered joint representations. We further find that probabilistic modeling of sample covariates enables post-hoc data integration without the need for fine-tuning. Additionally, we show that multiDGD can detect statistical associations between genes and regulatory regions conditioned on the learned representations. multiDGD is available as an scverse-compatible package on GitHub., Competing Interests: Competing interests: S.A.T. is a Scientific Advisory Board member of ForeSite Labs, Qiagen, OMass, and is a co-founder and equity holder of TransitionBio and EnsoCell Therapeutics, a non-executive director of 10x Genomics, as well as a part-time employee of GlaxoSmithKline. A.K., E.D., and V.S. declare no competing interests., (© 2024. The Author(s).)

The human vascular system, comprising endothelial (EC) and mural cells, covers a vast surface area in the body, providing a critical interface between blood and tissue environments. Functional differences exist across specific vascular beds, but their molecular determinants across tissues remain largely unknown. Here, we integrated single-cell transcriptomics data from 19 human organs and tissues, and defined 42 vascular cell states from ~67,000 cells (62 donors), including angiotypic transitional signatures along the arterial endothelial axis from large to small calibre vessels. We also characterised organotypic populations, including splenic littoral ECs and blood-brain barrier cells, thus clarifying the molecular profiles of these important cell states. Interrogating endothelial-mural cell molecular crosstalk revealed angiotypic and organotypic communication pathways related to Notch, Wnt, retinoic acid, prostaglandin, and cell adhesion signalling. Transcription factor network analysis revealed differential regulation of downstream target genes in tissue-specific modules, such as FOXF1 target genes across multiple lung vascular subpopulations. Additionally, we make mechanistic inferences of vascular drug targets within different vascular beds. This open access resource enhances our understanding of angiodiversity and organotypic molecular signatures in human vascular cells and has therapeutic implications for vascular diseases across tissues., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Cutaneous T cell lymphoma (CTCL) is a potentially fatal clonal malignancy of T cells primarily affecting the skin. The most common form of CTCL, mycosis fungoides, can be difficult to diagnose, resulting in treatment delay. We performed single-cell and spatial transcriptomics analysis of skin from patients with mycosis fungoides-type CTCL and an integrated comparative analysis with human skin cell atlas datasets from healthy and inflamed skin. We revealed the co-optation of T helper 2 (T H 2) cell-immune gene programs by malignant CTCL cells and modeling of the tumor microenvironment to support their survival. We identified MHC-II + fibroblasts and dendritic cells that can maintain T H 2 cell-like tumor cells. CTCL tumor cells are spatially associated with B cells, forming tertiary lymphoid structure-like aggregates. Finally, we validated the enrichment of B cells in CTCL and its association with disease progression across three independent patient cohorts. Our findings provide diagnostic aids, potential biomarkers for disease staging and therapeutic strategies for CTCL., Competing Interests: Competing interests In the past 3 years, S.A.T. has received remuneration for scientific advisory board membership from Sanofi, GlaxoSmithKline, Foresite Labs and Qiagen. She is a co-founder of and holds equity in Transition Bio and Ensocell. From 8 January 2024, she has been a part-time employee of GlaxoSmithKline. The remaining authors declare no competing interests., (© 2024. The Author(s).)

T cells develop from circulating precursor cells, which enter the thymus and migrate through specialized subcompartments that support their maturation and selection 1 . In humans, this process starts in early fetal development and is highly active until thymic involution in adolescence. To map the microanatomical underpinnings of this process in pre- and early postnatal stages, we established a quantitative morphological framework for the thymus-the Cortico-Medullary Axis-and used it to perform a spatially resolved analysis. Here, by applying this framework to a curated multimodal single-cell atlas, spatial transcriptomics and high-resolution multiplex imaging data, we demonstrate establishment of the lobular cytokine network, canonical thymocyte trajectories and thymic epithelial cell distributions by the beginning of the the second trimester of fetal development. We pinpoint tissue niches of thymic epithelial cell progenitors and distinct subtypes associated with Hassall's corpuscles and identify divergence in the timing of medullary entry between CD4 and CD8 T cell lineages. These findings provide a basis for a detailed understanding of T lymphocyte development and are complemented with a holistic toolkit for cross-platform imaging data analysis, annotation and OrganAxis construction (TissueTag), which can be applied to any tissue., Competing Interests: Competing interests: J.C.M. has been an employee of Genentech since September 2022. S.A.T. is a scientific advisory board member of ForeSite Labs, OMass Therapeutics, Qiagen, a co-founder and equity holder of TransitionBio and EnsoCell Therapeutics, a non-executive director of 10x Genomics and a part-time employee of GlaxoSmithKline. The other authors declare no competing interests., (© 2024. The Author(s).)

The gastrointestinal tract is a multi-organ system crucial for efficient nutrient uptake and barrier immunity. Advances in genomics and a surge in gastrointestinal diseases 1,2 has fuelled efforts to catalogue cells constituting gastrointestinal tissues in health and disease 3 . Here we present systematic integration of 25 single-cell RNA sequencing datasets spanning the entire healthy gastrointestinal tract in development and in adulthood. We uniformly processed 385 samples from 189 healthy controls using a newly developed automated quality control approach (scAutoQC), leading to a healthy reference atlas with approximately 1.1 million cells and 136 fine-grained cell states. We anchor 12 gastrointestinal disease datasets spanning gastrointestinal cancers, coeliac disease, ulcerative colitis and Crohn's disease to this reference. Utilizing this 1.6 million cell resource (gutcellatlas.org), we discover epithelial cell metaplasia originating from stem cells in intestinal inflammatory diseases with transcriptional similarity to cells found in pyloric and Brunner's glands. Although previously linked to mucosal healing 4 , we now implicate pyloric gland metaplastic cells in inflammation through recruitment of immune cells including T cells and neutrophils. Overall, we describe inflammation-induced changes in stem cells that alter mucosal tissue architecture and promote further inflammation, a concept applicable to other tissues and diseases., Competing Interests: Competing interests: S.A.T. is a scientific advisory board member of ForeSite Labs, OMass Therapeutics, a co-founder and equity holder of TransitionBio and EnsoCell Therapeutics, a non-executive director of 10x Genomics and a part-time employee of GlaxoSmithKline. R.E. is an equity holder in EnsoCell. P.K. has consulted for AstraZeneca, UCB, Biomunex and Infinitopes. N.M.P reports consulting fees from Infinitopes. J.S.-R. reports funding from GSK, Pfizer and Sanofi and fees/honoraria from Travere Therapeutics, Stadapharm, Astex, Owkin, Pfizer, Moderna and Grunenthal. A.S. is the recipient of research grants from Roche-Genentech, Abbvie, GSK, Scipher Medicine, Pfizer, Alimentiv, Boehringer Ingelheim and Agomab and has received consulting fees from Genentech, GSK, Pfizer, HotSpot Therapeutics, Alimentiv, Agomab, Goodgut and Orikine. R.E. and S.A.T are inventors on the patent GB2412853.0 filed in the UK, some components of which are related to this work. All other authors declare no competing interests., (© 2024. The Author(s).)

Human embryonic bone and joint formation is determined by coordinated differentiation of progenitors in the nascent skeleton. The cell states, epigenetic processes and key regulatory factors that underlie lineage commitment of these cells remain elusive. Here we applied paired transcriptional and epigenetic profiling of approximately 336,000 nucleus droplets and spatial transcriptomics to establish a multi-omic atlas of human embryonic joint and cranium development between 5 and 11 weeks after conception. Using combined modelling of transcriptional and epigenetic data, we characterized regionally distinct limb and cranial osteoprogenitor trajectories across the embryonic skeleton and further described regulatory networks that govern intramembranous and endochondral ossification. Spatial localization of cell clusters in our in situ sequencing data using a new tool, ISS-Patcher, revealed mechanisms of progenitor zonation during bone and joint formation. Through trajectory analysis, we predicted potential non-canonical cellular origins for human chondrocytes from Schwann cells. We also introduce SNP2Cell, a tool to link cell-type-specific regulatory networks to polygenic traits such as osteoarthritis. Using osteolineage trajectories characterized here, we simulated in silico perturbations of genes that cause monogenic craniosynostosis and implicate potential cell states and disease mechanisms. This work forms a detailed and dynamic regulatory atlas of bone and cartilage maturation and advances our fundamental understanding of cell-fate determination in human skeletal development., Competing Interests: Competing interests: C.D.B. is a founder of Mestag Therapeutics. S.A.T. is a scientific advisory board member of ForeSite Labs, OMass Therapeutics, Qiagen, a co-founder and equity holder of TransitionBio and EnsoCell Therapeutics, a non-executive director of 10x Genomics and a part-time employee of GlaxoSmithKline. Sarah also holds equity in TransitionBio and options for equity in Element Biosciences, where she is a former SAB member.  All other authors declare no competing interests., (© 2024. The Author(s).)

Human prenatal skin is populated by innate immune cells, including macrophages, but whether they act solely in immunity or have additional functions in morphogenesis is unclear. Here we assembled a comprehensive multi-omics reference atlas of prenatal human skin (7-17 post-conception weeks), combining single-cell and spatial transcriptomics data, to characterize the microanatomical tissue niches of the skin. This atlas revealed that crosstalk between non-immune and immune cells underpins the formation of hair follicles, is implicated in scarless wound healing and is crucial for skin angiogenesis. We systematically compared a hair-bearing skin organoid (SkO) model derived from human embryonic stem cells and induced pluripotent stem cells to prenatal and adult skin 1 . The SkO model closely recapitulated in vivo skin epidermal and dermal cell types during hair follicle development and expression of genes implicated in the pathogenesis of genetic hair and skin disorders. However, the SkO model lacked immune cells and had markedly reduced endothelial cell heterogeneity and quantity. Our in vivo prenatal skin cell atlas indicated that macrophages and macrophage-derived growth factors have a role in driving endothelial development. Indeed, vascular network remodelling was enhanced following transfer of autologous macrophages derived from induced pluripotent stem cells into SkO cultures. Innate immune cells are therefore key players in skin morphogenesis beyond their conventional role in immunity, a function they achieve through crosstalk with non-immune cells., Competing Interests: Competing interests: J.L. and K.R.K., with the Indiana University Research and Technology Corporation, have a patent relating to the methodology and composition of SkOs (PCT/US2016/058174). K.R.K. is a consultant for StemCell Technologies. All other authors declare no competing interests., (© 2024. The Author(s).)

Human limbs emerge during the fourth post-conception week as mesenchymal buds, which develop into fully formed limbs over the subsequent months 1 . This process is orchestrated by numerous temporally and spatially restricted gene expression programmes, making congenital alterations in phenotype common 2 . Decades of work with model organisms have defined the fundamental mechanisms underlying vertebrate limb development, but an in-depth characterization of this process in humans has yet to be performed. Here we detail human embryonic limb development across space and time using single-cell and spatial transcriptomics. We demonstrate extensive diversification of cells from a few multipotent progenitors to myriad differentiated cell states, including several novel cell populations. We uncover two waves of human muscle development, each characterized by different cell states regulated by separate gene expression programmes, and identify musculin (MSC) as a key transcriptional repressor maintaining muscle stem cell identity. Through assembly of multiple anatomically continuous spatial transcriptomic samples using VisiumStitcher, we map cells across a sagittal section of a whole fetal hindlimb. We reveal a clear anatomical segregation between genes linked to brachydactyly and polysyndactyly, and uncover transcriptionally and spatially distinct populations of the mesenchyme in the autopod. Finally, we perform single-cell RNA sequencing on mouse embryonic limbs to facilitate cross-species developmental comparison, finding substantial homology between the two species., Competing Interests: Competing interests: In the past 3 years, S.A.T. has consulted for or been a member of scientific advisory boards at Qiagen, Sanofi, GlaxoSmithKline and ForeSite Labs. She is a consultant and equity holder for TransitionBio and EnsoCell. J.C.M has been an employee of Genentech, Inc. since September 2022. The remaining authors declare no competing interests., (© 2023. The Author(s).)

In this issue of Cancer Cell, Gao et al. map fibroblast diversity across tumors and chronic inflammatory tissues. The authors uncover universal fibroblast subtypes such as LRRC15 + and MMP1 + myofibroblasts along with specialized tissue-specific subtypes. They reveal cellular roles of fibroblasts in immunosuppression through stromal niches and cell-cell interactions., Competing Interests: Declaration of interests S.A.T. is a scientific advisory board member of ForeSite Labs, OMass Therapeutics, and Qiagen; a co-founder and equity holder of TransitionBio and EnsoCell Therapeutics; a non-executive director of 10x Genomics; and a part-time employee of GlaxoSmithKline., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Over the last decade, single-cell genomics has revealed remarkable heterogeneity and plasticity of cell types in the lungs and airways. The challenge now is to understand how these cell types interact in three-dimensional space to perform lung functions, facilitating airflow and gas exchange while simultaneously providing barrier function to avoid infection. An explosion in novel spatially resolved gene expression technologies, coupled with computational tools that harness machine learning and deep learning, now promise to address this challenge. Here, we review the most commonly used spatial analysis workflows, highlighting their advantages and limitations, and outline recent developments in machine learning and artificial intelligence that will augment how we interpret spatial data. Together these technologies have the potential to transform our understanding of the respiratory system in health and disease, and we showcase studies in lung development, COVID-19, lung cancer, and fibrosis where spatially resolved transcriptomics is already providing novel insights.

Embryogenesis is a vulnerable time. Mutations in developmental cells can result in the wide dissemination of cells predisposed to disease within mature organs. We characterised the evolutionary history of four synchronous renal tumours from a 14-year-old girl using whole genome sequencing alongside single cell and bulk transcriptomic sequencing. Phylogenetic reconstruction timed the origin of all tumours to a multipotent embryonic cell committed to the right kidney, around 4 weeks post-conception. Biochemical and structural analysis of their shared MTOR mutation, absent from normal tissues, demonstrates enhanced protein flexibility, enabling a FAT domain hinge to dramatically increase activity of mTORC1 and mTORC2. Developmental mutations, not usually detected in traditional genetic screening, have vital clinical importance in guiding prognosis, targeted treatment, and family screening decisions for paediatric tumours., (© 2024. The Author(s).)

Single-cell data analysis can infer dynamic changes in cell populations, for example across time, space or in response to perturbation, thus deriving pseudotime trajectories. Current approaches comparing trajectories often use dynamic programming but are limited by assumptions such as the existence of a definitive match. Here we describe Genes2Genes, a Bayesian information-theoretic dynamic programming framework for aligning single-cell trajectories. It is able to capture sequential matches and mismatches of individual genes between a reference and query trajectory, highlighting distinct clusters of alignment patterns. Across both real world and simulated datasets, it accurately inferred alignments and demonstrated its utility in disease cell-state trajectory analysis. In a proof-of-concept application, Genes2Genes revealed that T cells differentiated in vitro match an immature in vivo state while lacking expression of genes associated with TNF signaling. This demonstrates that precise trajectory alignment can pinpoint divergence from the in vivo system, thus guiding the optimization of in vitro culture conditions., (© 2024. The Author(s).)

Summary: Visium HD by 10X Genomics is the first commercially available platform capable of capturing full scale transcriptomic data paired with a reference morphology image from archived FFPE blocks at sub-cellular resolution. However, aggregation of capture regions to single cells poses challenges. Bin2cell reconstructs cells from the highest resolution data (2 μm bins) by leveraging morphology image segmentation and gene expression information. It is compatible with established Python single cell and spatial transcriptomics software, and operates efficiently in a matter of minutes without requiring a GPU. We demonstrate improvements in downstream analysis when using the reconstructed cells over default 8 μm bins on mouse brain and human colorectal cancer data., Availability and Implementation: Bin2cell is available at https://github.com/Teichlab/bin2cell, along with documentation and usage examples, and can be installed from pip. Probe design functionality is available at https://github.com/Teichlab/gene2probe., (© The Author(s) 2024. Published by Oxford University Press.)

Metabolic dysfunction-associated steatotic liver disease (MASLD, formerly known as non-alcoholic fatty liver disease) is a leading cause of chronic liver disease worldwide. MASLD can progress to metabolic dysfunction-associated steatohepatitis (MASH, formerly known as non-alcoholic steatohepatitis) with subsequent liver cirrhosis and hepatocellular carcinoma formation. The advent of current technologies such as single-cell and single-nuclei RNA sequencing have transformed our understanding of the liver in homeostasis and disease. The next frontier is contextualizing this single-cell information in its native spatial orientation. This understanding will markedly accelerate discovery science in hepatology, resulting in a further step-change in our knowledge of liver biology and pathobiology. In this Review, we discuss up-to-date knowledge of MASLD development and progression and how the burgeoning field of spatial genomics is driving exciting new developments in our understanding of human liver disease pathogenesis and therapeutic target identification., (© 2024. Springer Nature Limited.)

Here we propose "molecular connectomics" to link molecular and morphological cell features in three dimensions across scales, using machine learning and artificial intelligence to reveal emergent properties of complex biological systems., Competing Interests: In the past 3 years, S.A.T. has received remuneration for scientific advisory board membership from Sanofi, GlaxoSmithKline, Foresite Labs and Qiagen. S.A.T. is a co-founder and holds equity in Transition Bio and Ensocell. From 8 January 2024, S.A.T. is a part-time employee of GlaxoSmithKline., (Copyright: © 2024 Megas et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The Human BioMolecular Atlas Program (HuBMAP) aims to construct a reference 3D structural, cellular, and molecular atlas of the healthy adult human body. The HuBMAP Data Portal (https://portal.hubmapconsortium.org) serves experimental datasets and supports data processing, search, filtering, and visualization. The Human Reference Atlas (HRA) Portal (https://humanatlas.io) provides open access to atlas data, code, procedures, and instructional materials. Experts from more than 20 consortia are collaborating to construct the HRA's Common Coordinate Framework (CCF), knowledge graphs, and tools that describe the multiscale structure of the human body (from organs and tissues down to cells, genes, and biomarkers) and to use the HRA to understand changes that occur at each of these levels with aging, disease, and other perturbations. The 6th release of the HRA v2.0 covers 36 organs with 4,499 unique anatomical structures, 1,195 cell types, and 2,089 biomarkers (e.g., genes, proteins, lipids) linked to ontologies and 2D/3D reference objects. New experimental data can be mapped into the HRA using (1) three cell type annotation tools (e.g., Azimuth) or (2) validated antibody panels (OMAPs), or (3) by registering tissue data spatially. This paper describes the HRA user stories, terminology, data formats, ontology validation, unified analysis workflows, user interfaces, instructional materials, application programming interface (APIs), flexible hybrid cloud infrastructure, and previews atlas usage applications.

Objective: Epigenetic mechanisms, including DNA methylation (DNAm), have been proposed to play a key role in Crohn's disease (CD) pathogenesis. However, the specific cell types and pathways affected as well as their potential impact on disease phenotype and outcome remain unknown. We set out to investigate the role of intestinal epithelial DNAm in CD pathogenesis., Design: We generated 312 intestinal epithelial organoids (IEOs) from mucosal biopsies of 168 patients with CD (n=72), UC (n=23) and healthy controls (n=73). We performed genome-wide molecular profiling including DNAm, bulk as well as single-cell RNA sequencing. Organoids were subjected to gene editing and the functional consequences of DNAm changes evaluated using an organoid-lymphocyte coculture and a nucleotide-binding oligomerisation domain, leucine-rich repeat and CARD domain containing 5 (NLRC5) dextran sulphate sodium (DSS) colitis knock-out mouse model., Results: We identified highly stable, CD-associated loss of DNAm at major histocompatibility complex (MHC) class 1 loci including NLRC5 and cognate gene upregulation. Single-cell RNA sequencing of primary mucosal tissue and IEOs confirmed the role of NLRC5 as transcriptional transactivator in the intestinal epithelium. Increased mucosal MHC-I and NLRC5 expression in adult and paediatric patients with CD was validated in additional cohorts and the functional role of MHC-I highlighted by demonstrating a relative protection from DSS-mediated mucosal inflammation in NLRC5-deficient mice. MHC-I DNAm in IEOs showed a significant correlation with CD disease phenotype and outcomes. Application of machine learning approaches enabled the development of a disease prognostic epigenetic molecular signature., Conclusions: Our study has identified epigenetically regulated intestinal epithelial MHC-I as a novel mechanism in CD pathogenesis., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2024. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

The functional diversity of natural killer (NK) cell repertoires stems from differentiation, homeostatic, receptor-ligand interactions and adaptive-like responses to viral infections. In the present study, we generated a single-cell transcriptional reference map of healthy human blood- and tissue-derived NK cells, with temporal resolution and fate-specific expression of gene-regulatory networks defining NK cell differentiation. Transfer learning facilitated incorporation of tumor-infiltrating NK cell transcriptomes (39 datasets, 7 solid tumors, 427 patients) into the reference map to analyze tumor microenvironment (TME)-induced perturbations. Of the six functionally distinct NK cell states identified, a dysfunctional stressed CD56 bright state susceptible to TME-induced immunosuppression and a cytotoxic TME-resistant effector CD56 dim state were commonly enriched across tumor types, the ratio of which was predictive of patient outcome in malignant melanoma and osteosarcoma. This resource may inform the design of new NK cell therapies and can be extended through transfer learning to interrogate new datasets from experimental perturbations or disease conditions., (© 2024. The Author(s).)