Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma.

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

Combined protein and nucleic acid imaging reveals virus-dependent B cell and macrophage immunosuppression of tissue microenvironments.

Understanding the mechanisms of HIV tissue persistence necessitates the ability to visualize tissue microenvironments where infected cells reside; however, technological barriers limit our ability to dissect the cellular components of these HIV reservoirs. Here, we developed protein and nucleic acid in situ imaging (PANINI) to simultaneously quantify DNA, RNA, and protein levels within these tissue compartments. By coupling PANINI with multiplexed ion beam imaging (MIBI), we measured over 30 parameters simultaneously across archival lymphoid tissues from healthy or simian immunodeficiency virus (SIV)-infected nonhuman primates. PANINI enabled the spatial dissection of cellular phenotypes, functional markers, and viral events resulting from infection. SIV infection induced IL-10 expression in lymphoid B cells, which correlated with local macrophage M2 polarization. This highlights a potential viral mechanism for conditioning an immunosuppressive tissue environment for virion production. The spatial multimodal framework here can be extended to decipher tissue responses in other infectious diseases and tumor biology.

Asprosin, a Fasting-Induced Glucogenic Protein Hormone.

Hepatic glucose release into the circulation is vital for brain function and survival during periods of fasting and is modulated by an array of hormones that precisely regulate plasma glucose levels. We have identified a fasting-induced protein hormone that modulates hepatic glucose release. It is the C-terminal cleavage product of profibrillin, and we name it Asprosin. Asprosin is secreted by white adipose, circulates at nanomolar levels, and is recruited to the liver, where it activates the G protein-cAMP-PKA pathway, resulting in rapid glucose release into the circulation. Humans and mice with insulin resistance show pathologically elevated plasma asprosin, and its loss of function via immunologic or genetic means has a profound glucose- and insulin-lowering effect secondary to reduced hepatic glucose release. Asprosin represents a glucogenic protein hormone, and therapeutically targeting it may be beneficial in type II diabetes and metabolic syndrome.

Organization of the human intestine at single-cell resolution.

The intestine is a complex organ that promotes digestion, extracts nutrients, participates in immune surveillance, maintains critical symbiotic relationships with microbiota and affects overall health1. The intesting has a length of over nine metres, along which there are differences in structure and function2. The localization of individual cell types, cell type development trajectories and detailed cell transcriptional programs probably drive these differences in function. Here, to better understand these differences, we evaluated the organization of single cells using multiplexed imaging and single-nucleus RNA and open chromatin assays across eight different intestinal sites from nine donors. Through systematic analyses, we find cell compositions that differ substantially across regions of the intestine and demonstrate the complexity of epithelial subtypes, and find that the same cell types are organized into distinct neighbourhoods and communities, highlighting distinct immunological niches that are present in the intestine. We also map gene regulatory differences in these cells that are suggestive of a regulatory differentiation cascade, and associate intestinal disease heritability with specific cell types. These results describe the complexity of the cell composition, regulation and organization for this organ, and serve as an important reference map for understanding human biology and disease.

Robust single-cell matching and multimodal analysis using shared and distinct features.

The ability to align individual cellular information from multiple experimental sources is fundamental for a systems-level understanding of biological processes. However, currently available tools are mainly designed for single-cell transcriptomics matching and integration, and generally rely on a large number of shared features across datasets for cell matching. This approach underperforms when applied to single-cell proteomic datasets due to the limited number of parameters simultaneously accessed and lack of shared markers across these experiments. Here, we introduce a cell-matching algorithm, matching with partial overlap (MARIO) that accounts for both shared and distinct features, while consisting of vital filtering steps to avoid suboptimal matching. MARIO accurately matches and integrates data from different single-cell proteomic and multimodal methods, including spatial techniques and has cross-species capabilities. MARIO robustly matched tissue macrophages identified from COVID-19 lung autopsies via codetection by indexing imaging to macrophages recovered from COVID-19 bronchoalveolar lavage fluid by cellular indexing of transcriptomes and epitopes by sequencing, revealing unique immune responses within the lung microenvironment of patients with COVID.

Twelve-hour rhythms in transcript expression within the human dorsolateral prefrontal cortex are altered in schizophrenia.

Twelve-hour (12 h) ultradian rhythms are a well-known phenomenon in coastal marine organisms. While 12 h cycles are observed in human behavior and physiology, no study has measured 12 h rhythms in the human brain. Here, we identify 12 h rhythms in transcripts that either peak at sleep/wake transitions (approximately 9 AM/PM) or static times (approximately 3 PM/AM) in the dorsolateral prefrontal cortex, a region involved in cognition. Subjects with schizophrenia (SZ) lose 12 h rhythms in genes associated with the unfolded protein response and neuronal structural maintenance. Moreover, genes involved in mitochondrial function and protein translation, which normally peak at sleep/wake transitions, peak instead at static times in SZ, suggesting suboptimal timing of these essential processes.

TRAF4-mediated nonproteolytic ubiquitination of androgen receptor promotes castration-resistant prostate cancer.

Castration-resistant prostate cancer (CRPC) poses a major clinical challenge with the androgen receptor (AR) remaining to be a critical oncogenic player. Several lines of evidence indicate that AR induces a distinct transcriptional program after androgen deprivation in CRPCs. However, the mechanism triggering AR binding to a distinct set of genomic loci in CRPC and how it promotes CRPC development remain unclear. We demonstrate here that atypical ubiquitination of AR mediated by an E3 ubiquitin ligase TRAF4 plays an important role in this process. TRAF4 is highly expressed in CRPCs and promotes CRPC development. It mediates K27-linked ubiquitination at the C-terminal tail of AR and increases its association with the pioneer factor FOXA1. Consequently, AR binds to a distinct set of genomic loci enriched with FOXA1- and HOXB13-binding motifs to drive different transcriptional programs including an olfactory transduction pathway. Through the surprising upregulation of olfactory receptor gene transcription, TRAF4 increases intracellular cAMP levels and boosts E2F transcription factor activity to promote cell proliferation under androgen deprivation conditions. Altogether, these findings reveal a posttranslational mechanism driving AR-regulated transcriptional reprogramming to provide survival advantages for prostate cancer cells under castration conditions.

Beyond circadian rhythms: emerging roles of ultradian rhythms in control of liver functions.

The mammalian liver must cope with various metabolic and physiological changes that normally recur every day and primarily stem from daily cycles of rest-activity and fasting-feeding. Although a large body of evidence supports the reciprocal regulation of circadian rhythms and liver function, the research on the hepatic ultradian rhythms have largely been lagging behind. However, with the advent of more cost-effective high-throughput omics technologies, high-resolution time-lapse imaging, and more robust and powerful mathematical tools, several recent studies have shed new light on the presence and functions of hepatic ultradian rhythms. In this review, we will first very briefly discuss the basic principles of circadian rhythms, and then cover in greater details the recent literature related to ultradian rhythms. Specifically, we will highlight the prevalence and mechanisms of hepatic 12-h rhythms, and 8-h rhythms, which cycle at the second and third harmonics of circadian frequency. Finally, we also refer to ultradian rhythms with other frequencies and examine the limitations of the current approaches as well as the challenges related to identifying ultradian rhythm and addressing their molecular underpinnings.

Metabolic enzyme PFKFB4 activates transcriptional coactivator SRC-3 to drive breast cancer.

Alterations in both cell metabolism and transcriptional programs are hallmarks of cancer that sustain rapid proliferation and metastasis 1 . However, the mechanisms that control the interaction between metabolic reprogramming and transcriptional regulation remain unclear. Here we show that the metabolic enzyme 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 4 (PFKFB4) regulates transcriptional reprogramming by activating the oncogenic steroid receptor coactivator-3 (SRC-3). We used a kinome-wide RNA interference-based screening method to identify potential kinases that modulate the intrinsic SRC-3 transcriptional response. PFKFB4, a regulatory enzyme that synthesizes a potent stimulator of glycolysis 2 , is found to be a robust stimulator of SRC-3 that coregulates oestrogen receptor. PFKFB4 phosphorylates SRC-3 at serine 857 and enhances its transcriptional activity, whereas either suppression of PFKFB4 or ectopic expression of a phosphorylation-deficient Ser857Ala mutant SRC-3 abolishes the SRC-3-mediated transcriptional output. Functionally, PFKFB4-driven SRC-3 activation drives glucose flux towards the pentose phosphate pathway and enables purine synthesis by transcriptionally upregulating the expression of the enzyme transketolase. In addition, the two enzymes adenosine monophosphate deaminase-1 (AMPD1) and xanthine dehydrogenase (XDH), which are involved in purine metabolism, were identified as SRC-3 targets that may or may not be directly involved in purine synthesis. Mechanistically, phosphorylation of SRC-3 at Ser857 increases its interaction with the transcription factor ATF4 by stabilizing the recruitment of SRC-3 and ATF4 to target gene promoters. Ablation of SRC-3 or PFKFB4 suppresses breast tumour growth in mice and prevents metastasis to the lung from an orthotopic setting, as does Ser857Ala-mutant SRC-3. PFKFB4 and phosphorylated SRC-3 levels are increased and correlate in oestrogen receptor-positive tumours, whereas, in patients with the basal subtype, PFKFB4 and SRC-3 drive a common protein signature that correlates with the poor survival of patients with breast cancer. These findings suggest that the Warburg pathway enzyme PFKFB4 acts as a molecular fulcrum that couples sugar metabolism to transcriptional activation by stimulating SRC-3 to promote aggressive metastatic tumours.

Logic of the Temporal Compartmentalization of the Hepatic Metabolic Cycle.

The mammalian liver must cope with various metabolic and physiological changes that normally recur every day and result primarily from rest-activity and fasting-feeding cycles. In this article, I present evidence supporting a temporal compartmentalization of rhythmic hepatic metabolic processes into four main clusters: regulation of energy homeostasis, maintenance of information integrity, immune response, and genetic information flow. I further review literatures and discuss how both the circadian and the newly discovered 12-h ultradian clock work together to regulate these four temporally separated processes in mouse liver, which, interestingly, is largely uncoupled from the liver zonation regulation.

12-h clock regulation of genetic information flow by XBP1s.

Our group recently characterized a cell-autonomous mammalian 12-h clock independent from the circadian clock, but its function and mechanism of regulation remain poorly understood. Here, we show that in mouse liver, transcriptional regulation significantly contributes to the establishment of 12-h rhythms of mRNA expression in a manner dependent on Spliced Form of X-box Binding Protein 1 (XBP1s). Mechanistically, the motif stringency of XBP1s promoter binding sites dictates XBP1s's ability to drive 12-h rhythms of nascent mRNA transcription at dawn and dusk, which are enriched for basal transcription regulation, mRNA processing and export, ribosome biogenesis, translation initiation, and protein processing/sorting in the Endoplasmic Reticulum (ER)-Golgi in a temporal order consistent with the progressive molecular processing sequence described by the central dogma information flow (CEDIF). We further identified GA-binding proteins (GABPs) as putative novel transcriptional regulators driving 12-h rhythms of gene expression with more diverse phases. These 12-h rhythms of gene expression are cell autonomous and evolutionarily conserved in marine animals possessing a circatidal clock. Our results demonstrate an evolutionarily conserved, intricate network of transcriptional control of the mammalian 12-h clock that mediates diverse biological pathways. We speculate that the 12-h clock is coopted to accommodate elevated gene expression and processing in mammals at the two rush hours, with the particular genes processed at each rush hour regulated by the circadian and/or tissue-specific pathways.

Coactivator-Dependent Oscillation of Chromatin Accessibility Dictates Circadian Gene Amplitude via REV-ERB Loading.

A central mechanism for controlling circadian gene amplitude remains elusive. We present evidence for a "facilitated repression (FR)" model that functions as an amplitude rheostat for circadian gene oscillation. We demonstrate that ROR and/or BMAL1 promote global chromatin decondensation during the activation phase of the circadian cycle to actively facilitate REV-ERB loading for repression of circadian gene expression. Mechanistically, we found that SRC-2 dictates global circadian chromatin remodeling through spatial and temporal recruitment of PBAF members of the SWI/SNF complex to facilitate loading of REV-ERB in the hepatic genome. Mathematical modeling highlights how the FR model sustains proper circadian rhythm despite fluctuations of REV-ERB levels. Our study not only reveals a mechanism for active communication between the positive and negative limbs of the circadian transcriptional loop but also establishes the concept that clock transcription factor binding dynamics is perhaps a central tenet for fine-tuning circadian rhythm.

Revealing the hidden reality of the mammalian 12-h ultradian rhythms.

Biological oscillations often cycle at different harmonics of the 24-h circadian rhythms, a phenomenon we coined "Musica Universalis" in 2017. Like the circadian rhythm, the 12-h oscillation is also evolutionarily conserved, robust, and has recently gained new traction in the field of chronobiology. Originally thought to be regulated by the circadian clock and/or environmental cues, recent new evidences support the notion that the majority of 12-h rhythms are regulated by a distinct and cell-autonomous pacemaker that includes the unfolded protein response (UPR) transcription factor spliced form of XBP1 (XBP1s). 12-h cycle of XBP1s level in turn transcriptionally generates robust 12-h rhythms of gene expression enriched in the central dogma information flow (CEDIF) pathway. Given the regulatory and functional separation of the 12-h and circadian clocks, in this review, we will focus our attention on the mammalian 12-h pacemaker, and discuss our current understanding of its prevalence, evolutionary origin, regulation, and functional roles in both physiological and pathological processes.

Correction: SRC-2-mediated coactivation of anti-tumorigenic target genes suppresses MYC-induced liver cancer.

[This corrects the article DOI: 10.1371/journal.pgen.1006650.].

SRC-2-mediated coactivation of anti-tumorigenic target genes suppresses MYC-induced liver cancer.

Hepatocellular carcinoma (HCC) is the fifth most common solid tumor in the world and the third leading cause of cancer-associated deaths. A Sleeping Beauty-mediated transposon mutagenesis screen previously identified mutations that cooperate with MYC to accelerate liver tumorigenesis. This revealed a tumor suppressor role for Steroid Receptor Coactivator 2/Nuclear Receptor Coactivator 2 (Src-2/Ncoa2) in liver cancer. In contrast, SRC-2 promotes survival and metastasis in prostate cancer cells, suggesting a tissue-specific and context-dependent role for SRC-2 in tumorigenesis. To determine if genetic loss of SRC-2 is sufficient to accelerate MYC-mediated liver tumorigenesis, we bred Src-2-/- mice with a MYC-induced liver tumor model and observed a significant increase in liver tumor burden. RNA sequencing of liver tumors and in vivo chromatin immunoprecipitation assays revealed a set of direct target genes that are bound by SRC-2 and exhibit downregulated expression in Src-2-/- liver tumors. We demonstrate that activation of SHP (Small Heterodimer Partner), DKK4 (Dickkopf-4), and CADM4 (Cell Adhesion Molecule 4) by SRC-2 suppresses tumorigenesis in vitro and in vivo. These studies suggest that SRC-2 may exhibit oncogenic or tumor suppressor activity depending on the target genes and nuclear receptors that are expressed in distinct tissues and illuminate the mechanisms of tumor suppression by SRC-2 in liver.

Regulatory mechanisms mediated by peroxisome proliferator-activated receptor-ß/δ in skin cancer.

Considerable progress has been made during the past 20 years towards elucidating the role of peroxisome proliferator-activated receptor-ß/δ (PPARß/δ) in skin cancer. In 1999, the original notion that PPARß/δ was involved with epithelial cell function was postulated based on a correlation between PPARß/δ expression and the induction of messenger RNAs encoding proteins that mediate terminal differentiation in keratinocytes. Subsequent studies definitively revealed that PPARß/δ could induce terminal differentiation and inhibit proliferation of keratinocytes. Molecular mechanisms have since been discovered to explain how this nuclear receptor can be targeted for preventing and treating skin cancer. This includes the regulation of terminal differentiation, mitotic signaling, endoplasmic reticulum stress, and cellular senescence. Interestingly, the effects of activating PPARß/δ can preferentially target keratinocytes with genetic mutations associated with skin cancer. This review provides the history and current understanding of how PPARß/δ can be targeted for both nonmelanoma skin cancer and melanoma and postulates how future approaches that modulate PPARß/δ signaling may be developed for the prevention and treatment of these diseases.

Acetylation on histone H3 lysine 9 mediates a switch from transcription initiation to elongation.

The transition from transcription initiation to elongation is a key regulatory step in gene expression, which requires RNA polymerase II (pol II) to escape promoter proximal pausing on chromatin. Although elongation factors promote pause release leading to transcription elongation, the role of epigenetic modifications during this critical transition step is poorly understood. Two histone marks on histone H3, lysine 4 trimethylation (H3K4me3) and lysine 9 acetylation (H3K9ac), co-localize on active gene promoters and are associated with active transcription. H3K4me3 can promote transcription initiation, yet the functional role of H3K9ac is much less understood. We hypothesized that H3K9ac may function downstream of transcription initiation by recruiting proteins important for the next step of transcription. Here, we describe a functional role for H3K9ac in promoting pol II pause release by directly recruiting the super elongation complex (SEC) to chromatin. H3K9ac serves as a substrate for direct binding of the SEC, as does acetylation of histone H4 lysine 5 to a lesser extent. Furthermore, lysine 9 on histone H3 is necessary for maximal pol II pause release through SEC action, and loss of H3K9ac increases the pol II pausing index on a subset of genes in HeLa cells. At select gene promoters, H3K9ac loss or SEC depletion reduces gene expression and increases paused pol II occupancy. We therefore propose that an ordered histone code can promote progression through the transcription cycle, providing new mechanistic insight indicating that SEC recruitment to certain acetylated histones on a subset of genes stimulates the subsequent release of paused pol II needed for transcription elongation.

SRC-2 orchestrates polygenic inputs for fine-tuning glucose homeostasis.

Despite extensive efforts to understand the monogenic contributions to perturbed glucose homeostasis, the complexity of genetic events that fractionally contribute to the spectrum of this pathology remain poorly understood. Proper maintenance of glucose homeostasis is the central feature of a constellation of comorbidities that define the metabolic syndrome. The ability of the liver to balance carbohydrate uptake and release during the feeding-to-fasting transition is essential to the regulation of peripheral glucose availability. The liver coordinates the expression of gene programs that control glucose absorption, storage, and secretion. Herein, we demonstrate that Steroid Receptor Coactivator 2 (SRC-2) orchestrates a hierarchy of nutritionally responsive transcriptional complexes to precisely modulate plasma glucose availability. Using DNA pull-down technology coupled with mass spectrometry, we have identified SRC-2 as an indispensable integrator of transcriptional complexes that control the rate-limiting steps of hepatic glucose release and accretion. Collectively, these findings position SRC-2 as a major regulator of polygenic inputs to metabolic gene regulation and perhaps identify a previously unappreciated model that helps to explain the clinical spectrum of glucose dysregulation.

Peroxisome proliferator-activated receptor-ß/δ inhibits human neuroblastoma cell tumorigenesis by inducing p53- and SOX2-mediated cell differentiation.

Neuroblastoma is a common childhood cancer typically treated by inducing differentiation with retinoic acid (RA). Peroxisome proliferator-activated receptor-ß/δ, (PPARß/δ) is known to promote terminal differentiation of many cell types. In the present study, PPARß/δ was over-expressed in three human neuroblastoma cell lines, NGP, SK-N-BE(2), and IMR-32, that exhibit high, medium, and low sensitivity, respectively, to retinoic acid-induced differentiation to determine if PPARß/δ and retinoic acid receptors (RARs) could be jointly targeted to increase the efficacy of treatment. All-trans-RA (atRA) decreased expression of SRY (sex determining region Y)-box 2 (SOX2), a stem cell regulator and marker of de-differentiation, in NGP and SK-N-BE(2) cells with inactive or mutant tumor suppressor p53, respectively. However, atRA did not suppress SOX2 expression in IMR-32 cells carrying wild-type p53. Over-expression and/or ligand activation of PPARß/δ reduced the average volume and weight of ectopic tumor xenografts from NGP, SK-N-BE(2), or IMR-32 cells compared to controls. Compared with that found with atRA, PPARß/δ suppressed SOX2 expression in NGP and SK-N-BE(2) cells and ectopic xenografts, and was also effective in suppressing SOX2 expression in IMR-32 cells that exhibit higher p53 expression compared to the former cell lines. Combined, these observations demonstrate that activating or over-expressing PPARß/δ induces cell differentiation through p53- and SOX2-dependent signaling pathways in neuroblastoma cells and tumors. This suggests that combinatorial activation of both RARα and PPARß/δ may be suitable as an alternative therapeutic approach for RA-resistant neuroblastoma patients.