Search

Reproducible definition and identification of cell types is essential to enable investigations into their biological function and to understand their relevance in the context of development, disease and evolution. Current approaches model variability in data as continuous latent factors, followed by clustering as a separate step, or immediately apply clustering on the data. We show that such approaches can suffer from qualitative mistakes in identifying cell types robustly, particularly when the number of such cell types is in the hundreds or even thousands. Here we propose an unsupervised method, Mixture Model Inference with Discrete-coupled AutoencoderS (MMIDAS), which combines a generalized mixture model with a multi-armed deep neural network to jointly infer the discrete type and continuous type-specific variability. Using four recent datasets of brain cells spanning different technologies, species and conditions, we demonstrate that MMIDAS can identify reproducible cell types and infer cell type-dependent continuous variability in both unimodal and multimodal datasets., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

The isocortex and hippocampal formation (HPF) in the mammalian brain play critical roles in perception, cognition, emotion, and learning. We profiled ∼1.3 million cells covering the entire adult mouse isocortex and HPF and derived a transcriptomic cell-type taxonomy revealing a comprehensive repertoire of glutamatergic and GABAergic neuron types. Contrary to the traditional view of HPF as having a simpler cellular organization, we discover a complete set of glutamatergic types in HPF homologous to all major subclasses found in the six-layered isocortex, suggesting that HPF and the isocortex share a common circuit organization. We also identify large-scale continuous and graded variations of cell types along isocortical depth, across the isocortical sheet, and in multiple dimensions in hippocampus and subiculum. Overall, our study establishes a molecular architecture of the mammalian isocortex and hippocampal formation and begins to shed light on its underlying relationship with the development, evolution, connectivity, and function of these two brain structures., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Consistent identification of neurons in different experimental modalities is a key problem in neuroscience. Although methods to perform multimodal measurements in the same set of single neurons have become available, parsing complex relationships across different modalities to uncover neuronal identity is a growing challenge. Here we present an optimization framework to learn coordinated representations of multimodal data and apply it to a large multimodal dataset profiling mouse cortical interneurons. Our approach reveals strong alignment between transcriptomic and electrophysiological characterizations, enables accurate cross-modal data prediction, and identifies cell types that are consistent across modalities.

Neurons are frequently classified into distinct types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 4,200 mouse visual cortical GABAergic interneurons and reconstructed the local morphologies of 517 of those neurons. We find that most transcriptomic types (t-types) occupy specific laminar positions within visual cortex, and, for most types, the cells mapping to a t-type exhibit consistent electrophysiological and morphological properties. These properties display both discrete and continuous variation among t-types. Through multimodal integrated analysis, we define 28 met-types that have congruent morphological, electrophysiological, and transcriptomic properties and robust mutual predictability. We identify layer-specific axon innervation pattern as a defining feature distinguishing different met-types. These met-types represent a unified definition of cortical GABAergic interneuron types, providing a systematic framework to capture existing knowledge and bridge future analyses across different modalities., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

There is growing interest in the description of short-lived patterns in the spatiotemporal cortical activity monitored via neuroimaging. Most traditional analysis methods, designed to estimate relatively long-term brain dynamics, are not always appropriate to capture these patterns. Here we introduce a novel data-driven approach for detecting short-lived fMRI brain activity patterns. Exploiting Density Peak Clustering (Rodriguez and Laio [2014]), our approach reveals well localized clusters by identifying and grouping together voxels whose time-series are similar, irrespective of their brain location, even when very short time windows (∼10 volumes) are used. The method, which we call Coherence Density Peak Clustering (CDPC), is first tested on simulated data and compared with a standard unsupervised approach for fMRI analysis, independent component analysis (ICA). CDPC identifies activated voxels with essentially no false-positives and proves more reliable than ICA, which is troubled by a number of false positives comparable to that of true positives. The reliability of the method is demonstrated on real fMRI data from a simple motor task, containing brief iterations of the same movement. The clusters identified are found in regions expected to be involved in the task, and repeat synchronously with the paradigm. The methodology proposed is especially suitable for the study of short-time brain dynamics and single trial experiments, where the event or task of interest cannot be repeated for the same subject, as happens, for instance, in problem-solving, learning and decision-making. A GUI implementation of our method is available for download at https://github.com/micheleallegra/CDPC. Hum Brain Mapp 38:1421-1437, 2017. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

The free energy landscape theory has been very successful in rationalizing the folding behaviour of globular proteins, as this representation provides intuitive information on the number of states involved in the folding process, their populations and pathways of interconversion. We extend here this formalism to the case of the Aβ40 peptide, a 40-residue intrinsically disordered protein fragment associated with Alzheimer's disease. By using an advanced sampling technique that enables free energy calculations to reach convergence also in the case of highly disordered states of proteins, we provide a precise structural characterization of the free energy landscape of this peptide. We find that such landscape has inverted features with respect to those typical of folded proteins. While the global free energy minimum consists of highly disordered structures, higher free energy regions correspond to a large variety of transiently structured conformations with secondary structure elements arranged in several different manners, and are not separated from each other by sizeable free energy barriers. From this peculiar structure of the free energy landscape we predict that this peptide should become more structured and not only more compact, with increasing temperatures, and we show that this is the case through a series of biophysical measurements.

By extended atomistic simulations in explicit solvent and bias-exchange metadynamics, we study the aggregation process of 18 chains of the C-terminal segment of amyloid-β, an intrinsically disordered protein involved in Alzheimer's disease and prone to form fibrils. Starting from a disordered aggregate, we are able to observe the formation of an ordered nucleus rich in beta sheets. The rate limiting step in the nucleation pathway involves crossing a barrier of approximately 40 kcal/mol and is associated with the formation of a very specific interdigitation of the side chains belonging to different sheets. This structural pattern is different from the one observed experimentally in a microcrystal of the same system, indicating that the structure of a "nascent" fibril may differ from the one of an "extended" fibril.

Starting from a disordered aggregate, we have simulated the formation of ordered amyloid-like beta structures in a system formed by 18 polyvaline chains in explicit solvent, employing molecular dynamics accelerated by bias-exchange metadynamics. We exploited 8 different collective variables to compute the free energy of hundreds of putative aggregate structures, with variable content of parallel and antiparallel β-sheets and different packing among the sheets. This allowed characterizing in detail a possible nucleation pathway for the formation of amyloid fibrils: first the system forms a relatively large ordered nucleus of antiparallel β-sheets, and then a few parallel sheets start appearing. The relevant nucleation process culminates at this point: when a sufficient number of parallel sheets is formed, the free energy starts to decrease toward a new minimum in which this structure is predominant. The complex nucleation pathway we found cannot be described within classical nucleation theory, namely employing a unique simple reaction coordinate like the total content of β-sheets., (© 2012 American Chemical Society)