Spatial Transcriptomics and In Situ Sequencing to Study Alzheimer's Disease.

Although complex inflammatory-like alterations are observed around the amyloid plaques of Alzheimer's disease (AD), little is known about the molecular changes and cellular interactions that characterize this response. We investigate here, in an AD mouse model, the transcriptional changes occurring in tissue domains in a 100-µm diameter around amyloid plaques using spatial transcriptomics. We demonstrate early alterations in a gene co-expression network enriched for myelin and oligodendrocyte genes (OLIGs), whereas a multicellular gene co-expression network of plaque-induced genes (PIGs) involving the complement system, oxidative stress, lysosomes, and inflammation is prominent in the later phase of the disease. We confirm the majority of the observed alterations at the cellular level using in situ sequencing on mouse and human brain sections. Genome-wide spatial transcriptomics analysis provides an unprecedented approach to untangle the dysregulated cellular network in the vicinity of pathogenic hallmarks of AD and other brain diseases.

A Single-Cell Transcriptome Atlas of the Aging Drosophila Brain.

The diversity of cell types and regulatory states in the brain, and how these change during aging, remains largely unknown. We present a single-cell transcriptome atlas of the entire adult Drosophila melanogaster brain sampled across its lifespan. Cell clustering identified 87 initial cell clusters that are further subclustered and validated by targeted cell-sorting. Our data show high granularity and identify a wide range of cell types. Gene network analyses using SCENIC revealed regulatory heterogeneity linked to energy consumption. During aging, RNA content declines exponentially without affecting neuronal identity in old brains. This single-cell brain atlas covers nearly all cells in the normal brain and provides the tools to study cellular diversity alongside other Drosophila and mammalian single-cell datasets in our unique single-cell analysis platform: SCope (http://scope.aertslab.org). These results, together with SCope, allow comprehensive exploration of all transcriptional states of an entire aging brain.

The γ-secretase substrate proteome and its role in cell signaling regulation.

γ-Secretases mediate the regulated intramembrane proteolysis (RIP) of more than 150 integral membrane proteins. We developed an unbiased γ-secretase substrate identification (G-SECSI) method to study to what extent these proteins are processed in parallel. We demonstrate here parallel processing of at least 85 membrane proteins in human microglia in steady-state cell culture conditions. Pharmacological inhibition of γ-secretase caused substantial changes of human microglial transcriptomes, including the expression of genes related to the disease-associated microglia (DAM) response described in Alzheimer disease (AD). While the overall effects of γ-secretase deficiency on transcriptomic cell states remained limited in control conditions, exposure of mouse microglia to AD-inducing amyloid plaques strongly blocked their capacity to mount this putatively protective DAM cell state. We conclude that γ-secretase serves as a critical signaling hub integrating the effects of multiple extracellular stimuli into the overall transcriptome of the cell.

Melanoma addiction to the long non-coding RNA SAMMSON.

Focal amplifications of chromosome 3p13-3p14 occur in about 10% of melanomas and are associated with a poor prognosis. The melanoma-specific oncogene MITF resides at the epicentre of this amplicon. However, whether other loci present in this amplicon also contribute to melanomagenesis is unknown. Here we show that the recently annotated long non-coding RNA (lncRNA) gene SAMMSON is consistently co-gained with MITF. In addition, SAMMSON is a target of the lineage-specific transcription factor SOX10 and its expression is detectable in more than 90% of human melanomas. Whereas exogenous SAMMSON increases the clonogenic potential in trans, SAMMSON knockdown drastically decreases the viability of melanoma cells irrespective of their transcriptional cell state and BRAF, NRAS or TP53 mutational status. Moreover, SAMMSON targeting sensitizes melanoma to MAPK-targeting therapeutics both in vitro and in patient-derived xenograft models. Mechanistically, SAMMSON interacts with p32, a master regulator of mitochondrial homeostasis and metabolism, to increase its mitochondrial targeting and pro-oncogenic function. Our results indicate that silencing of the lineage addiction oncogene SAMMSON disrupts vital mitochondrial functions in a cancer-cell-specific manner; this silencing is therefore expected to deliver highly effective and tissue-restricted anti-melanoma therapeutic responses.

Multiplex enhancer-reporter assays uncover unsophisticated TP53 enhancer logic.

Transcription factors regulate their target genes by binding to regulatory regions in the genome. Although the binding preferences of TP53 are known, it remains unclear what distinguishes functional enhancers from nonfunctional binding. In addition, the genome is scattered with recognition sequences that remain unoccupied. Using two complementary techniques of multiplex enhancer-reporter assays, we discovered that functional enhancers could be discriminated from nonfunctional binding events by the occurrence of a single TP53 canonical motif. By combining machine learning with a meta-analysis of TP53 ChIP-seq data sets, we identified a core set of more than 1000 responsive enhancers in the human genome. This TP53 cistrome is invariably used between cell types and experimental conditions, whereas differences among experiments can be attributed to indirect nonfunctional binding events. Our data suggest that TP53 enhancers represent a class of unsophisticated cell-autonomous enhancers containing a single TP53 binding site, distinct from complex developmental enhancers that integrate signals from multiple transcription factors.

Human stem cell-derived monocytes and microglia-like cells reveal impaired amyloid plaque clearance upon heterozygous or homozygous loss of TREM2.

INTRODUCTION: Murine microglia expressing the Alzheimer's disease-linked TREM2R47H mutation display variable decrease in phagocytosis, while impaired phagocytosis is reported following loss of TREM2. However, no data exist on TREM2+/R47H human microglia. Therefore, we created human pluripotent stem cell (hPSC) monocytes and transdifferentiated microglia-like cells (tMGs) to examine the effect of the TREM2+/R47H mutation and loss of TREM2 on phagocytosis. METHODS: We generated isogenic TREM2+/R47H, TREM2+/-, and TREM2-/- hPSCs using CRISPR/Cas9. Following differentiation to monocytes and tMGs, we studied the uptake of Escherichia coli fragments and analyzed amyloid plaque clearance from cryosections of APP/PS1+/- mouse brains. RESULTS: We demonstrated that tMGs resemble cultured human microglia. TREM2+/- and TREM2-/- hPSC monocytes and tMGs phagocytosed significantly less E. coli fragments and cleared less amyloid plaques than wild-type hPSC progeny, with no difference for TREM2+/R47H progeny. DISCUSSION: In vitro phagocytosis of hPSC monocytes and tMGs was not affected by the TREM2+/R47H mutation but was significantly impaired in TREM2+/- and TREM2-/- progeny.

An immunoaffinity-based method for isolating ultrapure adult astrocytes based on ATP1B2 targeting by the ACSA-2 antibody.

Astrocytes are a major cell type in the mammalian CNS. Astrocytes are now known to play a number of essential roles in processes including synapse formation and function, as well as blood-brain barrier formation and control of cerebral blood flow. However, our understanding of the molecular mechanisms underlying astrocyte development and function is still rudimentary. This lack of knowledge is at least partly due to the lack of tools currently available for astrocyte biology. ACSA-2 is a commercially available antibody originally developed for the isolation of astrocytes from young postnatal mouse brain, using magnetic cell-sorting methods, but its utility in isolating cells from adult tissue has not yet been published. Using a modified protocol, we now show that this tool can also be used to isolate ultrapure astrocytes from the adult brain. Furthermore, using a variety of techniques (including single-cell sequencing, overexpression and knockdown assays, immunoblotting, and immunohistochemistry), we identify the ACSA-2 epitope for the first time as ATP1B2 and characterize its distribution in the CNS. Finally, we show that ATP1B2 is stably expressed in multiple models of CNS injury and disease. Hence, we show that the ACSA-2 antibody possesses the potential to be an extremely valuable tool for astrocyte research, allowing the purification and characterization of astrocytes (potentially including injury and disease models) without the need for any specialized and expensive equipment. In fact, our results suggest that ACSA-2 should be a first-choice method for astrocyte isolation and characterization.

FXR2P Exerts a Positive Translational Control and Is Required for the Activity-Dependent Increase of PSD95 Expression.

In brain, specific RNA-binding proteins (RBPs) associate with localized mRNAs and function as regulators of protein synthesis at synapses exerting an indirect control on neuronal activity. Thus, the Fragile X Mental Retardation protein (FMRP) regulates expression of the scaffolding postsynaptic density protein PSD95, but the mode of control appears to be different from other FMRP target mRNAs. Here, we show that the fragile X mental retardation-related protein 2 (FXR2P) cooperates with FMRP in binding to the 3'-UTR of mouse PSD95/Dlg4 mRNA. Absence of FXR2P leads to decreased translation of PSD95/Dlg4 mRNA in the hippocampus, implying a role for FXR2P as translation activator. Remarkably, mGluR-dependent increase of PSD95 synthesis is abolished in neurons lacking Fxr2. Together, these findings show a coordinated regulation of PSD95/Dlg4 mRNA by FMRP and FXR2P that ultimately affects its fine-tuning during synaptic activity.

JAK3 mutants transform hematopoietic cells through JAK1 activation, causing T-cell acute lymphoblastic leukemia in a mouse model.

JAK3 is a tyrosine kinase that associates with the common γ chain of cytokine receptors and is recurrently mutated in T-cell acute lymphoblastic leukemia (T-ALL). We tested the transforming properties of JAK3 pseudokinase and kinase domain mutants using in vitro and in vivo assays. Most, but not all, JAK3 mutants transformed cytokine-dependent Ba/F3 or MOHITO cell lines to cytokine-independent proliferation. JAK3 pseudokinase mutants were dependent on Jak1 kinase activity for cellular transformation, whereas the JAK3 kinase domain mutant could transform cells in a Jak1 kinase-independent manner. Reconstitution of the IL7 receptor signaling complex in 293T cells showed that JAK3 mutants required receptor binding to mediate downstream STAT5 phosphorylation. Mice transplanted with bone marrow progenitor cells expressing JAK3 mutants developed a long-latency transplantable T-ALL-like disease, characterized by an accumulation of immature CD8(+) T cells. In vivo treatment of leukemic mice with the JAK3 selective inhibitor tofacitinib reduced the white blood cell count and caused leukemic cell apoptosis. Our data show that JAK3 mutations are drivers of T-ALL and require the cytokine receptor complex for transformation. These results warrant further investigation of JAK1/JAK3 inhibitors for the treatment of T-ALL.

Identification of High-Impact cis-Regulatory Mutations Using Transcription Factor Specific Random Forest Models.

Cancer genomes contain vast amounts of somatic mutations, many of which are passenger mutations not involved in oncogenesis. Whereas driver mutations in protein-coding genes can be distinguished from passenger mutations based on their recurrence, non-coding mutations are usually not recurrent at the same position. Therefore, it is still unclear how to identify cis-regulatory driver mutations, particularly when chromatin data from the same patient is not available, thus relying only on sequence and expression information. Here we use machine-learning methods to predict functional regulatory regions using sequence information alone, and compare the predicted activity of the mutated region with the reference sequence. This way we define the Predicted Regulatory Impact of a Mutation in an Enhancer (PRIME). We find that the recently identified driver mutation in the TAL1 enhancer has a high PRIME score, representing a "gain-of-target" for MYB, whereas the highly recurrent TERT promoter mutation has a surprisingly low PRIME score. We trained Random Forest models for 45 cancer-related transcription factors, and used these to score variations in the HeLa genome and somatic mutations across more than five hundred cancer genomes. Each model predicts only a small fraction of non-coding mutations with a potential impact on the function of the encompassing regulatory region. Nevertheless, as these few candidate driver mutations are often linked to gains in chromatin activity and gene expression, they may contribute to the oncogenic program by altering the expression levels of specific oncogenes and tumor suppressor genes.

Genomic analysis of the Kiwifruit pathogen Pseudomonas syringae pv. actinidiae provides insight into the origins of an emergent plant disease.

The origins of crop diseases are linked to domestication of plants. Most crops were domesticated centuries--even millennia--ago, thus limiting opportunity to understand the concomitant emergence of disease. Kiwifruit (Actinidia spp.) is an exception: domestication began in the 1930s with outbreaks of canker disease caused by P. syringae pv. actinidiae (Psa) first recorded in the 1980s. Based on SNP analyses of two circularized and 34 draft genomes, we show that Psa is comprised of distinct clades exhibiting negligible within-clade diversity, consistent with disease arising by independent samplings from a source population. Three clades correspond to their geographical source of isolation; a fourth, encompassing the Psa-V lineage responsible for the 2008 outbreak, is now globally distributed. Psa has an overall clonal population structure, however, genomes carry a marked signature of within-pathovar recombination. SNP analysis of Psa-V reveals hundreds of polymorphisms; however, most reside within PPHGI-1-like conjugative elements whose evolution is unlinked to the core genome. Removal of SNPs due to recombination yields an uninformative (star-like) phylogeny consistent with diversification of Psa-V from a single clone within the last ten years. Growth assays provide evidence of cultivar specificity, with rapid systemic movement of Psa-V in Actinidia chinensis. Genomic comparisons show a dynamic genome with evidence of positive selection on type III effectors and other candidate virulence genes. Each clade has highly varied complements of accessory genes encoding effectors and toxins with evidence of gain and loss via multiple genetic routes. Genes with orthologs in vascular pathogens were found exclusively within Psa-V. Our analyses capture a pathogen in the early stages of emergence from a predicted source population associated with wild Actinidia species. In addition to candidate genes as targets for resistance breeding programs, our findings highlight the importance of the source population as a reservoir of new disease.

iRegulon: from a gene list to a gene regulatory network using large motif and track collections.

Identifying master regulators of biological processes and mapping their downstream gene networks are key challenges in systems biology. We developed a computational method, called iRegulon, to reverse-engineer the transcriptional regulatory network underlying a co-expressed gene set using cis-regulatory sequence analysis. iRegulon implements a genome-wide ranking-and-recovery approach to detect enriched transcription factor motifs and their optimal sets of direct targets. We increase the accuracy of network inference by using very large motif collections of up to ten thousand position weight matrices collected from various species, and linking these to candidate human TFs via a motif2TF procedure. We validate iRegulon on gene sets derived from ENCODE ChIP-seq data with increasing levels of noise, and we compare iRegulon with existing motif discovery methods. Next, we use iRegulon on more challenging types of gene lists, including microRNA target sets, protein-protein interaction networks, and genetic perturbation data. In particular, we over-activate p53 in breast cancer cells, followed by RNA-seq and ChIP-seq, and could identify an extensive up-regulated network controlled directly by p53. Similarly we map a repressive network with no indication of direct p53 regulation but rather an indirect effect via E2F and NFY. Finally, we generalize our computational framework to include regulatory tracks such as ChIP-seq data and show how motif and track discovery can be combined to map functional regulatory interactions among co-expressed genes. iRegulon is available as a Cytoscape plugin from http://iregulon.aertslab.org.

Genomes of 'Candidatus Liberibacter solanacearum' Haplotype A from New Zealand and the United States Suggest Significant Genome Plasticity in the Species.

'Candidatus Liberibacter solanacearum' contains two solanaceous crop-infecting haplotypes, A and B. Two haplotype A draft genomes were assembled and compared with ZC1 (haplotype B), revealing inversion and relocation genomic rearrangements, numerous single-nucleotide polymorphisms, and differences in phage-related regions. Differences in prophage location and sequence were seen both within and between haplotype comparisons. OrthoMCL and BLAST analyses identified 46 putative coding sequences present in haplotype A that were not present in haplotype B. Thirty-eight of these loci were not found in sequences from other Liberibacter spp. Quantitative polymerase chain reaction (qPCR) assays designed to amplify sequences from 15 of these loci were screened against a panel of 'Ca. L. solanacearum'-positive samples to investigate genetic diversity. Seven of the assays demonstrated within-haplotype diversity; five failed to amplify loci in at least one haplotype A sample while three assays produced amplicons from some haplotype B samples. Eight of the loci assays showed consistent A-B differentiation. Differences in genome arrangements, prophage, and qPCR results suggesting locus diversity within the haplotypes provide more evidence for genetic complexity in this emerging bacterial species.

Structure and expression of GSL1 and GSL2 genes encoding gibberellin stimulated-like proteins in diploid and highly heterozygous tetraploid potato reveals their highly conserved and essential status.

BACKGROUND: GSL1 and GSL2, Gibberellin Stimulated-Like proteins (also known as Snakin-1 and Snakin-2), are cysteine-rich peptides from potato (Solanum tuberosum L.) with antimicrobial properties. Similar peptides in other species have been implicated in diverse biological processes and are hypothesised to play a role in several aspects of plant development, plant responses to biotic or abiotic stress through their participation in hormone crosstalk, and redox homeostasis. To help resolve the biological roles of GSL1 and GSL2 peptides we have undertaken an in depth analysis of the structure and expression of these genes in potato. RESULTS: We have characterised the full length genes for both GSL1 (chromosome 4) and GSL2 (chromosome 1) from diploid and tetraploid potato using the reference genome sequence of potato, coupled with further next generation sequencing of four highly heterozygous tetraploid cultivars. The frequency of SNPs in GSL1 and GSL2 were very low with only one SNP every 67 and 53 nucleotides in exon regions of GSL1 and GSL2, respectively. Analysis of comprehensive RNA-seq data substantiated the role of specific promoter motifs in transcriptional control of gene expression. Expression analysis based on the frequency of next generation sequence reads established that GSL2 was expressed at a higher level than GSL1 in 30 out of 32 tissue and treatment libraries. Furthermore, both the GSL1 and GSL2 genes exhibited constitutive expression that was not up regulated in response to biotic or abiotic stresses, hormone treatments or wounding. Potato transformation with antisense knock-down expression cassettes failed to recover viable plants. CONCLUSIONS: The potato GSL1 and GSL2 genes are very highly conserved suggesting they contribute to an important biological function. The known antimicrobial activity of the GSL proteins, coupled with the FPKM analysis from RNA-seq data, implies that both genes contribute to the constitutive defence barriers in potatoes. The lethality of antisense knock-down expression of GSL1 and GSL2, coupled with the rare incidence of SNPs in these genes, suggests an essential role for this gene family. These features are consistent with the GSL protein family playing a role in several aspects of plant development in addition to plant defence against biotic stresses.

Regulation of human microglial gene expression and function via RNAase-H active antisense oligonucleotides in vivo in Alzheimer's disease.

BACKGROUND: Microglia play important roles in maintaining brain homeostasis and neurodegeneration. The discovery of genetic variants in genes predominately or exclusively expressed in myeloid cells, such as Apolipoprotein E (APOE) and triggering receptor expressed on myeloid cells 2 (TREM2), as the strongest risk factors for Alzheimer's disease (AD) highlights the importance of microglial biology in the brain. The sequence, structure and function of several microglial proteins are poorly conserved across species, which has hampered the development of strategies aiming to modulate the expression of specific microglial genes. One way to target APOE and TREM2 is to modulate their expression using antisense oligonucleotides (ASOs). METHODS: In this study, we identified, produced, and tested novel, selective and potent ASOs for human APOE and TREM2. We used a combination of in vitro iPSC-microglia models, as well as microglial xenotransplanted mice to provide proof of activity in human microglial in vivo. RESULTS: We proved their efficacy in human iPSC microglia in vitro, as well as their pharmacological activity in vivo in a xenografted microglia model. We demonstrate ASOs targeting human microglia can modify their transcriptional profile and their response to amyloid-ß plaques in vivo in a model of AD. CONCLUSIONS: This study is the first proof-of-concept that human microglial can be modulated using ASOs in a dose-dependent manner to manipulate microglia phenotypes and response to neurodegeneration in vivo.

Xenografted human microglia display diverse transcriptomic states in response to Alzheimer's disease-related amyloid-ß pathology.

Microglia are central players in Alzheimer's disease pathology but analyzing microglial states in human brain samples is challenging due to genetic diversity, postmortem delay and admixture of pathologies. To circumvent these issues, here we generated 138,577 single-cell expression profiles of human stem cell-derived microglia xenotransplanted in the brain of the AppNL-G-F model of amyloid pathology and wild-type controls. Xenografted human microglia adopt a disease-associated profile similar to that seen in mouse microglia, but display a more pronounced human leukocyte antigen or HLA state, likely related to antigen presentation in response to amyloid plaques. The human microglial response also involves a pro-inflammatory cytokine/chemokine cytokine response microglia or CRM response to oligomeric Aß oligomers. Genetic deletion of TREM2 or APOE as well as APOE polymorphisms and TREM2R47H expression in the transplanted microglia modulate these responses differentially. The expression of other Alzheimer's disease risk genes is differentially regulated across the distinct cell states elicited in response to amyloid pathology. Thus, we have identified multiple transcriptomic cell states adopted by human microglia in a multipronged response to Alzheimer's disease-related pathology, which should be taken into account in translational studies.

Microglia-astrocyte crosstalk in the amyloid plaque niche of an Alzheimer's disease mouse model, as revealed by spatial transcriptomics.

The amyloid plaque niche is a pivotal hallmark of Alzheimer's disease (AD). Here, we employ two high-resolution spatial transcriptomics (ST) platforms, CosMx and Spatial Enhanced Resolution Omics-sequencing (Stereo-seq), to characterize the transcriptomic alterations, cellular compositions, and signaling perturbations in the amyloid plaque niche in an AD mouse model. We discover heterogeneity in the cellular composition of plaque niches, marked by an increase in microglial accumulation. We profile the transcriptomic alterations of glial cells in the vicinity of plaques and conclude that the microglial response to plaques is consistent across different brain regions, while the astrocytic response is more heterogeneous. Meanwhile, as the microglial density of plaque niches increases, astrocytes acquire a more neurotoxic phenotype and play a key role in inducing GABAergic signaling and decreasing glutamatergic signaling in hippocampal neurons. We thus show that the accumulation of microglia around hippocampal plaques disrupts astrocytic signaling, in turn inducing an imbalance in neuronal synaptic signaling.

Arabidopsis thaliana transcription factors bZIP19 and bZIP23 regulate the adaptation to zinc deficiency.

Zinc is an essential micronutrient for all living organisms. When facing a shortage in zinc supply, plants adapt by enhancing the zinc uptake capacity. The molecular regulators controlling this adaptation are not known. We present the identification of two closely related members of the Arabidopsis thaliana basic-region leucine-zipper (bZIP) transcription factor gene family, bZIP19 and bZIP23, that regulate the adaptation to low zinc supply. They were identified, in a yeast-one-hybrid screening, to associate to promoter regions of the zinc deficiency-induced ZIP4 gene of the Zrt- and Irt-related protein (ZIP) family of metal transporters. Although mutation of only one of the bZIP genes hardly affects plants, we show that the bzip19 bzip23 double mutant is hypersensitive to zinc deficiency. Unlike the wild type, the bzip19 bzip23 mutant is unable to induce the expression of a small set of genes that constitutes the primary response to zinc deficiency, comprising additional ZIP metal transporter genes. This set of target genes is characterized by the presence of one or more copies of a 10-bp imperfect palindrome in their promoter region, to which both bZIP proteins can bind. The bZIP19 and bZIP23 transcription factors, their target genes, and the characteristic cis zinc deficiency response elements they can bind to are conserved in higher plants. These findings are a significant step forward to unravel the molecular mechanism of zinc homeostasis in plants, allowing the improvement of zinc bio-fortification to alleviate human nutrition problems and phytoremediation strategies to clean contaminated soils.

MEG3 activates necroptosis in human neuron xenografts modeling Alzheimer's disease.

Neuronal cell loss is a defining feature of Alzheimer's disease (AD), but the underlying mechanisms remain unclear. We xenografted human or mouse neurons into the brain of a mouse model of AD. Only human neurons displayed tangles, Gallyas silver staining, granulovacuolar neurodegeneration (GVD), phosphorylated tau blood biomarkers, and considerable neuronal cell loss. The long noncoding RNA MEG3 was strongly up-regulated in human neurons. This neuron-specific long noncoding RNA is also up-regulated in AD patients. MEG3 expression alone was sufficient to induce necroptosis in human neurons in vitro. Down-regulation of MEG3 and inhibition of necroptosis using pharmacological or genetic manipulation of receptor-interacting protein kinase 1 (RIPK1), RIPK3, or mixed lineage kinase domain-like protein (MLKL) rescued neuronal cell loss in xenografted human neurons. This model suggests potential therapeutic approaches for AD and reveals a human-specific vulnerability to AD.