An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases.

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

Progressive, multi-organ, and multi-layered nature of cancer cachexia.

Cancer cachexia is a complex, multifaceted condition that negatively impacts the health, treatment efficacy, and economic status of cancer patients. The management of cancer cachexia is an essential clinical need. Cancer cachexia is currently defined mainly according to the severity of weight loss and sarcopenia (i.e., macrosymptoms). However, such macrosymptoms may be insufficient to give clinicians clues on how to manage this condition as these symptoms appear at the late stage of cancer. We need to understand earlier events during the progression of cancer cachexia so as not to miss a clinical opportunity to control this complex syndrome. Recent research indicates that cancer-induced changes in the host are much wider than previously recognized, including disruption of liver function and the immune system. Furthermore, such changes are observed before the occurrence of visible distant metastases (i.e., in early, localized cancers). In light of these findings, we propose to expand the definition of cancer cachexia to include all cancer-induced changes to host physiology, including changes caused by early, localized cancers. This new definition of cancer cachexia can provide a new perspective on this topic, which can stimulate the research and development of novel cancer cachexia therapies.

Time-series blood cytokine profiles correlate with treatment responses in triple-negative breast cancer patients.

BACKGROUND: Triple-negative breast cancer (TNBC) is the most heterogeneous breast cancer subtype. Partly due to its heterogeneity, it is currently challenging to stratify TNBC patients and predict treatment outcomes. METHODS: In this study, we examined blood cytokine profiles of TNBC patients throughout treatments (pre-treatment, during chemotherapy, pre-surgery, and 1 year after the surgery in a total of 294 samples). We analyzed the obtained cytokine datasets using weighted correlation network analyses, protein-protein interaction analyses, and logistic regression analyses. RESULTS: We identified five cytokines that correlate with good clinical outcomes: interleukin (IL)-1α, TNF-related apoptosis-inducing ligand (TRAIL), Stem Cell Factor (SCF), Chemokine ligand 5 (CCL5 also known as RANTES), and IL-16. The expression of these cytokines was decreased during chemotherapy and then restored after the treatment. Importantly, patients with good clinical outcomes had constitutively high expression of these cytokines during treatments. Protein-protein interaction analyses implicated that these five cytokines promote an immune response. Logistic regression analyses revealed that IL-1α and TRAIL expression levels at pre-treatment could predict treatment outcomes in our cohort. CONCLUSION: We concluded that time-series cytokine profiles in breast cancer patients may be useful for understanding immune cell activity during treatment and for predicting treatment outcomes, supporting precision medicine. TRIAL REGISTRATION: The study has been registered with the University Hospital Medical Information Network Clinical Trials Registry ( http://www.umin.ac.jp/ctr/index-j.htm ) with the unique trial number UMIN000023162. The association Japan Breast Cancer Research Group trial number is JBCRG-22. The clinical outcome of the JBCRG-22 study was published in Breast Cancer Research and Treatment on 25 March 2021. https://doi.org/10.1007/s10549-021-06184-w .

Nicotinamide-N-methyltransferase regulates lipid metabolism via SAM and 1-methylnicotinamide in the AML12 hepatocyte cell line.

Nicotinamide-N-methyltransferase (NNMT) is an enzyme that consumes S-adenosyl-methionine (SAM) and nicotinamide (NAM) to produce S-adenosyl-homocysteine (SAH) and 1-methylnicotinamide (MNAM). How much NNMT contributes to the quantity regulation of these four metabolites depends on whether NNMT is a major consumer or producer of these metabolites, which varies among various cellular contexts. Yet, whether NNMT critically regulates these metabolites in the AML12 hepatocyte cell line has been unexplored. To address this, we knockdown Nnmt in AML12 cells and investigate the effects of Nnmt RNAi on metabolism and gene expression. We find that Nnmt RNAi accumulates SAM and SAH, whereas it reduces MNAM with NAM being unaltered. These results indicate that NNMT is a significant consumer of SAM and critical for MNAM production in this cell line. Moreover, transcriptome analyses reveal that altered SAM and MNAM homeostasis is accompanied by various detrimental molecular phenotypes, as exemplified by the down-regulations of lipogenic genes, such as Srebf1. Consistent with this, oil-red O-staining experiments demonstrate the decrease of total neutral lipids upon Nnmt RNAi. Treating Nnmt RNAi AML12 cells with cycloleucine, an inhibitor of SAM biogenesis suppresses SAM accumulation and rescues the decrease of neutral lipids. MNAM also shows activity to elevate neutral lipids. These results suggest that NNMT contributes to lipid metabolism by maintaining proper SAM and MNAM homeostasis. This study provides an additional example where NNMT plays a critical role in regulating SAM and MNAM metabolism.

Serum amyloid alpha 1-2 are not required for liver inflammation in the 4T1 murine breast cancer model.

Cancers induce the production of acute phase proteins such as serum amyloid alpha (SAA) in the liver and cause inflammation in various host organs. Despite the well-known coincidence of acute phase response and inflammation, the direct roles of SAA proteins in inflammation in the cancer context remains incompletely characterized, particularly in vivo. Here, we investigate the in vivo significance of SAA proteins in liver inflammation in the 4T1 murine breast cancer model. 4T1 cancers elevate the expression of SAA1 and SAA2, the two major murine acute phase proteins in the liver. The elevation of Saa1-2 correlates with the up-regulation of immune cell-related genes including neutrophil markers. To examine this correlation in detail, we generate mice that lack Saa1-2 and investigate immune-cell phenotypes. RNA-seq experiments reveal that deletion of Saa1-2 does not strongly affect 4T1-induced activation of immune cell-related genes in the liver. Flow cytometry experiments demonstrate the dispensable roles of SAA1-2 in cancer-dependent neutrophil infiltration to the liver. Consistently, 4T1-induced gene expression changes in bone marrow do not require Saa1-2. This study clarifies the negligible contribution of SAA1-2 proteins in liver inflammation in the 4T1 breast cancer model.

Murine breast cancers disorganize the liver transcriptome in a zonated manner.

The spatially organized gene expression program within the liver specifies hepatocyte functions according to their relative distances to the bloodstream (i.e., zonation), contributing to liver homeostasis. Despite the knowledge that solid cancers remotely disrupt liver homeostasis, it remains unexplored whether solid cancers affect liver zonation. Here, using spatial transcriptomics, we thoroughly investigate the abundance and zonation of hepatic genes in cancer-bearing mice. We find that breast cancers affect liver zonation in various distinct manners depending on biological pathways. Aspartate metabolism and triglyceride catabolic processes retain relatively intact zonation patterns, but the zonation of xenobiotic catabolic process genes exhibits a strong disruption. The acute phase response is induced in zonated manners. Furthermore, we demonstrate that breast cancers activate innate immune cells in particular neutrophils in distinct zonated manners, rather than in a uniform fashion within the liver. Collectively, breast cancers disorganize hepatic transcriptomes in zonated manners, thereby disrupting zonated functions of the liver.

Tracing the evolutionary history of blood cells to the unicellular ancestor of animals.

Blood cells are thought to have emerged as phagocytes in the common ancestor of animals followed by the appearance of novel blood cell lineages such as thrombocytes, erythrocytes, and lymphocytes, during evolution. However, this speculation is not based on genetic evidence and it is still possible to argue that phagocytes in different species have different origins. It also remains to be clarified how the initial blood cells evolved; whether ancient animals have solely developed de novo programs for phagocytes or they have inherited a key program from ancestral unicellular organisms. Here, we traced the evolutionary history of blood cells, and cross-species comparison of gene expression profiles revealed that phagocytes in various animal species and Capsaspora (C.) owczarzaki, a unicellular organism, are transcriptionally similar to each other. We also found that both phagocytes and C. owczarzaki share a common phagocytic program, and that CEBPα is the sole transcription factor highly expressed in both phagocytes and C. owczarzaki. We further showed that the function of CEBPα to drive phagocyte program in nonphagocytic blood cells has been conserved in tunicate, sponge, and C. owczarzaki. We finally showed that, in murine hematopoiesis, repression of CEBPα to maintain nonphagocytic lineages is commonly achieved by polycomb complexes. These findings indicate that the initial blood cells emerged inheriting a unicellular organism program driven by CEBPα and that the program has also been seamlessly inherited in phagocytes of various animal species throughout evolution.

Remote solid cancers rewire hepatic nitrogen metabolism via host nicotinamide-N-methyltransferase.

Cancers disrupt host homeostasis in various manners but the identity of host factors underlying such disruption remains largely unknown. Here we show that nicotinamide-N-methyltransferase (NNMT) is a host factor that mediates metabolic dysfunction in the livers of cancer-bearing mice. Multiple solid cancers distantly increase expression of Nnmt and its product 1-methylnicotinamide (MNAM) in the liver. Multi-omics analyses reveal suppression of the urea cycle accompanied by accumulation of amino acids, and enhancement of uracil biogenesis in the livers of cancer-bearing mice. Importantly, genetic deletion of Nnmt leads to alleviation of these metabolic abnormalities, and buffers cancer-dependent weight loss and reduction of the voluntary wheel-running activity. Our data also demonstrate that MNAM is capable of affecting urea cycle metabolites in the liver. These results suggest that cancers up-regulate the hepatic NNMT pathway to rewire liver metabolism towards uracil biogenesis rather than nitrogen disposal via the urea cycle, thereby disrupting host homeostasis.

Retraction: Transcriptome analysis reveals a role for the endothelial ANP-GC-A signaling in interfering with pre-metastatic niche formation by solid cancers.

[This retracts the article DOI: 10.18632/oncotarget.18032.].

Identification of a genomic enhancer that enforces proper apoptosis induction in thymic negative selection.

During thymic negative selection, autoreactive thymocytes carrying T cell receptor (TCR) with overtly strong affinity to self-MHC/self-peptide are removed by Bim-dependent apoptosis, but how Bim is specifically regulated to link TCR activation and apoptosis induction is unclear. Here we identify a murine T cell-specific genomic enhancer EBAB (Bub1-Acoxl-Bim), whose deletion leads to accumulation of thymocytes expressing high affinity TCRs. Consistently, EBAB knockout mice have defective negative selection and fail to delete autoreactive thymocytes in various settings, with this defect accompanied by reduced Bim expression and apoptosis induction. By contrast, EBAB is dispensable for maintaining peripheral T cell homeostasis via Bim-dependent pathways. Our data thus implicate EBAB as an important, developmental stage-specific regulator of Bim expression and apoptosis induction to enforce thymic negative selection and suppress autoimmunity. Our study unravels a part of genomic enhancer codes that underlie complex and context-dependent gene regulation in TCR signaling.

A novel zebrafish intestinal tumor model reveals a role for cyp7a1-dependent tumor-liver crosstalk in causing adverse effects on the host.

The nature of host organs and genes that underlie tumor-induced physiological disruption on the host remains ill-defined. Here, we establish a novel zebrafish intestinal tumor model that is suitable for addressing this issue, and find that hepatic cyp7a1, the rate-limiting factor for synthesizing bile acids [or, in the case of zebrafish, bile alcohol (BA)], is such a host gene. Inducing krasG12D by Gal4 specifically expressed in the posterior intestine resulted in the formation of an intestinal tumor. The local intestinal tumor caused systemic detrimental effects on the host, including liver inflammation, hepatomegaly, growth defects and organismal death. Whole-organism-level gene expression analysis and metabolite measurements revealed that the intestinal tumor reduced total BA levels, possibly via altered expression of hepatic cyp7a1 Genetically overexpressing cyp7a1 in the liver restored BA synthesis and ameliorated tumor-induced liver inflammation, but not other tumor-dependent phenotypes. Thus, we found a previously unknown role of cyp7a1 as the host gene that links the intestinal tumor, hepatic cholesterol-BA metabolism and liver inflammation in tumor-bearing zebrafish larvae. Our model provides an important basis to discover host genes responsible for tumor-induced phenotypes and to uncover mechanisms underlying how tumors adversely affect host organisms.

Transcriptome analysis reveals a role for the endothelial ANP-GC-A signaling in interfering with pre-metastatic niche formation by solid cancers.

Cancer establishes a microenvironment called the pre-metastatic niche in distant organs where disseminated cancer cells can efficiently metastasize. Pre-metastatic niche formation requires various genetic factors. Previous studies suggest that inhibiting a single niche-factor is insufficient to completely block pre-metastatic niche formation especially in human patients. Here we show that the atrial natriuretic peptide (ANP), an endogenous hormone produced by the heart, inhibits pre-metastatic niche formation and metastasis of murine solid cancer models when pharmacologically supplied in vivo. On the basis of a wealth of comprehensive RNA-seq data, we demonstrated that ANP globally suppressed expression of cancer-induced genes including known niche-factors in the lung. The lungs of mice overexpressing GC-A, a receptor for ANP in endothelial cells, were conferred resistance against pre-metastatic niche formation. Importantly, neither ANP administration nor GC-A overexpression had a detrimental effect on lung gene expression in a cancer-free condition. The current study establishes endothelial ANP-GC-A signaling as a therapeutic target to control the pre-metastatic niche.

Remote reprogramming of hepatic circadian transcriptome by breast cancer.

Cancers adversely affect organismal physiology. To date, the genes within a patient responsible for systemically spreading cancer-induced physiological disruption remain elusive. To identify host genes responsible for transmitting disruptive, cancer-driven signals, we thoroughly analyzed the transcriptome of a suite of host organs from mice bearing 4T1 breast cancer, and discovered complexly rewired patterns of circadian gene expression in the liver. Our data revealed that 7 core clock transcription factors, represented by Rev-erba and Rorg, exhibited abnormal daily expression rhythm in the liver of 4T1-bearing mice. Accordingly, expression patterns of specific set of downstream circadian genes were compromised. Osgin1, a marker for oxidative stress, was an example. Specific downstream genes, including E2f8, a transcriptional repressor that controls cellular polyploidy, displayed a striking pattern of disruption, "day-night reversal." Meanwhile, we found that the liver of 4T1-bearing mice suffered from increased oxidative stress. The tetraploid hepatocytes population was concomitantly increased in 4T1-bearing mice, which has not been previously appreciated as a cancer-induced phenotype. In summary, the current study provides a comprehensive characterization of the 4T1-affected hepatic circadian transcriptome that possibly underlies cancer-induced physiological alteration in the liver.

The transcriptional cofactor TRIM33 prevents apoptosis in B lymphoblastic leukemia by deactivating a single enhancer.

Most mammalian transcription factors (TFs) and cofactors occupy thousands of genomic sites and modulate the expression of large gene networks to implement their biological functions. In this study, we describe an exception to this paradigm. TRIM33 is identified here as a lineage dependency in B cell neoplasms and is shown to perform this essential function by associating with a single cis element. ChIP-seq analysis of TRIM33 in murine B cell leukemia revealed a preferential association with two lineage-specific enhancers that harbor an exceptional density of motifs recognized by the PU.1 TF. TRIM33 is recruited to these elements by PU.1, yet acts to antagonize PU.1 function. One of the PU.1/TRIM33 co-occupied enhancers is upstream of the pro-apoptotic gene Bim, and deleting this enhancer renders TRIM33 dispensable for leukemia cell survival. These findings reveal an essential role for TRIM33 in preventing apoptosis in B lymphoblastic leukemia by interfering with enhancer-mediated Bim activation.

Silkworm HP1a transcriptionally enhances highly expressed euchromatic genes via association with their transcription start sites.

Heterochromatin protein 1 (HP1) is an evolutionarily conserved protein across different eukaryotic species and is crucial for heterochromatin establishment and maintenance. The silkworm, Bombyx mori, encodes two HP1 proteins, BmHP1a and BmHP1b. In order to investigate the role of BmHP1a in transcriptional regulation, we performed genome-wide analyses of the transcriptome, transcription start sites (TSSs), chromatin modification states and BmHP1a-binding sites of the silkworm ovary-derived BmN4 cell line. We identified a number of BmHP1a-binding loci throughout the silkworm genome and found that these loci included TSSs and frequently co-occurred with neighboring euchromatic histone modifications. In addition, we observed that genes with BmHP1a-associated TSSs were relatively highly expressed in BmN4 cells. RNA interference-mediated BmHP1a depletion resulted in the transcriptional repression of highly expressed genes with BmHP1a-associated TSSs, whereas genes not coupled with BmHP1a-binding regions were less affected by the treatment. These results demonstrate that BmHP1a binds near TSSs of highly expressed euchromatic genes and positively regulates their expression. Our study revealed a novel mode of transcriptional regulation mediated by HP1 proteins.

A single female-specific piRNA is the primary determiner of sex in the silkworm.

The silkworm Bombyx mori uses a WZ sex determination system that is analogous to the one found in birds and some reptiles. In this system, males have two Z sex chromosomes, whereas females have Z and W sex chromosomes. The silkworm W chromosome has a dominant role in female determination, suggesting the existence of a dominant feminizing gene in this chromosome. However, the W chromosome is almost fully occupied by transposable element sequences, and no functional protein-coding gene has been identified so far. Female-enriched PIWI-interacting RNAs (piRNAs) are the only known transcripts that are produced from the sex-determining region of the W chromosome, but the function(s) of these piRNAs are unknown. Here we show that a W-chromosome-derived, female-specific piRNA is the feminizing factor of B. mori. This piRNA is produced from a piRNA precursor which we named Fem. Fem sequences were arranged in tandem in the sex-determining region of the W chromosome. Inhibition of Fem-derived piRNA-mediated signalling in female embryos led to the production of the male-specific splice variants of B. mori doublesex (Bmdsx), a gene which acts at the downstream end of the sex differentiation cascade. A target gene of Fem-derived piRNA was identified on the Z chromosome of B. mori. This gene, which we named Masc, encoded a CCCH-type zinc finger protein. We show that the silencing of Masc messenger RNA by Fem piRNA is required for the production of female-specific isoforms of Bmdsx in female embryos, and that Masc protein controls both dosage compensation and masculinization in male embryos. Our study characterizes a single small RNA that is responsible for primary sex determination in the WZ sex determination system.

Making piRNAs in vitro.

Classical biochemical approaches have made great contribution to our current understanding of how small interfering RNA (siRNA) and microRNA (miRNA) are produced and function. However, it has been challenging to establish in vitro systems that can dissect the mechanism of PIWI-interacting RNAs (piRNAs), largely due to the lack of suitable cultured cell lines possessing the piRNA pathway endogenously. The silkworm ovary-derived BmN4 cell line is an emerging model for studying piRNAs, because this cell line harbors a fully functional germline piRNA pathway. We recently reported in vitro recapitulation of a part of piRNA biogenesis-loading of precursor piRNAs and 3'-end maturation-by using lysate prepared from BmN4 cells. Here we describe the methods for maintenance of BmN4 cells, preparation of cell lysate, and in vitro assays to dissect piRNA biogenesis.

Role of SWI/SNF in acute leukemia maintenance and enhancer-mediated Myc regulation.

Cancer cells frequently depend on chromatin regulatory activities to maintain a malignant phenotype. Here, we show that leukemia cells require the mammalian SWI/SNF chromatin remodeling complex for their survival and aberrant self-renewal potential. While Brg1, an ATPase subunit of SWI/SNF, is known to suppress tumor formation in several cell types, we found that leukemia cells instead rely on Brg1 to support their oncogenic transcriptional program, which includes Myc as one of its key targets. To account for this context-specific function, we identify a cluster of lineage-specific enhancers located 1.7 Mb downstream from Myc that are occupied by SWI/SNF as well as the BET protein Brd4. Brg1 is required at these distal elements to maintain transcription factor occupancy and for long-range chromatin looping interactions with the Myc promoter. Notably, these distal Myc enhancers coincide with a region that is focally amplified in ∼3% of acute myeloid leukemias. Together, these findings define a leukemia maintenance function for SWI/SNF that is linked to enhancer-mediated gene regulation, providing general insights into how cancer cells exploit transcriptional coactivators to maintain oncogenic gene expression programs.

Large scale full-length cDNA sequencing reveals a unique genomic landscape in a lepidopteran model insect, Bombyx mori.

The establishment of a complete genomic sequence of silkworm, the model species of Lepidoptera, laid a foundation for its functional genomics. A more complete annotation of the genome will benefit functional and comparative studies and accelerate extensive industrial applications for this insect. To realize these goals, we embarked upon a large-scale full-length cDNA collection from 21 full-length cDNA libraries derived from 14 tissues of the domesticated silkworm and performed full sequencing by primer walking for 11,104 full-length cDNAs. The large average intron size was 1904 bp, resulting from a high accumulation of transposons. Using gene models predicted by GLEAN and published mRNAs, we identified 16,823 gene loci on the silkworm genome assembly. Orthology analysis of 153 species, including 11 insects, revealed that among three Lepidoptera including Monarch and Heliconius butterflies, the 403 largest silkworm-specific genes were composed mainly of protective immunity, hormone-related, and characteristic structural proteins. Analysis of testis-/ovary-specific genes revealed distinctive features of sexual dimorphism, including depletion of ovary-specific genes on the Z chromosome in contrast to an enrichment of testis-specific genes. More than 40% of genes expressed in specific tissues mapped in tissue-specific chromosomal clusters. The newly obtained FL-cDNA sequences enabled us to annotate the genome of this lepidopteran model insect more accurately, enhancing genomic and functional studies of Lepidoptera and comparative analyses with other insect orders, and yielding new insights into the evolution and organization of lepidopteran-specific genes.