Cnidarian Cell Type Diversity and Regulation Revealed by Whole-Organism Single-Cell RNA-Seq.

The emergence and diversification of cell types is a leading factor in animal evolution. So far, systematic characterization of the gene regulatory programs associated with cell type specificity was limited to few cell types and few species. Here, we perform whole-organism single-cell transcriptomics to map adult and larval cell types in the cnidarian Nematostella vectensis, a non-bilaterian animal with complex tissue-level body-plan organization. We uncover eight broad cell classes in Nematostella, including neurons, cnidocytes, and digestive cells. Each class comprises different subtypes defined by the expression of multiple specific markers. In particular, we characterize a surprisingly diverse repertoire of neurons, which comparative analysis suggests are the result of lineage-specific diversification. By integrating transcription factor expression, chromatin profiling, and sequence motif analysis, we identify the regulatory codes that underlie Nematostella cell-specific expression. Our study reveals cnidarian cell type complexity and provides insights into the evolution of animal cell-specific genomic regulation.

Two independent modes of chromatin organization revealed by cohesin removal.

Imaging and chromosome conformation capture studies have revealed several layers of chromosome organization, including segregation into megabase-sized active and inactive compartments, and partitioning into sub-megabase domains (TADs). It remains unclear, however, how these layers of organization form, interact with one another and influence genome function. Here we show that deletion of the cohesin-loading factor Nipbl in mouse liver leads to a marked reorganization of chromosomal folding. TADs and associated Hi-C peaks vanish globally, even in the absence of transcriptional changes. By contrast, compartmental segregation is preserved and even reinforced. Strikingly, the disappearance of TADs unmasks a finer compartment structure that accurately reflects the underlying epigenetic landscape. These observations demonstrate that the three-dimensional organization of the genome results from the interplay of two independent mechanisms: cohesin-independent segregation of the genome into fine-scale compartments, defined by chromatin state; and cohesin-dependent formation of TADs, possibly by loop extrusion, which helps to guide distant enhancers to their target genes.

De Novo Variants Found in Three Distinct Schizophrenia Populations Hit a Common Core Gene Network Related to Microtubule and Actin Cytoskeleton Gene Ontology Classes.

Schizophrenia (SZ) is a heterogeneous and debilitating psychiatric disorder with a strong genetic component. To elucidate functional networks perturbed in schizophrenia, we analysed a large dataset of whole-genome studies that identified SNVs, CNVs, and a multi-stage schizophrenia genome-wide association study. Our analysis identified three subclusters that are interrelated and with small overlaps: GO:0007017~Microtubule-Based Process, GO:00015629~Actin Cytoskeleton, and GO:0007268~SynapticTransmission. We next analysed three distinct trio cohorts of 75 SZ Algerian, 45 SZ French, and 61 SZ Japanese patients. We performed Illumina HiSeq whole-exome sequencing and identified de novo mutations using a Bayesian approach. We validated 88 de novo mutations by Sanger sequencing: 35 in French, 21 in Algerian, and 32 in Japanese SZ patients. These 88 de novo mutations exhibited an enrichment in genes encoding proteins related to GO:0051015~actin filament binding (p = 0.0011) using David, and enrichments in GO: 0003774~transport (p = 0.019) and GO:0003729~mRNA binding (p = 0.010) using Amigo. One of these de novo variant was found in CORO1C coding sequence. We studied Coro1c haploinsufficiency in a Coro1c+/- mouse and found defects in the corpus callosum. These results could motivate future studies of the mechanisms surrounding genes encoding proteins involved in transport and the cytoskeleton, with the goal of developing therapeutic intervention strategies for a subset of SZ cases.

Single cell atlas of Xenoturbella bocki highlights limited cell-type complexity.

Phylogenetic analyses over the last two decades have united a few small, and previously orphan clades, the nematodermatids, acoels and xenoturbelids, into the phylum Xenacoelomorpha. Some phylogenetic analyses support a sister relationship between Xenacoelomorpha and Ambulacraria (Xenambulacraria), while others suggest that Xenacoelomorpha may be sister to the rest of the Bilateria (Nephrozoa). An understanding of the cell type complements of Xenacoelomorphs is essential to assessing these alternatives as well as to our broader understanding of bilaterian cell type evolution. Employing whole organism single-cell RNA-seq in the marine xenacoelomorph worm Xenoturbella bocki, we show that Xenambulacrarian nerve nets share regulatory features and a peptidergic identity with those found in cnidarians and protostomes and more broadly share muscle and gland cell similarities with other metazoans. Taken together, these data are consistent with broad homologies of animal gland, muscle, and neurons as well as more specific affinities between Xenoturbella and acoel gut and epidermal tissues, consistent with the monophyly of Xenacoelomorpha.

A major role for CD4+ T cells in driving cytokine release syndrome during CAR T cell therapy.

Anti-CD19 chimeric antigen receptor (CAR) T cell therapy represents a breakthrough for the treatment of B cell malignancies. Yet, it can lead to severe adverse events, including cytokine release syndrome (CRS), which may require urgent clinical management. Whether interpatient variability in CAR T cell subsets contributes to CRS is unclear. Here, we show that CD4+ CAR T cells are the main drivers of CRS. Using an immunocompetent model of anti-CD19 CAR T cell therapy, we report that CD4+, but not CD8+, CAR T cells elicit physiological CRS-like manifestations associated with the release of inflammatory cytokines. In CAR T cell-treated patients, CRS occurrence and severity are significantly associated with high absolute values of CD4+ CAR T cells in the blood. CRS in mice occurs independently of CAR T cell-derived interferon γ (IFN-γ) but requires elevated tumor burden. Thus, adjusting the CD4:CD8 CAR T cell ratio to patient tumor load may help mitigate CAR T cell-associated toxicities.

Transient suppression of SUMOylation in embryonic stem cells generates embryo-like structures.

Recent advances in synthetic embryology have opened new avenues for understanding the complex events controlling mammalian peri-implantation development. Here, we show that mouse embryonic stem cells (ESCs) solely exposed to chemical inhibition of SUMOylation generate embryo-like structures comprising anterior neural and trunk-associated regions. HypoSUMOylation-instructed ESCs give rise to spheroids that self-organize into gastrulating structures containing cell types spatially and functionally related to embryonic and extraembryonic compartments. Alternatively, spheroids cultured in a droplet microfluidic device form elongated structures that undergo axial organization reminiscent of natural embryo morphogenesis. Single-cell transcriptomics reveals various cellular lineages, including properly positioned anterior neuronal cell types and paraxial mesoderm segmented into somite-like structures. Transient SUMOylation suppression gradually increases DNA methylation genome wide and repressive mark deposition at Nanog. Interestingly, cell-to-cell variations in SUMOylation levels occur during early embryogenesis. Our approach provides a proof of principle for potentially powerful strategies to explore early embryogenesis by targeting chromatin roadblocks of cell fate change.

SMARCA2 and other genome-wide supported schizophrenia-associated genes: regulation by REST/NRSF, network organization and primate-specific evolution.

The SMARCA2 gene, which encodes BRM in the SWI/SNF chromatin-remodeling complex, was recently identified as being associated with schizophrenia (SZ) in a genome-wide approach. Polymorphisms in SMARCA2, associated with the disease, produce changes in the expression of the gene and/or in the encoded amino acid sequence. We show here that an SWI/SNF-centered network including the Smarca2 gene is modified by the down-regulation of REST/NRSF in a mouse neuronal cell line. REST/NRSF down-regulation also modifies the levels of Smarce1, Smarcd3 and SWI/SNF interactors (Hdac1, RcoR1 and Mecp2). Smarca2 down-regulation generates an abnormal dendritic spine morphology that is an intermediate phenotype of SZ. We further found that 8 (CSF2RA, HIST1H2BJ, NOTCH4, NRGN, SHOX, SMARCA2, TCF4 and ZNF804A) out of 10 genome-wide supported SZ-associated genes are part of an interacting network (including SMARCA2), 5 members of which encode transcription regulators. The expression of 3 (TCF4, SMARCA2 and CSF2RA) of the 10 genome-wide supported SZ-associated genes is modified when the REST/NRSF-SWI/SNF chromatin-remodeling complex is experimentally manipulated in mouse cell lines and in transgenic mouse models. The REST/NRSF-SWI/SNF deregulation also results in the differential expression of genes that are clustered in chromosomes suggesting the induction of genome-wide epigenetic changes. Finally, we found that SMARCA2 interactors and the genome-wide supported SZ-associated genes are considerably enriched in genes displaying positive selection in primates and in the human lineage which suggests the occurrence of novel protein interactions in primates. Altogether, these data identify the SWI/SNF chromatin-remodeling complex as a key component of the genetic architecture of SZ.

TATTOO-seq delineates spatial and cell type-specific regulatory programs in the developing limb.

The coordinated differentiation of progenitor cells into specialized cell types and their spatial organization into distinct domains is central to embryogenesis. Here, we developed and applied an unbiased spatially resolved single-cell transcriptomics method to identify the genetic programs underlying the emergence of specialized cell types during mouse limb development and their spatial integration. We identify multiple transcription factors whose expression patterns are predominantly associated with cell type specification or spatial position, suggesting two parallel yet highly interconnected regulatory systems. We demonstrate that the embryonic limb undergoes a complex multiscale reorganization upon perturbation of one of its spatial organizing centers, including the loss of specific cell populations, alterations of preexisting cell states' molecular identities, and changes in their relative spatial distribution. Our study shows how multidimensional single-cell, spatially resolved molecular atlases can allow the deconvolution of spatial identity and cell fate and reveal the interconnected genetic networks that regulate organogenesis and its reorganization upon genetic alterations.

Chr21 protein-protein interactions: enrichment in proteins involved in intellectual disability, autism, and late-onset Alzheimer's disease.

Down syndrome (DS) is caused by human chromosome 21 (HSA21) trisomy. It is characterized by a poorly understood intellectual disability (ID). We studied two mouse models of DS, one with an extra copy of the <i>Dyrk1A</i> gene (189N3) and the other with an extra copy of the mouse Chr16 syntenic region (Dp(16)1Yey). RNA-seq analysis of the transcripts deregulated in the embryonic hippocampus revealed an enrichment in genes associated with chromatin for the 189N3 model, and synapses for the Dp(16)1Yey model. A large-scale yeast two-hybrid screen (82 different screens, including 72 HSA21 baits and 10 rebounds) of a human brain library containing at least 10⁷ independent fragments identified 1,949 novel protein-protein interactions. The direct interactors of HSA21 baits and rebounds were significantly enriched in ID-related genes (<i>P</i>-value &lt; 2.29 × 10^-8). Proximity ligation assays showed that some of the proteins encoded by HSA21 were located at the dendritic spine postsynaptic density, in a protein network at the dendritic spine postsynapse. We located HSA21 DYRK1A and DSCAM, mutations of which increase the risk of autism spectrum disorder (ASD) 20-fold, in this postsynaptic network. We found that an intracellular domain of DSCAM bound either DLGs, which are multimeric scaffolds comprising receptors, ion channels and associated signaling proteins, or DYRK1A. The DYRK1A-DSCAM interaction domain is conserved in <i>Drosophila</i> and humans. The postsynaptic network was found to be enriched in proteins associated with ARC-related synaptic plasticity, ASD, and late-onset Alzheimer's disease. These results highlight links between DS and brain diseases with a complex genetic basis.

DYRK1A interacts with the REST/NRSF-SWI/SNF chromatin remodelling complex to deregulate gene clusters involved in the neuronal phenotypic traits of Down syndrome.

The molecular mechanisms that lead to the cognitive defects characteristic of Down syndrome (DS), the most frequent cause of mental retardation, have remained elusive. Here we use a transgenic DS mouse model (152F7 line) to show that DYRK1A gene dosage imbalance deregulates chromosomal clusters of genes located near neuron-restrictive silencer factor (REST/NRSF) binding sites. We found that Dyrk1a binds the SWI/SNF complex known to interact with REST/NRSF. The mutation of a REST/NRSF binding site in the promoter of the REST/NRSF target gene L1cam modifies the transcriptional effect of Dyrk1a-dosage imbalance on L1cam. Dyrk1a dosage imbalance perturbs Rest/Nrsf levels with decreased Rest/Nrsf expression in embryonic neurons and increased expression in adult neurons. Using transcriptome analysis of embryonic brain subregions of transgenic 152F7 mouse line, we identified a coordinated deregulation of multiple genes that are responsible for dendritic growth impairment present in DS. Similarly, Dyrk1a overexpression in primary mouse cortical neurons induced severe reduction of the dendritic growth and dendritic complexity. We propose that DYRK1A overexpression-related neuronal gene deregulation via disturbance of REST/NRSF levels, and the REST/NRSF-SWI/SNF chromatin remodelling complex, significantly contributes to the neural phenotypic changes that characterize DS.

A cross-talk between CAR T cell subsets and the tumor microenvironment is essential for sustained cytotoxic activity.

Chimeric antigen receptor (CAR) T cell therapy relies on the activity of a large pool of tumor-targeting cytotoxic effectors. Whether CAR T cells act autonomously or require interactions with the tumor microenvironment (TME) remains incompletely understood. Here, we report an essential cross-talk between CAR T cell subsets and the TME for tumor control in an immunocompetent mouse B cell lymphoma model of anti-CD19 CAR T cell therapy. Using single-cell RNA sequencing, we revealed substantial modification of the TME during CAR T cell therapy. Interferon-γ (IFN-γ) produced by CAR T cells not only enhanced endogenous T and natural killer cell activity but was also essential for sustaining CAR T cell cytotoxicity, as revealed by intravital imaging. CAR T cell-derived IFN-γ facilitated host interleukin-12 production that supported host immune and CAR T cell responses. Compared with CD8+ CAR T cells, CD4+ CAR T cells were more efficient at host immune activation but less capable of direct tumor killing. In summary, CAR T cells do not act independently in vivo but rely instead on cytokine-mediated cross-talk with the TME for optimal activity. Invigorating CAR T cell interplay with the host represents an attractive strategy to prevent relapses after therapy.

Convergent evidence identifying MAP/microtubule affinity-regulating kinase 1 (MARK1) as a susceptibility gene for autism.

Autism spectrum disorders (ASDs) are common, heritable, but genetically heterogeneous neurodevelopmental conditions. We recently defined a susceptibility locus for ASDs on chromosome 1q41-q42. High-resolution single-nucleotide polymorphisms (126 SNPs) genotyping across the chromosome 1q41-q42 region, followed by a MARK1 (microtubule affinity-regulating kinase 1)-tagged-SNP association study in 276 families with autism from the Autism Genetic Research Exchange, showed that several SNPs within the MARK1 gene were significantly associated with ASDs by transmission disequilibrium tests. Haplotype rs12740310*C-rs3737296*G-rs12410279*A was overtransmitted (P(corrected)= 0.0016), with a relative risk for autism of 1.8 in homozygous carriers. Furthermore, ASD-associated SNP rs12410279 modulates the level of transcription of MARK1. We found that MARK1 was overexpressed in the prefrontal cortex (BA46) but not in cerebellar granule cells, on postmortem brain tissues from patients. MARK1 displayed an accelerated evolution along the lineage leading to humans, suggesting possible involvement of this gene in cognition. MARK1 encodes a kinase-regulating microtubule-dependent transport in axons and dendrites. Both overexpression and silencing of MARK1 resulted in significantly shorter dendrite length in mouse neocortical neurons and modified dendritic transport speed. As expected for a gene encoding a key polarity determinant Par-1 protein kinase, MARK1 is involved in axon-dendrite specification. Thus, MARK1 overexpression in humans may be responsible for subtle changes in dendritic functioning.

HENA, heterogeneous network-based data set for Alzheimer's disease.

Alzheimer's disease and other types of dementia are the top cause for disabilities in later life and various types of experiments have been performed to understand the underlying mechanisms of the disease with the aim of coming up with potential drug targets. These experiments have been carried out by scientists working in different domains such as proteomics, molecular biology, clinical diagnostics and genomics. The results of such experiments are stored in the databases designed for collecting data of similar types. However, in order to get a systematic view of the disease from these independent but complementary data sets, it is necessary to combine them. In this study we describe a heterogeneous network-based data set for Alzheimer's disease (HENA). Additionally, we demonstrate the application of state-of-the-art graph convolutional networks, i.e. deep learning methods for the analysis of such large heterogeneous biological data sets. We expect HENA to allow scientists to explore and analyze their own results in the broader context of Alzheimer's disease research.

Somatostatin-IRES-Cre Mice: Between Knockout and Wild-Type?

The neuropeptide somatostatin (SOM) is widely expressed in rodent brain and somatostatin-IRES-Cre (SOM-cre) mouse strains are increasingly used to unravel the physiology of SOM-containing neurons. However, while knock-in targeting strategy greatly improves Cre-Lox system accuracy, recent reports have shown that genomic insertion of Cre construct per se can markedly affect physiological function. We show that Cre transgene insertion into the 3'UTR of the somatostatin gene leads to the selective and massive depletion of endogenous SOM in all tested brain regions. It also strongly impacts SOM-related neuroendocrine responses in a similar manner to what has been reported for SST KO mice: increased corticosterone levels after 30-min restraint stress, decreased amplitude and regularity of ultradian growth hormone secretory patterns accompanied by changes in sexually dimorphic liver gene expression (serpina1, Cyp2b9, Cyp2a4, Cyp2d9, and Cyp7b1). In addition to demonstrating the need for examination of the consequences of Cre transgenesis, these results also reveal how this SOM-cre strain may be a useful tool in studying the functional consequences of moderate to low SOM levels as reported in neurological and psychiatric disorders.

Fluorescent nanodiamond tracking reveals intraneuronal transport abnormalities induced by brain-disease-related genetic risk factors.

Brain diseases such as autism and Alzheimer's disease (each inflicting >1% of the world population) involve a large network of genes displaying subtle changes in their expression. Abnormalities in intraneuronal transport have been linked to genetic risk factors found in patients, suggesting the relevance of measuring this key biological process. However, current techniques are not sensitive enough to detect minor abnormalities. Here we report a sensitive method to measure the changes in intraneuronal transport induced by brain-disease-related genetic risk factors using fluorescent nanodiamonds (FNDs). We show that the high brightness, photostability and absence of cytotoxicity allow FNDs to be tracked inside the branches of dissociated neurons with a spatial resolution of 12 nm and a temporal resolution of 50 ms. As proof of principle, we applied the FND tracking assay on two transgenic mouse lines that mimic the slight changes in protein concentration (∼30%) found in the brains of patients. In both cases, we show that the FND assay is sufficiently sensitive to detect these changes.