Transcriptional Regulatory Circuits Research Articles

Optimal CXCR5 Expression during Tfh Maturation Involves the Bhlhe40-Pou2af1 Axis.

The pair of transcription factors Bcl6-Blimp1 is well-known for follicular T helper (Tfh) cell fate determination, however, the mechanism(s) for Bcl6-independent regulation of CXCR5 during Tfh migration into germinal center (GC) is still unclear. In this study, we uncovered another pair of transcription factors, Bhlhe40-Pou2af1, that regulates CXCR5 expression. Pou2af1 was specifically expressed in Tfh cells whereas Bhlhe40 expression was found high in non-Tfh cells. Pou2af1 promoted Tfh formation and migration into GC by upregulating CXCR5 but not Bcl6, while Bhlhe40 repressed this process by inhibiting Pou2af1 expression. RNA-Seq analysis of antigen-specific Tfh cells generated in vivo confirmed the role of Bhlhe40-Pou2af1 axis in regulating optimal CXCR5 expression. Thus, the regulation of CXCR5 expression and migration of Tfh cells into GC involves a transcriptional regulatory circuit consisting of Bhlhe40 and Pou2af1, which operates independent of the Bcl6-Blimp1 circuit that determines the Tfh cell fate.

Abstract 6270: Anti-metastatic effects of estrogen receptor β in inflammatory breast cancer

Abstract Inflammatory breast cancer (IBC) is a rare but highly aggressive form of breast cancer with a 5-year survival rate of less than 50%. Although 20-30% of women with IBC have distant metastasis at the time of diagnosis, the biology underlying the aggressive behavior of the disease is poorly understood and molecular markers that distinguish IBC from non-IBC have not yet been identified. In response to the urgent need of discovering new targeted therapies in IBC, we have studied the second estrogen receptor (ERβ) since prior work from our lab correlated higher tumor expression of ERβ in samples from IBC patients with longer metastasis free survival. We showed that ERβ exerted anti-metastatic activity in IBC by inhibiting actin-based cell migration through the suppression of RhoC GTPase signaling. However, since estrogen receptors elicit their effects by inducing global cellular alterations, we have now adopted an integrated genomic approach to fully understand the anti-metastatic effects of ERβ in IBC and delineate the underlying transcriptional regulatory circuits. We found ERβ binding sites in IBC cells and combined this information with gene expression data to identify downstream direct target genes and associated pathways. Our findings suggest the regulation of tumor micro-environment and metabolism as the primary targets of ERβ and we are currently working to pinpoint the crucial regulatory components in these functions. Our findings thus offer an opportunity to better understand the mechanism of ERβ action in aggressive breast cancer, serving as a critical step in validating its tumor repressive role. Citation Format: Harika Nagandla, Spyros Tastsoglou, Kim Cuong Cap, Aireana Phillips, Wei Qian, Jianying Zhou, Jenny Chang, Savitri Krishnamurthy, Naoto Ueno, Artemis Hatzigeorgiou, Christoforos Thomas. Anti-metastatic effects of estrogen receptor β in inflammatory breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6270.

Engineering and application of LacI mutants with stringent expressions.

Optimal transcriptional regulatory circuits are expected to exhibit stringent control, maintaining silence in the absence of inducers while exhibiting a broad induction dynamic range upon the addition of effectors. In the Plac /LacI pair, the promoter of the lac operon in Escherichia coli is characterized by its leakiness, attributed to the moderate affinity of LacI for its operator target. In response to this limitation, the LacI regulatory protein underwent engineering to enhance its regulatory properties. The M7 mutant, carrying I79T and N246S mutations, resulted in the lac promoter displaying approximately 95% less leaky expression and a broader induction dynamic range compared to the wild-type LacI. An in-depth analysis of each mutation revealed distinct regulatory profiles. In contrast to the wild-type LacI, the M7 mutant exhibited a tighter binding to the operator sequence, as evidenced by surface plasmon resonance studies. Leveraging the capabilities of the M7 mutant, a high-value sugar biosensor was constructed. This biosensor facilitated the selection of mutant galactosidases with approximately a seven-fold improvement in specific activity for transgalactosylation. Consequently, this advancement enabled enhanced biosynthesis of galacto-oligosaccharides (GOS).

Non-coding RNAs and angiogenesis in cardiovascular diseases: a comprehensive review.

Non-coding RNAs (ncRNAs) have key roles in the etiology of many illnesses, including heart failure, myocardial infarction, stroke, and in physiological processes like angiogenesis. In transcriptional regulatory circuits that control heart growth, signaling, and stress response, as well as remodeling in cardiac disease, ncRNAs have become important players. Studies on ncRNAs and cardiovascular disease have made great progress recently. Here, we go through the functions of non-coding RNAs (ncRNAs) like circular RNAs (circRNAs), and microRNAs (miRNAs) as well as long non-coding RNAs (lncRNAs) in modulating cardiovascular disorders.

Regulus infers signed regulatory relations from few samples' information using discretization and likelihood constraints.

Transcriptional regulation is performed by transcription factors (TF) binding to DNA in context-dependent regulatory regions and determines the activation or inhibition of gene expression. Current methods of transcriptional regulatory circuits inference, based on one or all of TF, regions and genes activity measurements require a large number of samples for ranking the candidate TF-gene regulation relations and rarely predict whether they are activations or inhibitions. We hypothesize that transcriptional regulatory circuits can be inferred from fewer samples by (1) fully integrating information on TF binding, gene expression and regulatory regions accessibility, (2) reducing data complexity and (3) using biology-based likelihood constraints to determine the global consistency between a candidate TF-gene relation and patterns of genes expressions and region activations, as well as qualify regulations as activations or inhibitions. We introduce Regulus, a method which computes TF-gene relations from gene expressions, regulatory region activities and TF binding sites data, together with the genomic locations of all entities. After aggregating gene expressions and region activities into patterns, data are integrated into a RDF (Resource Description Framework) endpoint. A dedicated SPARQL (SPARQL Protocol and RDF Query Language) query retrieves all potential relations between expressed TF and genes involving active regulatory regions. These TF-region-gene relations are then filtered using biological likelihood constraints allowing to qualify them as activation or inhibition. Regulus provides signed relations consistent with public databases and, when applied to biological data, identifies both known and potential new regulators. Regulus is devoted to context-specific transcriptional circuits inference in human settings where samples are scarce and cell populations are closely related, using discretization into patterns and likelihood reasoning to decipher the most robust regulatory relations.

Landscape and perturbation of enhancer-driven core transcription regulatory circuits in cancer

Landscape and perturbation of enhancer-driven core transcription regulatory circuits in cancer

Molecular Evidence for Olfactory Neuroblastoma as a Tumor of Malignant Globose Basal Cells

Olfactory neuroblastoma (ONB, esthesioneuroblastoma) is a sinonasal cancer with an underdeveloped diagnostic toolkit, and is the subject of many incidents of tumor misclassification throughout the literature. Despite its name, connections between the cancer and normal cells of the olfactory epithelium have not been systematically explored and markers of olfactory epithelial cell types are not deployed in clinical practice. Here, we utilize an integrated human-mouse single-cell atlas of the nasal mucosa, including the olfactory epithelium, to identify transcriptomic programs that link ONB to a specific population of stem/progenitor cells known as olfactory epithelial globose basal cells (GBCs). Expression of a GBC transcription factor NEUROD1 distinguishes both low- and high-grade ONB from sinonasal undifferentiated carcinoma, a potential histologic mimic with a distinctly unfavorable prognosis. Furthermore, we identify a reproducible subpopulation of highly proliferative ONB cells expressing the GBC stemness marker EZH2, suggesting that EZH2 inhibition may play a role in the targeted treatment of ONB. Finally, we study the cellular states comprising ONB parenchyma using single-cell transcriptomics and identify evidence of a conserved GBC transcriptional regulatory circuit that governs divergent neuronal-versus-sustentacular differentiation. These results link ONB to a specific cell type for the first time and identify conserved developmental pathways within ONB that inform diagnostic, prognostic, and mechanistic investigation.

A quantitative evaluation of topological motifs and their coupling in gene circuit state distributions

One of the major challenges in biology is to understand how gene interactions collaborate to determine overall functions of biological systems. Here, we present a new computational framework that enables systematic, high-throughput, and quantitative evaluation of how small transcriptional regulatory circuit motifs, and their coupling, contribute to functions of a dynamical biological system. We illustrate how this approach can be applied to identify four-node gene circuits, circuit motifs, and motif coupling responsible for various gene expression state distributions, including those derived from single-cell RNA sequencing data. We also identify seven major classes of four-node circuits from clustering analysis of state distributions. The method is applied to establish phenomenological models of gene circuits driving human neuron differentiation, revealing important biologically relevant regulatory interactions. Our study will shed light on a better understanding of gene regulatory mechanisms in creating and maintaining cellular states.

OsMADS58 Stabilizes Gene Regulatory Circuits during Rice Stamen Development.

Rice (Oryza sativa) OsMADS58 is a C-class MADS box protein, and characterization of a transposon insertion mutant osmads58 suggested that OsMADS58 plays a role in stamen development. However, as no null mutation has been obtained, its role has remained unclear. Here, we report that the CRISPR knockout mutant osmads58 exhibits complex altered phenotypes, including anomalous diploid germ cells, aberrant meiosis, and delayed tapetum degeneration. This CRISPR mutant line exhibited stronger changes in expression of OsMADS58 target genes compared with the osmads58 dSpm (transposon insertion) line, along with changes in multiple pathways related to early stamen development. Notably, transcriptional regulatory circuits in young panicles covering the stamen at stages 4-6 were substantially altered in the CRISPR line compared to the dSpm line. These findings strongly suggest that the pleiotropic effects of OsMADS58 on stamen development derive from a potential role in stabilizing gene regulatory circuits during early stamen development. Thus, this work opens new avenues for viewing and deciphering the regulatory mechanisms of early stamen development from a network perspective.

IRF4 as an Oncogenic Master Transcription Factor.

Simple SummaryMaster transcription factors regulate essential developmental processes and cellular maintenance that characterize cell identity. Many of them also serve as oncogenes when aberrantly expressed or activated. IRF4 is one of prime examples of oncogenic master transcription factors that has been implicated in various mature lymphoid neoplasms. IRF4 forms unique regulatory circuits and induces oncogenic transcription programs through the interactions with upstream pathways and binding partners. IRF4 is a transcription factor in the interferon regulatory factor (IRF) family. Since the discovery of this gene, various research fields including immunology and oncology have highlighted the unique characteristics and the importance of IRF4 in several biological processes that distinguish it from other IRF family members. In normal lymphocyte development and immunity, IRF4 mediates critical immune responses via interactions with upstream signaling pathways, such as the T-cell receptor and B-cell receptor pathways, as well as their binding partners, which are uniquely expressed in each cell type. On the other hand, IRF4 acts as an oncogene in various mature lymphoid neoplasms when abnormally expressed. IRF4 induces several oncogenes, such as MYC, as well as genes that characterize each cell type by utilizing its ability as a master regulator of immunity. IRF4 and its upstream factor NF-κB form a transcriptional regulatory circuit, including feedback and feedforward loops, to maintain the oncogenic transcriptional program in malignant lymphoid cells. In this review article, we provide an overview of the molecular functions of IRF4 in mature lymphoid neoplasms and highlight its upstream and downstream pathways, as well as the regulatory circuits mediated by IRF4.

Auxotrophic and prototrophic conditional genetic networks reveal the rewiring of transcription factors in Escherichia coli

Bacterial transcription factors (TFs) are widely studied in Escherichia coli. Yet it remains unclear how individual genes in the underlying pathways of TF machinery operate together during environmental challenge. Here, we address this by applying an unbiased, quantitative synthetic genetic interaction (GI) approach to measure pairwise GIs among all TF genes in E. coli under auxotrophic (rich medium) and prototrophic (minimal medium) static growth conditions. The resulting static and differential GI networks reveal condition-dependent GIs, widespread changes among TF genes in metabolism, and new roles for uncharacterized TFs (yjdC, yneJ, ydiP) as regulators of cell division, putrescine utilization pathway, and cold shock adaptation. Pan-bacterial conservation suggests TF genes with GIs are co-conserved in evolution. Together, our results illuminate the global organization of E. coli TFs, and remodeling of genetic backup systems for TFs under environmental change, which is essential for controlling the bacterial transcriptional regulatory circuits.

Engineering eukaryote-like regulatory circuits to expand artificial control mechanisms for metabolic engineering in Saccharomyces cerevisiae

Temporal control of heterologous pathway expression is critical to achieve optimal efficiency in microbial metabolic engineering. The broadly-used GAL promoter system for engineered yeast (Saccharomyces cerevisiae) suffers from several drawbacks; specifically, unintended induction during laboratory development, and unintended repression in industrial production applications, which decreases overall production capacity. Eukaryotic synthetic circuits have not been well examined to address these problems. Here, we explore a modularised engineering method to deploy new genetic circuits applicable for expanding the control of GAL promoter-driven heterologous pathways in S. cerevisiae. Trans- and cis- modules, including eukaryotic trans-activating-and-repressing mechanisms, were characterised to provide new and better tools for circuit design. A eukaryote-like tetracycline-mediated circuit that delivers stringent repression was engineered to minimise metabolic burden during strain development and maintenance. This was combined with a novel 37 °C induction circuit to relief glucose-mediated repression on the GAL promoter during the bioprocess. This delivered a 44% increase in production of the terpenoid nerolidol, to 2.54 g L−1 in flask cultivation. These negative/positive transcriptional regulatory circuits expand global strategies of metabolic control to facilitate laboratory maintenance and for industry applications.

Deciphering Transcriptional Regulatory Circuits: Transcription Factor Binding and Regulatory Variants Identification

Deciphering Transcriptional Regulatory Circuits: Transcription Factor Binding and Regulatory Variants Identification

Engineering of Promoters for Gene Expression in Pichia pastoris.

The availability of exceptionally strong and tightly regulated promoters is a key feature of Komagataella phaffii (syn. Pichia pastoris), a widely applied yeast expression system for heterologous protein production. Most commonly, the methanol-inducible promoter of the alcohol oxidase 1 gene (PAOX1) and the constitutive promoter of the glyceraldehyde 3 phosphate dehydrogenase gene (PGAP) have been used. Recently, also promising novel constitutive (PGCW14), regulated (PGTH1, PCAT1), and bidirectional promoters (histone promoters and synthetic hybrid variants) have been reported.As natural promoters showed so far limited tunability of expression levels and regulatory profiles, various promoter engineering efforts have been undertaken for P. pastoris . PAOX1, PDAS2, PGAP, and PGCW14 have been engineered by systematic deletion studies or random mutagenesis of upstream regulatory sequences. New engineering strategies have focused on PAOX1 core promoter modifications by random or rational approaches and transcriptional regulatory circuits to render PAOX1 independent of methanol induction. These promoter engineering efforts in P. pastoris have resulted in improved, sequence-diversified synthetic promoter variants allowing coordinated fine-tuning of gene expression for a multitude of biotechnological applications.

Cancer CRC: A Comprehensive Cancer Core Transcriptional Regulatory Circuit Resource and Analysis Platform.

A core transcriptional regulatory circuit (CRC) is a group of interconnected auto-regulating transcription factors (TFs) that form loops and can be identified by super-enhancers (SEs). Studies have indicated that CRCs play an important role in defining cellular identity and determining cellular fate. Additionally, core TFs in CRCs are regulators of cell-type-specific transcriptional regulation. However, a global view of CRC properties across various cancer types has not been generated. Thus, we integrated paired cancer ATAC-seq and H3K27ac ChIP-seq data for specific cell lines to develop the Cancer CRC (http://bio.liclab.net/Cancer_crc/index.html). This platform documented 94,108 cancer CRCs, including 325 core TFs. The cancer CRC also provided the “SE active core TFs analysis” and “TF enrichment analysis” tools to identify potentially key TFs in cancer. In addition, we performed a comprehensive analysis of core TFs in various cancer types to reveal conserved and cancer-specific TFs.

A comprehensive integrated post-GWAS analysis of Type 1 diabetes reveals enhancer-based immune dysregulation.

Type 1 diabetes (T1D) is an organ-specific autoimmune disease, whereby immune cell-mediated killing leads to loss of the insulin-producing β cells in the pancreas. Genome-wide association studies (GWAS) have identified over 200 genetic variants associated with risk for T1D. The majority of the GWAS risk variants reside in the non-coding regions of the genome, suggesting that gene regulatory changes substantially contribute to T1D. However, identification of causal regulatory variants associated with T1D risk and their affected genes is challenging due to incomplete knowledge of non-coding regulatory elements and the cellular states and processes in which they function. Here, we performed a comprehensive integrated post-GWAS analysis of T1D to identify functional regulatory variants in enhancers and their cognate target genes. Starting with 1,817 candidate T1D SNPs defined from the GWAS catalog and LDlink databases, we conducted functional annotation analysis using genomic data from various public databases. These include 1) Roadmap Epigenomics, ENCODE, and RegulomeDB for epigenome data; 2) GTEx for tissue-specific gene expression and expression quantitative trait loci data; and 3) lncRNASNP2 for long non-coding RNA data. Our results indicated a prevalent enhancer-based immune dysregulation in T1D pathogenesis. We identified 26 high-probability causal enhancer SNPs associated with T1D, and 64 predicted target genes. The majority of the target genes play major roles in antigen presentation and immune response and are regulated through complex transcriptional regulatory circuits, including those in HLA (6p21) and non-HLA (16p11.2) loci. These candidate causal enhancer SNPs are supported by strong evidence and warrant functional follow-up studies.

Poacic acid, a β‐1,3‐glucan–binding antifungal agent, inhibits cell‐wall remodeling and activates transcriptional responses regulated by the cell‐wall integrity and high‐osmolarity glycerol pathways in yeast

As a result of the relatively few available antifungals and the increasing frequency of resistance to them, the development of novel antifungals is increasingly important. The plant natural product poacic acid (PA) inhibits β-1,3-glucan synthesis in Saccharomyces cerevisiae and has antifungal activity against a wide range of plant pathogens. However, the mode of action of PA is unclear. Here, we reveal that PA specifically binds to β-1,3-glucan, its affinity for which is ~30-fold that for chitin. Besides its effect on β-1,3-glucan synthase activity, PA inhibited the yeast glucan-elongating activity of Gas1 and Gas2 and the chitin-glucan transglycosylase activity of Crh1. Regarding the cellular response to PA, transcriptional co-regulation was mediated by parallel activation of the cell-wall integrity (CWI) and high-osmolarity glycerol signaling pathways. Despite targeting β-1,3-glucan remodeling, the transcriptional profiles and regulatory circuits activated by caspofungin, zymolyase, and PA differed, indicating that their effects on CWI have different mechanisms. The effects of PA on the growth of yeast strains indicated that it has a mode of action distinct from that of echinocandins, suggesting it is a unique antifungal agent.

Precise transcriptional control of cellular quiescence by BRAVO/WOX5 complex in Arabidopsis roots

Understanding stem cell regulatory circuits is the next challenge in plant biology, as these cells are essential for tissue growth and organ regeneration in response to stress. In the Arabidopsis primary root apex, stem cell‐specific transcription factors BRAVO and WOX5 co‐localize in the quiescent centre (QC) cells, where they commonly repress cell division so that these cells can act as a reservoir to replenish surrounding stem cells, yet their molecular connection remains unknown. Genetic and biochemical analysis indicates that BRAVO and WOX5 form a transcription factor complex that modulates gene expression in the QC cells to preserve overall root growth and architecture. Furthermore, by using mathematical modelling we establish that BRAVO uses the WOX5/BRAVO complex to promote WOX5 activity in the stem cells. Our results unveil the importance of transcriptional regulatory circuits in plant stem cell development.

Bone Morphogenic Protein Signaling and Melanoma.

Malignant melanoma is a deadly form of skin cancer caused by neoplastic transformation of melanocytic cells. Despite recent progress in melanoma therapy, by inhibition of activated oncogenes or immunotherapy, survival rate for metastatic melanoma patients remains low. The remarkable phenotypic plasticity of melanoma cells allows for rapid development of invasive properties and metastatic tumors, the main cause of mortality in melanoma patients. Phenotypic and molecular analyses of developing tumors revealed that epithelial-mesenchymal transition (EMT), a cellular and molecular mechanism, controls transition from mature melanocyte to less differentiated melanocyte lineage progenitor cells forming melanoma tumors. This transition is facilitated by persistence of transcriptional regulatory circuit characteristic of embryonic stage in mature melanocytes. Switching of the developmental program of mature melanocyte to EMT is induced by accumulated mutations, especially targeting BRAF, N-RAS, or MEK1/2 signaling pathways, and further promoted by dynamic stimuli from local environment including hypoxia, interactions with extracellular matrix and growth factors or cytokines. Recent reports demonstrate that signaling mediated by transforming growth factor-β (TGF-β) and bone morphogenic proteins (BMPs) play critical roles in inducing EMT by controlling expression of critical transcription factors. BMPs are essential modulators of differentiation, proliferation, apoptosis, invasiveness, and metastases in developing melanoma tumors. They control transcription and epigenetic landscape of melanoma cells. Better understanding of the role of BMPs may lead to new strategies to control EMT processes in melanocyte cell lineage and to achieve clinical benefits for the patients.

Transfer RNA-derived non-coding RNAs (tncRNAs): Hidden regulation of plants' transcriptional regulatory circuits

The emergence of distinct classes of non-coding RNAs has led to better insights into the eukaryotic gene regulatory networks. Amongst them, the existence of transfer RNA (tRNA)-derived non-coding RNAs (tncRNAs) demands exploration in the plant kingdom. We have designed a methodology to uncover the entire perspective of tncRNAome in plants. Using this pipeline, we have identified diverse tncRNAs with a size ranging from 14 to 50 nucleotides (nt) by utilizing 2448 small RNA-seq samples from six angiosperms, and studied their various features, including length, codon-usage, cleavage pattern, and modified tRNA nucleosides. Codon-dependent generation of tncRNAs suggests that the tRNA cleavage is highly specific rather than random tRNA degradation. The nucleotide composition analysis of tncRNA cleavage positions indicates that they are generated through precise endoribonucleolytic cleavage machinery. Certain nucleoside modifications detected on tncRNAs were found to be conserved across the plants, and hence may influence tRNA cleavage, as well as tncRNA functions. Pathway enrichment analysis revealed that common tncRNA targets are majorly enriched during metabolic and developmental processes. Further distinct tissue-specific tncRNA clusters highlight their role in plant development. Significant number of tncRNAs differentially expressed under abiotic and biotic stresses highlights their potential role in stress resistance. In summary, this study has developed a platform that will help in the understanding of tncRNAs and their involvement in growth, development, and response to various stresses. The workflow, software package, and results are freely available at http://nipgr.ac.in/tncRNA.