Neuronal Computation Underlying Inferential Reasoning in Humans and Mice.

Every day we make decisions critical for adaptation and survival. We repeat actions with known consequences. But we also draw on loosely related events to infer and imagine the outcome of entirely novel choices. These inferential decisions are thought to engage a number of brain regions; however, the underlying neuronal computation remains unknown. Here, we use a multi-day cross-species approach in humans and mice to report the functional anatomy and neuronal computation underlying inferential decisions. We show that during successful inference, the mammalian brain uses a hippocampal prospective code to forecast temporally structured learned associations. Moreover, during resting behavior, coactivation of hippocampal cells in sharp-wave/ripples represent inferred relationships that include reward, thereby "joining-the-dots" between events that have not been observed together but lead to profitable outcomes. Computing mnemonic links in this manner may provide an important mechanism to build a cognitive map that stretches beyond direct experience, thus supporting flexible behavior.

Profiling of Pluripotency Factors in Single Cells and Early Embryos.

Cell fate decisions are governed by sequence-specific transcription factors (TFs) that act in small populations of cells within developing embryos. To understand their functions in vivo, it is important to identify TF binding sites in these cells. However, current methods cannot profile TFs genome-wide at or near the single-cell level. Here we adapt the cleavage under targets and release using nuclease (CUT&RUN) method to profile TFs in low cell numbers, including single cells and individual pre-implantation embryos. Single-cell experiments suggest that only a fraction of TF binding sites are occupied in most cells, in a manner broadly consistent with measurements of peak intensity from multi-cell studies. We further show that chromatin binding by the pluripotency TF NANOG is highly dependent on the SWI/SNF chromatin remodeling complex in individual blastocysts but not in cultured cells. Ultra-low input CUT&RUN (uliCUT&RUN) therefore enables interrogation of TF binding from rare cell populations of particular importance in development or disease.

Direct Promoter Repression by BCL11A Controls the Fetal to Adult Hemoglobin Switch.

Fetal hemoglobin (HbF, α2γ2) level is genetically controlled and modifies severity of adult hemoglobin (HbA, α2ß2) disorders, sickle cell disease, and ß-thalassemia. Common genetic variation affects expression of BCL11A, a regulator of HbF silencing. To uncover how BCL11A supports the developmental switch from γ- to ß- globin, we use a functional assay and protein binding microarray to establish a requirement for a zinc-finger cluster in BCL11A in repression and identify a preferred DNA recognition sequence. This motif appears in embryonic and fetal-expressed globin promoters and is duplicated in γ-globin promoters. The more distal of the duplicated motifs is mutated in individuals with hereditary persistence of HbF. Using the CUT&RUN approach to map protein binding sites in erythroid cells, we demonstrate BCL11A occupancy preferentially at the distal motif, which can be disrupted by editing the promoter. Our findings reveal that direct γ-globin gene promoter repression by BCL11A underlies hemoglobin switching.

Retinoid X receptor gamma dictates the activation threshold of group 2 innate lymphoid cells and limits type 2 inflammation in the small intestine.

Group 2 innate lymphoid cells (ILC2s) are crucial in promoting type 2 inflammation that contributes to both anti-parasite immunity and allergic diseases. However, the molecular checkpoints in ILC2s that determine whether to immediately launch a proinflammatory response are unknown. Here, we found that retinoid X receptor gamma (Rxrg) was highly expressed in small intestinal ILC2s and rapidly suppressed by alarmin cytokines. Genetic deletion of Rxrg did not impact ILC2 development but facilitated ILC2 responses and the tissue inflammation induced by alarmins. Mechanistically, RXRγ maintained the expression of its target genes that support intracellular cholesterol efflux, which in turn reduce ILC2 proliferation. Furthermore, RXRγ expression prevented ILC2 response to mild stimulations, including low doses of alarmin cytokine and mechanical skin injury. Together, we propose that RXRγ expression and its mediated lipid metabolic states function as a cell-intrinsic checkpoint that confers the threshold of ILC2 activation in the small intestine.

Integrated multi-omics analysis of zinc-finger proteins uncovers roles in RNA regulation.

RNA interactome studies have revealed that hundreds of zinc-finger proteins (ZFPs) are candidate RNA-binding proteins (RBPs), yet their RNA substrates and functional significance remain largely uncharacterized. Here, we present a systematic multi-omics analysis of the DNA- and RNA-binding targets and regulatory roles of more than 100 ZFPs representing 37 zinc-finger families. We show that multiple ZFPs are previously unknown regulators of RNA splicing, alternative polyadenylation, stability, or translation. The examined ZFPs show widespread sequence-specific RNA binding and preferentially bind proximal to transcription start sites. Additionally, several ZFPs associate with their targets at both the DNA and RNA levels. We highlight ZNF277, a C2H2 ZFP that binds thousands of RNA targets and acts as a multi-functional RBP. We also show that ZNF473 is a DNA/RNA-associated protein that regulates the expression and splicing of cell cycle genes. Our results reveal diverse roles for ZFPs in transcriptional and post-transcriptional gene regulation.

Single-nucleus multiomic mapping of m6A methylomes and transcriptomes in native populations of cells with sn-m6A-CT.

N6-methyladenosine (m6A) RNA modification plays important roles in the governance of gene expression and is temporally regulated in different cell states. In contrast to global m6A profiling in bulk sequencing, single-cell technologies for analyzing m6A heterogeneity are not extensively established. Here, we developed single-nucleus m6A-CUT&Tag (sn-m6A-CT) for simultaneous profiling of m6A methylomes and transcriptomes within a single nucleus using mouse embryonic stem cells (mESCs). m6A-CT is capable of enriching m6A-marked RNA molecules in situ, without isolating RNAs from cells. We adapted m6A-CT to the droplet-based single-cell omics platform and demonstrated high-throughput performance in analyzing nuclei isolated from thousands of cells from various cell types. We show that sn-m6A-CT profiling is sufficient to determine cell identity and allows the generation of cell-type-specific m6A methylome landscapes from heterogeneous populations. These indicate that sn-m6A-CT provides additional dimensions to multimodal datasets and insights into epitranscriptomic landscape in defining cell fate identity and states.

Stepwise activities of mSWI/SNF family chromatin remodeling complexes direct T cell activation and exhaustion.

Highly coordinated changes in gene expression underlie T cell activation and exhaustion. However, the mechanisms by which such programs are regulated and how these may be targeted for therapeutic benefit remain poorly understood. Here, we comprehensively profile the genomic occupancy of mSWI/SNF chromatin remodeling complexes throughout acute and chronic T cell stimulation, finding that stepwise changes in localization over transcription factor binding sites direct site-specific chromatin accessibility and gene activation leading to distinct phenotypes. Notably, perturbation of mSWI/SNF complexes using genetic and clinically relevant chemical strategies enhances the persistence of T cells with attenuated exhaustion hallmarks and increased memory features in vitro and in vivo. Finally, pharmacologic mSWI/SNF inhibition improves CAR-T expansion and results in improved anti-tumor control in vivo. These findings reveal the central role of mSWI/SNF complexes in the coordination of T cell activation and exhaustion and nominate small-molecule-based strategies for the improvement of current immunotherapy protocols.

PALI1 promotes tumor growth through competitive recruitment of PRC2 to G9A-target chromatin for dual epigenetic silencing.

PALI1 is a newly identified accessory protein of the Polycomb repressive complex 2 (PRC2) that catalyzes H3K27 methylation. However, the roles of PALI1 in cancer are yet to be defined. Here, we report that PALI1 is upregulated in advanced prostate cancer (PCa) and competes with JARID2 for binding to the PRC2 core subunit SUZ12. PALI1 further interacts with the H3K9 methyltransferase G9A, bridging the formation of a unique G9A-PALI1-PRC2 super-complex that occupies a subset of G9A-target genes to mediate dual H3K9/K27 methylation and gene repression. Many of these genes are developmental regulators required for cell differentiation, and their loss in PCa predicts poor prognosis. Accordingly, PALI1 and G9A drive PCa cell proliferation and invasion in vitro and xenograft tumor growth in vivo. Collectively, our study shows that PALI1 harnesses two central epigenetic mechanisms to suppress cellular differentiation and promote tumorigenesis, which can be targeted by dual EZH2 and G9A inhibition.

A unified atlas of CD8 T cell dysfunctional states in cancer and infection.

CD8 T cells play an essential role in defense against viral and bacterial infections and in tumor immunity. Deciphering T cell loss of functionality is complicated by the conspicuous heterogeneity of CD8 T cell states described across experimental and clinical settings. By carrying out a unified analysis of over 300 assay for transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-seq) experiments from 12 studies of CD8 T cells in cancer and infection, we defined a shared differentiation trajectory toward dysfunction and its underlying transcriptional drivers and revealed a universal early bifurcation of functional and dysfunctional T cell states across models. Experimental dissection of acute and chronic viral infection using single-cell ATAC (scATAC)-seq and allele-specific single-cell RNA (scRNA)-seq identified state-specific drivers and captured the emergence of similar TCF1+ progenitor-like populations at an early branch point, at which functional and dysfunctional T cells diverge. Our atlas of CD8 T cell states will facilitate mechanistic studies of T cell immunity and translational efforts.

Conserved Gsx2/Ind homeodomain monomer versus homodimer DNA binding defines regulatory outcomes in flies and mice.

How homeodomain proteins gain sufficient specificity to control different cell fates has been a long-standing problem in developmental biology. The conserved Gsx homeodomain proteins regulate specific aspects of neural development in animals from flies to mammals, and yet they belong to a large transcription factor family that bind nearly identical DNA sequences in vitro. Here, we show that the mouse and fly Gsx factors unexpectedly gain DNA binding specificity by forming cooperative homodimers on precisely spaced and oriented DNA sites. High-resolution genomic binding assays revealed that Gsx2 binds both monomer and homodimer sites in the developing mouse ventral telencephalon. Importantly, reporter assays showed that Gsx2 mediates opposing outcomes in a DNA binding site-dependent manner: Monomer Gsx2 binding represses transcription, whereas homodimer binding stimulates gene expression. In Drosophila, the Gsx homolog, Ind, similarly represses or stimulates transcription in a site-dependent manner via an autoregulatory enhancer containing a combination of monomer and homodimer sites. Integrating these findings, we test a model showing how the homodimer to monomer site ratio and the Gsx protein levels defines gene up-regulation versus down-regulation. Altogether, these data serve as a new paradigm for how cooperative homeodomain transcription factor binding can increase target specificity and alter regulatory outcomes.

Opposing Functions of BRD4 Isoforms in Breast Cancer.

Bromodomain-containing protein 4 (BRD4) is a cancer therapeutic target in ongoing clinical trials disrupting primarily BRD4-regulated transcription programs. The role of BRD4 in cancer has been attributed mainly to the abundant long isoform (BRD4-L). Here we show, by isoform-specific knockdown and endogenous protein detection, along with transgene expression, the less abundant BRD4 short isoform (BRD4-S) is oncogenic while BRD4-L is tumor-suppressive in breast cancer cell proliferation and migration, as well as mammary tumor formation and metastasis. Through integrated RNA-seq, genome-wide ChIP-seq, and CUT&RUN association profiling, we identify the Engrailed-1 (EN1) homeobox transcription factor as a key BRD4-S coregulator, particularly in triple-negative breast cancer. BRD4-S and EN1 comodulate the extracellular matrix (ECM)-associated matrisome network, including type II cystatin gene cluster, mucin 5, and cathepsin loci, via enhancer regulation of cancer-associated genes and pathways. Our work highlights the importance of targeted therapies for the oncogenic, but not tumor-suppressive, activity of BRD4.

RSC-Associated Subnucleosomes Define MNase-Sensitive Promoters in Yeast.

The classic view of nucleosome organization at active promoters is that two well-positioned nucleosomes flank a nucleosome-depleted region (NDR). However, this view has been recently disputed by contradictory reports as to whether wider (≳150 bp) NDRs instead contain unstable, micrococcal nuclease-sensitive ("fragile") nucleosomal particles. To determine the composition of fragile particles, we introduce CUT&RUN.ChIP, in which targeted nuclease cleavage and release is followed by chromatin immunoprecipitation. We find that fragile particles represent the occupancy of the RSC (remodeling the structure of chromatin) nucleosome remodeling complex and RSC-bound, partially unwrapped nucleosomal intermediates. We also find that general regulatory factors (GRFs) bind to partially unwrapped nucleosomes at these promoters. We propose that RSC binding and its action cause nucleosomes to unravel, facilitate subsequent binding of GRFs, and constitute a dynamic cycle of nucleosome deposition and clearance at the subset of wide Pol II promoter NDRs.

Nitrogen availability is important for preventing catastrophic mitosis in fission yeast.

Mitosis is a crucial stage in the cell cycle, controlled by a vast network of regulators responding to multiple internal and external factors. The fission yeast Schizosaccharomyces pombe demonstrates catastrophic mitotic phenotypes due to mutations or drug treatments. One of the factors provoking catastrophic mitosis is a disturbed lipid metabolism, resulting from, for example, mutations in the acetyl-CoA/biotin carboxylase (cut6), fatty acid synthase (fas2, also known as lsd1) or transcriptional regulator of lipid metabolism (cbf11) genes, as well as treatment with inhibitors of fatty acid synthesis. It has been previously shown that mitotic fidelity in lipid metabolism mutants can be partially rescued by ammonium chloride supplementation. In this study, we demonstrate that mitotic fidelity can be improved by multiple nitrogen sources. Moreover, this improvement is not limited to lipid metabolism disturbances but also applies to a number of unrelated mitotic mutants. Interestingly, the partial rescue is not achieved by restoring the lipid metabolism state, but rather indirectly. Our results highlight a novel role for nitrogen availability in mitotic fidelity.

The CUT&RUN greenlist: genomic regions of consistent noise are effective normalizing factors for quantitative epigenome mapping.

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a recent development for epigenome mapping, but its unique methodology can hamper proper quantitative analyses. As traditional normalization approaches have been shown to be inaccurate, we sought to determine endogenous normalization factors based on the human genome regions of constant nonspecific signal. This constancy was determined by applying Shannon's information entropy, and the set of normalizer regions, which we named the 'Greenlist', was extensively validated using publicly available datasets. We demonstrate here that the greenlist normalization outperforms the current top standards, and remains consistent across different experimental setups, cell lines and antibodies; the approach can even be applied to different species or to CUT&Tag. Requiring no additional experimental steps and no added cost, this approach can be universally applied to CUT&RUN experiments to greatly minimize the interference of technical variation over the biological epigenome changes of interest.

Review and meta-analysis of the genetic Minimal Cut Set approach for gene essentiality prediction in cancer metabolism.

Cancer metabolism is a marvellously complex topic, in part, due to the reprogramming of its pathways to self-sustain the malignant phenotype in the disease, to the detriment of its healthy counterpart. Understanding these adjustments can provide novel targeted therapies that could disrupt and impair proliferation of cancerous cells. For this very purpose, genome-scale metabolic models (GEMs) have been developed, with Human1 being the most recent reconstruction of the human metabolism. Based on GEMs, we introduced the genetic Minimal Cut Set (gMCS) approach, an uncontextualized methodology that exploits the concepts of synthetic lethality to predict metabolic vulnerabilities in cancer. gMCSs define a set of genes whose knockout would render the cell unviable by disrupting an essential metabolic task in GEMs, thus, making cellular proliferation impossible. Here, we summarize the gMCS approach and review the current state of the methodology by performing a systematic meta-analysis based on two datasets of gene essentiality in cancer. First, we assess several thresholds and distinct methodologies for discerning highly and lowly expressed genes. Then, we address the premise that gMCSs of distinct length should have the same predictive power. Finally, we question the importance of a gene partaking in multiple gMCSs and analyze the importance of all the essential metabolic tasks defined in Human1. Our meta-analysis resulted in parameter evaluation to increase the predictive power for the gMCS approach, as well as a significant reduction of computation times by only selecting the crucial gMCS lengths, proposing the pertinency of particular parameters for the peak processing of gMCS.

A single-cell multiomic analysis of kidney organoid differentiation.

Kidney organoids differentiated from pluripotent stem cells are powerful models of kidney development and disease but are characterized by cell immaturity and off-target cell fates. Comparing the cell-specific gene regulatory landscape during organoid differentiation with human adult kidney can serve to benchmark progress in differentiation at the epigenome and transcriptome level for individual organoid cell types. Using single-cell multiome and histone modification analysis, we report more broadly open chromatin in organoid cell types compared to the human adult kidney. We infer enhancer dynamics by cis-coaccessibility analysis and validate an enhancer driving transcription of HNF1B by CRISPR interference both in cultured proximal tubule cells and also during organoid differentiation. Our approach provides an experimental framework to judge the cell-specific maturation state of human kidney organoids and shows that kidney organoids can be used to validate individual gene regulatory networks that regulate differentiation.

The kinesin-5 protein Cut7 moves bidirectionally on fission yeast spindles with activity that increases in anaphase.

Kinesin-5 motors are essential to separate mitotic spindle poles and assemble a bipolar spindle in many organisms. These motors crosslink and slide apart antiparallel microtubules via microtubule plus-end-directed motility. However, kinesin-5 localization is enhanced away from antiparallel overlaps. Increasing evidence suggests this localization occurs due to bidirectional motility or trafficking. The purified fission-yeast kinesin-5 protein Cut7 moves bidirectionally, but bidirectionality has not been shown in cells, and the function of the minus-end-directed movement is unknown. Here, we characterized the motility of Cut7 on bipolar and monopolar spindles and observed movement toward both plus- and minus-ends of microtubules. Notably, the activity of the motor increased at anaphase B onset. Perturbations to microtubule dynamics only modestly changed Cut7 movement, whereas Cut7 mutation reduced movement. These results suggest that the directed motility of Cut7 contributes to the movement of the motor. Comparison of the Cut7 mutant and human Eg5 (also known as KIF11) localization suggest a new hypothesis for the function of minus-end-directed motility and spindle-pole localization of kinesin-5s.

A new CUT&RUN low volume-urea (LoV-U) protocol optimized for transcriptional co-factors uncovers Wnt/ß-catenin tissue-specific genomic targets.

Upon WNT/ß-catenin pathway activation, stabilized ß-catenin travels to the nucleus where it associates with the TCF/LEF transcription factors, constitutively bound to genomic Wnt-responsive elements (WREs), to activate target gene transcription. Discovering the binding profile of ß-catenin is therefore required to unambiguously assign direct targets of WNT signaling. Cleavage under targets and release using nuclease (CUT&RUN) has emerged as prime technique for mapping the binding profile of DNA-interacting proteins. Here, we present a modified version of CUT&RUN, named LoV-U (low volume and urea), that enables the robust and reproducible generation of ß-catenin binding profiles, uncovering direct WNT/ß-catenin target genes in human cells, as well as in cells isolated from developing mouse tissues. CUT&RUN-LoV-U outperforms original CUT&RUN when targeting co-factors that do not bind the DNA, can profile all classes of chromatin regulators and is well suited for simultaneous processing of several samples. We believe that the application of our protocol will allow the detection of the complex system of tissue-specific WNT/ß-catenin target genes, together with other non-DNA-binding transcriptional regulators that act downstream of ontogenetically fundamental signaling cascades.

Discovery and characterization of sgRNA-sequence-independent DNA cleavage from CRISPR/Cas9 in mouse embryos.

The CRISPR/Cas9 system can induce off-target effects in programmed gene editing, but there have been few reports on cleavage detection and their affection in embryo development. To study these events, sgRNAs with different off-target rates were designed and compared after micro-injected into mouse zygotes, and γH2AX was used for DNA cleavage sites analysis by immunostaining and CUT&Tag. Although the low off-target sgRNA were usually selected for production gene editing animals, γH2AX immunofluorescence indicated that there was a relative DSB peak at 15 h after Cas9 system injection, and the number of γH2AX foci at the peak was significantly higher in the low off-target sgRNA-injected group than in the control group. Further, the result of CUT&Tag sequencing analysis showed more double-strand breaks (DSBs) related sequences were detected in low off-target sgRNA-injected group than control and the distribution of DSB related sequences had no chromosome specificity. Gene Ontology (GO) annotation analysis of the DSB related sequences showed that these sequences were mainly concentrated at genes associated with some important biological processes, molecular functions, and cell components. In a conclusion, there are many sgRNA-sequence-independent DSBs in early mouse embryos when the Cas9 system is used for gene editing and the DSB related sequence could be detected and characterized in the genome. These results and method should also be considered in using or optimizing the Cas9 system.

Wild-type distributions of minimum inhibitory concentrations and epidemiological cut-off values-laboratory and clinical utility.

The characterization of wild-type minimum inhibitory concentration (MIC) and zone diameter distributions with the setting of epidemiological cut-off values (ECOFFs or ECVs) provides a reference for the otherwise relative MIC values in the international system for antimicrobial susceptibility testing. Distributions of MIC values for a species and an agent follow a log-normal distribution, which in the absence of resistance mechanisms is monomodal and designated wild type (WT). The upper end of the WT distribution, the ECOFF, can be identified with statistical methods. In the presence of phenotypically detectable resistance, the distribution has at least one more mode (the non-WT), but despite this, the WT is most often identifiable using the same methods. The ECOFF provides the most sensitive measure of resistance development in a species against an agent. The WT and non-WT modes are independent of the organism´s response to treatment, but when the European Committee on Antimicrobial Susceptibility Testing (EUCAST) determines the clinical breakpoints, the committee avoids breakpoints that split WT distributions of target species. This is to avoid the poorer reproducibility of susceptibility categorization when breakpoints split major populations but also because the EUCAST has failed to identify different clinical outcomes for isolates with different MIC values inside the wild-type distribution. In laboratory practice, the ECOFF is used to screen for and exclude resistance and allows the comparison of resistance between systems with different breakpoints from different breakpoint organizations, breakpoints evolving over time, and different breakpoints between human and animal medicine. The EUCAST actively encourages colleagues to question MIC distributions as presented on the website (https://www.eucast.org/mic_and_zone_distributions_and_ecoffs) and to contribute MIC and inhibition zone diameter data.