Common and Rare Variants in TMEM175 Gene Concur to the Pathogenesis of Parkinson's Disease in Italian Patients.

Parkinson's disease (PD) represents the most common neurodegenerative movement disorder. We recently identified 16 novel genes associated with PD. In this study, we focused the attention on the common and rare variants identified in the lysosomal K+ channel TMEM175. The study includes a detailed clinical and genetic analysis of 400 cases and 300 controls. Molecular studies were performed on patient-derived fibroblasts. The functional properties of the mutant channels were assessed by patch-clamp technique and co-immunoprecipitation. We have found that TMEM175 was highly expressed in dopaminergic neurons of the substantia nigra pars compacta and in microglia of the cerebral cortex of the human brain. Four common variants were associated with PD, including two novel variants rs2290402 (c.-10C > T) and rs80114247 (c.T1022C, p.M341T), located in the Kozak consensus sequence and TM3II domain, respectively. We also disclosed 13 novel highly penetrant detrimental mutations in the TMEM175 gene associated with PD. At least nine of these mutations (p.R35C, p. R183X, p.A270T, p.P308L, p.S348L, p. L405V, p.R414W, p.P427fs, p.R481W) may be sufficient to cause the disease, and the presence of mutations of other genes correlated with an earlier disease onset. In vitro functional analysis of the ion channel encoded by the mutated TMEM175 gene revealed a loss of the K+ conductance and a reduced channel affinity for Akt. Moreover, we observed an impaired autophagic/lysosomal proteolytic flux and an increase expression of unfolded protein response markers in patient-derived fibroblasts. These data suggest that mutations in TMEM175 gene may contribute to the pathophysiology of PD.

Integration of functional genomics data to uncover cell type-specific pathways affected in Parkinson's disease.

Parkinson's disease (PD) is the second most prevalent late-onset neurodegenerative disorder worldwide after Alzheimer's disease for which available drugs only deliver temporary symptomatic relief. Loss of dopaminergic neurons (DaNs) in the substantia nigra and intracellular alpha-synuclein inclusions are the main hallmarks of the disease but the events that cause this degeneration remain uncertain. Despite cell types other than DaNs such as astrocytes, microglia and oligodendrocytes have been recently associated with the pathogenesis of PD, we still lack an in-depth characterisation of PD-affected brain regions at cell-type resolution that could help our understanding of the disease mechanisms. Nevertheless, publicly available large-scale brain-specific genomic, transcriptomic and epigenomic datasets can be further exploited to extract different layers of cell type-specific biological information for the reconstruction of cell type-specific transcriptional regulatory networks. By intersecting disease risk variants within the networks, it may be possible to study the functional role of these risk variants and their combined effects at cell type- and pathway levels, that, in turn, can facilitate the identification of key regulators involved in disease progression, which are often potential therapeutic targets.

Combining multiomics and drug perturbation profiles to identify muscle-specific treatments for spinal muscular atrophy.

Spinal muscular atrophy (SMA) is a neuromuscular disorder caused by loss of survival motor neuron (SMN) protein. While SMN restoration therapies are beneficial, they are not a cure. We aimed to identify potentially novel treatments to alleviate muscle pathology combining transcriptomics, proteomics, and perturbational data sets. This revealed potential drug candidates for repurposing in SMA. One of the candidates, harmine, was further investigated in cell and animal models, improving multiple disease phenotypes, including lifespan, weight, and key molecular networks in skeletal muscle. Our work highlights the potential of multiple and parallel data-driven approaches for the development of potentially novel treatments for use in combination with SMN restoration therapies.

Targeted single-cell RNA sequencing of transcription factors enhances the identification of cell types and trajectories.

Single-cell RNA sequencing (scRNA-seq) is a widely used method for identifying cell types and trajectories in biologically heterogeneous samples, but it is limited in its detection and quantification of lowly expressed genes. This results in missing important biological signals, such as the expression of key transcription factors (TFs) driving cellular differentiation. We show that targeted sequencing of ∼1000 TFs (scCapture-seq) in iPSC-derived neuronal cultures greatly improves the biological information garnered from scRNA-seq. Increased TF resolution enhanced cell type identification, developmental trajectories, and gene regulatory networks. This allowed us to resolve differences among neuronal populations, which were generated in two different laboratories using the same differentiation protocol. ScCapture-seq improved TF-gene regulatory network inference and thus identified divergent patterns of neurogenesis into either excitatory cortical neurons or inhibitory interneurons. Furthermore, scCapture-seq revealed a role for of retinoic acid signaling in the developmental divergence between these different neuronal populations. Our results show that TF targeting improves the characterization of human cellular models and allows identification of the essential differences between cellular populations, which would otherwise be missed in traditional scRNA-seq. scCapture-seq TF targeting represents a cost-effective enhancement of scRNA-seq, which could be broadly applied to improve scRNA-seq resolution.

Dystonia genes functionally converge in specific neurons and share neurobiology with psychiatric disorders.

Dystonia is a neurological disorder characterized by sustained or intermittent muscle contractions causing abnormal movements and postures, often occurring in absence of any structural brain abnormality. Psychiatric comorbidities, including anxiety, depression, obsessive-compulsive disorder and schizophrenia, are frequent in patients with dystonia. While mutations in a fast-growing number of genes have been linked to Mendelian forms of dystonia, the cellular, anatomical, and molecular basis remains unknown for most genetic forms of dystonia, as does its genetic and biological relationship to neuropsychiatric disorders. Here we applied an unbiased systems-biology approach to explore the cellular specificity of all currently known dystonia-associated genes, predict their functional relationships, and test whether dystonia and neuropsychiatric disorders share a genetic relationship. To determine the cellular specificity of dystonia-associated genes in the brain, single-nuclear transcriptomic data derived from mouse brain was used together with expression-weighted cell-type enrichment. To identify functional relationships among dystonia-associated genes, we determined the enrichment of these genes in co-expression networks constructed from 10 human brain regions. Stratified linkage-disequilibrium score regression was used to test whether co-expression modules enriched for dystonia-associated genes significantly contribute to the heritability of anxiety, major depressive disorder, obsessive-compulsive disorder, schizophrenia, and Parkinson's disease. Dystonia-associated genes were significantly enriched in adult nigral dopaminergic neurons and striatal medium spiny neurons. Furthermore, 4 of 220 gene co-expression modules tested were significantly enriched for the dystonia-associated genes. The identified modules were derived from the substantia nigra, putamen, frontal cortex, and white matter, and were all significantly enriched for genes associated with synaptic function. Finally, we demonstrate significant enrichments of the heritability of major depressive disorder, obsessive-compulsive disorder and schizophrenia within the putamen and white matter modules, and a significant enrichment of the heritability of Parkinson's disease within the substantia nigra module. In conclusion, multiple dystonia-associated genes interact and contribute to pathogenesis likely through dysregulation of synaptic signalling in striatal medium spiny neurons, adult nigral dopaminergic neurons and frontal cortical neurons. Furthermore, the enrichment of the heritability of psychiatric disorders in the co-expression modules enriched for dystonia-associated genes indicates that psychiatric symptoms associated with dystonia are likely to be intrinsic to its pathophysiology.

A single-cell atlas of the human substantia nigra reveals cell-specific pathways associated with neurological disorders.

We describe a human single-nuclei transcriptomic atlas for the substantia nigra (SN), generated by sequencing approximately 17,000 nuclei from matched cortical and SN samples. We show that the common genetic risk for Parkinson's disease (PD) is associated with dopaminergic neuron (DaN)-specific gene expression, including mitochondrial functioning, protein folding and ubiquitination pathways. We identify a distinct cell type association between PD risk and oligodendrocyte-specific gene expression. Unlike Alzheimer's disease (AD), we find no association between PD risk and microglia or astrocytes, suggesting that neuroinflammation plays a less causal role in PD than AD. Beyond PD, we find associations between SN DaNs and GABAergic neuron gene expression and multiple neuropsychiatric disorders. Conditional analysis reveals that distinct neuropsychiatric disorders associate with distinct sets of neuron-specific genes but converge onto shared loci within oligodendrocytes and oligodendrocyte precursors. This atlas guides our aetiological understanding by associating SN cell type expression profiles with specific disease risk.

Addressing variability in iPSC-derived models of human disease: guidelines to promote reproducibility.

Induced pluripotent stem cell (iPSC) technologies have provided in vitro models of inaccessible human cell types, yielding new insights into disease mechanisms especially for neurological disorders. However, without due consideration, the thousands of new human iPSC lines generated in the past decade will inevitably affect the reproducibility of iPSC-based experiments. Differences between donor individuals, genetic stability and experimental variability contribute to iPSC model variation by impacting differentiation potency, cellular heterogeneity, morphology, and transcript and protein abundance. Such effects will confound reproducible disease modelling in the absence of appropriate strategies. In this Review, we explore the causes and effects of iPSC heterogeneity, and propose approaches to detect and account for experimental variation between studies, or even exploit it for deeper biological insight.

Cellular α-synuclein pathology is associated with bioenergetic dysfunction in Parkinson's iPSC-derived dopamine neurons.

Parkinson's disease (PD) is the second most common neurodegenerative disorder and a central role for α-synuclein (αSyn; SNCA) in disease aetiology has been proposed based on genetics and neuropathology. To better understand the pathological mechanisms of αSyn, we generated induced pluripotent stem cells (iPSCs) from healthy individuals and PD patients carrying the A53T SNCA mutation or a triplication of the SNCA locus and differentiated them into dopaminergic neurons (DAns). iPSC-derived DAn from PD patients carrying either mutation showed increased intracellular αSyn accumulation, and DAns from patients carrying the SNCA triplication displayed oligomeric αSyn pathology and elevated αSyn extracellular release. Transcriptomic analysis of purified DAns revealed perturbations in expression of genes linked to mitochondrial function, consistent with observed reduction in mitochondrial respiration, impairment in mitochondrial membrane potential, aberrant mitochondrial morphology and decreased levels of phosphorylated DRP1Ser616. Parkinson's iPSC-derived DAns showed increased endoplasmic reticulum stress and impairments in cholesterol and lipid homeostasis. Together, these data show a correlation between αSyn cellular pathology and deficits in metabolic and cellular bioenergetics in the pathology of PD.

Reproducibility of Molecular Phenotypes after Long-Term Differentiation to Human iPSC-Derived Neurons: A Multi-Site Omics Study.

Reproducibility in molecular and cellular studies is fundamental to scientific discovery. To establish the reproducibility of a well-defined long-term neuronal differentiation protocol, we repeated the cellular and molecular comparison of the same two iPSC lines across five distinct laboratories. Despite uncovering acceptable variability within individual laboratories, we detect poor cross-site reproducibility of the differential gene expression signature between these two lines. Factor analysis identifies the laboratory as the largest source of variation along with several variation-inflating confounders such as passaging effects and progenitor storage. Single-cell transcriptomics shows substantial cellular heterogeneity underlying inter-laboratory variability and being responsible for biases in differential gene expression inference. Factor analysis-based normalization of the combined dataset can remove the nuisance technical effects, enabling the execution of robust hypothesis-generating studies. Our study shows that multi-center collaborations can expose systematic biases and identify critical factors to be standardized when publishing novel protocols, contributing to increased cross-site reproducibility.

Accurate Ab Initio and Template-Based Prediction of Short Intrinsically-Disordered Regions by Bidirectional Recurrent Neural Networks Trained on Large-Scale Datasets.

Intrinsically-disordered regions lack a well-defined 3D structure, but play key roles in determining the function of many proteins. Although predictors of disorder have been shown to achieve relatively high rates of correct classification of these segments, improvements over the the years have been slow, and accurate methods are needed that are capable of accommodating the ever-increasing amount of structurally-determined protein sequences to try to boost predictive performances. In this paper, we propose a predictor for short disordered regions based on bidirectional recurrent neural networks and tested by rigorous five-fold cross-validation on a large, non-redundant dataset collected from MobiDB, a new comprehensive source of protein disorder annotations. The system exploits sequence and structural information in the forms of frequency profiles, predicted secondary structure and solvent accessibility and direct disorder annotations from homologous protein structures (templates) deposited in the Protein Data Bank. The contributions of sequence, structure and homology information result in large improvements in predictive accuracy. Additionally, the large scale of the training set leads to low false positive rates, making our systems a robust and efficient way to address high-throughput disorder prediction.

SCLpredT: Ab initio and homology-based prediction of subcellular localization by N-to-1 neural networks.

The prediction of protein subcellular localization is a important step towards the prediction of protein function, and considerable effort has gone over the last decade into the development of computational predictors of protein localization. In this article we design a new predictor of protein subcellular localization, based on a Machine Learning model (N-to-1 Neural Networks) which we have recently developed. This system, in three versions specialised, respectively, on Plants, Fungi and Animals, has a rich output which incorporates the class "organelle" alongside cytoplasm, nucleus, mitochondria and extracellular, and, additionally, chloroplast in the case of Plants. We investigate the information gain of introducing additional inputs, including predicted secondary structure, and localization information from homologous sequences. To accommodate the latter we design a new algorithm which we present here for the first time. While we do not observe any improvement when including predicted secondary structure, we measure significant overall gains when adding homology information. The final predictor including homology information correctly predicts 74%, 79% and 60% of all proteins in the case of Fungi, Animals and Plants, respectively, and outperforms our previous, state-of-the-art predictor SCLpred, and the popular predictor BaCelLo. We also observe that the contribution of homology information becomes dominant over sequence information for sequence identity values exceeding 50% for Animals and Fungi, and 60% for Plants, confirming that subcellular localization is less conserved than structure. SCLpredT is publicly available at http://distillf.ucd.ie/sclpredt/. Sequence- or template-based predictions can be obtained, and up to 32kbytes of input can be processed in a single submission.

Accurate prediction of protein enzymatic class by N-to-1 Neural Networks.

We present a novel ab initio predictor of protein enzymatic class. The predictor can classify proteins, solely based on their sequences, into one of six classes extracted from the enzyme commission (EC) classification scheme and is trained on a large, curated database of over 6,000 non-redundant proteins which we have assembled in this work. The predictor is powered by an ensemble of N-to-1 Neural Network, a novel architecture which we have recently developed. N-to-1 Neural Networks operate on the full sequence and not on predefined features. All motifs of a predefined length (31 residues in this work) are considered and are compressed by an N-to-1 Neural Network into a feature vector which is automatically determined during training. We test our predictor in 10-fold cross-validation and obtain state of the art results, with a 96% correct classification and 86% generalized correlation. All six classes are predicted with a specificity of at least 80% and false positive rates never exceeding 7%. We are currently investigating enhanced input encoding schemes which include structural information, and are analyzing trained networks to mine motifs that are most informative for the prediction, hence, likely, functionally relevant.

Classification of dendritic cell phenotypes from gene expression data.

BACKGROUND: The selection of relevant genes for sample classification is a common task in many gene expression studies. Although a number of tools have been developed to identify optimal gene expression signatures, they often generate gene lists that are too long to be exploited clinically. Consequently, researchers in the field try to identify the smallest set of genes that provide good sample classification. We investigated the genome-wide expression of the inflammatory phenotype in dendritic cells. Dendritic cells are a complex group of cells that play a critical role in vertebrate immunity. Therefore, the prediction of the inflammatory phenotype in these cells may help with the selection of immune-modulating compounds. RESULTS: A data mining protocol was applied to microarray data for murine cell lines treated with various inflammatory stimuli. The learning and validation data sets consisted of 155 and 49 samples, respectively. The data mining protocol reduced the number of probe sets from 5,802 to 10, then from 10 to 6 and finally from 6 to 3. The performances of a set of supervised classification models were compared. The best accuracy, when using the six following genes --Il12b, Cd40, Socs3, Irgm1, Plin2 and Lgals3bp-- was obtained by Tree Augmented Naïve Bayes and Nearest Neighbour (91.8%). Using the smallest set of three genes --Il12b, Cd40 and Socs3-- the performance remained satisfactory and the best accuracy was with Support Vector Machine (95.9%). These data mining models, using data for the genes Il12b, Cd40 and Socs3, were validated with a human data set consisting of 27 samples. Support Vector Machines (71.4%) and Nearest Neighbour (92.6%) gave the worst performances, but the remaining models correctly classified all the 27 samples. CONCLUSIONS: The genes selected by the data mining protocol proposed were shown to be informative for discriminating between inflammatory and steady-state phenotypes in dendritic cells. The robustness of the data mining protocol was confirmed by the accuracy for a human data set, when using only the following three genes: Il12b, Cd40 and Socs3. In summary, we analysed the longitudinal pattern of expression in dendritic cells stimulated with activating agents with the aim of identifying signatures that would predict or explain the dentritic cell response to an inflammatory agent.