RORγt is well recognized as the lineage-defining transcription factor for T helper 17 (TH17) cell development. However, the cell-intrinsic mechanisms that negatively regulate TH17 cell development and autoimmunity remain poorly understood. Here, we demonstrate that the transcriptional repressor REV-ERBα is exclusively expressed in TH17 cells, competes with RORγt for their shared DNA consensus sequence, and negatively regulates TH17 cell development via repression of genes traditionally characterized as RORγt dependent, including Il17a. Deletion of REV-ERBα enhanced TH17-mediated pro-inflammatory cytokine expression, exacerbating experimental autoimmune encephalomyelitis (EAE) and colitis. Treatment with REV-ERB-specific synthetic ligands, which have similar phenotypic properties as RORγ modulators, suppressed TH17 cell development, was effective in colitis intervention studies, and significantly decreased the onset, severity, and relapse rate in several models of EAE without affecting thymic cellularity. Our results establish that REV-ERBα negatively regulates pro-inflammatory TH17 responses in vivo and identifies the REV-ERBs as potential targets for the treatment of TH17-mediated autoimmune diseases.

A variety of mechanosensory neurons are involved in touch, proprioception, and pain. Many molecular components of the mechanotransduction machinery subserving these sensory modalities remain to be discovered. Here, we combine recordings of mechanosensitive (MS) currents in mechanosensory neurons with single-cell RNA sequencing. Transcriptional profiles are mapped onto previously identified sensory neuron types to identify cell-type correlates between datasets. Correlation of current signatures with single-cell transcriptomes provides a one-to-one correspondence between mechanoelectric properties and transcriptomically defined neuronal populations. Moreover, a gene-expression differential comparison provides a set of candidate genes for mechanotransduction complexes. Piezo2 is expectedly found to be enriched in rapidly adapting MS current-expressing neurons, whereas Tmem120a and Tmem150c, thought to mediate slow-type MS currents, are uniformly expressed in all mechanosensory neuron subtypes. Further knockdown experiments disqualify them as mediating MS currents in sensory neurons. This dataset constitutes an open resource to explore further the cell-type-specific determinants of mechanosensory properties.

Despite the clinical success of cancer immunotherapies, the majority of patients fail to respond or develop resistance through disruption of pathways that promote neo-antigen presentation on MHC I molecules. Here, we conducted a series of unbiased, genome-wide CRISPR/Cas9 screens to identify genes that limit natural killer (NK) cell anti-tumor activity. We identified that genes associated with antigen presentation and/or interferon-γ (IFN-γ) signaling protect tumor cells from NK cell killing. Indeed, Jak1-deficient melanoma cells were sensitized to NK cell killing through attenuated NK cell-derived IFN-γ-driven transcriptional events that regulate MHC I expression. Importantly, tumor cells that became resistant to T cell killing through enrichment of MHC I-deficient clones were highly sensitive to NK cell killing. Taken together, we reveal the genes targeted by tumor cells to drive checkpoint blockade resistance but simultaneously increase their vulnerability to NK cells, unveiling NK cell-based immunotherapies as a strategy to antagonize tumor immune escape.

Vitamin-D-binding protein (DBP) or group-specific component of serum (GC-globulin) carries vitamin D metabolites from the circulation to target tissues. DBP is highly localized to the liver and pancreatic α cells. Although DBP serum levels, gene polymorphisms, and autoantigens have all been associated with diabetes risk, the underlying mechanisms remain unknown. Here, we show that DBP regulates α cell morphology, α cell function, and glucagon secretion. Deletion of DBP leads to smaller and hyperplastic α cells, altered Na+ channel conductance, impaired α cell activation by low glucose, and reduced rates of glucagon secretion both in vivo and in vitro. Mechanistically, this involves reversible changes in islet microfilament abundance and density, as well as changes in glucagon granule distribution. Defects are also seen in β cell and δ cell function. Immunostaining of human pancreata reveals generalized loss of DBP expression as a feature of late-onset and long-standing, but not early-onset, type 1 diabetes. Thus, DBP regulates α cell phenotype, with implications for diabetes pathogenesis.

Translation initiation is often attributed as the rate-determining step of eukaryotic protein synthesis and key to gene expression control. Despite this centrality, the series of steps involved in this process is poorly understood. Here, we capture the transcriptome-wide occupancy of ribosomes across all stages of translation initiation, enabling us to characterize the transcriptome-wide dynamics of ribosome recruitment to mRNAs, scanning across 5' UTRs and stop codon recognition, in a higher eukaryote. We provide mechanistic evidence for ribosomes attaching to the mRNA by threading the mRNA through the small subunit. Moreover, we identify features that regulate the recruitment and processivity of scanning ribosomes and redefine optimal initiation contexts. Our approach enables deconvoluting translation initiation into separate stages and identifying regulators at each step.

Cryptosporidiosis is a leading cause of life-threatening diarrhea in young children in resource-poor settings. To explore microbial influences on susceptibility, we screened 85 microbiota-associated metabolites for their effects on Cryptosporidium parvum growth in vitro. We identify eight inhibitory metabolites in three main classes: secondary bile salts/acids, a vitamin B6 precursor, and indoles. Growth restriction of C. parvum by indoles does not depend on the host aryl hydrocarbon receptor (AhR) pathway. Instead, treatment impairs host mitochondrial function and reduces total cellular ATP, as well as directly reducing the membrane potential in the parasite mitosome, a degenerate mitochondria. Oral administration of indoles, or reconstitution of the gut microbiota with indole-producing bacteria, delays life cycle progression of the parasite in vitro and reduces the severity of C. parvum infection in mice. Collectively, these findings indicate that microbiota metabolites impair mitochondrial function and contribute to colonization resistance to Cryptosporidium infection.

The discovery of cancer-associated alterations has primarily focused on genetic variants. Nonetheless, altered epigenomes contribute to deregulate transcription and promote oncogenic pathways. Here, we designed an algorithmic approach (RESET) to identify aberrant DNA methylation and associated cis-transcriptional changes across >6,000 human tumors. Tumors exhibiting mutations of chromatin remodeling factors and Wnt signaling displayed DNA methylation instability, characterized by numerous hyper- and hypo-methylated loci. Most silenced and enhanced genes coalesced in specific pathways including apoptosis, DNA repair, and cell metabolism. Cancer-germline antigens (CG) were frequently epigenomically enhanced and their expression correlated with response to anti-PD-1, but not anti-CTLA4, in skin melanoma. Finally, we demonstrated the potential of our approach to explore DNA methylation changes in pediatric tumors, which frequently lack genetic drivers and exhibit epigenomic modifications. Our results provide a pan-cancer map of aberrant DNA methylation to inform functional and therapeutic studies.

The circadian clock operates as intrinsic time-keeping machinery to preserve homeostasis in response to the changing environment. While food is a known zeitgeber for clocks in peripheral tissues, it remains unclear how lack of food influences clock function. We demonstrate that the transcriptional response to fasting operates through molecular mechanisms that are distinct from time-restricted feeding regimens. First, fasting affects core clock genes and proteins, resulting in blunted rhythmicity of BMAL1 and REV-ERBα both in liver and skeletal muscle. Second, fasting induces a switch in temporal gene expression through dedicated fasting-sensitive transcription factors such as GR, CREB, FOXO, TFEB, and PPARs. Third, the rhythmic genomic response to fasting is sustainable by prolonged fasting and reversible by refeeding. Thus, fasting imposes specialized dynamics of transcriptional coordination between the clock and nutrient-sensitive pathways, thereby achieving a switch to fasting-specific temporal gene regulation.

The upper gastrointestinal tract, consisting of the esophagus, stomach, and duodenum, controls food transport, digestion, nutrient uptake, and hormone production. By single-cell analysis of healthy epithelia of these human organs, we molecularly define their distinct cell types. We identify a quiescent COL17A1high KRT15high stem/progenitor cell population in the most basal cell layer of the esophagus and detect substantial gene expression differences between identical cell types of the human and mouse stomach. Selective expression of BEST4, CFTR, guanylin, and uroguanylin identifies a rare duodenal cell type, referred to as BCHE cell, which likely mediates high-volume fluid secretion because of continual activation of the CFTR channel by guanylin/uroguanylin-mediated autocrine signaling. Serotonin-producing enterochromaffin cells in the antral stomach significantly differ in gene expression from duodenal enterochromaffin cells. We, furthermore, discover that the histamine-producing enterochromaffin-like cells in the oxyntic stomach express the luteinizing hormone, yet another member of the enteroendocrine hormone family.

Craniofacial abnormalities often involve sutures, the growth centers of the skull. To characterize the organization and processes governing their development, we profile the murine frontal suture, a model for sutural growth and fusion, at the tissue- and single-cell level on embryonic days (E)16.5 and E18.5. For the wild-type suture, bulk RNA sequencing (RNA-seq) analysis identifies mesenchyme-, osteogenic front-, and stage-enriched genes and biological processes, as well as alternative splicing events modifying the extracellular matrix. Single-cell RNA-seq analysis distinguishes multiple subpopulations, of which five define a mesenchyme-osteoblast differentiation trajectory and show variation along the anteroposterior axis. Similar analyses of in vivo mouse models of impaired frontal suturogenesis in Saethre-Chotzen and Apert syndromes, Twist1+/- and Fgfr2+/S252W, demonstrate distinct transcriptional changes involving angiogenesis and ribogenesis, respectively. Co-expression network analysis reveals gene expression modules from which we validate key driver genes regulating osteoblast differentiation. Our study provides a global approach to gain insights into suturogenesis.

Genome-wide association studies (GWASs) have identified hundreds of loci associated with psychiatric diseases, yet there is a lack of understanding of disease pathophysiology. Common risk variants can shed light on the underlying molecular mechanisms; however, identifying causal variants remains challenging. We map cis-regulatory elements in human neurons derived from pluripotent stem cells. This system allows us to determine enhancers that activate the transcription of neuronal activity-regulated gene programs, which are thought to be critical for synaptic plasticity and are not possible to identify from postmortem tissues. Using the activity-by-contact model, we create variant-to-gene maps to interpret the function of GWAS variants. Our work nominates a subset of variants to elucidate the molecular mechanisms involving GWAS-significant loci. It also highlights that in vitro human cellular models are a powerful platform for identifying and mechanistic studies of human trait-associated genetic variants in cell states that are inaccessible from other types of human samples.

Extensive remodeling of the female mammary epithelium during development and pregnancy has been linked to cancer susceptibility. The faithful response of mammary epithelial cells (MECs) to hormone signaling is key to avoiding breast cancer development. Here, we show that lactogenic differentiation of murine MECs requires silencing of genes encoding ribosomal RNA (rRNA) by the antisense transcript PAPAS. Accordingly, knockdown of PAPAS derepresses rRNA genes, attenuates the response to lactogenic hormones, and induces malignant transformation. Restoring PAPAS levels in breast cancer cells reduces tumorigenicity and lung invasion and activates many interferon-regulated genes previously linked to metastasis suppression. Mechanistically, PAPAS transcription depends on R-loop formation at the 3' end of rRNA genes, which is repressed by RNase H1 and replication protein A (RPA) overexpression in breast cancer cells. Depletion of PAPAS and upregulation of RNase H1 and RPA in human breast cancer underpin the clinical relevance of our findings.

The liver is a key regulator of systemic energy homeostasis whose proper function is dependent on the circadian clock. Here, we show that livers deficient in the oscillator component JARID1a exhibit a dysregulation of genes involved in energy metabolism. Importantly, we find that mice that lack hepatic JARID1a have decreased lean body mass, decreased respiratory exchange ratios, faster production of ketones, and increased glucose production in response to fasting. Finally, we find that JARID1a loss compromises the response of the hepatic transcriptome to nutrient availability. In all, ablation of hepatic JARID1a disrupts the coordination of hepatic metabolic programs with whole-body consequences.

EWSR1-FLI1, the chimeric oncogene specific for Ewing sarcoma (EwS), induces a cascade of signaling events leading to cell transformation. However, it remains elusive how genetically homogeneous EwS cells can drive the heterogeneity of transcriptional programs. Here, we combine independent component analysis of single-cell RNA sequencing data from diverse cell types and model systems with time-resolved mapping of EWSR1-FLI1 binding sites and of open chromatin regions to characterize dynamic cellular processes associated with EWSR1-FLI1 activity. We thus define an exquisitely specific and direct enhancer-driven EWSR1-FLI1 program. In EwS tumors, cell proliferation and strong oxidative phosphorylation metabolism are associated with a well-defined range of EWSR1-FLI1 activity. In contrast, a subpopulation of cells from below and above the intermediary EWSR1-FLI1 activity is characterized by increased hypoxia. Overall, our study reveals sources of intratumoral heterogeneity within EwS tumors.

Mitochondrial oxidative phosphorylation (OXPHOS) and substrate utilization critically regulate the function of hypothalamic proopiomelanocortin (POMC)-expressing neurons. Here, we demonstrate that inactivation of apoptosis-inducing factor (AIF) in POMC neurons mildly impairs mitochondrial respiration and decreases firing of POMC neurons in lean mice. In contrast, under diet-induced obese conditions, POMC-Cre-specific inactivation of AIF prevents obesity-induced silencing of POMC neurons, translating into improved glucose metabolism, improved leptin, and insulin sensitivity, as well as increased energy expenditure in AIFΔPOMC mice. On a cellular level, AIF deficiency improves mitochondrial morphology, facilitates the utilization of fatty acids for mitochondrial respiration, and increases reactive oxygen species (ROS) formation in POMC neurons from obese mice, ultimately leading to restored POMC firing upon HFD feeding. Collectively, partial impairment of mitochondrial function shifts substrate utilization of POMC neurons from glucose to fatty acid metabolism and restores their firing properties, resulting in improved systemic glucose and energy metabolism in obesity.

Many proteins that are needed for progression through S-phase are produced from transcripts that peak in the S-phase, linking temporal expression of those proteins to the time that they are required in cell cycle. Here, we explore the potential roles of long non-coding RNAs in cell cycle progression. We use a sensitive click-chemistry approach to isolate nascent RNAs in a human cell line, and we identify more than 900 long non-coding RNAs (lncRNAs) whose synthesis peaks during the S-phase. More than 200 of these are long intergenic non-coding RNAs (lincRNAs) with S-phase-specific expression. We characterize three of these lincRNAs by knockdown and find that all three lincRNAs are required for appropriate S-phase progression. We infer that non-coding RNAs are key regulatory effectors during the cell cycle, acting on distinct regulatory networks, and herein, we provide a large catalog of candidate cell-cycle regulatory RNAs.

Itch is an unpleasant skin sensation that can be triggered by exposure to many chemicals, including those released by mast cells. The natriuretic polypeptide b (Nppb)-expressing class of sensory neurons, when activated, elicits scratching responses in mice, but it is unclear which itch-inducing agents stimulate these cells and the receptors involved. Here, we identify receptors expressed by Nppb neurons and demonstrate the functional importance of these receptors as sensors of endogenous pruritogens released by mast cells. Our search for receptors in Nppb neurons reveals that they express leukotriene, serotonin, and sphingosine-1-phosphate receptors. Targeted cell ablation, calcium imaging of primary sensory neurons, and conditional receptor knockout studies demonstrate that these receptors induce itch by the direct stimulation of Nppb neurons and neurotransmission through the canonical gastrin-releasing peptide (GRP)-dependent spinal cord itch pathway. Together, our results define a molecular and cellular pathway for mast cell-induced itch.

Cerebral small vessel disease and brain white matter injury are worsened by cardiovascular risk factors including obesity. Molecular pathways in cerebral endothelial cells activated by chronic cerebrovascular risk factors alter cell-cell signaling, blocking endogenous and post-ischemic white matter repair. Using cell-specific translating ribosome affinity purification (RiboTag) in white matter endothelia and oligodendrocyte progenitor cells (OPCs), we identify a coordinated interleukin-chemokine signaling cascade within the oligovascular niche of subcortical white matter that is triggered by diet-induced obesity (DIO). DIO induces interleukin-17B (IL-17B) signaling that acts on the cerebral endothelia through IL-17Rb to increase both circulating and local endothelial expression of CXCL5. In white matter endothelia, CXCL5 promotes the association of OPCs with the vasculature and triggers OPC gene expression programs regulating cell migration through chemokine signaling. Targeted blockade of IL-17B reduced vessel-associated OPCs by reducing endothelial CXCL5 expression. In multiple human cohorts, blood levels of CXCL5 function as a diagnostic and prognostic biomarker of vascular cognitive impairment.

A major goal of current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine efforts is to elicit antibody responses that confer protection. Mapping the epitope targets of the SARS-CoV-2 antibody response is critical for vaccine design, diagnostics, and development of therapeutics. Here, we develop a pan-coronavirus phage display library to map antibody binding sites at high resolution within the complete viral proteomes of all known human-infecting coronaviruses in patients with mild or moderate/severe coronavirus disease 2019 (COVID-19). We find that the majority of immune responses to SARS-CoV-2 are targeted to the spike protein, nucleocapsid, and ORF1ab and include sites of mutation in current variants of concern. Some epitopes are identified in the majority of samples, while others are rare, and we find variation in the number of epitopes targeted between individuals. We find low levels of SARS-CoV-2 cross-reactivity in individuals with no exposure to the virus and significant cross-reactivity with endemic human coronaviruses (CoVs) in convalescent sera from patients with COVID-19.

Somatic DNA copy number variations (CNVs) are prevalent in cancer and can drive cancer progression, albeit with often uncharacterized roles in altering cell signaling states. Here, we integrate genomic and proteomic data for 5,598 tumor samples to identify CNVs leading to aberrant signal transduction. The resulting associations recapitulate known kinase-substrate relationships, and further network analysis prioritizes likely causal genes. Of the 303 significant associations we identify from the pan-tumor analysis, 43% are replicated in cancer cell lines, including 44 robust gene-phosphosite associations identified across multiple tumor types. Several predicted regulators of hippo signaling are experimentally validated. Using RNAi, CRISPR, and drug screening data, we find evidence of kinase addiction in cancer cell lines, identifying inhibitors for targeting of kinase-dependent cell lines. We propose copy number status of genes as a useful predictor of differential impact of kinase inhibition, a strategy that may be of use in the future for anticancer therapies.