Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data.

MOTIVATION: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts. Results: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays.

ON THE ROLE OF PERIOD-2 IN THE CIRCADIAN AND HOMEOSTATIC REGULATION OF SLEEP

Humans spend one third of their life sleeping, then we could raise the basic question: Why do we sleep? Despite the fact that we still don't fully understand its function, we made much progress in understanding at different levels how sleep is regulated. One model suggests that sleep is regulated by two processes: a homeostatic process that tracks the need for sleep and by a circadian rhythm that determines the preferred time-of-day sleep occurs. At the molecular level circadian rhythms are a property of interlocking transcriptional regula-tors referred to as clock genes. The heterodimeric transcription factors BMAL1::CLOCK/NPAS2 drive the transcription of many target genes including the clock genes Cryptochome1 (Cry1), Cry2, Period1 (Per1), and Per2. The encoded CRY/PER proteins are transcriptional inhibitors of BMAL1::CLOCK/NPAS2 thereby providing negative feedback to their own transcription. These genes seem, however, also involved in sleep homeostasis because the brain expression of clock genes, es-pecially that of Per2, increase as a function of time-spent-awake and because mice lacking clock genes display altered sleep homeostasis. The aim of first part of my doctoral work has been to advance our understanding the link that exists between sleep homeostasis and circadian rhythms investigating a possible mechanism by which sleep deprivation could alter clock gene expression by quantifying DNA-binding of the core-clock genes BMAL1, CLOCK and NPAS2 to their target chromatin loci including the E-box enhancers of the Per2 promoter. We made use of chromatin immunoprecipitation (ChIP) and quantitative poly-merase chain reaction (qPCR) to show that DNA-binding of CLOCK and BMAL1 to their target genes changes as a function of time-of-day in both liver and cerebral cortex. We then performed a 6h sleep deprivation (SD) and observed a significant decrease in DNA-binding of CLOCK and BMAL1 to Dbp. This is consistent with a decrease in Dbp mRNA levels after SD. The DNA-binding of NPAS2 and BMAL1 to Per2 was similarly decreased following SD. However, SD has been previously shown to in-crease Per2 expression in the cortex which seems paradoxical. Our results demonstrate that sleep-wake history can affect the molecular clock machinery directly at the level of the chromatin thereby altering the cortical expression of Dbp and Per2, and likely other targets. However, the precise dy-namic relationship between DNA-binding and mRNA expression, especially for Per2, remains elusive. The second aim of my doctoral work has been to perform an in depth characterization of cir-cadian rhythmicity, sleep architecture, analyze the response to SD in full null-Per2 knock-out (Per2-/-) mice, and Per1-/- mice, as well as their double knock-out offspring (Per1,2-/-) and littermate wildtype (Wt) mice. The techniques used include locomotor activity recording by passive infrared (PIR) sen-sors, EEG/EMG surgery, recording, and analysis, and cerebral cortex extraction and quantification of mRNA levels by qPCR. Under standard LD12:12 conditions, we found that wakefulness onset, as well as the time courses of clock gene expression in the brain and corticosterone plasma levels were ad-vanced by about 2h in Per2-/- mice compared to Wt mice. When released under constant dark condi-tions almost all Per2-/- mice (97%) became arrhythmic immediately. From these observations, we conclude that while Per2-/- mice seem to be able to anticipate dark onset, this does not result from a self-sustained circadian clock. Our results suggest instead that the earlier onset of activity results from a labile, not-self sustained 22h rhythm linked to light onset suggesting the existence of a light-driven rhythm. Analyses of sleep under LD12:12 conditions revealed that in both Per2-/- and Per1,2-/- mice the same sleep phenotypes are observed compared to Wt mice: increased NREM sleep frag-mentation and inability to adequately compensate the loss of NREM sleep. That suggests a possible role of PER2 in sleep consolidation and recovery.

Residual thermomechanical stresses in thinned-chip assemblies

A new technology for the three-dimensional (3-D) stacking of very thin chips on a substrate is currently under development within the ultrathin chip stacking (UTCS) Esprit Project 24910. In this work, we present the first-level UTCS structure and the analysis of the thermomechanical stresses produced by the manufacturing process. Chips are thinned up to 10 or 15 m. We discuss potentially critical points at the edges of the chips, the suppression of delamination problems of the peripheral dielectric matrix and produce a comparative study of several technological choices for the design of metallic interconnect structures. The purpose of these calculations is to give inputs for the definition of design rules for this technology. We have therefore undertaken a programme that analyzes the influence of sundry design parameters and alternative development options. Numerical analyses are based on the finite element method.

Thermal modeling and management in ultrathin chip stack technology

This paper presents a thermal modeling for power management of a new three-dimensional (3-D) thinned dies stacking process. Besides the high concentration of power dissipating sources, which is the direct consequence of the very interesting integration efficiency increase, this new ultra-compact packaging technology can suffer of the poor thermal conductivity (about 700 times smaller than silicon one) of the benzocyclobutene (BCB) used as both adhesive and planarization layers in each level of the stack. Thermal simulation was conducted using three-dimensional (3-D) FEM tool to analyze the specific behaviors in such stacked structure and to optimize the design rules. This study first describes the heat transfer limitation through the vertical path by examining particularly the case of the high dissipating sources under small area. First results of characterization in transient regime by means of dedicated test device mounted in single level structure are presented. For the design optimization, the thermal draining capabilities of a copper grid or full copper plate embedded in the intermediate layer of stacked structure are evaluated as a function of the technological parameters and the physical properties. It is shown an interest for the transverse heat extraction under the buffer devices dissipating most the power and generally localized in the peripheral zone, and for the temperature uniformization, by heat spreading mechanism, in the localized regions where the attachment of the thin die is altered. Finally, all conclusions of this analysis are used for the quantitative projections of the thermal performance of a first demonstrator based on a three-levels stacking structure for space application.

Dynamic compact thermal models with multiple power sources: application to an ultrathin chip stacking technology

Whereas numerical modeling using finite-element methods (FEM) can provide transient temperature distribution in the component with enough accuracy, it is of the most importance the development of compact dynamic thermal models that can be used for electrothermal simulation. While in most cases single power sources are considered, here we focus on the simultaneous presence of multiple sources. The thermal model will be in the form of a thermal impedance matrix containing the thermal impedance transfer functions between two arbitrary ports. Eachindividual transfer function element ( ) is obtained from the analysis of the thermal temperature transient at node ¿ ¿ after a power step at node ¿ .¿ Different options for multiexponential transient analysis are detailed and compared. Among the options explored, small thermal models can be obtained by constrained nonlinear least squares (NLSQ) methods if the order is selected properly using validation signals. The methods are applied to the extraction of dynamic compact thermal models for a new ultrathin chip stack technology (UTCS).

Phytogeographical origin of Madeiran common beans based on phaseolin patterns

The objective of this work was to determine the geographic origin of the Madeiran common bean (Phaseolus vulgaris) gene pool. Phaseolin patterns of 50 accessions representing the diversity of common bean collected in Madeira, Portugal, and conserved in the ISOPlexis Germplasm Bank, were analysed using the Experion automated electrophoresis system, based on lab-on-a-chip technology. Five common bean standard varieties with typical phaseolin patterns were used to determine the phytogeographical origin of the Madeiran common bean accessions. Ninety two percent of the accessions exhibited a phaseolin pattern consistent with the one of common bean types belonging to the Andean gene pool, while the origin of the remaining 8% of the accessions was indistinguishable. The application of a similarity coefficient of 85%, based on Pearson correlations, increases the number of accessions with uncertain pattern. The analytical approach used permitted the determination of the origin of the common bean gene pool, which is Andean in 98% of the cases, and clustering of the observed variability among the Madeiran common beans.

On the Condition Monitoring of Induction machines

The research of condition monitoring of electric motors has been wide for several decades. The research and development at universities and in industry has provided means for the predictive condition monitoring. Many different devices and systems are developed and are widely used in industry, transportation and in civil engineering. In addition, many methods are developed and reported in scientific arenas in order to improve existing methods for the automatic analysis of faults. The methods, however, are not widely used as a part of condition monitoring systems. The main reasons are, firstly, that many methods are presented in scientific papers but their performance in different conditions is not evaluated, secondly, the methods include parameters that are so case specific that the implementation of a systemusing such methods would be far from straightforward. In this thesis, some of these methods are evaluated theoretically and tested with simulations and with a drive in a laboratory. A new automatic analysis method for the bearing fault detection is introduced. In the first part of this work the generation of the bearing fault originating signal is explained and its influence into the stator current is concerned with qualitative and quantitative estimation. The verification of the feasibility of the stator current measurement as a bearing fault indicatoris experimentally tested with the running 15 kW induction motor. The second part of this work concentrates on the bearing fault analysis using the vibration measurement signal. The performance of the micromachined silicon accelerometer chip in conjunction with the envelope spectrum analysis of the cyclic bearing faultis experimentally tested. Furthermore, different methods for the creation of feature extractors for the bearing fault classification are researched and an automatic fault classifier using multivariate statistical discrimination and fuzzy logic is introduced. It is often important that the on-line condition monitoring system is integrated with the industrial communications infrastructure. Two types of a sensor solutions are tested in the thesis: the first one is a sensor withcalculation capacity for example for the production of the envelope spectra; the other one can collect the measurement data in memory and another device can read the data via field bus. The data communications requirements highly depend onthe type of the sensor solution selected. If the data is already analysed in the sensor the data communications are needed only for the results but in the other case, all measurement data need to be transferred. The complexity of the classification method can be great if the data is analysed at the management level computer, but if the analysis is made in sensor itself, the analyses must be simple due to the restricted calculation and memory capacity.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study.

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHRadjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHRadjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHRadjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

On the numerical simulation of machining processes

Numerical simulation of machining processes can be traced back to the early seventies when finite element models for continuous chip formation were proposed. The advent of fast computers and development of new techniques to model large plastic deformations have favoured machining simulation. Relevant aspects of finite element simulation of machining processes are discussed in this paper, such as solution methods, material models, thermo-mechanical coupling, friction models, chip separation and breakage strategies and meshing/re-meshing strategies.

Effects of the partial removal of wood hemicelluloses on the properties of kraft pulp

The objective of this work was to study the effects of partial removal of wood hemicelluloses on the properties of kraft pulp.The work was conducted by extracting hemicelluloses (1) by a softwood chip pretreatment process prior to kraft pulping, (2) by alkaline extraction from bleached birch kraft pulp, and (3) by enzymatic treatment, xylanase treatment in particular, of bleached birch kraft pulp. The qualitative and quantitative changes in fibers and paper properties were evaluated. In addition, the applicability of the extraction concepts and hemicellulose-extracted birch kraft pulp as a raw material in papermaking was evaluated in a pilot-scale papermaking environment. The results showed that each examined hemicellulose extraction method has its characteristic effects on fiber properties, seen as differences in both the physical and chemical nature of the fibers. A prehydrolysis process prior to the kraft pulping process offered reductions in cooking time, bleaching chemical consumption and produced fibers with low hemicellulose content that are more susceptible to mechanically induced damages and dislocations. Softwood chip pretreatment for hemicellulose recovery prior to cooking, whether acidic or alkaline, had an impact on the physical properties of the non-refined and refined pulp. In addition, all the pretreated pulps exhibited slower beating response than the unhydrolyzed reference pulp. Both alkaline extraction and enzymatic (xylanase) treatment of bleached birch kraft pulp fibers indicated very selective hemicellulose removal, particularly xylan removal. Furthermore, these two hemicellulose-extracted birch kraft pulps were utilized in a pilot-scale papermaking environment in order to evaluate the upscalability of the extraction concepts. Investigations made using pilot paper machine trials revealed that some amount of alkalineextracted birch kraft pulp, with a 24.9% reduction in the total amount of xylan, could be used in the papermaking stock as a mixture with non-extracted pulp when producing 75 g/m2 paper. For xylanase-treated fibers there were no reductions in the mechanical properties of the 180 g/m2 paper produced compared to paper made from the control pulp, although there was a 14.2% reduction in the total amount of xylan in the xylanase-treated pulp compared to the control birch kraft pulp. This work emphasized the importance of the hemicellulose extraction method in providing new solutions to create functional fibers and in providing a valuable hemicellulose co-product stream. The hemicellulose removal concept therefore plays an important role in the integrated forest biorefinery scenario, where the target is to the co-production of hemicellulose-extracted pulp and hemicellulose-based chemicals or fuels.

The euchromatic and heterochromatic landscapes are shaped by antagonizing effects of transcription on H2A.Z deposition

A role for variant histone H2A.Z in gene expression is now well established but little is known about the mechanisms by which it operates. Using a combination of ChIP-chip, knockdown and expression profiling experiments, we show that upon gene induction, human H2A.Z associates with gene promoters and helps in recruiting the transcriptional machinery. Surprisingly, we also found that H2A.Z is randomly incorporated in the genome at low levels and that active transcription antagonizes this incorporation in transcribed regions. After cessation of transcription, random H2A.Z quickly reappears on genes, demonstrating that this incorporation utilizes an active mechanism. Within facultative heterochromatin, we observe a hyper accumulation of the variant histone, which might be due to the lack of transcription in these regions. These results show how chromatin structure and transcription can antagonize each other, therefore shaping chromatin and controlling gene expression.

A mechanism for co-transcriptional recruitment of mRNA localization factor on nascent mRNAs in budding yeast

Le transport et la localisation des ARN messagers permettent de réguler l’expression spatiale et temporelle de facteurs spécifiques impliqués dans la détermination du destin cellulaire, la plasticité synaptique, la polarité cellulaire et la division asymétrique des cellules. Chez S.cerevisiæ, plus de trente transcrits sont transportés activement vers le bourgeon cellulaire. Parmi ces transcrits, l’ARNm ASH1 (asymetric synthesis of HO) est localisé à l’extrémité du bourgeon pendant l’anaphase. Ce processus va entrainer une localisation asymétrique de la protéine Ash1p, qui sera importée uniquement dans le noyau de la cellule fille, où elle entraine le changement de type sexuel. La localisation asymétrique de l’ARNm ASH1, et donc de Ash1p, implique la présence de différents facteurs de localisation. Parmi ces facteurs, les protéines She (She1p/Myo4p, She2p et She3p) et les répresseurs traductionnels (Puf6p, Loc1p et Khd1p) participent à ce mécanisme. La protéine navette She2p est capable de lier l’ARNm ASH1 et va entrainer le ciblage de cet ARNm vers l’extrémité du bourgeon en recrutant le complexe She3p-Myo4p. Des répresseurs traductionnels régulent la traduction de cet ARNm et évitent l’expression ectopique de la protéine Ash1p pendant son transport. Alors que la fonction cytoplasmique de She2p sur la localisation des ARNm est connue, sa fonction nucléaire est encore inconnue. Nous avons montré que She2p contient une séquence de localisation nucléaire non classique qui est essentielle à son import nucléaire médié par l’importine α (Srp1p). L’exclusion de She2p du noyau par mutation de son NLS empêche la liaison de Loc1p et Puf6p sur l’ARNm ASH1, entrainant un défaut de localisation de l’ARNm et de la protéine. Pour étudier plus en détail l’assemblage de la machinerie de localisation des ARNm dans le noyau, nous avons utilisé des techniques d’immunoprécipitation de chromatine afin de suivre le recrutement des facteurs de localisation et des répresseurs traductionnels sur les ARNm naissants. Nous avons montré que She2p est recruté sur le gène ASH1 pendant sa transcription, via son interaction avec l’ARNm ASH1 naissant. Puf6p est également recruté sur ASH1, mais d’une manière dépendante de la présence de She2p. De façon intéressante, nous avons détecté une interaction entre She2p et la plus grande sous-unité de l’ARN polymérase II (Rpb1p). Cette interaction est détectée avec la forme active en élongation de l’ARN polymérase II. Nous avons également démontré que She2p interagit avec le complexe d’élongation de la transcription Spt4p/Spt5p. Une délétion de SPT4 ou une mutation dans SPT5 (Ts spt5) à température restrictive empêche l’interaction entre She2p et Rpb1p, et diminue le recrutement de She2p au gène ASH1, entrainant un défaut de localisation de l’ARNm et un défaut de localisation asymétrique de la protéine Ash1p. De manière globale, nos résultats montrent que les facteurs impliqués dans la localisation cytoplasmique des ARNm et dans leur contrôle traductionnel sont recrutés de façon co-transcriptionnelle sur les ARNm naissants via leur interaction avec la machinerie de transcription, suggèrant un rôle important de la machinerie transcriptionelle dans la localisation des ARNm.

Développement d’outils pour l’analyse de données de ChIP-seq et l’identification des facteurs de transcription

La méthode ChIP-seq est une technologie combinant la technique de chromatine immunoprecipitation avec le séquençage haut-débit et permettant l’analyse in vivo des facteurs de transcription à grande échelle. Le traitement des grandes quantités de données ainsi générées nécessite des moyens informatiques performants et de nombreux outils ont vu le jour récemment. Reste cependant que cette multiplication des logiciels réalisant chacun une étape de l’analyse engendre des problèmes de compatibilité et complique les analyses. Il existe ainsi un besoin important pour une suite de logiciels performante et flexible permettant l’identification des motifs. Nous proposons ici un ensemble complet d’analyse de données ChIP-seq disponible librement dans R et composé de trois modules PICS, rGADEM et MotIV. A travers l’analyse de quatre jeux de données des facteurs de transcription CTCF, STAT1, FOXA1 et ER nous avons démontré l’efficacité de notre ensemble d’analyse et mis en avant les fonctionnalités novatrices de celui-ci, notamment concernant le traitement des résultats par MotIV conduisant à la découverte de motifs non détectés par les autres algorithmes.

Une nouvelle approche computationnelle pour la découverte des sites de fixation de facteurs de transcription à l’ADN, adaptée aux données de ChIP-chip et de ChIP-séquençage

Les facteurs de transcription sont des protéines spécialisées qui jouent un rôle important dans différents processus biologiques tel que la différenciation, le cycle cellulaire et la tumorigenèse. Ils régulent la transcription des gènes en se fixant sur des séquences d’ADN spécifiques (éléments cis-régulateurs). L’identification de ces éléments est une étape cruciale dans la compréhension des réseaux de régulation des gènes. Avec l’avènement des technologies de séquençage à haut débit, l’identification de tout les éléments fonctionnels dans les génomes, incluant gènes et éléments cis-régulateurs a connu une avancée considérable. Alors qu’on est arrivé à estimer le nombre de gènes chez différentes espèces, l’information sur les éléments qui contrôlent et orchestrent la régulation de ces gènes est encore mal définie. Grace aux techniques de ChIP-chip et de ChIP-séquençage il est possible d’identifier toutes les régions du génome qui sont liées par un facteur de transcription d’intérêt. Plusieurs approches computationnelles ont été développées pour prédire les sites fixés par les facteurs de transcription. Ces approches sont classées en deux catégories principales: les algorithmes énumératifs et probabilistes. Toutefois, plusieurs études ont montré que ces approches génèrent des taux élevés de faux négatifs et de faux positifs ce qui rend difficile l’interprétation des résultats et par conséquent leur validation expérimentale. Dans cette thèse, nous avons ciblé deux objectifs. Le premier objectif a été de développer une nouvelle approche pour la découverte des sites de fixation des facteurs de transcription à l’ADN (SAMD-ChIP) adaptée aux données de ChIP-chip et de ChIP-séquençage. Notre approche implémente un algorithme hybride qui combine les deux stratégies énumérative et probabiliste, afin d’exploiter les performances de chacune d’entre elles. Notre approche a montré ses performances, comparée aux outils de découvertes de motifs existants sur des jeux de données simulées et des jeux de données de ChIP-chip et de ChIP-séquençage. SAMD-ChIP présente aussi l’avantage d’exploiter les propriétés de distributions des sites liés par les facteurs de transcription autour du centre des régions liées afin de limiter la prédiction aux motifs qui sont enrichis dans une fenêtre de longueur fixe autour du centre de ces régions. Les facteurs de transcription agissent rarement seuls. Ils forment souvent des complexes pour interagir avec l’ADN pour réguler leurs gènes cibles. Ces interactions impliquent des facteurs de transcription dont les sites de fixation à l’ADN sont localisés proches les uns des autres ou bien médier par des boucles de chromatine. Notre deuxième objectif a été d’exploiter la proximité spatiale des sites liés par les facteurs de transcription dans les régions de ChIP-chip et de ChIP-séquençage pour développer une approche pour la prédiction des motifs composites (motifs composés par deux sites et séparés par un espacement de taille fixe). Nous avons testé ce module pour prédire la co-localisation entre les deux demi-sites ERE qui forment le site ERE, lié par le récepteur des œstrogènes ERα. Ce module a été incorporé à notre outil de découverte de motifs SAMD-ChIP.

Affectation de composantes basée sur des contraintes énergétiques dans une architecture multiprocesseurs en trois dimensions

La lithographie et la loi de Moore ont permis des avancées extraordinaires dans la fabrication des circuits intégrés. De nos jours, plusieurs systèmes très complexes peuvent être embarqués sur la même puce électronique. Les contraintes de développement de ces systèmes sont tellement grandes qu’une bonne planification dès le début de leur cycle de développement est incontournable. Ainsi, la planification de la gestion énergétique au début du cycle de développement est devenue une phase importante dans la conception de ces systèmes. Pendant plusieurs années, l’idée était de réduire la consommation énergétique en ajoutant un mécanisme physique une fois le circuit créé, comme par exemple un dissipateur de chaleur. La stratégie actuelle est d’intégrer les contraintes énergétiques dès les premières phases de la conception des circuits. Il est donc essentiel de bien connaître la dissipation d’énergie avant l’intégration des composantes dans une architecture d’un système multiprocesseurs de façon à ce que chaque composante puisse fonctionner efficacement dans les limites de ses contraintes thermiques. Lorsqu’une composante fonctionne, elle consomme de l’énergie électrique qui est transformée en dégagement de chaleur. Le but de ce mémoire est de trouver une affectation efficace des composantes dans une architecture de multiprocesseurs en trois dimensions en tenant compte des limites des facteurs thermiques de ce système.