Cellular Diversity in Human Subgenual Anterior Cingulate and Dorsolateral Prefrontal Cortex by Single-Nucleus RNA-Sequencing.

Regional cellular heterogeneity is a fundamental feature of the human neocortex; however, details of this heterogeneity are still undefined. We used single-nucleus RNA-sequencing to examine cell-specific transcriptional features in the dorsolateral PFC (DLPFC) and the subgenual anterior cingulate cortex (sgACC), regions implicated in major psychiatric disorders. Droplet-based nuclei-capture and library preparation were performed on replicate samples from 8 male donors without history of psychiatric or neurologic disorder. Unsupervised clustering identified major neural cell classes. Subsequent iterative clustering of neurons further revealed 20 excitatory and 22 inhibitory subclasses. Inhibitory cells were consistently more abundant in the sgACC and excitatory neuron subclusters exhibited considerable variability across brain regions. Excitatory cell subclasses also exhibited greater within-class transcriptional differences between the two regions. We used these molecular definitions to determine which cell classes might be enriched in loci carrying a genetic signal in genome-wide association studies or for differentially expressed genes in mental illness. We found that the heritable signals of psychiatric disorders were enriched in neurons and that, while the gene expression changes detected in bulk-RNA-sequencing studies were dominated by glial cells, some alterations could be identified in specific classes of excitatory and inhibitory neurons. Intriguingly, only two excitatory cell classes exhibited concomitant region-specific enrichment for both genome-wide association study loci and transcriptional dysregulation. In sum, by detailing the molecular and cellular diversity of the DLPFC and sgACC, we were able to generate hypotheses on regional and cell-specific dysfunctions that may contribute to the development of mental illness.SIGNIFICANCE STATEMENT Dysfunction of the subgenual anterior cingulate cortex has been implicated in mood disorders, particularly major depressive disorder, and the dorsolateral PFC, a subsection of the PFC involved in executive functioning, has been implicated in schizophrenia. Understanding the cellular composition of these regions is critical to elucidating the neurobiology underlying psychiatric and neurologic disorders. We studied cell type diversity of the subgenual anterior cingulate cortex and dorsolateral PFC of humans with no neuropsychiatric illness using a clustering analysis of single-nuclei RNA-sequencing data. Defining the transcriptomic profile of cellular subpopulations in these cortical regions is a first step to demystifying the cellular and molecular pathways involved in psychiatric disorders.

Antipsychotic drug use complicates assessment of gene expression changes associated with schizophrenia.

Recent postmortem transcriptomic studies of schizophrenia (SCZ) have shown hundreds of differentially expressed genes. However, the extent to which these gene expression changes reflect antipsychotic drug (APD) exposure remains uncertain. We compared differential gene expression in the prefrontal cortex of SCZ patients who tested positive for APDs at the time of death with SCZ patients who did not. APD exposure was associated with numerous changes in the brain transcriptome, especially among SCZ patients on atypical APDs. Brain transcriptome data from macaques chronically treated with APDs showed that APDs affect the expression of many functionally relevant genes, some of which show expression changes in the same directions as those observed in SCZ. Co-expression modules enriched for synaptic function showed convergent patterns between SCZ and some of the APD effects, while those associated with inflammation and glucose metabolism exhibited predominantly divergent patterns between SCZ and APD effects. In contrast, major cell-type shifts inferred in SCZ were primarily unaffected by APD use. These results show that APDs may confound SCZ-associated gene expression changes in postmortem brain tissue. Disentangling these effects will help identify causal genes and improve our neurobiological understanding of SCZ.

Polygenic scores for psychiatric disorders in a diverse postmortem brain tissue cohort.

A new era of human postmortem tissue research has emerged thanks to the development of 'omics technologies that measure genes, proteins, and spatial parameters in unprecedented detail. Also newly possible is the ability to construct polygenic scores, individual-level metrics of genetic risk (also known as polygenic risk scores/PRS), based on genome-wide association studies, GWAS. Here, we report on clinical, educational, and brain gene expression correlates of polygenic scores in ancestrally diverse samples from the Human Brain Collection Core (HBCC). Genotypes from 1418 donors were subjected to quality control filters, imputed, and used to construct polygenic scores. Polygenic scores for schizophrenia predicted schizophrenia status in donors of European ancestry (p = 4.7 × 10-8, 17.2%) and in donors with African ancestry (p = 1.6 × 10-5, 10.4% of phenotypic variance explained). This pattern of higher variance explained among European ancestry samples was also observed for other psychiatric disorders (depression, bipolar disorder, substance use disorders, anxiety disorders) and for height, body mass index, and years of education. For a subset of 223 samples, gene expression from dorsolateral prefrontal cortex (DLPFC) was available through the CommonMind Consortium. In this subgroup, schizophrenia polygenic scores also predicted an aggregate gene expression score for schizophrenia (European ancestry: p = 0.0032, African ancestry: p = 0.15). Overall, polygenic scores performed as expected in ancestrally diverse samples, given historical biases toward use of European ancestry samples and variable predictive power of polygenic scores across phenotypes. The transcriptomic results reported here suggest that inherited schizophrenia genetic risk influences gene expression, even in adulthood. For future research, these and additional polygenic scores are being made available for analyses, and for selecting samples, using postmortem tissue from the Human Brain Collection Core.

Natural Killer Cell Recruitment and Activation Are Regulated by CD47 Expression in the Tumor Microenvironment.

Elevated CD47 expression in some cancers is associated with decreased survival and limited clearance by phagocytes expressing the CD47 counterreceptor SIRPα. In contrast, elevated CD47 mRNA expression in human melanomas was associated with improved survival. Gene-expression data were analyzed to determine a potential mechanism for this apparent protective function and suggested that high CD47 expression increases recruitment of natural killer (NK) cells into the tumor microenvironment. The CD47 ligand thrombospondin-1 inhibited NK cell proliferation and CD69 expression in vitro Cd47 -/- NK cells correspondingly displayed augmented effector phenotypes, indicating an inhibitory function of CD47 on NK cells. Treating human NK cells with a CD47 antibody that blocks thrombospondin-1 binding abrogated its inhibitory effect on NK cell proliferation. Similarly, treating wild-type mice with a CD47 antibody that blocks thrombospondin-1 binding delayed B16 melanoma growth, associating with increased NK cell recruitment and increased granzyme B and interferon-γ levels in intratumoral NK but not CD8+ T cells. However, B16 melanomas grew faster in Cd47 -/- than in wild-type mice. Melanoma-bearing Cd47 -/- mice exhibited decreased splenic NK cell numbers, with impaired effector protein expression and elevated exhaustion markers. Proapoptotic gene expression in Cd47-/- NK cells was associated with stress-mediated increases in mitochondrial proton leak, reactive oxygen species, and apoptosis. Global gene-expression profiling in NK cells from tumor-bearing mice identified CD47-dependent transcriptional responses that regulate systemic NK activation and exhaustion. Therefore, CD47 positively and negatively regulates NK cell function, and therapeutic antibodies that block inhibitory CD47 signaling can enhance NK immune surveillance of melanomas.

Sterol uptake and sterol biosynthesis act coordinately to mediate antifungal resistance in Candida glabrata under azole and hypoxic stress.

Pathogenic fungi, including Candida glabrata, develop strategies to grow and survive both in vitro and in vivo under azole stress. However, the mechanisms by which yeast cells counteract the inhibitory effects of azoles are not completely understood. In the current study, it was demonstrated that the expression of the ergosterol biosynthetic genes ERG2, ERG3, ERG4, ERG10, and ERG11 was significantly upregulated in C. glabrata following fluconazole treatment. Inhibiting ergosterol biosynthesis using fluconazole also increased the expression of the sterol influx transporter AUS1 and the sterol metabolism regulators SUT1 and UPC2 in fungal cells. The microarray study quantified 35 genes with elevated mRNA levels, including AUS1, TIR3, UPC2, and 8 ERG genes, in a C. glabrata mutant strain lacking ERG1, indicating that sterol importing activity is increased to compensate for defective sterol biosynthesis in cells. Bioinformatic analyses further revealed that those differentially expressed genes were involved in multiple cellular processes and biological functions, such as sterol biosynthesis, lipid localization, and sterol transport. Finally, to assess whether sterol uptake affects yeast susceptibility to azoles, we generated a C. glabrata aus1∆ mutant strain. It was shown that loss of Aus1p in C. glabrata sensitized the pathogen to azoles and enhanced the efficacy of drug exposure under low oxygen tension. In contrast, the presence of exogenous cholesterol or ergosterol in medium rendered the C. glabrata AUS1 wild­type strain highly resistant to fluconazole and voriconazole, suggesting that the sterol importing mechanism is augmented when ergosterol biosynthesis is suppressed in the cell, thus allowing C. glabrata to survive under azole pressure. On the basis of these results, it was concluded that sterol uptake and sterol biosynthesis may act coordinately and collaboratively to sustain growth and to mediate antifungal resistance in C. glabrata through dynamic gene expression in response to azole stress and environmental challenges.

CD47 Expression in Natural Killer Cells Regulates Homeostasis and Modulates Immune Response to Lymphocytic Choriomeningitis Virus.

CD47 is a ubiquitous cell surface receptor that directly regulates T cell immunity by interacting with its inhibitory ligand thrombospondin-1 and limits clearance of cells by phagocytes that express its counter-receptor signal-regulatory protein-α. Murine natural killer (NK) cells express higher levels of CD47 than other lymphocytes, but the role of CD47 in regulating NK cell homeostasis and immune function remains unclear. Cd47-/- mice exhibited depletion of NK precursors in bone marrow, consistent with the antiphagocytic function of CD47. In contrast, antisense CD47 knockdown or gene disruption resulted in a dose dependent accumulation of immature and mature NK cells in spleen. Mature Cd47-/- NK cells exhibited increased expression of NK effector and interferon gene signatures and an increased proliferative response to interleukin-15 in vitro. Cd47-/- mice showed no defect in their early response to acute Armstrong lymphocytic choriomeningitis virus (LCMV) infection but were moderately impaired in controlling chronic Clone-13 LCMV infection, which was associated with depletion of splenic NK cells and loss of effector cytokine and interferon response gene expression in Cd47-/- NK cells. Broad CD47-dependent differences in NK activation, survival, and exhaustion pathways were observed in NK cell transcriptional signatures in LCMV infected mice. These data identify CD47 as a cell-intrinsic and systemic regulator of NK cell homeostasis and NK cell function in responding to a viral infection.

A new non-radioactive deoxyhypusine synthase assay adaptable to high throughput screening.

Deoxyhypusine synthase (DHS) catalyzes the post-translational modification of eukaryotic translation factor 5A (eIF5A) by the polyamine, spermidine, that converts one specific lysine residue to deoxyhypusine [N ε -4-aminobutyl(lysine)], which is subsequently hydroxylated to hypusine [N ε -4-amino-2-hydroxybutyl(lysine)]. Hypusine synthesis represents the most critical function of polyamine. As eIF5A has been implicated in various human diseases, identification of specific inhibitors of hypusine modification is of vital importance. DHS catalyzes a complex reaction that occurs in two stages, first, the NAD-dependent cleavage of spermidine to form an enzyme-butylimine intermediate and enzyme-bound NADH, and second, the transfer of the butylimine moiety from the enzyme intermediate to the eIF5A precursor and subsequent reduction of the eIF5A-butylimine intermediate by enzyme-bound NADH to form deoxyhypusine [N ε -4-aminobutyl(lysine)]. Our data demonstrate that there is a measurable release of enzyme-bound NADH in the absence of eIF5A precursor and that the DHS activity can be determined by coupling the first phase reaction with the NADH-Glo assay in which the generation of luminescence is dependent on NADH derived from the DHS partial reaction. The conventional DHS assay that measures the incorporation of radioactivity from [1,8-3H]spermidine into the eIF5A precursor in the complete reaction cannot be readily adapted for high throughput screening (HTS). In contrast, the non-radioactive DHS/NADH-Glo coupled assay is highly specific, sensitive and reproducible and could be configured for HTS of small molecule libraries for the identification of new inhibitors of DHS. Furthermore, the coupled assay provides new insights into the dynamics of the DHS reaction especially regarding the fate of NADH.

Global quantitative proteomics reveal up-regulation of endoplasmic reticulum stress response proteins upon depletion of eIF5A in HeLa cells.

The eukaryotic translation factor, eIF5A, is a translation factor essential for protein synthesis, cell growth and animal development. By use of a adenoviral eIF5A shRNA, we have achieved an effective depletion of eIF5A in HeLa cells and undertook in vivo comprehensive proteomic analyses to examine the effects of eIF5A depletion on the total proteome and to identify cellular pathways influenced by eIF5A. The proteome of HeLa cells transduced with eIF5A shRNA was compared with that of scramble shRNA-transduced counterpart by the iTRAQ method. We identified 972 proteins consistently detected in three iTRAQ experiments and 104 proteins with significantly altered levels (protein ratio ≥1.5 or ≤0.66, p-value ≤0.05) at 72 h and/or 96 h of Ad-eIF5A-shRNA transduction. The altered expression levels of key pathway proteins were validated by western blotting. Integration of functional ontology with expression data of the 104 proteins revealed specific biological processes that are prominently up- or down-regulated. Heatmap analysis and Cytoscape visualization of biological networks identified protein folding as the major cellular process affected by depletion of eIF5A. Our unbiased, quantitative, proteomic data demonstrate that the depletion of eIF5A leads to endoplasmic reticulum stress, an unfolded protein response and up-regulation of chaperone expression in HeLa cells.

Genome-scale RNA interference screen identifies antizyme 1 (OAZ1) as a target for improvement of recombinant protein production in mammalian cells.

For the purpose of improving recombinant protein production from mammalian cells, an unbiased, high-throughput whole-genome RNA interference screen was conducted using human embryonic kidney 293 (HEK 293) cells expressing firefly luciferase. A 21,585 human genes were individually silenced with three different siRNAs for each gene. The screen identified 56 genes that led to the greatest improvement in luciferase expression. These genes were found to be included in several pathways involved in spliceosome formation and mRNA processing, transcription, metabolic processes, transport, and protein folding. The 10 genes that most enhanced protein expression when downregulated, were further confirmed by measuring the effect of their silencing on the expression of three additional recombinant proteins. Among the confirmed genes, OAZ1-the gene encoding the ornithine decarboxylase antizyme1-was selected for detailed investigation, since its silencing improved the reporter protein production without affecting cell viability. Silencing OAZ1 caused an increase of the ornithine decarboxylase enzyme and the cellular levels of putrescine and spermidine; an indication that increased cellular polyamines enhances luciferase expression without affecting its transcription. The study shows that OAZ1 is a novel target for improving expression of recombinant proteins. The genome-scale screening performed in this work can establish the foundation for targeted design of an efficient mammalian cell platform for various biotechnological applications. Biotechnol. Bioeng. 2016;113: 2403-2415. © 2016 Wiley Periodicals, Inc.

Depletion of the polyamines spermidine and spermine by overexpression of spermidine/spermine N¹-acetyltransferase 1 (SAT1) leads to mitochondria-mediated apoptosis in mammalian cells.

The polyamines putrescine, spermidine and spermine are intimately involved in the regulation of cellular growth and viability. Transduction of human embryonic kidney (HEK) 293T cells with an adenovirus encoding a key polyamine catabolic enzyme, spermidine N¹-acetyltransferase 1 (SSAT1)/SAT1 (AdSAT1), leads to a rapid depletion of spermidine and spermine, arrest in cell growth and a decline in cell viability. Annexin V/propidium iodide FACS analyses, terminal uridine nucleotide end-labelling (TUNEL) and caspase 3 assays showed a clear indication of apoptosis in AdSAT1-transduced cells (at 24-72 h), but not in cells transduced with GFP-encoding adenovirus (AdGFP). Apoptosis in the polyamine-depleted cells occurs by the mitochondrial intrinsic pathway, as evidenced by loss of mitochondrial membrane potential, increase in pro-apoptotic Bax, decrease in anti-apoptotic Bcl-xl, Bcl2 and Mcl-1 and release of cytochrome c from mitochondria, upon transduction with AdSAT1. Moreover, TEM images of AdSAT1-transduced cells revealed morphological changes commonly associated with apoptosis, including cell shrinkage, nuclear fragmentation, mitochondrial alteration, vacuolization and membrane blebbing. The apoptosis appears to result largely from depletion of the polyamines spermidine and spermine, as the polyamine analogues α-methylspermidine (α-MeSpd) and N¹,N¹²-dimethylspermine (Me2Spm) that are not substrates for SAT1 could partially restore growth and prevent apoptosis of AdSAT1-transduced cells. Inhibition of polyamine oxidases did not restore the growth of AdSAT1-transduced cells or block apoptosis, suggesting that the growth arrest and apoptosis were not induced by oxidative stress resulting from accelerated polyamine catabolism. Taken together, these data provide strong evidence that the depletion of the polyamines spermidine and spermine leads to mitochondria-mediated apoptosis.

Genome-wide analyses and functional classification of proline repeat-rich proteins: potential role of eIF5A in eukaryotic evolution.

The eukaryotic translation factor, eIF5A has been recently reported as a sequence-specific elongation factor that facilitates peptide bond formation at consecutive prolines in Saccharomyces cerevisiae, as its ortholog elongation factor P (EF-P) does in bacteria. We have searched the genome databases of 35 representative organisms from six kingdoms of life for PPP (Pro-Pro-Pro) and/or PPG (Pro-Pro-Gly)-encoding genes whose expression is expected to depend on eIF5A. We have made detailed analyses of proteome data of 5 selected species, Escherichia coli, Saccharomyces cerevisiae, Drosophila melanogaster, Mus musculus and Homo sapiens. The PPP and PPG motifs are low in the prokaryotic proteomes. However, their frequencies markedly increase with the biological complexity of eukaryotic organisms, and are higher in newly derived proteins than in those orthologous proteins commonly shared in all species. Ontology classifications of S. cerevisiae and human genes encoding the highest level of polyprolines reveal their strong association with several specific biological processes, including actin/cytoskeletal associated functions, RNA splicing/turnover, DNA binding/transcription and cell signaling. Previously reported phenotypic defects in actin polarity and mRNA decay of eIF5A mutant strains are consistent with the proposed role for eIF5A in the translation of the polyproline-containing proteins. Of all the amino acid tandem repeats (≥3 amino acids), only the proline repeat frequency correlates with functional complexity of the five organisms examined. Taken together, these findings suggest the importance of proline repeat-rich proteins and a potential role for eIF5A and its hypusine modification pathway in the course of eukaryotic evolution.

Carboxylation of cytosine (5caC) in the CG dinucleotide in the E-box motif (CGCAG|GTG) increases binding of the Tcf3|Ascl1 helix-loop-helix heterodimer 10-fold.

Three oxidative products of 5-methylcytosine (5mC) occur in mammalian genomes. We evaluated if these cytosine modifications in a CG dinucleotide altered DNA binding of four B-HLH homodimers and three heterodimers to the E-Box motif CGCAG|GTG. We examined 25 DNA probes containing all combinations of cytosine in a CG dinucleotide and none changed binding except for carboxylation of cytosine (5caC) in the strand CGCAG|GTG. 5caC enhanced binding of all examined B-HLH homodimers and heterodimers, particularly the Tcf3|Ascl1 heterodimer which increased binding ~10-fold. These results highlight a potential function of the oxidative products of 5mC, changing the DNA binding of sequence-specific transcription factors.

Depletion of cellular polyamines, spermidine and spermine, causes a total arrest in translation and growth in mammalian cells.

The polyamines, putrescine, spermidine, and spermine, are essential polycations, intimately involved in the regulation of cellular proliferation. Although polyamines exert dynamic effects on the conformation of nucleic acids and macromolecular synthesis in vitro, their specific functions in vivo are poorly understood. We investigated the cellular function of polyamines by overexpression of a key catabolic enzyme, spermidine/spermine N(1)-acetyltransferase 1 (SAT1) in mammalian cells. Transient cotransfection of HeLa cells with GFP and SAT1 vectors suppressed GFP protein expression without lowering its mRNA level, an indication that the block in GFP expression was not at transcription, but at translation. Fluorescence single-cell imaging also revealed specific inhibition of endogenous protein synthesis in the SAT1 overexpressing cells, without any inhibition of synthesis of DNA or RNA. Overexpression of SAT1 using a SAT1 adenovirus led to rapid depletion of cellular spermidine and spermine, total inhibition of protein synthesis, and growth arrest within 24 h. The SAT1 effect is most likely due to depletion of spermidine and spermine, because stable polyamine analogs that are not substrates for SAT1 restored GFP and endogenous protein synthesis. Loss of polysomes with increased 80S monosomes in the polyamine-depleted cells suggests a direct role for polyamines in translation initiation. Our data provide strong evidence for a primary function of polyamines, spermidine and spermine, in translation in mammalian cells.

A key structural domain of the Candida albicans Mdr1 protein.

A major multidrug transporter, MDR1 (multidrug resistance 1), a member of the MFS (major facilitator superfamily), invariably contributes to an increased efflux of commonly used azoles and thus corroborates their direct involvement in MDR in Candida albicans. The Mdr1 protein has two transmembrane domains, each comprising six transmembrane helices, interconnected with extracellular loops and ICLs (intracellular loops). The introduction of deletions and insertions through mutagenesis was used to address the role of the largest interdomain ICL3 of the MDR1 protein. Most of the progressive deletants, when overexpressed, eliminated the drug resistance. Notably, restoration of the length of the ICL3 by insertional mutagenesis did not restore the functionality of the protein. Interestingly, most of the insertion and deletion variants of ICL3 became amenable to trypsinization, yielding peptide fragments. The homology model of the Mdr1 protein showed that the molecular surface-charge distribution was perturbed in most of the ICL3 mutant variants. Taken together, these results provide the first evidence that the CCL (central cytoplasmic loop) of the fungal MFS transporter of the DHA1 (drug/proton antiporter) family is critical for the function of MDR. Unlike other homologous proteins, ICL3 has no apparent role in imparting substrate specificity or in the recruitment of the transporter protein.

Phospholipidome of Candida: each species of Candida has distinctive phospholipid molecular species.

By employing electrospray ionization tandem mass spectrometry (ESI-MS/MS), the phospholipidomes of eight hemiascomycetous human pathogenic Candida species have been characterized. Over 200 phospholipid molecular species were identified and quantified. There were no large differences among Candida species in phosphoglyceride class composition; however, differences in phosphoglycerides components (i.e., fatty acyl chains) were identified. In contrast, differences in sphingolipid class composition as well as in molecular species were quite evident. The phospholipid compositions of C. albicans, C. glabrata, C. parapsilosis, C. kefyr, C. tropicalis, C. dubliniensis, C. krusei, and C. utilis could be further discriminated by principal component analysis. Notwithstanding that a single strain of each species was analyzed, our data do point to a typical molecular species imprint of Candida strains.

Divergent signature motifs of nucleotide binding domains of ABC multidrug transporter, CaCdr1p of pathogenic Candida albicans, are functionally asymmetric and noninterchangeable.

Nucleotide binding domains (NBDs) of the multidrug transporter of Candida albicans, CaCdr1p, possess unique divergent amino acids in their conserved motifs. For example, NBD1 (N-terminal-NBD) possesses conserved signature motifs, while the same motif is divergent in NBD2 (C-terminal-NBD). In this study, we have evaluated the contribution of these conserved and divergent signature motifs of CaCdr1p in ATP catalysis and drug transport. By employing site-directed mutagenesis, we made three categories of mutant variants. These included mutants where all the signature motif residues were replaced with either alanines or mutants with exchanged equipositional residues to mimic the conservancy and degeneracy in opposite domain. In addition, a set of mutants where signature motifs were swapped to have variants with either both the conserved or degenerated entire signature motif. We observed that conserved and equipositional residues of NBD1 and NBD2 and swapped signature motif mutants showed high susceptibility to all the tested drugs with simultaneous abrogation in ATPase and R6G efflux activities. However, some of the mutants displayed a selective increase in susceptibility to the drugs. Notably, none of the mutant variants and WT-CaCdr1p showed any difference in drug and nucleotide binding. Our mutational analyses show not only that certain conserved residues of NBD1 signature sequence (S304, G306, and E307) are important in ATP hydrolysis and R6G efflux but also that a few divergent residues (N1002 and E1004) of NBD2 signature motif have evolved to be functionally relevant and are not interchangeable. Taken together, our data suggest that the signature motifs of CaCdr1p, whether it is divergent or conserved, are nonexchangeable and are functionally critical for ATP hydrolysis.