Extensible visualization and analysis for multidimensional images using Vaa3D

Imaging, Visualization, and Computation in Developmental Biology

Embryonic development is highly complex and dynamic, requiring the coordination of numerous molecular and cellular events at precise times and places. Advances in imaging technology have made it possible to follow developmental processes at cellular, tissue, and organ levels over time as they take place in the intact embryo. Parallel innovations of in vivo probes permit imaging to report on molecular, physiological, and anatomical events of embryogenesis, but the resulting multidimensional data sets pose significant challenges for extracting knowledge. In this review, we discuss recent and emerging advances in imaging technologies, in vivo labeling, and data processing that offer the greatest potential for jointly deciphering the intricate cellular dynamics and the underlying molecular mechanisms. Our discussion of the emerging area of “image-omics” highlights both the challenges of data analysis and the promise of more fully embracing computation and data science for rapidly advancing our understanding of biology.

Classification of electrophysiological and morphological neuron types in the mouse visual cortex

Understanding the diversity of cell types in the brain has been an enduring challenge and requires detailed characterization of individual neurons in multiple dimensions. To systematically profile morpho-electric properties of mammalian neurons, we established a single-cell characterization pipeline using standardized patch-clamp recordings in brain slices and biocytin-based neuronal reconstructions. We built a publicly accessible online database, the Allen Cell Types Database, to display these datasets. Intrinsic physiological properties were measured from 1,938 neurons from the adult laboratory mouse visual cortex, morphological properties were measured from 461 reconstructed neurons, and 452 neurons had both measurements available. Quantitative features were used to classify neurons into distinct types using unsupervised methods. We established a taxonomy of morphologically and electrophysiologically defined cell types for this region of the cortex, with 17 electrophysiological types, 38 morphological types and 46 morpho-electric types. There was good correspondence with previously defined transcriptomic cell types and subclasses using the same transgenic mouse lines.

A Novel Type of Neuron Within the Dorsal Striatum

The striatum is predominantly composed of medium spiny projection neurons, with the remaining neurons consisting of several types of interneurons. Among the interneurons are a group of cells that express tyrosine hydroxylase (TH). Although the intrinsic electrical properties of these TH-expressing interneurons have been characterized, there is no agreement on the number of TH-expressing cell types and their electrical properties. Here, we have used transgenic mice in which YFP-tagged channelrhodopsin-2 (ChR2) was expressed in potential TH-expressing cells in a Cre-dependent manner. We found that the YFP+ neurons in the striatum were heterogeneous in their intrinsic electrical properties; unbiased clustering indicated that there are three main neuronal subtypes. One population of neurons had aspiny dendrites with high-frequency action potential (AP) firing and plateau potentials, resembling the TH interneurons (THIN) described previously. A second, very small population of labeled neurons resembled medium-sized spiny neurons (MSN). The third population of neurons had dendrites with an intermediate density of spines, showed substantial AP adaptation and generated prolonged spikes. This type of striatal neuron has not been previously identified in the adult mouse and we have named it the Frequency-Adapting Neuron with Spines (FANS). Because of their distinctive properties, FANS may play a unique role in striatal information processing.

Extreme Compartmentalization in a Drosophila Amacrine Cell

A neuron is conventionally regarded as a single processing unit. It receives input from one or several presynaptic cells, transforms these signals, and transmits one output signal to its postsynaptic partners. Exceptions exist: amacrine cells in the mammalian retina [1-3] or interneurons in the locust mesothoracic ganglion [4] are thought to represent many electrically isolated microcircuits within one neuron. An extreme case of such an amacrine cell has recently been described in the Drosophila visual system. This cell, called CT1, reaches into two neuropils of the optic lobe, where it visits each of 700 repetitive columns, thereby covering the whole visual field [5, 6]. Due to its unusual morphology, CT1 has been suspected to perform local computations [6, 7], but this has never been proven. Using 2-photon calcium imaging and visual stimulation, we find highly compartmentalized retinotopic response properties in neighboring terminals of CT1, with each terminal acting as an independent functional unit. Model simulations demonstrate that this extreme case of compartmentalization is at the biophysical limit of neural computation.

Diverse synaptic and dendritic mechanisms of complex spike burst generation in hippocampal CA3 pyramidal cells

Complex spike bursts (CSBs) represent a characteristic firing pattern of hippocampal pyramidal cells (PCs). In CA1PCs, CSBs are driven by regenerative dendritic plateau potentials, produced by correlated entorhinal cortical and CA3 inputs that simultaneously depolarize distal and proximal dendritic domains. However, in CA3PCs neither the generation mechanisms nor the computational role of CSBs are well elucidated. We show that CSBs are induced by dendritic Ca²⁺ spikes in CA3PCs. Surprisingly, the ability of CA3PCs to produce CSBs is heterogeneous, with non-uniform synaptic input-output transformation rules triggering CSBs. The heterogeneity is partly related to the topographic position of CA3PCs; we identify two ion channel types, HCN and Kv2 channels, whose proximodistal activity gradients contribute to subregion-specific modulation of CSB propensity. Our results suggest that heterogeneous dendritic integrative properties, along with previously reported synaptic connectivity gradients, define functional subpopulations of CA3PCs that may support CA3 network computations underlying associative memory processes.

Diffusion-weighted magnetic resonance spectroscopy enables cell-specific monitoring of astrocyte reactivity in vivo

Reactive astrocytes exhibit hypertrophic morphology and altered metabolism. Deciphering astrocytic status would be of great importance to understand their role and dysregulation in pathologies, but most analytical methods remain highly invasive or destructive. The diffusion of brain metabolites, as non-invasively measured using diffusion-weighted magnetic resonance spectroscopy (DW-MRS) in vivo, depends on the structure of their micro-environment. Here we perform advanced DW-MRS in a mouse model of reactive astrocytes to determine how cellular compartments confining metabolite diffusion are changing. This reveals myo-inositol as a specific intra-astrocytic marker whose diffusion closely reflects astrocytic morphology, enabling non-invasive detection of astrocyte hypertrophy (subsequently confirmed by confocal microscopy ex vivo). Furthermore, we measure massive variations of lactate diffusion properties, suggesting that intracellular lactate is predominantly astrocytic under control conditions, but predominantly neuronal in case of astrocyte reactivity. This indicates massive remodeling of lactate metabolism, as lactate compartmentation is tightly linked to the astrocyte-to-neuron lactate shuttle mechanism.

PAT-Probabilistic Axon Tracking for Densely Labeled Neurons in Large 3-D Micrographs

A major goal of contemporary neuroscience research is to map the structural connectivity of mammalian brain using microscopy imaging data. In this context, the reconstruction of densely labeled axons from two-photon microscopy images is a challenging and important task. The visually overlapping, crossing, and often strongly distorted images of the axons allow many ambiguous interpretations to be made. We address the problem of tracking axons in densely labeled samples of neurons in large image data sets acquired from marmoset brains. Our high-resolution images were acquired using two-photon microscopy and they provided whole brain coverage, occupying terabytes of memory. Both the image distortions and the large data set size frequently make it impractical to apply present-day neuron tracing algorithms to such data due to the optimization of such algorithms to the precise tracing of either single or sparse sets of neurons. Thus, new tracking techniques are needed. We propose a probabilistic axon tracking algorithm (PAT). PAT tackles the tracking of axons in two steps: locally (L-PAT) and globally (G-PAT). L-PAT is a probabilistic tracking algorithm that can tackle distorted, cluttered images of densely labeled axons. L-PAT divides a large micrograph into smaller image stacks. It then processes each image stack independently before mapping the axons in each image to a sparse model of axon trajectories. G-PAT merges the sparse L-PAT models into a single global model of axon trajectories by minimizing a global objective function using a probabilistic optimization method. We demonstrate the superior performance of PAT over standard approaches on synthetic data. Furthermore, we successfully apply PAT to densely labeled axons in large images acquired from marmoset brains.

Automated 3-D Neuron Tracing With Precise Branch Erasing and Confidence Controlled Back Tracking

The automatic reconstruction of single neurons from microscopic images is essential to enable large-scale data-driven investigations in neuron morphology research. However, few previous methods were able to generate satisfactory results automatically from 3-D microscopic images without human intervention. In this paper, we developed a new algorithm for automatic 3-D neuron reconstruction. The main idea of the proposed algorithm is to iteratively track backward from the potential neuronal termini to the soma centre. An online confidence score is computed to decide if a tracing iteration should be stopped and discarded from the final reconstruction. The performance improvements comparing with the previous methods are mainly introduced by a more accurate estimation of the traced area and the confidence controlled back-tracking algorithm. The proposed algorithm supports large-scale batch-processing by requiring only one user specified parameter for background segmentation. We bench tested the proposed algorithm on the images obtained from both the DIADEM challenge and the BigNeuron challenge. Our proposed algorithm achieved the state-of-the-art results.

Systematic Analysis of Transmitter Coexpression Reveals Organizing Principles of Local Interneuron Heterogeneity

Broad neuronal classes are surprisingly heterogeneous across many parameters, and subclasses often exhibit partially overlapping traits including transmitter coexpression. However, the extent to which transmitter coexpression occurs in predictable, consistent patterns is unknown. Here, we demonstrate that pairwise coexpression of GABA and multiple neuropeptide families by olfactory local interneurons (LNs) of the moth Manduca sexta is highly heterogeneous, with a single LN capable of expressing neuropeptides from at least four peptide families and few instances in which neuropeptides are consistently coexpressed. Using computational modeling, we demonstrate that observed coexpression patterns cannot be explained by independent probabilities of expression of each neuropeptide. Our analyses point to three organizing principles that, once taken into consideration, allow replication of overall coexpression structure: (1) peptidergic neurons are highly likely to coexpress GABA; (2) expression probability of allatotropin depends on myoinhibitory peptide expression; and (3) the all-or-none coexpression patterns of tachykinin neurons with several other neuropeptides. For other peptide pairs, the presence of one peptide was not predictive of the presence of the other, and coexpression probability could be replicated by independent probabilities. The stochastic nature of these coexpression patterns highlights the heterogeneity of transmitter content among LNs and argues against clear-cut definition of subpopulation types based on the presence of single neuropeptides. Furthermore, the receptors for all neuropeptides and GABA were expressed within each population of principal neuron type in the antennal lobe (AL). Thus, activation of any given LN results in a dynamic cocktail of modulators that have the potential to influence every level of olfactory processing within the AL.

3D Shape Modeling for Cell Nuclear Morphological Analysis and Classification

Quantitative analysis of morphological changes in a cell nucleus is important for the understanding of nuclear architecture and its relationship with pathological conditions such as cancer. However, dimensionality of imaging data, together with a great variability of nuclear shapes, presents challenges for 3D morphological analysis. Thus, there is a compelling need for robust 3D nuclear morphometric techniques to carry out population-wide analysis. We propose a new approach that combines modeling, analysis, and interpretation of morphometric characteristics of cell nuclei and nucleoli in 3D. We used robust surface reconstruction that allows accurate approximation of 3D object boundary. Then, we computed geometric morphological measures characterizing the form of cell nuclei and nucleoli. Using these features, we compared over 450 nuclei with about 1,000 nucleoli of epithelial and mesenchymal prostate cancer cells, as well as 1,000 nuclei with over 2,000 nucleoli from serum-starved and proliferating fibroblast cells. Classification of sets of 9 and 15 cells achieved accuracy of 95.4% and 98%, respectively, for prostate cancer cells, and 95% and 98% for fibroblast cells. To our knowledge, this is the first attempt to combine these methods for 3D nuclear shape modeling and morphometry into a highly parallel pipeline workflow for morphometric analysis of thousands of nuclei and nucleoli in 3D.

Electrophysiological properties and projections of lateral hypothalamic parvalbumin positive neurons

Cracking the cytoarchitectural organization, activity patterns, and neurotransmitter nature of genetically-distinct cell types in the lateral hypothalamus (LH) is fundamental to develop a mechanistic understanding of how activity dynamics within this brain region are generated and operate together through synaptic connections to regulate circuit function. However, the precise mechanisms through which LH circuits orchestrate such dynamics have remained elusive due to the heterogeneity of the intermingled and functionally distinct cell types in this brain region. Here we reveal that a cell type in the mouse LH identified by the expression of the calcium-binding protein parvalbumin (PVALB; LHPV) is fast-spiking, releases the excitatory neurotransmitter glutamate, and sends long range projections throughout the brain. Thus, our findings challenge long-standing concepts that define neurons with a fast-spiking phenotype as exclusively GABAergic. Furthermore, we provide for the first time a detailed characterization of the electrophysiological properties of these neurons. Our work identifies LHPV neurons as a novel functional component within the LH glutamatergic circuitry.

DeepNeuron: an open deep learning toolbox for neuron tracing

Reconstructing three-dimensional (3D) morphology of neurons is essential for understanding brain structures and functions. Over the past decades, a number of neuron tracing tools including manual, semiautomatic, and fully automatic approaches have been developed to extract and analyze 3D neuronal structures. Nevertheless, most of them were developed based on coding certain rules to extract and connect structural components of a neuron, showing limited performance on complicated neuron morphology. Recently, deep learning outperforms many other machine learning methods in a wide range of image analysis and computer vision tasks. Here we developed a new Open Source toolbox, DeepNeuron, which uses deep learning networks to learn features and rules from data and trace neuron morphology in light microscopy images. DeepNeuron provides a family of modules to solve basic yet challenging problems in neuron tracing. These problems include but not limited to: (1) detecting neuron signal under different image conditions, (2) connecting neuronal signals into tree(s), (3) pruning and refining tree morphology, (4) quantifying the quality of morphology, and (5) classifying dendrites and axons in real time. We have tested DeepNeuron using light microscopy images including bright-field and confocal images of human and mouse brain, on which DeepNeuron demonstrates robustness and accuracy in neuron tracing.

Nontoxic, double-deletion-mutant rabies viral vectors for retrograde targeting of projection neurons

Recombinant rabies viral vectors have proven useful for applications including retrograde targeting of projection neurons and monosynaptic tracing, but their cytotoxicity has limited their use to short-term experiments. Here we introduce a new class of double-deletion-mutant rabies viral vectors that left transduced cells alive and healthy indefinitely. Deletion of the viral polymerase gene abolished cytotoxicity and reduced transgene expression to trace levels but left vectors still able to retrogradely infect projection neurons and express recombinases, allowing downstream expression of other transgene products such as fluorophores and calcium indicators. The morphology of retrogradely targeted cells appeared unperturbed at 1 year postinjection. Whole-cell patch-clamp recordings showed no physiological abnormalities at 8 weeks. Longitudinal two-photon structural and functional imaging in vivo, tracking thousands of individual neurons for up to 4 months, showed that transduced neurons did not die but retained stable visual response properties even at the longest time points imaged.

Radiation therapy primes tumors for nanotherapeutic delivery via macrophage-mediated vascular bursts

Efficient delivery of therapeutic nanoparticles (TNPs) to tumors is critical in improving efficacy, yet strategies that universally maximize tumoral targeting by TNP modification have been difficult to achieve in the clinic. Instead of focusing on TNP optimization, we show that the tumor microenvironment itself can be therapeutically primed to facilitate accumulation of multiple clinically relevant TNPs. Building on the recent finding that tumor-associated macrophages (TAM) can serve as nanoparticle drug depots, we demonstrate that local tumor irradiation substantially increases TAM relative to tumor cells and, thus, TNP delivery. High-resolution intravital imaging reveals that after radiation, TAM primarily accumulate adjacent to microvasculature, elicit dynamic bursts of extravasation, and subsequently enhance drug uptake in neighboring tumor cells. TAM depletion eliminates otherwise beneficial radiation effects on TNP accumulation and efficacy, and controls with unencapsulated drug show that radiation effects are more pronounced with TNPs. Priming with combined radiation and cyclophosphamide enhances vascular bursting and tumoral TNP concentration, in some cases leading to a sixfold increase of TNP accumulation in the tumor, reaching 6% of the injected dose per gram of tissue. Radiation therapy alters tumors for enhanced TNP delivery in a TAM-dependent fashion, and these observations have implications for the design of next-generation tumor-targeted nanomaterials and clinical trials for adjuvant strategies.

Neuroanatomical details of the lateral neurons of Drosophila melanogaster support their functional role in the circadian system

Drosophila melanogaster is a long‐standing model organism in the circadian clock research. A major advantage is the relative small number of about 150 neurons, which built the circadian clock in Drosophila. In our recent work, we focused on the neuroanatomical properties of the lateral neurons of the clock network. By applying the multicolor‐labeling technique Flybow we were able to identify the anatomical similarity of the previously described E2 subunit of the evening oscillator of the clock, which is built by the 5th small ventrolateral neuron (5th s‐LN_v) and one ITP positive dorsolateral neuron (LN_d). These two clock neurons share the same spatial and functional properties. We found both neurons innervating the same brain areas with similar pre‐ and postsynaptic sites in the brain. Here the anatomical findings support their shared function as a main evening oscillator in the clock network like also found in previous studies. A second quite surprising finding addresses the large lateral ventral PDF‐neurons (l‐LN_vs). We could show that the four hardly distinguishable l‐LN_vs consist of two subgroups with different innervation patterns. While three of the neurons reflect the well‐known branching pattern reproduced by PDF immunohistochemistry, one neuron per brain hemisphere has a distinguished innervation profile and is restricted only to the proximal part of the medulla‐surface. We named this neuron “extra” l‐LN_v (l‐LN_vx). We suggest the anatomical findings reflect different functional properties of the two l‐LN_v subgroups.

Mild myelin disruption elicits early alteration in behavior and proliferation in the subventricular zone

Myelin, the insulating sheath around axons, supports axon function. An important question is the impact of mild myelin disruption. In the absence of the myelin protein proteolipid protein (PLP1), myelin is generated but with age, axonal function/maintenance is disrupted. Axon disruption occurs in Plp1-null mice as early as 2 months in cortical projection neurons. High-volume cellular quantification techniques revealed a region-specific increase in oligodendrocyte density in the olfactory bulb and rostral corpus callosum that increased during adulthood. A distinct proliferative response of progenitor cells was observed in the subventricular zone (SVZ), while the number and proliferation of parenchymal oligodendrocyte progenitor cells was unchanged. This SVZ proliferative response occurred prior to evidence of axonal disruption. Thus, a novel SVZ response contributes to the region-specific increase in oligodendrocytes in Plp1-null mice. Young adult Plp1-null mice exhibited subtle but substantial behavioral alterations, indicative of an early impact of mild myelin disruption.

Automated Neuron Tracing Methods: An Updated Account

The reconstruction of neuron morphology allows to investigate how the brain works, which is one of the foremost challenges in neuroscience. This process aims at extracting the neuronal structures from microscopic imaging data. The great advances in microscopic technologies have made a huge amount of data available at the micro-, or even lower, resolution where manual inspection is time consuming, prone to error and utterly impractical. This has motivated the development of methods to automatically trace the neuronal structures, a task also known as neuron tracing. This paper surveys the latest neuron tracing methods available in the scientific literature as well as a selection of significant older papers to better place these proposals into context. They are categorized into global processing, local processing and meta-algorithm approaches. Furthermore, we point out the algorithmic components used to design each method and we report information on the datasets and the performance metrics used.

Rivulet: 3D Neuron Morphology Tracing with Iterative Back-Tracking

The digital reconstruction of single neurons from 3D confocal microscopic images is an important tool for understanding the neuron morphology and function. However the accurate automatic neuron reconstruction remains a challenging task due to the varying image quality and the complexity in the neuronal arborisation. Targeting the common challenges of neuron tracing, we propose a novel automatic 3D neuron reconstruction algorithm, named Rivulet, which is based on the multi-stencils fast-marching and iterative back-tracking. The proposed Rivulet algorithm is capable of tracing discontinuous areas without being interrupted by densely distributed noises. By evaluating the proposed pipeline with the data provided by the Diadem challenge and the recent BigNeuron project, Rivulet is shown to be robust to challenging microscopic imagestacks. We discussed the algorithm design in technical details regarding the relationships between the proposed algorithm and the other state-of-the-art neuron tracing algorithms.

Design and implementation of multi-signal and time-varying neural reconstructions

Several efficient procedures exist to digitally trace neuronal structure from light microscopy, and mature community resources have emerged to store, share, and analyze these datasets. In contrast, the quantification of intracellular distributions and morphological dynamics is not yet standardized. Current widespread descriptions of neuron morphology are static and inadequate for subcellular characterizations. We introduce a new file format to represent multichannel information as well as an open-source Vaa3D plugin to acquire this type of data. Next we define a novel data structure to capture morphological dynamics, and demonstrate its application to different time-lapse experiments. Importantly, we designed both innovations as judicious extensions of the classic SWC format, thus ensuring full back-compatibility with popular visualization and modeling tools. We then deploy the combined multichannel/time-varying reconstruction system on developing neurons in live Drosophila larvae by digitally tracing fluorescently labeled cytoskeletal components along with overall dendritic morphology as they changed over time. This same design is also suitable for quantifying dendritic calcium dynamics and tracking arbor-wide movement of any subcellular substrate of interest.

Dendritic Cytoskeletal Architecture Is Modulated by Combinatorial Transcriptional Regulation in Drosophila melanogaster

Transcription factors (TFs) have emerged as essential cell autonomous mediators of subtype specific dendritogenesis; however, the downstream effectors of these TFs remain largely unknown, as are the cellular events that TFs control to direct morphological change. As dendritic morphology is largely dictated by the organization of the actin and microtubule (MT) cytoskeletons, elucidating TF-mediated cytoskeletal regulatory programs is key to understanding molecular control of diverse dendritic morphologies. Previous studies in Drosophila melanogaster have demonstrated that the conserved TFs Cut and Knot exert combinatorial control over aspects of dendritic cytoskeleton development, promoting actin and MT-based arbor morphology, respectively. To investigate transcriptional targets of Cut and/or Knot regulation, we conducted systematic neurogenomic studies, coupled with in vivo genetic screens utilizing multi-fluor cytoskeletal and membrane marker reporters. These analyses identified a host of putative Cut and/or Knot effector molecules, and a subset of these putative TF targets converge on modulating dendritic cytoskeletal architecture, which are grouped into three major phenotypic categories, based upon neuromorphometric analyses: complexity enhancer, complexity shifter, and complexity suppressor. Complexity enhancer genes normally function to promote higher order dendritic growth and branching with variable effects on MT stabilization and F-actin organization, whereas complexity shifter and complexity suppressor genes normally function in regulating proximal-distal branching distribution or in restricting higher order branching complexity, respectively, with spatially restricted impacts on the dendritic cytoskeleton. Collectively, we implicate novel genes and cellular programs by which TFs distinctly and combinatorially govern dendritogenesis via cytoskeletal modulation.

ImageJ2: ImageJ for the next generation of scientific image data

BACKGROUND

ImageJ is an image analysis program extensively used in the biological sciences and beyond. Due to its ease of use, recordable macro language, and extensible plug-in architecture, ImageJ enjoys contributions from non-programmers, amateur programmers, and professional developers alike. Enabling such a diversity of contributors has resulted in a large community that spans the biological and physical sciences. However, a rapidly growing user base, diverging plugin suites, and technical limitations have revealed a clear need for a concerted software engineering effort to support emerging imaging paradigms, to ensure the software's ability to handle the requirements of modern science.

RESULTS

We rewrote the entire ImageJ codebase, engineering a redesigned plugin mechanism intended to facilitate extensibility at every level, with the goal of creating a more powerful tool that continues to serve the existing community while addressing a wider range of scientific requirements. This next-generation ImageJ, called "ImageJ2" in places where the distinction matters, provides a host of new functionality. It separates concerns, fully decoupling the data model from the user interface. It emphasizes integration with external applications to maximize interoperability. Its robust new plugin framework allows everything from image formats, to scripting languages, to visualization to be extended by the community. The redesigned data model supports arbitrarily large, N-dimensional datasets, which are increasingly common in modern image acquisition. Despite the scope of these changes, backwards compatibility is maintained such that this new functionality can be seamlessly integrated with the classic ImageJ interface, allowing users and developers to migrate to these new methods at their own pace.

CONCLUSIONS

Scientific imaging benefits from open-source programs that advance new method development and deployment to a diverse audience. ImageJ has continuously evolved with this idea in mind; however, new and emerging scientific requirements have posed corresponding challenges for ImageJ's development. The described improvements provide a framework engineered for flexibility, intended to support these requirements as well as accommodate future needs. Future efforts will focus on implementing new algorithms in this framework and expanding collaborations with other popular scientific software suites.

ImageJ2: ImageJ for the next generation of scientific image data

BACKGROUND

ImageJ is an image analysis program extensively used in the biological sciences and beyond. Due to its ease of use, recordable macro language, and extensible plug-in architecture, ImageJ enjoys contributions from non-programmers, amateur programmers, and professional developers alike. Enabling such a diversity of contributors has resulted in a large community that spans the biological and physical sciences. However, a rapidly growing user base, diverging plugin suites, and technical limitations have revealed a clear need for a concerted software engineering effort to support emerging imaging paradigms, to ensure the software's ability to handle the requirements of modern science.

RESULTS

We rewrote the entire ImageJ codebase, engineering a redesigned plugin mechanism intended to facilitate extensibility at every level, with the goal of creating a more powerful tool that continues to serve the existing community while addressing a wider range of scientific requirements. This next-generation ImageJ, called "ImageJ2" in places where the distinction matters, provides a host of new functionality. It separates concerns, fully decoupling the data model from the user interface. It emphasizes integration with external applications to maximize interoperability. Its robust new plugin framework allows everything from image formats, to scripting languages, to visualization to be extended by the community. The redesigned data model supports arbitrarily large, N-dimensional datasets, which are increasingly common in modern image acquisition. Despite the scope of these changes, backwards compatibility is maintained such that this new functionality can be seamlessly integrated with the classic ImageJ interface, allowing users and developers to migrate to these new methods at their own pace.

CONCLUSIONS

Scientific imaging benefits from open-source programs that advance new method development and deployment to a diverse audience. ImageJ has continuously evolved with this idea in mind; however, new and emerging scientific requirements have posed corresponding challenges for ImageJ's development. The described improvements provide a framework engineered for flexibility, intended to support these requirements as well as accommodate future needs. Future efforts will focus on implementing new algorithms in this framework and expanding collaborations with other popular scientific software suites.

Fuzzy-based propagation of prior knowledge to improve large-scale image analysis pipelines

Many automatically analyzable scientific questions are well-posed and a variety of information about expected outcomes is available a priori. Although often neglected, this prior knowledge can be systematically exploited to make automated analysis operations sensitive to a desired phenomenon or to evaluate extracted content with respect to this prior knowledge. For instance, the performance of processing operators can be greatly enhanced by a more focused detection strategy and by direct information about the ambiguity inherent in the extracted data. We present a new concept that increases the result quality awareness of image analysis operators by estimating and distributing the degree of uncertainty involved in their output based on prior knowledge. This allows the use of simple processing operators that are suitable for analyzing large-scale spatiotemporal (3D+t) microscopy images without compromising result quality. On the foundation of fuzzy set theory, we transform available prior knowledge into a mathematical representation and extensively use it to enhance the result quality of various processing operators. These concepts are illustrated on a typical bioimage analysis pipeline comprised of seed point detection, segmentation, multiview fusion and tracking. The functionality of the proposed approach is further validated on a comprehensive simulated 3D+t benchmark data set that mimics embryonic development and on large-scale light-sheet microscopy data of a zebrafish embryo. The general concept introduced in this contribution represents a new approach to efficiently exploit prior knowledge to improve the result quality of image analysis pipelines. The generality of the concept makes it applicable to practically any field with processing strategies that are arranged as linear pipelines. The automated analysis of terabyte-scale microscopy data will especially benefit from sophisticated and efficient algorithms that enable a quantitative and fast readout.

Morphological determinants of dendritic arborization neurons in Drosophila larva

Pairing in vivo imaging and computational modeling of dendritic arborization (da) neurons from the fruit fly larva provides a unique window into neuronal growth and underlying molecular processes. We image, reconstruct, and analyze the morphology of wild-type, RNAi-silenced, and mutant da neurons. We then use local and global rule-based stochastic simulations to generate artificial arbors, and identify the parameters that statistically best approximate the real data. We observe structural homeostasis in all da classes, where an increase in size of one dendritic stem is compensated by a reduction in the other stems of the same neuron. Local rule models show that bifurcation probability is determined by branch order, while branch length depends on path distance from the soma. Global rule simulations suggest that most complex morphologies tend to be constrained by resource optimization, while simpler neuron classes privilege path distance conservation. Genetic manipulations affect both the local and global optimal parameters, demonstrating functional perturbations in growth mechanisms.

Large-scale 3-dimensional quantitative imaging of tissues: state-of-the-art and translational implications

Recent developments in automated optical sectioning microscope systems have enabled researchers to conduct high resolution, three-dimensional (3D) microscopy at the scale of millimeters in various types of tissues. This powerful technology allows the exploration of tissues at an unprecedented level of detail, while preserving the spatial context. By doing so, such technology will also enable researchers to explore cellular and molecular signatures within tissue and correlate with disease course. This will allow an improved understanding of pathophysiology and facilitate a precision medicine approach to assess the response to treatment. The ability to perform large-scale imaging in 3D cannot be realized without the widespread availability of accessible quantitative analysis. In this review, we will outline recent advances in large-scale 3D imaging and discuss the available methodologies to perform meaningful analysis and potential applications in translational research.

Hedgehog signaling regulates ciliary localization of mouse odorant receptors

The ciliary localization of odorant receptors (ORs) is evolutionary conserved and essential for olfactory transduction. However, how the transport of ORs is regulated in mammalian olfactory sensory neurons is poorly understood. Here we demonstrate that odorant responsiveness and OR transport is regulated by the Hedgehog pathway. OR transport is inhibited by conditional gene inactivation of the Hedgehog signal mediator Smoothened (Smo) as well as by systemic administration of the Smo inhibitor vismodegib, a clinically used anticancer drug reported to distort smell perception in patients. The ciliary phenotype of Smo inhibition is haploinsufficient, cell autonomous, and correlates with the accumulation of OR-containing putative transport vesicles in the cytosol. The Smo-dependent OR transport route works in parallel with a low basal transport of vesicle containing both ORs and other olfactory transduction components. These findings both define a physiological function of Hedgehog signaling in olfaction and provide an important evolutionary link between olfaction and the requirement of a ciliary compartment for Hedgehog signaling.

Automatic and adaptive heterogeneous refractive index compensation for light-sheet microscopy

Optical tissue clearing has revolutionized researchers’ ability to perform fluorescent measurements of molecules, cells, and structures within intact tissue. One common complication to all optically cleared tissue is a spatially heterogeneous refractive index, leading to light scattering and first-order defocus. We designed C-DSLM (cleared tissue digital scanned light-sheet microscopy) as a low-cost method intended to automatically generate in-focus images of cleared tissue. We demonstrate the flexibility and power of C-DSLM by quantifying fluorescent features in tissue from multiple animal models using refractive index matched and mismatched microscope objectives. This includes a unique measurement of myelin tracks within intact tissue using an endogenous fluorescent reporter where typical clearing approaches render such structures difficult to image. For all measurements, we provide independent verification using standard serial tissue sectioning and quantification methods. Paired with advancements in volumetric image processing, C-DSLM provides a robust methodology to quantify sub-micron features within large tissue sections.

Optical clearing of tissue has enabled optical imaging deeper into tissue due to significantly reduced light scattering. Here, Ryan et al. tackle first-order defocus, an artefact of a non-uniform refractive index, extending light-sheet microscopy to partially cleared samples.

N3DFix: an Algorithm for Automatic Removal of Swelling Artifacts in Neuronal Reconstructions

It is well established that not only electrophysiology but also morphology plays an important role in shaping the functional properties of neurons. In order to properly quantify morphological features it is first necessary to translate observational histological data into 3-dimensional geometric reconstructions of the neuronal structures. This reconstruction process, independently of being manual or (semi-)automatic, requires several preparation steps (e.g. histological processing) before data acquisition using specialized software. Unfortunately these processing steps likely produce artifacts which are then carried to the reconstruction, such as tissue shrinkage and formation of swellings. If not accounted for and corrected, these artifacts can change significantly the results from morphometric analysis and computer simulations. Here we present N3DFix, an open-source software which uses a correction algorithm to automatically find and fix swelling artifacts in neuronal reconstructions. N3DFix works as a post-processing tool and therefore can be used in either manual or (semi-)automatic reconstructions. The algorithm's internal parameters have been defined using a "ground truth" dataset produced by a neuroanatomist, involving two complementary manual reconstruction procedures: in the first, neuronal topology was faithfully reconstructed, including all swelling artifacts; in the second procedure a meticulous correction of the artifacts was manually performed directly during neuronal tracing. The internal parameters of N3DFix were set to minimize the differences between manual amendments and the algorithm's corrections. It is shown that the performance of N3DFix is comparable to careful manual correction of the swelling artifacts. To promote easy access and wide adoption, N3DFix is available in NEURON, Vaa3D and Py3DN.

Three-dimensional Imaging and Scanning: Current and Future Applications for Pathology

Imaging is vital for the assessment of physiologic and phenotypic details. In the past, biomedical imaging was heavily reliant on analog, low-throughput methods, which would produce two-dimensional images. However, newer, digital, and high-throughput three-dimensional (3D) imaging methods, which rely on computer vision and computer graphics, are transforming the way biomedical professionals practice. 3D imaging has been useful in diagnostic, prognostic, and therapeutic decision-making for the medical and biomedical professions. Herein, we summarize current imaging methods that enable optimal 3D histopathologic reconstruction: Scanning, 3D scanning, and whole slide imaging. Briefly mentioned are emerging platforms, which combine robotics, sectioning, and imaging in their pursuit to digitize and automate the entire microscopy workflow. Finally, both current and emerging 3D imaging methods are discussed in relation to current and future applications within the context of pathology.

SmartScope2: Simultaneous Imaging and Reconstruction of Neuronal Morphology

Quantitative analysis of neuronal morphology is critical in cell type classification and for deciphering how structure gives rise to function in the brain. Most current approaches to imaging and tracing neuronal 3D morphology are data intensive. We introduce SmartScope2, the first open source, automated neuron reconstruction machine integrating online image analysis with automated multiphoton imaging. SmartScope2 takes advantage of a neuron’s sparse morphology to improve imaging speed and reduce image data stored, transferred and analyzed. We show that SmartScope2 is able to produce the complex 3D morphology of human and mouse cortical neurons with six-fold reduction in image data requirements and three times the imaging speed compared to conventional methods.

TDat: An Efficient Platform for Processing Petabyte-Scale Whole-Brain Volumetric Images

Three-dimensional imaging of whole mammalian brains at single-neuron resolution has generated terabyte (TB)- and even petabyte (PB)-sized datasets. Due to their size, processing these massive image datasets can be hindered by the computer hardware and software typically found in biological laboratories. To fill this gap, we have developed an efficient platform named TDat, which adopts a novel data reformatting strategy by reading cuboid data and employing parallel computing. In data reformatting, TDat is more efficient than any other software. In data accessing, we adopted parallelization to fully explore the capability for data transmission in computers. We applied TDat in large-volume data rigid registration and neuron tracing in whole-brain data with single-neuron resolution, which has never been demonstrated in other studies. We also showed its compatibility with various computing platforms, image processing software and imaging systems.

Optical clearing and fluorescence deep-tissue imaging for 3D quantitative analysis of the brain tumor microenvironment

BACKGROUND

Three-dimensional visualization of the brain vasculature and its interactions with surrounding cells may shed light on diseases where aberrant microvascular organization is involved, including glioblastoma (GBM). Intravital confocal imaging allows 3D visualization of microvascular structures and migration of cells in the brain of mice, however, with limited imaging depth. To enable comprehensive analysis of GBM and the brain microenvironment, in-depth 3D imaging methods are needed. Here, we employed methods for optical tissue clearing prior to 3D microscopy to visualize the brain microvasculature and routes of invasion of GBM cells.

METHODS

We present a workflow for ex vivo imaging of optically cleared brain tumor tissues and subsequent computational modeling. This workflow was used for quantification of the microvasculature in relation to nuclear or cellular density in healthy mouse brain tissues and in human orthotopic, infiltrative GBM8 and E98 glioblastoma models.

RESULTS

Ex vivo cleared mouse brain tissues had a >10-fold imaging depth as compared to intravital imaging of mouse brain in vivo. Imaging of optically cleared brain tissue allowed quantification of the 3D microvascular characteristics in healthy mouse brains and in tissues with diffuse, infiltrative growing GBM8 brain tumors. Detailed 3D visualization revealed the organization of tumor cells relative to the vasculature, in both gray matter and white matter regions, and patterns of multicellular GBM networks collectively invading the brain parenchyma.

CONCLUSIONS

Optical tissue clearing opens new avenues for combined quantitative and 3D microscopic analysis of the topographical relationship between GBM cells and their microenvironment.

Quantifying three-dimensional rodent retina vascular development using optical tissue clearing and light-sheet microscopy

Retinal vasculature develops in a highly orchestrated three-dimensional (3-D) sequence. The stages of retinal vascularization are highly susceptible to oxygen perturbations. We demonstrate that optical tissue clearing of intact rat retinas and light-sheet microscopy provides rapid 3-D characterization of vascular complexity during retinal development. Compared with flat mount preparations that dissect the retina and primarily image the outermost vascular layers, intact cleared retinas imaged using light-sheet fluorescence microscopy display changes in the 3-D retinal vasculature rapidly without the need for point scanning techniques. Using a severe model of retinal vascular disruption, we demonstrate that a simple metric based on Sholl analysis captures the vascular changes observed during retinal development in 3-D. Taken together, these results provide a methodology for rapidly quantifying the 3-D development of the entire rodent retinal vasculature.

Smooth 2D manifold extraction from 3D image stack

Three-dimensional fluorescence microscopy followed by image processing is routinely used to study biological objects at various scales such as cells and tissue. However, maximum intensity projection, the most broadly used rendering tool, extracts a discontinuous layer of voxels, obliviously creating important artifacts and possibly misleading interpretation. Here we propose smooth manifold extraction, an algorithm that produces a continuous focused 2D extraction from a 3D volume, hence preserving local spatial relationships. We demonstrate the usefulness of our approach by applying it to various biological applications using confocal and wide-field microscopy 3D image stacks. We provide a parameter-free ImageJ/Fiji plugin that allows 2D visualization and interpretation of 3D image stacks with maximum accuracy.

Maximum Intensity Projection is a common tool to represent 3D biological imaging data in a 2D space, but it creates artefacts. Here the authors develop Smooth Manifold Extraction, an ImageJ/Fiji plugin, to preserve local spatial relationships when extracting the content of a 3D volume to a 2D space.

FluoRender: joint freehand segmentation and visualization for many-channel fluorescence data analysis

Background

Image segmentation and registration techniques have enabled biologists to place large amounts of volume data from fluorescence microscopy, morphed three-dimensionally, onto a common spatial frame. Existing tools built on volume visualization pipelines for single channel or red-green-blue (RGB) channels have become inadequate for the new challenges of fluorescence microscopy. For a three-dimensional atlas of the insect nervous system, hundreds of volume channels are rendered simultaneously, whereas fluorescence intensity values from each channel need to be preserved for versatile adjustment and analysis. Although several existing tools have incorporated support of multichannel data using various strategies, the lack of a flexible design has made true many-channel visualization and analysis unavailable. The most common practice for many-channel volume data presentation is still converting and rendering pseudosurfaces, which are inaccurate for both qualitative and quantitative evaluations.

Results

Here, we present an alternative design strategy that accommodates the visualization and analysis of about 100 volume channels, each of which can be interactively adjusted, selected, and segmented using freehand tools. Our multichannel visualization includes a multilevel streaming pipeline plus a triple-buffer compositing technique. Our method also preserves original fluorescence intensity values on graphics hardware, a crucial feature that allows graphics-processing-unit (GPU)-based processing for interactive data analysis, such as freehand segmentation. We have implemented the design strategies as a thorough restructuring of our original tool, FluoRender.

Conclusion

The redesign of FluoRender not only maintains the existing multichannel capabilities for a greatly extended number of volume channels, but also enables new analysis functions for many-channel data from emerging biomedical-imaging techniques.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1694-9) contains supplementary material, which is available to authorized users.

Fast assembling of neuron fragments in serial 3D sections

Reconstructing neurons from 3D image-stacks of serial sections of thick brain tissue is very time-consuming and often becomes a bottleneck in high-throughput brain mapping projects. We developed NeuronStitcher, a software suite for stitching non-overlapping neuron fragments reconstructed in serial 3D image sections. With its efficient algorithm and user-friendly interface, NeuronStitcher has been used successfully to reconstruct very large and complex human and mouse neurons.

Imaging the pharmacology of nanomaterials by intravital microscopy: Toward understanding their biological behavior

Therapeutic nanoparticles (NPs) can deliver cytotoxic chemotherapeutics and other drugs more safely and efficiently to patients; furthermore, selective delivery to target tissues can theoretically be accomplished actively through coating NPs with molecular ligands, and passively through exploiting physiological "enhanced permeability and retention" features. However, clinical trial results have been mixed in showing improved efficacy with drug nanoencapsulation, largely due to heterogeneous NP accumulation at target sites across patients. Thus, a clear need exists to better understand why many NP strategies fail in vivo and not result in significantly improved tumor uptake or therapeutic response. Multicolor in vivo confocal fluorescence imaging (intravital microscopy; IVM) enables integrated pharmacokinetic and pharmacodynamic (PK/PD) measurement at the single-cell level, and has helped answer key questions regarding the biological mechanisms of in vivo NP behavior. This review summarizes progress to date and also describes useful technical strategies for successful IVM experimentation.

Breathing control center neurons that promote arousal in mice

Slow, controlled breathing has been used for centuries to promote mental calming, and it is used clinically to suppress excessive arousal such as panic attacks. However, the physiological and neural basis of the relationship between breathing and higher-order brain activity is unknown. We found a neuronal subpopulation in the mouse preBötzinger complex (preBötC), the primary breathing rhythm generator, which regulates the balance between calm and arousal behaviors. Conditional, bilateral genetic ablation of the ~175 Cdh9/Dbx1 double-positive preBötC neurons in adult mice left breathing intact but increased calm behaviors and decreased time in aroused states. These neurons project to, synapse on, and positively regulate noradrenergic neurons in the locus coeruleus, a brain center implicated in attention, arousal, and panic that projects throughout the brain.

BFPTool: a software tool for analysis of Biomembrane Force Probe experiments

Background

The Biomembrane Force Probe is an approachable experimental technique commonly used for single-molecule force spectroscopy and experiments on biological interfaces. The technique operates in the range of forces from 0.1 pN to 1000 pN. Experiments are typically repeated many times, conditions are often not optimal, the captured video can be unstable and lose focus; this makes efficient analysis challenging, while out-of-the-box non-proprietary solutions are not freely available.

Results

This dedicated tool was developed to integrate and simplify the image processing and analysis of videomicroscopy recordings from BFP experiments. A novel processing feature, allowing the tracking of the pipette, was incorporated to address a limitation of preceding methods. Emphasis was placed on versatility and comprehensible user interface implemented in a graphical form.

Conclusions

An integrated analytical tool was implemented to provide a faster, simpler and more convenient way to process and analyse BFP experiments.

Electronic supplementary material

The online version of this article (doi:10.1186/s13628-016-0033-2) contains supplementary material, which is available to authorized users.

Optimized Protocol for Imaging Cleared Neural Tissues Using Light Microscopy

Understanding physical and chemical processes at an organismal scale is a fundamental goal in biology. While science is adept at explaining biological phenomena at both molecular and cellular levels, understanding how these processes translate to organismal functions remains a challenging problem. This issue is particularly significant for the nervous system where cell signaling and synaptic activities function in the context of broad neural networks. Recent progress in tissue clearing technologies lessens the barriers that previously prevented the study of large tissue samples while maintaining molecular and cellular resolution. While these new methods open vast opportunities and exciting new questions, the logistics of analyzing cellular processes in intact tissue have to be carefully considered. In this protocol, we outline a procedure to rapidly image intact brain tissue up to thousands of cubic millimeters. This experimental pipeline involves three steps: tissue clearing, tissue imaging, and data analysis. In an attempt to streamline the process for researchers entering this field, we address important considerations for each of these stages and describe an integrated solution to image intact biological tissues. Hopefully, this optimized protocol will lower the barrier of implementing high-resolution tissue imaging and facilitate the investigations of mesoscale questions at molecular and cellular resolution.

Fuzzy-Logic Based Detection and Characterization of Junctions and Terminations in Fluorescence Microscopy Images of Neurons

Digital reconstruction of neuronal cell morphology is an important step toward understanding the functionality of neuronal networks. Neurons are tree-like structures whose description depends critically on the junctions and terminations, collectively called critical points, making the correct localization and identification of these points a crucial task in the reconstruction process. Here we present a fully automatic method for the integrated detection and characterization of both types of critical points in fluorescence microscopy images of neurons. In view of the majority of our current studies, which are based on cultured neurons, we describe and evaluate the method for application to two-dimensional (2D) images. The method relies on directional filtering and angular profile analysis to extract essential features about the main streamlines at any location in an image, and employs fuzzy logic with carefully designed rules to reason about the feature values in order to make well-informed decisions about the presence of a critical point and its type. Experiments on simulated as well as real images of neurons demonstrate the detection performance of our method. A comparison with the output of two existing neuron reconstruction methods reveals that our method achieves substantially higher detection rates and could provide beneficial information to the reconstruction process.

Intrinsic Cornu Ammonis Area 1 Theta-Nested Gamma Oscillations Induced by Optogenetic Theta Frequency Stimulation

Gamma oscillations (30-120 Hz) are thought to be important for various cognitive functions, including perception and working memory, and disruption of these oscillations has been implicated in brain disorders, such as schizophrenia and Alzheimer's disease. The cornu ammonis area 1 (CA1) of the hippocampus receives gamma frequency inputs from upstream regions (cornu ammonis area 3 and medial entorhinal cortex) and generates itself a faster gamma oscillation. The exact nature and origin of the intrinsic CA1 gamma oscillation is still under debate. Here, we expressed channel rhodopsin-2 under the CaMKIIα promoter in mice and prepared hippocampal slices to produce a model of intrinsic CA1 gamma oscillations. Sinusoidal optical stimulation of CA1 at theta frequency was found to induce robust theta-nested gamma oscillations with a temporal and spatial profile similar to CA1 gamma in vivo The results suggest the presence of a single gamma rhythm generator with a frequency range of 65-75 Hz at 32 °C. Pharmacological analysis found that the oscillations depended on both AMPA and GABAA receptors. Cell-attached and whole-cell recordings revealed that excitatory neuron firing slightly preceded interneuron firing within each gamma cycle, suggesting that this intrinsic CA1 gamma oscillation is generated with a pyramidal-interneuron circuit mechanism.

SIGNIFICANCE STATEMENT

This study demonstrates that the cornu ammonis area 1 (CA1) is capable of generating intrinsic gamma oscillations in response to theta input. This gamma generator is independent of activity in the upstream regions, highlighting that CA1 can produce its own gamma oscillation in addition to inheriting activity from the upstream regions. This supports the theory that gamma oscillations predominantly function to achieve local synchrony, and that a local gamma generated in each area conducts the signal to the downstream region.

Identification, visualization and clonal analysis of intestinal stem cells in fish

Recently, a stochastic model of symmetrical stem cell division followed by neutral drift has been proposed for intestinal stem cells (ISCs), which has been suggested to represent the predominant mode of stem cell progression in mammals. In contrast, stem cells in the retina of teleost fish show an asymmetric division mode. To address whether the mode of stem cell division follows phylogenetic or ontogenetic routes, we analysed the entire gastrointestinal tract of the teleost medaka (Oryzias latipes). X-ray microcomputed tomography shows a correlation of 3D topography with the functional domains. Analysis of ISCs in proliferation assays and via genetically encoded lineage tracing highlights a stem cell niche in the furrow between the long intestinal folds that is functionally equivalent to mammalian intestinal crypts. Stem cells in this compartment are characterized by the expression of homologs of mammalian ISC markers – sox9, axin2 and lgr5 – emphasizing the evolutionary conservation of the Wnt pathway components in the stem cell niche of the intestine. The stochastic, sparse initial labelling of ISCs ultimately resulted in extended labelled or unlabelled domains originating from single stem cells in the furrow niche, contributing to both homeostasis and growth. Thus, different modes of stem cell division co-evolved within one organism, and in the absence of physical isolation in crypts, ISCs contribute to homeostatic growth.

Summary: Adult medaka intestinal stem cells (ISCs) proliferate within a niche functionally equivalent to that in the mammal. Like mammalian ISCs, but unlike medaka retinal stem cells, their mode of division is largely symmetric.

BlastNeuron for Automated Comparison, Retrieval and Clustering of 3D Neuron Morphologies

Characterizing the identity and types of neurons in the brain, as well as their associated function, requires a means of quantifying and comparing 3D neuron morphology. Presently, neuron comparison methods are based on statistics from neuronal morphology such as size and number of branches, which are not fully suitable for detecting local similarities and differences in the detailed structure. We developed BlastNeuron to compare neurons in terms of their global appearance, detailed arborization patterns, and topological similarity. BlastNeuron first compares and clusters 3D neuron reconstructions based on global morphology features and moment invariants, independent of their orientations, sizes, level of reconstruction and other variations. Subsequently, BlastNeuron performs local alignment between any pair of retrieved neurons via a tree-topology driven dynamic programming method. A 3D correspondence map can thus be generated at the resolution of single reconstruction nodes. We applied BlastNeuron to three datasets: (1) 10,000+ neuron reconstructions from a public morphology database, (2) 681 newly and manually reconstructed neurons, and (3) neurons reconstructions produced using several independent reconstruction methods. Our approach was able to accurately and efficiently retrieve morphologically and functionally similar neuron structures from large morphology database, identify the local common structures, and find clusters of neurons that share similarities in both morphology and molecular profiles.

Keto-polymethines: a versatile class of dyes with outstanding spectroscopic properties for in cellulo and in vivo two-photon microscopy imaging

The synthesis of keto-heptamethine derivatives has been expanded to various new symmetrical and asymmetrical structures, including an unprecedented di-anionic keto-polymethine. The spectroscopic behavior of these new dyes has been systematically and thoroughly investigated, revealing that the formation of hydrogen bond interactions with protic solvents is responsible for a dramatic enhancement of the fluorescence quantum yield in the far-red spectral region. The existence of these strong hydrogen-bond interactions was further confirmed by molecular dynamics simulations. These bis-dipolar polymethines exhibit large two-photon absorption (TPA) cross-sections (σ² in GM) in the near-infrared, making them ideal candidates for NIR-to-NIR two-photon microscopy imaging applications. We demonstrate that the molecular engineering of the hydrophilic/hydrophobic balance enables targeting of different cellular components, such as cytoplasm or cell membranes. Addition of appropriate substituents provides the molecule with high-water-solubility, affording efficient two-photon probes for angiography.

Emerging Imaging and Genomic Tools for Developmental Systems Biology

Animal development is a complex and dynamic process orchestrated by exquisitely timed cell lineage commitment, divisions, migration, and morphological changes at the single-cell level. In the past decade, extensive genetic, stem cell, and genomic studies provided crucial insights into molecular underpinnings and the functional importance of genetic pathways governing various cellular differentiation processes. However, it is still largely unknown how the precise coordination of these pathways is achieved at the whole-organism level and how the highly regulated spatiotemporal choreography of development is established in turn. Here, we discuss the latest technological advances in imaging and single-cell genomics that hold great promise for advancing our understanding of this intricate process. We propose an integrated approach that combines such methods to quantitatively decipher in vivo cellular dynamic behaviors and their underlying molecular mechanisms at the systems level with single-cell, single-molecule resolution.