Automated deep-phenotyping of the vertebrate brain

Here, we describe an automated platform suitable for large-scale deep-phenotyping of zebrafish mutant lines, which uses optical projection tomography to rapidly image brain-specific gene expression patterns in 3D at cellular resolution. Registration algorithms and correlation analysis are then used to compare 3D expression patterns, to automatically detect all statistically significant alterations in mutants, and to map them onto a brain atlas. Automated deep-phenotyping of a mutation in the master transcriptional regulator fezf2 not only detects all known phenotypes but also uncovers important novel neural deficits that were overlooked in previous studies. In the telencephalon, we show for the first time that fezf2 mutant zebrafish have significant patterning deficits, particularly in glutamatergic populations. Our findings reveal unexpected parallels between fezf2 function in zebrafish and mice, where mutations cause deficits in glutamatergic neurons of the telencephalon-derived neocortex.National Institutes of Health (U.S.) (Director’s Pioneer Award DP1-NS082101)David & Lucile Packard Foundation. Award in Science and EngineeringBroad Institute of MIT and Harvard (SPARC Award)Epilepsy Foundation of America (Postdoctoral Fellowship

Visualizing anatomically registered data with Brainrender

Three-dimensional (3D) digital brain atlases and high-throughput brain wide imaging techniques generate large multidimensional datasets that can be registered to a common reference frame. Generating insights from such datasets depends critically on visualization and interactive data exploration, but this a challenging task. Currently available software is dedicated to single atlases, model species or data types, and generating 3D renderings that merge anatomically registered data from diverse sources requires extensive development and programming skills. Here, we present brainrender: an open-source Python package for interactive visualization of multidimensional datasets registered to brain atlases. Brainrender facilitates the creation of complex renderings with different data types in the same visualization and enables seamless use of different atlas sources. High-quality visualizations can be used interactively and exported as high-resolution figures and animated videos. By facilitating the visualization of anatomically registered data, brainrender should accelerate the analysis, interpretation, and dissemination of brain-wide multidimensional data

Doctor of Philosophy

dissertationConfocal microscopy has become a popular imaging technique in biology research in recent years. It is often used to study three-dimensional (3D) structures of biological samples. Confocal data are commonly multichannel, with each channel resulting from a different fluorescent staining. This technique also results in finely detailed structures in 3D, such as neuron fibers. Despite the plethora of volume rendering techniques that have been available for many years, there is a demand from biologists for a flexible tool that allows interactive visualization and analysis of multichannel confocal data. Together with biologists, we have designed and developed FluoRender. It incorporates volume rendering techniques such as a two-dimensional (2D) transfer function and multichannel intermixing. Rendering results can be enhanced through tone-mappings and overlays. To facilitate analyses of confocal data, FluoRender provides interactive operations for extracting complex structures. Furthermore, we developed the Synthetic Brainbow technique, which takes advantage of the asynchronous behavior in Graphics Processing Unit (GPU) framebuffer loops and generates random colorizations for different structures in single-channel confocal data. The results from our Synthetic Brainbows, when applied to a sequence of developing cells, can then be used for tracking the movements of these cells. Finally, we present an application of FluoRender in the workflow of constructing anatomical atlases

Abnormal Brain Connectivity in the Primary Visual Pathway in Human Albinism

In albinism, the ipsilateral projection of retinal axons is significantly reduced, and most fibres project contralaterally. The retina and optic chiasm have been proposed as sites for misrouting. The number of lateral geniculate nucleus (LGN) relay neurons has been linked to LGN volume, suggesting a correlation between LGN size and the number of tracts traveling through the optic radiation (OR) to the primary visual cortex (V1). Using diffusion data and both deterministic and probabilistic tractography, we studied differences in OR between albinism and controls. Statistical analyses measured white matter integrity in areas corresponding to the OR, as well as LGN to V1 connectivity. Results revealed reduced white matter integrity and connectivity in the OR region in albinism compared to controls, suggesting altered structural development. Previous reports of smaller LGN and the altered thalamo-cortical connectivity reported here demonstrate the effect of misrouting on structural organization of the visual pathway in albinism

Geometric representation of neuroanatomical data observed in mouse brain at cellular and gross levels

This dissertation studies two problems related to geometric representation of neuroanatomical data: (i) spatial representation and organization of individual neurons, and (ii) reconstruction of three-dimensional neuroanatomical regions from sparse two-dimensional drawings. This work has been motivated by nearby development of new technology, Knife-Edge Scanning Microscopy (KESM), that images a whole mouse brain at cellular level in less than a month. A method is introduced to represent neuronal data observed in the mammalian brain at the cellular level using geometric primitives and spatial indexing. A data representation scheme is defined that captures the geometry of individual neurons using traditional geometric primitives, points and cross-sectional areas along a trajectory. This representation captures inferred synapses as directed links between primitives and spatially indexes observed neurons based on the locations of their cell bodies. This method provides a set of rules for acquisition, representation, and indexing of KESMgenerated data. Neuroanatomical data observed at the gross level provides the underlying regional framework for neuronal circuits. Accumulated expert knowledge on neuroanatomical organization is usually given as a series of sparse two-dimensional contours. A data structure and an algorithm are described to reconstruct separating surfaces among multiple regions from these sparse cross-sectional contours. A topology graph is defined for each region that describes the topological skeleton of the region’s boundary surface and that shows between which contours the surface patches should be generated. A graph-directed triangulation algorithm is provided to reconstruct surface patches between contours. This graph-directed triangulation algorithm combined together with a piecewise parametric curve fitting technique ensures that abutting or shared surface patches are precisely coincident. This method overcomes limitations in i) traditional surfaces-from-contours algorithms that assume binary, not multiple, regionalization of space, and in ii) few existing separating surfaces algorithms that assume conversion of input into a regular volumetric grid, which is not possible with sparse inter-planar resolution

A Pipeline for Automated Assessment of Cell Location in 3D Mouse Brain Image Sets

Mapping the neuronal connectivity of the mouse brain has long been hampered by the laborious and time-consuming process of slicing, staining and imaging the brain tissue. Recent developments in automated 3D fluorescence microscopy, such as serial two- photon tomography (STP) and light sheet fluorescence microscopy, now allow for automated rapid 3D imaging of a complete mouse brain at cellular resolution. In combination with transsynaptic viral tracers, this paves the way for high-throughput brain mapping studies, which could greatly advance our understanding of the function of the brain. Because transsynaptic tracers label synaptically connected cells, the analysis of these whole-brain scans requires detection of fluorescently labelled cells and anatomical segmentation of the data, which are very labour- and time-intensive manual tasks and prevent high-throughput analysis. This thesis presents and validates two software tools to automate anatomical segmentation and cell detection in serial two photon (STP) scans of the mouse brain. Automated mouse atlas propagation (aMAP) segments the scans into anatomical regions by matching a 3D reference atlas to the data using affine and free-form image registration. The fast automated cell counting tool (FACCT) then detects fluorescently labelled cells in the dataset with a novel approach of stepwise data reduction combined with object detection using artificial neuronal networks. The tools are optimised for large datasets and are capable of processing a 2.5TB STP scan in under two days. The performance of aMAP and FACCT is evaluated on STP scans from retrograde connectivity tracing experiments using rabies virus injections in the primary visual corte

Development and application of quantitative image analysis for preclinical MRI research

The aim of this thesis is to develop quantitative analysis methods to validate MRI and improve the detection of tumour infiltration. The major components include a description of the development the quantitative methods to better validate imaging biomarkers and detect of infiltration of tumour cells into normal tissue, which were then applied to a mouse model of glioblastoma invasion. To do this, a new histology model, called Stacked In-plane Histology (SIH), was developed to allow a quantitative analysis of MRI. Validating imaging biomarkers for glioblastoma infiltration Cancer can be defined as a disease in which a group of abnormal cells grow uncontrollably, often with fatal outcomes. According to (Cancer research UK, 2019), there are more than 363,000 new cancer cases in the UK every year, an increase from the 990 cases reported daily in 2014-2016, with only half of all patients recovering. Glioblastoma (GB) is the most frequent and malignant form of primary brain tumours with a very poor prognosis. Even with the development of modern diagnostic strategies and new therapies, the five-year survival rate is just 5%, with the median survival time only 14 months. Unfortunately, glioblastoma can affect patients at any age, including young children, but has a peak occurrence between the ages of 65 and 75 years. The standard treatment for GB consists of surgical resection, followed by radiotherapy and chemotherapy. However, the infiltration of GB cells into healthy adjacent brain tissue is a major obstacle for successful treatment, making complete removal of a tumour by surgery a difficult task, with the potential for tumour recurrence. Magnetic Resonance Imaging (MRI) is a non-invasive, multipurpose imaging tool used for the diagnosis and monitoring of cancerous tumours. It can provide morphological, physiological, and metabolic information about the tumour. Currently, MRI is the standard diagnostic tool for GB before the pathological examination of tissue from surgical resection or biopsy specimens. The standard MRI sequences used for diagnosis of GB include T2-Weighted (T2W), T1-Weighted (T1W), Fluid-Attenuated Inversion Recovery (FLAIR), and Contrast Enhance T1 gadolinium (CE-T1) scans. These conventional scans are used to localize the tumour, in addition to associated oedema and necrosis. Although these scans can identify the bulk of the tumour, it is known that they do not detect regions infiltrated by GB cells. The MRI signal depends upon many physical parameters including water content, local structure, tumbling rates, diffusion, and hypoxia (Dominietto, 2014). There has been considerable interest in identifying whether such biologically indirect image contrasts can be used as non-invasive imaging biomarkers, either for normal biological functions, pathogenic processes or pharmacological responses to therapeutic interventions (Atkinson et al., 2001). In fact, when new MRI methods are proposed as imaging biomarkers of particular diseases, it is crucial that they are validated against histopathology. In humans, such validation is limited to a biopsy, which is the gold standard of diagnosis for most types of cancer. Some types of biopsies can take an image-guided approach using MRI, Computed Tomography (CT) or Ultrasound (US). However, a biopsy may miss the most malignant region of the tumour and is difficult to repeat. Biomarker validation can be performed in preclinical disease models, where the animal can be terminated immediately after imaging for histological analysis. Here, in principle, co-registration of the biomarker images with the histopathology would allow for direct validation. However, in practice, most preclinical validation studies have been limited to using simple visual comparisons to assess the correlation between the imaging biomarker and underlying histopathology. First objective (Chapter 5): Histopathology is the gold standard for assessing non-invasive imaging biomarkers, with most validation approaches involving a qualitative visual inspection. To allow a more quantitative analysis, previous studies have attempted to co-register MRI with histology. However, these studies have focused on developing better algorithms to deal with the distortions common in histology sections. By contrast, we have taken an approach to improve the quality of the histological processing and analysis, for example, by taking into account the imaging slice orientation and thickness. Multiple histology sections were cut in the MR imaging plane to produce a Stacked In-plane Histology (SIH) map. This approach, which is applied to the next two objectives, creates a histopathology map which can be used as the gold standard to quantitatively validate imaging biomarkers. Second objective (Chapter 6): Glioblastoma is the most malignant form of primary brain tumour and recurrence following treatment is common. Non-invasive MR imaging is an important component of brain tumour diagnosis and treatment planning. Unfortunately, clinic MRI (T1W, T2W, CE-T1, and FLAIR) fails to detect regions of glioblastoma cell infiltration beyond the solid tumour region identified by contrast enhanced T1 scans. However, advanced MRI techniques such as Arterial Spin Labelling (ASL) could provide us with extra information (perfusion) which may allow better detection of infiltration. In order to assess whether local perfusion perturbation could provide a useful biomarker for glioblastoma cell infiltration, we quantitatively analysed the correlation between perfusion MRI (ASL) and stacked in-plane histology. This work used a mouse model of glioblastoma that mimics the infiltrative behaviour found in human patients. The results demonstrate the ability of perfusion imaging to probe regions of low tumour cell infiltration, while confirming the sensitivity limitations of clinic imaging modalities. Third objective (Chapter 7): It is widely hypothesised that Multiparametric MRI (mpMRI), can extract more information than is obtained from the constituent individual MR images, by reconstructing a single map that contains complementary information. Using the MRI and histology dataset from objective 2, we used a multi-regression algorithm to reconstruct a single map which was highly correlated (r>0.6) with histology. The results are promising, showing that mpMRI can better predict the whole tumour region, including the region of tumour cell infiltration