Mutually exclusive teams-like patterns of gene regulation characterize phenotypic heterogeneity along the noradrenergic-mesenchymal axis in neuroblastoma

Differential Regulation of Opsin Gene Expression in Response to Internal and External Stimuli

Determining how internal and external stimuli interact to determine developmental trajectories of traits is a challenge that requires the integration of different subfields of biology. Internal stimuli, such as hormones, control developmental patterns of phenotypic changes, which might be modified by external environmental cues (e.g. plasticity). Thyroid hormone (TH) modulates the timing of opsin gene expression in developing Midas cichlid fish (Amphilophus citrinellus). Moreover, fish reared in red light accelerate this developmental timing compared to fish reared in white light. Hence, we hypothesized that plasticity caused by variation in light conditions has coopted the TH signaling pathway to induce changes in opsin gene expression. We treated Midas cichlids with TH and crossed this treatment with two light conditions, white and red. We observed that not only opsin expression responded similarly to TH and red light but also that, at high TH levels, there is limited capacity for light-induced plasticity. Transcriptomic analysis of the eye showed that genes in the TH pathway were affected by TH, but not by light treatments. Coexpression network analyses further suggested that response to light was independent of the response to TH manipulations. Taken together, our results suggest independent mechanisms mediating development and plasticity during development of opsin gene expression, and that responses to environmental stimuli may vary depending on internal stimuli. This conditional developmental response to external factors depending on internal ones (e.g. hormones) might play a fundamental role in the patterns of phenotypic divergence observed in Midas cichlids and potentially other organisms.

Proneural-mesenchymal antagonism dominates the patterns of phenotypic heterogeneity in glioblastoma

The aggressive nature of glioblastoma (GBM) - one of the deadliest forms of brain tumors - is majorly attributed to underlying phenotypic heterogeneity. Early attempts to classify this heterogeneity at a transcriptomic level in TCGA GBM cohort proposed the existence of four distinct molecular subtypes: Proneural, Neural, Classical, and Mesenchymal. Further, a single-cell RNA sequencing (scRNA-seq) analysis of primary tumors also reported similar four subtypes mimicking neurodevelopmental lineages. However, it remains unclear whether these four subtypes identified via bulk and single-cell transcriptomics are mutually exclusive or not. Here, we perform pairwise correlations among individual genes and gene signatures corresponding to these proposed subtypes and show that the subtypes are not distinctly mutually antagonistic in either TCGA or scRNA-seq data. We observed that the proneural (or neural progenitor-like)-mesenchymal axis is the most prominent antagonistic pair, with the other two subtypes lying on this spectrum. These results are reinforced through a meta-analysis of over 100 single-cell and bulk transcriptomic datasets as well as in terms of functional association with metabolic switching, cell cycle, and immune evasion pathways. Finally, this proneural-mesenchymal antagonistic trend percolates to the association of relevant transcription factors with patient survival. These results suggest rethinking GBM phenotypic characterization for more effective therapeutic targeting efforts.