PNPO-PLP axis senses prolonged hypoxia in macrophages by regulating lysosomal activity.

Oxygen is critical for all metazoan organisms on the earth and impacts various biological processes in physiological and pathological conditions. While oxygen-sensing systems inducing acute hypoxic responses, including the hypoxia-inducible factor pathway, have been identified, those operating in prolonged hypoxia remain to be elucidated. Here we show that pyridoxine 5'-phosphate oxidase (PNPO), which catalyses bioactivation of vitamin B6, serves as an oxygen sensor and regulates lysosomal activity in macrophages. Decreased PNPO activity under prolonged hypoxia reduced an active form of vitamin B6, pyridoxal 5'-phosphate (PLP), and inhibited lysosomal acidification, which in macrophages led to iron dysregulation, TET2 protein loss and delayed resolution of the inflammatory response. Among PLP-dependent metabolism, supersulfide synthesis was suppressed in prolonged hypoxia, resulting in the lysosomal inhibition and consequent proinflammatory phenotypes of macrophages. The PNPO-PLP axis creates a distinct layer of oxygen sensing that gradually shuts down PLP-dependent metabolism in response to prolonged oxygen deprivation.

Whole blood transcriptome analysis for age- and gender-specific gene expression profiling in Japanese individuals.

Whole blood transcriptome analysis is a valuable approachin medical research, primarily due to the ease of sample collection and the richness of the information obtained. Since the expression profile of individual genes in the analysis is influenced by medical traits and demographic attributes such as age and gender, there has been a growing demand for a comprehensive database for blood transcriptome analysis. Here, we performed whole blood RNA sequencing (RNA-seq) analysis on 576 participants stratified by age (20-30s and 60-70s) and gender from cohorts of the Tohoku Medical Megabank (TMM). A part of female segment included pregnant women. We did not exclude the globin gene family in our RNA-seq study, which enabled us to identify instances of hereditary persistence of fetal hemoglobin based on the HBG1 and HBG2 expression information. Comparing stratified populations allowed us to identify groups of genes associated with age-related changes and gender differences. We also found that the immune response status, particularly measured by neutrophil-to-lymphocyte ratio (NLR), strongly influences the diversity of individual gene expression profiles in whole blood transcriptome analysis. This stratification has resulted in a data set that will be highly beneficial for future whole blood transcriptome analysis in the Japanese population.

C4S DB: Comprehensive Collection and Comparison for ChIP-Seq Database.

Combining multiple binding profiles, such as transcription factors and histone modifications, is a crucial step in revealing the functions of complex biological systems. Although a massive amount of chromatin immunoprecipitation followed by sequencing (ChIP-seq) data is available, existing ChIP-seq databases or repositories focus on individual experiments, and it is difficult to elucidate orchestrated regulation by DNA-binding elements. We developed the Comprehensive Collection and Comparison for ChIP-Seq Database (C4S DB) to provide researchers with insights into the combination of DNA binding elements based on quality-assessed public ChIP-seq data. The C4S DB is based on >16,000 human ChIP-seq experiments and provides two main web interfaces to discover the relationships between ChIP-seq data. "Gene browser" illustrates the landscape of distributions of binding elements around a specified gene, and "global similarity," a hierarchical clustering heatmap based on a similarity between two ChIP-seq experiments, gives an overview of genome-wide relations of regulatory elements. These functions promote the identification or evaluation of both gene-specific and genome-wide colocalization or mutually exclusive localization. Modern web technologies allow users to search for and aggregate large-scale experimental data through interactive web interfaces with quick responses. The C4S DB is available at https://c4s.site.

Deficiency of CHAMP1, a gene related to intellectual disability, causes impaired neuronal development and a mild behavioural phenotype.

CHAMP1 is a gene associated with intellectual disability, which was originally identified as being involved in the maintenance of kinetochore-microtubule attachment. To explore the neuronal defects caused by CHAMP1 deficiency, we established mice that lack CHAMP1. Mice that are homozygous knockout for CHAMP1 were slightly smaller than wild-type mice and died soon after birth on pure C57BL/6J background. Although gross anatomical defects were not found in CHAMP1 -/- mouse brains, mitotic cells were increased in the cerebral cortex. Neuronal differentiation was delayed in CHAMP1 -/- neural stem cells in vitro, which was also suggested in vivo by CHAMP1 knockdown. In a behavioural test battery, adult CHAMP1 heterozygous knockout mice showed mild memory defects, altered social interaction, and depression-like behaviours. In transcriptomic analysis, genes related to neurotransmitter transport and neurodevelopmental disorder were downregulated in embryonic CHAMP1 -/- brains. These results suggest that CHAMP1 plays a role in neuronal development, and CHAMP1-deficient mice resemble some aspects of individuals with CHAMP1 mutations.

High levels of chromosomal instability facilitate the tumor growth and sphere formation.

Most cancer cells show chromosomal instability (CIN), a condition in which chromosome missegregation occurs at high rates. Growing evidence suggests that CIN is not just a consequence of, but a driving force for, oncogenic transformation, although the relationship between CIN and tumorigenesis has not been fully elucidated. Here we found that conventional two-dimensional (2D) culture of HeLa cells, a cervical cancer-derived cell line, was a heterogenous population containing cells with different CIN levels. Although cells with high-CIN levels (high-CIN cells) grew more slowly compared with cells with low-CIN levels (low-CIN cells) in 2D monolayer culture, they formed tumors in nude mice and larger spheres in three-dimensional (3D) culture, which was more representative of the in vivo environment. The duration of mitosis was longer in high-CIN cells, reflecting their higher mitotic defects. Single-cell genome sequencing revealed that high-CIN cells exhibited a higher karyotype heterogeneity compared with low-CIN cells. Intriguingly, the karyotype heterogeneity was reduced in the spheres formed by high-CIN cells, suggesting that cells with growth advantages were selected, although genomic copy number changes specific for spheres were not identified. When we examined gene expression profiles, genes related to the K-ras signaling were upregulated, while those related to the unfolded protein response were downregulated in high-CIN cells in 3D culture compared with 2D culture, suggesting the relevance of these genes for their survival. Our data suggested that, although CIN is disadvantageous in monolayer culture, it promotes the selection of cells with growth advantages under in vivo environments, which may lead to tumorigenesis.

CEBPB is required for NRF2-mediated drug resistance in NRF2-activated non-small cell lung cancer cells.

NRF2 is a transcription activator that plays a key role in cytoprotection against oxidative stress. Although increased NRF2 activity is principally beneficial for our health, NRF2 activation in cancer cells is detrimental, as it drives their malignant progression. We previously found that CCAAT/enhancer-binding protein B (CEBPB) cooperates with NRF2 in NRF2-activated lung cancer and enhances tumour-initiating activity by promoting NOTCH3 expression. However, the general contribution of CEBPB in lung cancer is rather controversial, probably because the role of CEBPB depends on cooperating transcription factors in each cellular context. To understand how NRF2 shapes the function of CEBPB in NRF2-activated lung cancers and its biological consequence, we comprehensively explored NRF2-CEBPB-coregulated genes and found that genes involved in drug metabolism and detoxification were characteristically enriched. Indeed, CEBPB and NRF2 cooperatively contribute to the drug resistance. We also found that CEBPB is directly regulated by NRF2, which is likely to be advantageous for the coexpression and cooperative function of NRF2 and CEBPB. These results suggest that drug resistance of NRF2-activated lung cancers is achieved by the cooperative function of NRF2 and CEBPB.

Enhancer remodeling promotes tumor-initiating activity in NRF2-activated non-small cell lung cancers.

Transcriptional dysregulation, which can be caused by genetic and epigenetic alterations, is a fundamental feature of many cancers. A key cytoprotective transcriptional activator, NRF2, is often aberrantly activated in non-small cell lung cancers (NSCLCs) and supports both aggressive tumorigenesis and therapeutic resistance. Herein, we find that persistently activated NRF2 in NSCLCs generates enhancers at gene loci that are not normally regulated by transiently activated NRF2 under physiological conditions. Elevated accumulation of CEBPB in NRF2-activated NSCLCs is found to be one of the prerequisites for establishment of the unique NRF2-dependent enhancers, among which the NOTCH3 enhancer is shown to be critical for promotion of tumor-initiating activity. Enhancer remodeling mediated by NRF2-CEBPB cooperativity promotes tumor-initiating activity and drives malignancy of NRF2-activated NSCLCs via establishment of the NRF2-NOTCH3 regulatory axis.

Landscape of electrophilic and inflammatory stress-mediated gene regulation in human lymphoblastoid cell lines.

Human lymphoblastoid cell lines (LCLs) are valuable for the functional analyses of diseases. We have established more than 4200 LCLs as one of the resources of an integrated biobank. While oxidative and inflammatory stresses play critical roles in the onset and progression of various diseases, the responsiveness of LCLs, especially that of biobank-made LCLs, to these stresses has not been established. To address how LCLs respond to these stresses, in this study, we performed RNA sequencing of eleven human LCLs that were treated with an electrophile, diethyl maleate (DEM) and/or an inflammatory mediator, lipopolysaccharide (LPS). We found that over two thousand genes, including those regulated by a master regulator of the electrophilic/oxidative stress response, NRF2, were upregulated in LCLs treated with DEM, while approximately three hundred genes, including inflammation-related genes, were upregulated in LPS-treated LCLs. Of the LPS-induced genes, a subset of proinflammatory genes was repressed by DEM, supporting the notion that DEM suppresses the expression of proinflammatory genes through NRF2 activation. Conversely, a part of DEM-induced gene was repressed by LPS, suggesting reciprocal interference between electrophilic and inflammatory stress-mediated pathways. These data clearly demonstrate that LCLs maintain, by and large, responsive pathways against oxidative and inflammatory stresses and further endorse the usefulness of the LCL supply from the biobank.

Theoretical characterisation of strand cross-correlation in ChIP-seq.

BACKGROUND: Strand cross-correlation profiles are used for both peak calling pre-analysis and quality control (QC) in chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis. Despite its potential for robust and accurate assessments of signal-to-noise ratio (S/N) because of its peak calling independence, it remains unclear what aspects of quality such strand cross-correlation profiles actually measure. RESULTS: We introduced a simple model to simulate the mapped read-density of ChIP-seq and then derived the theoretical maximum and minimum of cross-correlation coefficients between strands. The results suggest that the maximum coefficient of typical ChIP-seq samples is directly proportional to the number of total mapped reads and the square of the ratio of signal reads, and inversely proportional to the number of peaks and the length of read-enriched regions. Simulation analysis supported our results and evaluation using 790 ChIP-seq data obtained from the public database demonstrated high consistency between calculated cross-correlation coefficients and estimated coefficients based on the theoretical relations and peak calling results. In addition, we found that the mappability-bias-correction improved sensitivity, enabling differentiation of maximum coefficients from the noise level. Based on these insights, we proposed virtual S/N (VSN), a novel peak call-free metric for S/N assessment. We also developed PyMaSC, a tool to calculate strand cross-correlation and VSN efficiently. VSN achieved most consistent S/N estimation for various ChIP targets and sequencing read depths. Furthermore, we demonstrated that a combination of VSN and pre-existing peak calling results enable the estimation of the numbers of detectable peaks for posterior experiments and assess peak calling results. CONCLUSIONS: We present the first theoretical insights into the strand cross-correlation, and the results reveal the potential and the limitations of strand cross-correlation analysis. Our quality assessment framework using VSN provides peak call-independent QC and will help in the evaluation of peak call analysis in ChIP-seq experiments.