Complementary roles of the neuron-enriched endosomal proteins NEEP21 and calcyon in neuronal vesicle trafficking

Cell-type deconvolution of bulk-blood RNA-seq reveals biological insights into neuropsychiatric disorders

Genome-wide association studies (GWASs) have uncovered susceptibility loci associated with psychiatric disorders such as bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome, and the causal mechanisms of the link between genetic variation and disease risk is unknown. Expression quantitative trait locus (eQTL) analysis of bulk tissue is a common approach used for deciphering underlying mechanisms, although this can obscure cell-type-specific signals and thus mask trait-relevant mechanisms. Although single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell-type proportions and cell-type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-seq from 1,730 samples derived from whole blood in a cohort ascertained from individuals with BP and SCZ, this study estimated cell-type proportions and their relation with disease status and medication. For each cell type, we found between 2,875 and 4,629 eGenes (genes with an associated eQTL), including 1,211 that are not found on the basis of bulk expression alone. We performed a colocalization test between cell-type eQTLs and various traits and identified hundreds of associations that occur between cell-type eQTLs and GWASs but that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on the regulation of cell-type expression loci and found examples of genes that are differentially regulated according to lithium use. Our study suggests that applying computational methods to large bulk RNA-seq datasets of non-brain tissue can identify disease-relevant, cell-type-specific biology of psychiatric disorders and psychiatric medication.

Neuron-specific gene NSG1 binds to and positively regulates sortilin ectodomain shedding via a metalloproteinase-dependent mechanism

Increasing evidence suggests that aberrant regulation of sortilin ectodomain shedding can contribute to amyloid-β pathology and frontotemporal dementia, although the mechanism by which this occurs has not been elucidated. Here, we probed for novel binding partners of sortilin using multiple and complementary approaches and identified two proteins of the neuron-specific gene (NSG) family, NSG1 and NSG2, that physically interact and colocalize with sortilin. We show both NSG1 and NSG2 induce subcellular redistribution of sortilin to NSG1- and NSG2-enriched compartments. However, using cell surface biotinylation, we found only NSG1 reduced sortilin cell surface expression, which caused significant reductions in uptake of progranulin, a molecular determinant for frontotemporal dementia. In contrast, we demonstrate NSG2 has no effect on sortilin cell surface abundance or progranulin uptake, suggesting specificity for NSG1 in the regulation of sortilin cell surface expression. Using metalloproteinase inhibitors and A disintegrin and metalloproteinase 10 KO cells, we further show that NSG1-dependent reduction of cell surface sortilin occurred via proteolytic processing by A disintegrin and metalloproteinase 10 with a concomitant increase in shedding of sortilin ectodomain to the extracellular space. This represents a novel regulatory mechanism for sortilin ectodomain shedding that is regulated in a neuron-specific manner. Furthermore, this finding has implications for the development of strategies for brain-specific regulation of sortilin and possibly sortilin-driven pathologies.

Predicting Key Genes and Therapeutic Molecular Modelling to Explain the Association between Porphyromonas gingivalis (P. gingivalis) and Alzheimer’s Disease (AD)

The association between Porphyromonas gingivalis (P. gingivalis) and Alzheimer’s disease (AD) remains unclear. The major aim of this study was to elucidate the role of genes and molecular targets in P. gingivalis-associated AD. Two Gene Expression Omnibus (GEO) datasets, GSE5281 for AD (n = 84 Alzheimer’s, n = 74 control) and GSE9723 (n = 4 P. gingivalis, n = 4 control), were downloaded from the GEO database. Differentially expressed genes (DEGs) were obtained, and genes common to both diseases were drawn. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis was performed from the top 100 genes (50 upregulated and 50 downregulated genes). We then proceeded with CMap analysis to screen for possible small drug molecules targeting these genes. Subsequently, we performed molecular dynamics simulations. A total of 10 common genes (CALD1, HES1, ID3, PLK2, PPP2R2D, RASGRF1, SUN1, VPS33B, WTH3DI/RAB6A, and ZFP36L1) were identified with a p-value < 0.05. The PPI network of the top 100 genes showed UCHL1, SST, CHGB, CALY, and INA to be common in the MCC, DMNC, and MNC domains. Out of the 10 common genes identified, only 1 was mapped in CMap. We found three candidate small drug molecules to be a fit for PLK2, namely PubChem ID: 24971422, 11364421, and 49792852. We then performed molecular docking of PLK2 with PubChem ID: 24971422, 11364421, and 49792852. The best target, 11364421, was used to conduct the molecular dynamics simulations. The results of this study unravel novel genes to P. gingivalis-associated AD that warrant further validation.

Detecting Fear-Memory-Related Genes from Neuronal scRNA-seq Data by Diverse Distributions and Bhattacharyya Distance

The detection of differentially expressed genes (DEGs) is one of most important computational challenges in the analysis of single-cell RNA sequencing (scRNA-seq) data. However, due to the high heterogeneity and dropout noise inherent in scRNAseq data, challenges in detecting DEGs exist when using a single distribution of gene expression levels, leaving much room to improve the precision and robustness of current DEG detection methods. Here, we propose the use of a new method, DEGman, which utilizes several possible diverse distributions in combination with Bhattacharyya distance. DEGman can automatically select the best-fitting distributions of gene expression levels, and then detect DEGs by permutation testing of Bhattacharyya distances of the selected distributions from two cell groups. Compared with several popular DEG analysis tools on both large-scale simulation data and real scRNA-seq data, DEGman shows an overall improvement in the balance of sensitivity and precision. We applied DEGman to scRNA-seq data of TRAP; Ai14 mouse neurons to detect fear-memory-related genes that are significantly differentially expressed in neurons with and without fear memory. DEGman detected well-known fear-memory-related genes and many novel candidates. Interestingly, we found 25 DEGs in common in five neuron clusters that are functionally enriched for synaptic vesicles, indicating that the coupled dynamics of synaptic vesicles across in neurons plays a critical role in remote memory formation. The proposed method leverages the advantage of the use of diverse distributions in DEG analysis, exhibiting better performance in analyzing composite scRNA-seq datasets in real applications.

Heat-induced increases in body temperature in lactating dairy cows: impact on the cumulus and granulosa cell transcriptome of the periovulatory follicle

Cows acutely heat stressed after a pharmacologically induced luteinizing hormone (LH) surge had periovulatory changes in the follicular fluid proteome that may potentiate ovulation and impact oocyte developmental competence. Because the cellular origins of differentially abundant proteins were not known, we have examined the cumulus and granulosa cell transcriptomes from the periovulatory follicle in cows exhibiting varying levels of hyperthermia when occurring after the LH surge. After pharmacological induction of a dominant follicle, lactating dairy cows were administered gonadotropin releasing hormone (GnRH) and maintained in thermoneutral conditions (~67 temperature-humidity index [THI]) or heat stress conditions where THI was steadily increased for ~12 h (71 to 86 THI) and was sufficient to steadily elevate rectal temperatures. Cumulus-oocyte complexes and mural granulosa cells were recovered by transvaginal aspiration of dominant follicle content ~16 h after GnRH. Rectal temperature was used as a continuous, independent variable to identify differentially expressed genes (DEGs) increased or decreased per each 1 °C change in temperature. Cumulus (n = 9 samples) and granulosa (n = 8 samples) cells differentially expressed (false discovery rate [FDR] < 0.05) 25 and 87 genes, respectively. The majority of DEGs were upregulated by hyperthermia. Steady increases in THI are more like the "turning of a dial" than the "flipping of a switch." The moderate but impactful increases in rectal temperature induced modest fold changes in gene expression (<2-fold per 1 °C change in rectal temperature). Identification of cumulus DEGs involved in cell junctions, plasma membrane rafts, and cell-cycle regulation are consistent with marked changes in the interconnectedness and function of cumulus after the LH surge. Depending on the extent to which impacts may be occurring at the junctional level, cumulus changes may have indirect but impactful consequences on the oocyte as it undergoes meiotic maturation. Two granulosa cell DEGs have been reported by others to promote ovulation. Based on what is known, several other DEGs are suggestive of impacts on collagen formation or angiogenesis. Collectively these and other findings provide important insight regarding the extent to which the transcriptomes of the components of the periovulatory follicle (cumulus and mural granulosa cells) are affected by varying degrees of hyperthermia.

Transcriptional Profiling of Hippocampus Identifies Network Alterations in Alzheimer’s Disease

Alzheimer’s disease (AD) is a neurodegenerative disease characterized by rapid brain cell degeneration affecting different areas of the brain. Hippocampus is one of the earliest involved brain regions in the disease. Modern technologies based on high-throughput data have identified transcriptional profiling of several neurological diseases, including AD, for a better comprehension of genetic mechanisms of the disease. In this study, we investigated differentially expressed genes (DEGs) from six Gene Expression Omnibus (GEO) datasets of hippocampus of AD patients. The identified DEGs were submitted to Weighted correlation network analysis (WGCNA) and ClueGo to explore genes with a higher degree centrality and to comprehend their biological role. Subsequently, MCODE was used to identify subnetworks of interconnected DEGs. Our study found 40 down-regulated genes and 36 up-regulated genes as consensus DEGs. Analysis of the co-expression network revealed ACOT7, ATP8A2, CDC42, GAD1, GOT1, INA, NCALD, and WWTR1 to be genes with a higher degree centrality. ClueGO revealed the pathways that were mainly enriched, such as clathrin coat assembly, synaptic vesicle endocytosis, and DNA damage response signal transduction by p53 class mediator. In addition, we found a subnetwork of 12 interconnected genes (AMPH, CA10, CALY, NEFL, SNAP25, SNAP91, SNCB, STMN2, SV2B, SYN2, SYT1, and SYT13). Only CA10 and CALY are targets of known drugs while the others could be potential novel drug targets.

Global loss of Neuron-specific gene 1 causes alterations in motor coordination, increased anxiety, and diurnal hyperactivity in male mice

The Neuron-specific gene family (NSG1-3) consists of small endolysosomal proteins that are critical for trafficking multiple receptors and signaling molecules in neurons. NSG1 has been shown to play a critical role in AMPAR recycling from endosomes to plasma membrane during synaptic plasticity. However, to date nothing is known about whether NSG1 is required for normal behavior at an organismal level. Here we performed a battery of behavioral tests to determine whether loss of NSG1 would affect motor, cognitive, and/or affective behaviors, as well as circadian-related activity. Consistent with unique cerebellar expression of NSG1 among family members, we found that NSG1 was obligatory for motor coordination but not for gross motor function or learning. NSG1 knockout (KO) also altered performance across other behavioral modalities including anxiety-related and diurnal activity paradigms. Surprisingly, NSG1 KO did not cause significant impairments across all tasks within a given modality, but had specific effects within each modality. For instance, we found increases in anxiety-related behaviors in tasks with multiple stressors (e.g., elevation and exposure), but not those with a single main stressor (e.g., exposure). Interestingly, NSG1 KO animals displayed a significant increase in locomotor activity during subjective daytime, suggesting a possible impact on diurnal activity rhythms or vigilance. Surprisingly, loss of NSG1 had no effect on hippocampal-dependent learning despite previous studies showing deficits in CA1 long-term potentiation. Together, these findings do not support a role of NSG1 in hippocampal-dependent learning, but support a role in mediating proper neuronal function across amygdalar and cerebellar circuits.

Artificial intelligence and leukocyte epigenomics: Evaluation and prediction of late-onset Alzheimer's disease

We evaluated the utility of leucocyte epigenomic-biomarkers for Alzheimer’s Disease (AD) detection and elucidates its molecular pathogeneses. Genome-wide DNA methylation analysis was performed using the Infinium MethylationEPIC BeadChip array in 24 late-onset AD (LOAD) and 24 cognitively healthy subjects. Data were analyzed using six Artificial Intelligence (AI) methodologies including Deep Learning (DL) followed by Ingenuity Pathway Analysis (IPA) was used for AD prediction. We identified 152 significantly (FDR p<0.05) differentially methylated intragenic CpGs in 171 distinct genes in AD patients compared to controls. All AI platforms accurately predicted AD with AUCs ≥0.93 using 283,143 intragenic and 244,246 intergenic/extragenic CpGs. DL had an AUC = 0.99 using intragenic CpGs, with both sensitivity and specificity being 97%. High AD prediction was also achieved using intergenic/extragenic CpG sites (DL significance value being AUC = 0.99 with 97% sensitivity and specificity). Epigenetically altered genes included CR1L & CTSV (abnormal morphology of cerebral cortex), S1PR1 (CNS inflammation), and LTB4R (inflammatory response). These genes have been previously linked with AD and dementia. The differentially methylated genes CTSV & PRMT5 (ventricular hypertrophy and dilation) are linked to cardiovascular disease and of interest given the known association between impaired cerebral blood flow, cardiovascular disease, and AD. We report a novel, minimally invasive approach using peripheral blood leucocyte epigenomics, and AI analysis to detect AD and elucidate its pathogenesis.

Increased arterial pressure in mice with overexpression of the ADHD candidate gene calcyon in forebrain

The link between blood pressure (BP) and cerebral function is well established. However, it is not clear whether a common mechanism could underlie the relationship between elevated BP and cognitive deficits. The expression of calcyon, a gene abundant in catecholaminergic and hypothalamic nuclei along with other forebrain regions, is increased in the brain of the spontaneously hypertensive rat (SHR) which is a widely accepted animal model of essential hypertension and attention deficit hyperactivity disorder (ADHD). Previous studies demonstrated that mice with up-regulation of calcyon in forebrain (CalOE) exhibit deficits in working memory. To date, there is no evidence directly connecting calcyon to BP regulation. Here, we investigated whether forebrain up-regulation of calcyon alters BP using radiotelemetry. We found that CalOE mice exhibited higher mean arterial pressure (MAP) compared to tTA controls. Plasma norepinephrine levels were significantly higher in CalOE mice compared to tTA controls. Silencing the transgene with doxycycline normalized BP in CalOE mice, whereas challenging the mice with 4% high salt diet for 12 days exacerbated the MAP differences between CalOE and tTA mice. High salt diet challenge also increased proteinuria and urinary thiobarbituric acid reactive substances (TBARs) in tTA and CalOE; and the increases were more prominent in CalOE mice. Taken together, our data suggest that upregulation of calcyon in forebrain could increase BP via alterations in noradrenergic transmission and increased oxidative stress during high salt challenge. Overall, this study reveals that calcyon could be a novel neural regulator of BP raising the possibility that it could play a role in the development of vascular abnormalities.

Neuron-Specific Gene 2 (NSG2) Encodes an AMPA Receptor Interacting Protein That Modulates Excitatory Neurotransmission

Neurons have evolved a number of unique protein-coding genes that regulate trafficking of protein complexes within small organelles throughout dendrites and axons. Neuron-specific gene 2 (NSG2) encodes for one of the most abundant proteins in the nervous system during perinatal development. NSG2 belongs to a family of small neuronal endosomal proteins but its function has remained uncharacterized to date. Here, we show that NSG2 is found in discrete punctae restricted to the somatodendritic arbors of developing mouse and human neurons, and a significant proportion of NSG2 punctae colocalize with postsynaptic HOMER1 and surface-expressed AMPA-type glutamate receptors (AMPARs) at excitatory synapses. Immunoprecipitation revealed that NSG2 physically interacts with both the GluA1 and GluA2 AMPAR subunits in mouse brain. Knock-out of NSG2 in mouse hippocampal neurons selectively impaired the frequency of miniature EPSCs (mEPSCs) and caused alterations in PSD95 expression at postsynaptic densities (PSDs). In contrast, NSG2 overexpression caused a significant increase in the amplitude of mEPSCs as well as GluA2 surface expression. Thus, NSG2 functions as an AMPAR-binding protein that is required for normal synapse formation and/or maintenance, and has unique functions compared with other NSG family members.

Astroprincin (FAM171A1, C10orf38): A Regulator of Human Cell Shape and Invasive Growth

Our group originally found and cloned cDNA for a 98-kDa type 1 transmembrane glycoprotein of unknown function. Because of its abundant expression in astrocytes, it was called the protein astroprincin (APCN). Two thirds of the evolutionarily conserved protein is intracytoplasmic, whereas the extracellular domain carries two N-glycosidic side chains. APCN is physiologically expressed in placental trophoblasts, skeletal and hearth muscle, and kidney and pancreas. Overexpression of APCN (cDNA) in various cell lines induced sprouting of slender projections, whereas knockdown of APCN expression by siRNA caused disappearance of actin stress fibers. Immunohistochemical staining of human cancers for endogenous APCN showed elevated expression in invasive tumor cells compared with intratumoral cells. Human melanoma cells (SK-MEL-28) transfected with APCN cDNA acquired the ability of invasive growth in semisolid medium (Matrigel) not seen with control cells. A conserved carboxyterminal stretch of 21 amino acids was found to be essential for APCN to induce cell sprouting and invasive growth. Yeast two-hybrid screening revealed several interactive partners, of which ornithine decarboxylase antizyme-1, NEEP21 (NSG1), and ADAM10 were validated by coimmunoprecipitation. This is the first functional description of APCN. These data show that APCN regulates the dynamics of the actin cytoskeletal and, thereby, the cell shape and invasive growth potential of tumor cells.

The endosomal neuronal proteins Nsg1/NEEP21 and Nsg2/P19 are itinerant, not resident proteins of dendritic endosomes

Membrane traffic critically regulates most aspects of neuronal function. Neurons express many neuronal-specific proteins that regulate membrane traffic, including the poorly understood small transmembrane proteins neural-specific gene 1 and 2 (Nsg1/NEEP21 and Nsg2/P19). Nsg1 has been implicated in regulating endosomal recycling and sorting of several important neuronal receptors. Nsg2 is largely unstudied. At steady-state, Nsg1 and Nsg2 only partially co-localize with known endosomal compartments, and it was suggested that they mark a neuronal-specific endosome. Since Nsg1 localizes to and functions in the dendritic endosome, we set out to discover how Nsg1 and Nsg2 localization to endosomes is regulated in primary rat hippocampal neurons, using quadruple immunolocalization against endogenous proteins, live imaging of dendritic endosomes, and interference approaches against the endosomal regulators Rab5 and Rab7. In contrast to previous conclusions, we now show that Nsg1 and Nsg2 are not resident endosomal proteins, but traffic rapidly from the cell surface to lysosomes and have a half-life of less than two hours. Their partial co-localization with canonical endosomal markers thus reflects their rapid flux towards degradation rather than specific targeting to a singular compartment. These findings will require rethinking of how this class of endosomal proteins regulates trafficking of much longer-lived receptors.

Dynein binds and stimulates axonal motility of the endosome adaptor and NEEP21 family member, calcyon

The neuron-enriched, endosomal protein Calcyon (Caly) regulates endocytosis and vesicle sorting, and is important for synaptic plasticity and brain development. In the current investigation of Caly interacting proteins in brain, the microtubule retrograde motor subunit, cytoplasmic dynein 1 heavy chain (DYNC1H), and microtubule structural proteins, α and β tubulin, were identified as Caly associated proteins by MALDI-ToF/ToF. Direct interaction of the Caly-C terminus with dynein and tubulin was further confirmed in in vitro studies. In Cos-7 cells, mCherry-Caly moved along the microtubule network in organelles largely labeled by the late endosome marker Rab7. Expression of the dynein inhibitor CC1, produced striking alterations in Caly distribution, consistent with retrograde motors playing a prominent role in Caly localization and movement. In axons of cultured adult rat sensory neurons, Caly-positive organelles co-localized with dynein intermediate chain (DYNC1I1-isoform IC-1B) and the dynein regulator, lissencephaly 1 (LIS1), both of which co-precipitated from brain with the Caly C-terminus. Manipulation of dynein function in axons altered the motile properties of Caly indicating that Caly vesicles utilize the retrograde motor. Altogether, the current evidence for association with dynein motors raises the possibility that the endocytic and cargo sorting functions of Caly in neurons could be regulated by interaction with the microtubule transport system.

Axonal transport deficits in multiple sclerosis: spiraling into the abyss

The transport of mitochondria and other cellular components along the axonal microtubule cytoskeleton plays an essential role in neuronal survival. Defects in this system have been linked to a large number of neurological disorders. In multiple sclerosis (MS) and associated models such as experimental autoimmune encephalomyelitis (EAE), alterations in axonal transport have been shown to exist before neurodegeneration occurs. Genome-wide association (GWA) studies have linked several motor proteins to MS susceptibility, while neuropathological studies have shown accumulations of proteins and organelles suggestive for transport deficits. A reduced effectiveness of axonal transport can lead to neurodegeneration through inhibition of mitochondrial motility, disruption of axoglial interaction or prevention of remyelination. In MS, demyelination leads to dysregulation of axonal transport, aggravated by the effects of TNF-alpha, nitric oxide and glutamate on the cytoskeleton. The combined effect of all these pathways is a vicious cycle in which a defective axonal transport system leads to an increase in ATP consumption through loss of membrane organization and a reduction in available ATP through inhibition of mitochondrial transport, resulting in even further inhibition of transport. The persistent activity of this positive feedback loop contributes to neurodegeneration in MS.

Ctip2-, Satb2-, Prox1-, and GAD65-Expressing Neurons in Rat Cultures: Preponderance of Single- and Double-Positive Cells, and Cell Type-Specific Expression of Neuron-Specific Gene Family Members, Nsg-1 (NEEP21) and Nsg-2 (P19)

The brain consists of many distinct neuronal cell types, but which cell types are present in widely used primary cultures of embryonic rodent brain is often not known. We characterized how abundantly four cell type markers (Ctip2, Satb2, Prox1, GAD65) were represented in cultured rat neurons, how easily neurons expressing different markers can be transfected with commonly used plasmids, and whether neuronal-enriched endosomal proteins Nsg-1 (NEEP21) and Nsg-2 (P19) are ubiquitously expressed in all types of cultured neurons. We found that cultured neurons stably maintain cell type identities that are reflective of cell types in vivo. This includes neurons maintaining simultaneous expression of two transcription factors, such as Ctip2+/Satb2+ or Prox1+/Ctip2+ double-positive cells, which have also been described in vivo. Secondly, we established the superior efficiency of CAG promoters for both Lipofectamine-mediated transfection as well as for electroporation. Thirdly, we discovered that Nsg-1 and Nsg-2 were not expressed equally in all neurons: whereas high levels of both Nsg-1 and Nsg-2 were found in Satb2-, Ctip2-, and GAD65-positive neurons, Prox1-positive neurons in hippocampal cultures expressed low levels of both. Our findings thus highlight the importance of identifying neuronal cell types for doing cell biology in cultured neurons: Keeping track of neuronal cell type might uncover effects in assays that might otherwise be masked by the mixture of responsive and non-responsive neurons in the dish.