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.

Chromatin assembly in a yeast whole-cell extract

A simple in vitro system that supports chromatin assembly was developed for Saccharomyces cerevisiae. The assembly reaction is ATP-dependent, uses soluble histones and assembly factors, and generates physiologically spaced nucleosomes. We analyze the pathway of histone recruitment into nucleosomes, using this system in combination with genetic methods for the manipulation of yeast. This analysis supports the model of sequential recruitment of H3/H4 tetramers and H2A/H2B dimers into nucleosomes. Using a similar approach, we show that DNA ligase I can play an important role in template repair during assembly. These studies demonstrate the utility of this system for the combined biochemical and genetic analysis of chromatin assembly in yeast.

Structure of the thrombin complex with triabin, a lipocalin-like exosite-binding inhibitor derived from a triatomine bug

Triabin, a 142-residue protein from the saliva of the blood-sucking triatomine bug Triatoma pallidipennis, is a potent and selective thrombin inhibitor. Its stoichiometric complex with bovine α-thrombin was crystallized, and its crystal structure was solved by Patterson search methods and refined at 2.6-Å resolution to an R value of 0.184. The analysis revealed that triabin is a compact one-domain molecule essentially consisting of an eight-stranded β-barrel. The eight strands A to H are arranged in the order A-C-B-D-E-F-G-H, with the first four strands exhibiting a hitherto unobserved up-up-down-down topology. Except for the B-C inversion, the triabin fold exhibits the regular up-and-down topology of lipocalins. In contrast to the typical ligand-binding lipocalins, however, the triabin barrel encloses a hydrophobic core intersected by a unique salt-bridge cluster. Triabin interacts with thrombin exclusively via its fibrinogen-recognition exosite. Surprisingly, most of the interface interactions are hydrophobic. A prominent exception represents thrombin’s Arg-77A side chain, which extends into a hydrophobic triabin pocket forming partially buried salt bridges with Glu-128 and Asp-135 of the inhibitor. The fully accessible active site of thrombin in this complex is in agreement with its retained hydrolytic activity toward small chromogenic substrates. Impairment of thrombin’s fibrinogen converting activity or of its thrombomodulin-mediated protein C activation capacity upon triabin binding is explained by usage of overlapping interaction sites of fibrinogen, thrombomodulin, and triabin on thrombin. These data demonstrate that triabin inhibits thrombin via a novel and unique mechanism that might be of interest in the context of potential therapeutic applications.

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.

Histone H1 Reduces the Frequency of Initiation in Xenopus Egg Extract by Limiting the Assembly of Prereplication Complexes on Sperm Chromatin

Somatic histone H1 reduces both the rate and extent of DNA replication in Xenopus egg extract. We show here that H1 inhibits replication directly by reducing the number of replication forks, but not the rate of fork progression, in Xenopus sperm nuclei. Density substitution experiments demonstrate that those forks that are active in H1 nuclei elongate to form large tracts of fully replicated DNA, indicating that inhibition is due to a reduction in the frequency of initiation and not the rate or extent of elongation. The observation that H1 dramatically reduces the number of replication foci in sperm nuclei supports this view. The establishment of replication competent DNA in egg extract requires the assembly of prereplication complexes (pre-RCs) on sperm chromatin. H1 reduces binding of the pre-RC proteins, XOrc2, XCdc6, and XMcm3, to chromatin. Replication competence can be restored in these nuclei, however, only under conditions that promote the loss of H1 from chromatin and licensing of the DNA. Thus, H1 inhibits replication in egg extract by preventing the assembly of pre-RCs on sperm chromatin, thereby reducing the frequency of initiation. These data raise the interesting possibility that H1 plays a role in regulating replication origin use during Xenopus development.

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.