BATF epigenetically and transcriptionally controls the activation program of regulatory T cells in human tumors.

Regulatory T (T) cells suppress effective antitumor immunity in tumor-bearing hosts, thereby becoming promising targets in cancer immunotherapy. Despite the importance of T cells in tumor immunity, little is known about their differentiation process and epigenetic profiles in the tumor microenvironment (TME). Here, we showed that T cells in the TME of human lung cancers harbored a completely different open chromatin profile compared with CD8 T cells, conventional CD4 T cells in the TME, and peripheral T cells.

Regulatory T-cell development in the tumor microenvironment.

Regulatory T (Treg) cells are required for maintaining self-tolerance and preventing the development of autoimmune diseases. However, Treg cells are abundant in tumors and suppress antitumor immunity, contributing to tumor development and growth. Thus, the selective deletion of tumor-infiltrating Treg cells is important for successful Treg cell-targeted therapies, providing effective antitumor immunity without inducing deleterious autoimmune disorders.

Bronchoalveolar lavage fluid reveals factors contributing to the efficacy of PD-1 blockade in lung cancer.

Bronchoalveolar lavage is commonly performed to assess inflammation and identify responsible pathogens in lung diseases. Findings from bronchoalveolar lavage might be used to evaluate the immune profile of the lung tumor microenvironment (TME). To investigate whether bronchoalveolar lavage fluid (BALF) analysis can help identify patients with non-small cell lung cancer (NSCLC) who respond to immune checkpoint inhibitors (ICIs), BALF and blood were prospectively collected before initiating nivolumab.

PD-1 blockade therapy promotes infiltration of tumor-attacking exhausted T cell clonotypes.

PD-1 blockade exerts clinical efficacy against various types of cancer by reinvigorating T cells that directly attack tumor cells (tumor-specific T cells) in the tumor microenvironment (TME), and tumor-infiltrating lymphocytes (TILs) also comprise nonspecific bystander T cells. Here, using single-cell sequencing, we show that TILs include skewed T cell clonotypes, which are characterized by exhaustion (T) or nonexhaustion signatures (T). Among skewed clonotypes, those in the T, but not those in the T, cluster respond to autologous tumor cell lines.

Lactic acid promotes PD-1 expression in regulatory T cells in highly glycolytic tumor microenvironments.

The balance of programmed death-1 (PD-1)-expressing CD8 T cells and regulatory T (Treg) cells in the tumor microenvironment (TME) determines the clinical efficacy of PD-1 blockade therapy through the competition of their reactivation. However, factors that determine this balance remain unknown. Here, we show that Treg cells gain higher PD-1 expression than effector T cells in highly glycolytic tumors, including MYC-amplified tumors and liver tumors.

Depletion of central memory CD8 T cells might impede the antitumor therapeutic effect of Mogamulizumab.

Regulatory T (Treg) cells are important negative regulators of immune homeostasis, but in cancers they tone down the anti-tumor immune response. They are distinguished by high expression levels of the chemokine receptor CCR4, hence their targeting by the anti-CCR4 monoclonal antibody mogamulizumab holds therapeutic promise. Here we show that despite a significant reduction in peripheral effector Treg cells, clinical responses are minimal in a cohort of patients with advanced CCR4-negative solid cancer in a phase Ib study (NCT01929486).

Importance of lymph node immune responses in MSI-H/dMMR colorectal cancer.

Patients with colorectal cancers (CRCs) generally exhibit improved survival through intensive lymph node (LN) dissection. However, recent progress in cancer immunotherapy revisits the potential importance of regional LNs, where T cells are primed to attack tumor cells. To elucidate the role of regional LN, we investigated the immunological status of nonmetastatic regional LN lymphocytes (LNLs) in comparison with those of the tumor microenvironment (tumor-infiltrating lymphocytes; TILs) using flow cytometry and next-generation sequencing.

Potentiality of multiple modalities for single-cell analyses to evaluate the tumor microenvironment in clinical specimens.

Single-cell level analysis is powerful tool to assess the heterogeneity of cellular components in tumor microenvironments (TME). In this study, we investigated immune-profiles using the single-cell analyses of endoscopically- or surgically-resected tumors, and peripheral blood mononuclear cells from gastric cancer patients. Furthermore, we technically characterized two distinct platforms of the single-cell analysis; RNA-seq-based analysis (scRNA-seq), and mass cytometry-based analysis (CyTOF), both of which are broadly embraced technologies.

The PD-1 expression balance between effector and regulatory T cells predicts the clinical efficacy of PD-1 blockade therapies.

Immune checkpoint blockade has provided a paradigm shift in cancer therapy, but the success of this approach is very variable; therefore, biomarkers predictive of clinical efficacy are urgently required. Here, we show that the frequency of PD-1CD8 T cells relative to that of PD-1 regulatory T (T) cells in the tumor microenvironment can predict the clinical efficacy of programmed cell death protein 1 (PD-1) blockade therapies and is superior to other predictors, including PD ligand 1 (PD-L1) expression or tumor mutational burden. PD-1 expression by CD8 T cells and T cells negatively impacts effector and immunosuppressive functions, respectively.

Immune Suppression by PD-L2 against Spontaneous and Treatment-Related Antitumor Immunity.

Purpose: To evaluate the detailed immunosuppressive role(s) of PD-L2 given that its detailed role(s) remains unclear in PD-1 signal blockade therapy in animal models and humans.

Experimental Design: We generated mouse cell lines harboring various status of PD-L1/PD-L2 and evaluated the tumor growth and phenotypes of tumor-infiltrated lymphocytes using several PD-1 signal blockades in animal models. In humans, the correlation between immune-related gene expression and (encoding PD-L1) or (encoding PD-L2) was investigated using The Cancer Genome Atlas (TCGA) datasets.

A new computational method to predict transcriptional activity of a DNA sequence from diverse datasets of massively parallel reporter assays.

In recent years, the dramatic increase in the number of applications for massively parallel reporter assay (MPRA) technology has produced a large body of data for various purposes. However, a computational model that can be applied to decipher regulatory codes for diverse MPRAs does not exist yet. Here, we propose a new computational method to predict the transcriptional activity of MPRAs, as well as luciferase reporter assays, based on the TRANScription FACtor database.

Analysis of RNA decay factor mediated RNA stability contributions on RNA abundance.

Background: Histone epigenome data determined by chromatin immunoprecipitation sequencing (ChIP-seq) is used in identifying transcript regions and estimating expression levels. However, this estimation does not always correlate with eventual RNA expression levels measured by RNA sequencing (RNA-seq). Part of the inconsistency may arise from the variance in RNA stability, where the transcripts that are more or less abundant than predicted RNA expression from histone epigenome data are inferred to be more or less stable.

BRIC-seq: a genome-wide approach for determining RNA stability in mammalian cells.

We recently developed a novel transcriptome analysis method, termed 5'-bromo-uridine (BrU) immunoprecipitation chase-deep sequencing analysis (BRIC-seq). BRIC-seq enables the determination of genome-wide RNA stability by chasing chronological decreases of BrU-labeled RNAs under physiologically undisturbed conditions. The RNA half-life of each transcript is calculated from the decreasing number of BrU-labeled RNA sequence tags measured by deep sequencing of BrU-labeled RNAs.

Identification of hundreds of novel UPF1 target transcripts by direct determination of whole transcriptome stability.

UPF1 eliminates aberrant mRNAs harboring premature termination codons, and regulates the steady-state levels of normal physiological mRNAs. Although genome-wide studies of UPF1 targets performed, previous studies did not distinguish indirect UPF1 targets because they could not determine UPF1-dependent altered RNA stabilities. Here, we measured the decay rates of the whole transcriptome in UPF1-depleted HeLa cells using BRIC-seq, an inhibitor-free method for directly measuring RNA stability.

Linear regression models predicting strength of transcriptional activity of promoters.

We developed linear regression models which predict strength of transcriptional activity of promoters from their sequences. Intrinsic transcriptional strength data of 451 human promoter sequences in three cell lines (HEK293, MCF7 and 3T3), which were measured by systematic luciferase reporter gene assays, were used to build the models. The models sum up contributions of CG dinucleotide content and transcription factor binding sites (TFBSs) to transcriptional strength.

Predicting promoter activities of primary human DNA sequences.

We developed a computer program that can predict the intrinsic promoter activities of primary human DNA sequences. We observed promoter activity using a quantitative luciferase assay and generated a prediction model using multiple linear regression. Our program achieved a prediction accuracy correlation coefficient of 0.

Massive transcriptional start site analysis of human genes in hypoxia cells.

Combining our full-length cDNA method and the massively parallel sequencing technology, we developed a simple method to collect precise positional information of transcriptional start sites (TSSs) together with digital information of the gene-expression levels in a high throughput manner. We applied this method to observe gene-expression changes in a colon cancer cell line cultured in normoxic and hypoxic conditions. We generated more than 100 million 36-base TSS-tag sequences and revealed comprehensive features of hypoxia responsive alterations in the transcriptional landscape of the human genome.

Low conservation and species-specific evolution of alternative splicing in humans and mice: comparative genomics analysis using well-annotated full-length cDNAs.

Using full-length cDNA sequences, we compared alternative splicing (AS) in humans and mice. The alignment of the human and mouse genomes showed that 86% of 199 426 total exons in human AS variants were conserved in the mouse genome. Of the 20 392 total human AS variants, however, 59% consisted of all conserved exons.

Distinct class of putative "non-conserved" promoters in humans: comparative studies of alternative promoters of human and mouse genes.

Although recent studies have revealed that the majority of human genes are subject to regulation of alternative promoters, the biological relevance of this phenomenon remains unclear. We have also demonstrated that roughly half of the human RefSeq genes examined contain putative alternative promoters (PAPs). Here we report large-scale comparative studies of PAPs between human and mouse counterpart genes.

Intrinsic promoter activities of primary DNA sequences in the human genome.

In order to understand an overview of promoter activities intrinsic to primary DNA sequences in the human genome within a particular cell type, we carried out systematic quantitative luciferase assays of DNA fragments corresponding to putative promoters for 472 human genes which are expressed in HEK (human embryonic kidney epithelial) 293 cells. We observed the promoter activities of them were distributed in a bimodal manner; putative promoters belonging to the first group (with strong promoter activities) were designated as P1 and the latter (with weak promoter activities) as P2. The frequencies of the TATA-boxes, the CpG islands, and the overall G + C-contents were significantly different between these two populations, indicating there are two separate groups of promoters.

Diverse DNA methylation statuses at alternative promoters of human genes in various tissues.

We characterized the DNA methylation status at 144 tissue-biased and 37 non-tissue-biased alternative promoters of 61 human genes in five normal tissues. Analysis of the collected data revealed that (i) DNA methylation status differed greatly among alternative promoters belonging to the same gene; (ii) DNA methylation status differed between tissues for the majority of the individual promoters, and (iii) 80-90% of CpG-island-containing promoters were not methylated on either allele throughout the five tissues examined. Furthermore, although the statistical significance was not as clear as for the above features, we also found that (iv) the DNA methylation patterns of tissue-biased promoters changed more drastically than those of non-tissue-biased promoters; (v) tissue-biased promoters tended to be less methylated than their respective alternative promoters in the tissues where they were preferentially expressed, and (vi) the 'null' methylation pattern of a given promoter was enriched in the tissues where the transcription was most active.