Single-cell multi-ome regression models identify functional and disease-associated enhancers and enable chromatin potential analysis

We present a gene-level regulatory model, single-cell ATAC + RNA linking (SCARlink), which predicts single-cell gene expression and links enhancers to target genes using multi-ome (scRNA-seq and scATAC-seq co-assay) sequencing data. The approach uses regularized Poisson regression on tile-level accessibility data to jointly model all regulatory effects at a gene locus, avoiding the limitations of pairwise gene-peak correlations and dependence on peak calling. SCARlink outperformed existing gene scoring methods for imputing gene expression from chromatin accessibility across high-coverage multi-ome datasets while giving comparable to improved performance on low-coverage datasets. Shapley value analysis on trained models identified cell-type-specific gene enhancers that are validated by promoter capture Hi-C and are 11× to 15× and 5× to 12× enriched in fine-mapped eQTLs and fine-mapped genome-wide association study (GWAS) variants, respectively. We further show that SCARlink-predicted and observed gene expression vectors provide a robust way to compute a chromatin potential vector field to enable developmental trajectory analysis.

Multicenter integrated analysis of noncoding CRISPRi screens

The ENCODE Consortium's efforts to annotate noncoding cis-regulatory elements (CREs) have advanced our understanding of gene regulatory landscapes. Pooled, noncoding CRISPR screens offer a systematic approach to investigate cis-regulatory mechanisms. The ENCODE4 Functional Characterization Centers conducted 108 screens in human cell lines, comprising >540,000 perturbations across 24.85 megabases of the genome. Using 332 functionally confirmed CRE-gene links in K562 cells, we established guidelines for screening endogenous noncoding elements with CRISPR interference (CRISPRi), including accurate detection of CREs that exhibit variable, often low, transcriptional effects. Benchmarking five screen analysis tools, we find that CASA produces the most conservative CRE calls and is robust to artifacts of low-specificity single guide RNAs. We uncover a subtle DNA strand bias for CRISPRi in transcribed regions with implications for screen design and analysis. Together, we provide an accessible data resource, predesigned single guide RNAs for targeting 3,275,697 ENCODE SCREEN candidate CREs with CRISPRi and screening guidelines to accelerate functional characterization of the noncoding genome.

Universal recording of immune cell interactions in vivo

Immune cells rely on transient physical interactions with other immune and non-immune populations to regulate their function1. To study these 'kiss-and-run' interactions directly in vivo, we previously developed LIPSTIC (labelling immune partnerships by SorTagging intercellular contacts)2, an approach that uses enzymatic transfer of a labelled substrate between the molecular partners CD40L and CD40 to label interacting cells. Reliance on this pathway limited the use of LIPSTIC to measuring interactions between CD4+ T helper cells and antigen-presenting cells, however. Here we report the development of a universal version of LIPSTIC (uLIPSTIC), which can record physical interactions both among immune cells and between immune and non-immune populations irrespective of the receptors and ligands involved. We show that uLIPSTIC can be used, among other things, to monitor the priming of CD8+ T cells by dendritic cells, reveal the steady-state cellular partners of regulatory T cells and identify germinal centre-resident T follicular helper cells on the basis of their ability to interact cognately with germinal centre B cells. By coupling uLIPSTIC with single-cell transcriptomics, we build a catalogue of the immune populations that physically interact with intestinal epithelial cells at the steady state and profile the evolution of the interactome of lymphocytic choriomeningitis virus-specific CD8+ T cells in multiple organs following systemic infection. Thus, uLIPSTIC provides a broadly useful technology for measuring and understanding cell-cell interactions across multiple biological systems.

Lessons for the Next Generation of Scientists from the Second Annual Arthur and Sandra Irving Cancer Immunology Symposium

The Arthur and Sandra Irving Cancer Immunology Symposium has been created as a platform for established cancer immunologists to mentor trainees and young investigators as they launch their research career in the field. By sharing their different paths to success, the senior faculty mentors provide an invaluable resource to support the development of the next generation of leaders in the cancer immunology community. This Commentary describes some of the key topics that were discussed during the 2022 symposium: scientific and career trajectory, leadership, mentoring, collaborations, and publishing. For each of these topics, established investigators discussed the elements that facilitate success in these areas as well as mistakes that can hinder progress. Herein, we outline the critical points raised in these discussions for establishing a successful independent research career. These points are highly relevant for the broader scientific community.

Universal recording of cell-cell contacts in vivo for interaction-based transcriptomics

Cellular interactions are essential for tissue organization and functionality. In particular, immune cells rely on direct and usually transient interactions with other immune and non-immune populations to specify and regulate their function. To study these "kiss-and-run" interactions directly in vivo, we previously developed LIPSTIC (Labeling Immune Partnerships by SorTagging Intercellular Contacts), an approach that uses enzymatic transfer of a labeled substrate between the molecular partners CD40L and CD40 to label interacting cells. Reliance on this pathway limited the use of LIPSTIC to measuring interactions between CD4+ helper T cells and antigen presenting cells, however. Here, we report the development of a universal version of LIPSTIC (uLIPSTIC), which can record physical interactions both among immune cells and between immune and non-immune populations irrespective of the receptors and ligands involved. We show that uLIPSTIC can be used, among other things, to monitor the priming of CD8+ T cells by dendritic cells, reveal the cellular partners of regulatory T cells in steady state, and identify germinal center (GC)-resident T follicular helper (Tfh) cells based on their ability to interact cognately with GC B cells. By coupling uLIPSTIC with single-cell transcriptomics, we build a catalog of the immune populations that physically interact with intestinal epithelial cells (IECs) and find evidence of stepwise acquisition of the ability to interact with IECs as CD4+ T cells adapt to residence in the intestinal tissue. Thus, uLIPSTIC provides a broadly useful technology for measuring and understanding cell-cell interactions across multiple biological systems.

Genomic Analysis of Foxp3 Function in Regulatory T Cells

Regulatory T (Treg) cells are critical for tolerance to self-antigens and for preventing autoimmunity. Foxp3 has been identified as a Treg cell lineage-defining transcription factor controlling Treg cell differentiation and function. In this article, we review the current mechanistic and systemic understanding of Foxp3 function enabled by experimental and computational advances in high-throughput genomics.

Genetic tracing reveals transcription factor Foxp3-dependent and Foxp3-independent functionality of peripherally induced Treg cells

Regulatory T (Treg) cells expressing the transcription factor Foxp3 are an essential suppressive T cell lineage of dual origin: Foxp3 induction in thymocytes and mature CD4⁺ T cells gives rise to thymic (tTreg) and peripheral (pTreg) Treg cells, respectively. While tTreg cells suppress autoimmunity, pTreg cells enforce tolerance to food and commensal microbiota. However, the role of Foxp3 in pTreg cells and the mechanisms supporting their differentiation remain poorly understood. Here, we used genetic tracing to identify microbiota-induced pTreg cells and found that many of their distinguishing features were Foxp3 independent. Lineage-committed, microbiota-dependent pTreg-like cells persisted in the colon in the absence of Foxp3. While Foxp3 was critical for the suppression of a Th17 cell program, colitis, and mastocytosis, pTreg cells suppressed colonic effector T cell expansion in a Foxp3-independent manner. Thus, Foxp3 and the tolerogenic signals that precede and promote its expression independently confer distinct facets of pTreg functionality.

Hierarchical regulation of the resting and activated T cell epigenome by major transcription factor families

T cell activation, a key early event in the adaptive immune response, is subject to elaborate transcriptional control. Here, we examined how the activities of eight major transcription factor (TF) families are integrated to shape the epigenome of naïve and activated CD4 and CD8 T cells. By leveraging extensive polymorphisms in evolutionarily divergent mice, we identified the “heavy lifters” positively influencing chromatin accessibility. Members of Ets, Runx, and TCF/Lef TF families occupied the vast majority of accessible chromatin regions, acting as “housekeepers”, “universal amplifiers”, and “placeholders”, respectively, at sites that maintained or gained accessibility upon T cell activation. Additionally, a small subset of strongly induced immune response genes displayed a non-canonical TF recruitment pattern. Our study provides a key resource and foundation for the understanding of transcriptional and epigenetic regulation in T cells and offers a new perspective on the hierarchical interactions between critical TFs.

Ectopic activation of the miR-200c-EpCAM axis enhances antitumor T cell responses in models of adoptive cell therapy

Adoptive T cell therapy (ACT) is a promising strategy for treating cancer, but it often fails because of cell intrinsic regulatory programs that limit the degree or duration of T cell function. In this study, we found that ectopic expression of microRNA-200c (miR-200c) markedly enhanced the antitumor activity of CD8+ cytotoxic T lymphocytes (CTLs) during ACT in multiple mouse models. CTLs transduced with miR-200c exhibited reduced apoptosis during engraftment and enhanced in vivo persistence, accompanied by up-regulation of the transcriptional regulator T cell factor 1 (TCF1) and the inflammatory cytokine tumor necrosis factor (TNF). miR-200c elicited these changes by suppressing the transcription factor Zeb1 and thereby inducing genes characteristic of epithelial cells. Overexpression of one of these genes, Epcam, was sufficient to augment therapeutic T cell responses against both solid and liquid tumors. These results identify the miR-200c–EpCAM axis as an avenue for improving ACT and demonstrate that select genetic perturbations can produce phenotypically distinct T cells with advantageous therapeutic properties.

A unified atlas of CD8 T cell dysfunctional states in cancer and infection

CD8 T cells play an essential role in defense against viral and bacterial infections and in tumor immunity. Deciphering T cell loss of functionality is complicated by the conspicuous heterogeneity of CD8 T cell states described across experimental and clinical settings. By carrying out a unified analysis of over 300 assay for transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-seq) experiments from 12 studies of CD8 T cells in cancer and infection, we defined a shared differentiation trajectory toward dysfunction and its underlying transcriptional drivers and revealed a universal early bifurcation of functional and dysfunctional T cell states across models. Experimental dissection of acute and chronic viral infection using single-cell ATAC (scATAC)-seq and allele-specific single-cell RNA (scRNA)-seq identified state-specific drivers and captured the emergence of similar TCF1+ progenitor-like populations at an early branch point, at which functional and dysfunctional T cells diverge. Our atlas of CD8 T cell states will facilitate mechanistic studies of T cell immunity and translational efforts.

FXR mediates T cell-intrinsic responses to reduced feeding during infection

Reduced nutrient intake is a widely conserved manifestation of sickness behavior with poorly characterized effects on adaptive immune responses. During infectious challenges, naive T cells encountering their cognate antigen become activated and differentiate into highly proliferative effector T cells. Despite their evident metabolic shift upon activation, it remains unclear how effector T cells respond to changes in nutrient availability in vivo. Here, we show that spontaneous or imposed feeding reduction during infection decreases the numbers of splenic lymphocytes. Effector T cells showed cell-intrinsic responses dependent on the nuclear receptor Farnesoid X Receptor (FXR). Deletion of FXR in T cells prevented starvation-induced loss of lymphocytes and increased effector T cell fitness in nutrient-limiting conditions, but imparted greater weight loss to the host. FXR deficiency increased the contribution of glutamine and fatty acids toward respiration and enhanced cell survival under low-glucose conditions. Provision of glucose during anorexia of infection rescued effector T cells, suggesting that this sugar is a limiting nutrient for activated lymphocytes and that alternative fuel usage may affect cell survival in starved animals. Altogether, we identified a mechanism by which the host scales immune responses according to food intake, featuring FXR as a T cell-intrinsic sensor.

High-Resolution In Vivo Identification of miRNA Targets by Halo-Enhanced Ago2 Pull-Down

The identification of microRNA (miRNA) targets by Ago2 crosslinking-immunoprecipitation (CLIP) methods has provided major insights into the biology of this important class of non-coding RNAs. However, these methods are technically challenging and not easily applicable to an in vivo setting. To overcome these limitations and facilitate the investigation of miRNA functions in vivo, we have developed a method based on a genetically engineered mouse harboring a conditional Halo-Ago2 allele expressed from the endogenous Ago2 locus. By using a resin conjugated to the HaloTag ligand, Ago2-miRNA-mRNA complexes can be purified from cells and tissues expressing the endogenous Halo-Ago2 allele. We demonstrate the reproducibility and sensitivity of this method in mouse embryonic stem cells, developing embryos, adult tissues, and autochthonous mouse models of human brain and lung cancers. This method and the datasets we have generated will facilitate the characterization of miRNA-mRNA networks in vivo under physiological and pathological conditions.

Non-redundant ISGF3 Components Promote NK Cell Survival in an Auto-regulatory Manner during Viral Infection

Natural killer (NK) cells are innate lymphocytes that possess adaptive features, including antigen-specific clonal expansion and long-lived memory responses. Although previous work demonstrated that type I interferon (IFN) signaling is crucial for NK cell expansion and memory cell formation following mouse cytomegalovirus (MCMV) infection, the global transcriptional mechanisms underlying type I IFN-mediated responses remained to be determined. Here, we demonstrate that among the suite of transcripts induced in activated NK cells, IFN-α is necessary and sufficient to promote expression of its downstream transcription factors STAT1, STAT2, and IRF9, via an auto-regulatory, feedforward loop. Similar to STAT1 deficiency, we show that STAT2- or IRF9-deficient NK cells are defective in their ability to expand following MCMV infection, in part because of diminished survival rather than an inability to proliferate. Thus, our findings demonstrate that individual ISGF3 components are crucial cell-autonomous and non-redundant regulators of the NK cell response to viral infection.

Using RNA-seq and ChIP-seq, Geary et al. investigate the impacts of type I interferon on NK cells during MCMV infection and demonstrate crucial and non-redundant roles for STAT1, STAT2, and IRF9 in promoting cytotoxicity and survival of antiviral NK cells.

Transcriptional Basis of Mouse and Human Dendritic Cell Heterogeneity

Dendritic cells (DCs) play a critical role in orchestrating adaptive immune responses due to their unique ability to initiate T cell responses and direct their differentiation into effector lineages. Classical DCs have been divided into two subsets, cDC1 and cDC2, based on phenotypic markers and their distinct abilities to prime CD8 and CD4 T cells. While the transcriptional regulation of the cDC1 subset has been well characterized, cDC2 development and function remain poorly understood. By combining transcriptional and chromatin analyses with genetic reporter expression, we identified two principal cDC2 lineages defined by distinct developmental pathways and transcriptional regulators, including T-bet and RORγt, two key transcription factors known to define innate and adaptive lymphocyte subsets. These novel cDC2 lineages were characterized by distinct metabolic and functional programs. Extending our findings to humans revealed conserved DC heterogeneity and the presence of the newly defined cDC2 subsets in human cancer.

•

Single-cell analyses reveal novel dendritic cell subsets

•

Major cDC2 subsets differentially express T-bet and RORγt

•

Distinct pro- and anti-inflammatory potential of T-bet⁺ and Tbet^– cDC2s

•

Transcriptional basis for cDC2 heterogeneity conserved across mouse and human

Abstract 3384: Genome-wide epigenetic and transcriptional comparison of CD8 T cell functional states between mouse models of cancer and chronic viral infection

Recent studies in mouse models of cancer and chronic viral infection applied genomic assays to profile various states of impaired T cell function and demonstrated that T cell dysfunction is epigenetically imprinted. However, comprehensive characterization of T cell dysfunction across models based on their epigenetic and transcriptional profiles is lacking. We reanalyzed a large collection of recently published chromatin accessibility (ATAC-seq) and gene expression (RNA-seq) data sets. After batch effect correction, we observed that epigenetic profiles of dysfunctional tumor-infiltrating T cells and exhausted T cells in chronic viral infection were surprisingly extremely similar. Furthermore, we observed that temporal progression towards dysfunction in chronic infection resembles that in tumorigenesis, and found that T cells committed to becoming dysfunctional early after activation. Motif analysis using generalized linear modeling allowed us to find candidate transcription factors associated with development of dysfunction in both immune settings. This analysis provides a better systematic understanding of cell-intrinsic mechanisms driving different functional states of CD8 T cells and demonstrates the power of systems biology approaches in cancer immunology.

Citation Format: Yuri Pritykin, Christina Leslie. Genome-wide epigenetic and transcriptional comparison of CD8 T cell functional states between mouse models of cancer and chronic viral infection [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3384.

Abstract B038: A unified genome-wide analysis of dysfunctional T-cell states in cancer and chronic viral infection

Tumor-specific T-cells that have differentiated into a terminal dysfunctional state exist in the tumor microenvironment. A systematic understanding of the requirements of immunotherapeutic rescue of these cells is critically needed to improve clinical results in patients. Mouse models of chronic infection and cancer have been studied to elucidate biologic mechanisms of persistent antigen stimulation resulting in T-cell dysfunction, or “exhaustion.” Recently, chromatin accessibility imprinting has been associated with T-cells falling back into the dysfunctional state after temporary rescue by checkpoint blockade, suggesting epigenetic mechanisms in control of T-cell dysfunction. However, comprehensive characterization of T-cell dysfunction across models based on their epigenetic and transcriptional profiles is lacking.We collected 106 chromatin accessibility (ATAC-seq) samples and 87 gene expression (RNA-seq) samples from seven recent publications. We analyzed these data by first applying batch effect correction using generalized linear modeling. This enabled mapping profiles of chromatin accessibility peaks in gene promoters and enhancers from different studies into the same space. We observed that epigenetic profiles of dysfunctional tumor-infiltrating T-cells and dysfunctional T-cells in chronic viral infection were, surprisingly, extremely similar. Furthermore, a recently characterized discrete distinction between epigenetic profiles of early (day 7-8) and late (day 28-35) dysfunction in the tumor was recapitulated in the model of chronic infection. Overall we observed across mouse models that T-cells committed to becoming dysfunctional early after an immune challenge, rather than first mounting and then loosing an effector response. These observations were also largely recapitulated in gene expression analysis. Differentially expressed genes with massive differential accessibility of their promoter and enhancer peaks during development of dysfunction, observed consistently across models, including transcription factors (TF) well studied in immunity such as Tcf7, Lef1, Satb1, Ikzf2, Tox, are good candidates for further targeted analysis.We then turned to TF binding analysis. We associated absolute levels of chromatin accessibility in peaks of each sample with TF binding (predicted by motif analysis) using regularized negative binomial regression with cross-validation. We estimated the effect of each TF in each sample, which allowed us to map chromatin accessibility profiles into the TF activity space of much lower dimensionality. This mapping largely preserved the hierarchy of relative similarities between samples. We identified key TFs whose binding was associated with open or closed chromatin in functional and dysfunctional cell states. For example, not surprisingly, binding of well known effector factors Eomes and Batf was associated with closed chromatin in naive cells and open chromatin in effector cells. Strikingly, the strongest association with closing chromatin in dysfunction, consistently across mouse models, was observed for Tcf7/Lef1 binding, further suggesting the role of these TFs in establishing the terminal CD8 T-cell dysfunctional state. Notably, we found two large groups of TFs anti-correlated to each other whose predicted binding sites had higher accessibility in either functional or dysfunctional cells. This suggested that coordinated activity of a broad range of TFs (not necessarily binding at the same sites) might be responsible for establishing and maintaining T-cell functional states, and focusing on one or a handful of TFs is not sufficient to explain them. This analysis provides a better systematic understanding of cell-intrinsic mechanisms driving different functional states of CD8 T-cells, and the developed computational methods are broadly applicable in other experimental setups where diverse cell states are profiled by high-throughput genomic assays.

Citation Format: Yuri Pritykin, Christina Leslie. A unified genome-wide analysis of dysfunctional T-cell states in cancer and chronic viral infection [abstract]. In: Proceedings of the Fourth CRI-CIMT-EATI-AACR International Cancer Immunotherapy Conference: Translating Science into Survival; Sept 30-Oct 3, 2018; New York, NY. Philadelphia (PA): AACR; Cancer Immunol Res 2019;7(2 Suppl):Abstract nr B038.

Transcription factor Foxp1 regulates Foxp3 chromatin binding and coordinates regulatory T cell function

Regulatory T cells (T_reg cells), whose differentiation and function are controlled by transcription factor Foxp3, express the closely related family member Foxp1. Here we explored Foxp1 function in T_reg cells. We found that a large number of Foxp3-bound genomic sites in T_reg cells were occupied by Foxp1 in both T_reg cells and conventional T cells (T_conv cells). In T_reg cells, Foxp1 markedly increased Foxp3 binding to these sites. Foxp1 deficiency in T_reg cells resulted in their impaired function and competitive fitness, associated with markedly reduced CD25 expression and interleukin-2 (IL-2) responsiveness, diminished CTLA-4 expression and increased SATB1 expression. The characteristic expression patterns of CD25, Foxp3 and CTLA-4 in T_reg cells were fully or partially rescued by strong IL-2 signaling. Our studies suggest that Foxp1 serves an essential non-redundant function in T_reg cells by enforcing Foxp3-mediated regulation of gene expression and enabling efficient IL-2 signaling in these cells.

ILC1 Confer Early Host Protection at Initial Sites of Viral Infection

Infection is restrained by the concerted activation of tissue-resident and circulating immune cells. Whether tissue-resident lymphocytes confer early antiviral immunity at local sites of primary infection prior to the initiation of circulating responses is not well understood. Furthermore, the kinetics of initial antiviral responses at sites of infection remain unclear. Here, we show that tissue-resident type 1 innate lymphoid cells (ILC1) serve an essential early role in host immunity through rapid production of interferon (IFN)-γ following viral infection. Ablation of Zfp683-dependent liver ILC1 lead to increased viral load in the presence of intact adaptive and innate immune cells critical for mouse cytomegalovirus (MCMV) clearance. Swift production of interleukin (IL)-12 by tissue-resident XCR1⁺ conventional dendritic cells (cDC1) promoted ILC1 production of IFN-γ in a STAT4-dependent manner to limit early viral burden. Thus, ILC1 contribute an essential role in viral immunosurveillance at sites of initial infection in response to local cDC1-derived proinflammatory cytokines.

Integrative analysis unveils new functions for the Drosophila Cutoff protein in noncoding RNA biogenesis and gene regulation

Piwi-interacting RNAs (piRNAs) are central components of the piRNA pathway, which directs transposon silencing and guarantees genome integrity in the germ cells of several metazoans. In Drosophila, piRNAs are produced from discrete regions of the genome termed piRNA clusters, whose expression relies on the RDC complex comprised of the core proteins Rhino, Deadlock, and Cutoff. To date, the RDC complex has been exclusively implicated in the regulation of the piRNA loci. Here we further elucidate the function of Cutoff and the RDC complex by performing genome-wide ChIP-seq and RNA-seq assays in the Drosophila ovaries and analyzing these data together with other publicly available data sets. In agreement with previous studies, we confirm that Cutoff is involved in the transcriptional regulation of piRNA clusters and in the repression of transposable elements in germ cells. Surprisingly, however, we find that Cutoff is enriched at and affects the expression of other noncoding RNAs, including spliceosomal RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). At least in some instances, Cutoff appears to act at a transcriptional level in concert with Rhino and perhaps Deadlock. Finally, we show that mutations in Cutoff result in the deregulation of hundreds of protein-coding genes in germ cells. Our study uncovers a broader function for the RDC complex in the Drosophila germline development.

GuideScan software for improved single and paired CRISPR guide RNA design

We present GuideScan software for the design of CRISPR guide RNA libraries that can be used to edit coding and noncoding genomic regions. GuideScan produces high-density sets of guide RNAs (gRNAs) for single- and paired-gRNA genome-wide screens. We also show that the trie data structure of GuideScan enables the design of gRNAs that are more specific than those designed by existing tools.

Cancer Immunosurveillance by Tissue-resident Innate Lymphoid Cells and Innate-like T Cells

Malignancy can be suppressed by the immune system in a process termed immunosurveillance. However, to what extent immunosurveillance occurs in spontaneous cancers and the composition of participating cell types remain obscure. Here we show that cell transformation triggers a tissue-resident lymphocyte response in oncogene-induced murine cancer models. Non-circulating cytotoxic lymphocytes, derived from innate, TCRαβ and TCRγδ lineages, expand in early tumors. Characterized by high expression of NK1.1, CD49a and CD103, these cells share a gene expression signature distinct from those of conventional NK cells, T cells and invariant NKT cells. Generation of these lymphocytes is dependent on the cytokine IL-15, but not the transcription factor Nfil3 that is required for the differentiation of tumor-infiltrating NK cells, and IL-15, but not Nfil3, deficiency results in accelerated tumor growth. These findings reveal a tumor-elicited immunosurveillance mechanism that engages unconventional type 1-like innate lymphoid cells and type 1 innate-like T cells.

Cancer Immunosurveillance by Tissue-Resident Innate Lymphoid Cells and Innate-like T Cells

Malignancy can be suppressed by the immune system in a process termed immunosurveillance. However, to what extent immunosurveillance occurs in spontaneous cancers and the composition of participating cell types remains obscure. Here, we show that cell transformation triggers a tissue-resident lymphocyte response in oncogene-induced murine cancer models. Non-circulating cytotoxic lymphocytes, derived from innate, T cell receptor (TCR)αβ, and TCRγδ lineages, expand in early tumors. Characterized by high expression of NK1.1, CD49a, and CD103, these cells share a gene-expression signature distinct from those of conventional NK cells, T cells, and invariant NKT cells. Generation of these lymphocytes is dependent on the cytokine IL-15, but not the transcription factor Nfil3 that is required for the differentiation of tumor-infiltrating NK cells, and IL-15 deficiency, but not Nfil3 deficiency, results in accelerated tumor growth. These findings reveal a tumor-elicited immunosurveillance mechanism that engages unconventional type-1-like innate lymphoid cells and type 1 innate-like T cells.

Genome-Wide Detection and Analysis of Multifunctional Genes

Many genes can play a role in multiple biological processes or molecular functions. Identifying multifunctional genes at the genome-wide level and studying their properties can shed light upon the complexity of molecular events that underpin cellular functioning, thereby leading to a better understanding of the functional landscape of the cell. However, to date, genome-wide analysis of multifunctional genes (and the proteins they encode) has been limited. Here we introduce a computational approach that uses known functional annotations to extract genes playing a role in at least two distinct biological processes. We leverage functional genomics data sets for three organisms—H. sapiens, D. melanogaster, and S. cerevisiae—and show that, as compared to other annotated genes, genes involved in multiple biological processes possess distinct physicochemical properties, are more broadly expressed, tend to be more central in protein interaction networks, tend to be more evolutionarily conserved, and are more likely to be essential. We also find that multifunctional genes are significantly more likely to be involved in human disorders. These same features also hold when multifunctionality is defined with respect to molecular functions instead of biological processes. Our analysis uncovers key features about multifunctional genes, and is a step towards a better genome-wide understanding of gene multifunctionality.

Almost every aspect of cellular function depends on protein activity. In spite of being fine-tuned to carry out highly specific functions, proteins can also multitask. Experimental studies have identified genes and proteins endowed with more than one molecular function, or participating in very different biological processes. These studies suggest that the degree of functional plasticity exhibited by proteins might go well beyond a simple “one protein—one function” relationship. However, systematic studies of the properties of multifunctional genes (and their encoded proteins) have been limited. Here we present a computational framework to identify putative multifunctional genes, and compare their properties with those of other genes. We find that multifunctional genes are significantly different from other genes with respect to their physicochemical properties, expression profiles, and interaction properties. We also observe that multifunctional genes tend to be more conserved, and that a greater fraction of them are associated with human disorders. Taken together, these results represent a step towards a more complete understanding of the role multifunctional genes play in the functional organization of the cell.