Spatiotemporal analysis of human intestinal development at single-cell resolution.

Development of the human intestine is not well understood. Here, we link single-cell RNA sequencing and spatial transcriptomics to characterize intestinal morphogenesis through time. We identify 101 cell states including epithelial and mesenchymal progenitor populations and programs linked to key morphogenetic milestones. We describe principles of crypt-villus axis formation; neural, vascular, mesenchymal morphogenesis, and immune population of the developing gut. We identify the differentiation hierarchies of developing fibroblast and myofibroblast subtypes and describe diverse functions for these including as vascular niche cells. We pinpoint the origins of Peyer's patches and gut-associated lymphoid tissue (GALT) and describe location-specific immune programs. We use our resource to present an unbiased analysis of morphogen gradients that direct sequential waves of cellular differentiation and define cells and locations linked to rare developmental intestinal disorders. We compile a publicly available online resource, spatio-temporal analysis resource of fetal intestinal development (STAR-FINDer), to facilitate further work.

Structural Remodeling of the Human Colonic Mesenchyme in Inflammatory Bowel Disease.

Intestinal mesenchymal cells play essential roles in epithelial homeostasis, matrix remodeling, immunity, and inflammation. But the extent of heterogeneity within the colonic mesenchyme in these processes remains unknown. Using unbiased single-cell profiling of over 16,500 colonic mesenchymal cells, we reveal four subsets of fibroblasts expressing divergent transcriptional regulators and functional pathways, in addition to pericytes and myofibroblasts. We identified a niche population located in proximity to epithelial crypts expressing SOX6, F3 (CD142), and WNT genes essential for colonic epithelial stem cell function. In colitis, we observed dysregulation of this niche and emergence of an activated mesenchymal population. This subset expressed TNF superfamily member 14 (TNFSF14), fibroblastic reticular cell-associated genes, IL-33, and Lysyl oxidases. Further, it induced factors that impaired epithelial proliferation and maturation and contributed to oxidative stress and disease severity in vivo. Our work defines how the colonic mesenchyme remodels to fuel inflammation and barrier dysfunction in IBD.

Can AlphaFold's breakthrough in protein structure help decode the fundamental principles of adaptive cellular immunity?

T cells are essential immune cells responsible for identifying and eliminating pathogens. Through interactions between their T-cell antigen receptors (TCRs) and antigens presented by major histocompatibility complex molecules (MHCs) or MHC-like molecules, T cells discriminate foreign and self peptides. Determining the fundamental principles that govern these interactions has important implications in numerous medical contexts. However, reconstructing a map between T cells and their antagonist antigens remains an open challenge for the field of immunology, and success of in silico reconstructions of this relationship has remained incremental. In this Perspective, we discuss the role that new state-of-the-art deep-learning models for predicting protein structure may play in resolving some of the unanswered questions the field faces linking TCR and peptide-MHC properties to T-cell specificity. We provide a comprehensive overview of structural databases and the evolution of predictive models, and highlight the breakthrough AlphaFold provided the field.

Colonic epithelial cell diversity in health and inflammatory bowel disease.

The colonic epithelium facilitates host-microorganism interactions to control mucosal immunity, coordinate nutrient recycling and form a mucus barrier. Breakdown of the epithelial barrier underpins inflammatory bowel disease (IBD). However, the specific contributions of each epithelial-cell subtype to this process are unknown. Here we profile single colonic epithelial cells from patients with IBD and unaffected controls. We identify previously unknown cellular subtypes, including gradients of progenitor cells, colonocytes and goblet cells within intestinal crypts. At the top of the crypts, we find a previously unknown absorptive cell, expressing the proton channel OTOP2 and the satiety peptide uroguanylin, that senses pH and is dysregulated in inflammation and cancer. In IBD, we observe a positional remodelling of goblet cells that coincides with downregulation of WFDC2-an antiprotease molecule that we find to be expressed by goblet cells and that inhibits bacterial growth. In vivo, WFDC2 preserves the integrity of tight junctions between epithelial cells and prevents invasion by commensal bacteria and mucosal inflammation. We delineate markers and transcriptional states, identify a colonic epithelial cell and uncover fundamental determinants of barrier breakdown in IBD.

Evaluating performance of existing computational models in predicting CD8+ T cell pathogenic epitopes and cancer neoantigens.

T cell recognition of a cognate peptide-major histocompatibility complex (pMHC) presented on the surface of infected or malignant cells is of the utmost importance for mediating robust and long-term immune responses. Accurate predictions of cognate pMHC targets for T cell receptors would greatly facilitate identification of vaccine targets for both pathogenic diseases and personalized cancer immunotherapies. Predicting immunogenic peptides therefore has been at the center of intensive research for the past decades but has proven challenging. Although numerous models have been proposed, performance of these models has not been systematically evaluated and their success rate in predicting epitopes in the context of human pathology has not been measured and compared. In this study, we evaluated the performance of several publicly available models, in identifying immunogenic CD8+ T cell targets in the context of pathogens and cancers. We found that for predicting immunogenic peptides from an emerging virus such as severe acute respiratory syndrome coronavirus 2, none of the models perform substantially better than random or offer considerable improvement beyond HLA ligand prediction. We also observed suboptimal performance for predicting cancer neoantigens. Through investigation of potential factors associated with ill performance of models, we highlight several data- and model-associated issues. In particular, we observed that cross-HLA variation in the distribution of immunogenic and non-immunogenic peptides in the training data of the models seems to substantially confound the predictions. We additionally compared key parameters associated with immunogenicity between pathogenic peptides and cancer neoantigens and observed evidence for differences in the thresholds of binding affinity and stability, which suggested the need to modulate different features in identifying immunogenic pathogen versus cancer peptides. Overall, we demonstrate that accurate and reliable predictions of immunogenic CD8+ T cell targets remain unsolved; thus, we hope our work will guide users and model developers regarding potential pitfalls and unsettled questions in existing immunogenicity predictors.

HLA-dependent variation in SARS-CoV-2 CD8 + T cell cross-reactivity with human coronaviruses.

The conditions and extent of cross-protective immunity between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and common-cold human coronaviruses (HCoVs) remain open despite several reports of pre-existing T cell immunity to SARS-CoV-2 in individuals without prior exposure. Using a pool of functionally evaluated SARS-CoV-2 peptides, we report a map of 126 immunogenic peptides with high similarity to 285 MHC-presented peptides from at least one HCoV. Employing this map of SARS-CoV-2-non-homologous and homologous immunogenic peptides, we observe several immunogenic peptides with high similarity to human proteins, some of which have been reported to have elevated expression in severe COVID-19 patients. After combining our map with SARS-CoV-2-specific TCR repertoire data from COVID-19 patients and healthy controls, we show that public repertoires for the majority of convalescent patients are dominated by TCRs cognate to non-homologous SARS-CoV-2 peptides. We find that for a subset of patients, >50% of their public SARS-CoV-2-specific repertoires consist of TCRs cognate to homologous SARS-CoV-2-HCoV peptides. Further analysis suggests that this skewed distribution of TCRs cognate to homologous or non-homologous peptides in COVID-19 patients is likely to be HLA-dependent. Finally, we provide 10 SARS-CoV-2 peptides with known cognate TCRs that are conserved across multiple coronaviruses and are predicted to be recognized by a high proportion of the global population. These findings may have important implications for COVID-19 heterogeneity, vaccine-induced immune responses, and robustness of immunity to SARS-CoV-2 and its variants.

SARS-CoV-2 peptides bind to NKG2D and increase NK cell activity.

Immune dysregulation is commonly observed in patients with coronavirus disease 2019 (COVID-19). Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) induces severe lung inflammation and innate immune cell dysregulation. However, the precise interaction between SARS-CoV-2 and the innate immune system is currently unknown. To understand the interaction between SARS-CoV-2 and natural killer (NK) cells, several SARS-CoV-2 S protein peptides capable of binding to the NKG2D receptor were screened by in silico analysis. Among them, two peptides, cov1 and cov2, bound to NK cells and NKG2D receptors. These cov peptides increased NK cytotoxicity toward lung cancer cells, stimulated interferon gamma (IFN-γ) production by NK cells, and likely mediated these responses through the phosphorylation of Vav1, a key downstream-signaling molecule of NKG2D and NK activation genes. The direct interaction between SARS-CoV-2 and NK cells is a novel finding, and modulation of this interaction has potential clinical application as a therapeutic target for COVID-19.

Discovery of a CD10-negative B-progenitor in human fetal life identifies unique ontogeny-related developmental programs.

Human lymphopoiesis is a dynamic lifelong process that starts in utero 6 weeks postconception. Although fetal B-lymphopoiesis remains poorly defined, it is key to understanding leukemia initiation in early life. Here, we provide a comprehensive analysis of the human fetal B-cell developmental hierarchy. We report the presence in fetal tissues of 2 distinct CD19+ B-progenitors, an adult-type CD10+ve ProB-progenitor and a new CD10-ve PreProB-progenitor, and describe their molecular and functional characteristics. PreProB-progenitors and ProB-progenitors appear early in the first trimester in embryonic liver, followed by a sustained second wave of B-progenitor development in fetal bone marrow (BM), where together they form >40% of the total hematopoietic stem cell/progenitor pool. Almost one-third of fetal B-progenitors are CD10-ve PreProB-progenitors, whereas, by contrast, PreProB-progenitors are almost undetectable (0.53% ± 0.24%) in adult BM. Single-cell transcriptomics and functional assays place fetal PreProB-progenitors upstream of ProB-progenitors, identifying them as the first B-lymphoid-restricted progenitor in human fetal life. Although fetal BM PreProB-progenitors and ProB-progenitors both give rise solely to B-lineage cells, they are transcriptionally distinct. As with their fetal counterparts, adult BM PreProB-progenitors give rise only to B-lineage cells in vitro and express the expected B-lineage gene expression program. However, fetal PreProB-progenitors display a distinct, ontogeny-related gene expression pattern that is not seen in adult PreProB-progenitors, and they share transcriptomic signatures with CD10-ve B-progenitor infant acute lymphoblastic leukemia blast cells. These data identify PreProB-progenitors as the earliest B-lymphoid-restricted progenitor in human fetal life and suggest that this fetal-restricted committed B-progenitor might provide a permissive cellular context for prenatal B-progenitor leukemia initiation.

A new phosphate-starvation response in fission yeast requires the endocytic function of myosin I.

Endocytosis is essential for uptake of many substances into the cell, but how it links to nutritional signalling is poorly understood. Here, we show a new role for endocytosis in regulating the response to low phosphate in Schizosaccharomyces pombe. Loss of function of myosin I (Myo1), Sla2/End4 or Arp2, proteins involved in the early steps of endocytosis, led to increased proliferation in low-phosphate medium compared to controls. We show that once cells are deprived of phosphate they undergo a quiescence response that is dependent on the endocytic function of Myo1. Transcriptomic analysis revealed a wide perturbation of gene expression with induction of stress-regulated genes upon phosphate starvation in wild-type but not Δmyo1 cells. Thus, endocytosis plays a pivotal role in mediating the cellular response to nutrients, bridging the external environment and internal molecular functions of the cell.

A comparison of clustering models for inference of T cell receptor antigen specificity.

The vast potential sequence diversity of TCRs and their ligands has presented an historic barrier to computational prediction of TCR epitope specificity, a holy grail of quantitative immunology. One common approach is to cluster sequences together, on the assumption that similar receptors bind similar epitopes. Here, we provide the first independent evaluation of widely used clustering algorithms for TCR specificity inference, observing some variability in predictive performance between models, and marked differences in scalability. Despite these differences, we find that different algorithms produce clusters with high degrees of similarity for receptors recognising the same epitope. Our analysis strengthens the case for use of clustering models to identify signals of common specificity from large repertoires, whilst highlighting scope for improvement of complex models over simple comparators.

Tracking in situ checkpoint inhibitor-bound target T cells in patients with checkpoint-induced colitis.

The success of checkpoint inhibitors (CPIs) for cancer has been tempered by immune-related adverse effects including colitis. CPI-induced colitis is hallmarked by expansion of resident mucosal IFNγ cytotoxic CD8+ T cells, but how these arise is unclear. Here, we track CPI-bound T cells in intestinal tissue using multimodal single-cell and subcellular spatial transcriptomics (ST). Target occupancy was increased in inflamed tissue, with drug-bound T cells located in distinct microdomains distinguished by specific intercellular signaling and transcriptional gradients. CPI-bound cells were largely CD4+ T cells, including enrichment in CPI-bound peripheral helper, follicular helper, and regulatory T cells. IFNγ CD8+ T cells emerged from both tissue-resident memory (TRM) and peripheral populations, displayed more restricted target occupancy profiles, and co-localized with damaged epithelial microdomains lacking effective regulatory cues. Our multimodal analysis identifies causal pathways and constitutes a resource to inform novel preventive strategies.

Can we predict T cell specificity with digital biology and machine learning?

Recent advances in machine learning and experimental biology have offered breakthrough solutions to problems such as protein structure prediction that were long thought to be intractable. However, despite the pivotal role of the T cell receptor (TCR) in orchestrating cellular immunity in health and disease, computational reconstruction of a reliable map from a TCR to its cognate antigens remains a holy grail of systems immunology. Current data sets are limited to a negligible fraction of the universe of possible TCR-ligand pairs, and performance of state-of-the-art predictive models wanes when applied beyond these known binders. In this Perspective article, we make the case for renewed and coordinated interdisciplinary effort to tackle the problem of predicting TCR-antigen specificity. We set out the general requirements of predictive models of antigen binding, highlight critical challenges and discuss how recent advances in digital biology such as single-cell technology and machine learning may provide possible solutions. Finally, we describe how predicting TCR specificity might contribute to our understanding of the broader puzzle of antigen immunogenicity.

A systems approach evaluating the impact of SARS-CoV-2 variant of concern mutations on CD8+ T cell responses.

T cell recognition of SARS-CoV-2 antigens after vaccination and/or natural infection has played a central role in resolving SARS-CoV-2 infections and generating adaptive immune memory. However, the clinical impact of SARS-CoV-2-specific T cell responses is variable and the mechanisms underlying T cell interaction with target antigens are not fully understood. This is especially true given the virus' rapid evolution, which leads to new variants with immune escape capacity. In this study, we used the Omicron variant as a model organism and took a systems approach to evaluate the impact of mutations on CD8+ T cell immunogenicity. We computed an immunogenicity potential score for each SARS-CoV-2 peptide antigen from the ancestral strain and Omicron, capturing both antigen presentation and T cell recognition probabilities. By comparing ancestral vs. Omicron immunogenicity scores, we reveal a divergent and heterogeneous landscape of impact for CD8+ T cell recognition of mutated targets in Omicron variants. While T cell recognition of Omicron peptides is broadly preserved, we observed mutated peptides with deteriorated immunogenicity that may assist breakthrough infection in some individuals. We then combined our scoring scheme with an in silico mutagenesis, to characterise the position- and residue-specific theoretical mutational impact on immunogenicity. While we predict many escape trajectories from the theoretical landscape of substitutions, our study suggests that Omicron mutations in T cell epitopes did not develop under cell-mediated pressure. Our study provides a generalisable platform for fostering a deeper understanding of existing and novel variant impact on antigen-specific vaccine- and/or infection-induced T cell immunity.

A robust deep learning workflow to predict CD8 + T-cell epitopes.

BACKGROUND: T-cells play a crucial role in the adaptive immune system by triggering responses against cancer cells and pathogens, while maintaining tolerance against self-antigens, which has sparked interest in the development of various T-cell-focused immunotherapies. However, the identification of antigens recognised by T-cells is low-throughput and laborious. To overcome some of these limitations, computational methods for predicting CD8 + T-cell epitopes have emerged. Despite recent developments, most immunogenicity algorithms struggle to learn features of peptide immunogenicity from small datasets, suffer from HLA bias and are unable to reliably predict pathology-specific CD8 + T-cell epitopes. METHODS: We developed TRAP (T-cell recognition potential of HLA-I presented peptides), a robust deep learning workflow for predicting CD8 + T-cell epitopes from MHC-I presented pathogenic and self-peptides. TRAP uses transfer learning, deep learning architecture and MHC binding information to make context-specific predictions of CD8 + T-cell epitopes. TRAP also detects low-confidence predictions for peptides that differ significantly from those in the training datasets to abstain from making incorrect predictions. To estimate the immunogenicity of pathogenic peptides with low-confidence predictions, we further developed a novel metric, RSAT (relative similarity to autoantigens and tumour-associated antigens), as a complementary to 'dissimilarity to self' from cancer studies. RESULTS: TRAP was used to identify epitopes from glioblastoma patients as well as SARS-CoV-2 peptides, and it outperformed other algorithms in both cancer and pathogenic settings. TRAP was especially effective at extracting immunogenicity-associated properties from restricted data of emerging pathogens and translating them onto related species, as well as minimising the loss of likely epitopes in imbalanced datasets. We also demonstrated that the novel metric termed RSAT was able to estimate immunogenic of pathogenic peptides of various lengths and species. TRAP implementation is available at: https://github.com/ChloeHJ/TRAP . CONCLUSIONS: This study presents a novel computational workflow for accurately predicting CD8 + T-cell epitopes to foster a better understanding of antigen-specific T-cell response and the development of effective clinical therapeutics.

Group A Streptococcus induces CD1a-autoreactive T cells and promotes psoriatic inflammation.

Group A Streptococcus (GAS) infection is associated with multiple clinical sequelae, including different subtypes of psoriasis. Such post-streptococcal disorders have been long known but are largely unexplained. CD1a is expressed at constitutively high levels by Langerhans cells and presents lipid antigens to T cells, but the potential relevance to GAS infection has not been studied. Here, we investigated whether GAS-responsive CD1a-restricted T cells contribute to the pathogenesis of psoriasis. Healthy individuals had high frequencies of circulating and cutaneous GAS-responsive CD4+ and CD8+ T cells with rapid effector functions, including the production of interleukin-22 (IL-22). Human skin and blood single-cell CITE-seq analyses of IL-22-producing T cells showed a type 17 signature with proliferative potential, whereas IFN-γ-producing T cells displayed cytotoxic T lymphocyte characteristics. Furthermore, individuals with psoriasis had significantly higher frequencies of circulating GAS-reactive T cells, enriched for markers of activation, cytolytic potential, and tissue association. In addition to responding to GAS, subsets of expanded GAS-reactive T cell clones/lines were found to be autoreactive, which included the recognition of the self-lipid antigen lysophosphatidylcholine. CD8+ T cell clones/lines produced cytolytic mediators and lysed infected CD1a-expressing cells. Furthermore, we established cutaneous models of GAS infection in a humanized CD1a transgenic mouse model and identified enhanced and prolonged local and systemic inflammation, with resolution through a psoriasis-like phenotype. Together, these findings link GAS infection to the CD1a pathway and show that GAS infection promotes the proliferation and activation of CD1a-autoreactive T cells, with relevance to post-streptococcal disease, including the pathogenesis and treatment of psoriasis.

Dual RNA sequencing reveals dendritic cell reprogramming in response to typhoidal Salmonella invasion.

Salmonella enterica represent a major disease burden worldwide. S. enterica serovar Typhi (S. Typhi) is responsible for potentially life-threatening Typhoid fever affecting 10.9 million people annually. While non-typhoidal Salmonella (NTS) serovars usually trigger self-limiting diarrhoea, invasive NTS bacteraemia is a growing public health challenge. Dendritic cells (DCs) are key professional antigen presenting cells of the human immune system. The ability of pathogenic bacteria to subvert DC functions and prevent T cell recognition contributes to their survival and dissemination within the host. Here, we adapted dual RNA-sequencing to define how different Salmonella pathovariants remodel their gene expression in tandem with that of infected DCs. We find DCs harness iron handling pathways to defend against invading Salmonellas, which S. Typhi is able to circumvent by mounting a robust response to nitrosative stress. In parallel, we uncover the alternative strategies invasive NTS employ to impair DC functions.

HLA variants have different preferences to present proteins with specific molecular functions which are complemented in frequent haplotypes.

Human leukocyte antigen (HLA) genes are the most polymorphic loci in the human genome and code for proteins that play a key role in guiding adaptive immune responses by presenting foreign and self peptides (ligands) to T cells. Each person carries up to 6 HLA class I variants (maternal and paternal copies of HLA-A, HLA-B and HLA-C genes) and also multiple HLA class II variants, which cumulatively define the landscape of peptides presented to T cells. Each HLA variant has its own repertoire of presented peptides with a certain sequence motif which is mainly defined by peptide anchor residues (typically the second and the last positions for HLA class I ligands) forming key interactions with the peptide-binding groove of HLA. In this study, we aimed to characterize HLA binding preferences in terms of molecular functions of presented proteins. To focus on the ligand presentation bias introduced specifically by HLA-peptide interaction we performed large-scale in silico predictions of binding of all peptides from human proteome for a wide range of HLA variants and established which functions are characteristic for proteins that are more or less preferentially presented by different HLA variants using statistical calculations and gene ontology (GO) analysis. We demonstrated marked distinctions between HLA variants in molecular functions of preferentially presented proteins (e.g. some HLA variants preferentially present membrane and receptor proteins, while others - ribosomal and DNA-binding proteins) and reduced presentation of extracellular matrix and collagen proteins by the majority of HLA variants. To explain these observations we demonstrated that HLA preferentially presents proteins enriched in amino acids which are required as anchor residues for the particular HLA variant. Our observations can be extrapolated to explain the protective effect of certain HLA alleles in infectious diseases, and we hypothesize that they can also explain susceptibility to certain autoimmune diseases and cancers. We demonstrate that these differences lead to differential presentation of HIV, influenza virus, SARS-CoV-1 and SARS-CoV-2 proteins by various HLA alleles. Taking into consideration that HLA alleles are inherited in haplotypes, we hypothesized that haplotypes composed of a combination of HLA variants with different presentation preferences should be more advantageous as they allow presenting a larger repertoire of peptides and avoiding holes in immunopeptidome. Indeed, we demonstrated that HLA-A/HLA-B and HLA-A/HLA-C haplotypes which have a high frequency in the human population are comprised of HLA variants that are more distinct in terms of functions of preferentially presented proteins than the control pairs.

Decitabine increases neoantigen and cancer testis antigen expression to enhance T-cell-mediated toxicity against glioblastoma.

BACKGROUND: Glioblastoma (GBM) is the most common and malignant primary brain tumor in adults. Despite maximal treatment, median survival remains dismal at 14-24 months. Immunotherapies, such as checkpoint inhibition, have revolutionized management of some cancers but have little benefit for GBM patients. This is, in part, due to the low mutational and neoantigen burden in this immunogenically "cold" tumor. METHODS: U87MG and patient-derived cell lines were treated with 5-aza-2'-deoxycytidine (DAC) and underwent whole-exome and transcriptome sequencing. Cell lines were then subjected to cellular assays with neoantigen and cancer testis antigen (CTA) specific T cells. RESULTS: We demonstrate that DAC increases neoantigen and CTA mRNA expression through DNA hypomethylation. This results in increased neoantigen presentation by MHC class I in tumor cells, leading to increased neoantigen- and CTA-specific T-cell activation and killing of DAC-treated cancer cells. In addition, we show that patients have endogenous cancer-specific T cells in both tumor and blood, which show increased tumor-specific activation in the presence of DAC-treated cells. CONCLUSIONS: Our work shows that DAC increases GBM immunogenicity and consequent susceptibility to T-cell responses in vitro. Our results support a potential use of DAC as a sensitizing agent for immunotherapy.

An alignment-free model for comparison of regulatory sequences.

MOTIVATION: Some recent comparative studies have revealed that regulatory regions can retain function over large evolutionary distances, even though the DNA sequences are divergent and difficult to align. It is also known that such enhancers can drive very similar expression patterns. This poses a challenge for the in silico detection of biologically related sequences, as they can only be discovered using alignment-free methods. RESULTS: Here, we present a new computational framework called Regulatory Region Scoring (RRS) model for the detection of functional conservation of regulatory sequences using predicted occupancy levels of transcription factors of interest. We demonstrate that our model can detect the functional and/or evolutionary links between some non-alignable enhancers with a strong statistical significance. We also identify groups of enhancers that are likely to be similarly regulated. Our model is motivated by previous work on prediction of expression patterns and it can capture similarity by strong binding sites, weak binding sites and even the statistically significant absence of sites. Our results support the hypothesis that weak binding sites contribute to the functional similarity of sequences. Our model fills a gap between two families of models: detailed, data-intensive models for the prediction of precise spatio-temporal expression patterns on the one side, and crude, generally applicable models on the other side. Our model borrows some of the strengths of each group and addresses their drawbacks. AVAILABILITY: The RRS source code is freely available upon publication of this manuscript: http://www2.warwick.ac.uk/fac/sci/systemsbiology/staff/ott/tools_and_software/rrs.