An adaptive multiscale method for density-driven instabilities

Accurate modeling of flow instabilities requires computational tools able to deal with several interacting scales, from the scale at which fingers are triggered up to the scale at which their effects need to be described. The Multiscale Finite Volume (MsFV) method offers a framework to couple fine-and coarse-scale features by solving a set of localized problems which are used both to define a coarse-scale problem and to reconstruct the fine-scale details of the flow. The MsFV method can be seen as an upscaling-downscaling technique, which is computationally more efficient than standard discretization schemes and more accurate than traditional upscaling techniques. We show that, although the method has proven accurate in modeling density-driven flow under stable conditions, the accuracy of the MsFV method deteriorates in case of unstable flow and an iterative scheme is required to control the localization error. To avoid large computational overhead due to the iterative scheme, we suggest several adaptive strategies both for flow and transport. In particular, the concentration gradient is used to identify a front region where instabilities are triggered and an accurate (iteratively improved) solution is required. Outside the front region the problem is upscaled and both flow and transport are solved only at the coarse scale. This adaptive strategy leads to very accurate solutions at roughly the same computational cost as the non-iterative MsFV method. In many circumstances, however, an accurate description of flow instabilities requires a refinement of the computational grid rather than a coarsening. For these problems, we propose a modified iterative MsFV, which can be used as downscaling method (DMsFV). Compared to other grid refinement techniques the DMsFV clearly separates the computational domain into refined and non-refined regions, which can be treated separately and matched later. This gives great flexibility to employ different physical descriptions in different regions, where different equations could be solved, offering an excellent framework to construct hybrid methods.

Discriminant ECOC: a heuristic method for application dependent design of error correcting output codes

We present a heuristic method for learning error correcting output codes matrices based on a hierarchical partition of the class space that maximizes a discriminative criterion. To achieve this goal, the optimal codeword separation is sacrificed in favor of a maximum class discrimination in the partitions. The creation of the hierarchical partition set is performed using a binary tree. As a result, a compact matrix with high discrimination power is obtained. Our method is validated using the UCI database and applied to a real problem, the classification of traffic sign images.

Extrapolation theory for the real interpolation method

We develop an abstract extrapolation theory for the real interpolation method that covers and improves the most recent versions of the celebrated theorems of Yano and Zygmund. As a consequence of our method, we give new endpoint estimates of the embedding Sobolev theorem for an arbitrary domain Omega

High-throughput SELEX SAGE method for quantitative modeling of transcription-factor binding sites.

The ability to determine the location and relative strength of all transcription-factor binding sites in a genome is important both for a comprehensive understanding of gene regulation and for effective promoter engineering in biotechnological applications. Here we present a bioinformatically driven experimental method to accurately define the DNA-binding sequence specificity of transcription factors. A generalized profile was used as a predictive quantitative model for binding sites, and its parameters were estimated from in vitro-selected ligands using standard hidden Markov model training algorithms. Computer simulations showed that several thousand low- to medium-affinity sequences are required to generate a profile of desired accuracy. To produce data on this scale, we applied high-throughput genomics methods to the biochemical problem addressed here. A method combining systematic evolution of ligands by exponential enrichment (SELEX) and serial analysis of gene expression (SAGE) protocols was coupled to an automated quality-controlled sequence extraction procedure based on Phred quality scores. This allowed the sequencing of a database of more than 10,000 potential DNA ligands for the CTF/NFI transcription factor. The resulting binding-site model defines the sequence specificity of this protein with a high degree of accuracy not achieved earlier and thereby makes it possible to identify previously unknown regulatory sequences in genomic DNA. A covariance analysis of the selected sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism.

Critical roles of PPAR beta/delta in keratinocyte response to inflammation.

The immediate response to skin injury is the release of inflammatory signals. It is shown here, by use of cultures of primary keratinocytes from wild-type and PPAR beta/delta(-/-) mice, that such signals including TNF-alpha and IFN-gamma, induce keratinocyte differentiation. This cytokine-dependent cell differentiation pathway requires up-regulation of the PPAR beta/delta gene via the stress-associated kinase cascade, which targets an AP-1 site in the PPAR beta/delta promoter. In addition, the pro-inflammatory cytokines also initiate the production of endogenous PPAR beta/delta ligands, which are essential for PPAR beta/delta activation and action. Activated PPAR beta/delta regulates the expression of genes associated with apoptosis resulting in an increased resistance of cultured keratinocytes to cell death. This effect is also observed in vivo during wound healing after an injury, as shown in dorsal skin of PPAR beta/delta(+/+) and PPAR beta/delta(+/-) mice.

Semiquantitative RT-PCR measurement of gene expression in rat tissues including a correction for varying cell size and number

Background: Current methodology of gene expression analysis limits the possibilities of comparison between cells/tissues of organs in which cell size and/or number changes as a consequence of the study (e.g. starvation). A method relating the abundance of specific mRNA copies per cell may allow direct comparison or different organs and/or changing physiological conditions. Methods: With a number of selected genes, we analysed the relationship of the number of bases and the fluorescence recorded at a present level using cDNA standards. A lineal relationship was found between the final number of bases and the length of the transcript. The constants of this equation and those of the relationship between fluorescence and number of bases in cDNA were determined and a general equation linking the length of the transcript and the initial number of copies of mRNA was deduced for a given pre-established fluorescence setting. This allowed the calculation of the concentration of the corresponding mRNAs per g of tissue. The inclusion of tissue RNA and the DNA content per cell, allowed the calculation of the mRNA copies per cell. Results: The application of this procedure to six genes: Arbp, cyclophilin, ChREBP, T4 deiodinase 2, acetyl-CoA carboxylase 1 and IRS-1, in liver and retroperitoneal adipose tissue of food-restricted rats allowed precise measures of their changes irrespective of the shrinking of the tissue, the loss of cells or changes in cell size, factors that deeply complicate the comparison between changing tissue conditions. The percentage results obtained with the present methods were essentially the same obtained with the delta-delta procedure and with individual cDNA standard curve quantitative RT-PCR estimation. Conclusion: The method presented allows the comparison (i.e. as copies of mRNA per cell) between different genes and tissues, establishing the degree of abundance of the different molecular species tested.

Role of β-catenin and γ-catenin in haematopoiesis : Generation of floxed Delta like 4 gene targeted mice

Résumé La voie de signalisation de Wnt est extrêmement conservée au cours de l'évolution. Les protéines Wnt sont des molécules sécrétées qui se lient à la famille de récepteurs Frizzled. Cette interaction mène à la stabilisation de la protéine β-caténine, qui va s'accumuler dans le cytoplasme puis migrer dans le noyau où elle peut s'hétérodimériser avec les facteurs de transcription de la famille TCF/LEF. Il a été démontré que cette voie de signalisation joue un rôle important durant la lymphopoïèse et de récents résultats suggèrent un rôle clé de cette voie dans le renouvellement des Cellules Souches Hématopoïétique (CSH). Des études se basant sur un système de surexpression de protéines montrent clairement que la voie Wnt peut influencer l'hématopoïèse. Cependant, le rôle de la protéine β-caténine dans le système hématopoïétique n'a jamais été testé directement. Ce projet de thèse se propose d'étudier la fonction de la protéine β-caténine par sa délétion inductible via le système Cre-loxP. De façon surprenante, nous avons pu démontrer que les progéniteurs de la moelle osseuse, déficients en β-caténine, ne montrent aucune altération dans leur capacité à s'auto-renouveler et/ou à reconstituer toutes les lignées hématopoïétiques (myéloïde, érythroïde et lymphoïde) dans les souris-chimères. De plus, le développement, la survie des thymocytes ainsi que la prolifération des cellules T périphériques induite par un antigène, sont indépendants de β-caténine. Ces résultats suggèrent soit que la protéine β-caténine ne joue pas un rôle primordial dans le système hématopoiétique, soit que son absence pourrait être compensée par une autre protéine. Un candidat privilégié susceptible de se substituer à β-caténine, serait plakoglobine, aussi connu sous le nom de γ-caténine. En effet, ces deux protéines partagent de multiples caractéristiques structurelles. Afin de démontrer que la protéine γ-caténine peut compenser l'absence de β-caténine, nous avons généré des souris dans lesquelles, le système hématopoïétique est déficient pour ces deux protéines. Cette déficience combinée de β- caténine et γ-caténine ne perturbe pas la capacité des Cellules Souche Hématopoïétique-Long Terme (CSH-LT) de se renouveler, par contre elle agit sur un progéniteur précoce déjà différencié de la moelle osseuse. Ces résultats mettent en évidence que la protéine γ-caténine est capable de compenser l'absence de protéine β-caténine dans le système hématopoïétique. Par conséquent, ce travail contribue à une meilleure connaissance de la cascade Wnt dans l'hématopoïèse. Summary The canonical Wnt signal transduction pathway is a developmentally highly conserved. Wnts are secreted molecules which bind to the family of Frizzled receptors in a complex with the low density lipoprotein receptor related protein (LRP-5/6). This initial activation step leads to the stabilization and accumulation of β-catenin, first in the cytoplasm and subsequently in the nucleus where it forms heterodimers with TCF/LEF transcription factor family members. Wnt signalling has been shown to be important during early lymphopoiesis and has more recently, been suggested to be a key player in self-renewal of haematopoietic stem cells (HSCs). Although mostly gain of function studies indicate that components of the Wnt signalling pathway can influence the haematopoietic system, the role of β-catenin has never been directly investigated. The aim of this thesis project is to investigate the putatively critical role of β-catenin in vivo using the Cre-loxP mediated conditional loss of function approach. Surprisingly, β-catenin deficient bone marrow (BM) progenitors arc not impaired in their ability to self-renew and/or to reconstitute all haematopoietic lineages (myeloid, erythroid and lymphoid) in both mixed and straight bone marrow chimeras. In addition, both thymocyte development and survival, and antigen-induced proliferation of peripheral T cells are β- catenin independent. Our results do not necessarily exclude the possibility of an important function for β-catenin mediated Wnt signalling in the haematopoietic system, it rather raises the question that β-catenin is compensated for by another protein. A prime candidate that may take over the function of β-catenin in its absence, is the close relative plakoglobin, also know as γ-catenin. This protein shares multiple structural features with β-catenin. In order to investigate whether γ-catenin can compensate for the loss of β-catenin we have generated mice in which the haematopoietic compartment is deficient for both proteins. Combined deficiency of β-catenin and γ-catenin does not perturb Long Term-Haematopoietic Stem Cells (LT-HSC) self renewal, but affects an already lineage committed progenitor population within the BM. Our results demonstrate that y-catenin can indeed compensate for the loss of β-catenin within the haematopoietie system.

A novel method of dynamic culture surface expansion improves mesenchymal stem cell proliferation and phenotype.

Repeated passaging in conventional cell culture reduces pluripotency and proliferation capacity of human mesenchymal stem cells (MSC). We introduce an innovative cell culture method whereby the culture surface is dynamically enlarged during cell proliferation. This approach maintains constantly high cell density while preventing contact inhibition of growth. A highly elastic culture surface was enlarged in steps of 5% over the course of a 20-day culture period to 800% of the initial surface area. Nine weeks of dynamic expansion culture produced 10-fold more MSC compared with conventional culture, with one-third the number of trypsin passages. After 9 weeks, MSC continued to proliferate under dynamic expansion but ceased to grow in conventional culture. Dynamic expansion culture fully retained the multipotent character of MSC, which could be induced to differentiate into adipogenic, chondrogenic, osteogenic, and myogenic lineages. Development of an undesired fibrogenic myofibroblast phenotype was suppressed. Hence, our novel method can rapidly provide the high number of autologous, multipotent, and nonfibrogenic MSC needed for successful regenerative medicine.

A multilevel multiscale finite volume method

The Multiscale Finite Volume (MsFV) method has been developed to efficiently solve reservoir-scale problems while conserving fine-scale details. The method employs two grid levels: a fine grid and a coarse grid. The latter is used to calculate a coarse solution to the original problem, which is interpolated to the fine mesh. The coarse system is constructed from the fine-scale problem using restriction and prolongation operators that are obtained by introducing appropriate localization assumptions. Through a successive reconstruction step, the MsFV method is able to provide an approximate, but fully conservative fine-scale velocity field. For very large problems (e.g. one billion cell model), a two-level algorithm can remain computational expensive. Depending on the upscaling factor, the computational expense comes either from the costs associated with the solution of the coarse problem or from the construction of the local interpolators (basis functions). To ensure numerical efficiency in the former case, the MsFV concept can be reapplied to the coarse problem, leading to a new, coarser level of discretization. One challenge in the use of a multilevel MsFV technique is to find an efficient reconstruction step to obtain a conservative fine-scale velocity field. In this work, we introduce a three-level Multiscale Finite Volume method (MlMsFV) and give a detailed description of the reconstruction step. Complexity analyses of the original MsFV method and the new MlMsFV method are discussed, and their performances in terms of accuracy and efficiency are compared.

Development of a two-step screening ESI-TOF-MS method for rapid determination of significant stress-induced metabolome modifications in plant leaf extracts: the wound response in Arabidopsis thaliana as a case study.

To study the stress-induced effects caused by wounding under a new perspective, a metabolomic strategy based on HPLC-MS has been devised for the model plant Arabidopsis thaliana. To detect induced metabolites and precisely localise these compounds among the numerous constitutive metabolites, HPLC-MS analyses were performed in a two-step strategy. In a first step, rapid direct TOF-MS measurements of the crude leaf extract were performed with a ballistic gradient on a short LC-column. The HPLC-MS data were investigated by multivariate analysis as total mass spectra (TMS). Principal components analysis (PCA) and hierarchical cluster analysis (HCA) on principal coordinates were combined for data treatment. PCA and HCA demonstrated a clear clustering of plant specimens selecting the highest discriminating ions given by the complete data analysis, leading to the specific detection of discrete-induced ions (m/z values). Furthermore, pool constitution with plants of homogeneous behaviour was achieved for confirmatory analysis. In this second step, long high-resolution LC profilings on an UPLC-TOF-MS system were used on pooled samples. This allowed to precisely localise the putative biological marker induced by wounding and by specific extraction of accurate m/z values detected in the screening procedure with the TMS spectra.

The systematic profiling of false identity documents: method validation and performance evaluation using seizures known to originate from common and different sources

False identity documents constitute a potential powerful source of forensic intelligence because they are essential elements of transnational crime and provide cover for organized crime. In previous work, a systematic profiling method using false documents' visual features has been built within a forensic intelligence model. In the current study, the comparison process and metrics lying at the heart of this profiling method are described and evaluated. This evaluation takes advantage of 347 false identity documents of four different types seized in two countries whose sources were known to be common or different (following police investigations and dismantling of counterfeit factories). Intra-source and inter-sources variations were evaluated through the computation of more than 7500 similarity scores. The profiling method could thus be validated and its performance assessed using two complementary approaches to measuring type I and type II error rates: a binary classification and the computation of likelihood ratios. Very low error rates were measured across the four document types, demonstrating the validity and robustness of the method to link documents to a common source or to differentiate them. These results pave the way for an operational implementation of a systematic profiling process integrated in a developed forensic intelligence model.