An atlas of combinatorial transcriptional regulation in mouse and man

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Smoking-associated gene expression alterations in nasal epithelium reveal immune impairment linked to lung cancer risk

BACKGROUND

Lung cancer is the leading cause of cancer-related death in the world. In contrast to many other cancers, a direct connection to modifiable lifestyle risk in the form of tobacco smoke has long been established. More than 50% of all smoking-related lung cancers occur in former smokers, 40% of which occur more than 15 years after smoking cessation. Despite extensive research, the molecular processes for persistent lung cancer risk remain unclear. We thus set out to examine whether risk stratification in the clinic and in the general population can be improved upon by the addition of genetic data and to explore the mechanisms of the persisting risk in former smokers.

METHODS

We analysed transcriptomic data from accessible airway tissues of 487 subjects, including healthy volunteers and clinic patients of different smoking statuses. We developed a computational model to assess smoking-associated gene expression changes and their reversibility after smoking is stopped, comparing healthy subjects to clinic patients with and without lung cancer.

RESULTS

We find persistent smoking-associated immune alterations to be a hallmark of the clinic patients. Integrating previous GWAS data using a transcriptional network approach, we demonstrate that the same immune- and interferon-related pathways are strongly enriched for genes linked to known genetic risk factors, demonstrating a causal relationship between immune alteration and lung cancer risk. Finally, we used accessible airway transcriptomic data to derive a non-invasive lung cancer risk classifier.

CONCLUSIONS

Our results provide initial evidence for germline-mediated personalized smoke injury response and risk in the general population, with potential implications for managing long-term lung cancer incidence and mortality.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced H3K27ac. Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data supports clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.

Intermittent fasting shifts the diurnal transcriptome atlas of transcription factors

Intermittent fasting remains a safe and effective strategy to ameliorate various age-related diseases, but its specific mechanisms are not fully understood. Considering that transcription factors (TFs) determine the response to environmental signals, here, we profiled the diurnal expression of 600 samples across four metabolic tissues sampled every 4 over 24 h from mice placed on five different feeding regimens to provide an atlas of TFs in biological space, time, and feeding regimen. Results showed that 1218 TFs exhibited tissue-specific and temporal expression profiles in ad libitum mice, of which 974 displayed significant oscillations at least in one tissue. Intermittent fasting triggered more than 90% (1161 in 1234) of TFs to oscillate somewhere in the body and repartitioned their tissue-specific expression. A single round of fasting generally promoted TF expression, especially in skeletal muscle and adipose tissues, while intermittent fasting mainly suppressed TF expression. Intermittent fasting down-regulated aging pathway and upregulated the pathway responsible for the inhibition of mammalian target of rapamycin (mTOR). Intermittent fasting shifts the diurnal transcriptome atlas of TFs, and mTOR inhibition may orchestrate intermittent fasting-induced health improvements. This atlas offers a reference and resource to understand how TFs and intermittent fasting may contribute to diurnal rhythm oscillation and bring about specific health benefits.

Cell-state dependent regulation of PPAR γ signaling by ZBTB9 in adipocytes

Adipocytes play a critical role in metabolic homeostasis. Peroxisome proliferator-activated receptor- γ (PPAR γ ) is a nuclear hormone receptor that is a master regulator of adipocyte differentiation and function. ZBTB9 was predicted to interact with PPAR γ based on large-scale protein interaction experiments. In addition, GWAS studies in the type 2 diabetes (T2D) Knowledge Portal revealed associations between Z btb9 and both BMI and T2D risk. Here we show that ZBTB9 positively regulates PPAR γ activity in mature adipocytes. Surprisingly Z btb9 knockdown (KD) also increased adipogenesis in 3T3-L1 cells and human preadipocytes. E2F activity was increased and E2F downstream target genes were upregulated in Zbtb9 -KD preadipocytes. Accordingly, RB phosphorylation, which regulates E2F activity, was enhanced in Zbtb9 -KD preadipocytes. Critically, an E2F1 inhibitor blocked the effects of Zbtb9 deficiency on adipogenic gene expression and lipid accumulation. Collectively, these results demonstrate that Zbtb9 inhibits adipogenesis as a negative regulator of Pparg expression via altered RB-E2F1 signaling. Our findings reveal complex cell-state dependent roles of ZBTB9 in adipocytes, identifying a new molecule that regulates adipogenesis and adipocyte biology as both a positive and negative regulator of PPAR γ signaling depending on the cellular context, and thus may be important in the pathogenesis and treatment of obesity and T2D.

Impaired bone morphogenetic protein (BMP) signaling pathways disrupt decidualization in endometriosis

Endometriosis is linked to increased infertility and pregnancy complications due to defective endometrial decidualization. We hypothesized that identification of altered signaling pathways during decidualization could identify the underlying cause of infertility and pregnancy complications. Our study reveals that transforming growth factor β (TGFβ) pathways are impaired in the endometrium of individuals with endometriosis, leading to defective decidualization. Through detailed transcriptomic analyses, we discovered abnormalities in TGFβ signaling pathways and key regulators, such as SMAD4, in the endometrium of affected individuals. We also observed compromised activity of bone morphogenetic proteins (BMP), a subset of the TGFβ family, that control endometrial receptivity. Using 3-dimensional models of endometrial stromal and epithelial assembloids, we showed that exogenous BMP2 improved decidual marker expression in individuals with endometriosis. Our findings reveal dysfunction of BMP/SMAD signaling in the endometrium of individuals with endometriosis, explaining decidualization defects and subsequent pregnancy complications in these individuals.

A time- and single-cell-resolved model of murine bone marrow hematopoiesis

The paradigmatic hematopoietic tree model is increasingly recognized to be limited, as it is based on heterogeneous populations largely defined by non-homeostatic assays testing cell fate potentials. Here, we combine persistent labeling with time-series single-cell RNA sequencing to build a real-time, quantitative model of in vivo tissue dynamics for murine bone marrow hematopoiesis. We couple cascading single-cell expression patterns with dynamic changes in differentiation and growth speeds. The resulting explicit linkage between molecular states and cellular behavior reveals widely varying self-renewal and differentiation properties across distinct lineages. Transplanted stem cells show strong acceleration of differentiation at specific stages of erythroid and neutrophil production, illustrating how the model can quantify the impact of perturbations. Our reconstruction of dynamic behavior from snapshot measurements is akin to how a kinetoscope allows sequential images to merge into a movie. We posit that this approach is generally applicable to understanding tissue-scale dynamics at high resolution.

A modular dCas9-based recruitment platform for combinatorial epigenome editing

Targeted epigenome editing tools allow precise manipulation and investigation of genome modifications, however they often display high context dependency and variable efficacy between target genes and cell types. While systems that simultaneously recruit multiple distinct 'effector' chromatin regulators can improve efficacy, they generally lack control over effector composition and spatial organisation. To overcome this we have created a modular combinatorial epigenome editing platform, called SSSavi. This system is an interchangeable and reconfigurable docking platform fused to dCas9 that enables simultaneous recruitment of up to four different effectors, allowing precise control of effector composition and spatial ordering. We demonstrate the activity and specificity of the SSSavi system and, by testing it against existing multi-effector targeting systems, demonstrate its comparable efficacy. Furthermore, we demonstrate the importance of the spatial ordering of the recruited effectors for effective transcriptional regulation. Together, the SSSavi system enables exploration of combinatorial effector co-recruitment to enhance manipulation of chromatin contexts previously resistant to targeted editing.

KnockTF 2.0: a comprehensive gene expression profile database with knockdown/knockout of transcription (co-)factors in multiple species

Transcription factors (TFs), transcription co-factors (TcoFs) and their target genes perform essential functions in diseases and biological processes. KnockTF 2.0 (http://www.licpathway.net/KnockTF/index.html) aims to provide comprehensive gene expression profile datasets before/after T(co)F knockdown/knockout across multiple tissue/cell types of different species. Compared with KnockTF 1.0, KnockTF 2.0 has the following improvements: (i) Newly added T(co)F knockdown/knockout datasets in mice, Arabidopsis thaliana and Zea mays and also an expanded scale of datasets in humans. Currently, KnockTF 2.0 stores 1468 manually curated RNA-seq and microarray datasets associated with 612 TFs and 172 TcoFs disrupted by different knockdown/knockout techniques, which are 2.5 times larger than those of KnockTF 1.0. (ii) Newly added (epi)genetic annotations for T(co)F target genes in humans and mice, such as super-enhancers, common SNPs, methylation sites and chromatin interactions. (iii) Newly embedded and updated search and analysis tools, including T(co)F Enrichment (GSEA), Pathway Downstream Analysis and Search by Target Gene (BLAST). KnockTF 2.0 is a comprehensive update of KnockTF 1.0, which provides more T(co)F knockdown/knockout datasets and (epi)genetic annotations across multiple species than KnockTF 1.0. KnockTF 2.0 facilitates not only the identification of functional T(co)Fs and target genes but also the investigation of their roles in the physiological and pathological processes.

TFCheckpoint database update, a cross-referencing system for transcription factors from human, mouse and rat

Prior knowledge about DNA-binding transcription factors (dbTFs), transcription co-regulators (coTFs) and general transcriptional factors (GTFs) is crucial for the study and understanding of the regulation of transcription. This is reflected by the many publications and database resources describing knowledge about TFs. We previously launched the TFCheckpoint database, an integrated resource focused on human, mouse and rat dbTFs, providing users access to a comprehensive overview of these proteins. Here, we describe TFCheckpoint 2.0 (https://www.tfcheckpoint.org/index.php), comprising 13 collections of dbTFs, coTFs and GTFs. TFCheckpoint 2.0 provides an easy and versatile cross-referencing system for users to view and download collections that may otherwise be cumbersome to find, compare and retrieve.

Essential transcription factors for induced neuron differentiation

Neurogenins are proneural transcription factors required to specify neuronal identity. Their overexpression in human pluripotent stem cells rapidly produces cortical-like neurons with spiking activity and, because of this, they have been widely adopted for human neuron disease models. However, we do not fully understand the key downstream regulatory effectors responsible for driving neural differentiation. Here, using inducible expression of NEUROG1 and NEUROG2, we identify transcription factors (TFs) required for directed neuronal differentiation by combining expression and chromatin accessibility analyses with a pooled in vitro CRISPR-Cas9 screen targeting all ~1900 TFs in the human genome. The loss of one of these essential TFs (ZBTB18) yields few MAP2-positive neurons. Differentiated ZBTB18-null cells have radically altered gene expression, leading to cytoskeletal defects and stunted neurites and spines. In addition to identifying key downstream TFs for neuronal differentiation, our work develops an integrative multi-omics and TFome-wide perturbation platform to rapidly characterize essential TFs for the differentiation of any human cell type.

Single-cell RNA sequencing unveils Lrg1's role in cerebral ischemia‒reperfusion injury by modulating various cells

BACKGROUND AND PURPOSE

Cerebral ischemia‒reperfusion injury causes significant harm to human health and is a major contributor to stroke-related deaths worldwide. Current treatments are limited, and new, more effective prevention and treatment strategies that target multiple cell components are urgently needed. Leucine-rich alpha-2 glycoprotein 1 (Lrg1) appears to be associated with the progression of cerebral ischemia‒reperfusion injury, but the exact mechanism of it is unknown.

METHODS

Wild-type (WT) and Lrg1 knockout (Lrg1-/-) mice were used to investigate the role of Lrg1 after cerebral ischemia‒reperfusion injury. The effects of Lrg1 knockout on brain infarct volume, blood‒brain barrier permeability, and neurological score (based on 2,3,5-triphenyl tetrazolium chloride, evans blue dye, hematoxylin, and eosin staining) were assessed. Single-cell RNA sequencing (scRNA-seq), immunofluorescence, and microvascular albumin leakage tests were utilized to investigate alterations in various cell components in brain tissue after Lrg1 knockout.

RESULTS

Lrg1 expression was increased in various cell types of brain tissue after cerebral ischemia‒reperfusion injury. Lrg1 knockout reduced cerebral edema and infarct size and improved neurological function after cerebral ischemia‒reperfusion injury. Single-cell RNA sequencing analysis of WT and Lrg1-/- mouse brain tissues after cerebral ischemia‒reperfusion injury revealed that Lrg1 knockout enhances blood‒brain barrier (BBB) by upregulating claudin 11, integrin β5, protocadherin 9, and annexin A2. Lrg1 knockout also promoted an anti-inflammatory and tissue-repairing phenotype in microglia and macrophages while reducing neuron and oligodendrocyte cell death.

CONCLUSIONS

Our results has shown that Lrg1 mediates numerous pathological processes involved in cerebral ischemia‒reperfusion injury by altering the functional states of various cell types, thereby rendering it a promising therapeutic target for cerebral ischemia‒reperfusion injury.

Impaired bone morphogenetic protein (BMP) signaling pathways disrupt decidualization in endometriosis

Endometriosis is linked to increased infertility and pregnancy complications due to defective endometrial decidualization. We hypothesized that identification of altered signaling pathways during decidualization could identify the underlying cause of infertility and pregnancy complications. Our study reveals that transforming growth factor β (TGFβ) pathways are impaired in the endometrium of individuals with endometriosis, leading to defective decidualization. Through detailed transcriptomic analyses, we discovered abnormalities in TGFβ signaling pathways and key regulators, such as SMAD4, in the endometrium of affected individuals. We also observed compromised activity of bone morphogenetic proteins (BMP), a subset of the TGFβ family, that control endometrial receptivity. Using 3-dimensional models of endometrial stromal and epithelial assembloids, we showed that exogenous BMP2 improved decidual marker expression in individuals with endometriosis. Our findings unveil a previously unidentified dysfunction in BMP/SMAD signaling in the endometrium of individuals with endometriosis, explaining decidualization defects and subsequent pregnancy complications in these individuals.

Psychiatric risk gene Transcription Factor 4 (TCF4) regulates the density and connectivity of distinct inhibitory interneuron subtypes

Transcription factor 4 (TCF4) is a basic helix-loop-helix transcription factor that is implicated in a variety of psychiatric disorders including autism spectrum disorder (ASD), major depression, and schizophrenia. Autosomal dominant mutations in TCF4 are causal for a specific ASD called Pitt-Hopkins Syndrome (PTHS). However, our understanding of etiological and pathophysiological mechanisms downstream of TCF4 mutations is incomplete. Single cell sequencing indicates TCF4 is highly expressed in GABAergic interneurons (INs). Here, we performed cell-type specific expression analysis (CSEA) and cellular deconvolution (CD) on bulk RNA sequencing data from 5 different PTHS mouse models. Using CSEA we observed differentially expressed genes (DEGs) were enriched in parvalbumin expressing (PV+) INs and CD predicted a reduction in the PV+ INs population. Therefore, we investigated the role of TCF4 in regulating the development and function of INs in the Tcf4+/tr mouse model of PTHS. In Tcf4+/tr mice, immunohistochemical (IHC) analysis of subtype-specific IN markers and reporter mice identified reductions in PV+, vasoactive intestinal peptide (VIP+), and cortistatin (CST+) expressing INs in the cortex and cholinergic (ChAT+) INs in the striatum, with the somatostatin (SST+) IN population being spared. The reduction of these specific IN populations led to cell-type specific alterations in the balance of excitatory and inhibitory inputs onto PV+ and VIP+ INs and excitatory pyramidal neurons within the cortex. These data indicate TCF4 is a critical regulator of the development of specific subsets of INs and highlight the inhibitory network as an important source of pathophysiology in PTHS.

Paired yeast one-hybrid assays to detect DNA-binding cooperativity and antagonism across transcription factors

Cooperativity and antagonism between transcription factors (TFs) can drastically modify their binding to regulatory DNA elements. While mapping these relationships between TFs is important for understanding their context-specific functions, existing approaches either rely on DNA binding motif predictions, interrogate one TF at a time, or study individual TFs in parallel. Here, we introduce paired yeast one-hybrid (pY1H) assays to detect cooperativity and antagonism across hundreds of TF-pairs at DNA regions of interest. We provide evidence that a wide variety of TFs are subject to modulation by other TFs in a DNA region-specific manner. We also demonstrate that TF-TF relationships are often affected by alternative isoform usage and identify cooperativity and antagonism between human TFs and viral proteins from human papillomaviruses, Epstein-Barr virus, and other viruses. Altogether, pY1H assays provide a broadly applicable framework to study how different functional relationships affect protein occupancy at regulatory DNA regions.

Characterization of allosteric modulators that disrupt androgen receptor co-activator protein-protein interactions to alter transactivation-Drug leads for metastatic castration resistant prostate cancer

Three series of compounds were prioritized from a high content screening campaign that identified molecules that blocked dihydrotestosterone (DHT) induced formation of Androgen Receptor (AR) protein-protein interactions (PPIs) with the Transcriptional Intermediary Factor 2 (TIF2) coactivator and also disrupted preformed AR-TIF2 PPI complexes; the hydrobenzo-oxazepins (S1), thiadiazol-5-piperidine-carboxamides (S2), and phenyl-methyl-indoles (S3). Compounds from these series inhibited AR PPIs with TIF2 and SRC-1, another p160 coactivator, in mammalian 2-hybrid assays and blocked transcriptional activation in reporter assays driven by full length AR or AR-V7 splice variants. Compounds inhibited the growth of five prostate cancer cell lines, with many exhibiting differential cytotoxicity towards AR positive cell lines. Representative compounds from the 3 series substantially reduced both endogenous and DHT-enhanced expression and secretion of the prostate specific antigen (PSA) cancer biomarker in the C4-2 castration resistant prostate cancer (CRPC) cell line. The comparatively weak activities of series compounds in the H3-DHT and/or TIF2 box 3 LXXLL-peptide binding assays to the recombinant ligand binding domain of AR suggest that direct antagonism at the orthosteric ligand binding site or AF-2 surface respectively are unlikely mechanisms of action. Cellular enhanced thermal stability assays (CETSA) indicated that compounds engaged AR and reduced the maximum efficacy and right shifted the EC50 of DHT-enhanced AR thermal stabilization consistent with the effects of negative allosteric modulators. Molecular docking of potent representative hits from each series to AR structures suggest that S1-1 and S2-6 engage a novel binding pocket (BP-1) adjacent to the orthosteric ligand binding site, while S3-11 occupies the AR binding function 3 (BF-3) allosteric pocket. Hit binding poses indicate spaces and residues adjacent to the BP-1 and BF-3 pockets that will be exploited in future medicinal chemistry optimization studies. Small molecule allosteric modulators that prevent/disrupt AR PPIs with coactivators like TIF2 to alter transcriptional activation in the presence of orthosteric agonists might evade the resistance mechanisms to existing prostate cancer drugs and provide novel starting points for medicinal chemistry lead optimization and future development into therapies for metastatic CRPC.

Impaired bone morphogenetic protein signaling pathways disrupt decidualization in endometriosis

It is hypothesized that impaired endometrial decidualization contributes to decreased fertility in individuals with endometriosis. To identify the molecular defects that underpin defective decidualization in endometriosis, we subjected endometrial stromal cells from individuals with or without endometriosis to time course in vitro decidualization with estradiol, progesterone, and 8-bromo-cyclic-AMP (EPC) for 2, 4, 6, or 8 days. Transcriptomic profiling identified differences in key pathways between the two groups, including defective bone morphogenetic protein (BMP)/SMAD4 signaling (ID2, ID3, FST), oxidate stress response (NFE2L2, ALOX15, SLC40A1), and retinoic acid signaling pathways (RARRES, RARB, ALDH1B1). Genome-wide binding analyses identified an altered genomic distribution of SMAD4 and H3K27Ac in the decidualized stromal cells from individuals without endometriosis relative to those with endometriosis, with target genes enriched in pathways related to signaling by transforming growth factor β (TGFβ), neurotrophic tyrosine kinase receptors (NTRK), and nerve growth factor (NGF)-stimulated transcription. We found that direct SMAD1/5/4 target genes control FOXO, PI3K/AKT, and progesterone-mediated signaling in decidualizing cells and that BMP2 supplementation in endometriosis patient-derived assembloids elevated the expression of decidualization markers. In summary, transcriptomic and genome-wide binding analyses of patient-derived endometrial cells and assembloids identified that a functional BMP/SMAD1/5/4 signaling program is crucial for engaging decidualization.

High-throughput functional characterization of combinations of transcriptional activators and repressors

Despite growing knowledge of the functions of individual human transcriptional effector domains, much less is understood about how multiple effector domains within the same protein combine to regulate gene expression. Here, we measure transcriptional activity for 8,400 effector domain combinations by recruiting them to reporter genes in human cells. In our assay, weak and moderate activation domains synergize to drive strong gene expression, whereas combining strong activators often results in weaker activation. In contrast, repressors combine linearly and produce full gene silencing, and repressor domains often overpower activation domains. We use this information to build a synthetic transcription factor whose function can be tuned between repression and activation independent of recruitment to target genes by using a small-molecule drug. Altogether, we outline the basic principles of how effector domains combine to regulate gene expression and demonstrate their value in building precise and flexible synthetic biology tools. A record of this paper's transparent peer review process is included in the supplemental information.

Using i-GONAD for Cell-Type-Specific and Systematic Analysis of Developmental Transcription Factors In Vivo

Transcription factors (TFs) regulate gene expression via direct DNA binding together with cofactors and in chromatin remodeling complexes. Their function is thus regulated in a spatiotemporal and cell-type-specific manner. To analyze the functions of TFs in a cell-type-specific context, genome-wide DNA binding, as well as the identification of interacting proteins, is required. We used i-GONAD (improved genome editing via oviductal nucleic acids delivery) in mice to genetically modify TFs by adding fluorescent reporter and affinity tags that can be exploited for the imaging and enrichment of target cells as well as chromatin immunoprecipitation and pull-down assays. As proof-of-principle, we showed the functional genetic modification of the closely related developmental TFs, Bcl11a and Bcl11b, in defined cell types of newborn mice. i-GONAD is a highly efficient procedure for modifying TF-encoding genes via the integration of small insertions, such as reporter and affinity tags. The novel Bcl11a and Bcl11b mouse lines, described in this study, will be used to improve our understanding of the Bcl11 family's function in neurodevelopment and associated disease.

Zinc Ions Modulate YY1 Activity: Relevance in Carcinogenesis

YY1 is widely recognized as an intrinsically disordered transcription factor that plays a role in development of many cancers. In most cases, its overexpression is correlated with tumor progression and unfavorable patient outcomes. Our latest research focusing on the role of zinc ions in modulating YY1's interaction with DNA demonstrated that zinc enhances the protein's multimeric state and affinity to its operator. In light of these findings, changes in protein concentration appear to be just one element relevant to modulating YY1-dependent processes. Thus, alterations in zinc ion concentration can directly and specifically impact the regulation of gene expression by YY1, in line with reports indicating a correlation between zinc ion levels and advancement of certain tumors. This review concentrates on other potential consequences of YY1 interaction with zinc ions that may act by altering charge distribution, conformational state distribution, or oligomerization to influence its interactions with molecular partners that can disrupt gene expression patterns.

RgnTX: Colocalization analysis of transcriptome elements in the presence of isoform heterogeneity and ambiguity

Colocalization analysis of genomic region sets has been widely adopted to unveil potential functional interactions between corresponding biological attributes, which often serves as the basis for further investigation. A number of methods have been developed for colocalization analysis of genomic elements. However, none of them explicitly considered the transcriptome heterogeneity and isoform ambiguity, making them less appropriate for analyzing transcriptome elements. Here, we developed RgnTX, an R/Bioconductor tool for the colocalization analysis of transcriptome elements with permutation tests. Different from existing approaches, RgnTX directly takes advantage of transcriptome annotation, and offers high flexibility in the null model to simulate realistic transcriptome-wide background, such as the complex alternative splicing patterns. Importantly, it supports the testing of transcriptome elements without clear isoform association, which is often the real scenario due to technical limitations. Proposed package offers a wide selection of pre-defined functions, easy to be utilized by users for visualizing permutation results, calculating shifted z-scores and conducting multiple hypothesis testing under Benjamini-Hochberg correction. Moreover, with synthetic and real datasets, we show that RgnTX novel testing modes return distinct and more significant results compared to existing genome-based methods. We believe RgnTX should make a useful tool to characterize the randomness of the transcriptome, and for conducting statistical association analysis for genomic region sets within the heterogeneous transcriptome. The package now has been accepted by Bioconductor and is freely available at: https://bioconductor.org/packages/RgnTX.

Transcriptomic Effects of Low-Dose Inorganic Arsenic Exposure on Murine Bone Marrow-Derived Macrophages

Both tissue-resident macrophages and monocytes recruited from the bone marrow that transform into tissue-resident cells play critical roles in mediating homeostasis as well as in the pathology of inflammatory diseases. Inorganic arsenic (iAs) is the most common drinking water contaminant worldwide and represents a major public health concern. Several diseases that macrophages have implicated involvement in are caused by iAs exposure, including cardiovascular disease, cancer, and increased risk of infectious disease. Therefore, understanding the effects of iAs exposure on macrophages can help us better grasp the full range of arsenic immunotoxicity and better design therapeutic targets for iAs-induced diseases particularly in exposed populations. In this study, we analyzed the transcriptome of low dose iAs-exposed male and female murine bone marrow-derived macrophages (BMDMs) with either M0, M1, or M2 stimulation. We identified differentially expressed genes by iAs in a sex- and stimulation-dependent manner and used bioinformatics tools to predict protein-protein interactions, transcriptional regulatory networks, and associated biological processes. Overall, our data suggest that M1-stimulated, especially female-derived, BMDMs are most susceptible to iAs exposure. Most notably, we observed significant downregulation of major proinflammatory transcription factors, like IRF8, and its downstream targets, as well as genes encoding proteins involved in pattern recognition and antigen presentation, such as TLR7, TLR8, and H2-D1, potentially providing causal insight regarding arsenic's role in perturbing immune responses to infectious diseases. We also observed significant downregulation of genes involved in processes crucial to coordinating a proinflammatory response including leukocyte migration, differentiation, and cytokine and chemokine production and response. Finally, we discovered that 24 X-linked genes were dysregulated in iAs-exposed female stimulation groups compared to only 3 across the iAs-exposed male stimulation groups. These findings elucidate the potential mechanisms underlying the sex-differential iAs-associated immune-related disease risk.

MYC determines lineage commitment in KRAS-driven primary liver cancer development

BACKGROUND & AIMS

Primary liver cancer (PLC) comprises hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA), two frequent and lethal tumour types that differ regarding their tumour biology and responses to cancer therapies. Liver cells harbour a high degree of cellular plasticity and can give rise to either HCC or iCCA. However, little is known about the cell-intrinsic mechanisms directing an oncogenically transformed liver cell to either HCC or iCCA. The scope of this study was to identify cell-intrinsic factors determining lineage commitment in PLC.

METHODS

Cross-species transcriptomic and epigenetic profiling was applied to murine HCCs and iCCAs and to two human PLC cohorts. Integrative data analysis comprised epigenetic Landscape In Silico deletion Analysis (LISA) of transcriptomic data and Hypergeometric Optimization of Motif EnRichment (HOMER) analysis of chromatin accessibility data. Identified candidate genes were subjected to functional genetic testing in non-germline genetically engineered PLC mouse models (shRNAmir knockdown or overexpression of full-length cDNAs).

RESULTS

Integrative bioinformatic analyses of transcriptomic and epigenetic data pinpointed the Forkhead-family transcription factors FOXA1 and FOXA2 as MYC-dependent determination factors of the HCC lineage. Conversely, the ETS family transcription factor ETS1 was identified as a determinant of the iCCA lineage, which was found to be suppressed by MYC during HCC development. Strikingly, shRNA-mediated suppression of FOXA1 and FOXA2 with concomitant ETS1 expression fully switched HCC to iCCA development in PLC mouse models.

CONCLUSIONS

The herein reported data establish MYC as a key determinant of lineage commitment in PLC and provide a molecular explanation why common liver-damaging risk factors such as alcoholic or non-alcoholic steatohepatitis can lead to either HCC or iCCA.

IMPACT AND IMPLICATIONS

Liver cancer is a major health problem and comprises hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA), two frequent and lethal tumour types that differ regarding their morphology, tumour biology, and responses to cancer therapies. We identified the transcription factor and oncogenic master regulator MYC as a switch between HCC and iCCA development. When MYC levels are high at the time point when a hepatocyte becomes a tumour cell, an HCC is growing out. Conversely, if MYC levels are low at this time point, the result is the outgrowth of an iCCA. Our study provides a molecular explanation why common liver-damaging risk factors such as alcoholic or non-alcoholic steatohepatitis can lead to either HCC or iCCA. Furthermore, our data harbour potential for the development of better PLC therapies.

Combinatorial transcription factor binding encodes cis-regulatory wiring of forebrain GABAergic neurogenesis

Transcription factors (TFs) bind combinatorially to genomic cis-regulatory elements (cREs), orchestrating transcription programs. While studies of chromatin state and chromosomal interactions have revealed dynamic neurodevelopmental cRE landscapes, parallel understanding of the underlying TF binding lags. To elucidate the combinatorial TF-cRE interactions driving mouse basal ganglia development, we integrated ChIP-seq for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and transcriptional state, and transgenic enhancer assays. We identified TF-cREs modules with distinct chromatin features and enhancer activity that have complementary roles driving GABAergic neurogenesis and suppressing other developmental fates. While the majority of distal cREs were bound by one or two TFs, a small proportion were extensively bound, and these enhancers also exhibited exceptional evolutionary conservation, motif density, and complex chromosomal interactions. Our results provide new insights into how modules of combinatorial TF-cRE interactions activate and repress developmental expression programs and demonstrate the value of TF binding data in modeling gene regulatory wiring.

NFAT and NF-κB dynamically co-regulate TCR and CAR signaling responses in human T cells

While it has been established that the responses of T cells to antigens are combinatorially regulated by multiple signaling pathways, it remains elusive what mechanisms cells utilize to quantitatively modulate T cell responses during pathway integration. Here, we show that two key pathways in T cell signaling, calcium/nuclear factor of activated T cells (NFAT) and protein kinase C (PKC)/nuclear factor κB (NF-κB), integrate through a dynamic and combinatorial strategy to fine-tune T cell response genes. At the cis-regulatory level, the two pathways integrate through co-binding of NFAT and NF-κB to immune response genes. Pathway integration is further regulated temporally, where T cell receptor (TCR) and chimeric antigen receptor (CAR) activation signals modulate the temporal relationships between the nuclear localization dynamics of NFAT and NF-κB. Such physical and temporal integrations together contribute to distinct modes of expression modulation for genes. Thus, the temporal relationships between regulators can be modulated to affect their co-targets during immune responses, underscoring the importance of dynamic combinatorial regulation in cellular signaling.

Zinc finger and SCAN domain-containing 18 suppresses the proliferation, self-renewal, and drug resistance of glioblastoma cells

Elucidation of cellular and molecular mechanisms key to glioblastoma growth, self-renewal, survival, and metastasis is important for developing novel therapeutic strategies. In this study, the expression and function of zinc finger and SCAN domain-containing 18 (ZSCAN18) in human glioblastoma cell lines were characterized. Compared with normal astrocytes, ZSCAN18 was significantly down-regulated in all tested glioblastoma cell lines, with the LN-229 cell line having the lowest ZSCAN18 expression. Lentivirus-mediated ZSCAN18 overexpression suppressed glioblastoma cell proliferation, sphere formation, and SOX2 and OCT4 expression, implying the negative role of ZSCAN18 in glioblastoma development. ZSCAN18 overexpression enhanced the sensitivity of glioblastoma cells to Temozolomide. The glioblastoma implantation model showed a consistent inhibitory effect of ZSCAN18 on the proliferation and self-renewal of glioblastoma cells in vivo. Notably, ZSCAN18 overexpression resulted in the down-regulation of glioma-associated oncogene homolog 1 (GLI1) which is the terminal component of the Hedgehog signaling. Lentivirus-mediated GLI1 overexpression restored the proliferation and promoted the resistance of glioblastoma cells to Temozolomide. However, GLI1 overexpression did not affect the self-renewal of ZSCAN18-overexpressing glioblastoma cells. Taken together, this research uncovers the role of ZSCAN18 in regulating glioblastoma cell growth and maintenance. ZSCAN18 could be a potential glioblastoma biomarker.

Complementary strategies for directing in vivo transcription factor binding through DNA binding domains and intrinsically disordered regions

DNA binding domains (DBDs) of transcription factors (TFs) recognize DNA sequence motifs that are highly abundant in genomes. Within cells, TFs bind a subset of motif-containing sites as directed by either their DBDs or DBD-external (nonDBD) sequences. To define the relative roles of DBDs and nonDBDs in directing binding preferences, we compared the genome-wide binding of 48 (∼30%) budding yeast TFs with their DBD-only, nonDBD-truncated, and nonDBD-only mutants. With a few exceptions, binding locations differed between DBDs and TFs, resulting from the cumulative action of multiple determinants mapped mostly to disordered nonDBD regions. Furthermore, TFs' preferences for promoters of the fuzzy nucleosome architecture were lost in DBD-only mutants, whose binding spread across promoters, implicating nonDBDs' preferences in this hallmark of budding yeast regulatory design. We conclude that DBDs and nonDBDs employ complementary DNA-targeting strategies, whose balance defines TF binding specificity along genomes.

Transcriptomic analysis of the Non-Obstructive Azoospermia (NOA) to address gene expression regulation in human testis

There is a need to understand the molecular basis of testes under Non-Obstructive Azoospermia (NOA), a state of failed spermatogenesis. There has been a lack of attention to the transcriptome at the level of alternatively spliced mRNAs (iso-mRNAs) and the mechanism of gene expression regulation. Hence, we aimed to establish a reliable iso-mRNA profile of NOA-testes, and explore molecular mechanisms - especially those related to gene expression regulation. We sequenced mRNAs from testicular samples of donors with complete spermatogenesis (control samples) and a failure of spermatogenesis (NOA samples). We identified differentially expressed genes and their iso-mRNAs via standard NGS data analyses. We then listed these iso-mRNAs hierarchically based on the extent of consistency of differential quantities across samples and groups, and validated the lists via RT-qPCRs (for 80 iso-mRNAs). In addition, we performed extensive bioinformatic analysis of the splicing features, domains, interactions, and functions of differentially expressed genes and iso-mRNAs. Many top-ranking down-regulated genes and iso-mRNAs, i.e., those down-regulated more consistently across the NOA samples, are associated with mitosis, replication, meiosis, cilium, RNA regulation, and post-translational modifications such as ubiquitination and phosphorylation. Most down-regulated iso-mRNAs correspond to full-length proteins that include all expected domains. The predominance of alternative promoters and termination sites in these iso-mRNAs indicate their gene expression regulation via promoters and UTRs. We compiled a new, comprehensive list of human transcription factors (TFs) and used it to identify TF-'TF gene' interactions with potential significance in down-regulating genes under the NOA condition. The results indicate that RAD51 suppression by HSF4 prevents SP1-activation, and SP1, in turn, could regulate multiple TF genes. This potential regulatory axis and other TF interactions identified in this study could explain the down-regulation of multiple genes in NOA-testes. Such molecular interactions may also have key regulatory roles during normal human spermatogenesis.

Growing Directed Acyclic Graphs: Optimization Functions for Pathway Reconstruction Algorithms

A major challenge in molecular systems biology is to understand how proteins work to transmit external signals to changes in gene expression. Computationally reconstructing these signaling pathways from protein interaction networks can help understand what is missing from existing pathway databases. We formulate a new pathway reconstruction problem, one that iteratively grows directed acyclic graphs (DAGs) from a set of starting proteins in a protein interaction network. We present an algorithm that provably returns the optimal DAGs for two different cost functions and evaluate the pathway reconstructions when applied to six diverse signaling pathways from the NetPath database. The optimal DAGs outperform an existing k-shortest paths method for pathway reconstruction, and the new reconstructions are enriched for different biological processes. Growing DAGs is a promising step toward reconstructing pathways that provably optimize a specific cost function.

A stay of execution: ATF4 regulation and potential outcomes for the integrated stress response

ATF4 is a cellular stress induced bZIP transcription factor that is a hallmark effector of the integrated stress response. The integrated stress response is triggered by phosphorylation of the alpha subunit of the eukaryotic initiation factor 2 complex that can be carried out by the cellular stress responsive kinases; GCN2, PERK, PKR, and HRI. eIF2α phosphorylation downregulates mRNA translation initiation en masse, however ATF4 translation is upregulated. The integrated stress response can output two contradicting outcomes in cells; pro-survival or apoptosis. The mechanism for choice between these outcomes is unknown, however combinations of ATF4 heterodimerisation partners and post-translational modifications have been linked to this regulation. This semi-systematic review article covers ATF4 target genes, heterodimerisation partners and post-translational modifications. Together, this review aims to be a useful resource to elucidate the mechanisms controlling the effects of the integrated stress response. Additional putative roles of the ATF4 protein in cell division and synaptic plasticity are outlined.

Effects of Positive Fighting Experience and Its Subsequent Deprivation on the Expression Profile of Mouse Hippocampal Genes Associated with Neurogenesis

The hippocampus is known as the brain region implicated in visuospatial processes and processes associated with learning and short- and long-term memory. An important functional characteristic of the hippocampus is lifelong neurogenesis. A decrease or increase in adult hippocampal neurogenesis is associated with a wide range of neurological diseases. We have previously shown that in adult male mice with a chronic positive fighting experience in daily agonistic interactions, there is an increase in the proliferation of progenitor neurons and the production of young neurons in the dentate gyrus (in hippocampus), and these neurogenesis parameters remain modified during 2 weeks of deprivation of further fights. The aim of the present work was to identify hippocampal genes associated with neurogenesis and involved in the formation of behavioral features in mice with the chronic experience of wins in aggressive confrontations, as well as during the subsequent 2-week deprivation of agonistic interactions. Hippocampal gene expression profiles were compared among three groups of adult male mice: chronically winning for 20 days in the agonistic interactions, chronically victorious for 20 days followed by the 2-week deprivation of fights, and intact (control) mice. Neurogenesis-associated genes were identified whose transcription levels changed during the social confrontations and in the subsequent period of deprivation of fights. In the experimental males, some of these genes are associated with behavioral traits, including abnormal aggression-related behavior, an abnormal anxiety-related response, and others. Two genes encoding transcription factors (Nr1d1 and Fmr1) were likely to contribute the most to the between-group differences. It can be concluded that the chronic experience of wins in agonistic interactions alters hippocampal levels of transcription of multiple genes in adult male mice. The transcriptome changes get reversed only partially after the 2-week period of deprivation of fights. The identified differentially expressed genes associated with neurogenesis and involved in the control of a behavior/neurological phenotype can be used in further studies to identify targets for therapeutic correction of the neurological disturbances that develop in winners under the conditions of chronic social confrontations.

Discovering Entities Similarities in Biological Networks Using a Hybrid Immune Algorithm

Disease phenotypes are generally caused by the failure of gene modules which often have similar biological roles. Through the study of biological networks, it is possible to identify the intrinsic structure of molecular interactions in order to identify the so-called “disease modules”. Community detection is an interesting and valuable approach to discovering the structure of the community in a complex network, revealing the internal organization of the nodes, and has become a leading research topic in the analysis of complex networks. This work investigates the link between biological modules and network communities in test-case biological networks that are commonly used as a reference point and which include Protein–Protein Interaction Networks, Metabolic Networks and Transcriptional Regulation Networks. In order to identify small and structurally well-defined communities in the biological context, a hybrid immune metaheuristic algorithm Hybrid-IA is proposed and compared with several metaheuristics, hyper-heuristics, and the well-known greedy algorithm Louvain, with respect to modularity maximization. Considering the limitation of modularity optimization, which can fail to identify smaller communities, the reliability of Hybrid-IA was also analyzed with respect to three well-known sensitivity analysis measures (NMI, ARI and NVI) that assess how similar the detected communities are to real ones. By inspecting all outcomes and the performed comparisons, we will see that on one hand Hybrid-IA finds slightly lower modularity values than Louvain, but outperforms all other metaheuristics, while on the other hand, it can detect communities more similar to the real ones when compared to those detected by Louvain.

A transcription factor atlas of directed differentiation

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

Targeting Transcription Factors in Cancer: From "Undruggable" to "Druggable"

Deregulation of transcription factors is critical to hallmarks of cancer. Genetic mutations, gene fusions, amplifications or deletions, epigenetic alternations, and aberrant post-transcriptional modification of transcription factors are involved in the regulation of various stages of carcinogenesis, including cancer initiation, progression, and metastasis. Thus, targeting the dysfunctional transcription factors may lead to new cancer therapeutic strategies. However, transcription factors are conventionally considered as "undruggable." Here, we summarize the recent progresses in understanding the regulation of transcription factors in cancers and strategies to target transcription factors and co-factors for preclinical and clinical drug development, particularly focusing on c-Myc, YAP/TAZ, and β-catenin due to their significance and interplays in cancer.

Detecting clusters of transcription factors based on a nonhomogeneous poisson process model

BACKGROUND

Rapidly growing genome-wide ChIP-seq data have provided unprecedented opportunities to explore transcription factor (TF) binding under various cellular conditions. Despite the rich resources, development of analytical methods for studying the interaction among TFs in gene regulation still lags behind.

RESULTS

In order to address cooperative TF binding and detect TF clusters with coordinative functions, we have developed novel computational methods based on clustering the sample paths of nonhomogeneous Poisson processes. Simulation studies demonstrated the capability of these methods to accurately detect TF clusters and uncover the hierarchy of TF interactions. A further application to the multiple-TF ChIP-seq data in mouse embryonic stem cells (ESCs) showed that our methods identified the cluster of core ESC regulators reported in the literature and provided new insights on functional implications of transcrisptional regulatory modules.

CONCLUSIONS

Effective analytical tools are essential for studying protein-DNA relations. Information derived from this research will help us better understand the orchestration of transcription factors in gene regulation processes.

Integrated proteomic and transcriptomic landscape of macrophages in mouse tissues

Macrophages are involved in tissue homeostasis and are critical for innate immune responses, yet distinct macrophage populations in different tissues exhibit diverse gene expression patterns and biological processes. While tissue-specific macrophage epigenomic and transcriptomic profiles have been reported, proteomes of different macrophage populations remain poorly characterized. Here we use mass spectrometry and bulk RNA sequencing to assess the proteomic and transcriptomic patterns, respectively, of 10 primary macrophage populations from seven mouse tissues, bone marrow-derived macrophages and the cell line RAW264.7. The results show distinct proteomic landscape and protein copy numbers between tissue-resident and recruited macrophages. Construction of a hierarchical regulatory network finds cell-type-specific transcription factors of macrophages serving as hubs for denoting tissue and functional identity of individual macrophage subsets. Finally, Il18 is validated to be essential in distinguishing molecular signatures and cellular function features between tissue-resident and recruited macrophages in the lung and liver. In summary, these deposited datasets and our open proteome server ( http://macrophage.mouseprotein.cn ) integrating all information will provide a valuable resource for future functional and mechanistic studies of mouse macrophages.

scHumanNet: a single-cell network analysis platform for the study of cell-type specificity of disease genes

A major challenge in single-cell biology is identifying cell-type-specific gene functions, which may substantially improve precision medicine. Differential expression analysis of genes is a popular, yet insufficient approach, and complementary methods that associate function with cell type are required. Here, we describe scHumanNet (https://github.com/netbiolab/scHumanNet), a single-cell network analysis platform for resolving cellular heterogeneity across gene functions in humans. Based on cell-type-specific gene networks (CGNs) constructed under the guidance of the HumanNet reference interactome, scHumanNet displayed higher functional relevance to the cellular context than CGNs built by other methods on single-cell transcriptome data. Cellular deconvolution of gene signatures based on network compactness across cell types revealed breast cancer prognostic markers associated with T cells. scHumanNet could also prioritize genes associated with particular cell types using CGN centrality and identified the differential hubness of CGNs between disease and healthy conditions. We demonstrated the usefulness of scHumanNet by uncovering T-cell-specific functional effects of GITR, a prognostic gene for breast cancer, and functional defects in autism spectrum disorder genes specific for inhibitory neurons. These results suggest that scHumanNet will advance our understanding of cell-type specificity across human disease genes.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

Transcription factor NFE2L3 promotes the proliferation of esophageal squamous cell carcinoma cells and causes radiotherapy resistance by regulating IL-6

OBJECTIVE

To scrutinize the impact of overexpression and interference of NFE2L3 on radiosensitivity of esophageal squamous cell carcinoma cells (ESCC) and its downstream mechanism and to assess whether NFE2L3 expression alters in vivo radiosensitivity of ESCC by developing a subcutaneous tumor model in mice.

METHODS

Through RNA-Seq, we compared the differentially expressed genes between the ECA-109R cell line and its parent ECA-109 cell line. The differentially expressed genes were selected and verified by qRT-PCR. Transfection of ESCC cell lines with NFE2L3 inhibitor or mimic lentivirus constructs was done to study the activity of NFE2L3. To assess the effect of NFE2L3 on cellular growth and proliferation, clonogenic survival assay, EdU incorporation assay, and CCK-8 assay were done after irradiation. To probe how many irradiated DNA double-strand breaks were produced, the corresponding intensity of γ-H2AX foci were detected by immunofluorescence. Apoptotic cells were assayed by flow cytometry assay after irradiation; To investigate the downstream genes of NFE2L3, we knocked NFE2L3, and RNA-Seq was used to find out the downstream genes. qRT-PCR and western blot ensued to score associated protein profiles. The in vivo ESCC cell radiosensitivity was scrutinized by nude mouse xenograft models.

RESULTS

The differential genes between ECA-109R cells and its parent ECA-109 cells were compared by qRT-PCR to unveil a significant increase in NFE2L3 expression. Functional analysis indicated that NFE2L3 increased radioresistance in ESCC cells. Then, through high-throughput sequencing and bioinformatics analysis, IL-6 was found to be a hub gene that played a role downstream of NFE2L3 and was verified by qRT-PCR, western blot, and double luciferase reporter gene experiment. NFE2L3 could regulate ESCC cell radiosensitivity via the IL-6-STAT3 signaling pathway, and downregulation of IL-6 expression could reverse the effects of highly expressed NFE2L3. In vivo tumor xenograft experiments confirmed that NFE2L3 affects the sensitivity to radiation therapy.

CONCLUSION

NFE2L3 can affect the radiosensitivity of ESCC cells through IL-6 transcription and IL-6/STAT3 signaling pathway. This makes NFE2L3 a putative target to regulate ESCC cell radiosensitivity.

A map of cis-regulatory modules and constituent transcription factor binding sites in 80% of the mouse genome

BACKGROUND

Mouse is probably the most important model organism to study mammal biology and human diseases. A better understanding of the mouse genome will help understand the human genome, biology and diseases. However, despite the recent progress, the characterization of the regulatory sequences in the mouse genome is still far from complete, limiting its use to understand the regulatory sequences in the human genome.

RESULTS

Here, by integrating binding peaks in ~ 9,000 transcription factor (TF) ChIP-seq datasets that cover 79.9% of the mouse mappable genome using an efficient pipeline, we were able to partition these binding peak-covered genome regions into a cis-regulatory module (CRM) candidate (CRMC) set and a non-CRMC set. The CRMCs contain 912,197 putative CRMs and 38,554,729 TF binding sites (TFBSs) islands, covering 55.5% and 24.4% of the mappable genome, respectively. The CRMCs tend to be under strong evolutionary constraints, indicating that they are likely cis-regulatory; while the non-CRMCs are largely selectively neutral, indicating that they are unlikely cis-regulatory. Based on evolutionary profiles of the genome positions, we further estimated that 63.8% and 27.4% of the mouse genome might code for CRMs and TFBSs, respectively.

CONCLUSIONS

Validation using experimental data suggests that at least most of the CRMCs are authentic. Thus, this unprecedentedly comprehensive map of CRMs and TFBSs can be a good resource to guide experimental studies of regulatory genomes in mice and humans.

The Transcription Factor YY-1 Is an Essential Regulator of T Follicular Helper Cell Differentiation

T follicular helper (TFH) cells are a specialized subset of CD4 T cells that deliver critical help signals to B cells for the production of high-affinity Abs. Understanding the genetic program regulating TFH differentiation is critical if one wants to manipulate TFH cells during vaccination. A large number of transcription factor (TFs) involved in the regulation of TFH differentiation have been characterized. However, there are likely additional unknown TFs required for this process. To identify new TFs, we screened a large short hairpin RNA library targeting 353 TFs in mice using an in vivo RNA interference screen. Yin Yang 1 (YY-1) was identified as a novel positive regulator of TFH differentiation. Ablation of YY-1 severely impaired TFH differentiation following acute viral infection and protein immunization. We found that the zinc fingers of YY-1 are critical to support TFH differentiation. Thus, we discovered a novel TF involved in the regulation of TFH cells.

Identification of Hypothalamic Long Noncoding RNAs Associated with Hypertension and the Behavior/Neurological Phenotype of Hypertensive ISIAH Rats

Long noncoding RNAs (lncRNAs) play an important role in the control of many physiological and pathophysiological processes, including the development of hypertension and other cardiovascular diseases. Nonetheless, the understanding of the regulatory function of many lncRNAs is still incomplete. This work is a continuation of our earlier study on the sequencing of hypothalamic transcriptomes of hypertensive ISIAH rats and control normotensive WAG rats. It aims to identify lncRNAs that may be involved in the formation of the hypertensive state and the associated behavioral features of ISIAH rats. Interstrain differences in the expression of seven lncRNAs were validated by quantitative PCR. Differential hypothalamic expression of lncRNAs LOC100910237 and RGD1562890 between hypertensive and normotensive rats was shown for the first time. Expression of four lncRNAs (Snhg4, LOC100910237, RGD1562890, and Tnxa-ps1) correlated with transcription levels of many hypothalamic genes differentially expressed between ISIAH and WAG rats (DEGs), including genes associated with the behavior/neurological phenotype and hypertension. After functional annotation of these DEGs, it was concluded that lncRNAs Snhg4, LOC100910237, RGD1562890, and Tnxa-ps1 may be involved in the hypothalamic processes related to immune-system functioning and in the response to various exogenous and endogenous factors, including hormonal stimuli. Based on the functional enrichment analysis of the networks, an association of lncRNAs LOC100910237 and Tnxa-ps1 with retinol metabolism and an association of lncRNAs RGD1562890 and Tnxa-ps1 with type 1 diabetes mellitus are proposed for the first time. Based on a discussion, it is hypothesized that previously functionally uncharacterized lncRNA LOC100910237 is implicated in the regulation of hypothalamic processes associated with dopaminergic synaptic signaling, which may contribute to the formation of the behavioral/neurological phenotype and hypertensive state of ISIAH rats.

Exploration of Target Spaces in the Human Genome for Protein and Peptide Drugs

After decades of development, protein and peptide drugs have now grown into a major drug class in the marketplace. Target identification and validation are crucial for the discovery of protein and peptide drugs, and bioinformatics prediction of targets based on the characteristics of known target proteins will help improve the efficiency and success rate of target selection. However, owing to the developmental history in the pharmaceutical industry, previous systematic exploration of the target spaces has mainly focused on traditional small-molecule drugs, while studies related to protein and peptide drugs are lacking. Here, we systematically explore the target spaces in the human genome specifically for protein and peptide drugs. Compared with other proteins, both successful protein and peptide drug targets have many special characteristics, and are also significantly different from those of small-molecule drugs in many aspects. Based on these features, we develop separate effective genome-wide target prediction models for protein and peptide drugs. Finally, a user-friendly web server, Predictor Of Protein and PeptIde drugs' therapeutic Targets (POPPIT) (http://poppit.ncpsb.org.cn/), is established, which provides not only target prediction specifically for protein and peptide drugs but also abundant annotations for predicted targets.

Temporally divergent regulatory mechanisms govern neuronal diversification and maturation in the mouse and marmoset neocortex

Mammalian neocortical neurons span one of the most diverse cell type spectra of any tissue. Cortical neurons are born during embryonic development, and their maturation extends into postnatal life. The regulatory strategies underlying progressive neuronal development and maturation remain unclear. Here we present an integrated single-cell epigenomic and transcriptional analysis of individual mouse and marmoset cortical neuron classes, spanning both early postmitotic stages of identity acquisition and later stages of neuronal plasticity and circuit integration. We found that, in both species, the regulatory strategies controlling early and late stages of pan-neuronal development diverge. Early postmitotic neurons use more widely shared and evolutionarily conserved molecular regulatory programs. In contrast, programs active during later neuronal maturation are more brain- and neuron-specific and more evolutionarily divergent. Our work uncovers a temporal shift in regulatory choices during neuronal diversification and maturation in both mice and marmosets, which likely reflects unique evolutionary constraints on distinct events of neuronal development in the neocortex.

Foxh1 engages in chromatin regulation revealed by protein interactome analyses

Early embryonic cell fates are specified through coordinated integration of transcription factor activities and epigenetic states of the genome. Foxh1 is a key maternal transcription factor controlling the mesendodermal gene regulatory program. Proteomic interactome analyses using FOXH1 as a bait in mouse embryonic stem cells revealed that FOXH1 interacts with PRC2 subunits and HDAC1. Foxh1 physically interacts with Hdac1, and confers transcriptional repression of mesendodermal genes in Xenopus ectoderm. Our findings reveal a central role of Foxh1 in coordinating the chromatin states of the Xenopus embryonic genome.

Dynamic regulatory module networks for inference of cell type-specific transcriptional networks

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type–specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type–specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type–specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

TGFβ-induced FOXS1 controls epithelial-mesenchymal transition and predicts a poor prognosis in liver cancer

Transforming growth factor beta (TGF-β) plays a key role in tumor progression, notably as a potent inducer of epithelial-mesenchymal transition (EMT). However, all of the molecular effectors driving TGFβ-induced EMT are not fully characterized. Here, we report that forkhead box S1 (FOXS1) is a SMAD (mothers against decapentaplegic)-dependent TGFβ-induced transcription factor, which regulates the expression of genes required for the initial steps of EMT (e.g., snail family transcription repressor 1) and to maintain a mesenchymal phenotype in hepatocellular carcinoma (HCC) cells. In human HCC, we report that FOXS1 is a biomarker of poorly differentiated and aggressive tumor subtypes. Importantly, FOXS1 expression level and activity are associated with a poor prognosis (e.g., reduced patient survival), not only in HCC but also in colon, stomach, and kidney cancers. Conclusion: FOXS1 constitutes a clinically relevant biomarker for tumors in which the pro-metastatic arm of TGF-β is active (i.e., patients who may benefit from targeted therapies using inhibitors of the TGF-β pathway).

Single-cell analysis of transcription factor regulatory networks reveals molecular basis for subtype-specific dysregulation in acute myeloid leukemia

Highly heterogeneous acute myeloid leukemia (AML) exhibits dysregulated transcriptional programs. Transcription factor (TF) regulatory networks underlying AML subtypes have not been elucidated at single-cell resolution. Here, we comprehensively mapped malignancy-related TFs activated in different AML subtypes by analyzing single-cell RNA sequencing data from AMLs and healthy donors. We first identified six modules of regulatory networks which were prevalently dysregulated in all AML patients. AML subtypes featured with different malignant cellular composition possessed subtype-specific regulatory TFs associated with differentiation suppression or immune modulation. At last, we validated that ERF was crucial for the development of hematopoietic stem/progenitor cells by performing loss- and gain-of-function experiments in zebrafish embryos. Collectively, our work thoroughly documents an abnormal spectrum of transcriptional regulatory networks in AML and reveals subtype-specific dysregulation basis, which provides a prospective view to AML pathogenesis and potential targets for both diagnosis and therapy.

The BTB transcription factors ZBTB11 and ZFP131 maintain pluripotency by repressing pro-differentiation genes

In pluripotent cells, a delicate activation-repression balance maintains pro-differentiation genes ready for rapid activation. The identity of transcription factors (TFs) that specifically repress pro-differentiation genes remains obscure. By targeting ∼1,700 TFs with CRISPR loss-of-function screen, we found that ZBTB11 and ZFP131 are required for embryonic stem cell (ESC) pluripotency. ESCs without ZBTB11 or ZFP131 lose colony morphology, reduce proliferation rate, and upregulate transcription of genes associated with three germ layers. ZBTB11 and ZFP131 bind proximally to pro-differentiation genes. ZBTB11 or ZFP131 loss leads to an increase in H3K4me3, negative elongation factor (NELF) complex release, and concomitant transcription at associated genes. Together, our results suggest that ZBTB11 and ZFP131 maintain pluripotency by preventing premature expression of pro-differentiation genes and present a generalizable framework to maintain cellular potency.

The Fra-1/AP-1 Oncoprotein: From the "Undruggable" Transcription Factor to Therapeutic Targeting

The genetic and epigenetic changes affecting transcription factors, coactivators, and chromatin modifiers are key determinants of the hallmarks of cancer. The acquired dependence on oncogenic transcriptional regulators, representing a major determinant of cancer cell vulnerability, points to transcription factors as ideal therapeutic targets. However, given the unavailability of catalytic activities or binding pockets for small-molecule inhibitors, transcription factors are generally regarded as undruggable proteins. Among components of the AP-1 complex, the FOS-family transcription factor Fra-1, encoded by FOSL1, has emerged as a prominent therapeutic target. Fra-1 is overexpressed in most solid tumors, in response to the BRAF-MAPK, Wnt-beta-catenin, Hippo-YAP, IL-6-Stat3, and other major oncogenic pathways. In vitro functional analyses, validated in onco-mouse models and corroborated by prognostic correlations, show that Fra-1-containing dimers control tumor growth and disease progression. Fra-1 participates in key mechanisms of cancer cell invasion, Epithelial-to-Mesenchymal Transition, and metastatic spreading, by driving the expression of EMT-inducing transcription factors, cytokines, and microRNAs. Here we survey various strategies aimed at inhibiting tumor growth, metastatic dissemination, and drug resistance by interfering with Fra-1 expression, stability, and transcriptional activity. We summarize several tools aimed at the design and tumor-specific delivery of Fra-1/AP-1-specific drugs. Along with RNA-based therapeutics targeting the FOSL1 gene, its mRNA, or cognate regulatory circRNAs, we will examine the exploitation of blocking peptides, small molecule inhibitors, and innovative Fra-1 protein degraders. We also consider the possible caveats concerning Fra-1 inhibition in specific therapeutic contexts. Finally, we discuss a recent suicide gene therapy-based approach, aimed at selectively killing the Fra-1-overexpressing neoplastic cells.