Challenges of Decoding Transcription Factor Dynamics in Terms of Gene Regulation

Co-imaging of RelA and c-Rel reveals features of NF-κB signaling for ligand discrimination

Individual cell sensing of external cues has evolved through the temporal patterns in signaling. Since nuclear factor κB (NF-κB) signaling dynamics have been examined using a single subunit, RelA, it remains unclear whether more information might be transmitted via other subunits. Using NF-κB double-knockin reporter mice, we monitored both canonical NF-κB subunits, RelA and c-Rel, simultaneously in single macrophages by quantitative live-cell imaging. We show that signaling features of RelA and c-Rel convey more information about the stimuli than those of either subunit alone. Machine learning is used to predict the ligand identity accurately based on RelA and c-Rel signaling features without considering the co-activated factors. Ligand discrimination is achieved through selective non-redundancy of RelA and c-Rel signaling dynamics, as well as their temporal coordination. These results suggest a potential role of c-Rel in fine-tuning immune responses and highlight the need for approaches that will elucidate the mechanisms regulating NF-κB subunit specificity.

Gene Regulatory Networks in Coronary Artery Disease

PURPOSE OF REVIEW

Coronary artery disease is a complex disorder and the leading cause of mortality worldwide. As technologies for the generation of high-throughput multiomics data have advanced, gene regulatory network modeling has become an increasingly powerful tool in understanding coronary artery disease. This review summarizes recent and novel gene regulatory network tools for bulk tissue and single cell data, existing databases for network construction, and applications of gene regulatory networks in coronary artery disease.

RECENT FINDINGS

New gene regulatory network tools can integrate multiomics data to elucidate complex disease mechanisms at unprecedented cellular and spatial resolutions. At the same time, updates to coronary artery disease expression data in existing databases have enabled researchers to build gene regulatory networks to study novel disease mechanisms. Gene regulatory networks have proven extremely useful in understanding CAD heritability beyond what is explained by GWAS loci and in identifying mechanisms and key driver genes underlying disease onset and progression. Gene regulatory networks can holistically and comprehensively address the complex nature of coronary artery disease. In this review, we discuss key algorithmic approaches to construct gene regulatory networks and highlight state-of-the-art methods that model specific modes of gene regulation. We also explore recent applications of these tools in coronary artery disease patient data repositories to understand disease heritability and shared and distinct disease mechanisms and key driver genes across tissues, between sexes, and between species.

Tuning Plant Promoters Using a Simple Split Luciferase Method to Assess Transcription Factor-DNA Interactions

Sequence features, including the affinity of binding motifs for their cognate transcription factors, are important contributors to promoter behavior. The ability to predictably recode affinity enables the development of synthetic promoters with varying levels of response to known cellular signals. Here we describe a luminescence-based microplate assay for comparing the interactions of transcription factors with short DNA probes. We then demonstrate how these data can be used to design synthetic plant promoters of varying strengths that respond to the same transcription factor.

SCING: Inference of robust, interpretable gene regulatory networks from single cell and spatial transcriptomics

Gene regulatory network (GRN) inference is an integral part of understanding physiology and disease. Single cell/nuclei RNA-seq (scRNA-seq/snRNA-seq) data has been used to elucidate cell-type GRNs; however, the accuracy and speed of current scRNAseq-based GRN approaches are suboptimal. Here, we present Single Cell INtegrative Gene regulatory network inference (SCING), a gradient boosting and mutual information-based approach for identifying robust GRNs from scRNA-seq, snRNA-seq, and spatial transcriptomics data. Performance evaluation using Perturb-seq datasets, held-out data, and the mouse cell atlas combined with the DisGeNET database demonstrates the improved accuracy and biological interpretability of SCING compared to existing methods. We applied SCING to the entire mouse single cell atlas, human Alzheimer's disease (AD), and mouse AD spatial transcriptomics. SCING GRNs reveal unique disease subnetwork modeling capabilities, have intrinsic capacity to correct for batch effects, retrieve disease relevant genes and pathways, and are informative on spatial specificity of disease pathogenesis.

De-novo reconstruction and identification of transcriptional gene regulatory network modules differentiating single-cell clusters

Single-cell RNA sequencing (scRNA-seq) technology provides an unprecedented opportunity to understand gene functions and interactions at single-cell resolution. While computational tools for scRNA-seq data analysis to decipher differential gene expression profiles and differential pathway expression exist, we still lack methods to learn differential regulatory disease mechanisms directly from the single-cell data. Here, we provide a new methodology, named DiNiro, to unravel such mechanisms de novo and report them as small, easily interpretable transcriptional regulatory network modules. We demonstrate that DiNiro is able to uncover novel, relevant, and deep mechanistic models that not just predict but explain differential cellular gene expression programs. DiNiro is available at https://exbio.wzw.tum.de/diniro/.

NF-κB dynamics in the language of immune cells

Biological discovery has been driven by advances in throughput and resolution of analysis technologies. They have also created an indelible bias for snapshot-based knowledge. Even though recent methods such as multi-omics single-cell assays have empowered immunological investigations, they still provide snapshots of cellular behaviors and thus, have inherent limitations in reconstructing unsynchronized dynamic events across individual cells. Here, we present a rationale for how NF-κB may convey specificity of contextual information through subtle quantitative features of its signaling dynamics. The next frontier of predictive understanding should involve functional characterization of NF-κB signaling dynamics and their immunological implications. This may help solve the apparent paradox that a ubiquitously activated transcription factor can shape accurate responses to different immune challenges.

Transcriptome Analysis of Goat Mammary Gland Tissue Reveals the Adaptive Strategies and Molecular Mechanisms of Lactation and Involution

To understand how genes precisely regulate lactation physiological activity and the molecular genetic mechanisms underlying mammary gland involution, this study investigated the transcriptome characteristics of goat mammary gland tissues at the late gestation (LG), early lactation (EL), peak lactation (PL), late lactation (LL), dry period (DP), and involution (IN) stages. A total of 13,083 differentially expressed transcripts were identified by mutual comparison of mammary gland tissues at six developmental stages. Genes related to cell growth, apoptosis, immunity, nutrient transport, synthesis, and metabolism make adaptive transcriptional changes to meet the needs of mammary lactation. Notably, platelet derived growth factor receptor beta (PDGFRB) was screened as a hub gene of the mammary gland developmental network, which is highly expressed during the DP and IN. Overexpression of PDGFRB in vitro could slow down the G1/S phase arrest of goat mammary epithelial cell cycle and promote cell proliferation by regulating the PI3K/Akt signaling pathway. In addition, PDGFRB overexpression can also affect the expression of genes related to apoptosis, matrix metalloproteinase family, and vascular development, which is beneficial to the remodeling of mammary gland tissue during involution. These findings provide new insights into the molecular mechanisms involved in lactation and mammary gland involution.

Bilateral Feedback in Oscillator Model Is Required to Explain the Coupling Dynamics of Hes1 with the Cell Cycle

Biological processes are governed by the expression of proteins, and for some proteins, their level of expression can fluctuate periodically over time (i.e., they oscillate). Many oscillatory proteins (e.g., cell cycle proteins and those from the HES family of transcription factors) are connected in complex ways, often within large networks. This complexity can be elucidated by developing intuitive mathematical models that describe the underlying critical aspects of the relationships between these processes. Here, we provide a mathematical explanation of a recently discovered biological phenomenon: the phasic position of the gene Hes1’s oscillatory expression at the beginning of the cell cycle of an individual human breast cancer stem cell can have a predictive value on how long that cell will take to complete a cell cycle. We use a two-component model of coupled oscillators to represent Hes1 and the cell cycle in the same cell with minimal assumptions. Inputting only the initial phase angles, we show that this model is capable of predicting the dynamic mitosis to mitosis behaviour of Hes1 and predicting cell cycle length patterns as found in real-world experimental data. Moreover, we discover that bidirectional coupling between Hes1 and the cell cycle is critical within the system for the data to be reproduced and that nonfixed asymmetry in the interactions between the oscillators is required. The phase dynamics we present here capture the complex interplay between Hes1 and the cell cycle, helping to explain nongenetic cell cycle variability, which has critical implications in cancer treatment contexts.

The bone marrow microenvironment of pre-B acute lymphoblastic leukemia at single-cell resolution

The bone marrow microenvironment (BMM) plays a key role in leukemia progression, but its molecular complexity in pre-B cell acute lymphoblastic leukemia (B-ALL), the most common cancer in children, remains poorly understood. To gain further insight, we used single-cell RNA sequencing to characterize the kinetics of the murine BMM during B-ALL progression. Normal pro- and pre-B cells were found to be the most affected at the earliest stages of disease and this was associated with changes in expression of genes regulated by the AP1-transcription factor complex and regulatory factors NELFE, MYC and BCL11A. Granulocyte–macrophage progenitors show reduced expression of the tumor suppressor long non-coding RNA Neat1 and disruptions in the rate of transcription. Intercellular communication networks revealed monocyte-dendritic precursors to be consistently active during B-ALL progression, with enriched processes including cytokine-mediated signaling pathway, neutrophil-mediated immunity and regulation of cell migration and proliferation. In addition, we confirmed that the hematopoietic stem and progenitor cell compartment was perturbed during leukemogenesis. These findings extend our understanding of the complexity of changes and molecular interactions among the normal cells of the BMM during B-ALL progression.

Quantitative relationships between SMAD dynamics and target gene activation kinetics in single live cells

The transduction of extracellular signals through signaling pathways that culminate in a transcriptional response is central to many biological processes. However, quantitative relationships between activities of signaling pathway components and transcriptional output of target genes remain poorly explored. Here we developed a dual bioluminescence imaging strategy allowing simultaneous monitoring of nuclear translocation of the SMAD4 and SMAD2 transcriptional activators upon TGF-β stimulation, and the transcriptional response of the endogenous connective tissue growth factor (ctgf) gene. Using cell lines allowing to vary exogenous SMAD4/2 expression levels, we performed quantitative measurements of the temporal profiles of SMAD4/2 translocation and ctgf transcription kinetics in hundreds of individual cells at high temporal resolution. We found that while nuclear translocation efficiency had little impact on initial ctgf transcriptional activation, high total cellular SMAD4 but not SMAD2 levels increased the probability of cells to exhibit a sustained ctgf transcriptional response. The approach we present here allows time-resolved single cell quantification of transcription factor dynamics and transcriptional responses and thereby sheds light on the quantitative relationship between SMADs and target gene responses.