Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Mesurer la stigmatisation perçue chez les personnes souffrant de troubles psychiques : traduction française, validation et adaptation de la Stigma Scale [French translation, validation and adaptation of the Stigma Scale].

AIM: People suffering from mental illness are exposed to stigma. However, only few tools are available to assess stigmatization as perceived from the patient's perspective. The aim of this study is to adapt and validate a French version of the Stigma Scale (King et al., 2007 [8]). This self-report questionnaire has a three-factor structure: discrimination, disclosure and positive aspects of mental illness. Discrimination subscale refers to perceived negative reactions of others. Disclosure subscale refers mainly to managing disclosure to avoid discrimination and finally positive aspects subscale taps into how patients are becoming more accepting, more understanding toward their illness. METHOD: In the first step, internal consistency, convergent validity and test-retest reliability of the French adaptation of the 28-item scale were assessed in a sample of 183 patients. Results of confirmatory factor analyses (CFA) did not confirm the hypothesized structure. In the light of the failed attempts to validate the original version, an alternative 9-item short-form version of the Stigma Scale, maintaining the integrity of the original model, was developed based on results of exploratory factor analyses in the first sample and cross-validated in a new sample of 234 patients. RESULTS: Results of CFA did not confirm that the data fitted well to the three-factor model of the 28-item Stigma Scale (χ(2)/df=2.02, GFI=0.77, AGFI=0.73, RMSEA=0.07, CFI=0.77 and NNFI=0.75). Cronbach's α was excellent for discrimination (0.84) and disclosure (0.83) subscales but poor for potential positive aspects (0.46). External validity was satisfactory. Overall Stigma Scale total score was negatively correlated with the score on Rosenberg's Self-Esteem Scale (r=-0.49), and each subscale was significantly correlated with a visual analogue scale that referred to the specific aspect of stigma (0.43≤|r|≤0.60). Intraclass correlation coefficients between 0.68 and 0.89 indicated good test-retest reliability. The results of the CFA demonstrated that the items chosen for the short version of the Stigma Scale have the expected fit properties (χ(2)/df=1.02, GFI=0.98, AGFI=0.98, RMSEA=0.01, CFI=1.0 and NNFI=1.0). Considering the small number (three) of items in each subscale of the short version of the Stigma Scale, α coefficients for discrimination (0.57), disclosure (0.80) and potential positive aspects subscales (0.62) are considered as good. CONCLUSION: Our results suggest that the 9-item French short version of the Stigma Scale is a useful, reliable and valid self-report questionnaire to assess perceived stigmatization in people suffering from mental illness. The time of completion is really short and questions are well understood and accepted by the patients.

India in Translation through Hindi Literature. A Plurality of Voices

What role have translations from Hindi literary works played in shaping and transforming our knowledge about India? In this book, renowned scholars, translators and Hindi writers from India, Europe, and the United States offer their approaches to this question. Their articles deal with the political, cultural, and linguistic criteria germane to the selection and translation of Hindi works, the nature of the enduring links between India and Europe, and the reception of translated texts, particularly through the perspective of book history. More personal essays, both on the writing process itself or on the practice of translation, complete the volume and highlight the plurality of voices that are inherent to any translation. As the outcome of an international symposium held at the University of Lausanne, Switzerland, in 2008, India in Translation through Hindi Literature engages in the building of critical histories of the encounter between India and the «West», the use and impact of translations in this context, and Hindi literature and culture in connection to English (post)colonial power, literature and culture.

Low-cost MR-compatible moving heart phantom.

For the development and evaluation of cardiac magnetic resonance (MR) imaging sequences and methodologies, the availability of a periodically moving phantom to model respiratory and cardiac motion would be of substantial benefit. Given the specific physical boundary conditions in an MR environment, the choice of materials and power source of such phantoms is heavily restricted. Sophisticated commercial solutions are available; however, they are often relatively costly and user-specific modifications may not easily be implemented. We therefore sought to construct a low-cost MR-compatible motion phantom that could be easily reproduced and had design flexibility. A commercially available K'NEX construction set (Hyper Space Training Tower, K'NEX Industries, Inc., Hatfield, PA) was used to construct a periodically moving phantom head. The phantom head performs a translation with a superimposed rotation, driven by a motor over a 2-m rigid rod. To synchronize the MR data acquisition with phantom motion (without introducing radiofrequency-related image artifacts), a fiberoptic control unit generates periodic trigger pulses synchronized to the phantom motion. Total material costs of the phantom are US$ < 200.00, and a total of 80 man-hours were required to design and construct the original phantom. With schematics of the present solution, the phantom reproduction may be achieved in approximately 15 man-hours. The presented MR-compatible periodically moving phantom can easily be reproduced, and user-specific modifications may be implemented. Such an approach allows a detailed investigation of motion-related phenomena in MR images.

Ornament of Reason : an annoted english translation of rMa bya Byang chub brtson 'grus's : commentary to the Mūlamadhyamakakārikā

The present work contains an English translation of the Ornament of Reason (Tib. dBu ma'i rtsa ba shes rab kyi 'grel ba thad pa'i rgyan), an influential 12th century Tibetan commentary to a classic of Indian Buddhism, Nàgàrjuna's Root of the Middle Way (Skt. Mûlamadhyamaka-kârikâ Tib. dBu ma'i rtsa ba shes rab ces bya ba). The author, rMa bya byang chub brtson 'grus's, was among the first Tibetans to explicitly align themselves with the philosophical project of the 7th century Indian Madhyamaka teacher, CandrakTrti. With his stern revision and transformative re-employment of the epistemological framework that CandrakTrti otherwise scorns, rMa bya was however to develop an original and highly influential Madhyamaka interpretation.¦rMa bya's "topical outline" (Tib. sa bead) of his commentary is contained in the first appendix to the translation. The second appendix contains a digitalized version of the Rumtek edition of rMa bya's commentary, cross-referenced to the page numbers of the translation. The annotations contain concise explanatory remarks and suggestions for further reading, as well as identifications of the passages that rMa bya cites.

Translation and assembly of HLA-DR antigens in Xenopus oocytes injected with mRNA from a human B-cell line.

HLA-DR antigens are polymorphic cell surface glycoproteins, expressed primarily in B lymphocytes and macrophages, which are thought to play an important role in the immune response. Two polypeptide chains, alpha and beta, are associated at the cell surface, and a third chain associates with alpha and beta intracellularly. RNA isolated from the human B-cell line Raji was injected in Xenopus laevis oocytes. Immunoprecipitates of translation products with several monoclonal antibodies revealed the presence of HLA-DR antigens similar to those synthesized in Raji cells. One monoclonal antibody was able to bind the beta chain after dissociation of the three polypeptide chains with detergent. The presence of all three chains was confirmed by two-dimensional gel electrophoresis. The glycosylation pattern of the three chains was identical to that observed in vivo, as evidenced in studies using tunicamycin, an inhibitor of N-linked glycosylation. The presence of alpha chains assembled with beta chains in equimolar ratio was further demonstrated by amino-terminal sequencing. An RNA fraction enriched for the three mRNAs, encoding alpha, beta, and intracellular chains, was isolated. This translation-assembly system and the availability of monoclonal antibodies make it possible to assay for mRNA encoding specific molecules among the multiple human Ia-like antigens.

Translation of polarity cues into asymmetric spindle positioning in "Caenorhabditis elegans"

Abstract: Asymmetric cell division is important to generate tissue diversity. The Caenorhabditis elegans embryo is well suited to study the mechanisms of asymmetric cell division. In wild type one-cell stage embryos, the spindle sets up along the anterior-posterior axis (AP). During anaphase, the spindle elongates. While the anterior spindle pole is relatively immobile, the posterior spindle pole moves towards the posterior cortex during anaphase leading to an asymmetric spindle position. As a result, the first cleavage gives rise to a large anterior blastomere and a smaller posterior one, which differs also in cell fate determinants. This posterior spindle displacement occurs in response to polarity cues set up along the AP axis by the PAR proteins and is due to imbalanced pulling forces acting on the two spindle poles, with net forces acting on the posterior spindle pole being more extensive than those at the anterior one. The project of my thesis was to characterize the involvement of two new components, gpr-1 and gpr-2, in spindle positioning. These genes encode essentially identical proteins containing a GoLoco motif characteristic of proteins interacting with α subunits of heterotrimeric G protein (Gα). In gpr-1/2(RNAi) embryos and in embryos lacking simultaneously two α subunits, goa-1 and gpa-16, (Ga(RNAi) embryos), there is a minimal posterior displacement of the spindle during anaphase, and the first division is equal. I found that the pulling forces acting on the two spindle poles is weak and equal in gpr-1/2(RNAi) and Gα (RNAi) embryos. I found that GPR-1/2 acts downstream of polarity cues for generation of pulling forces. Furthermore, I showed that GPR-1/2 distribution was enriched at the posterior cortex during metaphase whereas GOA-1 and GPA-16 were uniformly distributed at the cell cortex throughout the cell cycle. Gα subunits oscillate between GDP- and GTP-bound forms. Gα signaling is turned on by GDP/GTP exchange catalyzed by guanine nucleotide exchange factors (GEFs) and turned off by hydrolysis of GTP catalyzed by GTPase activating proteins (GAPs). A third class of proteins, the guanine dissociation inhibitors (GDIs), binds the GDP-bound form of Gα subunits and inhibits nucleotide exchange. I found that GPR-1/2 acts as a GDI for GOA-1. Taken together, my findings suggest a model in which differential activation of Gα subunits along the AP axis may translate into generation of differential pulling forces on the anterior and posterior spindle poles, and, thus, asymmetric cell division. Résumé L'embryon du nématode Caenorhabditis elegans est un modèle approprié pour étudier les mécanismes de la division asymétrique. Chez l'embryon précoce, le fuseau mitotique se forme le long de l'axe antéro-postérieur (A/P) et au centre de l'embryon, le pôle antérieur restant relativement immobile alors que le pôle postérieur du fuseau se déplace vers le cortex postérieur au cours de l'anaphase conduisant à une position excentrée du fuseau. 11 en résulte une première division qui génère un blastomère antérieur et postérieur de grande et petite taille respectivement et qui diffèrent en facteurs développementaux. Ce déplacement postérieur se produit en réponse de la polarité établie par la distribution polarisée des protéines PAR et est le résultat de la génération de forces inégales tirant sur les deux pôles du fuseau, les forces agissant sur le pôle postérieur du fuseau étant plus grandes. Le projet de ma thèse était d'identifier la fonction de deux nouveaux constituants, gpr-1 et gpr-2 dans le positionnement asymétrique du fuseau. Ces gènes codent essentiellement pour la même protéine qui contient un motif GoLoco, caractéristique des protéines interagissant avec la sous-unité alpha des protéines G hétérotrimériques. Chez l'embryon gpr-1/2(RNAi) et chez les embryons dépourvus d'activité de deux sous-unités alpha, goa-1 et gpa-16, (Gα(RNAi)), j'ai montré qu'il y avait un déplacement minimal du fuseau vers le pôle postérieur au cours de l'anaphase et la première division est symétrique en raison de forces faibles et égales agissant sur les deux pôles du fuseau. J'ai également montré que gpr-1/2 était requis en aval des signaux établissant la polarité pour générer les forces responsables du positionnement asymétrique du fuseau. De plus, j'ai montré que GPR-1/2 était enrichi au pôle postérieur lors de la métaphase alors que GOA-1 et GPA-16 étaient localisés de façon uniforme au cortex de l'embryon précoce. Gas oscillent entre une forme liée au GDP et une forme liée au GTP. La signalisation des Gas est activée par l'échange GDP/GTP qui est catalysé par des protéines GEFs. La signalisation des Gas est désactivée par l'hydrolyse du GTP qui est catalysée par des protéines GAPs. Une troisième classe de protéines, GDIs lie la forme GDP et inhibe l'échange de nucléotides. J'ai montré que GPR-1/2 agissait comme un GDI pour GOA-1. Mes résultats suggèrent un modèle dans lequel une activation différentielle des Gα le long de l'axe A/P pourrait générer des forces différentielles sur le pôle antérieur et postérieur du fuseau.

Development and application of computational methodologies in qualitative modeling

La présente thèse s'intitule "Développent et Application des Méthodologies Computationnelles pour la Modélisation Qualitative". Elle comprend tous les différents projets que j'ai entrepris en tant que doctorante. Plutôt qu'une mise en oeuvre systématique d'un cadre défini a priori, cette thèse devrait être considérée comme une exploration des méthodes qui peuvent nous aider à déduire le plan de processus regulatoires et de signalisation. Cette exploration a été mue par des questions biologiques concrètes, plutôt que par des investigations théoriques. Bien que tous les projets aient inclus des systèmes divergents (réseaux régulateurs de gènes du cycle cellulaire, réseaux de signalisation de cellules pulmonaires) ainsi que des organismes (levure à fission, levure bourgeonnante, rat, humain), nos objectifs étaient complémentaires et cohérents. Le projet principal de la thèse est la modélisation du réseau de l'initiation de septation (SIN) du S.pombe. La cytokinèse dans la levure à fission est contrôlée par le SIN, un réseau signalant de protéines kinases qui utilise le corps à pôle-fuseau comme échafaudage. Afin de décrire le comportement qualitatif du système et prédire des comportements mutants inconnus, nous avons décidé d'adopter l'approche de la modélisation booléenne. Dans cette thèse, nous présentons la construction d'un modèle booléen étendu du SIN, comprenant la plupart des composantes et des régulateurs du SIN en tant que noeuds individuels et testable expérimentalement. Ce modèle utilise des niveaux d'activité du CDK comme noeuds de contrôle pour la simulation d'évènements du SIN à différents stades du cycle cellulaire. Ce modèle a été optimisé en utilisant des expériences d'un seul "knock-out" avec des effets phénotypiques connus comme set d'entraînement. Il a permis de prédire correctement un set d'évaluation de "knock-out" doubles. De plus, le modèle a fait des prédictions in silico qui ont été validées in vivo, permettant d'obtenir de nouvelles idées de la régulation et l'organisation hiérarchique du SIN. Un autre projet concernant le cycle cellulaire qui fait partie de cette thèse a été la construction d'un modèle qualitatif et minimal de la réciprocité des cyclines dans la S.cerevisiae. Les protéines Clb dans la levure bourgeonnante présentent une activation et une dégradation caractéristique et séquentielle durant le cycle cellulaire, qu'on appelle communément les vagues des Clbs. Cet évènement est coordonné avec la courbe d'activation inverse du Sic1, qui a un rôle inhibitoire dans le système. Pour l'identification des modèles qualitatifs minimaux qui peuvent expliquer ce phénomène, nous avons sélectionné des expériences bien définies et construit tous les modèles minimaux possibles qui, une fois simulés, reproduisent les résultats attendus. Les modèles ont été filtrés en utilisant des simulations ODE qualitatives et standardisées; seules celles qui reproduisaient le phénotype des vagues ont été gardées. L'ensemble des modèles minimaux peut être utilisé pour suggérer des relations regulatoires entre les molécules participant qui peuvent ensuite être testées expérimentalement. Enfin, durant mon doctorat, j'ai participé au SBV Improver Challenge. Le but était de déduire des réseaux spécifiques à des espèces (humain et rat) en utilisant des données de phosphoprotéines, d'expressions des gènes et des cytokines, ainsi qu'un réseau de référence, qui était mis à disposition comme donnée préalable. Notre solution pour ce concours a pris la troisième place. L'approche utilisée est expliquée en détail dans le dernier chapitre de la thèse. -- The present dissertation is entitled "Development and Application of Computational Methodologies in Qualitative Modeling". It encompasses the diverse projects that were undertaken during my time as a PhD student. Instead of a systematic implementation of a framework defined a priori, this thesis should be considered as an exploration of the methods that can help us infer the blueprint of regulatory and signaling processes. This exploration was driven by concrete biological questions, rather than theoretical investigation. Even though the projects involved divergent systems (gene regulatory networks of cell cycle, signaling networks in lung cells), as well as organisms (fission yeast, budding yeast, rat, human), our goals were complementary and coherent. The main project of the thesis is the modeling of the Septation Initiation Network (SIN) in S.pombe. Cytokinesis in fission yeast is controlled by the SIN, a protein kinase signaling network that uses the spindle pole body as scaffold. In order to describe the qualitative behavior of the system and predict unknown mutant behaviors we decided to adopt a Boolean modeling approach. In this thesis, we report the construction of an extended, Boolean model of the SIN, comprising most SIN components and regulators as individual, experimentally testable nodes. The model uses CDK activity levels as control nodes for the simulation of SIN related events in different stages of the cell cycle. The model was optimized using single knock-out experiments of known phenotypic effect as a training set, and was able to correctly predict a double knock-out test set. Moreover, the model has made in silico predictions that have been validated in vivo, providing new insights into the regulation and hierarchical organization of the SIN. Another cell cycle related project that is part of this thesis was to create a qualitative, minimal model of cyclin interplay in S.cerevisiae. CLB proteins in budding yeast present a characteristic, sequential activation and decay during the cell cycle, commonly referred to as Clb waves. This event is coordinated with the inverse activation curve of Sic1, which has an inhibitory role in the system. To generate minimal qualitative models that can explain this phenomenon, we selected well-defined experiments and constructed all possible minimal models that, when simulated, reproduce the expected results. The models were filtered using standardized qualitative ODE simulations; only the ones reproducing the wave-like phenotype were kept. The set of minimal models can be used to suggest regulatory relations among the participating molecules, which will subsequently be tested experimentally. Finally, during my PhD I participated in the SBV Improver Challenge. The goal was to infer species-specific (human and rat) networks, using phosphoprotein, gene expression and cytokine data and a reference network provided as prior knowledge. Our solution to the challenge was selected as in the final chapter of the thesis.

In vitro assays and clinical trials in red blood cell aging: Lost in translation.

The age of erythrocyte concentrates (EC) in transfusion medicine and the adverse outcomes when transfusing long-term-stored EC are highly controversial issues. Whereas the definition of a short-term-stored EC or a long-term-stored EC is unclear in clinical trials, data based on in vitro storage assays can help defining a limit in addition of the expiration date. The present review merges together these data in order to highlight an EC age cut-off and points out potential misleading consideration. The analysis of in vitro data highlights the presence of reversible and irreversible storage lesions and demonstrates that red blood cells (RBC) exhibit two limits during storage: one around 2 weeks and another one around 4 weeks of storage. Of particular importance, the first lesions to appear, i.e. the reversible ones, are per se reversible once transfused, whereas the irreversible lesions are not. In clinical trials, the EC age cut-off for short-term storage is in general fewer than 14 days (11 ± 4 days) and more disperse for long-term-stored EC (17 ± 13 days), regardless the clinical outcomes. Taking together, EC age cut-off in clinical trials does not totally fall into line of in vitro aging data, whereas it is the key criteria in clinical studies. Long-term-stored EC considered in clinical trials are not probably old enough to answer the question: "Does transfusion of long-term-stored EC (older than 4 weeks) result in worse clinical outcomes?" Depending on ethical concerns and clinical practices, older EC than currently assayed in clinical trials should have to be considered. These two worlds trying to understand the aging of erythrocytes and the impact on patients do not seem to speak the same language.

Lost in translation: dietary management of cardiovascular risk factors is seldom implemented.

OBJECTIVE: To assess dietary management of cardiovascular risk factors (CVRFs) in the general population. METHOD: Cross-sectional study conducted between 2009 and 2012 on 4811 participants (2567 women, 58±11years) living in Lausanne, Switzerland. RESULTS: Sixteen percent of participants diagnosed with overweight/obesity reported a slimming diet. Slimming diet was associated with diagnosis of hypertension: Odds ratio and (95% confidence interval): 0.61 (0.40-0.93); older age [0.84 (0.58-1.21), 0.79 (0.53-1.18) and 0.47 (0.27-0.81) for [50-60[, [60-70[ and [70+ years, respectively]; female gender [1.84 (1.36-2.48)] and diagnosis of diabetes [2.16 (1.13-4.12)]. Only 8% of participants diagnosed with hypertension reported a low-salt diet. Low-salt diet was associated with antihypertensive drug treatment [2.17 (1.28-3.68)] and diagnosis of diabetes [2.72 (1.26-5.86)]. One-third of participants diagnosed with dyslipidemia reported a low-fat diet. Low-fat diet was associated with female gender [1.47 (1.17-1.86)]; older age [1.29 (0.89-1.87), 1.71 (1.18-2.48) and 2.01 (1.33-3.03) for [50-60[, [60-70[ and [70+ years, respectively]; hypolipidemic drug treatment [OR=1.68 (1.29-2.18)]; current smoking [0.70 (0.51-0.96)] and obesity [0.67 (0.45-1.00)]. Approximately half of participants diagnosed with diabetes reported an antidiabetic diet. Antidiabetic diet was associated with current smoking [0.44 (0.22-0.88)] and antidiabetic drug treatment [OR=3.26 (1.81-5.86)]. CONCLUSION: Dietary management of CVRFs is seldom implemented in Switzerland.