Volumetric Ca2+ Imaging in the Mouse Brain Using Hybrid Multiplexed Sculpted Light Microscopy.

Calcium imaging using two-photon scanning microscopy has become an essential tool in neuroscience. However, in its typical implementation, the tradeoffs between fields of view, acquisition speeds, and depth restrictions in scattering brain tissue pose severe limitations. Here, using an integrated systems-wide optimization approach combined with multiple technical innovations, we introduce a new design paradigm for optical microscopy based on maximizing biological information while maintaining the fidelity of obtained neuron signals. Our modular design utilizes hybrid multi-photon acquisition and allows volumetric recording of neuroactivity at single-cell resolution within up to 1 × 1 × 1.22 mm volumes at up to 17 Hz in awake behaving mice. We establish the capabilities and potential of the different configurations of our imaging system at depth and across brain regions by applying it to in vivo recording of up to 12,000 neurons in mouse auditory cortex, posterior parietal cortex, and hippocampus.

Neuromodulatory Control of Long-Term Behavioral Patterns and Individuality across Development.

Animals generate complex patterns of behavior across development that may be shared or unique to individuals. Here, we examine the contributions of developmental programs and individual variation to behavior by monitoring single Caenorhabditis elegans nematodes over their complete developmental trajectories and quantifying their behavior at high spatiotemporal resolution. These measurements reveal reproducible trajectories of spontaneous foraging behaviors that are stereotyped within and between developmental stages. Dopamine, serotonin, the neuropeptide receptor NPR-1, and the TGF-ß peptide DAF-7 each have stage-specific effects on behavioral trajectories, implying the existence of a modular temporal program controlled by neuromodulators. In addition, a fraction of individuals within isogenic populations raised in controlled environments have consistent, non-genetic behavioral biases that persist across development. Several neuromodulatory systems increase or decrease the degree of non-genetic individuality to shape sustained patterns of behavior across the population.

High-Throughput, High-Resolution Mapping of Protein Localization in Mammalian Brain by In Vivo Genome Editing.

A scalable and high-throughput method to identify precise subcellular localization of endogenous proteins is essential for integrative understanding of a cell at the molecular level. Here, we developed a simple and generalizable technique to image endogenous proteins with high specificity, resolution, and contrast in single cells in mammalian brain tissue. The technique, single-cell labeling of endogenous proteins by clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9-mediated homology-directed repair (SLENDR), uses in vivo genome editing to insert a sequence encoding an epitope tag or a fluorescent protein to a gene of interest by CRISPR-Cas9-mediated homology-directed repair (HDR). Single-cell, HDR-mediated genome editing was achieved by delivering the editing machinery to dividing neuronal progenitors through in utero electroporation. We demonstrate that SLENDR allows rapid determination of the localization and dynamics of many endogenous proteins in various cell types, regions, and ages in the brain. Thus, SLENDR provides a high-throughput platform to map the subcellular localization of endogenous proteins with the resolution of micro- to nanometers in the brain.

Mapping of Brain Activity by Automated Volume Analysis of Immediate Early Genes.

Understanding how neural information is processed in physiological and pathological states would benefit from precise detection, localization, and quantification of the activity of all neurons across the entire brain, which has not, to date, been achieved in the mammalian brain. We introduce a pipeline for high-speed acquisition of brain activity at cellular resolution through profiling immediate early gene expression using immunostaining and light-sheet fluorescence imaging, followed by automated mapping and analysis of activity by an open-source software program we term ClearMap. We validate the pipeline first by analysis of brain regions activated in response to haloperidol. Next, we report new cortical regions downstream of whisker-evoked sensory processing during active exploration. Last, we combine activity mapping with axon tracing to uncover new brain regions differentially activated during parenting behavior. This pipeline is widely applicable to different experimental paradigms, including animal species for which transgenic activity reporters are not readily available.

Whole-brain imaging with single-cell resolution using chemical cocktails and computational analysis.

Systems-level identification and analysis of cellular circuits in the brain will require the development of whole-brain imaging with single-cell resolution. To this end, we performed comprehensive chemical screening to develop a whole-brain clearing and imaging method, termed CUBIC (clear, unobstructed brain imaging cocktails and computational analysis). CUBIC is a simple and efficient method involving the immersion of brain samples in chemical mixtures containing aminoalcohols, which enables rapid whole-brain imaging with single-photon excitation microscopy. CUBIC is applicable to multicolor imaging of fluorescent proteins or immunostained samples in adult brains and is scalable from a primate brain to subcellular structures. We also developed a whole-brain cell-nuclear counterstaining protocol and a computational image analysis pipeline that, together with CUBIC reagents, enable the visualization and quantification of neural activities induced by environmental stimulation. CUBIC enables time-course expression profiling of whole adult brains with single-cell resolution.

Molecular profiling of neurons based on connectivity.

The complexity and cellular heterogeneity of neural circuitry presents a major challenge to understanding the role of discrete neural populations in controlling behavior. While neuroanatomical methods enable high-resolution mapping of neural circuitry, these approaches do not allow systematic molecular profiling of neurons based on their connectivity. Here, we report the development of an approach for molecularly profiling projective neurons. We show that ribosomes can be tagged with a camelid nanobody raised against GFP and that this system can be engineered to selectively capture translating mRNAs from neurons retrogradely labeled with GFP. Using this system, we profiled neurons projecting to the nucleus accumbens. We then used an AAV to selectively profile midbrain dopamine neurons projecting to the nucleus accumbens. By comparing the captured mRNAs from each experiment, we identified a number of markers specific to VTA dopaminergic projection neurons. The current method provides a means for profiling neurons based on their projections.

Light-Sheet Microscopy in Neuroscience.

Light-sheet microscopy is an imaging approach that offers unique advantages for a diverse range of neuroscience applications. Unlike point-scanning techniques such as confocal and two-photon microscopy, light-sheet microscopes illuminate an entire plane of tissue, while imaging this plane onto a camera. Although early implementations of light sheet were optimized for longitudinal imaging of embryonic development in small specimens, emerging implementations are capable of capturing light-sheet images in freely moving, unconstrained specimens and even the intact in vivo mammalian brain. Meanwhile, the unique photobleaching and signal-to-noise benefits afforded by light-sheet microscopy's parallelized detection deliver the ability to perform volumetric imaging at much higher speeds than can be achieved using point scanning. This review describes the basic principles and evolution of light-sheet microscopy, followed by perspectives on emerging applications and opportunities for both imaging large, cleared, and expanded neural tissues and high-speed, functional imaging in vivo.

Neurodesk: an accessible, flexible and portable data analysis environment for reproducible neuroimaging.

Neuroimaging research requires purpose-built analysis software, which is challenging to install and may produce different results across computing environments. The community-oriented, open-source Neurodesk platform ( https://www.neurodesk.org/ ) harnesses a comprehensive and growing suite of neuroimaging software containers. Neurodesk includes a browser-accessible virtual desktop, command-line interface and computational notebook compatibility, allowing for accessible, flexible, portable and fully reproducible neuroimaging analysis on personal workstations, high-performance computers and the cloud.

brainlife.io: a decentralized and open-source cloud platform to support neuroscience research.

Neuroscience is advancing standardization and tool development to support rigor and transparency. Consequently, data pipeline complexity has increased, hindering FAIR (findable, accessible, interoperable and reusable) access. brainlife.io was developed to democratize neuroimaging research. The platform provides data standardization, management, visualization and processing and automatically tracks the provenance history of thousands of data objects. Here, brainlife.io is described and evaluated for validity, reliability, reproducibility, replicability and scientific utility using four data modalities and 3,200 participants.

A cortical surface template for human neuroscience.

Neuroimaging data analysis relies on normalization to standard anatomical templates to resolve macroanatomical differences across brains. Existing human cortical surface templates sample locations unevenly because of distortions introduced by inflation of the folded cortex into a standard shape. Here we present the onavg template, which affords uniform sampling of the cortex. We created the onavg template based on openly available high-quality structural scans of 1,031 brains-25 times more than existing cortical templates. We optimized the vertex locations based on cortical anatomy, achieving an even distribution. We observed consistently higher multivariate pattern classification accuracies and representational geometry inter-participant correlations based on onavg than on other templates, and onavg only needs three-quarters as much data to achieve the same performance compared with other templates. The optimized sampling also reduces CPU time across algorithms by 1.3-22.4% due to less variation in the number of vertices in each searchlight.

Causal mapping of human brain function.

Mapping human brain function is a long-standing goal of neuroscience that promises to inform the development of new treatments for brain disorders. Early maps of human brain function were based on locations of brain damage or brain stimulation that caused a functional change. Over time, this approach was largely replaced by technologies such as functional neuroimaging, which identify brain regions in which activity is correlated with behaviours or symptoms. Despite their advantages, these technologies reveal correlations, not causation. This creates challenges for interpreting the data generated from these tools and using them to develop treatments for brain disorders. A return to causal mapping of human brain function based on brain lesions and brain stimulation is underway. New approaches can combine these causal sources of information with modern neuroimaging and electrophysiology techniques to gain new insights into the functions of specific brain areas. In this Review, we provide a definition of causality for translational research, propose a continuum along which to assess the relative strength of causal information from human brain mapping studies and discuss recent advances in causal brain mapping and their relevance for developing treatments.

Longitudinal single-subject neuroimaging study reveals effects of daily environmental, physiological, and lifestyle factors on functional brain connectivity.

Our behavior and mental states are constantly shaped by our environment and experiences. However, little is known about the response of brain functional connectivity to environmental, physiological, and behavioral changes on different timescales, from days to months. This gives rise to an urgent need for longitudinal studies that collect high-frequency data. To this end, for a single subject, we collected 133 days of behavioral data with smartphones and wearables and performed 30 functional magnetic resonance imaging (fMRI) scans measuring attention, memory, resting state, and the effects of naturalistic stimuli. We find traces of past behavior and physiology in brain connectivity that extend up as far as 15 days. While sleep and physical activity relate to brain connectivity during cognitively demanding tasks, heart rate variability and respiration rate are more relevant for resting-state connectivity and movie-watching. This unique data set is openly accessible, offering an exceptional opportunity for further discoveries. Our results demonstrate that we should not study brain connectivity in isolation, but rather acknowledge its interdependence with the dynamics of the environment, changes in lifestyle, and short-term fluctuations such as transient illnesses or restless sleep. These results reflect a prolonged and sustained relationship between external factors and neural processes. Overall, precision mapping designs such as the one employed here can help to better understand intraindividual variability, which may explain some of the observed heterogeneity in fMRI findings. The integration of brain connectivity, physiology data and environmental cues will propel future environmental neuroscience research and support precision healthcare.

Integration of estimated regional gene expression with neuroimaging and clinical phenotypes at biobank scale.

An understanding of human brain individuality requires the integration of data on brain organization across people and brain regions, molecular and systems scales, as well as healthy and clinical states. Here, we help advance this understanding by leveraging methods from computational genomics to integrate large-scale genomic, transcriptomic, neuroimaging, and electronic-health record data sets. We estimated genetically regulated gene expression (gr-expression) of 18,647 genes, across 10 cortical and subcortical regions of 45,549 people from the UK Biobank. First, we showed that patterns of estimated gr-expression reflect known genetic-ancestry relationships, regional identities, as well as inter-regional correlation structure of directly assayed gene expression. Second, we performed transcriptome-wide association studies (TWAS) to discover 1,065 associations between individual variation in gr-expression and gray-matter volumes across people and brain regions. We benchmarked these associations against results from genome-wide association studies (GWAS) of the same sample and found hundreds of novel associations relative to these GWAS. Third, we integrated our results with clinical associations of gr-expression from the Vanderbilt Biobank. This integration allowed us to link genes, via gr-expression, to neuroimaging and clinical phenotypes. Fourth, we identified associations of polygenic gr-expression with structural and functional MRI phenotypes in the Human Connectome Project (HCP), a small neuroimaging-genomic data set with high-quality functional imaging data. Finally, we showed that estimates of gr-expression and magnitudes of TWAS were generally replicable and that the p-values of TWAS were replicable in large samples. Collectively, our results provide a powerful new resource for integrating gr-expression with population genetics of brain organization and disease.

Hemispheric functional organization, as revealed by naturalistic neuroimaging, in pediatric epilepsy patients with cortical resections.

Functional changes in the pediatric brain following neural injuries attest to remarkable feats of plasticity. Investigations of the neurobiological mechanisms that underlie this plasticity have largely focused on activation in the penumbra of the lesion or in contralesional, homotopic regions. Here, we adopt a whole-brain approach to evaluate the plasticity of the cortex in patients with large unilateral cortical resections due to drug-resistant childhood epilepsy. We compared the functional connectivity (FC) in patients' preserved hemisphere with the corresponding hemisphere of matched controls as they viewed and listened to a movie excerpt in a functional magnetic resonance imaging (fMRI) scanner. The preserved hemisphere was segmented into 180 and 200 parcels using two different anatomical atlases. We calculated all pairwise multivariate statistical dependencies between parcels, or parcel edges, and between 22 and 7 larger-scale functional networks, or network edges, aggregated from the smaller parcel edges. Both the left and right hemisphere-preserved patient groups had widespread reductions in FC relative to matched controls, particularly for within-network edges. A case series analysis further uncovered subclusters of patients with distinctive edgewise changes relative to controls, illustrating individual postoperative connectivity profiles. The large-scale differences in networks of the preserved hemisphere potentially reflect plasticity in the service of maintained and/or retained cognitive function.

Modeling genotype-protein interaction and correlation for Alzheimer's disease: a multi-omics imaging genetics study.

Integrating and analyzing multiple omics data sets, including genomics, proteomics and radiomics, can significantly advance researchers' comprehensive understanding of Alzheimer's disease (AD). However, current methodologies primarily focus on the main effects of genetic variation and protein, overlooking non-additive effects such as genotype-protein interaction (GPI) and correlation patterns in brain imaging genetics studies. Importantly, these non-additive effects could contribute to intermediate imaging phenotypes, finally leading to disease occurrence. In general, the interaction between genetic variations and proteins, and their correlations are two distinct biological effects, and thus disentangling the two effects for heritable imaging phenotypes is of great interest and need. Unfortunately, this issue has been largely unexploited. In this paper, to fill this gap, we propose $\textbf{M}$ulti-$\textbf{T}$ask $\textbf{G}$enotype-$\textbf{P}$rotein $\textbf{I}$nteraction and $\textbf{C}$orrelation disentangling method ($\textbf{MT-GPIC}$) to identify GPI and extract correlation patterns between them. To ensure stability and interpretability, we use novel and off-the-shelf penalties to identify meaningful genetic risk factors, as well as exploit the interconnectedness of different brain regions. Additionally, since computing GPI poses a high computational burden, we develop a fast optimization strategy for solving MT-GPIC, which is guaranteed to converge. Experimental results on the Alzheimer's Disease Neuroimaging Initiative data set show that MT-GPIC achieves higher correlation coefficients and classification accuracy than state-of-the-art methods. Moreover, our approach could effectively identify interpretable phenotype-related GPI and correlation patterns in high-dimensional omics data sets. These findings not only enhance the diagnostic accuracy but also contribute valuable insights into the underlying pathogenic mechanisms of AD.

Robust machine learning segmentation for large-scale analysis of heterogeneous clinical brain MRI datasets.

Every year, millions of brain MRI scans are acquired in hospitals, which is a figure considerably larger than the size of any research dataset. Therefore, the ability to analyze such scans could transform neuroimaging research. Yet, their potential remains untapped since no automated algorithm is robust enough to cope with the high variability in clinical acquisitions (MR contrasts, resolutions, orientations, artifacts, and subject populations). Here, we present SynthSeg+, an AI segmentation suite that enables robust analysis of heterogeneous clinical datasets. In addition to whole-brain segmentation, SynthSeg+ also performs cortical parcellation, intracranial volume estimation, and automated detection of faulty segmentations (mainly caused by scans of very low quality). We demonstrate SynthSeg+ in seven experiments, including an aging study on 14,000 scans, where it accurately replicates atrophy patterns observed on data of much higher quality. SynthSeg+ is publicly released as a ready-to-use tool to unlock the potential of quantitative morphometry.

Bias in machine learning models can be significantly mitigated by careful training: Evidence from neuroimaging studies.

Despite the great promise that machine learning has offered in many fields of medicine, it has also raised concerns about potential biases and poor generalization across genders, age distributions, races and ethnicities, hospitals, and data acquisition equipment and protocols. In the current study, and in the context of three brain diseases, we provide evidence which suggests that when properly trained, machine learning models can generalize well across diverse conditions and do not necessarily suffer from bias. Specifically, by using multistudy magnetic resonance imaging consortia for diagnosing Alzheimer's disease, schizophrenia, and autism spectrum disorder, we find that well-trained models have a high area-under-the-curve (AUC) on subjects across different subgroups pertaining to attributes such as gender, age, racial groups and different clinical studies and are unbiased under multiple fairness metrics such as demographic parity difference, equalized odds difference, equal opportunity difference, etc. We find that models that incorporate multisource data from demographic, clinical, genetic factors, and cognitive scores are also unbiased. These models have a better predictive AUC across subgroups than those trained only with imaging features, but there are also situations when these additional features do not help.

In vivo three-dimensional brain imaging with chemiluminescence probes in Alzheimer's disease models.

Optical three-dimensional (3D) molecular imaging is highly desirable for providing precise distribution of the target-of-interest in disease models. However, such 3D imaging is still far from wide applications in biomedical research; 3D brain optical molecular imaging, in particular, has rarely been reported. In this report, we designed chemiluminescence probes with high quantum yields, relatively long emission wavelengths, and high signal-to-noise ratios to fulfill the requirements for 3D brain imaging in vivo. With assistance from density-function theory (DFT) computation, we designed ADLumin-Xs by locking up the rotation of the double bond via fusing the furan ring to the phenyl ring. Our results showed that ADLumin-5 had a high quantum yield of chemiluminescence and could bind to amyloid beta (Aß). Remarkably, ADLumin-5's radiance intensity in brain areas could reach 4 × 107 photon/s/cm2/sr, which is probably 100-fold higher than most chemiluminescence probes for in vivo imaging. Because of its strong emission, we demonstrated that ADLumin-5 could be used for in vivo 3D brain imaging in transgenic mouse models of Alzheimer's disease.

A Unifying Model for Discordant and Concordant Results in Human Neuroimaging Studies of Facial Viewpoint Selectivity.

Recognizing faces regardless of their viewpoint is critical for social interactions. Traditional theories hold that view-selective early visual representations gradually become tolerant to viewpoint changes along the ventral visual hierarchy. Newer theories, based on single-neuron monkey electrophysiological recordings, suggest a three-stage architecture including an intermediate face-selective patch abruptly achieving invariance to mirror-symmetric face views. Human studies combining neuroimaging and multivariate pattern analysis (MVPA) have provided convergent evidence of view selectivity in early visual areas. However, contradictory conclusions have been reached concerning the existence in humans of a mirror-symmetric representation like that observed in macaques. We believe these contradictions arise from low-level stimulus confounds and data analysis choices. To probe for low-level confounds, we analyzed images from two face databases. Analyses of image luminance and contrast revealed biases across face views described by even polynomials-i.e., mirror-symmetric. To explain major trends across neuroimaging studies, we constructed a network model incorporating three constraints: cortical magnification, convergent feedforward projections, and interhemispheric connections. Given the identified low-level biases, we show that a gradual increase of interhemispheric connections across network-layers is sufficient to replicate view-tuning in early processing stages and mirror-symmetry in later stages. Data analysis decisions-pattern dissimilarity measure and data recentering-accounted for the inconsistent observation of mirror-symmetry across prior studies. Pattern analyses of human fMRI data (of either sex) revealed biases compatible with our model. The model provides a unifying explanation of MVPA studies of viewpoint selectivity and suggests observations of mirror-symmetry originate from ineffectively normalized signal imbalances across different face views.

Neuroimaging meta regression for coordinate based meta analysis data with a spatial model.

Coordinate-based meta-analysis combines evidence from a collection of neuroimaging studies to estimate brain activation. In such analyses, a key practical challenge is to find a computationally efficient approach with good statistical interpretability to model the locations of activation foci. In this article, we propose a generative coordinate-based meta-regression (CBMR) framework to approximate a smooth activation intensity function and investigate the effect of study-level covariates (e.g. year of publication, sample size). We employ a spline parameterization to model the spatial structure of brain activation and consider four stochastic models for modeling the random variation in foci. To examine the validity of CBMR, we estimate brain activation on 20 meta-analytic datasets, conduct spatial homogeneity tests at the voxel level, and compare the results to those generated by existing kernel-based and model-based approaches. Keywords: generalized linear models; meta-analysis; spatial statistics; statistical modeling.