Immune responses are different between individuals and personal health histories and unique environmental conditions should collectively determine the present state of immune cells. However, the molecular systems underlying such heterogeneity remain elusive. Here, we conducted a systematic time-lapse single-cell analysis, using 171 single-cell libraries and 30 mass cytometry datasets intensively for seven healthy individuals. We found substantial diversity in immune-cell profiles between different individuals. These patterns showed daily fluctuations even within the same individual. Similar diversities were also observed for the T-cell and B-cell receptor repertoires. Detailed immune-cell profiles at healthy statuses should give essential background information to understand their immune responses, when the individual is exposed to various environmental conditions. To demonstrate this idea, we conducted the similar analysis for the same individuals on the vaccination of influenza and SARS-CoV-2. In fact, we detected distinct responses to vaccines between individuals, although key responses are common. Single-cell immune-cell profile data should make fundamental data resource to understand variable immune responses, which are unique to each individual.

Targeted sequencing for somatic mutations across the hotspots of 50 cancer-related genes was performed using biopsy specimens to investigate whether clinicopathological factors and genomic alterations correlated with prognosis in locally advanced cervical cancer. Seventy patients diagnosed with International Federation of Obstetrics and Gynecology (FIGO) stage III to IVA cervical cancer underwent radiotherapy or concurrent chemoradiotherapy at the National Cancer Center Hospital between January 2008 and December 2017. Mutations were detected in 47 of 70 [67% of cases; frequency of genetic alterations was as follows: PIK3CA (51%), FBXW7 (10%), PTEN (7.1%), and TP53 (5.7%)]. The Cancer Genome Atlas (TCGA) datasets showed a similar distribution of somatic mutations, but PIK3CA mutation frequency was significantly higher in our cohort than in TCGA datasets (P = 0.028). Patients with TP53 mutation were significantly related to poor progression-free survival (PFS) (hazard ratio [HR] = 3.53, P = 0.042). Patients with tumor diameters > 70 mm were associated with poor prognosis (HR = 2.96, P = 0.0048). Patients with non-HPV16/18 genotypes had worse prognosis than those with HPV16/18 genotypes (HR = 2.15, P = 0.030). Hence, patients with locally advanced cervical cancer, TP53 mutation, large tumor diameter, and non-HPV16/18 genotype were independently correlated with poor PFS, despite concurrent chemoradiotherapy.

Activation of the KRAS oncogene is a source of replication stress, but how this stress is generated and how it is tolerated by cancer cells remain poorly understood. Here we show that induction of KRASG12V expression in untransformed cells triggers H3K27me3 and HP1-associated chromatin compaction in an RNA transcription dependent manner, resulting in replication fork slowing and cell death. Furthermore, elevated ATR expression is necessary and sufficient for tolerance of KRASG12V-induced replication stress to expand replication stress-tolerant cells (RSTCs). PrimPol is phosphorylated at Ser255, a potential Chk1 substrate site, under KRASG12V-induced replication stress and promotes repriming to maintain fork progression and cell survival in an ATR/Chk1-dependent manner. However, ssDNA gaps are generated at heterochromatin by PrimPol-dependent repriming, leading to genomic instability. These results reveal a role of ATR-PrimPol in enabling precancerous cells to survive KRAS-induced replication stress and expand clonally with accumulation of genomic instability.

Radiotherapy (RT) plus immunotherapy is a promising modality; however, the therapeutic effects are insufficient, and the molecular mechanism requires clarification to further develop combination therapies. Here, we found that the RNA virus sensor pathway dominantly regulates the cellular immune response in NSCLC and ESCC cell lines. Notably, transposable elements (TEs), especially long terminal repeats (LTRs), functioned as key ligands for the RNA virus sensor RIG-I, and the mTOR-LTR-RIG-I axis induced the cellular immune response and dendritic cell and macrophage infiltration after irradiation. Moreover, RIG-I-dependent immune activation was observed in ESCC patient tissue. scRNA sequencing and spatial transcriptome analysis revealed that radiotherapy induced the expression of LTRs, and the RNA virus sensor pathway in immune and cancer cells; this pathway was also found to mediate tumour conversion to an immunological hot state. Here, we report the upstream and ligand of the RNA virus sensor pathway functions in irradiated cancer tissues.

Radiotherapy (RT) combined with immunotherapy is promising; however, the immune response signature in the clinical setting after RT remains unclear. Here, by integrative spatial and single-cell analyses using multiplex immunostaining (CODEX), spatial transcriptome (VISIUM), and single-cell RNA sequencing, we substantiated the infiltration of immune cells into tumors with dynamic changes in immunostimulatory and immunosuppressive gene expression after RT. In addition, our comprehensive analysis uncovered time- and cell type-dependent alterations in the gene expression profile after RT. Furthermore, myeloid cells showed prominent up-regulation of immune response-associated genes after RT. Notably, a subset of infiltrating tumor-associated myeloid cells showing PD-L1 positivity exhibited significant up-regulation of immunostimulatory (HMGB1 and ISG15), immunosuppressive (SIRPA and IDO1), and protumor genes (CXCL8, CCL3, IL-6, and IL-1AB), which can be targets of immunotherapy in combination with PD-L1. These datasets will provide information on the RT-induced gene signature to seek an appropriate target for personalized immunotherapy combined with RT and guide the timing of combination therapy.

Small cell lung cancer (SCLC) is the most aggressive type of lung cancer with high mortality. One of the MYC family genes, MYC, MYCL or MYCN, is amplified in ~20% of the SCLCs; therefore, MYC proteins are potential therapeutic targets in SCLC patients. We investigated the therapeutic impact of Omomyc, a MYC dominant negative, in a panel of SCLC cell lines. Strikingly, Omomyc suppressed the growth of all tested cell lines by inducing cell cycle arrest and/or apoptosis. Induction of G1 arrest by Omomyc was found to be dependent on the activation of CDKN1A, in part, through the TP73 pathway. Our results strongly indicate that SCLC cells carrying amplification of MYC, MYCL or MYCN are addicted to MYC function, suggesting that MYC targeting would be an efficient therapeutic option for SCLC patients.

We aimed to elucidate the effect of JQ1, a BET inhibitor, on small cell lung cancers (SCLCs) with MYCL amplification and/or expression. Fourteen SCLC cell lines, including four with MYCL amplification, were examined for the effects of JQ1 on protein and gene expression by Western blot and mRNA microarray analyses. The sensitivity of SCLC cells to JQ1 was assessed by cell growth and apoptosis assays. MYCL was expressed in all the 14 cell lines, whereas MYC/MYCN expression was restricted mostly to cell lines with gene amplification. ASCL1, a transcription factor shown to play a role in SCLC, was also expressed in 11/14 cell lines. All SCLC cell lines were sensitive to JQ1 with GI50 values ≤1.23 μM, with six of them showing GI50 values <0.1 μM. Expression of MYCL as well as MYCN, ASCL1 and other driver oncogenes including CDK6 was reduced by JQ1 treatment, in particular in the cell lines with high expression of the respective genes; however, no association was observed between the sensitivity to JQ1 and the levels of MYCL, MYCN and ASCL1 expression. In contrast, levels of CDK6 expression and its reduction rates by JQ1 were associated with JQ1 sensitivity. Therefore, we concluded that CDK6 is a novel target of JQ1 and predictive marker for JQ1 sensitivity in SCLC cells.

Thyroid hormone (TH) is important for tissue repair because it regulates cellular differentiation. Intrahepatic TH activity is controlled by both serum TH levels and hepatic deiodinases. TH substrate (T4) is converted into active hormone (T3) by deiodinase 1 (D1) but into inactive hormone (rT3) by deiodinase 3 (D3). Although the relative expressions of D1 and D3 are known to change during liver injury, the cell types and signaling mechanisms involved are unclear. We evaluated the hypothesis that changes in hepatic deiodinases result from repair-related activation of the Hedgehog pathway in stromal cells. We localized deiodinase expression, assessed changes during injury, and determined how targeted manipulation of Hedgehog signaling in stromal cells impacted hepatic deiodinase expression, TH content, and TH action in rodents. Humans with chronic liver disease were also studied. In healthy liver, hepatocytes strongly expressed D1 and stromal cells weakly expressed D3. During injury, hepatocyte expression of D1 decreased, whereas stromal expression of D3 increased, particularly in myofibroblasts. Conditionally disrupting Hedgehog signaling in myofibroblasts normalized deiodinase expression. Repair-related changes in deiodinases were accompanied by reduced hepatic TH content and TH-regulated gene expression. In patients, this was reflected by increased serum rT3. Moreover, the decreases in the free T3 to rT3 and free T4 to rT3 ratios distinguished advanced from mild fibrosis, even in individuals with similar serum levels of TSH and free T4. In conclusion, the Hedgehog-dependent changes in liver stromal cells drive repair-related changes in hepatic deiodinase expression that promote intrahepatic hypothyroidism, thereby limiting exposure to T3, an important factor for cellular differentiation.

Dysfunction of the p53 network is a major cause of cancer development, and selective elimination of p53-inactivated cancer cells therefore represents an ideal therapeutic strategy. In this study, we performed a microRNA target screen that identified NEK9 (NIMA-related kinase 9) as a crucial regulator of cell-cycle progression in p53-inactivated cancer cells. NEK9 depletion selectively inhibited proliferation in p53-deficient cancer cells both in vitro and in vivo. The resultant cell-cycle arrest occurred predominantly in G1 phase, and exhibited senescence-like features. Furthermore, NEK9 repression affected expression of a broad range of genes encoding cell-cycle regulators and factors involved in mRNA processing, suggesting a novel role for NEK9 in p53-deficient cells. Lung adenocarcinoma patients with positive staining for NEK9 and mutant p53 proteins exhibited significantly poorer prognoses, suggesting that expression of both proteins promotes tumor growth. Our findings demonstrate that a novel NEK9 network regulates the growth of cancer cells lacking functional p53.

Here we conducted an integrative multi-omics analysis to understand how cancers harbor various types of aberrations at the genomic, epigenomic and transcriptional levels. In order to elucidate biological relevance of the aberrations and their mutual relations, we performed whole-genome sequencing, RNA-Seq, bisulfite sequencing and ChIP-Seq of 26 lung adenocarcinoma cell lines. The collected multi-omics data allowed us to associate an average of 536 coding mutations and 13,573 mutations in promoter or enhancer regions with aberrant transcriptional regulations. We detected the 385 splice site mutations and 552 chromosomal rearrangements, representative cases of which were validated to cause aberrant transcripts. Averages of 61, 217, 3687 and 3112 mutations are located in the regulatory regions which showed differential DNA methylation, H3K4me3, H3K4me1 and H3K27ac marks, respectively. We detected distinct patterns of aberrations in transcriptional regulations depending on genes. We found that the irregular histone marks were characteristic to EGFR and CDKN1A, while a large genomic deletion and hyper-DNA methylation were most frequent for CDKN2A. We also used the multi-omics data to classify the cell lines regarding their hallmarks of carcinogenesis. Our datasets should provide a valuable foundation for biological interpretations of interlaced genomic and epigenomic aberrations.

The serine threonine kinase Akt is a core survival factor that underlies a variety of human diseases. Although regulatory phosphorylation and dephosphorylation have been well documented, the other posttranslational mechanisms that modulate Akt activity remain unclear. We show here that tetratricopeptide repeat domain 3 (TTC3) is an E3 ligase that interacts with Akt. TTC3 contains a canonical RING finger motif, a pair of tetratricopeptide motifs, a putative Akt phosphorylation site, and nuclear localization signals, and is encoded by a gene within the Down syndrome (DS) critical region on chromosome 21. TTC3 is an Akt-specific E3 ligase that binds to phosphorylated Akt and facilitates its ubiquitination and degradation within the nucleus. Moreover, DS cells exhibit elevated TTC3 expression, reduced phosphorylated Akt, and accumulation in the G(2)M phase, which can be reversed by TTC3 siRNA or Myr-Akt. Thus, interaction between TTC3 and Akt may contribute to the clinical symptoms of DS.

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal malignancies. The genomic landscape of the PDAC genome features four frequently mutated genes (KRAS, CDKN2A, TP53, and SMAD4) and dozens of candidate driver genes altered at low frequency, including potential clinical targets. Circulating cell-free DNA (cfDNA) is a promising resource to detect and monitor molecular characteristics of tumors. In the present study, we determined the mutational status of KRAS in plasma cfDNA using multiplex picoliter-droplet digital PCR in 259 patients with PDAC. We constructed a novel modified SureSelect-KAPA-Illumina platform and an original panel of 60 genes. We then performed targeted deep sequencing of cfDNA and matched germline DNA samples in 48 patients who had ≥1% mutant allele frequencies of KRAS in plasma cfDNA. Importantly, potentially targetable somatic mutations were identified in 14 of 48 patients (29.2%) examined by targeted deep sequencing of cfDNA. We also analyzed somatic copy number alterations based on the targeted sequencing data using our in-house algorithm, and potentially targetable amplifications were detected. Assessment of mutations and copy number alterations in plasma cfDNA may provide a prognostic and diagnostic tool to assist decisions regarding optimal therapeutic strategies for PDAC patients.

To identify stage I lung adenocarcinoma patients with a poor prognosis who will benefit from adjuvant therapy.

Mice deficient for the fibulin-5 gene (Fbln5(-/-)) develop pelvic organ prolapse (POP) due to compromised elastic fibers and upregulation of matrix metalloprotease (MMP)-9. Here, we used casein zymography, inhibitor profiling, affinity pull-down, and mass spectrometry to discover additional protease upregulated in the vaginal wall of Fbln5(-/-) mice, herein named V1 (25 kDa). V1 was a serine protease with trypsin-like activity similar to protease, serine (PRSS) 3, a major extrapancreatic trypsinogen, was optimum at pH 8.0, and predominantly detected in estrogenized vaginal epithelium of Fbln5(-/-) mice. PRSS3 was (a) localized in epithelial secretions, (b) detected in media of vaginal organ culture from both Fbln5(-/-) and wild type mice, and (c) cleaved fibulin-5 in vitro. Expression of two serine protease inhibitors [Serpina1a (α1-antitrypsin) and Elafin] was dysregulated in Fbln5(-/-) epithelium. Finally, we confirmed that PRSS3 was expressed in human vaginal epithelium and that SERPINA1 and Elafin were downregulated in vaginal tissues from women with POP. These data collectively suggest that the balance between proteases and their inhibitors contributes to support of the pelvic organs in humans and mice.

ARID1A encodes an SWI/SNF chromatin-remodeling factor and is frequently mutated in various cancers. This study demonstrates that ARID1A-deficient cancer cells are specifically vulnerable to inhibition of the antioxidant glutathione (GSH) and the glutamate-cysteine ligase synthetase catalytic subunit (GCLC), a rate-limiting enzyme for GSH synthesis. Inhibition of GCLC markedly decreased GSH in ARID1A-deficient cancer cells, leading to apoptotic cell death triggered by excessive amounts of reactive oxygen species. The vulnerability of ARID1A-deficient cancer cells results from low basal levels of GSH due to impaired expression of SLC7A11. The SLC7A11-encoded cystine transporter supplies cells with cysteine, a key source of GSH, and its expression is enhanced by ARID1A-mediated chromatin remodeling. Thus, ARID1A-deficient cancers are susceptible to synthetic lethal targeting of GCLC.

Advanced cancer genomics technologies are now being employed in clinical sequencing, where next-generation sequencers are used to simultaneously identify multiple types of DNA alterations for prescription of molecularly targeted drugs. However, no computational tool is available to accurately detect DNA alterations in formalin-fixed paraffin-embedded (FFPE) samples commonly used in hospitals. Here, we developed a computational tool tailored to the detection of single nucleotide variations, indels, fusions, and copy number alterations in FFPE samples. Elaborated multilayer noise filters reduced the inherent noise while maintaining high sensitivity, as evaluated in tumor-unmatched normal samples using orthogonal technologies. This tool, cisCall, should facilitate clinical sequencing in everyday diagnostics. It is available at https://www.ciscall.org .

Histone acetylation at DNA double-strand break (DSB) sites by CBP and p300 histone acetyltransferases (HATs) is critical for the recruitment of DSB repair proteins to chromatin. Here, we show that CBP and p300 HATs also function in DSB repair by transcriptionally activating the BRCA1 and RAD51 genes, which are involved in homologous recombination (HR), a major DSB repair system. siRNA-mediated depletion of CBP and p300 impaired HR activity and downregulated BRCA1 and RAD51 at the protein and mRNA levels. Chromatin immunoprecipitation assays showed that CBP and p300 bind to the promoter regions of the BRCA1 and RAD51 genes, and that depletion of CBP and/or p300 reduces H3 and H4 acetylation and inhibits binding of the transcription factor E2F1 to these promoters. Depletion of CBP and p300 impaired DNA damage-induced phosphorylation and chromatin binding of the single-strand DNA-binding protein RPA following BRCA1-mediated DNA end resection. Consistent with this, subsequent phosphorylation of CHK1 and activation of the G2/M damage checkpoint were also impaired. These results indicate that the HATs CBP and p300 play multiple roles in the activation of the cellular response to DSBs.

We analyzed whole-exome sequencing data from 97 Japanese lung adenocarcinoma patients and identified several putative cancer-related genes and pathways. Particularly, we observed that cancer-related mutation patterns were significantly different between different ethnic groups. As previously reported, mutations in the EGFR gene were characteristic to Japanese, while those in the KRAS gene were more frequent in Caucasians. Furthermore, during the course of this analysis, we found that cancer-specific somatic mutations can be detected without sequencing normal tissue counterparts. 64% of the germline variants could be excluded using a total of 217 external Japanese exome datasets. We also show that a similar approach may be used for other three ethnic groups, although the discriminative power depends on the ethnic group. We demonstrate that the ATM gene and the PAPPA2 gene could be identified as cancer prognosis related genes. By bypassing the sequencing of normal tissue counterparts, this approach provides a useful means of not only reducing the time and cost of sequencing but also analyzing archive samples, for which normal tissue counterparts are not available.

Mortality amongst patients hospitalised for heart failure (HHF) in Western and Asian countries may differ, but this has not been investigated using individual patient-level data (IPLD). We sought to remedy this through rigorous statistical analysis of HHF registries and variable selection from a systematic literature review.

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. Our study revealed that the scRNA-seq analysis could cover a broader range of immune cells of TME in the biopsy-resected small samples of tumors, detecting even small subgroups of B cells or Treg cells in the tumors, although CyTOF could distinguish the specific populations in more depth. These findings demonstrate that scRNA-seq analysis is a highly-feasible platform for elucidating the complexity of TME in small biopsy tumors, which would provide a novel strategies to overcome a therapeutic difficulties against cancer heterogeneity in TME.