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.

Quantitative insight into proliferation and differentiation of oligodendrocyte type 2 astrocyte progenitor cells in vitro

As part of our attempts at understanding fundamental principles that underlie the generation of nondividing terminally differentiated progeny from dividing precursor cells, we have developed approaches to a quantitative analysis of proliferation and differentiation of oligodendrocyte type 2 astrocyte (O-2A) progenitor cells at the clonal level. Owing to extensive previous studies of clonal differentiation in this lineage, O-2A progenitor cells represent an excellent system for such an analysis. Previous studies have resulted in two competing hypotheses; one of them suggests that progenitor cell differentiation is symmetric, the other hypothesis introduces an asymmetric process of differentiation. We propose a general model that incorporates both such extreme hypotheses as special cases. Our analysis of experimental data has shown, however, that neither of these extreme cases completely explains the observed kinetics of O-2A progenitor cell proliferation and oligodendrocyte generation in vitro. Instead, our results indicate that O-2A progenitor cells become competent for differentiation after they complete a certain number of critical mitotic cycles that represent a period of symmetric development. This number varies from clone to clone and may be thought of as a random variable; its probability distribution was estimated from experimental data. Those O-2A cells that have undergone the critical divisions then may differentiate into an oligodendrocyte in each of the subsequent mitotic cycles with a certain probability, thereby exhibiting the asymmetric type of differentiation.

Generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification

We have developed a technique called the generation of longer cDNA fragments from serial analysis of gene expression (SAGE) tags for gene identification (GLGI), to convert SAGE tags of 10 bases into their corresponding 3′ cDNA fragments covering hundred bases. A primer containing the 10-base SAGE tag is used as the sense primer, and a single base anchored oligo(dT) primer is used as an antisense primer in PCR, together with Pfu DNA polymerase. By using this approach, a cDNA fragment extending from the SAGE tag toward the 3′ end of the corresponding sequence can be generated. Application of the GLGI technique can solve two critical issues in applying the SAGE technique: one is that a longer fragment corresponding to a SAGE tag, which has no match in databases, can be generated for further studies; the other is that the specific fragment corresponding to a SAGE tag can be identified from multiple sequences that match the same SAGE tag. The development of the GLGI method provides several potential applications. First, it provides a strategy for even wider application of the SAGE technique for quantitative analysis of global gene expression. Second, a combined application of SAGE/GLGI can be used to complete the catalogue of the expressed genes in human and in other eukaryotic species. Third, it can be used to identify the 3′ cDNA sequence from any exon within a gene. It can also be used to confirm the reality of exons predicted by bioinformatic tools in genomic sequences. Fourth, a combined application of SAGE/GLGI can be applied to define the 3′ boundary of expressed genes in the genomic sequences in human and in other eukaryotic genomes.

Appetite of an epiphyte: Quantitative monitoring of bacterial sugar consumption in the phyllosphere

We report here the construction, characterization, and application of a bacterial bioreporter for fructose and sucrose that was designed to monitor the availability of these sugars to microbial colonizers of the phyllosphere. Plasmid pPfruB-gfp[AAV] carries the Escherichia coli fruB promoter upstream from the gfp[AAV] allele that codes for an unstable variant of green fluorescent protein (GFP). In Erwinia herbicola, this plasmid brings about the accumulation of GFP fluorescence in response to both fructose and sucrose. Cells of E. herbicola (pPfruB-gfp[AAV]) were sprayed onto bean plants, recovered from leaves at various time intervals after inoculation, and analyzed individually for GFP content by quantitative analysis of digital microscope images. We observed a positive correlation between single-cell GFP accumulation and ribosomal content as determined by fluorescence in situ hybridization, indicating that foliar growth of E. herbicola occurred at the expense of fructose and/or sucrose. One hour after inoculation, nearly all bioreporter cells appeared to be actively engaged in fructose consumption. This fraction dropped to approximately 11% after 7 h and to ≈1% a day after inoculation. This pattern suggests a highly heterogeneous availability of fructose to individual E. herbicola cells as they colonize the phyllosphere. We estimated that individual cells were exposed to local initial fructose abundances ranging from less than 0.15 pg fructose to more than 4.6 pg.

Quantitative Intercellular Localization of NADH-Dependent Glutamate Synthase Protein in Different Types of Root Cells in Rice Plants1

The quantitative analysis with immunogold-electron microscopy using a single-affinity-purified anti-NADH-glutamate synthase (GOGAT) immunoglobulin G (IgG) as the primary antibody showed that the NADH-GOGAT protein was present in various forms of plastids in the cells of the epidermis and exodermis, in the cortex parenchyma, and in the vascular parenchyma of root tips (<10 mm) of rice (Oryza sativa) seedlings supplied with 1 mm NH4+ for 24 h. The values of the mean immunolabeling density of plastids were almost equal among these different cell types in the roots. However, the number of plastids per individual cell type was not identical, and some parts of the cells in the epidermis and exodermis contained large numbers of plastids that were heavily immunolabeled. Although there was an indication of labeling in the mitochondria using the single-affinity-purified anti-NADH-GOGAT IgG, this was not confirmed when a twice-affinity-purified IgG was used, indicating an exclusively plastidial location of the NADH-GOGAT protein in rice roots. These results, together with previous work from our laboratory (K. Ishiyama, T. Hayakawa, and T. Yamaya [1998] Planta 204: 288–294), suggest that the assimilation of exogeneously supplied NH4+ ions is primarily via the cytosolic glutamine synthetase/plastidial NADH-GOGAT cycle in specific regions of the epidermis and exodermis in rice roots. We also discuss the role of the NADH-GOGAT protein in vascular parenchyma cells.

Initiation sites of protein folding by NMR analysis.

Detailed characterization of denatured states of proteins is necessary to understand the interactions that funnel the large number of possible conformations along fast routes for folding. Nuclear magnetic resonance experiments based on the nuclear Overhauser effect (NOE) detect hydrogen atoms close in space and provide information about local structure. Here we present an NMR procedure that detects almost all sequential NOEs between amide hydrogen atoms (HN-HN NOE), including those in random coil regions in a protein, barnase, in urea solutions. A semi-quantitative analysis of these HN-HN NOEs identified partly structured regions that are in remarkable agreement with those found to form early on the reaction pathway. Our results strongly suggest that the folding of barnase initiates at the first helix and the beta-turn between the third and the fourth strands. This strategy of defining residual structure has also worked for cold-denatured barstar and guanidinium hydrochloride-denatured chymotrypsin inhibitor 2 and so should be generally applicable.

Serve analysis of professional players in beach volleyball

This study has been developed for the European Beach Volleyball Championship in 2005. The video-recorded analysis was held using Sportcode Pro v.8.5.2 software. The aim of study was to determine the types of serve used, depending on the time of the set in which they occur. Quantitative analysis with a sample of 10 players that make up 5 teams with a total of four meetings with a total of 327 serves analyzed. The serves were classified depending on the period which they occurred, the period being 1 (items 1 to 7), period 2 (from point 8 to 14) and period 3 (point 15 to 21). he statistical analysis was conducted using the statistical software SPSS 19, Chi-square test established significant differences between the different types of serve for period 1 and 2 (p<0.05), but no significant differences were established in the period 3 to floating serve and power jump (p>0.05). The results showed a decrease of using a serve with jump power at period 1 (89.7%) compared to the period 3 (27.3%), while the floating and floating serve jump respectively increase at period 1 (6.3% -4%) in the period 3 (23.4% -49.4%).

Quantitative stability of linear infinite inequality systems under block perturbations with applications to convex systems

The original motivation for this paper was to provide an efficient quantitative analysis of convex infinite (or semi-infinite) inequality systems whose decision variables run over general infinite-dimensional (resp. finite-dimensional) Banach spaces and that are indexed by an arbitrary fixed set J. Parameter perturbations on the right-hand side of the inequalities are required to be merely bounded, and thus the natural parameter space is l ∞(J). Our basic strategy consists of linearizing the parameterized convex system via splitting convex inequalities into linear ones by using the Fenchel–Legendre conjugate. This approach yields that arbitrary bounded right-hand side perturbations of the convex system turn on constant-by-blocks perturbations in the linearized system. Based on advanced variational analysis, we derive a precise formula for computing the exact Lipschitzian bound of the feasible solution map of block-perturbed linear systems, which involves only the system’s data, and then show that this exact bound agrees with the coderivative norm of the aforementioned mapping. In this way we extend to the convex setting the results of Cánovas et al. (SIAM J. Optim. 20, 1504–1526, 2009) developed for arbitrary perturbations with no block structure in the linear framework under the boundedness assumption on the system’s coefficients. The latter boundedness assumption is removed in this paper when the decision space is reflexive. The last section provides the aimed application to the convex case.

Involvement and detachment in David Cameron and Ed Miliband’s political discourse. An analysis of deictic shields on the pronoun ‘I’ in UK 2015 General Election news interviews

The goal of my study is to investigate the relationship between selected deictic shields on the pronoun ‘I’ and the involvement/detachment dichotomy in a sample of television news interviews. I focus on the use of personal pronouns in political discourse. Drawing upon Caffi’s (2007) classification of mitigating devices into bushes, hedges and shields, I focus on deictic shields on the pronoun ‘I’: I examine the way a selection of ‘I’-related deictic shields is employed in a collection of news interviews broadcast during the electoral campaign prior to the UK 2015 General Election. My purpose is to uncover the frequencies of each of the linguistic items selected and the pragmatic functions of those linguistic items in the involvement/detachment dichotomy. The research is structured as follows. Chapter 1 provides an account of previous studies on the three main areas of research: speech event analysis, institutional interaction and the news interview, and the UK 2015 General Election television programmes. Chapter 2 is centred on the involvement/detachment dichotomy: I provide an overview of nonlinguistic and linguistic features of involvement and detachment at all levels of sentence structure. Chapter 3 contains a detailed account of the data collection and data analysis process. Chapter 4 provides an accurate description of results in three steps: quantitative analysis, qualitative analysis and discussion of the pragmatic functions of the selected linguistic features of involvement and detachment. Chapter 5 includes a brief summary of the investigation, reviews the main findings, and indicates limitations of the study and possible inputs for further research. The results of the analysis confirm that, while some of the linguistic items examined point toward involvement, others have a detaching effect. I therefore conclude that deictic shields on the pronoun ‘I’ permit the realisation of the involvement/detachment dichotomy in the speech genre of the news interview.

Zebrafish Caudal Fin Angiogenesis Assay-Advanced Quantitative Assessment Including 3-Way Correlative Microscopy.

BACKGROUND Researchers evaluating angiomodulating compounds as a part of scientific projects or pre-clinical studies are often confronted with limitations of applied animal models. The rough and insufficient early-stage compound assessment without reliable quantification of the vascular response counts, at least partially, to the low transition rate to clinics. OBJECTIVE To establish an advanced, rapid and cost-effective angiogenesis assay for the precise and sensitive assessment of angiomodulating compounds using zebrafish caudal fin regeneration. It should provide information regarding the angiogenic mechanisms involved and should include qualitative and quantitative data of drug effects in a non-biased and time-efficient way. APPROACH & RESULTS Basic vascular parameters (total regenerated area, vascular projection area, contour length, vessel area density) were extracted from in vivo fluorescence microscopy images using a stereological approach. Skeletonization of the vasculature by our custom-made software Skelios provided additional parameters including "graph energy" and "distance to farthest node". The latter gave important insights into the complexity, connectivity and maturation status of the regenerating vascular network. The employment of a reference point (vascular parameters prior amputation) is unique for the model and crucial for a proper assessment. Additionally, the assay provides exceptional possibilities for correlative microscopy by combining in vivo-imaging and morphological investigation of the area of interest. The 3-way correlative microscopy links the dynamic changes in vivo with their structural substrate at the subcellular level. CONCLUSIONS The improved zebrafish fin regeneration model with advanced quantitative analysis and optional 3-way correlative morphology is a promising in vivo angiogenesis assay, well-suitable for basic research and preclinical investigations.

Quantitative photoluminescence of broad band absorbing melanins: a procedure to correct for inner filter and re-absorption effects

We report methods for correcting the photoluminescence emission and excitation spectra of highly absorbing samples for re-absorption and inner filter effects. We derive the general form of the correction, and investigate various methods for determining the parameters. Additionally, the correction methods are tested with highly absorbing fluorescein and melanin (broadband absorption) solutions; the expected linear relationships between absorption and emission are recovered upon application of the correction, indicating that the methods are valid. These procedures allow accurate quantitative analysis of the emission of low quantum yield samples (such as melanin) at concentrations where absorption is significant. (c) 2004 Elsevier B.V. All rights reserved.

Folate: Methods of analysis

Folates and its derivatives occur as polyglutamates in nature. The multiplicity of forms and the generally low levels in foods makes quantitative analysis of folate a difficult task. The assay of folates from foods generally involves three steps: liberation of folates from the cellular matrix; deconjugation from the polyglutamate to the mono and di-glutamate forms; and the detection of the biological activity or chemical concentration of the resulting folates. The detection methods used are the microbiological assay relying on the turbidimetric bacterial growth of Lactobacillus rhamnosus which is by far the most commonly used method; the HPLC and LC/MS techniques and bio-specific procedures. This review attempts to describe the methods along with the merits and demerits of using each of these methods.

Genome-wide analysis of mammalian promoter architecture and evolution

Mammalian promoters can be separated into two classes, conserved TATA box-enriched promoters, which initiate at a welldefined site, and more plastic, broad and evolvable CpG-rich promoters. We have sequenced tags corresponding to several hundred thousand transcription start sites (TSSs) in the mouse and human genomes, allowing precise analysis of the sequence architecture and evolution of distinct promoter classes. Different tissues and families of genes differentially use distinct types of promoters. Our tagging methods allow quantitative analysis of promoter usage in different tissues and show that differentially regulated alternative TSSs are a common feature in protein-coding genes and commonly generate alternative N termini. Among the TSSs, we identified new start sites associated with the majority of exons and with 3' UTRs. These data permit genome-scale identification of tissue-specific promoters and analysis of the cis-acting elements associated with them.

Conformational analysis and dielectric relaxation of polycarbosilanes and related materials

The conformational characteristics of poly(dimethylsilmethylene), poly(dimethylsilethene), poly(dimethylsilethane) and a related material, poly(2,2,5,5-tetramethyl-1-oxa-2,5-disilapentane), have been investigated using the method of molecular mechanics. In this method, a quantitative analysis of the factors affecting the nature and magnitude of the bond rotation potentials governing their conformational behaviour has been undertaken. Along with their structural data, the results obtained were employed to calculate a variety of conformationally-dependent properties for these polymers, including the characteristic ratio, the dipole moment ratio and the mean-square radius of gyration. In addition, the dielectric relaxation behaviour of two samples of poly(2,2,5,5-tetramethyl-1-oxa-2,5-disilapentane) with molar masses Mw = 28000 and Mw = 46000 respectively, have been studied as a function of temperature (179K-205K) and frequency (100-105Hz). Activation energies for the α-relaxation process and Davidson-Cole empirical distribution factors have been calculated.

Development of an on-demand beam pulsing system for Pixe-analysis of thick targets

This work concerns the developnent of a proton irduced X-ray emission (PIXE) analysis system and a multi-sample scattering chamber facility. The characteristics of the beam pulsing system and its counting rate capabilities were evaluated by observing the ion-induced X-ray emission from pure thick copper targets, with and without beam pulsing operation. The characteristic X-rays were detected with a high resolution Si(Li) detector coupled to a rrulti-channel analyser. The removal of the pile-up continuum by the use of the on-demand beam pulsing is clearly demonstrated in this work. This new on-demand pu1sirg system with its counting rate capability of 25, 18 and 10 kPPS corresponding to 2, 4 am 8 usec main amplifier time constant respectively enables thick targets to be analysed more readily. Reproducibility tests of the on-demard beam pulsing system operation were checked by repeated measurements of the system throughput curves, with and without beam pulsing. The reproducibility of the analysis performed using this system was also checked by repeated measurements of the intensity ratios from a number of standard binary alloys during the experimental work. A computer programme has been developed to evaluate the calculations of the X-ray yields from thick targets bornbarded by protons, taking into account the secondary X-ray yield production due to characteristic X-ray fluorescence from an element energetically higher than the absorption edge energy of the other element present in the target. This effect was studied on metallic binary alloys such as Fe/Ni and Cr/Fe. The quantitative analysis of Fe/Ni and Cr/Fe alloy samples to determine their elemental composition taking into account the enhancement has been demonstrated in this work. Furthermore, the usefulness of the Rutherford backscattering (R.B.S.) technique to obtain the depth profiles of the elements in the upper micron of the sample is discussed.