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.

Assembling models of embryo development: Image analysis and the construction of digital atlases

Digital atlases of animal development provide a quantitative description of morphogenesis, opening the path toward processes modeling. Prototypic atlases offer a data integration framework where to gather information from cohorts of individuals with phenotypic variability. Relevant information for further theoretical reconstruction includes measurements in time and space for cell behaviors and gene expression. The latter as well as data integration in a prototypic model, rely on image processing strategies. Developing the tools to integrate and analyze biological multidimensional data are highly relevant for assessing chemical toxicity or performing drugs preclinical testing. This article surveys some of the most prominent efforts to assemble these prototypes, categorizes them according to salient criteria and discusses the key questions in the field and the future challenges toward the reconstruction of multiscale dynamics in model organisms.

Reproductive and nutritional management on ovarian response and embryo quality on rabbit does

Rabbit does in modern rabbitries are under intensive reproductive rhythms. Females are high milk producers with high energetic expenses due to the extensive overlap between lactation and gestation. This situation leads to a negative energy balance with a mobilization of body fat especially in primiparous rabbit does. Poor body condition and poor health status severely affect the reproductive features (fertility rate and lifespan of the doe as well as ovarian physiology). This paper reviews some reproductive and nutritional approaches used in the last years to improve the reproductive performance of rabbit females, mainly focusing on the influence on ovarian response and embryo quality and with emphasis on epigenetic modifications in pre-implantation embryos and offspring consequences.

The identification of a 9-cis retinol dehydrogenase in the mouse embryo reveals a pathway for synthesis of 9-cis retinoic acid

The ligand-controlled retinoic acid (RA) receptors and retinoid X receptors are important for several physiological processes, including normal embryonic development, but little is known about how their ligands, all-trans and 9-cis RA, are generated. Here we report the identification of a stereo-specific 9-cis retinol dehydrogenase, which is abundantly expressed in embryonic tissues known to be targets in the retinoid signaling pathway. The membrane-bound enzyme is a member of the short-chain alcohol dehydrogenase/reductase superfamily, able to oxidize 9-cis retinol into 9-cis retinaldehyde, an intermediate in 9-cis RA biosynthesis. Analysis by nonradioactive in situ hybridization in mouse embryos shows that expression of the enzyme is temporally and spatially well controlled during embryogenesis with prominent expression in parts of the developing central nervous system, sensory organs, somites and myotomes, and several tissues of endodermal origin. The identification of this enzyme reveals a pathway in RA biosynthesis, where 9-cis retinol is generated for subsequent oxidation to 9-cis RA.

Postgastrulation Smad2-deficient embryos show defects in embryo turning and anterior morphogenesis

SMAD2 is a member of the transforming growth factor β and activin-signaling pathway. To examine the role of Smad2 in postgastrulation development, we independently generated mice with a null mutation in this gene. Smad2-deficient embryos die around day 7.5 of gestation because of failure of gastrulation and failure to establish an anterior–posterior (A-P) axis. Expression of the homeobox gene Hex (the earliest known marker of the A-P polarity and the prospective head organizer) was found to be missing in Smad2-deficient embryos. Homozygous mutant embryos and embryonic stem cells formed mesoderm derivatives revealing that mesoderm induction is SMAD2 independent. In the presence of wild-type extraembryonic tissues, Smad2-deficient embryos developed beyond 7.5 and up to 10.5 days postcoitum, demonstrating a requirement for SMAD2 in extraembryonic tissues for the generation of an A-P axis and gastrulation. The rescued postgastrulation embryos showed malformation of head structures, abnormal embryo turning, and cyclopia. Our results show that Smad2 expression is required at several stages during embryogenesis.

A serotonin transporter gene intron 2 polymorphic region, correlated with affective disorders, has allele-dependent differential enhancer-like properties in the mouse embryo

Polymorphic regions consisting of a variable number of tandem repeats within intron 2 of the gene coding for the serotonin transporter protein 5-HTT have been associated with susceptibility to affective disorders. We have cloned two of these intronic polymorphisms, Stin2.10 and Stin2.12, into an expression vector containing a heterologous minimal promoter and the bacterial LacZ reporter gene. These constructs were then used to produce transgenic mice. In embryonic day 10.5 embryos, both Stin2.10 and Stin2.12 produced consistent β-galactosidase expression in the embryonic midbrain, hindbrain, and spinal cord floor plate. However, we observed that the levels of β-galactosidase expression produced by both the Stin2.10 and Stin2.12 within the rostral hindbrain differed significantly at embryonic day 10.5. Our data suggest that these polymorphic variable number of tandem repeats regions act as transcriptional regulators and have allele-dependent differential enhancer-like properties within an area of the hindbrain where the 5-HTT gene is known to be transcribed at this stage of development.

TAFII mutations disrupt Dorsal activation in the Drosophila embryo

In this study, we present evidence that the Dorsal activator interacts with limiting amounts of the TFIID complex in the Drosophila embryo. In vitro transcription reactions and protein binding assays implicate the TAFII110 and TAFII60 subunits of the TFIID complex in contributing to Dorsal-mediated activation. Mutations in TAFII110 and TAFII60 result in altered patterns of snail and twist transcription in embryos derived from dl/+ females. These results suggest that TAFIIs contribute to the activation of transcription in vivo and support the hypothesis that subunits of TFIID may serve as targets of enhancer binding proteins.

5-Methylcytosine-DNA glycosylase activity is present in a cloned G/T mismatch DNA glycosylase associated with the chicken embryo DNA demethylation complex

We previously have shown that DNA demethylation by chicken embryo 5-methylcytosine DNA glycosylase (5-MCDG) needs both RNA and proteins. One of these proteins is a RNA helicase. Further peptides were sequenced, and three of them are identical to the mammalian G/T mismatch DNA glycosylase. A 3,233-bp cDNA coding for the chicken homologue of human G/T mismatch DNA glycosylase was isolated and sequenced. The derived amino acid sequence (408 aa) shows 80% identity with the human G/T mismatch DNA glycosylase, and both the C and N-terminal parts have about 50% identity. As for the highly purified chicken embryo DNA demethylation complex the recombinant protein expressed in Escherichia coli has both G/T mismatch and 5-MCDG activities. The recombinant protein has the same substrate specificity as the chicken embryo 5-MCDG where hemimethylated DNA is a better substrate than symmetrically methylated CpGs. The activity ratio of G/T mismatch and 5-MCDG is about 30:1 for the recombinant protein expressed in E. coli and 3:1 for the purified enzyme from chicken embryos. The incubation of a recombinant CpG-rich RNA isolated from the purified DNA demethylation complex with the recombinant enzyme strongly inhibits G/T mismatch glycosylase while slightly stimulating the activity of 5-MCDG. Deletion mutations indicate that G/T mismatch and 5-MCDG activities share the same areas of the N- and C-terminal parts of the protein. In reconstitution experiments RNA helicase in the presence of recombinant RNA and ATP potentiates the activity of 5-MCDG.

Expression and parent-of-origin effects for FIS2, MEA, and FIE in the endosperm and embryo of developing Arabidopsis seeds

The promoters of MEA (FIS1), FIS2, and FIE (FIS3), genes that repress seed development in the absence of pollination, were fused to β-glucuronidase (GUS) to study their activity pattern. The FIS2∷GUS product is found in the embryo sac, in each of the polar cell nuclei, and in the central cell nucleus. After pollination, the maternally derived FIS2∷GUS protein occurs in the nuclei of the cenocytic endosperm. Before cellularization of the endosperm, activity is terminated in the micropylar and central nuclei of the endosperm and subsequently in the nuclei of the chalazal cyst. MEA∷GUS has a pattern of activity similar to that of FIS2∷GUS, but FIE∷GUS protein is found in many tissues, including the prepollination embryo sac, and in embryo and endosperm postpollination. The similarity in mutant phenotypes; the activity of FIE, MEA, and FIS2 in the same cells in the embryo sac; and the fact that MEA and FIE proteins interact in a yeast two-hybrid system suggest that these proteins operate in the same system of control of seed development. Maternal and not paternal FIS2∷GUS, MEA∷GUS, and FIE∷GUS show activity in early endosperm, so these genes may be imprinted. When fis2, mea, and fie mutants are pollinated, seed development is arrested at the heart embryo stage. The seed arrest of mea and fis2 is avoided when they are fertilized by a low methylation parent. The wild-type alleles of MEA or FIS2 are not required. The parent-of-origin-determined differential activity of MEA, FIS2, and FIE is not dependent on DNA methylation, but methylation does control some gene(s) that have key roles in seed development.

A modifier screen in the eye reveals control genes for Krüppel activity in the Drosophila embryo

Irregular facets (If) is a dominant mutation of Drosophila that results in small eyes with fused ommatidia. Previous results showed that the gene Krüppel (Kr), which is best known for its early segmentation function, is expressed ectopically in If mutant eye discs. However, it was not known whether ectopic Kr activity is either the cause or the result of the If mutation. Here, we show that If is a gain-of-function allele of Kr. We then used the If mutation in a genetic screen to identify dominant enhancers and suppressors of Kr activity on the third chromosome. Of 30 identified Kr-interacting loci, two were cloned, and we examined whether they also represent components of a natural Kr-dependent developmental pathway of the embryo. We show that the two genes, eyelid (eld) and extramacrochaetae (emc), which encode a Bright family-type DNA binding protein and a helix-loop-helix factor, respectively, are necessary to achieve the singling-out of a unique Kr-expressing cell during the development of the Malpighian tubules, the excretory organs of the fly. The results indicate that the Kr gain-of-function mutation If provides a tool to identify genes that are active during eye development and that a number of them function also in the control of Kr-dependent developmental processes.

Segment-specific pattern of sympathetic preganglionic projections in the chicken embryo spinal cord is altered by retinoids

Sympathetic preganglionic neurons exhibit segment-specific projections. Preganglionic neurons located in rostral spinal segments project rostrally within the sympathetic chain, those located in caudal spinal segments project caudally, and those in midthoracic segments project either rostrally or caudally in segmentally graded proportions. Moreover, rostrally and caudally projecting preganglionic neurons are skewed toward the rostral and caudal regions, respectively, of each midthoracic segment. The mechanisms that establish these segment-specific projections are unknown. Here we show that experimental manipulation of retinoid signaling in the chicken embryo alters the segment-specific pattern of sympathetic preganglionic projections and that this effect is mediated by the somitic mesoderm. Application of exogenous retinoic acid to a single rostral thoracic somite decreases the number of rostrally projecting preganglionic neurons at that level. Conversely, disrupting endogenous synthesis of retinoic acid in a single caudal thoracic somite increases the number of rostrally projecting preganglionic neurons at that level. The number of caudally projecting neurons does not change in either case, indicating that the effect is specific for rostrally projecting preganglionic neurons. These results indicate that the sizes of the rostrally and caudally projecting populations may be independently regulated by different factors. Opposing gradients of such factors along the longitudinal axis of the thoracic region of the embryo could be sufficient, in combination, to determine the segment-specific identity of preganglionic projections.

Endogenous E2F-1 promotes timely G0 exit of resting mouse embryo fibroblasts

Much evidence strongly suggests a positive role for one or more E2F species in the control of exit from G0/G1. Results described here provide direct evidence that endogenous E2F-1, as predicted, contributes to progression from G0 to S. By contrast, cycling cells lacking an intact E2F-1 gene demonstrated normal cell cycle distribution. Therefore, E2F-1 exerts a unique function leading to timely G0 exit of resting cultured primary cells, while at the same time being unnecessary for normal G1 to S phase progression of cycling cells.

From cell death to embryo arrest: Mathematical models of human preimplantation embryo development

Human preimplantation embryos exhibit high levels of apoptotic cells and high rates of developmental arrest during the first week in vitro. The relation between the two is unclear and difficult to determine by conventional experimental approaches, partly because of limited numbers of embryos. We apply a mixture of experiment and mathematical modeling to show that observed levels of cell death can be reconciled with the high levels of embryo arrest seen in the human only if the developmental competence of embryos is already established at the zygote stage, and environmental factors merely modulate this. This suggests that research on improving in vitro fertilization success rates should move from its current concentration on optimizing culture media to focus more on the generation of a healthy zygote and on understanding the mechanisms that cause chromosomal and other abnormalities during early cleavage stages.

Production of zebrafish germ-line chimeras from embryo cell cultures

Although the zebrafish possesses many characteristics that make it a valuable model for genetic studies of vertebrate development, one deficiency of this model system is the absence of methods for cell-mediated gene transfer and targeted gene inactivation. In mice, embryonic stem cell cultures are routinely used for gene transfer and provide the advantage of in vitro selection for rare events such as homologous recombination and targeted mutation. Transgenic animals possessing a mutated copy of the targeted gene are generated when the selected cells contribute to the germ line of a chimeric embryo. Although zebrafish embryo cell cultures that exhibit characteristics of embryonic stem cells have been described, successful contribution of the cells to the germ-cell lineage of a host embryo has not been reported. In this study, we demonstrate that short-term zebrafish embryo cell cultures maintained in the presence of cells from a rainbow trout spleen cell line (RTS34st) are able to produce germ-line chimeras when introduced into a host embryo. Messenger RNA encoding the primordial germ-cell marker, vasa, was present for more than 30 days in embryo cells cocultured with RTS34st cells or their conditioned medium and disappeared by 5 days in the absence of the spleen cells. The RTS34st cells also inhibited melanocyte and neuronal cell differentiation in the embryo cell cultures. These results suggest that the RTS34st splenic–stromal cell line will be a valuable tool in the development of a cell-based gene transfer approach to targeted gene inactivation in zebrafish.

PEA3 sites within the progression elevated gene-3 (PEG-3) promoter and mitogen-activated protein kinase contribute to differential PEG-3 expression in Ha-ras and v-raf oncogene transformed rat embryo cells

Transformation of normal cloned rat embryo fibroblast (CREF) cells with cellular oncogenes results in acquisition of anchorage-independent growth and oncogenic potential in nude mice. These cellular changes correlate with an induction in the expression of a cancer progression-promoting gene, progression elevated gene-3 (PEG-3). To define the mechanism of activation of PEG-3 as a function of transformation by the Ha-ras and v-raf oncogenes, evaluations of the signaling and transcriptional regulation of the ~2.0 kb promoter region of the PEG-3 gene, PEG-Prom, was undertaken. The full-length and various mutated regions of the PEG-Prom were linked to a luciferase reporter construct and tested for promoter activity in CREF and oncogene-transformed CREF cells. An analysis was also performed using CREF cells doubly transformed with Ha-ras and the Ha-ras specific suppressor gene Krev-1, which inhibits the transformed phenotype in vitro. These assays document an association between expression of the transcription regulator PEA3 and PEG-3. The levels of PEA3 and PEG-3 RNA and proteins are elevated in the oncogenically transformed CREF cells, and reduced in transformation and tumorigenic suppressed Ha-ras/Krev-1 doubly transformed CREF cells. Enhanced tumorigenic behavior, PEG-3 promoter function and PEG-3 expression in Ha-ras transformed cells were all dependent upon increased activity within the mitogen-activated protein kinase (MAPK) pathway. Electrophoretic mobility shift assays and DNase I footprinting experiments indicate that PEA3 binds to sites within the PEG-Prom in transformed rodent cells in an area adjacent to the TATA box in a MAPK-dependent fashion. These findings demonstrate an association between Ha-ras and v-raf transformation of CREF cells with elevated PEA3 and PEG-3 expression, and they implicate MAPK signaling via PEA3 as a signaling cascade involved in activation of the PEG-Prom.