Epigenetic characterization of housekeeping core promoters and their importance in tumor suppression.

In this research, we elucidate the presence of around 11,000 housekeeping cis-regulatory elements (HK-CREs) and describe their main characteristics. Besides the trivial promoters of housekeeping genes, most HK-CREs reside in promoter regions and are involved in a broader role beyond housekeeping gene regulation. HK-CREs are conserved regions rich in unmethylated CpG sites. Their distribution highly correlates with that of protein-coding genes, and they interact with many genes over long distances. We observed reduced activity of a subset of HK-CREs in diverse cancer subtypes due to aberrant methylation, particularly those located in chromosome 19 and associated with zinc finger genes. Further analysis of samples from 17 cancer subtypes showed a significantly increased survival probability of patients with higher expression of these genes, suggesting them as housekeeping tumor suppressor genes. Overall, our work unravels the presence of housekeeping CREs indispensable for the maintenance and stability of cells.

Regulation of lymphoid-myeloid lineage bias through regnase-1/3-mediated control of Nfkbiz.

ABSTRACT: Regulation of lineage biases in hematopoietic stem and progenitor cells (HSPCs) is pivotal for balanced hematopoietic output. However, little is known about the mechanism behind lineage choice in HSPCs. Here, we show that messenger RNA (mRNA) decay factors regnase-1 (Reg1; Zc3h12a) and regnase-3 (Reg3; Zc3h12c) are essential for determining lymphoid fate and restricting myeloid differentiation in HSPCs. Loss of Reg1 and Reg3 resulted in severe impairment of lymphopoiesis and a mild increase in myelopoiesis in the bone marrow. Single-cell RNA sequencing analysis revealed that Reg1 and Reg3 regulate lineage directions in HSPCs via the control of a set of myeloid-related genes. Reg1- and Reg3-mediated control of mRNA encoding Nfkbiz, a transcriptional and epigenetic regulator, was essential for balancing lymphoid/myeloid lineage output in HSPCs in vivo. Furthermore, single-cell assay for transposase-accessible chromatin sequencing analysis revealed that Reg1 and Reg3 control the epigenetic landscape on myeloid-related gene loci in early stage HSPCs via Nfkbiz. Consistently, an antisense oligonucleotide designed to inhibit Reg1- and Reg3-mediated Nfkbiz mRNA degradation primed hematopoietic stem cells toward myeloid lineages by enhancing Nfkbiz expression. Collectively, the collaboration between posttranscriptional control and chromatin remodeling by the Reg1/Reg3-Nfkbiz axis governs HSPC lineage biases, ultimately dictating the fate of lymphoid vs myeloid differentiation.

Hematopoietic cell-derived IL-15 supports NK cell development in scattered and clustered localization within the bone marrow.

Natural killer (NK) cells are innate immune cells critical for protective immune responses against infection and cancer. Although NK cells differentiate in the bone marrow (BM) in an interleukin-15 (IL-15)-dependent manner, the cellular source of IL-15 remains elusive. Using NK cell reporter mice, we show that NK cells are localized in the BM in scattered and clustered manners. NK cell clusters overlap with monocyte and dendritic cell accumulations, whereas scattered NK cells require CXCR4 signaling. Using cell-specific IL-15-deficient mice, we show that hematopoietic cells, but not stromal cells, support NK cell development in the BM through IL-15. In particular, IL-15 produced by monocytes and dendritic cells appears to contribute to NK cell development. These results demonstrate that hematopoietic cells are the IL-15 niche for NK cell development in the BM and that BM NK cells are present in scattered and clustered compartments by different mechanisms, suggesting their distinct functions in the immune response.

A universal tool for predicting differentially active features in single-cell and spatial genomics data.

With the growing complexity of single-cell and spatial genomics data, there is an increasing importance of unbiased and efficient exploratory data analysis tools. One common exploratory data analysis step is the prediction of genes with different levels of activity in a subset of cells or locations inside a tissue. We previously developed singleCellHaystack, a method for predicting differentially expressed genes from single-cell transcriptome data, without relying on comparisons between clusters of cells. Here we present an update to singleCellHaystack, which is now a universally applicable method for predicting differentially active features: (1) singleCellHaystack now accepts continuous features that can be RNA or protein expression, chromatin accessibility or module scores from single-cell, spatial and even bulk genomics data, and (2) it can handle 1D trajectories, 2-3D spatial coordinates, as well as higher-dimensional latent spaces as input coordinates. Performance has been drastically improved, with up to ten times reduction in computational time and scalability to millions of cells, making singleCellHaystack a suitable tool for exploratory analysis of atlas level datasets. singleCellHaystack is available as packages in both R and Python.

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.

The N6-methyladenosine methyltransferase METTL16 enables erythropoiesis through safeguarding genome integrity.

During erythroid differentiation, the maintenance of genome integrity is key for the success of multiple rounds of cell division. However, molecular mechanisms coordinating the expression of DNA repair machinery in erythroid progenitors are poorly understood. Here, we discover that an RNA N6-methyladenosine (m6A) methyltransferase, METTL16, plays an essential role in proper erythropoiesis by safeguarding genome integrity via the control of DNA-repair-related genes. METTL16-deficient erythroblasts exhibit defective differentiation capacity, DNA damage and activation of the apoptotic program. Mechanistically, METTL16 controls m6A deposition at the structured motifs in DNA-repair-related transcripts including Brca2 and Fancm mRNAs, thereby upregulating their expression. Furthermore, a pairwise CRISPRi screen revealed that the MTR4-nuclear RNA exosome complex is involved in the regulation of METTL16 substrate mRNAs in erythroblasts. Collectively, our study uncovers that METTL16 and the MTR4-nuclear RNA exosome act as essential regulatory machinery to maintain genome integrity and erythropoiesis.

Regnase-1 Prevents Pulmonary Arterial Hypertension Through mRNA Degradation of Interleukin-6 and Platelet-Derived Growth Factor in Alveolar Macrophages.

BACKGROUND: Pulmonary arterial hypertension (PAH) is a type of pulmonary hypertension (PH) characterized by obliterative pulmonary vascular remodeling, resulting in right-sided heart failure. Although the pathogenesis of PAH is not fully understood, inflammatory responses and cytokines have been shown to be associated with PAH, in particular, with connective tissue disease-PAH. In this sense, Regnase-1, an RNase that regulates mRNAs encoding genes related to immune reactions, was investigated in relation to the pathogenesis of PH. METHODS: We first examined the expression levels of ZC3H12A (encoding Regnase-1) in peripheral blood mononuclear cells from patients with PH classified under various types of PH, searching for an association between the ZC3H12A expression and clinical features. We then generated mice lacking Regnase-1 in myeloid cells, including alveolar macrophages, and examined right ventricular systolic pressures and histological changes in the lung. We further performed a comprehensive analysis of the transcriptome of alveolar macrophages and pulmonary arteries to identify genes regulated by Regnase-1 in alveolar macrophages. RESULTS: ZC3H12A expression in peripheral blood mononuclear cells was inversely correlated with the prognosis and severity of disease in patients with PH, in particular, in connective tissue disease-PAH. The critical role of Regnase-1 in controlling PAH was also reinforced by the analysis of mice lacking Regnase-1 in alveolar macrophages. These mice spontaneously developed severe PAH, characterized by the elevated right ventricular systolic pressures and irreversible pulmonary vascular remodeling, which recapitulated the pathology of patients with PAH. Transcriptomic analysis of alveolar macrophages and pulmonary arteries of these PAH mice revealed that Il6, Il1b, and Pdgfa/b are potential targets of Regnase-1 in alveolar macrophages in the regulation of PAH. The inhibition of IL-6 (interleukin-6) by an anti-IL-6 receptor antibody or platelet-derived growth factor by imatinib but not IL-1ß (interleukin-1ß) by anakinra, ameliorated the pathogenesis of PAH. CONCLUSIONS: Regnase-1 maintains lung innate immune homeostasis through the control of IL-6 and platelet-derived growth factor in alveolar macrophages, thereby suppressing the development of PAH in mice. Furthermore, the decreased expression of Regnase-1 in various types of PH implies its involvement in PH pathogenesis and may serve as a disease biomarker, and a therapeutic target for PH as well.

Delamination of trophoblast-like syncytia from the amniotic ectodermal analogue in human primed embryonic stem cell-based differentiation model.

Human primed embryonic stem cells (ESCs) are known to be converted to cells with several trophoblast properties, but it has remained controversial whether this phenomenon represents the inherent differentiation competence of human primed ESCs to trophoblast lineages. In this study, we report that chemical blockage of ACTIVIN/NODAL and FGF signals is sufficient to steer human primed ESCs into GATA3-expressing cells that give rise to placental hormone-producing syncytia analogous to syncytiotrophoblasts of the post-implantation stage of the human embryo. Despite their cytological similarity to syncytiotrophoblasts, these syncytia arise from the non-trophoblastic differentiation trajectory that recapitulates amniogenesis. These results provide insights into the possible extraembryonic differentiation pathway that is unique in primate embryogenesis.

Enhancement of Regnase-1 expression with stem loop-targeting antisense oligonucleotides alleviates inflammatory diseases.

Regnase-1 is an ribonuclease that plays essential roles in restricting inflammation through degrading messenger RNAs (mRNAs) involved in immune reactions via the recognition of stem-loop (SL) structures in the 3' untranslated regions (3'UTRs). Dysregulated expression of Regnase-1 is associated with the pathogenesis of inflammatory and autoimmune diseases in mice and humans. Here, we developed a therapeutic strategy to suppress inflammatory responses by blocking Regnase-1 self-regulation, which was mediated by the simultaneous use of two antisense phosphorodiamidate morpholino oligonucleotides (MOs) to alter the binding of Regnase-1 toward the SL structures in its 3'UTR. Regnase-1-targeting MOs not only enhanced Regnase-1 expression by stabilizing mRNAs but also effectively reduced the expression of multiple proinflammatory transcripts that were controlled by Regnase-1 in macrophages. Intratracheal administration of Regnase-1-targeting MOs ameliorated acute respiratory distress syndrome and chronic fibrosis through suppression of inflammatory cascades. In addition, intracranial treatment with Regnase-1-targeting MOs attenuated the development of experimental autoimmune encephalomyelitis by promoting the expansion of homeostatic microglia and regulatory T cell populations. Regnase-1 expression was inversely correlated with disease severity in patients with multiple sclerosis, and MOs targeting human Regnase-1 SL structures were effective in mitigating cytokine production in human immune cells. Collectively, MO-mediated disruption of the Regnase-1 self-regulation pathway is a potential therapeutic strategy to enhance Regnase-1 abundance, which, in turn, provides therapeutic benefits for treating inflammatory diseases by suppressing inflammation.

Cyclin J-CDK complexes limit innate immune responses by reducing proinflammatory changes in macrophage metabolism.

Toll-like receptor (TLR) stimulation induces glycolysis and the production of mitochondrial reactive oxygen species (ROS), both of which are critical for inflammatory responses in macrophages. Here, we demonstrated that cyclin J, a TLR-inducible member of the cyclin family, reduced cytokine production in macrophages by coordinately controlling glycolysis and mitochondrial functions. Cyclin J interacted with cyclin-dependent kinases (CDKs), which increased the phosphorylation of a subset of CDK substrates, including the transcription factor FoxK1 and the GTPase Drp1. Cyclin J-dependent phosphorylation of FoxK1 decreased the transcription of glycolytic genes and Hif-1α activation, whereas hyperactivation of Drp1 by cyclin J-dependent phosphorylation promoted mitochondrial fragmentation and impaired the production of mitochondrial ROS. In mice, cyclin J in macrophages limited the growth of tumor xenografts and protected against LPS-induced shock but increased the susceptibility to bacterial infection. Collectively, our findings indicate that cyclin J-CDK signaling promotes antitumor immunity and the resolution of inflammation by opposing the metabolic changes that drive inflammatory responses in macrophages.

Evaluation of critical data processing steps for reliable prediction of gene co-expression from large collections of RNA-seq data.

MOTIVATION: Gene co-expression analysis is an attractive tool for leveraging enormous amounts of public RNA-seq datasets for the prediction of gene functions and regulatory mechanisms. However, the optimal data processing steps for the accurate prediction of gene co-expression from such large datasets remain unclear. Especially the importance of batch effect correction is understudied. RESULTS: We processed RNA-seq data of 68 human and 76 mouse cell types and tissues using 50 different workflows into 7,200 genome-wide gene co-expression networks. We then conducted a systematic analysis of the factors that result in high-quality co-expression predictions, focusing on normalization, batch effect correction, and measure of correlation. We confirmed the key importance of high sample counts for high-quality predictions. However, choosing a suitable normalization approach and applying batch effect correction can further improve the quality of co-expression estimates, equivalent to a >80% and >40% increase in samples. In larger datasets, batch effect removal was equivalent to a more than doubling of the sample size. Finally, Pearson correlation appears more suitable than Spearman correlation, except for smaller datasets. CONCLUSION: A key point for accurate prediction of gene co-expression is the collection of many samples. However, paying attention to data normalization, batch effects, and the measure of correlation can significantly improve the quality of co-expression estimates.

Profibrotic function of pulmonary group 2 innate lymphoid cells is controlled by regnase-1.

Regnase-1 is an RNase critical for post-transcriptional control of pulmonary immune homeostasis in mice by degrading immune-related mRNAs. However, little is known about the cell types Regnase-1 controls in the lung, and its relevance to human pulmonary diseases.Regnase-1-dependent changes in lung immune cell types were examined by a competitive bone marrow transfer mouse model, and group 2 innate lymphoid cells (ILC2s) were identified. Then the associations between Regnase-1 in ILC2s and human diseases were investigated by transcriptome analysis and a bleomycin-induced pulmonary fibrosis mouse model. The clinical significance of Regnase-1 in ILC2s was further assessed using patient-derived cells.Regnase-1-deficiency resulted in the spontaneous proliferation and activation of ILC2s in the lung. Intriguingly, genes associated with pulmonary fibrosis were highly upregulated in Regnase-1-deficient ILC2s compared with wild-type, and supplementation of Regnase-1-deficient ILC2s augmented bleomycin-induced pulmonary fibrosis in mice. Regnase-1 suppresses mRNAs encoding transcription factors Gata3 and Egr1, which are potent to regulate fibrosis-associated genes. Clinically, Regnase-1 protein levels in ILC2 negatively correlated with the ILC2 population in bronchoalveolar lavage fluid. Furthermore, idiopathic pulmonary fibrosis (IPF) patients with ILC2s >1500 cells·mL-1 peripheral blood exhibited poorer prognosis than patients with lower numbers, implying the contribution of Regnase-1 in ILC2s for the progression of IPF.Collectively, Regnase-1 was identified as a critical post-transcriptional regulator of the profibrotic function of ILC2s both in mouse and human, suggesting that Regnase-1 may be a novel therapeutic target for IPF.

A clustering-independent method for finding differentially expressed genes in single-cell transcriptome data.

A common analysis of single-cell sequencing data includes clustering of cells and identifying differentially expressed genes (DEGs). How cell clusters are defined has important consequences for downstream analyses and the interpretation of results, but is often not straightforward. To address this difficulty, we present singleCellHaystack, a method that enables the prediction of DEGs without relying on explicit clustering of cells. Our method uses Kullback-Leibler divergence to find genes that are expressed in subsets of cells that are non-randomly positioned in a multidimensional space. Comparisons with existing DEG prediction approaches on artificial datasets show that singleCellHaystack has higher accuracy. We illustrate the usage of singleCellHaystack through applications on 136 real transcriptome datasets and a spatial transcriptomics dataset. We demonstrate that our method is a fast and accurate approach for DEG prediction in single-cell data. singleCellHaystack is implemented as an R package and is available from CRAN and GitHub.

Zinc Finger Protein St18 Protects against Septic Death by Inhibiting VEGF-A from Macrophages.

Zinc finger protein St18 was initially reported as candidate tumor suppressor gene, and also suggested that fibroblast St18 positively regulates NF-κB activation. Despite the pleiotropic functions of St18, little is known about its roles in macrophages. Here, we report that myeloid St18 is a potent inhibitor of VEGF-A. Mice lacking St18 in myeloid lineages exhibit increased retinal vasculature with enhanced serum VEGF-A concentrations. Despite the normal activation of NF-κB target genes, these mice are highly susceptible to LPS-induced shock, polymicrobial sepsis, and experimental colitis, accompanied by enhanced vascular and intestinal leakage. Pharmacological inhibition of VEGF signaling rescued the high mortality rate of myeloid-specific St18-deficient mice in response to inflammation. Mechanistically, St18 directly binds to Sp1 and attenuates its activity, leading to the suppression of Sp1 target gene VEGF-A. Using mouse genetic and pharmacological models, we reveal myeloid St18 as a critical septic death protector.

Codon bias confers stability to human mRNAs.

Codon bias has been implicated as one of the major factors contributing to mRNA stability in several model organisms. However, the molecular mechanisms of codon bias on mRNA stability remain unclear in humans. Here, we show that human cells possess a mechanism to modulate RNA stability through a unique codon bias. Bioinformatics analysis showed that codons could be clustered into two distinct groups-codons with G or C at the third base position (GC3) and codons with either A or T at the third base position (AT3): the former stabilizing while the latter destabilizing mRNA. Quantification of codon bias showed that increased GC3-content entails proportionately higher GC-content. Through bioinformatics, ribosome profiling, and in vitro analysis, we show that decoupling the effects of codon bias reveals two modes of mRNA regulation, one GC3- and one GC-content dependent. Employing an immunoprecipitation-based strategy, we identify ILF2 and ILF3 as RNA-binding proteins that differentially regulate global mRNA abundances based on codon bias. Our results demonstrate that codon bias is a two-pronged system that governs mRNA abundance.

Waves of chromatin modifications in mouse dendritic cells in response to LPS stimulation.

BACKGROUND: The importance of transcription factors (TFs) and epigenetic modifications in the control of gene expression is widely accepted. However, causal relationships between changes in TF binding, histone modifications, and gene expression during the response to extracellular stimuli are not well understood. Here, we analyze the ordering of these events on a genome-wide scale in dendritic cells in response to lipopolysaccharide (LPS) stimulation. RESULTS: Using a ChIP-seq time series dataset, we find that the LPS-induced accumulation of different histone modifications follows clearly distinct patterns. Increases in H3K4me3 appear to coincide with transcriptional activation. In contrast, H3K9K14ac accumulates early after stimulation, and H3K36me3 at later time points. Integrative analysis with TF binding data reveals potential links between TF activation and dynamics in histone modifications. Especially, LPS-induced increases in H3K9K14ac and H3K4me3 are associated with binding by STAT1/2 and were severely impaired in Stat1-/- cells. CONCLUSIONS: While the timing of short-term changes of some histone modifications coincides with changes in transcriptional activity, this is not the case for others. In the latter case, dynamics in modifications more likely reflect strict regulation by stimulus-induced TFs and their interactions with chromatin modifiers.

Pulmonary Regnase-1 orchestrates the interplay of epithelium and adaptive immune systems to protect against pneumonia.

Inhaled pathogens including Pseudomonas aeruginosa initially encounter airway epithelial cells (AECs), which are poised to evoke cell-intrinsic innate defense, affecting second tier of hematopoietic cell-mediated immune reaction. However, it is largely unknown how pulmonary immune responses mediated by a variety of immune cells are coordinated. Here we show that Regnase-1, an endoribonuclease expressed in AECs and immune cells, plays an essential role in coordinating innate responses and adaptive immunity against P. aeruginosa infection. Intratracheal treatment of mice with heat-killed P. aeruginosa resulted in prolonged disappearance of Regnase-1 consistent with sustained expression of Regnase-1 target inflammatory genes, whereas the transcription factor NF-κB was only transiently activated. AEC-specific deletion of Regnase-1 not only augmented innate defenses against P. aeruginosa but also enhanced secretion of Pseudomonas-specific IgA and Th17 accumulation in the lung, culminating in conferring significant resistance against P. aeruginosa re-infection in vivo. Although Regnase-1 directly controls distinct sets of genes in each of AECs and T cells, degradation of Regnase-1 in both cell types is beneficial for maximizing acquired immune responses. Collectively, these results demonstrate that Regnase-1 orchestrates AEC-mediated and immune cell-mediated host defense against pulmonary bacterial infection.

Addendum: Guidance of regulatory T cell development by Satb1-dependent super-enhancer establishment.

This corrects the article DOI: 10.1038/ni.3646.

Modeling the cis-regulatory modules of genes expressed in developmental stages of Drosophila melanogaster.

Because transcription is the first step in the regulation of gene expression, understanding how transcription factors bind to their DNA binding motifs has become absolutely necessary. It has been shown that the promoters of genes with similar expression profiles share common structural patterns. This paper presents an extensive study of the regulatory regions of genes expressed in 24 developmental stages of Drosophila melanogaster. It proposes the use of a combination of structural features, such as positioning of individual motifs relative to the transcription start site, orientation, pairwise distance between motifs, and presence of motifs anywhere in the promoter for predicting gene expression from structural features of promoter sequences. RNA-sequencing data was utilized to create and validate the 24 models. When genes with high-scoring promoters were compared to those identified by RNA-seq samples, 19 (79.2%) statistically significant models, a number that exceeds previous studies, were obtained. Each model yielded a set of highly informative features, which were used to search for genes with similar biological functions.

Regnase-1 Maintains Iron Homeostasis via the Degradation of Transferrin Receptor 1 and Prolyl-Hydroxylase-Domain-Containing Protein 3 mRNAs.

Iron metabolism is regulated by transcriptional and post-transcriptional mechanisms. The mRNA of the iron-controlling gene, transferrin receptor 1 (TfR1), has long been believed to be negatively regulated by a yet-unidentified endonuclease. Here, we show that the endonuclease Regnase-1 is critical for the degradation of mRNAs involved in iron metabolism in vivo. First, we demonstrate that Regnase-1 promotes TfR1 mRNA decay. Next, we show that Regnase-1-/- mice suffer from severe iron deficiency anemia, although hepcidin expression is downregulated. The iron deficiency anemia is induced by a defect in duodenal iron uptake. We reveal that duodenal Regnase-1 controls the expression of PHD3, which impairs duodenal iron uptake via HIF2α suppression. Finally, we show that Regnase-1 is a HIF2α-inducible gene and thus provides a positive feedback loop for HIF2α activation via PHD3. Collectively, these results demonstrate that Regnase-1-mediated regulation of iron-related transcripts is essential for the maintenance of iron homeostasis.