Embryo size, specification, and homeostasis are regulated by a complex gene regulatory and signaling network. Here we used gene expression signatures of Wnt-activated mouse embryonic stem cell (mESC) clones to reverse engineer an mESC regulatory network. We identify NKX1-2 as a novel master regulator of preimplantation embryo development. We find that Nkx1-2 inhibition reduces nascent RNA synthesis, downregulates genes controlling ribosome biogenesis, RNA translation, and transport, and induces severe alteration of nucleolus structure, resulting in the exclusion of RNA polymerase I from nucleoli. In turn, NKX1-2 loss of function leads to chromosome missegregation in the 2- to 4-cell embryo stages, severe decrease in blastomere numbers, alterations of tight junctions (TJs), and impairment of microlumen coarsening. Overall, these changes impair the blastocoel expansion-collapse cycle and embryo cavitation, leading to altered lineage specification and developmental arrest.

UNASSIGNED: Sporadic and familial amyotrophic lateral sclerosis (ALS) is a fatal progressive neurodegenerative disease that results in loss of motor neurons and, in some patients, associates with frontotemporal dementia (FTD). Apart from the accumulation of proteinaceous deposits, emerging literature indicates that aberrant mitochondrial bioenergetics may contribute to the onset and progression of ALS/FTD. Here we sought to investigate the pathophysiological signatures of mitochondrial dysfunction associated with ALS/FTD.
UNASSIGNED: By means of label-free mass spectrometry (MS) and mRNA sequencing (mRNA-seq), we report pre-symptomatic changes in the cortices of TDP-43 and FUS mutant mouse models. Using tissues from transgenic mouse models of mitochondrial diseases as a reference, we performed comparative analyses and extracted unique and common mitochondrial signatures that revealed neuroprotective compensatory mechanisms in response to early damage.
UNASSIGNED: In this regard, upregulation of both Acyl-CoA Synthetase Long-Chain Family Member 3 (ACSL3) and mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) were the most representative change in pre-symptomatic ALS/FTD tissues, suggesting that fatty acid beta-oxidation and mitochondrial protein translation are mechanisms of adaptation in response to ALS/FTD pathology.
UNASSIGNED: Together, our unbiased integrative analyses unveil novel molecular components that may influence mitochondrial homeostasis in the earliest phase of ALS.

Due to their immunosuppressive role, tumor-infiltrating regulatory T cells (TI-Tregs) represent attractive immuno-oncology targets. Analysis of TI vs. peripheral Tregs (P-Tregs) from 36 patients, across four malignancies, identified 17 candidate master regulators (MRs) as mechanistic determinants of TI-Treg transcriptional state. Pooled CRISPR-Cas9 screening in vivo, using a chimeric hematopoietic stem cell transplant model, confirmed the essentiality of eight MRs in TI-Treg recruitment and/or retention without affecting other T cell subtypes, and targeting one of the most significant MRs (Trps1) by CRISPR KO significantly reduced ectopic tumor growth. Analysis of drugs capable of inverting TI-Treg MR activity identified low-dose gemcitabine as the top prediction. Indeed, gemcitabine treatment inhibited tumor growth in immunocompetent but not immunocompromised allografts, increased anti-PD-1 efficacy, and depleted MR-expressing TI-Tregs in vivo. This study provides key insight into Treg signaling, specifically in the context of cancer, and a generalizable strategy to systematically elucidate and target MR proteins in immunosuppressive subpopulations.

The expression and activity of transcription factors, which directly mediate gene transcription, are strictly regulated to control numerous normal cellular processes. In cancer, transcription factor activity is often dysregulated, resulting in abnormal expression of genes related to tumorigenesis and development. The carcinogenicity of transcription factors can be reduced through targeted therapy. However, most studies on the pathogenic and drug-resistant mechanisms of ovarian cancer have focused on the expression and signaling pathways of individual transcription factors. To improve the prognosis and treatment of patients with ovarian cancer, multiple transcription factors should be evaluated simultaneously to determine the effects of their protein activity on drug therapies. In this study, the transcription factor activity of ovarian cancer samples was inferred from virtual inference of protein activity by enriched regulon algorithm using mRNA expression data. Patients were clustered according to their transcription factor protein activities to investigate the association of transcription factor activities of different subtypes with prognosis and drug sensitivity for filtering subtype-specific drugs. Meanwhile, master regulator analysis was utilized to identify master regulators of differential protein activity between clustering subtypes, thereby identifying transcription factors associated with prognosis and assessing their potential as therapeutic targets. Master regulator risk scores were then constructed for guiding patients\' clinical treatment, providing new insights into the treatment of ovarian cancer at the level of transcriptional regulation.

The molecular mechanisms and gene regulatory networks sustaining cell proliferation in neuroblastoma (NBL) cells are still not fully understood. In this tumor context, it has been proposed that anti-proliferative drugs, such as the pan-HDAC inhibitor panobinostat, could be tested to mitigate tumor progression. Here, we set out to investigate the effects of panobinostat treatment at the unprecedented resolution offered by single-cell sequencing. We identified a global senescence signature paired with reduction in proliferation in treated Kelly cells and more isolated transcriptional responses compatible with early neuronal differentiation. Using master regulator analysis, we identified BAZ1A, HCFC1, MAZ, and ZNF146 as the transcriptional regulators most significantly repressed by panobinostat. Experimental silencing of these transcription factors (TFs) confirmed their role in sustaining NBL cell proliferation in vitro.

Retinoblastoma (Rb) is a rare intraocular tumour in early childhood, with an approximate incidence of 1 in 18 000 live births. Experimental studies for Rb are complex due to the challenges associated with obtaining a normal retina to contrast with diseased tissue. In this work, we reanalyse a dataset that contains normal retina samples. We identified the individual genes whose expression is different in Rb in contrast with normal tissue, determined the pathways whose global expression pattern is more distant from the global expression observed in normal tissue, and finally, we identified which transcription factors regulate the highest number of differentially expressed genes (DEGs) and proposed as transcriptional master regulators (TMRs). The enrichment of DEGs in the phototransduction and retrograde endocannabinoid signalling pathways could be associated with abnormal behaviour of the processes leading to cellular differentiation and cellular proliferation. On the other hand, the TMRs nuclear receptor subfamily 5 group A member 2 and hepatocyte nuclear factor 4 gamma are involved in hepatocyte differentiation. Therefore, the enrichment of aberrant expression in these transcription factors could suggest an abnormal retina development that could be involved in Rb origin and progression.

A cancer immune phenotype characterized by an active T-helper 1 (Th1)/cytotoxic response is associated with responsiveness to immunotherapy and favorable prognosis across different tumors. However, in some cancers, such an intratumoral immune activation does not confer protection from progression or relapse. Defining mechanisms associated with immune evasion is imperative to refine stratification algorithms, to guide treatment decisions and to identify candidates for immune-targeted therapy. Molecular alterations governing mechanisms for immune exclusion are still largely unknown. The availability of large genomic datasets offers an opportunity to ascertain key determinants of differential intratumoral immune response. We follow a network-based protocol to identify transcription regulators (TRs) associated with poor immunologic antitumor activity. We use a consensus of four different pipelines consisting of two state-of-the-art gene regulatory network inference techniques, regularized gradient boosting machines and ARACNE to determine TR regulons, and three separate enrichment techniques, including fast gene set enrichment analysis, gene set variation analysis and virtual inference of protein activity by enriched regulon analysis to identify the most important TRs affecting immunologic antitumor activity. These TRs, referred to as master regulators (MRs), are unique to immune-silent and immune-active tumors, respectively. We validated the MRs coherently associated with the immune-silent phenotype across cancers in The Cancer Genome Atlas and a series of additional datasets in the Prediction of Clinical Outcomes from Genomic Profiles repository. A downstream analysis of MRs specific to the immune-silent phenotype resulted in the identification of several enriched candidate pathways, including NOTCH1, TGF-$\\beta $, Interleukin-1 and TNF-$\\alpha $ signaling pathways. TGFB1I1 emerged as one of the main negative immune modulators preventing the favorable effects of a Th1/cytotoxic response.

Schisandra chinensis fruit extract (SCE) has been used as a traditional medicine for treating vascular diseases. However, little is known about how SCE and schisandrin B (SchB) affect transcriptional output-a crucial factor for shaping the fibrotic responses of the transforming growth factor β (TGFβ) signaling pathways in in vascular smooth muscle cells (VSMC). In this study, to assess the pharmacological effect of SCE and SchB on TGFβ-induced transcriptional output, we performed DNA microarray experiments in A7r5 VSMCs. We found that TGFβ induced distinctive changes in the gene expression profile and that these changes were considerably reversed by SCE and SchB. Gene Set Enrichment Analysis (GSEA) with Hallmark signature suggested that SCE or SchB inhibits a range of fibrosis-associated biological processes, including inflammation, cell proliferation and migration. With our VSMC-specific transcriptional interactome network, master regulator analysis identified crucial transcription factors that regulate the expression of SCE- and SchB-effective genes (i.e., TGFβ-reactive genes whose expression are reversed by SCE and SchB). Our results provide novel perspective and insight into understanding the pharmacological action of SCE and SchB at the transcriptome level and will support further investigations to develop multitargeted strategies for the treatment of vascular fibrosis.

Neuroblastoma (NBL) is a pediatric cancer responsible for more than 15% of cancer deaths in children, with 800 new cases each year in the United States alone. Genomic amplification of the MYC oncogene family member MYCN characterizes a subset of high-risk pediatric neuroblastomas. Several cellular models have been implemented to study this disease over the years. Two of these, SK-N-BE-2-C (BE2C) and Kelly, are amongst the most used worldwide as models of MYCN-Amplified human NBL. Here, we provide a transcriptome-wide quantitative measurement of gene expression and transcriptional network activity in BE2C and Kelly cell lines at an unprecedented single-cell resolution. We obtained 1105 Kelly and 962 BE2C unsynchronized cells, with an average number of mapped reads/cell of roughly 38,000. The single-cell data recapitulate gene expression signatures previously generated from bulk RNA-Seq. We highlight low variance for commonly used housekeeping genes between different cells (ACTB, B2M and GAPDH), while showing higher than expected variance for metallothionein transcripts in Kelly cells. The high number of samples, despite the relatively low read coverage of single cells, allowed for robust pathway enrichment analysis and master regulator analysis (MRA), both of which highlight the more mesenchymal nature of BE2C cells as compared to Kelly cells, and the upregulation of TWIST1 and DNAJC1 transcriptional networks. We further defined master regulators at the single cell level and showed that MYCN is not constantly active or expressed within Kelly and BE2C cells, independently of cell cycle phase. The dataset, alongside a detailed and commented programming protocol to analyze it, is fully shared and reusable.

Colorectal cancer (CRC) has been most extensively studied for characterizing genetic mutations along its development. However, we still have a poor understanding of CRC initiation due to limited measures of its observation and analysis. If we can unveil CRC initiation events, we might identify novel prognostic markers and therapeutic targets for early cancer detection and prevention. To tackle this problem, we establish the early CRC development model and perform transcriptome analysis of its single cell RNA-sequencing data. Interestingly, we find two subtypes, fast growing vs. slowly growing populations of distinct growth rate and gene signatures, and identify CCDC85B as a master regulator that can transform the cellular state of fast growing subtype cells into that of slowly growing subtype cells. We further validate this by in vitro experiments and suggest CCDC85B as a novel potential therapeutic target that may prevent malignant CRC development by suppressing stemness and uncontrolled cell proliferation.