The Hallmarks of Precancer

Research on precancers, as defined as at-risk tissues and early lesions, is of high significance given the effectiveness of early intervention. We discuss the need for risk stratification to prevent overtreatment, an emphasis on the role of genetic and epigenetic aging when considering risk, and the importance of integrating macroenvironmental risk factors with molecules and cells in lesions and at-risk normal tissues for developing effective intervention and health policy strategies.

Single-cell transcriptional landscape of temporal neutrophil response to burn wound in larval zebrafish

Neutrophils accumulate early in tissue injury. However, the cellular and functional heterogeneity of neutrophils during homeostasis and in response to tissue damage remains unclear. Here, we use larval zebrafish to understand neutrophil responses to thermal injury. Single-cell transcriptional mapping of myeloid cells during a 3-day time course in burn and control larvae revealed distinct neutrophil subsets and their cell-cell interactions with macrophages across time and conditions. The trajectory formed by three zebrafish neutrophil subsets resembles human neutrophil maturation, with varying transition patterns between conditions. Through ligand-receptor cell-cell interaction analysis, we found neutrophils communicate more in burns in a pathway and temporal manner. Finally, we identified the correlation between zebrafish myeloid signatures and human burn severity, establishing GPR84+ neutrophils as a potential marker of early innate immune response in burns. This work builds the molecular foundation and a comparative single-cell genomic framework to identify neutrophil markers of tissue damage using model organisms.

Fallopian tube single cell analysis reveals myeloid cell alterations in high-grade serous ovarian cancer

Most high-grade serous ovarian cancers (HGSCs) likely initiate from fallopian tube (FT) epithelia. While epithelial subtypes have been characterized using single-cell RNA-sequencing (scRNA-Seq), heterogeneity of other compartments and their involvement in tumor progression are poorly defined. Integrated analysis of human FT scRNA-Seq and HGSC-related tissues, including tumors, revealed greater immune and stromal transcriptional diversity than previously reported. We identified abundant monocytes in FTs across two independent cohorts. The ratio of macrophages to monocytes is similar between benign FTs, ovaries, and adjacent normal tissues but significantly greater in tumors. FT-defined monocyte and macrophage signatures, cell-cell communication, and gene set enrichment analyses identified monocyte- and macrophage-specific interactions and functional pathways in paired tumors and adjacent normal tissues. Further reanalysis of HGSC scRNA-Seq identified monocyte and macrophage subsets associated with neoadjuvant chemotherapy. Taken together, our work provides data that an altered FT myeloid cell composition could inform the discovery of early detection markers for HGSC.

PTP1B phosphatase dampens iPSC-derived neutrophil motility and antimicrobial function

Neutrophils are rapidly recruited to sites of infection and are critical for pathogen clearance. Therapeutic use of primary neutrophils has been limited as they have a short lifespan and are not amenable to genetic manipulation. Human induced pluripotent stem cells (iPSCs) can provide a robust source of neutrophils for infusion and are genetically tractable. However, current work has indicated that dampened intracellular signaling limits iPSC-derived neutrophil (iNeutrophil) cellular activation and antimicrobial response. Here, we show that protein tyrosine phosphatase 1B (PTP1B) inhibits intracellular signaling and dampens iNeutrophil effector function. Deletion of the PTP1B phosphatase increased PI3K and ERK signaling and was associated with increased F-actin polymerization, cell migration and phagocytosis. In contrast, other effector functions like NETosis and ROS production were reduced. PTP1B-deficient neutrophils were more responsive to A. fumigatus and displayed rapid recruitment and control of hyphal growth. Accordingly, depletion of PTP1B increased production of inflammatory factors including the neutrophil chemokine IL-8. Taken together, these findings suggest that PTP1B limits iNeutrophil motility and antimicrobial function.

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.

Human CD79b+ neutrophils in the blood are associated with early-stage melanoma

Purpose

Due to their abundance in the blood, low RNA content, and short lifespan, neutrophils have been classically considered to be one homogenous pool. However, recent work has found that mature neutrophils and neutrophil progenitors are composed of unique subsets exhibiting context-dependent functions. In this study, we ask if neutrophil heterogeneity is associated with melanoma incidence and/or disease stage.

Experimental design

Using mass cytometry, we profiled melanoma patient blood for unique cell surface markers among neutrophils. Markers were tested for their predictiveness using flow cytometry data and random forest machine learning.

Results

We identified CD79b+ neutrophils (CD3-CD56-CD19-Siglec8-CD203c-CD86LoCD66b+CD79b+) that are normally restricted to the bone marrow in healthy humans but appear in the blood of subjects with early-stage melanoma. Further, we found CD79b+ neutrophils present in tumors of subjects with head and neck cancer. AI-mediated machine learning analysis of neutrophils from subjects with melanoma confirmed that CD79b expression among peripheral blood neutrophils is highly important in identifying melanoma incidence. We noted that CD79b+ neutrophils possessed a neutrophilic appearance but have transcriptional and surface-marker phenotypes reminiscent of B cells. Compared to remaining blood neutrophils, CD79b+ neutrophils are primed for NETosis, express higher levels of antigen presentation-related proteins, and have an increased capacity for phagocytosis.

Conclusion

Our work suggests that CD79b+ neutrophils are associated with early-stage melanoma.

Stress Keratin 17 Is a Predictive Biomarker Inversely Associated with Response to Immune Check-Point Blockade in Head and Neck Squamous Cell Carcinomas and Beyond

Low response rates in immune check-point blockade (ICB)-treated head and neck squamous cell carcinoma (HNSCC) drive a critical need for robust, clinically validated predictive biomarkers. Our group previously showed that stress keratin 17 (CK17) suppresses macrophage-mediated CXCL9/CXCL10 chemokine signaling involved in attracting activated CD8+ T cells into tumors, correlating with decreased response rate to pembrolizumab-based therapy in a pilot cohort of ICB-treated HNSCC (n = 26). Here, we performed an expanded analysis of the predictive value of CK17 in ICB-treated HNSCC according to the REMARK criteria and investigated the gene expression profiles associated with high CK17 expression. Pretreatment samples from pembrolizumab-treated HNSCC patients were stained via immunohistochemistry using a CK17 monoclonal antibody (n = 48) and subjected to spatial transcriptomic profiling (n = 8). Our findings were validated in an independent retrospective cohort (n = 22). CK17 RNA expression in pembrolizumab-treated patients with various cancer types was investigated for predictive significance. Of the 48 patients (60% male, median age of 61.5 years), 21 (44%) were CK17 high, and 27 (56%) were CK17 low. A total of 17 patients (35%, 77% CK17 low) had disease control, while 31 patients (65%, 45% CK17 low) had progressive disease. High CK17 expression was associated with a lack of disease control (p = 0.037), shorter time to treatment failure (p = 0.025), and progression-free survival (PFS, p = 0.004), but not overall survival (OS, p = 0.06). A high CK17 expression was associated with lack of disease control in an independent validation cohort (p = 0.011). PD-L1 expression did not correlate with CK17 expression or clinical outcome. CK17 RNA expression was predictive of PFS and OS in 552 pembrolizumab-treated cancer patients. Our findings indicate that high CK17 expression may predict resistance to ICB in HNSCC patients and beyond.

Cell-autonomous regulation of complement C3 by factor H limits macrophage efferocytosis and exacerbates atherosclerosis

Complement factor H (CFH) negatively regulates consumption of complement component 3 (C3), thereby restricting complement activation. Genetic variants in CFH predispose to chronic inflammatory disease. Here, we examined the impact of CFH on atherosclerosis development. In a mouse model of atherosclerosis, CFH deficiency limited plaque necrosis in a C3-dependent manner. Deletion of CFH in monocyte-derived inflammatory macrophages propagated uncontrolled cell-autonomous C3 consumption without downstream C5 activation and heightened efferocytotic capacity. Among leukocytes, Cfh expression was restricted to monocytes and macrophages, increased during inflammation, and coincided with the accumulation of intracellular C3. Macrophage-derived CFH was sufficient to dampen resolution of inflammation, and hematopoietic deletion of CFH in atherosclerosis-prone mice promoted lesional efferocytosis and reduced plaque size. Furthermore, we identified monocyte-derived inflammatory macrophages expressing C3 and CFH in human atherosclerotic plaques. Our findings reveal a regulatory axis wherein CFH controls intracellular C3 levels of macrophages in a cell-autonomous manner, evidencing the importance of on-site complement regulation in the pathogenesis of inflammatory diseases.

Nonclassical monocytes potentiate anti-tumoral CD8+ T cell responses in the lungs

CD8⁺ T cells drive anti-cancer immunity in response to antigen-presenting cells such as dendritic cells and subpopulations of monocytes and macrophages. While CD14⁺ classical monocytes modulate CD8⁺ T cell responses, the contributions of CD16⁺ nonclassical monocytes to this process remain unclear. Herein we explored the role of nonclassical monocytes in CD8⁺ T cell activation by utilizing E2-deficient (E2^-/-) mice that lack nonclassical monocytes. During early metastatic seeding, modeled by B16F10-OVA cancer cells injected into E2^-/- mice, we noted lower CD8⁺ effector memory and effector T cell frequencies within the lungs as well as in lung-draining mediastinal lymph nodes in the E2^-/- mice. Analysis of the myeloid compartment revealed that these changes were associated with depletion of MHC-II^loLy6C^lo nonclassical monocytes within these tissues, with little change in other monocyte or macrophage populations. Additionally, nonclassical monocytes preferentially trafficked to primary tumor sites in the lungs, rather than to the lung-draining lymph nodes, and did not cross-present antigen to CD8⁺ T cells. Examination of the lung microenvironment in E2^-/- mice revealed reduced CCL21 expression in endothelial cells, which is chemokine involved in T cell trafficking. Our results highlight the previously unappreciated importance of nonclassical monocytes in shaping the tumor microenvironment via CCL21 production and CD8⁺ T cell recruitment.

A tessellated lymphoid network provides whole-body T cell surveillance in zebrafish

Homeostatic trafficking to lymph nodes allows T cells to efficiently survey the host for cognate antigen. Nonmammalian jawed vertebrates lack lymph nodes but maintain diverse T cell pools. Here, we exploit in vivo imaging of transparent zebrafish to investigate how T cells organize and survey for antigen in an animal devoid of lymph nodes. We find that naïve-like T cells in zebrafish organize into a previously undescribed whole-body lymphoid network that supports streaming migration and coordinated trafficking through the host. This network has the cellular hallmarks of a mammalian lymph node, including naïve T cells and CCR7-ligand expressing nonhematopoietic cells, and facilitates rapid collective migration. During infection, T cells transition to a random walk that supports antigen-presenting cell interactions and subsequent activation. Our results reveal that T cells can toggle between collective migration and individual random walks to prioritize either large-scale trafficking or antigen search in situ. This lymphoid network thus facilitates whole-body T cell trafficking and antigen surveillance in the absence of a lymph node system.

Single-cell transcriptomic analysis of endometriosis

Endometriosis is a common condition in women that causes chronic pain and infertility and is associated with an elevated risk of ovarian cancer. We profiled transcriptomes of >370,000 individual cells from endometriomas (n = 8), endometriosis (n = 28), eutopic endometrium (n = 10), unaffected ovary (n = 4) and endometriosis-free peritoneum (n = 4), generating a cellular atlas of endometrial-type epithelial cells, stromal cells and microenvironmental cell populations across tissue sites. Cellular and molecular signatures of endometrial-type epithelium and stroma differed across tissue types, suggesting a role for cellular restructuring and transcriptional reprogramming in the disease. Epithelium, stroma and proximal mesothelial cells of endometriomas showed dysregulation of pro-inflammatory pathways and upregulation of complement proteins. Somatic ARID1A mutation in epithelial cells was associated with upregulation of pro-angiogenic and pro-lymphangiogenic factors and remodeling of the endothelial cell compartment, with enrichment of lymphatic endothelial cells. Finally, signatures of ciliated epithelial cells were enriched in ovarian cancers, reinforcing epidemiologic associations between these two diseases.

A tessellated lymphoid network provides whole-body T cell surveillance in zebrafish

Homeostatic trafficking to lymph nodes allows T cells to efficiently survey the host for cognate antigen. Non-mammalian jawed vertebrates lack lymph nodes but maintain similarly diverse T cell pools. Here, we exploit in vivo imaging of transparent zebrafish to investigate how T cells organize and survey for antigen in an animal devoid of lymph nodes. We find that naïve-like T cells in zebrafish organize into a previously undescribed whole-body lymphoid network that supports streaming migration and coordinated trafficking through the host. This network has the cellular hallmarks of a mammalian lymph node, including naïve T cells and CCR7-ligand expressing non-hematopoietic cells, and facilitates rapid collective migration. During infection, T cells transition to a random walk that supports antigen presenting cell interactions and subsequent activation. Our results reveal that T cells can toggle between collective migration and individual random walks to prioritize either large-scale trafficking or antigen search in situ . This novel lymphoid network thus facilitates whole-body T cell trafficking and antigen surveillance in the absence of a lymph node system.

Significance Statement

In mammals, lymph nodes play a critical role in the initiation of adaptive immune responses by providing a dedicated place for T cells to scan antigen-presenting cells. Birds, reptiles, amphibians, and fish all maintain diverse repertoires of T cells but lack lymph nodes, raising questions about how adaptive immunity functions in lower jawed vertebrates. Here, we describe a novel network of lymphocytes in zebrafish that supports whole-body T cell trafficking and provides a site for antigen search, mirroring the function of mammalian lymph nodes. Within this network, T cells can prioritize large-scale trafficking or antigen scanning by toggling between two distinct modes of migration. This network provides valuable insights into the evolution of adaptive immunity.

Single-cell immune profiling reveals long-term changes in myeloid cells and identifies a novel subset of CD9+ monocytes associated with COVID-19 hospitalization

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection can result in severe immune dysfunction, hospitalization, and death. Many patients also develop long-COVID-19, experiencing symptoms months after infection. Although significant progress has been made in understanding the immune response to acute SARS-CoV-2 infection, gaps remain in our knowledge of how innate immunity influences disease kinetics and severity. We hypothesized that cytometry by time-of-flight analysis of PBMCs from healthy and infected subjects would identify novel cell surface markers and innate immune cell subsets associated with COVID-19 severity. In this pursuit, we identified monocyte and dendritic cell subsets that changed in frequency during acute SARS-CoV-2 infection and correlated with clinical parameters of disease severity. Subsets of nonclassical monocytes decreased in frequency in hospitalized subjects, yet increased in the most severe patients and positively correlated with clinical values associated with worse disease severity. CD9, CD163, PDL1, and PDL2 expression significantly increased in hospitalized subjects, and CD9 and 6-Sulfo LacNac emerged as the markers that best distinguished monocyte subsets amongst all subjects. CD9+ monocytes remained elevated, whereas nonclassical monocytes remained decreased, in the blood of hospitalized subjects at 3-4 months postinfection. Finally, we found that CD9+ monocytes functionally released more IL-8 and MCP-1 after LPS stimulation. This study identifies new monocyte subsets present in the blood of COVID-19 patients that correlate with disease severity, and links CD9+ monocytes to COVID-19 progression.

Stress Keratin 17 Expression in Head and Neck Cancer Contributes to Immune Evasion and Resistance to Immune-Checkpoint Blockade

PURPOSE

We investigated whether in human head and neck squamous cell carcinoma (HNSCC) high levels of expression of stress keratin 17 (K17) are associated with poor survival and resistance to immunotherapy.

EXPERIMENTAL DESIGN

We investigated the role of K17 in regulating both the tumor microenvironment and immune responsiveness of HNSCC using a syngeneic mouse HNSCC model, MOC2. MOC2 gives rise to immunologically cold tumors that are resistant to immune-checkpoint blockade (ICB). We engineered multiple, independent K17 knockout (KO) MOC2 cell lines and monitored their growth and response to ICB. We also measured K17 expression in human HNSCC of patients undergoing ICB.

RESULTS

MOC2 tumors were found to express K17 at high levels. When knocked out for K17 (K17KO MOC2), these cells formed tumors that grew slowly or spontaneously regressed and had a high CD8+ T-cell infiltrate in immunocompetent syngeneic C57BL/6 mice compared with parental MOC2 tumors. This phenotype was reversed when we depleted mice for T cells. Whereas parental MOC2 tumors were resistant to ICB treatment, K17KO MOC2 tumors that did not spontaneously regress were eliminated upon ICB treatment. In a cohort of patients with HNSCC receiving pembrolizumab, high K17 expression correlated with poor response. Single-cell RNA-sequencing analysis revealed broad differences in the immune landscape of K17KO MOC2 tumors compared with parental MOC2 tumors, including differences in multiple lymphoid and myeloid cell types.

CONCLUSIONS

We demonstrate that K17 expression in HNSCC contributes to immune evasion and resistance to ICB treatment by broadly altering immune landscapes of tumors.

Abstract 5148: Immune landscapes and interactions in human fallopian tubes and their implications for early high-grade serous ovarian cancers

The five-year survival rate of ovarian cancer (OC) dramatically rises from below 40% to 90% with patients who are diagnosed early (stage I), when the disease restricts to the ovary. Thus, the ultimate goals of our research are to learn how to detect ovarian cancer at its earliest stages of development and thereby improve patient survival and to learn how the immune microenvironment contributes to the development of OC. Achieving these goals requires us to better characterize the earliest steps in ovarian carcinogenesis. The human fallopian tube harbors the cell-of-origin for most high-grade serous 'ovarian' cancers (HGSOCs), but its cellular composition, particularly the immune component, is poorly characterized. Emerging evidence from several research groups using molecular sequencing technologies of fallopian tube epithelium strongly support this hypothesis. We recently described our use of single-cell transcriptome technology (scRNA-Seq) to profile the heterogeneity of non-malignant FT and define a novel transcriptional gene regulatory network of epithelial differentiation and tumorigenesis in the human FT. Here we present an in-depth characterization of the immune cell atlas for human fallopian tubes, which has not been reported. Our unpublished analyses of our FT dataset and the publicly available HGSOC scRNA-Seq dataset demonstrate a significant difference in the immune cell landscape. There are limited macrophages and a high abundance of monocytes in non-malignant FT and adjacent normal ovaries compared to those in HGSOC tumors. This finding leads us to hypothesize that the macrophage-to-monocyte ratio and their interactions with other immune cell types in human fallopian tubes create a permissive microenvironment allowing for the development of HGSOC and can provide a biomarker for early cancer detection. To test if those myeloid cells could be a marker for early HGSOC, we performed scRNA-seq on FT from women carrying mutant alleles of BRCA1 or BRCA2 who are at high risk of HGSOC. Intriguingly, the ratio of monocytes to macrophages in the FTs from half of these high-risk individuals is more like that found in HGSOC than in normal FT. Furthermore, computational analyses on myeloid cell types and their interactions with different epithelial cells, stromal fibroblast cells, and T and NK cells revealed the extended contribution of the non-epithelial microenvironment in FT, potentially driving early HGSOC development and progression. This cellular and molecular compendium of the human fallopian tube immune landscapes in cancer-free and cancer-high-risk women is expected to advance our understanding of the earliest stages of fallopian epithelial neoplasia.

Citation Format: Josh Brand, Marcela Harro, Forough Abbasi, Xianzhi Lin, Matthew Siedhoff, Li Andrew, Bobbie J. Rimel, Fabiola Medeiros, Kate Lawrenson, Huy Q. Dinh. Immune landscapes and interactions in human fallopian tubes and their implications for early high-grade serous ovarian cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5148.

Abstract 1994: Defining tumor infiltrating myeloid microenvironment landscapes in HPV+ and HPV- head and neck cancers

Head and neck cancer (HNC) has a <20% response rate to immunotherapy, making it imperative to learn what causes high treatment resistance and how to overcome these challenges. One way to generate a holistic view of both human papillomavirus positive (HPV+) and negative (HPV-) HNC is to study the tumor immune microenvironment (TIME). Since HPV+ HNC often results in more favorable outcomes than HPV- HNC, assessing differences between these two disease states may provide critical mechanistic insight into what causes discrepancies in treatment response. This study uses publicly available and in-house collected single-cell transcriptome (scRNA-seq) data from preclinical murine models and human HNC to assess the differences between HPV+ and HPV- HNC TIME. Recent scRNA-Seq data showed extensive characterization of infiltrating T and B cells between HPV+ and HPV- HNC, but little is known about the myeloid compartment. As myeloid cells with both pro- and anti-tumoral functions make up a bulk of the HNC TIME, we sought to identify the differences in myeloid cells between HPV+ vs. HPV- HNC and their roles in response to immune checkpoint inhibitor (ICI) therapy. We identified higher myeloid heterogeneity than previously reported from both published and our generated scRNA-Seq. In human HNC patients, we identified four macrophage, seven monocyte, and five dendritic cell subsets. One macrophage subset expressed high levels of CXCL9 in both humans and our murine HNC model. Interestingly, recent studies identified CXCL9 as a biomarker for immunotherapy response. Observing conserved Cxcl9-expressing macrophages between murine and human HNC could help assess treatment response in cancer patients. We will assess the translational value of murine HPV+ and HPV- HNC models compared to their human counterparts. To this end, we are generating scRNA-Seq data from different HNC murine models: MOC1, MOC2, mEERL. The MOC2 and MOC1 cell lines recapitulate a responsive and non-responsive HPV- HNC tumor phenotype, respectively. The mEERL cell line recapitulates an HPV+ HNC phenotype. Our analysis shows high concordance between mouse models with immunological cold TIME (MOC2) and a proportion of human HNCs. We also identified changes in myeloid landscapes between immunological cold-to-hot TIME mice receiving ICI (anti-PD1 + anti-CTL4) therapies. We will present an in-depth comparative analysis to provide mechanistic insights into observed differences in the TIME between mice and humans. We will also present a novel bioinformatics analysis to interrogate how cell-cell signaling interactions between myeloid and T cell subsets influence their downstream gene regulatory network and tumor immune phenotypes. This analysis provides a reference for future studies in HNC immunology.

Citation Format: Athena E. Golfinos, Wei Wang, Aisha Mergaert, Paul F. Lambert, Huy Q. Dinh. Defining tumor infiltrating myeloid microenvironment landscapes in HPV+ and HPV- head and neck cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1994.

CD33 Expression on Peripheral Blood Monocytes Predicts Efficacy of Anti-PD-1 Immunotherapy Against Non-Small Cell Lung Cancer

Non-small cell lung carcinoma (NSCLC) is the leading cause of cancer-related deaths globally. Immune checkpoint blockade (ICB) has transformed cancer medicine, with anti-programmed cell death protein 1 (anti-PD-1) therapy now well-utilized for treating NSCLC. Still, not all patients with NSCLC respond positively to anti-PD-1 therapy, and some patients acquire resistance to treatment. There remains an urgent need to find markers predictive of anti-PD-1 responsiveness. To this end, we performed mass cytometry on peripheral blood mononuclear cells from 26 patients with NSCLC during anti-PD-1 treatment. Patients who responded to anti-PD-1 ICB displayed significantly higher levels of antigen-presenting myeloid cells, including CD9+ nonclassical monocytes, and CD33hi classical monocytes. Using matched pre-post treatment samples, we found that the baseline pre-treatment frequencies of CD33hi monocytes predicted patient responsiveness to anti-PD-1 therapy. Moreover, some of these classical and nonclassical monocyte subsets were associated with reduced immunosuppression by T regulatory (CD4+FOXP3+CD25+) cells in the same patients. Our use of machine learning corroborated the association of specific monocyte markers with responsiveness to ICB. Our work provides a high-dimensional profile of monocytes in NSCLC and links CD33 expression on monocytes with anti-PD-1 effectiveness in patients with NSCLC.

Integrated single-cell transcriptome analysis reveals heterogeneity of esophageal squamous cell carcinoma microenvironment

The tumor microenvironment is a highly complex ecosystem of diverse cell types, which shape cancer biology and impact the responsiveness to therapy. Here, we analyze the microenvironment of esophageal squamous cell carcinoma (ESCC) using single-cell transcriptome sequencing in 62,161 cells from blood, adjacent nonmalignant and matched tumor samples from 11 ESCC patients. We uncover heterogeneity in most cell types of the ESCC stroma, particularly in the fibroblast and immune cell compartments. We identify a tumor-specific subset of CST1⁺ myofibroblasts with prognostic values and potential biological significance. CST1⁺ myofibroblasts are also highly tumor-specific in other cancer types. Additionally, a subset of antigen-presenting fibroblasts is revealed and validated. Analyses of myeloid and T lymphoid lineages highlight the immunosuppressive nature of the ESCC microenvironment, and identify cancer-specific expression of immune checkpoint inhibitors. This work establishes a rich resource of stromal cell types of the ESCC microenvironment for further understanding of ESCC biology.

Isoform-level quantification for single-cell RNA sequencing

MOTIVATION

RNA expression at isoform level is biologically more informative than at gene level and can potentially reveal cellular subsets and corresponding biomarkers that are not visible at gene level. However, due to the strong 3' bias sequencing protocol, mRNA quantification for high-throughput single-cell RNA sequencing such as Chromium Single Cell 3' 10× Genomics is currently performed at the gene level.

RESULTS

We have developed an isoform-level quantification method for high-throughput single-cell RNA sequencing by exploiting the concepts of transcription clusters and isoform paralogs. The method, called Scasa, compares well in simulations against competing approaches including Alevin, Cellranger, Kallisto, Salmon, Terminus and STARsolo at both isoform- and gene-level expression. The reanalysis of a CITE-Seq dataset with isoform-based Scasa reveals a subgroup of CD14 monocytes missed by gene-based methods.

AVAILABILITY AND IMPLEMENTATION

Implementation of Scasa including source code, documentation, tutorials and test data supporting this study is available at Github: https://github.com/eudoraleer/scasa and Zenodo: https://doi.org/10.5281/zenodo.5712503.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Atherosclerosis Impairs Naive CD4 T-Cell Responses via Disruption of Glycolysis

Objective

CD4 T cells are important regulators of atherosclerotic progression. The metabolic profile of CD4 T cells controls their signaling and function, but how atherosclerosis affects T-cell metabolism is unknown. Here, we sought to determine the impact of atherosclerosis on CD4 T-cell metabolism and the contribution of such metabolic alterations to atheroprogression.

Approach and Results

Using PCR arrays, we profiled the expression of metabolism genes in CD4 T cells from atherosclerotic apolipoprotein-E knockout mice fed a Western diet. These cells exhibited dysregulated expression of genes critically involved in glycolysis and fatty acid degradation, compared with those from animals fed a standard laboratory diet. We examined how T-cell metabolism was changed in either Western diet–fed apolipoprotein-E knockout mice or samples from patients with cardiovascular disease by measuring glucose uptake, activation, and proliferation in CD4 T cells. We found that naive CD4 T cells from Western diet–fed apolipoprotein-E knockout mice failed to uptake glucose and displayed impaired proliferation and activation, compared with CD4 T cells from standard laboratory diet–fed animals. Similarly, we observed that naive CD4 T-cell frequencies were reduced in the circulation of human subjects with high cardiovascular disease compared with low cardiovascular disease. Naive T cells from high cardiovascular disease subjects also showed reduced proliferative capacity.

Conclusions

These results highlight the dysfunction that occurs in CD4 T-cell metabolism and immune responses during atherosclerosis. Targeting metabolic pathways within naive CD4 T cells could thus yield novel therapeutic approaches for improving CD4 T-cell responses against atheroprogression.

Coexpression of CD71 and CD117 Identifies an Early Unipotent Neutrophil Progenitor Population in Human Bone Marrow

Neutrophils are the most abundant peripheral immune cells and thus, are continually replenished by bone marrow-derived progenitors. Still, how newly identified neutrophil subsets fit into the bone marrow neutrophil lineage remains unclear. Here, we use mass cytometry to show that two recently defined human neutrophil progenitor populations contain a homogeneous progenitor subset we term "early neutrophil progenitors" (eNePs) (Lin^-CD66b⁺CD117⁺CD71⁺). Surface marker- and RNA-expression analyses, together with in vitro colony formation and in vivo adoptive humanized mouse transfers, indicate that eNePs are the earliest human neutrophil progenitors. Furthermore, we identified CD71 as a marker associated with the earliest neutrophil developmental stages. Expression of CD71 marks proliferating neutrophils, which were expanded in the blood of melanoma patients and detectable in blood and tumors from lung cancer patients. In summary, we establish CD117⁺CD71⁺ eNeP as the inceptive human neutrophil progenitor and propose a refined model of the neutrophil developmental lineage in bone marrow.

Single-cell transcriptomics identifies gene expression networks driving differentiation and tumorigenesis in the human fallopian tube

The human fallopian tube harbors the cell of origin for the majority of high-grade serous "ovarian" cancers (HGSCs), but its cellular composition, particularly the epithelial component, is poorly characterized. We perform single-cell transcriptomic profiling of around 53,000 individual cells from 12 primary fallopian specimens to map their major cell types. We identify 10 epithelial subpopulations with diverse transcriptional programs. Based on transcriptional signatures, we reconstruct a trajectory whereby secretory cells differentiate into ciliated cells via a RUNX3^high intermediate. Computational deconvolution of advanced HGSCs identifies the "early secretory" population as a likely precursor state for the majority of HGSCs. Its signature comprises both epithelial and mesenchymal features and is enriched in mesenchymal-type HGSCs (p = 6.7 × 10^-27), a group known to have particularly poor prognoses. This cellular and molecular compendium of the human fallopian tube in cancer-free women is expected to advance our understanding of the earliest stages of fallopian epithelial neoplasia.

Naive CD8+ T Cells Expressing CD95 Increase Human Cardiovascular Disease Severity

OBJECTIVE

Cardiovascular disease (CVD) remains a significant global health concern with a high degree of mortality. While CD4⁺ T cells have been extensively studied in CVD, the importance of CD8⁺ T cells in this disease, despite their abundance and increased activation in human atherosclerotic plaques, remains largely unknown. Thus, the objective of this study was to compare peripheral T-cell signatures between humans with a high (severe) risk of CVD (including myocardial infarction or stroke) and those with a low risk of CVD. Approach and Results: Using mass cytometry, we uncovered a naive CD8⁺ T (T_N) cell population expressing CD95 (termed CD95⁺CD8⁺ stem cell memory T [CD8 T_SCM] cells) that was enriched in patients with high compared with low CVD. This T-cell subset enrichment within individuals with high CVD was a relative increase and resulted from the loss of CD95^lo cells within the T_N compartment. We found that CD8 T_SCM cells positively correlated with CVD risk in humans, while CD8⁺ T_N cells were inversely correlated. Atherosclerotic apolipoprotein E-deficient (ApoE^-/-) mice also displayed respective 7- and 2-fold increases in CD8⁺ T_SCM frequencies within the peripheral blood and aorta-draining paraaortic lymph nodes compared with C57BL/6J mice. CD8⁺ T_SCM cells were 1.7-fold increased in aortas from western diet fed ApoE^-/- mice compared with normal laboratory diet-fed ApoE^-/- mice. Importantly, transfer of T_SCM cells into immune-deficient Rag.Ldlr recipient mice that lacked T cells increased atherosclerosis, illustrating the importance of these cells in atherogenesis.

CONCLUSIONS

CD8⁺ T_SCM cells are increased in humans with high CVD. As these T_SCM cells promote atherosclerosis, targeting them may attenuate atherosclerotic plaque progression.

Meta-Analysis of Leukocyte Diversity in Atherosclerotic Mouse Aortas

The diverse leukocyte infiltrate in atherosclerotic mouse aortas was recently analyzed in 9 single-cell RNA sequencing and 2 mass cytometry studies. In a comprehensive meta-analysis, we confirm 4 known macrophage subsets-resident, inflammatory, interferon-inducible cell, and Trem2 (triggering receptor expressed on myeloid cells-2) foamy macrophages-and identify a new macrophage subset resembling cavity macrophages. We also find that monocytes, neutrophils, dendritic cells, natural killer cells, innate lymphoid cells-2, and CD (cluster of differentiation)-8 T cells form prominent and separate immune cell populations in atherosclerotic aortas. Many CD4 T cells express IL (interleukin)-17 and the chemokine receptor CXCR (C-X-C chemokine receptor)-6. A small number of regulatory T cells and T helper 1 cells is also identified. Immature and naive T cells are present in both healthy and atherosclerotic aortas. Our meta-analysis overcomes limitations of individual studies that, because of their experimental approach, over- or underrepresent certain cell populations. Mass cytometry studies demonstrate that cell surface phenotype provides valuable information beyond the cell transcriptomes. The present analysis helps resolve some long-standing controversies in the field. First, Trem2+ foamy macrophages are not proinflammatory but interferon-inducible cell and inflammatory macrophages are. Second, about half of all foam cells are smooth muscle cell-derived, retaining smooth muscle cell transcripts rather than transdifferentiating to macrophages. Third, Pf4, which had been considered specific for platelets and megakaryocytes, is also prominently expressed in the main population of resident vascular macrophages. Fourth, a new type of resident macrophage shares transcripts with cavity macrophages. Finally, the discovery of a prominent innate lymphoid cell-2 cluster links the single-cell RNA sequencing work to recent flow cytometry data suggesting a strong atheroprotective role of innate lymphoid cells-2. This resolves apparent discrepancies regarding the role of T helper 2 cells in atherosclerosis based on studies that predated the discovery of innate lymphoid cells-2 cells.

Frontline Science: Kindlin-3 is essential for patrolling and phagocytosis functions of nonclassical monocytes during metastatic cancer surveillance

Nonclassical monocytes maintain vascular homeostasis by patrolling the vascular endothelium, responding to inflammatory signals, and scavenging cellular debris. Nonclassical monocytes also prevent metastatic tumor cells from seeding new tissues, but whether the patrolling function of nonclassical monocytes is required for this process is unknown. To answer this question, we utilized an inducible-knockout mouse that exhibits loss of the integrin-adaptor protein Kindlin-3 specifically in nonclassical monocytes. We show that Kindlin-3-deficient nonclassical monocytes are unable to patrol the vascular endothelium in either the lungs or periphery. We also find that Kindlin-3-deficient nonclassical monocytes cannot firmly adhere to, and instead "slip" along, the vascular endothelium. Loss of patrolling activity by nonclassical monocytes was phenocopied by ablation of LFA-1, an integrin-binding partner of Kindlin-3. When B16F10 murine melanoma tumor cells were introduced into Kindlin-3-deficient mice, nonclassical monocytes showed defective patrolling towards tumor cells and failure to ingest tumor particles in vivo. Consequently, we observed a significant, 4-fold increase in lung tumor metastases in mice possessing Kindlin-3-deficient nonclassical monocytes. Thus, we conclude that the patrolling function of nonclassical monocytes is mediated by Kindlin-3 and essential for these cells to maintain vascular endothelial homeostasis and prevent tumor metastasis to the lung.

Macrophage Syk-PI3Kγ Inhibits Antitumor Immunity: SRX3207, a Novel Dual Syk-PI3K Inhibitory Chemotype Relieves Tumor Immunosuppression

Macrophages (MΦ) play a critical role in tumor growth, immunosuppression, and inhibition of adaptive immune responses in cancer. Hence, targeting signaling pathways in MΦs that promote tumor immunosuppression will provide therapeutic benefit. PI3Kγ has been recently established by our group and others as a novel immuno-oncology target. Herein, we report that an MΦ Syk-PI3K axis drives polarization of immunosuppressive MΦs that establish an immunosuppressive tumor microenvironment in in vivo syngeneic tumor models. Genetic or pharmacologic inhibition of Syk and/or PI3Kγ in MΦs promotes a proinflammatory MΦ phenotype, restores CD8+ T-cell activity, destabilizes HIF under hypoxia, and stimulates an antitumor immune response. Assay for transposase-accessible Chromatin using Sequencing (ATAC-seq) analyses on the bone marrow-derived macrophages (BMDM) show that inhibition of Syk kinase promotes activation and binding of NF-κB motif in SykMC-KO BMDMs, thus stimulating immunostimulatory transcriptional programming in MΦs to suppress tumor growth. Finally, we have developed in silico the "first-in-class" dual Syk/PI3K inhibitor, SRX3207, for the combinatorial inhibition of Syk and PI3K in one small molecule. This chemotype demonstrates efficacy in multiple tumor models and represents a novel combinatorial approach to activate antitumor immunity.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Identification of an Early Unipotent Neutrophil Progenitor with Pro-tumoral Activity in Mouse and Human Bone Marrow

Neutrophils are short-lived cells that play important roles in both health and disease. Neutrophils and monocytes originate from the granulocyte monocyte progenitor (GMP) in bone marrow; however, unipotent neutrophil progenitors are not well defined. Here, we use cytometry by time of flight (CyTOF) and single-cell RNA sequencing (scRNA-seq) methodologies to identify a committed unipotent early-stage neutrophil progenitor (NeP) in adult mouse bone marrow. Importantly, we found a similar unipotent NeP (hNeP) in human bone marrow. Both NeP and hNeP generate only neutrophils. NeP and hNeP both significantly increase tumor growth when transferred into murine cancer models, including a humanized mouse model. hNeP are present in the blood of treatment-naive melanoma patients but not of healthy subjects. hNeP can be readily identified by flow cytometry and could be used as a biomarker for early cancer discovery. Understanding the biology of hNeP should allow the development of new therapeutic targets for neutrophil-related diseases, including cancer.

Human Monocyte Heterogeneity as Revealed by High-Dimensional Mass Cytometry

Objective- Three distinct human monocyte subsets have been identified based on the surface marker expression of CD14 and CD16. We hypothesized that monocytes were likely more heterogeneous in composition. Approach and Results- We used the high dimensionality of mass cytometry together with the FlowSOM clustering algorithm to accurately identify and define monocyte subsets in blood of healthy human subjects and those with coronary artery disease (CAD). To study the behavior and functionality of the newly defined monocyte subsets, we performed RNA sequencing, transwell migration, and efferocytosis assays. Here, we identify 8 human monocyte subsets based on their surface marker phenotype. We found that 3 of these subsets fall within the CD16⁺ nonclassical monocyte population and 4 subsets belong to the CD14⁺ classical monocytes, illustrating significant monocyte heterogeneity in humans. As nonclassical monocytes are important in modulating atherosclerosis in mice, we studied the functions of our 3 newly identified nonclassical monocytes in subjects with CAD. We found a marked expansion of a Slan⁺CXCR6⁺ nonclassical monocyte subset in CAD subjects, which was positively correlated with CAD severity. This nonclassical subset can migrate towards CXCL16 and shows an increased efferocytosis capacity, indicating it may play an atheroprotective role. Conclusions- Our data demonstrate that human nonclassical monocytes are a heterogeneous population, existing of several subsets with functional differences. These subsets have changed frequencies in the setting of severe CAD. Understanding how these newly identified subsets modulate CAD will be important for CAD-based therapies that target myeloid cells.

Preparation of Whole Bone Marrow for Mass Cytometry Analysis of Neutrophil-lineage Cells

In this article, we present a protocol that is optimized to preserve neutrophil-lineage cells in fresh BM for whole BM CyTOF analysis. We utilized a myeloid-biased 39-antibody CyTOF panel to evaluate the hematopoietic system with a focus on the neutrophil-lineage cells by using this protocol. The CyTOF result was analyzed with an open-resource dimensional reduction algorithm, viSNE, and the data was presented to demonstrate the outcome of this protocol. We have discovered new neutrophil-lineage cell populations based on this protocol. This protocol of fresh whole BM preparation may be used for 1), CyTOF analysis to discover unidentified cell populations from whole BM, 2), investigating whole BM defects for patients with blood disorders such as leukemia, 3), assisting optimization of fluorescence-activated flow cytometry protocols that utilize fresh whole BM.

Monocyte heterogeneity and functions in cancer

Monocytes are innate immune cells of the mononuclear phagocyte system that have emerged as important regulators of cancer development and progression. Our understanding of monocytes has advanced from viewing these cells as a homogenous population to a heterogeneous system of cells that display diverse responses to different stimuli. During cancer, different monocyte subsets perform functions that contribute to both pro‐ and antitumoral immunity, including phagocytosis, secretion of tumoricidal mediators, promotion of angiogenesis, remodeling of the extracellular matrix, recruitment of lymphocytes, and differentiation into tumor‐associated macrophages and dendritic cells. The ability of cancer to evade immune recognition and clearance requires protumoral signals to outweigh ongoing attempts by the host immune system to prevent tumor growth. This review discusses current understanding of monocyte heterogeneity during homeostasis, highlights monocyte functions in cancer progression, and describes monocyte‐targeted therapeutic strategies for cancer treatment.

Identification of distinct mutational patterns and new driver genes in oesophageal squamous cell carcinomas and adenocarcinomas

OBJECTIVES

Oesophageal squamous cell carcinoma (OSCC) and adenocarcinoma (OAC) are distinct cancers in terms of a number of clinical and epidemiological characteristics, complicating the design of clinical trials and biomarker developments. We analysed 1048 oesophageal tumour-germline pairs from both subtypes, to characterise their genomic features, and biological and clinical significance.

DESIGN

Previously exome-sequenced samples were re-analysed to identify significantly mutated genes (SMGs) and mutational signatures. The biological functions of novel SMGs were investigated using cell line and xenograft models. We further performed whole-genome bisulfite sequencing and chromatin immunoprecipitation (ChIP)-seq to characterise epigenetic alterations.

RESULTS

OSCC and OAC displayed nearly mutually exclusive sets of driver genes, indicating that they follow independent developmental paths. The combined sample size allowed the statistical identification of a number of novel subtype-specific SMGs, mutational signatures and prognostic biomarkers. Particularly, we identified a novel mutational signature similar to Catalogue Of Somatic Mutations In Cancer (COSMIC)signature 16, which has prognostic value in OSCC. Two newly discovered SMGs, CUL3 and ZFP36L2, were validated as important tumour-suppressors specific to the OSCC subtype. We further identified their additional loss-of-function mechanisms. CUL3 was homozygously deleted specifically in OSCC and other squamous cell cancers (SCCs). Notably, ZFP36L2 is associated with super-enhancer in healthy oesophageal mucosa; DNA hypermethylation in its super-enhancer reduced active histone markers in squamous cancer cells, suggesting an epigenetic inactivation of a super-enhancer-associated SCC suppressor.

CONCLUSIONS

These data comprehensively contrast differences between OSCC and OAC at both genomic and epigenomic levels, and reveal novel molecular features for further delineating the pathophysiological mechanisms and treatment strategies for these cancers.

DNA methylation loss in late-replicating domains is linked to mitotic cell division

DNA methylation loss occurs frequently in cancer genomes, primarily within lamina-associated, late-replicating regions termed Partially Methylated Domains (PMDs). We profiled 39 diverse primary tumors and 8 matched adjacent tissues using Whole-Genome Bisulfite Sequencing (WGBS), and analyzed them alongside 343 additional human and 206 mouse WGBS datasets. We identified a local CpG sequence context associated with preferential hypomethylation in PMDs. Analysis of CpGs in this context (“Solo-WCGWs”) revealed previously undetected PMD hypomethylation in almost all healthy tissue types. PMD hypomethylation increased with age, beginning during fetal development, and appeared to track the accumulation of cell divisions. In cancer, PMD hypomethylation depth correlated with somatic mutation density and cell-cycle gene expression, consistent with its reflection of mitotic history, and suggesting its application as a mitotic clock. We propose that late replication leads to lifelong progressive methylation loss, which acts as a biomarker for cellular aging and which may contribute to oncogenesis.

Atlas of the Immune Cell Repertoire in Mouse Atherosclerosis Defined by Single-Cell RNA-Sequencing and Mass Cytometry

RATIONALE

Atherosclerosis is a chronic inflammatory disease that is driven by the interplay of pro- and anti-inflammatory leukocytes in the aorta. Yet, the phenotypic and transcriptional diversity of aortic leukocytes is poorly understood.

OBJECTIVE

We characterized leukocytes from healthy and atherosclerotic mouse aortas in-depth by single-cell RNA-sequencing and mass cytometry (cytometry by time of flight) to define an atlas of the immune cell landscape in atherosclerosis.

METHODS AND RESULTS

Using single-cell RNA-sequencing of aortic leukocytes from chow diet- and Western diet-fed Apoe^-/- and Ldlr^-/- mice, we detected 11 principal leukocyte clusters with distinct phenotypic and spatial characteristics while the cellular repertoire in healthy aortas was less diverse. Gene set enrichment analysis on the single-cell level established that multiple pathways, such as for lipid metabolism, proliferation, and cytokine secretion, were confined to particular leukocyte clusters. Leukocyte populations were differentially regulated in atherosclerotic Apoe^-/- and Ldlr^-/- mice. We confirmed the phenotypic diversity of these clusters with a novel mass cytometry 35-marker panel with metal-labeled antibodies and conventional flow cytometry. Cell populations retrieved by these protein-based approaches were highly correlated to transcriptionally defined clusters. In an integrated screening strategy of single-cell RNA-sequencing, mass cytometry, and fluorescence-activated cell sorting, we detected 3 principal B-cell subsets with alterations in surface markers, functional pathways, and in vitro cytokine secretion. Leukocyte cluster gene signatures revealed leukocyte frequencies in 126 human plaques by a genetic deconvolution strategy. This approach revealed that human carotid plaques and microdissected mouse plaques were mostly populated by macrophages, T-cells, and monocytes. In addition, the frequency of genetically defined leukocyte populations in carotid plaques predicted cardiovascular events in patients.

CONCLUSIONS

The definition of leukocyte diversity by high-dimensional analyses enables a fine-grained analysis of aortic leukocyte subsets, reveals new immunologic mechanisms and cell-type-specific pathways, and establishes a functional relevance for lesional leukocytes in human atherosclerosis.

Genomic and Epigenomic Heterogeneity of Hepatocellular Carcinoma

Understanding the intratumoral heterogeneity of hepatocellular carcinoma is instructive for developing personalized therapy and identifying molecular biomarkers. Here we applied whole-exome sequencing to 69 samples from 11 patients to resolve the genetic architecture of subclonal diversification. Spatial genomic diversity was found in all 11 hepatocellular carcinoma cases, with 29% of driver mutations being heterogeneous, including TERT, ARID1A, NOTCH2, and STAG2. Similar with other cancer types, TP53 mutations were always shared between all tumor regions, that is, located on the "trunk" of the evolutionary tree. In addition, we found that variants within several drug targets such as KIT, SYK, and PIK3CA were mutated in a fully clonal manner, indicating their therapeutic potentials for hepatocellular carcinoma. Temporal dissection of mutational signatures suggested that mutagenic processes associated with exposure to aristolochic acid and aflatoxin might play a more important role in early, as opposed to late, stages of hepatocellular carcinoma development. Moreover, we observed extensive intratumoral epigenetic heterogeneity in hepatocellular carcinoma based on multiple independent analytical methods and showed that intratumoral methylation heterogeneity might play important roles in the biology of hepatocellular carcinoma cells. Our results also demonstrated prominent heterogeneity of intratumoral methylation even in a stable hepatocellular carcinoma genome. Together, these findings highlight widespread intratumoral heterogeneity at both the genomic and epigenomic levels in hepatocellular carcinoma and provide an important molecular foundation for better understanding the pathogenesis of this malignancy. Cancer Res; 77(9); 2255-65. ©2017 AACR.

The RdDM pathway is required for basal heat tolerance in Arabidopsis

Heat stress affects epigenetic gene silencing in Arabidopsis. To test for a mechanistic involvement of epigenetic regulation in heat-stress responses, we analyzed the heat tolerance of mutants defective in DNA methylation, histone modifications, chromatin-remodeling, or siRNA-based silencing pathways. Plants deficient in NRPD2, the common second-largest subunit of RNA polymerases IV and V, and in the Rpd3-type histone deacetylase HDA6 were hypersensitive to heat exposure. Microarray analysis demonstrated that NRPD2 and HDA6 have independent roles in transcriptional reprogramming in response to temperature stress. The misexpression of protein-coding genes in nrpd2 mutants recovering from heat correlated with defective epigenetic regulation of adjacent transposon remnants which involved the loss of control of heat-stress-induced read-through transcription. We provide evidence that the transcriptional response to temperature stress, at least partially, relies on the integrity of the RNA-dependent DNA methylation pathway.

Advanced methylome analysis after bisulfite deep sequencing: an example in Arabidopsis

Deep sequencing after bisulfite conversion (BS-Seq) is the method of choice to generate whole genome maps of cytosine methylation at single base-pair resolution. Its application to genomic DNA of Arabidopsis flower bud tissue resulted in the first complete methylome, determining a methylation rate of 6.7% in this tissue. BS-Seq reads were mapped onto an in silico converted reference genome, applying the so-called 3-letter genome method. Here, we present BiSS (Bisufite Sequencing Scorer), a new method applying Smith-Waterman alignment to map bisulfite-converted reads to a reference genome. In addition, we introduce a comprehensive adaptive error estimate that accounts for sequencing errors, erroneous bisulfite conversion and also wrongly mapped reads. The re-analysis of the Arabidopsis methylome data with BiSS mapped substantially more reads to the genome. As a result, it determines the methylation status of an extra 10% of cytosines and estimates the methylation rate to be 7.7%. We validated the results by individual traditional bisulfite sequencing for selected genomic regions. In addition to predicting the methylation status of each cytosine, BiSS also provides an estimate of the methylation degree at each genomic site. Thus, BiSS explores BS-Seq data more extensively and provides more information for downstream analysis.

Uncovering cis-regulatory sequence requirements for context-specific transcription factor binding

The regulation of gene expression is mediated at the transcriptional level by enhancer regions that are bound by sequence-specific transcription factors (TFs). Recent studies have shown that the in vivo binding sites of single TFs differ between developmental or cellular contexts. How this context-specific binding is encoded in the cis-regulatory DNA sequence has, however, remained unclear. We computationally dissect context-specific TF binding sites in Drosophila, Caenorhabditis elegans, mouse, and human and find distinct combinations of sequence motifs for partner factors, which are predictive and reveal specific motif requirements of individual binding sites. We predict that TF binding in the early Drosophila embryo depends on motifs for the early zygotic TFs Vielfaltig (also known as Zelda) and Tramtrack. We validate experimentally that the activity of Twist-bound enhancers and Twist binding itself depend on Vielfaltig motifs, suggesting that Vielfaltig is more generally important for early transcription. Our finding that the motif content can predict context-specific binding and that the predictions work across different Drosophila species suggests that characteristic motif combinations are shared between sites, revealing context-specific motif codes (cis-regulatory signatures), which appear to be conserved during evolution. Taken together, this study establishes a novel approach to derive predictive cis-regulatory motif requirements for individual TF binding sites and enhancers. Importantly, the method is generally applicable across different cell types and organisms to elucidate cis-regulatory sequence determinants and the corresponding trans-acting factors from the increasing number of tissue- and cell-type-specific TF binding studies.

Alternative splicing and nonsense-mediated decay modulate expression of important regulatory genes in Arabidopsis

Alternative splicing (AS) coupled to nonsense-mediated decay (NMD) is a post-transcriptional mechanism for regulating gene expression. We have used a high-resolution AS RT-PCR panel to identify endogenous AS isoforms which increase in abundance when NMD is impaired in the Arabidopsis NMD factor mutants, upf1-5 and upf3-1. Of 270 AS genes (950 transcripts) on the panel, 102 transcripts from 97 genes (32%) were identified as NMD targets. Extrapolating from these data around 13% of intron-containing genes in the Arabidopsis genome are potentially regulated by AS/NMD. This cohort of naturally occurring NMD-sensitive AS transcripts also allowed the analysis of the signals for NMD in plants. We show the importance of AS in introns in 5' or 3'UTRs in modulating NMD-sensitivity of mRNA transcripts. In particular, we identified upstream open reading frames overlapping the main start codon as a new trigger for NMD in plants and determined that NMD is induced if 3'-UTRs were >350 nt. Unexpectedly, although many intron retention transcripts possess NMD features, they are not sensitive to NMD. Finally, we have shown that AS/NMD regulates the abundance of transcripts of many genes important for plant development and adaptation including transcription factors, RNA processing factors and stress response genes.

Genetic rearrangements can modify chromatin features at epialleles

Analogous to genetically distinct alleles, epialleles represent heritable states of different gene expression from sequence-identical genes. Alleles and epialleles both contribute to phenotypic heterogeneity. While alleles originate from mutation and recombination, the source of epialleles is less well understood. We analyze active and inactive epialleles that were found at a transgenic insert with a selectable marker gene in Arabidopsis. Both converse expression states are stably transmitted to progeny. The silent epiallele was previously shown to change its state upon loss-of-function of trans-acting regulators and drug treatments. We analyzed the composition of the epialleles, their chromatin features, their nuclear localization, transcripts, and homologous small RNA. After mutagenesis by T-DNA transformation of plants carrying the silent epiallele, we found new active alleles. These switches were associated with different, larger or smaller, and non-overlapping deletions or rearrangements in the 3' regions of the epiallele. These cis-mutations caused different degrees of gene expression stability depending on the nature of the sequence alteration, the consequences for transcription and transcripts, and the resulting chromatin organization upstream. This illustrates a tight dependence of epigenetic regulation on local structures and indicates that sequence alterations can cause epigenetic changes at some distance in regions not directly affected by the mutation. Similar effects may also be involved in gene expression and chromatin changes in the vicinity of transposon insertions or excisions, recombination events, or DNA repair processes and could contribute to the origin of new epialleles.

Pleiotropic effects of the wheat dehydrin DHN-5 on stress responses in Arabidopsis

We have previously reported that transgenic Arabidopsis plants overexpressing the wheat dehydrin DHN-5 show enhanced tolerance to osmotic stresses. In order to understand the mechanisms through which DHN-5 exerts this effect, we performed transcriptome profiling using the Affymetrix ATH1 microarray. Our data show an altered expression of 77 genes involved mainly in transcriptional regulation, cellular metabolism, stress tolerance and signaling. Among the up-regulated genes, we identified those which are known to be stress-related genes. Several late embryogenesis abundant (LEA) genes, ABA/stress-related genes (such as RD29B) and those involved in pathogen responses (PR genes) are among the most up-regulated genes. In addition, the MDHAR gene involved in the ascorbate biosynthetic pathway was also up-regulated. This up-regulation was correlated with higher ascorbate content in two dehydrin transgenic lines. In agreement with this result and as ascorbate is known to be an antioxidant, we found that both transgenic lines show enhanced tolerance to oxidative stress caused by H₂O₂. On the other hand, multiple types of transcription factors constitute the largest group of the down-regulated genes. Moreover, three members of the jasmonate-ZIM domain (JAZ) proteins which are negative regulators of jasmonate signaling were severely down-regulated. Interestingly, the dehydrin-overexpressing lines exhibit less sensitivity to jasmonate than wild-type plants and changes in regulation of jasmonate-responsive genes, in a manner similar to that in the jasmonate-insensitive jai3-1 mutant. Altogether, our data unravel the potential pleiotropic effects of DHN-5 on both abiotic and biotic stress responses in Arabidopsis.

Epigenetic regulation of repetitive elements is attenuated by prolonged heat stress in Arabidopsis

Epigenetic factors determine responses to internal and external stimuli in eukaryotic organisms. Whether and how environmental conditions feed back to the epigenetic landscape is more a matter of suggestion than of substantiation. Plants are suitable organisms with which to address this question due to their sessile lifestyle and diversification of epigenetic regulators. We show that several repetitive elements of Arabidopsis thaliana that are under epigenetic regulation by transcriptional gene silencing at ambient temperatures and upon short term heat exposure become activated by prolonged heat stress. Activation can occur without loss of DNA methylation and with only minor changes to histone modifications but is accompanied by loss of nucleosomes and by heterochromatin decondensation. Whereas decondensation persists, nucleosome loading and transcriptional silencing are restored upon recovery from heat stress but are delayed in mutants with impaired chromatin assembly functions. The results provide evidence that environmental conditions can override epigenetic regulation, at least transiently, which might open a window for more permanent epigenetic changes.

Cooperation of multiple chromatin modifications can generate unanticipated stability of epigenetic States in Arabidopsis

Epigenetic changes of gene expression can potentially be reversed by developmental programs, genetic manipulation, or pharmacological interference. However, a case of transcriptional gene silencing, originally observed in tetraploid Arabidopsis thaliana plants, created an epiallele resistant to many mutations or inhibitor treatments that activate many other suppressed genes. This raised the question about the molecular basis of this extreme stability. A combination of forward and reverse genetics and drug application provides evidence for an epigenetic double lock that is only alleviated upon the simultaneous removal of both DNA methylation and histone methylation. Therefore, the cooperation of multiple chromatin modifications can generate unanticipated stability of epigenetic states and contributes to heritable diversity of gene expression patterns.