Sex and cell-specific gene expression in corticolimbic brain regions associated with psychiatric disorders revealed by bulk and single-nuclei RNA sequencing.

BACKGROUND: There are sex-specific differences in the prevalence, symptomology and course of psychiatric disorders. However, preclinical models have primarily used males, such that the molecular mechanisms underlying sex-specific differences in psychiatric disorders are not well established. METHODS: In this study, we compared transcriptome-wide gene expression profiles in male and female rats within the corticolimbic system, including the cingulate cortex, nucleus accumbens medial shell (NAcS), ventral dentate gyrus and the basolateral amygdala (n = 22-24 per group/region). FINDINGS: We found over 3000 differentially expressed genes (DEGs) in the NAcS between males and females. Of these DEGs in the NAcS, 303 showed sex-dependent conservation DEGs in humans and were significantly enriched for gene ontology terms related to blood vessel morphogenesis and regulation of cell migration. Single nuclei RNA sequencing in the NAcS of male and female rats identified widespread sex-dependent expression, with genes upregulated in females showing a notable enrichment for synaptic function. Female upregulated genes in astrocytes, Drd3+MSNs and oligodendrocyte were also enriched in several psychiatric genome-wide association studies (GWAS). INTERPRETATION: Our data provide comprehensive evidence of sex- and cell-specific molecular profiles in the NAcS. Importantly these differences associate with anxiety, bipolar disorder, schizophrenia, and cross-disorder, suggesting an intrinsic molecular basis for sex-based differences in psychiatric disorders that strongly implicates the NAcS. FUNDING: This work was supported by funding from the Hope for Depression Research Foundation (MJM).

A single-cell atlas identifies pretreatment features of primary imatinib resistance in chronic myeloid leukemia.

Primary resistance to tyrosine kinase inhibitors (TKIs) is a significant barrier to optimal outcomes in chronic myeloid leukemia (CML), but factors contributing to response heterogeneity remain unclear. Using single-cell RNA (scRNA) sequencing, we identified 8 statistically significant features in pretreatment bone marrow, which correlated with either sensitivity (major molecular response or MMR) or extreme resistance to imatinib (eventual blast crisis [BC] transformation). Employing machine-learning, we identified leukemic stem cell (LSC) and natural killer (NK) cell gene expression profiles predicting imatinib response with >80% accuracy, including no false positives for predicting BC. A canonical erythroid-specifying (TAL1/KLF1/GATA1) regulon was a hallmark of LSCs from patients with MMR and was associated with erythroid progenitor [ERP] expansion in vivo (P < .05), and a 2- to 10-fold (6.3-fold in group A vs 1.09-fold in group C) erythroid over myeloid bias in vitro. Notably, ERPs demonstrated exquisite TKI sensitivity compared with myeloid progenitors (P < .001). These LSC features were lost with progressive resistance, and MYC- and IRF1-driven inflammatory regulons were evident in patients who progressed to transformation. Patients with MMR also exhibited a 56-fold expansion (P < .01) of a normally rare subset of hyperfunctional adaptive-like NK cells, which diminished with progressive resistance, whereas patients destined for BC accumulated inhibitory NKG2A+ NK cells favoring NK cell tolerance. Finally, we developed antibody panels to validate our scRNA-seq findings. These panels may be useful for prospective studies of primary resistance, and in assessing the contribution of predetermined vs acquired factors in TKI response heterogeneity.

Integrative multi-omics landscape of fluoxetine action across 27 brain regions reveals global increase in energy metabolism and region-specific chromatin remodelling.

Depression and anxiety are major global health burdens. Although SSRIs targeting the serotonergic system are prescribed over 200 million times annually, they have variable therapeutic efficacy and side effects, and mechanisms of action remain incompletely understood. Here, we comprehensively characterise the molecular landscape of gene regulatory changes associated with fluoxetine, a widely-used SSRI. We performed multimodal analysis of SSRI response in 27 mammalian brain regions using 310 bulk RNA-seq and H3K27ac ChIP-seq datasets, followed by in-depth characterisation of two hippocampal regions using single-cell RNA-seq (20 datasets). Remarkably, fluoxetine induced profound region-specific shifts in gene expression and chromatin state, including in the nucleus accumbens shell, locus coeruleus and septal areas, as well as in more well-studied regions such as the raphe and hippocampal dentate gyrus. Expression changes were strongly enriched at GWAS loci for depression and antidepressant drug response, stressing the relevance to human phenotypes. We observed differential expression at dozens of signalling receptors and pathways, many of which are previously unknown. Single-cell analysis revealed stark differences in fluoxetine response between the dorsal and ventral hippocampal dentate gyri, particularly in oligodendrocytes, mossy cells and inhibitory neurons. Across diverse brain regions, integrative omics analysis consistently suggested increased energy metabolism via oxidative phosphorylation and mitochondrial changes, which we corroborated in vitro; this may thus constitute a shared mechanism of action of fluoxetine. Similarly, we observed pervasive chromatin remodelling signatures across the brain. Our study reveals unexpected regional and cell type-specific heterogeneity in SSRI action, highlights under-studied brain regions that may play a major role in antidepressant response, and provides a rich resource of candidate cell types, genes, gene regulatory elements and pathways for mechanistic analysis and identifying new therapeutic targets for depression and anxiety.

Histone acetylome-wide associations in immune cells from individuals with active Mycobacterium tuberculosis infection.

Host cell chromatin changes are thought to play an important role in the pathogenesis of infectious diseases. Here we describe a histone acetylome-wide association study (HAWAS) of an infectious disease, on the basis of genome-wide H3K27 acetylation profiling of peripheral blood granulocytes and monocytes from persons with active Mycobacterium tuberculosis (Mtb) infection and healthy controls. We detected >2,000 differentially acetylated loci in either cell type in a Singapore Chinese discovery cohort (n = 46), which were validated in a subsequent multi-ethnic Singapore cohort (n = 29), as well as a longitudinal cohort from South Africa (n = 26), thus demonstrating that HAWAS can be independently corroborated. Acetylation changes were correlated with differential gene expression. Differential acetylation was enriched near potassium channel genes, including KCNJ15, which modulates apoptosis and promotes Mtb clearance in vitro. We performed histone acetylation quantitative trait locus (haQTL) analysis on the dataset and identified 69 candidate causal variants for immune phenotypes among granulocyte haQTLs and 83 among monocyte haQTLs. Our study provides proof-of-principle for HAWAS to infer mechanisms of host response to pathogens.

DUBStepR is a scalable correlation-based feature selection method for accurately clustering single-cell data.

Feature selection (marker gene selection) is widely believed to improve clustering accuracy, and is thus a key component of single cell clustering pipelines. Existing feature selection methods perform inconsistently across datasets, occasionally even resulting in poorer clustering accuracy than without feature selection. Moreover, existing methods ignore information contained in gene-gene correlations. Here, we introduce DUBStepR (Determining the Underlying Basis using Stepwise Regression), a feature selection algorithm that leverages gene-gene correlations with a novel measure of inhomogeneity in feature space, termed the Density Index (DI). Despite selecting a relatively small number of genes, DUBStepR substantially outperformed existing single-cell feature selection methods across diverse clustering benchmarks. Additionally, DUBStepR was the only method to robustly deconvolve T and NK heterogeneity by identifying disease-associated common and rare cell types and subtypes in PBMCs from rheumatoid arthritis patients. DUBStepR is scalable to over a million cells, and can be straightforwardly applied to other data types such as single-cell ATAC-seq. We propose DUBStepR as a general-purpose feature selection solution for accurately clustering single-cell data.

RCA2: a scalable supervised clustering algorithm that reduces batch effects in scRNA-seq data.

The transcriptomic diversity of cell types in the human body can be analysed in unprecedented detail using single cell (SC) technologies. Unsupervised clustering of SC transcriptomes, which is the default technique for defining cell types, is prone to group cells by technical, rather than biological, variation. Compared to de-novo (unsupervised) clustering, we demonstrate using multiple benchmarks that supervised clustering, which uses reference transcriptomes as a guide, is robust to batch effects and data quality artifacts. Here, we present RCA2, the first algorithm to combine reference projection (batch effect robustness) with graph-based clustering (scalability). In addition, RCA2 provides a user-friendly framework incorporating multiple commonly used downstream analysis modules. RCA2 also provides new reference panels for human and mouse and supports generation of custom panels. Furthermore, RCA2 facilitates cell type-specific QC, which is essential for accurate clustering of data from heterogeneous tissues. We demonstrate the advantages of RCA2 on SC data from human bone marrow, healthy PBMCs and PBMCs from COVID-19 patients. Scalable supervised clustering methods such as RCA2 will facilitate unified analysis of cohort-scale SC datasets.

Early Fate Defines Microglia and Non-parenchymal Brain Macrophage Development.

Central nervous system (CNS) macrophages comprise microglia and border-associated macrophages (BAMs) residing in the meninges, the choroid plexus, and the perivascular spaces. Most CNS macrophages emerge during development, with the exception of choroid plexus and dural macrophages, which are replaced by monocytes in adulthood. Whether microglia and BAMs share a developmental program or arise from separate lineages remains unknown. Here, we identified two phenotypically, transcriptionally, and locally distinct brain macrophages throughout development, giving rise to either microglia or BAMs. Two macrophage populations were already present in the yolk sac suggesting an early segregation. Fate-mapping models revealed that BAMs mostly derived from early erythro-myeloid progenitors in the yolk sac. The development of microglia was dependent on TGF-ß, whereas the genesis of BAMs occurred independently of this cytokine. Collectively, our data show that developing parenchymal and non-parenchymal brain macrophages are separate entities in terms of ontogeny, gene signature, and requirement for TGF-ß.

Intrinsic mutant HTT-mediated defects in oligodendroglia cause myelination deficits and behavioral abnormalities in Huntington disease.

White matter abnormalities are a nearly universal pathological feature of neurodegenerative disorders including Huntington disease (HD). A long-held assumption is that this white matter pathology is simply a secondary outcome of the progressive neuronal loss that manifests with advancing disease. Using a mouse model of HD, here we show that white matter and myelination abnormalities are an early disease feature appearing before the manifestation of any behavioral abnormalities or neuronal loss. We further show that selective inactivation of mutant huntingtin (mHTT) in the NG2+ oligodendrocyte progenitor cell population prevented myelin abnormalities and certain behavioral deficits in HD mice. Strikingly, the improvements in behavioral outcomes were seen despite the continued expression of mHTT in nonoligodendroglial cells including neurons, astrocytes, and microglia. Using RNA-seq and ChIP-seq analyses, we implicate a pathogenic mechanism that involves enhancement of polycomb repressive complex 2 (PRC2) activity by mHTT in the intrinsic oligodendroglial dysfunction and myelination deficits observed in HD. Our findings challenge the long-held dogma regarding the etiology of white matter pathology in HD and highlight the contribution of epigenetic mechanisms to the observed intrinsic oligodendroglial dysfunction. Our results further suggest that ameliorating white matter pathology and oligodendroglial dysfunction may be beneficial for HD.

Massive contribution of transposable elements to mammalian regulatory sequences.

Barbara McClintock discovered the existence of transposable elements (TEs) in the late 1940s and initially proposed that they contributed to the gene regulatory program of higher organisms. This controversial idea gained acceptance only much later in the 1990s, when the first examples of TE-derived promoter sequences were uncovered. It is now known that half of the human genome is recognizably derived from TEs. It is thus important to understand the scope and nature of their contribution to gene regulation. Here, we provide a timeline of major discoveries in this area and discuss how transposons have revolutionized our understanding of mammalian genomes, with a special emphasis on the massive contribution of TEs to primate evolution. Our analysis of primate-specific functional elements supports a simple model for the rate at which new functional elements arise in unique and TE-derived DNA. Finally, we discuss some of the challenges and unresolved questions in the field, which need to be addressed in order to fully characterize the impact of TEs on gene regulation, evolution and disease processes.

Comprehensive benchmarking reveals H2BK20 acetylation as a distinctive signature of cell-state-specific enhancers and promoters.

Although over 35 different histone acetylation marks have been described, the overwhelming majority of regulatory genomics studies focus exclusively on H3K27ac and H3K9ac. In order to identify novel epigenomic traits of regulatory elements, we constructed a benchmark set of validated enhancers by performing 140 enhancer assays in human T cells. We tested 40 chromatin signatures on this unbiased enhancer set and identified H2BK20ac, a little-studied histone modification, as the most predictive mark of active enhancers. Notably, we detected a novel class of functionally distinct enhancers enriched in H2BK20ac but lacking H3K27ac, which was present in all examined cell lines and also in embryonic forebrain tissue. H2BK20ac was also unique in highlighting cell-type-specific promoters. In contrast, other acetylation marks were present in all active promoters, regardless of cell-type specificity. In stimulated microglial cells, H2BK20ac was more correlated with cell-state-specific expression changes than H3K27ac, with TGF-beta signaling decoupling the two acetylation marks at a subset of regulatory elements. In summary, our study reveals a previously unknown connection between histone acetylation and cell-type-specific gene regulation and indicates that H2BK20ac profiling can be used to uncover new dimensions of gene regulation.

Noncoding origins of anthropoid traits and a new null model of transposon functionalization.

Little is known about novel genetic elements that drove the emergence of anthropoid primates. We exploited the sequencing of the marmoset genome to identify 23,849 anthropoid-specific constrained (ASC) regions and confirmed their robust functional signatures. Of the ASC base pairs, 99.7% were noncoding, suggesting that novel anthropoid functional elements were overwhelmingly cis-regulatory. ASCs were highly enriched in loci associated with fetal brain development, motor coordination, neurotransmission, and vision, thus providing a large set of candidate elements for exploring the molecular basis of hallmark primate traits. We validated ASC192 as a primate-specific enhancer in proliferative zones of the developing brain. Unexpectedly, transposable elements (TEs) contributed to >56% of ASCs, and almost all TE families showed functional potential similar to that of nonrepetitive DNA. Three L1PA repeat-derived ASCs displayed coherent eye-enhancer function, thus demonstrating that the "gene-battery" model of TE functionalization applies to enhancers in vivo. Our study provides fundamental insights into genome evolution and the origins of anthropoid phenotypes and supports an elegantly simple new null model of TE exaptation.

Uniform, optimal signal processing of mapped deep-sequencing data.

Despite their apparent diversity, many problems in the analysis of high-throughput sequencing data are merely special cases of two general problems, signal detection and signal estimation. Here we adapt formally optimal solutions from signal processing theory to analyze signals of DNA sequence reads mapped to a genome. We describe DFilter, a detection algorithm that identifies regulatory features in ChIP-seq, DNase-seq and FAIRE-seq data more accurately than assay-specific algorithms. We also describe EFilter, an estimation algorithm that accurately predicts mRNA levels from as few as 1-2 histone profiles (R ∼0.9). Notably, the presence of regulatory motifs in promoters correlates more with histone modifications than with mRNA levels, suggesting that histone profiles are more predictive of cis-regulatory mechanisms. We show by applying DFilter and EFilter to embryonic forebrain ChIP-seq data that regulatory protein identification and functional annotation are feasible despite tissue heterogeneity. The mathematical formalism underlying our tools facilitates integrative analysis of data from virtually any sequencing-based functional profile.

Costunolide inhibits proinflammatory cytokines and iNOS in activated murine BV2 microglia.

Costunolide, a sesquiterpene lactone present in Costus speciosus root exerts a variety of pharmacological activity but its effects on neuroinflammation have not been studied. Microglia, the resident phagocytic cells in the central nervous system respond to neuroinflammation and their overwhelming response in turn aggravate brain damage during infection, ischemia and neurodegenerative diseases. In this study, we report the effect of Costunolide on the production of proinflammatory mediators and mechanisms involved in BV2 microglial cells stimulated with LPS. Costunolide attenuated the expression of tumour necrosis factor-alpha, interleukin-1,6, inducible nitric oxide synthase, monocyte chemotactic protein 1 and cyclooxygenase 2 in activated microglia. This Costunolide-mediated inhibition was correspondent with the inhibition of NFkappaB activation. It has been further shown that Costunolide suppressed MAPK pathway activation by inducing MKP-1 production. Collectively our results suggest that Costunolide shows an ability to inhibit expression of multiple neuroinflammatory mediators and this is attributable to the compounds inhibition of NFkappaB and MAPK activation. This novel role of Costunolide upon investigation may aid in developing better therapeutic strategies for treatment of neuroinflammatory diseases.