Connectome-driven neural inventory of a complete visual system

Vision provides animals with detailed information about their surroundings, conveying diverse features such as color, form, and movement across the visual scene. Computing these parallel spatial features requires a large and diverse network of neurons, such that in animals as distant as flies and humans, visual regions comprise half the brain's volume. These visual brain regions often reveal remarkable structure-function relationships, with neurons organized along spatial maps with shapes that directly relate to their roles in visual processing. To unravel the stunning diversity of a complex visual system, a careful mapping of the neural architecture matched to tools for targeted exploration of that circuitry is essential. Here, we report a new connectome of the right optic lobe from a male Drosophila central nervous system FIB-SEM volume and a comprehensive inventory of the fly's visual neurons. We developed a computational framework to quantify the anatomy of visual neurons, establishing a basis for interpreting how their shapes relate to spatial vision. By integrating this analysis with connectivity information, neurotransmitter identity, and expert curation, we classified the ~53,000 neurons into 727 types, about half of which are systematically described and named for the first time. Finally, we share an extensive collection of split-GAL4 lines matched to our neuron type catalog. Together, this comprehensive set of tools and data unlock new possibilities for systematic investigations of vision in Drosophila, a foundation for a deeper understanding of sensory processing.

OME-Zarr: a cloud-optimized bioimaging file format with international community support

A growing community is constructing a next-generation file format (NGFF) for bioimaging to overcome problems of scalability and heterogeneity. Organized by the Open Microscopy Environment (OME), individuals and institutes across diverse modalities facing these problems have designed a format specification process (OME-NGFF) to address these needs. This paper brings together a wide range of those community members to describe the cloud-optimized format itself-OME-Zarr-along with tools and data resources available today to increase FAIR access and remove barriers in the scientific process. The current momentum offers an opportunity to unify a key component of the bioimaging domain-the file format that underlies so many personal, institutional, and global data management and analysis tasks.

A scalable implementation of the recursive least-squares algorithm for training spiking neural networks

Training spiking recurrent neural networks on neuronal recordings or behavioral tasks has become a popular way to study computations performed by the nervous system. As the size and complexity of neural recordings increase, there is a need for efficient algorithms that can train models in a short period of time using minimal resources. We present optimized CPU and GPU implementations of the recursive least-squares algorithm in spiking neural networks. The GPU implementation can train networks of one million neurons, with 100 million plastic synapses and a billion static synapses, about 1,000 times faster than an unoptimized reference CPU implementation. We demonstrate the code's utility by training a network, in less than an hour, to reproduce the activity of > 66, 000 recorded neurons of a mouse performing a decision-making task. The fast implementation enables a more interactive in-silico study of the dynamics and connectivity underlying multi-area computations. It also admits the possibility to train models as in-vivo experiments are being conducted, thus closing the loop between modeling and experiments.

Architecture and dynamics of a desmosome-endoplasmic reticulum complex

The endoplasmic reticulum (ER) forms a dynamic network that contacts other cellular membranes to regulate stress responses, calcium signalling and lipid transfer. Here, using high-resolution volume electron microscopy, we find that the ER forms a previously unknown association with keratin intermediate filaments and desmosomal cell-cell junctions. Peripheral ER assembles into mirror image-like arrangements at desmosomes and exhibits nanometre proximity to keratin filaments and the desmosome cytoplasmic plaque. ER tubules exhibit stable associations with desmosomes, and perturbation of desmosomes or keratin filaments alters ER organization, mobility and expression of ER stress transcripts. These findings indicate that desmosomes and the keratin cytoskeleton regulate the distribution, function and dynamics of the ER network. Overall, this study reveals a previously unknown subcellular architecture defined by the structural integration of ER tubules with an epithelial intercellular junction.

OME-Zarr: a cloud-optimized bioimaging file format with international community support

A growing community is constructing a next-generation file format (NGFF) for bioimaging to overcome problems of scalability and heterogeneity. Organized by the Open Microscopy Environment (OME), individuals and institutes across diverse modalities facing these problems have designed a format specification process (OME-NGFF) to address these needs. This paper brings together a wide range of those community members to describe the cloud-optimized format itself -- OME-Zarr -- along with tools and data resources available today to increase FAIR access and remove barriers in the scientific process. The current momentum offers an opportunity to unify a key component of the bioimaging domain -- the file format that underlies so many personal, institutional, and global data management and analysis tasks.

Automated reconstruction of whole-embryo cell lineages by learning from sparse annotations

We present a method to automatically identify and track nuclei in time-lapse microscopy recordings of entire developing embryos. The method combines deep learning and global optimization. On a mouse dataset, it reconstructs 75.8% of cell lineages spanning 1 h, as compared to 31.8% for the competing method. Our approach improves understanding of where and when cell fate decisions are made in developing embryos, tissues, and organs.

Long-term imaging reveals behavioral plasticity during C. elegans dauer exit

BACKGROUND

During their lifetime, animals must adapt their behavior to survive in changing environments. This ability requires the nervous system to undergo adjustments at distinct temporal scales, from short-term dynamic changes in expression of neurotransmitters and receptors to longer-term growth, spatial and connectivity reorganization, while integrating external stimuli. The nematode Caenorhabditis elegans provides a model of nervous system plasticity, in particular its dauer exit decision. Under unfavorable conditions, larvae will enter the non-feeding and non-reproductive stress-resistant dauer stage and adapt their behavior to cope with the harsh new environment, with active reversal under improved conditions leading to resumption of reproductive development. However, how different environmental stimuli regulate the exit decision mechanism and thereby drive the larva's behavioral change is unknown. To fill this gap and provide insights on behavioral changes over extended periods of time, we developed a new open hardware method for long-term imaging (12h) of C. elegans larvae.

RESULTS

Our WormObserver platform comprises open hardware and software components for video acquisition, automated processing of large image data (> 80k images/experiment) and data analysis. We identified dauer-specific behavioral motifs and characterized the behavioral trajectory of dauer exit in different environments and genetic backgrounds to identify key decision points and stimuli promoting dauer exit. Combining long-term behavioral imaging with transcriptomics data, we find that bacterial ingestion triggers a change in neuropeptide gene expression to establish post-dauer behavior.

CONCLUSIONS

Taken together, we show how a developing nervous system can robustly integrate environmental changes activate a developmental switch and adapt the organism's behavior to a new environment. WormObserver is generally applicable to other research questions within and beyond the C. elegans field, having a modular and customizable character and allowing assessment of behavioral plasticity over longer periods.

RS-FISH: precise, interactive, fast, and scalable FISH spot detection

Fluorescent in-situ hybridization (FISH)-based methods extract spatially resolved genetic and epigenetic information from biological samples by detecting fluorescent spots in microscopy images, an often challenging task. We present Radial Symmetry-FISH (RS-FISH), an accurate, fast, and user-friendly software for spot detection in two- and three-dimensional images. RS-FISH offers interactive parameter tuning and readily scales to large datasets and image volumes of cleared or expanded samples using distributed processing on workstations, clusters, or the cloud. RS-FISH maintains high detection accuracy and low localization error across a wide range of signal-to-noise ratios, a key feature for single-molecule FISH, spatial transcriptomics, or spatial genomics applications.

Cap-dependent translation initiation monitored in living cells

mRNA translation is tightly regulated to preserve cellular homeostasis. Despite extensive biochemical, genetic, and structural studies, a detailed understanding of mRNA translation regulation is lacking. Imaging methodologies able to resolve the binding dynamics of translation factors at single-cell and single-mRNA resolution were necessary to fully elucidate regulation of this paramount process. Here live-cell spectroscopy and single-particle tracking were combined to interrogate the binding dynamics of endogenous initiation factors to the 5'cap. The diffusion of initiation factors (IFs) changed markedly upon their association with mRNA. Quantifying their diffusion characteristics revealed the sequence of IFs assembly and disassembly in cell lines and the clustering of translation in neurons. This approach revealed translation regulation at high spatial and temporal resolution that can be applied to the formation of any endogenous complex that results in a measurable shift in diffusion.

Best practice standards for circular RNA research

Circular RNAs (circRNAs) are formed in all domains of life and via different mechanisms. There has been an explosion in the number of circRNA papers in recent years; however, as a relatively young field, circRNA biology has an urgent need for common experimental standards for isolating, analyzing, expressing and depleting circRNAs. Here we propose a set of guidelines for circRNA studies based on the authors' experience. This Perspective will specifically address the major class of circRNAs in Eukarya that are generated by a spliceosome-catalyzed back-splicing event. We hope that the implementation of best practice principles for circRNA research will help move the field forward and allow a better functional understanding of this fascinating group of RNAs.

The H4K20 demethylase DPY-21 regulates the dynamics of condensin DC binding

Condensin is a multi-subunit structural maintenance of chromosomes (SMC) complex that binds to and compacts chromosomes. Here, we addressed the regulation of condensin binding dynamics using Caenorhabditis elegans condensin DC, which represses X chromosomes in hermaphrodites for dosage compensation. We established fluorescence recovery after photobleaching (FRAP) using the SMC4 homolog DPY-27 and showed that a well-characterized ATPase mutation abolishes DPY-27 binding to X chromosomes. Next, we performed FRAP in the background of several chromatin modifier mutants that cause varying degrees of X chromosome derepression. The greatest effect was in a null mutant of the H4K20me2 demethylase DPY-21, where the mobile fraction of condensin DC reduced from ∼30% to 10%. In contrast, a catalytic mutant of dpy-21 did not regulate condensin DC mobility. Hi-C sequencing data from the dpy-21 null mutant showed little change compared to wild-type data, uncoupling Hi-C-measured long-range DNA contacts from transcriptional repression of the X chromosomes. Taken together, our results indicate that DPY-21 has a non-catalytic role in regulating the dynamics of condensin DC binding, which is important for transcription repression.

Summary: The histone demethylase DPY-21 has catalytic and non-catalytic roles in condensin DC-mediated X chromosome repression. The non-catalytic activity regulates dynamics of condensin DC binding to X chromosomes.

Dual-view light-sheet imaging through a tilted glass interface using a deformable mirror

Light-sheet microscopy has become indispensable for imaging developing organisms, and imaging from multiple directions (views) is essential to improve its spatial resolution. We combine multi-view light-sheet microscopy with microfluidics using adaptive optics (deformable mirror) which corrects aberrations introduced by the 45^o-tilted glass coverslip. The optimal shape of the deformable mirror is computed by an iterative algorithm that optimizes the point-spread function in two orthogonal views. Simultaneous correction in two optical arms is achieved via a knife-edge mirror that splits the excitation path and combines the detection paths. Our design allows multi-view light-sheet microscopy with microfluidic devices for precisely controlled experiments and high-content screening.

FRC-QE: a robust and comparable 3D microscopy image quality metric for cleared organoids

SUMMARY

Here, we propose Fourier ring correlation-based quality estimation (FRC-QE) as a new metric for automated image quality estimation in 3D fluorescence microscopy acquisitions of cleared organoids that yields comparable measurements across experimental replicates, clearing protocols and works for different microscopy modalities.

AVAILABILITY AND IMPLEMENTATION

FRC-QE is written in ImgLib2/Java and provided as an easy-to-use and macro-scriptable plugin for Fiji. Code, documentation, sample images and further information can be found under https://github.com/PreibischLab/FRC-QE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Regulatory encoding of quantitative variation in spatial activity of a Drosophila enhancer

Developmental enhancers control the expression of genes prefiguring morphological patterns. The activity of an enhancer varies among cells of a tissue, but collectively, expression levels in individual cells constitute a spatial pattern of gene expression. How the spatial and quantitative regulatory information is encoded in an enhancer sequence is elusive. To link spatial pattern and activity levels of an enhancer, we used systematic mutations of the yellow spot enhancer, active in developing Drosophila wings, and tested their effect in a reporter assay. Moreover, we developed an analytic framework based on the comprehensive quantification of spatial reporter activity. We show that the quantitative enhancer activity results from densely packed regulatory information along the sequence, and that a complex interplay between activators and multiple tiers of repressors carves the spatial pattern. Our results shed light on how an enhancer reads and integrates trans-regulatory landscape information to encode a spatial quantitative pattern.

Stochastic transcription in the p53-mediated response to DNA damage is modulated by burst frequency

Discontinuous transcription has been described for different mammalian cell lines and numerous promoters. However, our knowledge of how the activity of individual promoters is adjusted by dynamic signaling inputs from transcription factors is limited. To address this question, we characterized the activity of selected target genes that are regulated by pulsatile accumulation of the tumor suppressor p53 in response to ionizing radiation. We performed time-resolved measurements of gene expression at the single-cell level by smFISH and used the resulting data to inform a mathematical model of promoter activity. We found that p53 target promoters are regulated by frequency modulation of stochastic bursting and can be grouped along three archetypes of gene expression. The occurrence of these archetypes cannot solely be explained by nuclear p53 abundance or promoter binding of total p53. Instead, we provide evidence that the time-varying acetylation state of p53's C-terminal lysine residues is critical for gene-specific regulation of stochastic bursting.

Oscillations of MyoD and Hes1 proteins regulate the maintenance of activated muscle stem cells

Lahmann et al. show that Hes1 controls the balance between proliferation and differentiation of activated muscle stem cells in both developing and regenerating muscle. Hes1 is expressed in an oscillatory manner in activated stem cells, where it drives the oscillatory expression of MyoD.

The balance between proliferation and differentiation of muscle stem cells is tightly controlled, ensuring the maintenance of a cellular pool needed for muscle growth and repair. We demonstrate here that the transcriptional regulator Hes1 controls the balance between proliferation and differentiation of activated muscle stem cells in both developing and regenerating muscle. We observed that Hes1 is expressed in an oscillatory manner in activated stem cells where it drives the oscillatory expression of MyoD. MyoD expression oscillates in activated muscle stem cells from postnatal and adult muscle under various conditions: when the stem cells are dispersed in culture, when they remain associated with single muscle fibers, or when they reside in muscle biopsies. Unstable MyoD oscillations and long periods of sustained MyoD expression are observed in differentiating cells. Ablation of the Hes1 oscillator in stem cells interfered with stable MyoD oscillations and led to prolonged periods of sustained MyoD expression, resulting in increased differentiation propensity. This interfered with the maintenance of activated muscle stem cells, and impaired muscle growth and repair. We conclude that oscillatory MyoD expression allows the cells to remain in an undifferentiated and proliferative state and is required for amplification of the activated stem cell pool.

Multi-view light-sheet imaging and tracking with the MaMuT software reveals the cell lineage of a direct developing arthropod limb

During development, coordinated cell behaviors orchestrate tissue and organ morphogenesis. Detailed descriptions of cell lineages and behaviors provide a powerful framework to elucidate the mechanisms of morphogenesis. To study the cellular basis of limb development, we imaged transgenic fluorescently-labeled embryos from the crustacean Parhyale hawaiensis with multi-view light-sheet microscopy at high spatiotemporal resolution over several days of embryogenesis. The cell lineage of outgrowing thoracic limbs was reconstructed at single-cell resolution with new software called Massive Multi-view Tracker (MaMuT). In silico clonal analyses suggested that the early limb primordium becomes subdivided into anterior-posterior and dorsal-ventral compartments whose boundaries intersect at the distal tip of the growing limb. Limb-bud formation is associated with spatial modulation of cell proliferation, while limb elongation is also driven by preferential orientation of cell divisions along the proximal-distal growth axis. Cellular reconstructions were predictive of the expression patterns of limb development genes including the BMP morphogen Decapentaplegic.

During early life, animals develop from a single fertilized egg cell to hundreds, millions or even trillions of cells. These cells specialize to do different tasks; forming different tissues and organs like muscle, skin, lungs and liver. For more than a century, scientists have strived to understand the details of how animal cells become different and specialize, and have created many new techniques and technologies to help them achieve this goal.

Limbs – such as arms, legs and wings – form from small lumps of cells called limb buds. Scientists use the shrimp-like crustacean, Parhyale hawaiensis, to study development, including limb growth. This species is useful because it is easy to grow, manipulate and observe its developing young in the laboratory. Understanding how its limbs develop offers important new insights into how limbs develop in other animals too. Wolff, Tinevez, Pietzsch et al. have now combined advanced microscopy with custom computer software, called Massive Multi-view Tracker (MaMuT) to investigate this.

As limbs develop in Parhyale, the MaMuT software tracks how cells behave, and how they are organized. This analysis revealed that for cells to produce a limb bud, they need to split at an early stage into separate groups. These groups are organized along two body axes, one that goes from head to tail, and one that runs from back to belly. The limb grows perpendicular to these main body axes, along a new ‘proximal-distal’ axis that goes from nearest to furthest from the body. Wolff et al. found that the cells that contribute to the extremities of the limb divide faster than the ones that stay closer to the body. Finally, the results show that when cells in a limb divide, they mostly divide along the proximal-distal axis, producing one cell that is further from the body than the other. These cell activities may help limbs to get longer as they grow.

Notably, the groups of cells seen by Wolff et al. were expressing genes that had previously been identified in developing limbs. This helps to validate the new results and to identify which active genes control the behaviors of the analyzed cells.

These findings reveal new ways to study animal development. This approach could have many research uses and may help to link the mechanisms of cell biology to their effects. It could also contribute to new understanding of developmental and genetic conditions that affect human health.

Using Light Sheet Fluorescence Microscopy to Image Zebrafish Eye Development

Light sheet fluorescence microscopy (LSFM) is gaining more and more popularity as a method to image embryonic development. The main advantages of LSFM compared to confocal systems are its low phototoxicity, gentle mounting strategies, fast acquisition with high signal to noise ratio and the possibility of imaging samples from various angles (views) for long periods of time. Imaging from multiple views unleashes the full potential of LSFM, but at the same time it can create terabyte-sized datasets. Processing such datasets is the biggest challenge of using LSFM. In this protocol we outline some solutions to this problem. Until recently, LSFM was mostly performed in laboratories that had the expertise to build and operate their own light sheet microscopes. However, in the last three years several commercial implementations of LSFM became available, which are multipurpose and easy to use for any developmental biologist. This article is primarily directed to those researchers, who are not LSFM technology developers, but want to employ LSFM as a tool to answer specific developmental biology questions. Here, we use imaging of zebrafish eye development as an example to introduce the reader to LSFM technology and we demonstrate applications of LSFM across multiple spatial and temporal scales. This article describes a complete experimental protocol starting with the mounting of zebrafish embryos for LSFM. We then outline the options for imaging using the commercially available light sheet microscope. Importantly, we also explain a pipeline for subsequent registration and fusion of multiview datasets using an open source solution implemented as a Fiji plugin. While this protocol focuses on imaging the developing zebrafish eye and processing data from a particular imaging setup, most of the insights and troubleshooting suggestions presented here are of general use and the protocol can be adapted to a variety of light sheet microscopy experiments.

An automated workflow for parallel processing of large multiview SPIM recordings

Summary: Selective Plane Illumination Microscopy (SPIM) allows to image developing organisms in 3D at unprecedented temporal resolution over long periods of time. The resulting massive amounts of raw image data requires extensive processing interactively via dedicated graphical user interface (GUI) applications. The consecutive processing steps can be easily automated and the individual time points can be processed independently, which lends itself to trivial parallelization on a high performance computing (HPC) cluster. Here, we introduce an automated workflow for processing large multiview, multichannel, multiillumination time-lapse SPIM data on a single workstation or in parallel on a HPC cluster. The pipeline relies on snakemake to resolve dependencies among consecutive processing steps and can be easily adapted to any cluster environment for processing SPIM data in a fraction of the time required to collect it.

Availability and implementation: The code is distributed free and open source under the MIT license http://opensource.org/licenses/MIT. The source code can be downloaded from github: https://github.com/mpicbg-scicomp/snakemake-workflows. Documentation can be found here: http://fiji.sc/Automated_workflow_for_parallel_Multiview_Reconstruction.

Contact: schmied@mpi-cbg.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Nuclear accessibility of β-actin mRNA is measured by 3D single-molecule real-time tracking

Imaging single proteins or RNAs allows direct visualization of the inner workings of the cell. Typically, three-dimensional (3D) images are acquired by sequentially capturing a series of 2D sections. The time required to step through the sample often impedes imaging of large numbers of rapidly moving molecules. Here we applied multifocus microscopy (MFM) to instantaneously capture 3D single-molecule real-time images in live cells, visualizing cell nuclei at 10 volumes per second. We developed image analysis techniques to analyze messenger RNA (mRNA) diffusion in the entire volume of the nucleus. Combining MFM with precise registration between fluorescently labeled mRNA, nuclear pore complexes, and chromatin, we obtained globally optimal image alignment within 80-nm precision using transformation models. We show that β-actin mRNAs freely access the entire nucleus and fewer than 60% of mRNAs are more than 0.5 µm away from a nuclear pore, and we do so for the first time accounting for spatial inhomogeneity of nuclear organization.

Efficient Bayesian-based multiview deconvolution

Light-sheet fluorescence microscopy is able to image large specimens with high resolution by capturing the samples from multiple angles. Multiview deconvolution can substantially improve the resolution and contrast of the images, but its application has been limited owing to the large size of the data sets. Here we present a Bayesian-based derivation of multiview deconvolution that drastically improves the convergence time, and we provide a fast implementation using graphics hardware.

Bioimage Informatics in the context of Drosophila research

Modern biological research relies heavily on microscopic imaging. The advanced genetic toolkit of Drosophila makes it possible to label molecular and cellular components with unprecedented level of specificity necessitating the application of the most sophisticated imaging technologies. Imaging in Drosophila spans all scales from single molecules to the entire populations of adult organisms, from electron microscopy to live imaging of developmental processes. As the imaging approaches become more complex and ambitious, there is an increasing need for quantitative, computer-mediated image processing and analysis to make sense of the imagery. Bioimage Informatics is an emerging research field that covers all aspects of biological image analysis from data handling, through processing, to quantitative measurements, analysis and data presentation. Some of the most advanced, large scale projects, combining cutting edge imaging with complex bioimage informatics pipelines, are realized in the Drosophila research community. In this review, we discuss the current research in biological image analysis specifically relevant to the type of systems level image datasets that are uniquely available for the Drosophila model system. We focus on how state-of-the-art computer vision algorithms are impacting the ability of Drosophila researchers to analyze biological systems in space and time. We pay particular attention to how these algorithmic advances from computer science are made usable to practicing biologists through open source platforms and how biologists can themselves participate in their further development.

ImgLib2--generic image processing in Java

Summary: ImgLib2 is an open-source Java library for n-dimensional data representation and manipulation with focus on image processing. It aims at minimizing code duplication by cleanly separating pixel-algebra, data access and data representation in memory. Algorithms can be implemented for classes of pixel types and generic access patterns by which they become independent of the specific dimensionality, pixel type and data representation. ImgLib2 illustrates that an elegant high-level programming interface can be achieved without sacrificing performance. It provides efficient implementations of common data types, storage layouts and algorithms. It is the data model underlying ImageJ2, the KNIME Image Processing toolbox and an increasing number of Fiji-Plugins.

Availability: ImgLib2 is licensed under BSD. Documentation and source code are available at http://imglib2.net and in a public repository at https://github.com/imagej/imglib.

Supplementary Information:Supplementary data are available at Bioinformatics Online.

Contact: saalfeld@mpi-cbg.de

Fiji: an open-source platform for biological-image analysis

Fiji is a distribution of the popular open-source software ImageJ focused on biological-image analysis. Fiji uses modern software engineering practices to combine powerful software libraries with a broad range of scripting languages to enable rapid prototyping of image-processing algorithms. Fiji facilitates the transformation of new algorithms into ImageJ plugins that can be shared with end users through an integrated update system. We propose Fiji as a platform for productive collaboration between computer science and biology research communities.

Fiji: an open-source platform for biological-image analysis

Fiji is a distribution of the popular open-source software ImageJ focused on biological-image analysis. Fiji uses modern software engineering practices to combine powerful software libraries with a broad range of scripting languages to enable rapid prototyping of image-processing algorithms. Fiji facilitates the transformation of new algorithms into ImageJ plugins that can be shared with end users through an integrated update system. We propose Fiji as a platform for productive collaboration between computer science and biology research communities.

Fiji: an open-source platform for biological-image analysis

Fiji is a distribution of the popular open-source software ImageJ focused on biological-image analysis. Fiji uses modern software engineering practices to combine powerful software libraries with a broad range of scripting languages to enable rapid prototyping of image-processing algorithms. Fiji facilitates the transformation of new algorithms into ImageJ plugins that can be shared with end users through an integrated update system. We propose Fiji as a platform for productive collaboration between computer science and biology research communities.

TrakEM2 software for neural circuit reconstruction

A key challenge in neuroscience is the expeditious reconstruction of neuronal circuits. For model systems such as Drosophila and C. elegans, the limiting step is no longer the acquisition of imagery but the extraction of the circuit from images. For this purpose, we designed a software application, TrakEM2, that addresses the systematic reconstruction of neuronal circuits from large electron microscopical and optical image volumes. We address the challenges of image volume composition from individual, deformed images; of the reconstruction of neuronal arbors and annotation of synapses with fast manual and semi-automatic methods; and the management of large collections of both images and annotations. The output is a neural circuit of 3d arbors and synapses, encoded in NeuroML and other formats, ready for analysis.

Restoration of uneven illumination in light sheet microscopy images

Light microscopy images suffer from poor contrast due to light absorption and scattering by the media. The resulting decay in contrast varies exponentially across the image along the incident light path. Classical space invariant deconvolution approaches, while very effective in deblurring, are not designed for the restoration of uneven illumination in microscopy images. In this article, we present a modified radiative transfer theory approach to solve the contrast degradation problem of light sheet microscopy (LSM) images. We confirmed the effectiveness of our approach through simulation as well as real LSM images.

Global analysis of nascent RNA reveals transcriptional pausing in terminal exons

Pre-mRNA splicing is catalyzed by the spliceosome, which can assemble on pre-mRNA cotranscriptionally. However, whether splicing generally occurs during transcription has not been addressed. Indeed, splicing catalysis is expected to occur posttranscriptionally in yeast, where the shortness of terminal exons should leave insufficient time for splicing. Here, we isolate endogenous S. cerevisiae nascent RNA and determine gene-specific splicing efficiencies and transcription profiles, using high-density tiling microarrays. Surprisingly, we find that splicing occurs cotranscriptionally for the majority of intron-containing genes. Analysis of transcription profiles reveals Pol II pausing within the terminal exons of these genes. Intronless and inefficiently spliced genes lack this pause. In silico simulations of transcription and splicing kinetics confirm that this pausing event provides sufficient time for splicing before termination. The discovery of terminal exon pausing demonstrates functional coupling of transcription and splicing near gene ends.

An integrated micro- and macroarchitectural analysis of the Drosophila brain by computer-assisted serial section electron microscopy

The analysis of microcircuitry (the connectivity at the level of individual neuronal processes and synapses), which is indispensable for our understanding of brain function, is based on serial transmission electron microscopy (TEM) or one of its modern variants. Due to technical limitations, most previous studies that used serial TEM recorded relatively small stacks of individual neurons. As a result, our knowledge of microcircuitry in any nervous system is very limited. We applied the software package TrakEM2 to reconstruct neuronal microcircuitry from TEM sections of a small brain, the early larval brain of Drosophila melanogaster. TrakEM2 enables us to embed the analysis of the TEM image volumes at the microcircuit level into a light microscopically derived neuro-anatomical framework, by registering confocal stacks containing sparsely labeled neural structures with the TEM image volume. We imaged two sets of serial TEM sections of the Drosophila first instar larval brain neuropile and one ventral nerve cord segment, and here report our first results pertaining to Drosophila brain microcircuitry. Terminal neurites fall into a small number of generic classes termed globular, varicose, axiform, and dendritiform. Globular and varicose neurites have large diameter segments that carry almost exclusively presynaptic sites. Dendritiform neurites are thin, highly branched processes that are almost exclusively postsynaptic. Due to the high branching density of dendritiform fibers and the fact that synapses are polyadic, neurites are highly interconnected even within small neuropile volumes. We describe the network motifs most frequently encountered in the Drosophila neuropile. Our study introduces an approach towards a comprehensive anatomical reconstruction of neuronal microcircuitry and delivers microcircuitry comparisons between vertebrate and insect neuropile.

Planar cell polarity signalling regulates cell adhesion properties in progenitors of the zebrafish laterality organ

Organ formation requires the precise assembly of progenitor cells into a functional multicellular structure. Mechanical forces probably participate in this process but how they influence organ morphogenesis is still unclear. Here, we show that Wnt11- and Prickle1a-mediated planar cell polarity (PCP) signalling coordinates the formation of the zebrafish ciliated laterality organ (Kupffer's vesicle) by regulating adhesion properties between organ progenitor cells (the dorsal forerunner cells, DFCs). Combined inhibition of Wnt11 and Prickle1a reduces DFC cell-cell adhesion and impairs their compaction and arrangement during vesicle lumen formation. This leads to the formation of a mis-shapen vesicle with small fragmented lumina and shortened cilia, resulting in severely impaired organ function and, as a consequence, randomised laterality of both molecular and visceral asymmetries. Our results reveal a novel role for PCP-dependent cell adhesion in coordinating the supracellular organisation of progenitor cells during vertebrate laterality organ formation.

Globally optimal stitching of tiled 3D microscopic image acquisitions

Motivation: Modern anatomical and developmental studies often require high-resolution imaging of large specimens in three dimensions (3D). Confocal microscopy produces high-resolution 3D images, but is limited by a relatively small field of view compared with the size of large biological specimens. Therefore, motorized stages that move the sample are used to create a tiled scan of the whole specimen. The physical coordinates provided by the microscope stage are not precise enough to allow direct reconstruction (Stitching) of the whole image from individual image stacks.

Results: To optimally stitch a large collection of 3D confocal images, we developed a method that, based on the Fourier Shift Theorem, computes all possible translations between pairs of 3D images, yielding the best overlap in terms of the cross-correlation measure and subsequently finds the globally optimal configuration of the whole group of 3D images. This method avoids the propagation of errors by consecutive registration steps. Additionally, to compensate the brightness differences between tiles, we apply a smooth, non-linear intensity transition between the overlapping images. Our stitching approach is fast, works on 2D and 3D images, and for small image sets does not require prior knowledge about the tile configuration.

Availability: The implementation of this method is available as an ImageJ plugin distributed as a part of the Fiji project (FijiisjustImageJ: http://pacific.mpi-cbg.de/).

Contact:tomancak@mpi-cbg.de

"Hook"-calibration of GeneChip-microarrays: chip characteristics and expression measures

Background

Microarray experiments rely on several critical steps that may introduce biases and uncertainty in downstream analyses. These steps include mRNA sample extraction, amplification and labelling, hybridization, and scanning causing chip-specific systematic variations on the raw intensity level. Also the chosen array-type and the up-to-dateness of the genomic information probed on the chip affect the quality of the expression measures. In the accompanying publication we presented theory and algorithm of the so-called hook method which aims at correcting expression data for systematic biases using a series of new chip characteristics.

Results

In this publication we summarize the essential chip characteristics provided by this method, analyze special benchmark experiments to estimate transcript related expression measures and illustrate the potency of the method to detect and to quantify the quality of a particular hybridization. It is shown that our single-chip approach provides expression measures responding linearly on changes of the transcript concentration over three orders of magnitude. In addition, the method calculates a detection call judging the relation between the signal and the detection limit of the particular measurement. The performance of the method in the context of different chip generations and probe set assignments is illustrated. The hook method characterizes the RNA-quality in terms of the 3'/5'-amplification bias and the sample-specific calling rate. We show that the proper judgement of these effects requires the disentanglement of non-specific and specific hybridization which, otherwise, can lead to misinterpretations of expression changes. The consequences of modifying probe/target interactions by either changing the labelling protocol or by substituting RNA by DNA targets are demonstrated.

Conclusion

The single-chip based hook-method provides accurate expression estimates and chip-summary characteristics using the natural metrics given by the hybridization reaction with the potency to develop new standards for microarray quality control and calibration.

"Hook"-calibration of GeneChip-microarrays: theory and algorithm

BACKGROUND

: The improvement of microarray calibration methods is an essential prerequisite for quantitative expression analysis. This issue requires the formulation of an appropriate model describing the basic relationship between the probe intensity and the specific transcript concentration in a complex environment of competing interactions, the estimation of the magnitude these effects and their correction using the intensity information of a given chip and, finally the development of practicable algorithms which judge the quality of a particular hybridization and estimate the expression degree from the intensity values.

RESULTS

: We present the so-called hook-calibration method which co-processes the log-difference (delta) and -sum (sigma) of the perfect match (PM) and mismatch (MM) probe-intensities. The MM probes are utilized as an internal reference which is subjected to the same hybridization law as the PM, however with modified characteristics. After sequence-specific affinity correction the method fits the Langmuir-adsorption model to the smoothed delta-versus-sigma plot. The geometrical dimensions of this so-called hook-curve characterize the particular hybridization in terms of simple geometric parameters which provide information about the mean non-specific background intensity, the saturation value, the mean PM/MM-sensitivity gain and the fraction of absent probes. This graphical summary spans a metrics system for expression estimates in natural units such as the mean binding constants and the occupancy of the probe spots. The method is single-chip based, i.e. it separately uses the intensities for each selected chip.

CONCLUSION

: The hook-method corrects the raw intensities for the non-specific background hybridization in a sequence-specific manner, for the potential saturation of the probe-spots with bound transcripts and for the sequence-specific binding of specific transcripts. The obtained chip characteristics in combination with the sensitivity corrected probe-intensity values provide expression estimates scaled in natural units which are given by the binding constants of the particular hybridization.

Base pair interactions and hybridization isotherms of matched and mismatched oligonucleotide probes on microarrays

The microarray technology enables the expression degree of thousands of genes to be estimated at once by the measurement of the abundance of the respective messenger RNA. This method is based on the sequence specific binding of RNA to DNA probes and its detection using fluorescent labels. The raw intensity data are affected by the sequence-specific affinity of probe and RNA for duplex formation, by the background intensity due to nonspecific hybridization at small transcript concentrations and by the saturation of the probes at high transcript concentration owing to surface adsorption. We address these issues using a binding model which describes specific and nonspecific hybridization in terms of a competitive two-species Langmuir isotherm and DNA/RNA duplex formation in terms of sequence-specific, single-base related interactions. The GeneChip microarrays technology uses pairs of so-called perfect match (PM) and mismatch (MM) oligonucleotide probes to estimate the amount of nonspecific hybridization. The mean affinity of the probes decrease according to PM(specific) > MM(specific) >> PM(nonspecific) approximately MM(nonspecific). The stability of specific and nonspecific DNA/RNA duplexes is mainly determined by Watson Crick (WC) pairings. Mismatched self-complementary pairings in the middle of the MM sequence only weakly contribute to the duplex stability. The asymmetry of base pair interaction in the DNA/RNA hybrid duplexes gives rise to a duplet-like symmetry of the PM - MM intensity difference at dominating nonspecific hybridization and a triplet-like symmetry at specific hybridization. The signal intensities of the PM and MM probes and their difference are assessed in terms of sensitivity and specificity. The presented results imply the refinement of existing algorithms of probe level analysis to correct microarray data for nonspecific background intensities and saturation on the basis of the probe sequence.

Specific and nonspecific hybridization of oligonucleotide probes on microarrays

Gene expression analysis by means of microarrays is based on the sequence-specific binding of RNA to DNA oligonucleotide probes and its measurement using fluorescent labels. The binding of RNA fragments involving sequences other than the intended target is problematic because it adds a chemical background to the signal, which is not related to the expression degree of the target gene. The article presents a molecular signature of specific and nonspecific hybridization with potential consequences for gene expression analysis. We analyzed the signal intensities of perfect match (PM) and mismatch (MM) probes of GeneChip microarrays to specify the effect of specific and nonspecific hybridization. We found that these events give rise to different relations between the PM and MM intensities as function of the middle base of the PM, namely a triplet-like (C > G approximately T > A > 0) and a duplet-like (C approximately T > 0 > G approximately A) pattern of the PM-MM log-intensity difference upon binding of specific and nonspecific RNA fragments, respectively. The systematic behavior of the intensity difference can be rationalized on the level of basepairings of DNA/RNA oligonucleotide duplexes in the middle of the probe sequence. Nonspecific binding is characterized by the reversal of the central Watson-Crick (WC) pairing for each PM/MM probe pair, whereas specific binding refers to the combination of a WC and a self-complementary (SC) pairing in PM and MM probes, respectively. The Gibbs free energy contribution of WC pairs to duplex stability is asymmetric for purines and pyrimidines of the PM and decreases according to C > G approximately T > A. SC pairings on the average only weakly contribute to duplex stability. The intensity of complementary MM introduces a systematic source of variation which decreases the precision of expression measures based on the MM intensities.