NeuroPAL: A Multicolor Atlas for Whole-Brain Neuronal Identification in C. elegans.

Comprehensively resolving neuronal identities in whole-brain images is a major challenge. We achieve this in C. elegans by engineering a multicolor transgene called NeuroPAL (a neuronal polychromatic atlas of landmarks). NeuroPAL worms share a stereotypical multicolor fluorescence map for the entire hermaphrodite nervous system that resolves all neuronal identities. Neurons labeled with NeuroPAL do not exhibit fluorescence in the green, cyan, or yellow emission channels, allowing the transgene to be used with numerous reporters of gene expression or neuronal dynamics. We showcase three applications that leverage NeuroPAL for nervous-system-wide neuronal identification. First, we determine the brainwide expression patterns of all metabotropic receptors for acetylcholine, GABA, and glutamate, completing a map of this communication network. Second, we uncover changes in cell fate caused by transcription factor mutations. Third, we record brainwide activity in response to attractive and repulsive chemosensory cues, characterizing multimodal coding for these stimuli.

Molecular topography of an entire nervous system.

We have produced gene expression profiles of all 302 neurons of the C. elegans nervous system that match the single-cell resolution of its anatomy and wiring diagram. Our results suggest that individual neuron classes can be solely identified by combinatorial expression of specific gene families. For example, each neuron class expresses distinct codes of ∼23 neuropeptide genes and ∼36 neuropeptide receptors, delineating a complex and expansive "wireless" signaling network. To demonstrate the utility of this comprehensive gene expression catalog, we used computational approaches to (1) identify cis-regulatory elements for neuron-specific gene expression and (2) reveal adhesion proteins with potential roles in process placement and synaptic specificity. Our expression data are available at https://cengen.org and can be interrogated at the web application CengenApp. We expect that this neuron-specific directory of gene expression will spur investigations of underlying mechanisms that define anatomy, connectivity, and function throughout the C. elegans nervous system.

Unique homeobox codes delineate all the neuron classes of C. elegans.

It is not known at present whether neuronal cell-type diversity-defined by cell-type-specific anatomical, biophysical, functional and molecular signatures-can be reduced to relatively simple molecular descriptors of neuronal identity1. Here we show, through examination of the expression of all of the conserved homeodomain proteins encoded by the Caenorhabditis elegans genome2, that the complete set of 118 neuron classes of C. elegans can be described individually by unique combinations of the expression of homeodomain proteins, thereby providing-to our knowledge-the simplest currently known descriptor of neuronal diversity. Computational and genetic loss-of-function analyses corroborate the notion that homeodomain proteins not only provide unique descriptors of neuron type, but also have a critical role in specifying neuronal identity. We speculate that the pervasive use of homeobox genes in defining unique neuronal identities reflects the evolutionary history of neuronal cell-type specification.

Visualizing the organization and differentiation of the male-specific nervous system of C. elegans.

Sex differences in the brain are prevalent throughout the animal kingdom and particularly well appreciated in the nematode Caenorhabditis elegans, where male animals contain a little-studied set of 93 male-specific neurons. To make these neurons amenable for future study, we describe here how a multicolor reporter transgene, NeuroPAL, is capable of visualizing the distinct identities of all male-specific neurons. We used NeuroPAL to visualize and characterize a number of features of the male-specific nervous system. We provide several proofs of concept for using NeuroPAL to identify the sites of expression of gfp-tagged reporter genes and for cellular fate analysis by analyzing the effect of removal of several developmental patterning genes on neuronal identity acquisition. We use NeuroPAL and its intrinsic cohort of more than 40 distinct differentiation markers to show that, even though male-specific neurons are generated throughout all four larval stages, they execute their terminal differentiation program in a coordinated manner in the fourth larval stage. This coordinated wave of differentiation, which we call 'just-in-time' differentiation, couples neuronal maturation programs with the appearance of sexual organs.

Psychosis brain subtypes validated in first-episode cohorts and related to illness remission: results from the PHENOM consortium.

Using machine learning, we recently decomposed the neuroanatomical heterogeneity of established schizophrenia to discover two volumetric subgroups-a 'lower brain volume' subgroup (SG1) and an 'higher striatal volume' subgroup (SG2) with otherwise normal brain structure. In this study, we investigated whether the MRI signatures of these subgroups were also already present at the time of the first-episode of psychosis (FEP) and whether they were related to clinical presentation and clinical remission over 1-, 3-, and 5-years. We included 572 FEP and 424 healthy controls (HC) from 4 sites (Sao Paulo, Santander, London, Melbourne) of the PHENOM consortium. Our prior MRI subgrouping models (671 participants; USA, Germany, and China) were applied to both FEP and HC. Participants were assigned into 1 of 4 categories: subgroup 1 (SG1), subgroup 2 (SG2), no subgroup membership ('None'), and mixed SG1 + SG2 subgroups ('Mixed'). Voxel-wise analyses characterized SG1 and SG2 subgroups. Supervised machine learning analyses characterized baseline and remission signatures related to SG1 and SG2 membership. The two dominant patterns of 'lower brain volume' in SG1 and 'higher striatal volume' (with otherwise normal neuromorphology) in SG2 were identified already at the first episode of psychosis. SG1 had a significantly higher proportion of FEP (32%) vs. HC (19%) than SG2 (FEP, 21%; HC, 23%). Clinical multivariate signatures separated the SG1 and SG2 subgroups (balanced accuracy = 64%; p < 0.0001), with SG2 showing higher education but also greater positive psychosis symptoms at first presentation, and an association with symptom remission at 1-year, 5-year, and when timepoints were combined. Neuromorphological subtypes of schizophrenia are already evident at illness onset, separated by distinct clinical presentations, and differentially associated with subsequent remission. These results suggest that the subgroups may be underlying risk phenotypes that could be targeted in future treatment trials and are critical to consider when interpreting neuroimaging literature.

Interaction of the Joining Region in Junctophilin-2 With the L-Type Ca2+ Channel Is Pivotal for Cardiac Dyad Assembly and Intracellular Ca2+ Dynamics.

RATIONALE: Ca2+-induced Ca2+ release (CICR) in normal hearts requires close approximation of L-type calcium channels (LTCCs) within the transverse tubules (T-tubules) and RyR (ryanodine receptors) within the junctional sarcoplasmic reticulum. CICR is disrupted in cardiac hypertrophy and heart failure, which is associated with loss of T-tubules and disruption of cardiac dyads. In these conditions, LTCCs are redistributed from the T-tubules to disrupt CICR. The molecular mechanism responsible for LTCCs recruitment to and from the T-tubules is not well known. JPH (junctophilin) 2 enables close association between T-tubules and the junctional sarcoplasmic reticulum to ensure efficient CICR. JPH2 has a so-called joining region that is located near domains that interact with T-tubular plasma membrane, where LTCCs are housed. The idea that this joining region directly interacts with LTCCs and contributes to LTCC recruitment to T-tubules is unknown. OBJECTIVE: To determine if the joining region in JPH2 recruits LTCCs to T-tubules through direct molecular interaction in cardiomyocytes to enable efficient CICR. METHODS AND RESULTS: Modified abundance of JPH2 and redistribution of LTCC were studied in left ventricular hypertrophy in vivo and in cultured adult feline and rat ventricular myocytes. Protein-protein interaction studies showed that the joining region in JPH2 interacts with LTCC-α1C subunit and causes LTCCs distribution to the dyads, where they colocalize with RyRs. A JPH2 with induced mutations in the joining region (mutPG1JPH2) caused T-tubule remodeling and dyad loss, showing that an interaction between LTCC and JPH2 is crucial for T-tubule stabilization. mutPG1JPH2 caused asynchronous Ca2+-release with impaired excitation-contraction coupling after ß-adrenergic stimulation. The disturbed Ca2+ regulation in mutPG1JPH2 overexpressing myocytes caused calcium/calmodulin-dependent kinase II activation and altered myocyte bioenergetics. CONCLUSIONS: The interaction between LTCC and the joining region in JPH2 facilitates dyad assembly and maintains normal CICR in cardiomyocytes.

A Multidimensional Neural Maturation Index Reveals Reproducible Developmental Patterns in Children and Adolescents.

Adolescence is a time of extensive neural restructuring, leaving one susceptible to atypical development. Although neural maturation in humans can be measured using functional and structural MRI, the subtle patterns associated with the initial stages of abnormal change may be difficult to identify, particularly at an individual level. Brain age prediction models may have utility in assessing brain development in an individualized manner, as deviations between chronological age and predicted brain age could reflect one's divergence from typical development. Here, we built a support vector regression model to summarize high-dimensional neuroimaging as an index of brain age in both sexes. Using structural and functional MRI data from two large pediatric datasets and a third clinical dataset, we produced and validated a two-dimensional neural maturation index (NMI) that characterizes typical brain maturation patterns and identifies those who deviate from this trajectory. Examination of brain signatures associated with NMI scores revealed that elevated scores were related to significantly lower gray matter volume and significantly higher white matter volume, particularly in high-order regions such as the prefrontal cortex. Additionally, those with higher NMI scores exhibited enhanced connectivity in several functional brain networks, including the default mode network. Analysis of data from a sample of male and female patients with schizophrenia revealed an association between advanced NMI scores and schizophrenia diagnosis in participants aged 16-22, confirming the NMI's utility as a marker of atypicality. Altogether, our findings support the NMI as an individualized, interpretable measure by which neural development in adolescence may be assessed.SIGNIFICANCE STATEMENT The substantial neural restructuring that occurs during adolescence increases one's vulnerability to aberration. A brain index that is capable of capturing one's conformance with typical development will allow for individualized assessment and enhance our understanding of typical and atypical development. In this analysis, we produce a neural maturation index (NMI) using support vector regression and a large pediatric sample. This index generalizes across multiple cohorts and shows potential in the identification of clinical groups. We also implement a novel method for examining the developmental trajectory through data-driven analysis. The signatures identified by the NMI reflect key stages of the extensive neural development that occurs during adolescence and support its utility as a metric of typical brain development.

Two distinct neuroanatomical subtypes of schizophrenia revealed using machine learning.

Neurobiological heterogeneity in schizophrenia is poorly understood and confounds current analyses. We investigated neuroanatomical subtypes in a multi-institutional multi-ethnic cohort, using novel semi-supervised machine learning methods designed to discover patterns associated with disease rather than normal anatomical variation. Structural MRI and clinical measures in established schizophrenia (n = 307) and healthy controls (n = 364) were analysed across three sites of PHENOM (Psychosis Heterogeneity Evaluated via Dimensional Neuroimaging) consortium. Regional volumetric measures of grey matter, white matter, and CSF were used to identify distinct and reproducible neuroanatomical subtypes of schizophrenia. Two distinct neuroanatomical subtypes were found. Subtype 1 showed widespread lower grey matter volumes, most prominent in thalamus, nucleus accumbens, medial temporal, medial prefrontal/frontal and insular cortices. Subtype 2 showed increased volume in the basal ganglia and internal capsule, and otherwise normal brain volumes. Grey matter volume correlated negatively with illness duration in Subtype 1 (r = -0.201, P = 0.016) but not in Subtype 2 (r = -0.045, P = 0.652), potentially indicating different underlying neuropathological processes. The subtypes did not differ in age (t = -1.603, df = 305, P = 0.109), sex (chi-square = 0.013, df = 1, P = 0.910), illness duration (t = -0.167, df = 277, P = 0.868), antipsychotic dose (t = -0.439, df = 210, P = 0.521), age of illness onset (t = -1.355, df = 277, P = 0.177), positive symptoms (t = 0.249, df = 289, P = 0.803), negative symptoms (t = 0.151, df = 289, P = 0.879), or antipsychotic type (chi-square = 6.670, df = 3, P = 0.083). Subtype 1 had lower educational attainment than Subtype 2 (chi-square = 6.389, df = 2, P = 0.041). In conclusion, we discovered two distinct and highly reproducible neuroanatomical subtypes. Subtype 1 displayed widespread volume reduction correlating with illness duration, and worse premorbid functioning. Subtype 2 had normal and stable anatomy, except for larger basal ganglia and internal capsule, not explained by antipsychotic dose. These subtypes challenge the notion that brain volume loss is a general feature of schizophrenia and suggest differential aetiologies. They can facilitate strategies for clinical trial enrichment and stratification, and precision diagnostics.

MIDAS: Regionally linear multivariate discriminative statistical mapping.

Statistical parametric maps formed via voxel-wise mass-univariate tests, such as the general linear model, are commonly used to test hypotheses about regionally specific effects in neuroimaging cross-sectional studies where each subject is represented by a single image. Despite being informative, these techniques remain limited as they ignore multivariate relationships in the data. Most importantly, the commonly employed local Gaussian smoothing, which is important for accounting for registration errors and making the data follow Gaussian distributions, is usually chosen in an ad hoc fashion. Thus, it is often suboptimal for the task of detecting group differences and correlations with non-imaging variables. Information mapping techniques, such as searchlight, which use pattern classifiers to exploit multivariate information and obtain more powerful statistical maps, have become increasingly popular in recent years. However, existing methods may lead to important interpretation errors in practice (i.e., misidentifying a cluster as informative, or failing to detect truly informative voxels), while often being computationally expensive. To address these issues, we introduce a novel efficient multivariate statistical framework for cross-sectional studies, termed MIDAS, seeking highly sensitive and specific voxel-wise brain maps, while leveraging the power of regional discriminant analysis. In MIDAS, locally linear discriminative learning is applied to estimate the pattern that best discriminates between two groups, or predicts a variable of interest. This pattern is equivalent to local filtering by an optimal kernel whose coefficients are the weights of the linear discriminant. By composing information from all neighborhoods that contain a given voxel, MIDAS produces a statistic that collectively reflects the contribution of the voxel to the regional classifiers as well as the discriminative power of the classifiers. Critically, MIDAS efficiently assesses the statistical significance of the derived statistic by analytically approximating its null distribution without the need for computationally expensive permutation tests. The proposed framework was extensively validated using simulated atrophy in structural magnetic resonance imaging (MRI) and further tested using data from a task-based functional MRI study as well as a structural MRI study of cognitive performance. The performance of the proposed framework was evaluated against standard voxel-wise general linear models and other information mapping methods. The experimental results showed that MIDAS achieves relatively higher sensitivity and specificity in detecting group differences. Together, our results demonstrate the potential of the proposed approach to efficiently map effects of interest in both structural and functional data.

Heterogeneity of neuroanatomical patterns in prodromal Alzheimer's disease: links to cognition, progression and biomarkers.

See Coulthard and Knight (doi:10.1093/aww335) for a scientific commentary on this article.Individuals with mild cognitive impairment and Alzheimer's disease clinical diagnoses can display significant phenotypic heterogeneity. This variability likely reflects underlying genetic, environmental and neuropathological differences. Characterizing this heterogeneity is important for precision diagnostics, personalized predictions, and recruitment of relatively homogeneous sets of patients into clinical trials. In this study, we apply state-of-the-art semi-supervised machine learning methods to the Alzheimer's disease Neuroimaging cohort (ADNI) to elucidate the heterogeneity of neuroanatomical differences between subjects with mild cognitive impairment (n = 530) and Alzheimer's disease (n = 314) and cognitively normal individuals (n = 399), thereby adding to an increasing literature aiming to establish neuroanatomical and neuropathological (e.g. amyloid and tau deposition) dimensions in Alzheimer's disease and its prodromal stages. These dimensional approaches aim to provide surrogate measures of heterogeneous underlying pathologic processes leading to cognitive impairment. We relate these neuroimaging patterns to cerebrospinal fluid biomarkers, white matter hyperintensities, cognitive and clinical measures, and longitudinal trajectories. We identified four such atrophy patterns: (i) individuals with largely normal neuroanatomical profiles, who also turned out to have the least abnormal cognitive and cerebrospinal fluid biomarker profiles and the slowest clinical progression during follow-up; (ii) individuals with classical Alzheimer's disease neuroanatomical, cognitive, cerebrospinal fluid biomarkers and clinical profile, who presented the fastest clinical progression; (iii) individuals with a diffuse pattern of atrophy with relatively less pronounced involvement of the medial temporal lobe, abnormal cerebrospinal fluid amyloid-ß1-42 values, and proportionally greater executive impairment; and (iv) individuals with notably focal involvement of the medial temporal lobe and a slow steady progression, likely representing in early Alzheimer's disease stages. These four atrophy patterns effectively define a 4-dimensional categorization of neuroanatomical alterations in mild cognitive impairment and Alzheimer's disease that can complement existing dimensional approaches for staging Alzheimer's disease using a variety of biomarkers, which offer the potential for enabling precision diagnostics and prognostics, as well as targeted patient recruitment of relatively homogeneous subgroups of subjects for clinical trials.

HYDRA: Revealing heterogeneity of imaging and genetic patterns through a multiple max-margin discriminative analysis framework.

Multivariate pattern analysis techniques have been increasingly used over the past decade to derive highly sensitive and specific biomarkers of diseases on an individual basis. The driving assumption behind the vast majority of the existing methodologies is that a single imaging pattern can distinguish between healthy and diseased populations, or between two subgroups of patients (e.g., progressors vs. non-progressors). This assumption effectively ignores the ample evidence for the heterogeneous nature of brain diseases. Neurodegenerative, neuropsychiatric and neurodevelopmental disorders are largely characterized by high clinical heterogeneity, which likely stems in part from underlying neuroanatomical heterogeneity of various pathologies. Detecting and characterizing heterogeneity may deepen our understanding of disease mechanisms and lead to patient-specific treatments. However, few approaches tackle disease subtype discovery in a principled machine learning framework. To address this challenge, we present a novel non-linear learning algorithm for simultaneous binary classification and subtype identification, termed HYDRA (Heterogeneity through Discriminative Analysis). Neuroanatomical subtypes are effectively captured by multiple linear hyperplanes, which form a convex polytope that separates two groups (e.g., healthy controls from pathologic samples); each face of this polytope effectively defines a disease subtype. We validated HYDRA on simulated and clinical data. In the latter case, we applied the proposed method independently to the imaging and genetic datasets of the Alzheimer's Disease Neuroimaging Initiative (ADNI 1) study. The imaging dataset consisted of T1-weighted volumetric magnetic resonance images of 123 AD patients and 177 controls. The genetic dataset consisted of single nucleotide polymorphism information of 103 AD patients and 139 controls. We identified 3 reproducible subtypes of atrophy in AD relative to controls: (1) diffuse and extensive atrophy, (2) precuneus and extensive temporal lobe atrophy, as well some prefrontal atrophy, (3) atrophy pattern very much confined to the hippocampus and the medial temporal lobe. The genetics dataset yielded two subtypes of AD characterized mainly by the presence/absence of the apolipoprotein E (APOE) ε4 genotype, but also involving differential presence of risk alleles of CD2AP, SPON1 and LOC39095 SNPs that were associated with differences in the respective patterns of brain atrophy, especially in the precuneus. The results demonstrate the potential of the proposed approach to map disease heterogeneity in neuroimaging and genetic studies.

2P-NucTag: on-demand phototagging for molecular analysis of functionally identified cortical neurons.

Neural circuits are characterized by genetically and functionally diverse cell types. A mechanistic understanding of circuit function is predicated on linking the genetic and physiological properties of individual neurons. However, it remains highly challenging to map the functional properties of transcriptionally heterogeneous neuronal subtypes in mammalian cortical circuits in vivo. Here, we introduce a high-throughput two-photon nuclear phototagging (2P-NucTag) approach optimized for on-demand and indelible labeling of single neurons via a photoactivatable red fluorescent protein following in vivo functional characterization in behaving mice. We demonstrate the utility of this function-forward pipeline by selectively labeling and transcriptionally profiling previously inaccessible 'place' and 'silent' cells in the mouse hippocampus. Our results reveal unexpected differences in gene expression between these hippocampal pyramidal neurons with distinct spatial coding properties. Thus, 2P-NucTag opens a new way to uncover the molecular principles that govern the functional organization of neural circuits.

Learning Probabilistic Piecewise Rigid Atlases of Model Organisms via Generative Deep Networks.

Atlases are crucial to imaging statistics as they enable the standardization of inter-subject and inter-population analyses. While existing atlas estimation methods based on fluid/elastic/diffusion registration yield high-quality results for the human brain, these deformation models do not extend to a variety of other challenging areas of neuroscience such as the anatomy of C. elegans worms and fruit flies. To this end, this work presents a general probabilistic deep network-based framework for atlas estimation and registration which can flexibly incorporate various deformation models and levels of keypoint supervision that can be applied to a wide class of model organisms. Of particular relevance, it also develops a deformable piecewise rigid atlas model which is regularized to preserve inter-observation distances between neighbors. These modeling considerations are shown to improve atlas construction and key-point alignment across a diversity of datasets with small sample sizes including neuron positions in C. elegans hermaphrodites, fluorescence microscopy of male C. elegans, and images of fruit fly wings. Code is accessible at https://github.com/amin-nejat/Deformable-Atlas.

MULTIMODAL MICROSCOPY IMAGE ALIGNMENT USING SPATIAL AND SHAPE INFORMATION AND A BRANCH-AND-BOUND ALGORITHM.

Multimodal microscopy experiments that image the same population of cells under different experimental conditions have become a widely used approach in systems and molecular neuroscience. The main obstacle is to align the different imaging modalities to obtain complementary information about the observed cell population (e.g., gene expression and calcium signal). Traditional image registration methods perform poorly when only a small subset of cells are present in both images, as is common in multimodal experiments. We cast multimodal microscopy alignment as a cell subset matching problem. To solve this non-convex problem, we introduce an efficient and globally optimal branch-and-bound algorithm to find subsets of point clouds that are in rotational alignment with each other. In addition, we use complementary information about cell shape and location to compute the matching likelihood of cell pairs in two imaging modalities to further prune the optimization search tree. Finally, we use the maximal set of cells in rigid rotational alignment to seed image deformation fields to obtain a final registration result. Our framework performs better than the state-of-the-art histology alignment approaches regarding matching quality and is faster than manual alignment, providing a viable solution to improve the throughput of multimodal microscopy experiments.

ROBUST ONLINE MULTIBAND DRIFT ESTIMATION IN ELECTROPHYSIOLOGY DATA.

High-density electrophysiology probes have opened new possibilities for systems neuroscience in human and non-human animals, but probe motion poses a challenge for downstream analyses, particularly in human recordings. We improve on the state of the art for tracking this motion with four major contributions. First, we extend previous decentralized methods to use multiband information, leveraging the local field potential (LFP) in addition to spikes. Second, we show that the LFP-based approach enables registration at sub-second temporal resolution. Third, we introduce an efficient online motion tracking algorithm, enabling the method to scale up to longer and higher-resolution recordings, and possibly facilitating real-time applications. Finally, we improve the robustness of the approach by introducing a structure-aware objective and simple methods for adaptive parameter selection. Together, these advances enable fully automated scalable registration of challenging datasets from human and mouse.

Bypassing spike sorting: Density-based decoding using spike localization from dense multielectrode probes.

Neural decoding and its applications to brain computer interfaces (BCI) are essential for understanding the association between neural activity and behavior. A prerequisite for many decoding approaches is spike sorting, the assignment of action potentials (spikes) to individual neurons. Current spike sorting algorithms, however, can be inaccurate and do not properly model uncertainty of spike assignments, therefore discarding information that could potentially improve decoding performance. Recent advances in high-density probes (e.g., Neuropixels) and computational methods now allow for extracting a rich set of spike features from unsorted data; these features can in turn be used to directly decode behavioral correlates. To this end, we propose a spike sorting-free decoding method that directly models the distribution of extracted spike features using a mixture of Gaussians (MoG) encoding the uncertainty of spike assignments, without aiming to solve the spike clustering problem explicitly. We allow the mixing proportion of the MoG to change over time in response to the behavior and develop variational inference methods to fit the resulting model and to perform decoding. We benchmark our method with an extensive suite of recordings from different animals and probe geometries, demonstrating that our proposed decoder can consistently outperform current methods based on thresholding (i.e. multi-unit activity) and spike sorting. Open source code is available at https://github.com/yzhang511/density_decoding.

Neuroimaging-AI Endophenotypes of Brain Diseases in the General Population: Towards a Dimensional System of Vulnerability.

Disease heterogeneity poses a significant challenge for precision diagnostics in both clinical and sub-clinical stages. Recent work leveraging artificial intelligence (AI) has offered promise to dissect this heterogeneity by identifying complex intermediate phenotypes - herein called dimensional neuroimaging endophenotypes (DNEs) - which subtype various neurologic and neuropsychiatric diseases. We investigate the presence of nine such DNEs derived from independent yet harmonized studies on Alzheimer's disease (AD1-2)1, autism spectrum disorder (ASD1-3)2, late-life depression (LLD1-2)3, and schizophrenia (SCZ1-2)4, in the general population of 39,178 participants in the UK Biobank study. Phenome-wide associations revealed prominent associations between the nine DNEs and phenotypes related to the brain and other human organ systems. This phenotypic landscape aligns with the SNP-phenotype genome-wide associations, revealing 31 genomic loci associated with the nine DNEs (Bonferroni corrected P-value < 5×10-8/9). The DNEs exhibited significant genetic correlations, colocalization, and causal relationships with multiple human organ systems and chronic diseases. A causal effect (odds ratio=1.25 [1.11, 1.40], P-value=8.72×1-4) was established from AD2, characterized by focal medial temporal lobe atrophy, to AD. The nine DNEs and their polygenic risk scores significantly improved the prediction accuracy for 14 systemic disease categories and mortality. These findings underscore the potential of the nine DNEs to identify individuals at a high risk of developing the four brain diseases during preclinical stages for precision diagnostics. All results are publicly available at: http://labs.loni.usc.edu/medicine/.

Modified Neuropixels probes for recording human neurophysiology in the operating room.

Neuropixels are silicon-based electrophysiology-recording probes with high channel count and recording-site density. These probes offer a turnkey platform for measuring neural activity with single-cell resolution and at a scale that is beyond the capabilities of current clinically approved devices. Our team demonstrated the first-in-human use of these probes during resection surgery for epilepsy or tumors and deep brain stimulation electrode placement in patients with Parkinson's disease. Here, we provide a better understanding of the capabilities and challenges of using Neuropixels as a research tool to study human neurophysiology, with the hope that this information may inform future efforts toward regulatory approval of Neuropixels probes as research devices. In perioperative procedures, the major concerns are the initial sterility of the device, maintaining a sterile field during surgery, having multiple referencing and grounding schemes available to de-noise recordings (if necessary), protecting the silicon probe from accidental contact before insertion and obtaining high-quality action potential and local field potential recordings. The research team ensures that the device is fully operational while coordinating with the surgical team to remove sources of electrical noise that could otherwise substantially affect the signals recorded by the sensitive hardware. Prior preparation using the equipment and training in human clinical research and working in operating rooms maximize effective communication within and between the teams, ensuring high recording quality and minimizing the time added to the surgery. The perioperative procedure requires ~4 h, and the entire protocol requires multiple weeks.

Assessment of Neuroanatomical Endophenotypes of Autism Spectrum Disorder and Association With Characteristics of Individuals With Schizophrenia and the General Population.

Importance: Autism spectrum disorder (ASD) is associated with significant clinical, neuroanatomical, and genetic heterogeneity that limits precision diagnostics and treatment. Objective: To assess distinct neuroanatomical dimensions of ASD using novel semisupervised machine learning methods and to test whether the dimensions can serve as endophenotypes also in non-ASD populations. Design, Setting, and Participants: This cross-sectional study used imaging data from the publicly available Autism Brain Imaging Data Exchange (ABIDE) repositories as the discovery cohort. The ABIDE sample included individuals diagnosed with ASD aged between 16 and 64 years and age- and sex-match typically developing individuals. Validation cohorts included individuals with schizophrenia from the Psychosis Heterogeneity Evaluated via Dimensional Neuroimaging (PHENOM) consortium and individuals from the UK Biobank to represent the general population. The multisite discovery cohort included 16 internationally distributed imaging sites. Analyses were performed between March 2021 and March 2022. Main Outcomes and Measures: The trained semisupervised heterogeneity through discriminative analysis models were tested for reproducibility using extensive cross-validations. It was then applied to individuals from the PHENOM and the UK Biobank. It was hypothesized that neuroanatomical dimensions of ASD would display distinct clinical and genetic profiles and would be prominent also in non-ASD populations. Results: Heterogeneity through discriminative analysis models trained on T1-weighted brain magnetic resonance images of 307 individuals with ASD (mean [SD] age, 25.4 [9.8] years; 273 [88.9%] male) and 362 typically developing control individuals (mean [SD] age, 25.8 [8.9] years; 309 [85.4%] male) revealed that a 3-dimensional scheme was optimal to capture the ASD neuroanatomy. The first dimension (A1: aginglike) was associated with smaller brain volume, lower cognitive function, and aging-related genetic variants (FOXO3; Z = 4.65; P = 1.62 × 10-6). The second dimension (A2: schizophrenialike) was characterized by enlarged subcortical volumes, antipsychotic medication use (Cohen d = 0.65; false discovery rate-adjusted P = .048), partially overlapping genetic, neuroanatomical characteristics to schizophrenia (n = 307), and significant genetic heritability estimates in the general population (n = 14 786; mean [SD] h2, 0.71 [0.04]; P < 1 × 10-4). The third dimension (A3: typical ASD) was distinguished by enlarged cortical volumes, high nonverbal cognitive performance, and biological pathways implicating brain development and abnormal apoptosis (mean [SD] ß, 0.83 [0.02]; P = 4.22 × 10-6). Conclusions and Relevance: This cross-sectional study discovered 3-dimensional endophenotypic representation that may elucidate the heterogeneous neurobiological underpinnings of ASD to support precision diagnostics. The significant correspondence between A2 and schizophrenia indicates a possibility of identifying common biological mechanisms across the 2 mental health diagnoses.

DREDge: robust motion correction for high-density extracellular recordings across species.

High-density microelectrode arrays (MEAs) have opened new possibilities for systems neuroscience in human and non-human animals, but brain tissue motion relative to the array poses a challenge for downstream analyses, particularly in human recordings. We introduce DREDge (Decentralized Registration of Electrophysiology Data), a robust algorithm which is well suited for the registration of noisy, nonstationary extracellular electrophysiology recordings. In addition to estimating motion from spikes in the action potential (AP) frequency band, DREDge enables automated tracking of motion at high temporal resolution in the local field potential (LFP) frequency band. In human intraoperative recordings, which often feature fast (period <1s) motion, DREDge correction in the LFP band enabled reliable recovery of evoked potentials, and significantly reduced single-unit spike shape variability and spike sorting error. Applying DREDge to recordings made during deep probe insertions in nonhuman primates demonstrated the possibility of tracking probe motion of centimeters across several brain regions while simultaneously mapping single unit electrophysiological features. DREDge reliably delivered improved motion correction in acute mouse recordings, especially in those made with an recent ultra-high density probe. We also implemented a procedure for applying DREDge to recordings made across tens of days in chronic implantations in mice, reliably yielding stable motion tracking despite changes in neural activity across experimental sessions. Together, these advances enable automated, scalable registration of electrophysiological data across multiple species, probe types, and drift cases, providing a stable foundation for downstream scientific analyses of these rich datasets.