Computational methods to create and analyze a digital gene expression atlas of embryo development from microscopy images

Abstract The creation of atlases, or digital models where information from different subjects can be combined, is a field of increasing interest in biomedical imaging. When a single image does not contain enough information to appropriately describe the organism under study, it is then necessary to acquire images of several individuals, each of them containing complementary data with respect to the rest of the components in the cohort. This approach allows creating digital prototypes, ranging from anatomical atlases of human patients and organs, obtained for instance from Magnetic Resonance Imaging, to gene expression cartographies of embryo development, typically achieved from Light Microscopy. Within such context, in this PhD Thesis we propose, develop and validate new dedicated image processing methodologies that, based on image registration techniques, bring information from multiple individuals into alignment within a single digital atlas model. We also elaborate a dedicated software visualization platform to explore the resulting wealth of multi-dimensional data and novel analysis algo-rithms to automatically mine the generated resource in search of bio¬logical insights. In particular, this work focuses on gene expression data from developing zebrafish embryos imaged at the cellular resolution level with Two-Photon Laser Scanning Microscopy. Disposing of quantitative measurements relating multiple gene expressions to cell position and their evolution in time is a fundamental prerequisite to understand embryogenesis multi-scale processes. However, the number of gene expressions that can be simultaneously stained in one acquisition is limited due to optical and labeling constraints. These limitations motivate the implementation of atlasing strategies that can recreate a virtual gene expression multiplex. The developed computational tools have been tested in two different scenarios. The first one is the early zebrafish embryogenesis where the resulting atlas constitutes a link between the phenotype and the genotype at the cellular level. The second one is the late zebrafish brain where the resulting atlas allows studies relating gene expression to brain regionalization and neurogenesis. The proposed computational frameworks have been adapted to the requirements of both scenarios, such as the integration of partial views of the embryo into a whole embryo model with cellular resolution or the registration of anatom¬ical traits with deformable transformation models non-dependent on any specific labeling. The software implementation of the atlas generation tool (Match-IT) and the visualization platform (Atlas-IT) together with the gene expression atlas resources developed in this Thesis are to be made freely available to the scientific community. Lastly, a novel proof-of-concept experiment integrates for the first time 3D gene expression atlas resources with cell lineages extracted from live embryos, opening up the door to correlate genetic and cellular spatio-temporal dynamics. La creación de atlas, o modelos digitales, donde la información de distintos sujetos puede ser combinada, es un campo de creciente interés en imagen biomédica. Cuando una sola imagen no contiene suficientes datos como para describir apropiadamente el organismo objeto de estudio, se hace necesario adquirir imágenes de varios individuos, cada una de las cuales contiene información complementaria respecto al resto de componentes del grupo. De este modo, es posible crear prototipos digitales, que pueden ir desde atlas anatómicos de órganos y pacientes humanos, adquiridos por ejemplo mediante Resonancia Magnética, hasta cartografías de la expresión genética del desarrollo de embrionario, típicamente adquiridas mediante Microscopía Optica. Dentro de este contexto, en esta Tesis Doctoral se introducen, desarrollan y validan nuevos métodos de procesado de imagen que, basándose en técnicas de registro de imagen, son capaces de alinear imágenes y datos provenientes de múltiples individuos en un solo atlas digital. Además, se ha elaborado una plataforma de visualization específicamente diseñada para explorar la gran cantidad de datos, caracterizados por su multi-dimensionalidad, que resulta de estos métodos. Asimismo, se han propuesto novedosos algoritmos de análisis y minería de datos que permiten inspeccionar automáticamente los atlas generados en busca de conclusiones biológicas significativas. En particular, este trabajo se centra en datos de expresión genética del desarrollo embrionario del pez cebra, adquiridos mediante Microscopía dos fotones con resolución celular. Disponer de medidas cuantitativas que relacionen estas expresiones genéticas con las posiciones celulares y su evolución en el tiempo es un prerrequisito fundamental para comprender los procesos multi-escala característicos de la morfogénesis. Sin embargo, el número de expresiones genéticos que pueden ser simultáneamente etiquetados en una sola adquisición es reducido debido a limitaciones tanto ópticas como del etiquetado. Estas limitaciones requieren la implementación de estrategias de creación de atlas que puedan recrear un multiplexado virtual de expresiones genéticas. Las herramientas computacionales desarrolladas han sido validadas en dos escenarios distintos. El primer escenario es el desarrollo embrionario temprano del pez cebra, donde el atlas resultante permite constituir un vínculo, a nivel celular, entre el fenotipo y el genotipo de este organismo modelo. El segundo escenario corresponde a estadios tardíos del desarrollo del cerebro del pez cebra, donde el atlas resultante permite relacionar expresiones genéticas con la regionalización del cerebro y la formación de neuronas. La plataforma computacional desarrollada ha sido adaptada a los requisitos y retos planteados en ambos escenarios, como la integración, a resolución celular, de vistas parciales dentro de un modelo consistente en un embrión completo, o el alineamiento entre estructuras de referencia anatómica equivalentes, logrado mediante el uso de modelos de transformación deformables que no requieren ningún marcador específico. Está previsto poner a disposición de la comunidad científica tanto la herramienta de generación de atlas (Match-IT), como su plataforma de visualización (Atlas-IT), así como las bases de datos de expresión genética creadas a partir de estas herramientas. Por último, dentro de la presente Tesis Doctoral, se ha incluido una prueba conceptual innovadora que permite integrar los mencionados atlas de expresión genética tridimensionales dentro del linaje celular extraído de una adquisición in vivo de un embrión. Esta prueba conceptual abre la puerta a la posibilidad de correlar, por primera vez, las dinámicas espacio-temporales de genes y células.

Monte-Carlo dosimetry on a realistic cell monolayer geometry exposed to alpha particles

The energy and specific energy absorbed in the main cell compartments (nucleus and cytoplasm) in typical radiobiology experiments are usually estimated by calculations as they are not accessible for a direct measurement. In most of the work, the cell geometry is modelled using the combination of simple mathematical volumes. We propose a method based on high resolution confocal imaging and ion beam analysis (IBA) in order to import realistic cell nuclei geometries in Monte-Carlo simulations and thus take into account the variety of different geometries encountered in a typical cell population. Seventy-six cell nuclei have been imaged using confocal microscopy and their chemical composition has been measured using IBA. A cellular phantom was created from these data using the ImageJ image analysis software and imported in the Geant4 Monte-Carlo simulation toolkit. Total energy and specific energy distributions in the 76 cell nuclei have been calculated for two types of irradiation protocols: a 3 MeV alpha particle microbeam used for targeted irradiation and a 239Pu alpha source used for large angle random irradiation. Qualitative images of the energy deposited along the particle tracks have been produced and show good agreement with images of DNA double strand break signalling proteins obtained experimentally. The methodology presented in this paper provides microdosimetric quantities calculated from realistic cellular volumes. It is based on open-source oriented software that is publicly available.

Spatio-temporal filtering with morphological operators for robust cell migration estimation in "in-vivo" images

The understanding of the embryogenesis in living systems requires reliable quantitative analysis of the cell migration throughout all the stages of development. This is a major challenge of the "in-toto" reconstruction based on different modalities of "in-vivo" imaging techniques -spatio-temporal resolution and image artifacts and noise. Several methods for cell tracking are available, but expensive manual interaction -time and human resources- is always required to enforce coherence. Because of this limitation it is necessary to restrict the experiments or assume an uncontrolled error rate. Is it possible to obtain automated reliable measurements of migration? can we provide a seed for biologists to complete cell lineages efficiently? We propose a filtering technique that considers trajectories as spatio-temporal connected structures that prunes out those that might introduce noise and false positives by using multi-dimensional morphological operators.

Automated processing of zebrafish imaging data: a survey

Due to the relative transparency of its embryos and larvae, the zebrafish is an ideal model organism for bioimaging approaches in vertebrates. Novel microscope technologies allow the imaging of developmental processes in unprecedented detail, and they enable the use of complex image-based read-outs for high-throughput/high-content screening. Such applications can easily generate Terabytes of image data, the handling and analysis of which becomes a major bottleneck in extracting the targeted information. Here, we describe the current state of the art in computational image analysis in the zebrafish system. We discuss the challenges encountered when handling high-content image data, especially with regard to data quality, annotation, and storage. We survey methods for preprocessing image data for further analysis, and describe selected examples of automated image analysis, including the tracking of cells during embryogenesis, heartbeat detection, identification of dead embryos, recognition of tissues and anatomical landmarks, and quantification of behavioral patterns of adult fish. We review recent examples for applications using such methods, such as the comprehensive analysis of cell lineages during early development, the generation of a three-dimensional brain atlas of zebrafish larvae, and high-throughput drug screens based on movement patterns. Finally, we identify future challenges for the zebrafish image analysis community, notably those concerning the compatibility of algorithms and data formats for the assembly of modular analysis pipelines.

Methods for the analysis of multi-scale cell dynamics from fluorescence microscopy images

La embriogénesis es el proceso mediante el cual una célula se convierte en un ser un vivo. A lo largo de diferentes etapas de desarrollo, la población de células va proliferando a la vez que el embrión va tomando forma y se configura. Esto es posible gracias a la acción de varios procesos genéticos, bioquímicos y mecánicos que interaccionan y se regulan entre ellos formando un sistema complejo que se organiza a diferentes escalas espaciales y temporales. Este proceso ocurre de manera robusta y reproducible, pero también con cierta variabilidad que permite la diversidad de individuos de una misma especie. La aparición de la microscopía de fluorescencia, posible gracias a proteínas fluorescentes que pueden ser adheridas a las cadenas de expresión de las células, y los avances en la física óptica de los microscopios han permitido observar este proceso de embriogénesis in-vivo y generar secuencias de imágenes tridimensionales de alta resolución espacio-temporal. Estas imágenes permiten el estudio de los procesos de desarrollo embrionario con técnicas de análisis de imagen y de datos, reconstruyendo dichos procesos para crear la representación de un embrión digital. Una de las más actuales problemáticas en este campo es entender los procesos mecánicos, de manera aislada y en interacción con otros factores como la expresión genética, para que el embrión se desarrolle. Debido a la complejidad de estos procesos, estos problemas se afrontan mediante diferentes técnicas y escalas específicas donde, a través de experimentos, pueden hacerse y confrontarse hipótesis, obteniendo conclusiones sobre el funcionamiento de los mecanismos estudiados. Esta tesis doctoral se ha enfocado sobre esta problemática intentando mejorar las metodologías del estado del arte y con un objetivo específico: estudiar patrones de deformación que emergen del movimiento organizado de las células durante diferentes estados del desarrollo del embrión, de manera global o en tejidos concretos. Estudios se han centrado en la mecánica en relación con procesos de señalización o interacciones a nivel celular o de tejido. En este trabajo, se propone un esquema para generalizar el estudio del movimiento y las interacciones mecánicas que se desprenden del mismo a diferentes escalas espaciales y temporales. Esto permitiría no sólo estudios locales, si no estudios sistemáticos de las escalas de interacción mecánica dentro de un embrión. Por tanto, el esquema propuesto obvia las causas de generación de movimiento (fuerzas) y se centra en la cuantificación de la cinemática (deformación y esfuerzos) a partir de imágenes de forma no invasiva. Hoy en día las dificultades experimentales y metodológicas y la complejidad de los sistemas biológicos impiden una descripción mecánica completa de manera sistemática. Sin embargo, patrones de deformación muestran el resultado de diferentes factores mecánicos en interacción con otros elementos dando lugar a una organización mecánica, necesaria para el desarrollo, que puede ser cuantificado a partir de la metodología propuesta en esta tesis. La metodología asume un medio continuo descrito de forma Lagrangiana (en función de las trayectorias de puntos materiales que se mueven en el sistema en lugar de puntos espaciales) de la dinámica del movimiento, estimado a partir de las imágenes mediante métodos de seguimiento de células o de técnicas de registro de imagen. Gracias a este esquema es posible describir la deformación instantánea y acumulada respecto a un estado inicial para cualquier dominio del embrión. La aplicación de esta metodología a imágenes 3D + t del pez zebra sirvió para desvelar estructuras mecánicas que tienden a estabilizarse a lo largo del tiempo en dicho embrión, y que se organizan a una escala semejante al del mapa de diferenciación celular y con indicios de correlación con patrones de expresión genética. También se aplicó la metodología al estudio del tejido amnioserosa de la Drosophila (mosca de la fruta) durante el cierre dorsal, obteniendo indicios de un acoplamiento entre escalas subcelulares, celulares y supracelulares, que genera patrones complejos en respuesta a la fuerza generada por los esqueletos de acto-myosina. En definitiva, esta tesis doctoral propone una estrategia novedosa de análisis de la dinámica celular multi-escala que permite cuantificar patrones de manera inmediata y que además ofrece una representación que reconstruye la evolución de los procesos como los ven las células, en lugar de como son observados desde el microscopio. Esta metodología por tanto permite nuevas formas de análisis y comparación de embriones y tejidos durante la embriogénesis a partir de imágenes in-vivo. ABSTRACT The embryogenesis is the process from which a single cell turns into a living organism. Through several stages of development, the cell population proliferates at the same time the embryo shapes and the organs develop gaining their functionality. This is possible through genetic, biochemical and mechanical factors that are involved in a complex interaction of processes organized in different levels and in different spatio-temporal scales. The embryogenesis, through this complexity, develops in a robust and reproducible way, but allowing variability that makes possible the diversity of living specimens. The advances in physics of microscopes and the appearance of fluorescent proteins that can be attached to expression chains, reporting about structural and functional elements of the cell, have enabled for the in-vivo observation of embryogenesis. The imaging process results in sequences of high spatio-temporal resolution 3D+time data of the embryogenesis as a digital representation of the embryos that can be further analyzed, provided new image processing and data analysis techniques are developed. One of the most relevant and challenging lines of research in the field is the quantification of the mechanical factors and processes involved in the shaping process of the embryo and their interactions with other embryogenesis factors such as genetics. Due to the complexity of the processes, studies have focused on specific problems and scales controlled in the experiments, posing and testing hypothesis to gain new biological insight. However, methodologies are often difficult to be exported to study other biological phenomena or specimens. This PhD Thesis is framed within this paradigm of research and tries to propose a systematic methodology to quantify the emergent deformation patterns from the motion estimated in in-vivo images of embryogenesis. Thanks to this strategy it would be possible to quantify not only local mechanisms, but to discover and characterize the scales of mechanical organization within the embryo. The framework focuses on the quantification of the motion kinematics (deformation and strains), neglecting the causes of the motion (forces), from images in a non-invasive way. Experimental and methodological challenges hamper the quantification of exerted forces and the mechanical properties of tissues. However, a descriptive framework of deformation patterns provides valuable insight about the organization and scales of the mechanical interactions, along the embryo development. Such a characterization would help to improve mechanical models and progressively understand the complexity of embryogenesis. This framework relies on a Lagrangian representation of the cell dynamics system based on the trajectories of points moving along the deformation. This approach of analysis enables the reconstruction of the mechanical patterning as experienced by the cells and tissues. Thus, we can build temporal profiles of deformation along stages of development, comprising both the instantaneous events and the cumulative deformation history. The application of this framework to 3D + time data of zebrafish embryogenesis allowed us to discover mechanical profiles that stabilized through time forming structures that organize in a scale comparable to the map of cell differentiation (fate map), and also suggesting correlation with genetic patterns. The framework was also applied to the analysis of the amnioserosa tissue in the drosophila’s dorsal closure, revealing that the oscillatory contraction triggered by the acto-myosin network organized complexly coupling different scales: local force generation foci, cellular morphology control mechanisms and tissue geometrical constraints. In summary, this PhD Thesis proposes a theoretical framework for the analysis of multi-scale cell dynamics that enables to quantify automatically mechanical patterns and also offers a new representation of the embryo dynamics as experienced by cells instead of how the microscope captures instantaneously the processes. Therefore, this framework enables for new strategies of quantitative analysis and comparison between embryos and tissues during embryogenesis from in-vivo images.

PRIM: Proximity imaging of green fluorescent protein-tagged polypeptides

We report a serendipitous discovery that extends the impressive catalog of reporter functions performed by green fluorescent protein (GFP) or its derivatives. When two GFP molecules are brought into proximity, changes in the relative intensities of green fluorescence emitted upon excitation at 395 vs. 475 nm result. These spectral changes provide a sensitive ratiometric index of the extent of self-association that can be exploited to quantitatively image homo-oligomerization or clustering processes of GFP-tagged proteins in vivo. The method, which we term proximity imaging (PRIM), complements fluorescence resonance energy transfer between a blue fluorescent protein donor and a GFP acceptor, a powerful method for imaging proximity relationships between different proteins. However, unlike fluorescence resonance energy transfer (which is a spectral interaction), PRIM depends on direct contact between two GFP modules, which can lead to structural perturbations and concomitant spectral changes within a module. Moreover, the precise spatial arrangement of the GFP molecules within a given dimer determines the magnitude and direction of the spectral change. We have used PRIM to detect FK1012-induced dimerization of GFP fused to FK506-binding protein and clustering of glycosylphosphatidylinositol-anchored GFP at cell surfaces.

Cell surface expression of the Alzheimer disease-related presenilin proteins

The presenilin proteins PS-1 and PS-2 are crucially involved in Alzheimer disease (AD), but their molecular functions are not known. They are integral membrane proteins, but whether they can be expressed at the surface of cells has been in dispute. Here we show by immunofluorescence experiments, using anti-peptide antibodies specific for either PS-1 or PS-2, that live cultured DAMI cells and differentiated human NT2N neuronal cells are specifically immunolabeled for their endogenous as well as transfected presenilins, although the cells cannot be immunolabeled for their intracellular tubulin, unless they are first fixed and permeabilized. These and other results establish that portions of the presenilins are indeed expressed at the surfaces of these cells. These findings support our previous proposal that the presenilins on the surface of a cell engage in intercellular interactions with the β-amyloid precursor protein on the surface of a neighboring cell, as a critical step in the molecular and cellular mechanisms that lead to AD.

Production of identical sextuplet mice by transferring metaphase nuclei from four-cell embryos

Mouse clones were produced by serial nuclear transfer commencing with the transfer of four-cell nuclei at metaphase into unfertilized ooplasts. The donor four-cell-stage nuclei were synchronized in metaphase with nocodazole. The oocytes receiving a four-cell nucleus at metaphase formed two nuclei after artificial activation and inhibition of cytokinesis with cytochalasin B. To obtain embryos with diploid sets of chromosomes, nuclei from each reconstructed embryo were transferred individually into separate enucleated fertilized one-cell embryos, thus doubling the number of identical embryos. This procedure produced a high frequency of development of reconstructed embryos to the blastocyst stage. Of 11 sets of identical embryos produced by serial nuclear transplantation, 83% developed into blastocysts, including three sets of identical septuplet blastocysts. After transfer to recipient mice, a total of 25 (57%) live young were obtained, which included one set of identical sextuplet and two sets of identical quadruplet mice.

STEAP: A prostate-specific cell-surface antigen highly expressed in human prostate tumors

In search of novel genes expressed in metastatic prostate cancer, we subtracted cDNA isolated from benign prostatic hypertrophic tissue from cDNA isolated from a prostate cancer xenograft model that mimics advanced disease. One novel gene that is highly expressed in advanced prostate cancer encodes a 339-amino acid protein with six potential membrane-spanning regions flanked by hydrophilic amino- and carboxyl-terminal domains. This structure suggests a potential function as a channel or transporter protein. This gene, named STEAP for six-transmembrane epithelial antigen of the prostate, is expressed predominantly in human prostate tissue and is up-regulated in multiple cancer cell lines, including prostate, bladder, colon, ovarian, and Ewing sarcoma. Immunohistochemical analysis of clinical specimens demonstrates significant STEAP expression at the cell–cell junctions of the secretory epithelium of prostate and prostate cancer cells. Little to no staining was detected at the plasma membranes of normal, nonprostate human tissues, except for bladder tissue, which expressed low levels of STEAP at the cell membrane. Protein analysis located STEAP at the cell surface of prostate-cancer cell lines. Our results support STEAP as a cell-surface tumor-antigen target for prostate cancer therapy and diagnostic imaging.

Birefringence Imaging Directly Reveals Architectural Dynamics of Filamentous Actin in Living Growth ConesV⃞

We have investigated the dynamic behavior of cytoskeletal fine structure in the lamellipodium of nerve growth cones using a new type of polarized light microscope (the Pol-Scope). Pol-Scope images display with exquisite resolution and definition birefringent fine structures, such as filaments and membranes, without having to treat the cell with exogenous dyes or fluorescent labels. Furthermore, the measured birefringence of protein fibers in the thin lamellipodial region can be interpreted in terms of the number of filaments in the bundles. We confirmed that birefringent fibers are actin-based using conventional fluorescence-labeling methods. By recording movies of time-lapsed Pol-Scope images, we analyzed the creation and dynamic composition of radial fibers, filopodia, and intrapodia in advancing growth cones. The strictly quantitative information available in time-lapsed Pol-Scope images confirms previously deduced behavior and provides new insight into the architectural dynamics of filamentous actin.

Real-Time Imaging of the Axonal Transport of Granules Containing a Tissue Plasminogen Activator/Green Fluorescent Protein Hybrid

A hybrid protein, tPA/GFP, consisting of rat tissue plasminogen activator (tPA) and green fluorescent protein (GFP) was expressed in PC12 cells and used to study the distribution, secretory behavior, and dynamics of secretory granules containing tPA in living cells with a neuronal phenotype. High-resolution images demonstrate that tPA/GFP has a growth cone-biased distribution in differentiated cells and that tPA/GFP is transported in granules of the regulated secretory pathway that colocalize with granules containing secretogranin II. Time-lapse images of secretion reveal that secretagogues induce substantial loss of cellular tPA/GFP fluorescence, most importantly from growth cones. Time-lapse images of the axonal transport of granules containing tPA/GFP reveal a surprising complexity to granule dynamics. Some granules undergo canonical fast axonal transport; others move somewhat more slowly, especially in highly fluorescent neurites. Most strikingly, granules traffic bidirectionally along neurites to an extent that depends on granule accumulation, and individual granules can reverse their direction of motion. The retrograde component of this bidirectional transport may help to maintain cellular homeostasis by transporting excess tPA/GFP back toward the cell body. The results presented here provide a novel view of the axonal transport of secretory granules. In addition, the results suggest that tPA is targeted for regulated secretion from growth cones of differentiated cells, strategically positioning tPA to degrade extracellular barriers or to activate other barrier-degrading proteases during axonal elongation.

Myosin Heavy Chain Phosphorylation Sites Regulate Myosin Localization during Cytokinesis in Live Cells

Conventional myosin II plays a fundamental role in the process of cytokinesis where, in the form of bipolar thick filaments, it is thought to be the molecular motor that generates the force necessary to divide the cell. In Dictyostelium, the formation of thick filaments is regulated by the phosphorylation of three threonine residues in the tail region of the myosin heavy chain. We report here on the effects of this regulation on the localization of myosin in live cells undergoing cytokinesis. We imaged fusion proteins of the green-fluorescent protein with wild-type myosin and with myosins where the three critical threonines had been changed to either alanine or aspartic acid. We provide evidence that thick filament formation is required for the accumulation of myosin in the cleavage furrow and that if thick filaments are overproduced, this accumulation is markedly enhanced. This suggests that myosin localization in dividing cells is regulated by myosin heavy chain phosphorylation.

Imaging adenoviral-directed reporter gene expression in living animals with positron emission tomography

We are developing quantitative assays to repeatedly and noninvasively image expression of reporter genes in living animals, using positron emission tomography (PET). We synthesized positron-emitting 8-[18F]fluoroganciclovir (FGCV) and demonstrated that this compound is a substrate for the herpes simplex virus 1 thymidine kinase enzyme (HSV1-TK). Using positron-emitting FGCV as a PET reporter probe, we imaged adenovirus-directed hepatic expression of the HSV1-tk reporter gene in living mice. There is a significant positive correlation between the percent injected dose of FGCV retained per gram of liver and the levels of hepatic HSV1-tk reporter gene expression (r2 > 0.80). Over a similar range of HSV1-tk expression in vivo, the percent injected dose retained per gram of liver was 0–23% for ganciclovir and 0–3% for FGCV. Repeated, noninvasive, and quantitative imaging of PET reporter gene expression should be a valuable tool for studies of human gene therapy, of organ/cell transplantation, and of both environmental and behavioral modulation of gene expression in transgenic mice.