A multi-level multi-scale approach to study essential genes in Mycobacterium tuberculosis

Background: The set of indispensable genes that are required by an organism to grow and sustain life are termed as essential genes. There is a strong interest in identification of the set of essential genes, particularly in pathogens, not only for a better understanding of the pathogen biology, but also for identifying drug targets and the minimal gene set for the organism. Essentiality is inherently a systems property and requires consideration of the system as a whole for their identification. The available experimental approaches capture some aspects but each method comes with its own limitations. Moreover, they do not explain the basis for essentiality in most cases. A powerful prediction method to recognize this gene pool including rationalization of the known essential genes in a given organism would be very useful. Here we describe a multi-level multi-scale approach to identify the essential gene pool in a deadly pathogen, Mycobacterium tuberculosis. Results: The multi-level workflow analyses the bacterial cell by studying (a) genome-wide gene expression profiles to identify the set of genes which show consistent and significant levels of expression in multiple samples of the same condition, (b) indispensability for growth by using gene expression integrated flux balance analysis of a genome-scale metabolic model, (c) importance for maintaining the integrity and flow in a protein-protein interaction network and (d) evolutionary conservation in a set of genomes of the same ecological niche. In the gene pool identified, the functional basis for essentiality has been addressed by studying residue level conservation and the sub-structure at the ligand binding pockets, from which essential amino acid residues in that pocket have also been identified. 283 genes were identified as essential genes with high-confidence. An agreement of about 73.5% is observed with that obtained from the experimental transposon mutagenesis technique. A large proportion of the identified genes belong to the class of intermediary metabolism and respiration. Conclusions: The multi-scale, multi-level approach described can be generally applied to other pathogens as well. The essential gene pool identified form a basis for designing experiments to probe their finer functional roles and also serve as a ready shortlist for identifying drug targets.

Wide Screening of Phage-Displayed Libraries Identifies Immune Targets in Planta

Microbe-Associated Molecular Patterns and virulence effectors are recognized by plants as a first step to mount a defence response against potential pathogens. This recognition involves a large family of extracellular membrane receptors and other immune proteins located in different sub-cellular compartments. We have used phage-display technology to express and select for Arabidopsis proteins able to bind bacterial pathogens. To rapidly identify microbe-bound phage, we developed a monitoring method based on microarrays. This combined strategy allowed for a genome-wide screening of plant proteins involved in pathogen perception. Two phage libraries for high-throughput selection were constructed from cDNA of plants infected with Pseudomonas aeruginosa PA14, or from combined samples of the virulent isolate DC3000 of Pseudomonas syringae pv. tomato and its avirulent variant avrRpt2. These three pathosystems represent different degrees in the specificity of plant-microbe interactions. Libraries cover up to 26107 different plant transcripts that can be displayed as functional proteins on the surface of T7 bacteriophage. A number of these were selected in a bio-panning assay for binding to Pseudomonas cells. Among the selected clones we isolated the ethylene response factor ATERF-1, which was able to bind the three bacterial strains in competition assays. ATERF-1 was rapidly exported from the nucleus upon infiltration of either alive or heat-killed Pseudomonas. Moreover, aterf-1 mutants exhibited enhanced susceptibility to infection. These findings suggest that ATERF-1 contains a microbe-recognition domain with a role in plant defence. To identify other putative pathogen-binding proteins on a genome-wide scale, the copy number of selected-vs.-total clones was compared by hybridizing phage cDNAs with Arabidopsis microarrays. Microarray analysis revealed a set of 472 candidates with significant fold change. Within this set defence-related genes, including well-known targets of bacterial effectors, are over-represented. Other genes non-previously related to defence can be associated through this study with general or strain-specific recognition of Pseudomonas.

Functional genomic studies of the structure and regulation of eukaryotic transcriptomes

The main focus of this thesis is the use of high-throughput sequencing technologies in functional genomics (in particular in the form of ChIP-seq, chromatin immunoprecipitation coupled with sequencing, and RNA-seq) and the study of the structure and regulation of transcriptomes. Some parts of it are of a more methodological nature while others describe the application of these functional genomic tools to address various biological problems. A significant part of the research presented here was conducted as part of the ENCODE (ENCyclopedia Of DNA Elements) Project.

The first part of the thesis focuses on the structure and diversity of the human transcriptome. Chapter 1 contains an analysis of the diversity of the human polyadenylated transcriptome based on RNA-seq data generated for the ENCODE Project. Chapter 2 presents a simulation-based examination of the performance of some of the most popular computational tools used to assemble and quantify transcriptomes. Chapter 3 includes a study of variation in gene expression, alternative splicing and allelic expression bias on the single-cell level and on a genome-wide scale in human lymphoblastoid cells; it also brings forward a number of critical to the practice of single-cell RNA-seq measurements methodological considerations.

The second part presents several studies applying functional genomic tools to the study of the regulatory biology of organellar genomes, primarily in mammals but also in plants. Chapter 5 contains an analysis of the occupancy of the human mitochondrial genome by TFAM, an important structural and regulatory protein in mitochondria, using ChIP-seq. In Chapter 6, the mitochondrial DNA occupancy of the TFB2M transcriptional regulator, the MTERF termination factor, and the mitochondrial RNA and DNA polymerases is characterized. Chapter 7 consists of an investigation into the curious phenomenon of the physical association of nuclear transcription factors with mitochondrial DNA, based on the diverse collections of transcription factor ChIP-seq datasets generated by the ENCODE, mouseENCODE and modENCODE consortia. In Chapter 8 this line of research is further extended to existing publicly available ChIP-seq datasets in plants and their mitochondrial and plastid genomes.

The third part is dedicated to the analytical and experimental practice of ChIP-seq. As part of the ENCODE Project, a set of metrics for assessing the quality of ChIP-seq experiments was developed, and the results of this activity are presented in Chapter 9. These metrics were later used to carry out a global analysis of ChIP-seq quality in the published literature (Chapter 10). In Chapter 11, the development and initial application of an automated robotic ChIP-seq (in which these metrics also played a major role) is presented.

The fourth part presents the results of some additional projects the author has been involved in, including the study of the role of the Piwi protein in the transcriptional regulation of transposon expression in Drosophila (Chapter 12), and the use of single-cell RNA-seq to characterize the heterogeneity of gene expression during cellular reprogramming (Chapter 13).

The last part of the thesis provides a review of the results of the ENCODE Project and the interpretation of the complexity of the biochemical activity exhibited by mammalian genomes that they have revealed (Chapters 15 and 16), an overview of the expected in the near future technical developments and their impact on the field of functional genomics (Chapter 14), and a discussion of some so far insufficiently explored research areas, the future study of which will, in the opinion of the author, provide deep insights into many fundamental but not yet completely answered questions about the transcriptional biology of eukaryotes and its regulation.

Gene duplication in the genome of parasitic Giardia lamblia

Background: Giardia are a group of widespread intestinal protozoan parasites in a number of vertebrates. Much evidence from G. lamblia indicated they might be the most primitive extant eukaryotes. When and how such a group of the earliest branching unicellular eukaryotes developed the ability to successfully parasitize the latest branching higher eukaryotes (vertebrates) is an intriguing question. Gene duplication has long been thought to be the most common mechanism in the production of primary resources for the origin of evolutionary novelties. In order to parse the evolutionary trajectory of Giardia parasitic lifestyle, here we carried out a genome-wide analysis about gene duplication patterns in G. lamblia. Results: Although genomic comparison showed that in G. lamblia the contents of many fundamental biologic pathways are simplified and the whole genome is very compact, in our study 40% of its genes were identified as duplicated genes. Evolutionary distance analyses of these duplicated genes indicated two rounds of large scale duplication events had occurred in G. lamblia genome. Functional annotation of them further showed that the majority of recent duplicated genes are VSPs (Variant-specific Surface Proteins), which are essential for the successful parasitic life of Giardia in hosts. Based on evolutionary comparison with their hosts, it was found that the rapid expansion of VSPs in G. lamblia is consistent with the evolutionary radiation of placental mammals. Conclusions: Based on the genome-wide analysis of duplicated genes in G. lamblia, we found that gene duplication was essential for the origin and evolution of Giardia parasitic lifestyle. The recent expansion of VSPs uniquely occurring in G. lamblia is consistent with the increment of its hosts. Therefore we proposed a hypothesis that the increment of Giradia hosts might be the driving force for the rapid expansion of VSPs.

Sequence context analysis in the mouse genome: Single nucleotide polymorphisms and CpG island sequences

A genome-wide view of sequence mutability in mice is still limited, although biologists usually assume the same scenario for mice as for humans. In this study, we examined the sequence context in the local environment of 482,528 mouse single nucleotide po

Origin of avian genome size and structure in non-avian dinosaurs

Chris L. Organ, Andrew M. Shedlock, Andrew Meade, Mark Pagel and Scott V. Edwards (2007). Origin of avian genome size and structure in non-avian dinosaurs. Nature, 46(7132), 180-184. RAE2008

Mapping the gene for autosomal dominant restless legs syndrome in an Irish family

Restless Legs Syndrome (RLS) is a common neurological disorder affecting nearly 15% of the general population. Ironically, RLS can be described as the most common condition one has never heard of. It is usually characterised by uncomfortable, unpleasant sensations in the lower limbs inducing an uncontrollable desire to move the legs. RLS exhibits a circadian pattern with symptoms present predominantly in the evening or at night, thus leading to sleep disruption and daytime somnolence. RLS is generally classified into primary (idiopathic) and secondary (symptomatic) forms. Primary RLS includes sporadic and familial cases of which the age of onset is usually less than 45 years and progresses slowly with a female to male ratio of 2:1. Secondary forms often occur as a complication of another health condition, such as iron deficiency or thyroid dysfunction. The age of onset is usually over 45 years, with an equal male to female ratio and more rapid progression. Ekbom described the familial component of the disorder in 1945 and since then many studies have been published on the familial forms of the disorder. Molecular genetic studies have so far identified ten loci (5q, 12q, 14p, 9p, 20p, 16p, 19p, 4q, 17p). No specific gene within these loci has been identified thus far. Association mapping has highlighted a further five areas of interest. RLS6 has been found to be associated with SNPs in the BTBD9 gene. Four other variants were found within intronic and intergenic regions of MEIS1, MAP2K5/LBXCOR1, PTPRD and NOS1. The pathophysiology of RLS is complex and remains to be fully elucidated. Conditions associated with secondary RLS, such as pregnancy or end-stage renal disease, are characterised by iron deficiency, which suggests that disturbed iron homeostasis plays a role. Dopaminergic dysfunction in subcortical systems also appears to play a central role. An ongoing study within the Department of Pathology (University College Cork) is investigating the genetic characteristics of RLS in Irish families. A three generation RLS pedigree RLS3002 consisting of 11 affected and 7 unaffected living family members was recruited. The family had been examined for four of the known loci (5q, 12q, 14p and 9p) (Abdulrahim 2008). The aim of this study was to continue examining this Irish RLS pedigree for possible linkage to the previously described loci and associated regions. Using informative microsatellite markers linkage was excluded to the loci on 5q, 12q, 14p, 9p, 20p, 16p, 19p, 4q, 17p and also within the regions reported to be associated with RLS. This suggested the presence of a new unidentified locus. A genome-wide scan was performed using two microsatellite marker screening sets (Research Genetics Inc. Mapping set and the Applied Biosystems Linkage mapping set version 2.5). Linkage analysis was conducted under an autosomal dominant model with a penetrance of 95% and an allele frequency of 0.01. A maximum LOD score of 3.59 at θ=0.00 for marker D19S878 indicated significant linkage on chromosome 19p. Haplotype analysis defined a genetic region of 6.57 cM on chromosome 19p13.3, corresponding to 2.5 Mb. There are approximately 100 genes annotated within the critical region. Sequencing of two candidate genes, KLF16 and GAMT, selected on the assumed pathophysiology of RLS, did not identify any sequence variant. This study provides evidence of a novel RLS locus in an Irish pedigree, thus supporting the picture of RLS as a genetically heterogeneous trait.

Association Analysis of the PIP4K2A Gene on Chromosome 10p12 and Schizophrenia in the Irish Study of High Density Schizophrenia Families (ISHDSF) and the Irish Case-Control Study of Schizophrenia (ICCSS)

Molecular studies support pharmacological evidence that phosphoinositide signaling is perturbed in schizophrenia and bipolar disorder. The phosphatidylinositol-4-phosphate-5-kinase type-II alpha (PIP4K2A) gene is located on chromosome 10p12. This region has been implicated in both diseases by linkage, and PIP4K2A directly by association. Given linkage evidence in the Irish Study of High Density Schizophrenia Families (ISHDSF) to a region including 10p12, we performed an association study between genetic variants at PIP4K2A and disease. No association was detected through single-marker or haplotype analysis of the whole sample. However, stratification into families positive and negative for the ISHDSF schizophrenia high-risk haplotype (HRH) in the DTNBP1 gene and re-analysis for linkage showed reduced amplitude of the 10p12 linkage peak in the DTNBP1 HRH positive families. Association analysis of the stratified sample showed a trend toward association of PIP4K2A SNPs rs1417374 and rs1409395 with schizophrenia in the DTNBP1 HRH positive families. Despite this apparent paradox, our data may therefore suggest involvement of PIP4K2A in schizophrenia in those families for whom genetic variation in DTNBP1 appears also to be a risk factor. This trend appears to arise from under-transmission of common alleles to female cases. Follow-up association analysis in a large Irish schizophrenia case-control control sample (ICCSS) showed significant association with disease of a haplotype comprising these same SNPs rs1417374-rs1409395, again more so in affected females, and in cases with negative family history of the disease. This study supports a minor role for PIP4K2A in schizophrenia etiology in the Irish population. (C) 2009 Wiley-Liss, Inc.

DNA hypermethylation and DNA hypomethylation is present at different loci in chronic kidney disease

Genetic risk factors for chronic kidney disease (CKD) are being identified through international collaborations. By comparison, epigenetic risk factors for CKD have only recently been considered using population-based approaches. DNA methylation is a major epigenetic modification that is associated with complex diseases, so we investigated methylome-wide loci for association with CKD. A total of 485,577 unique features were evaluated in 255 individuals with CKD (cases) and 152 individuals without evidence of renal disease (controls). Following stringent quality control, raw data were quantile normalized and β values calculated to reflect the methylation status at each site. The difference in methylation status was evaluated between cases and controls with resultant P values adjusted for multiple testing. Genes with significantly increased and decreased levels of DNA methylation were considered for biological relevance by functional enrichment analysis using KEGG pathways in Partek Genomics Suite. Twenty-three genes, where more than one CpG per loci was identified with Padjusted < 10−8, demonstrated significant methylation changes associated with CKD and additional support for these associated loci was sought from published literature. Strong biological candidates for CKD that showed statistically significant differential methylation include CUX1, ELMO1, FKBP5, INHBA-AS1, PTPRN2, and PRKAG2 genes; several genes are differentially methylated in kidney tissue and RNA-seq supports a functional role for differential methylation in ELMO1 and PRKAG2 genes. This study reports the largest, most comprehensive, genome-wide quantitative evaluation of DNA methylation for association with CKD. Evidence confirming methylation sites influence development of CKD would stimulate research to identify epigenetic therapies that might be clinically useful for CKD.

Genomic markers for decision making:What is preventing us from using markers?

The advent of novel genomic technologies that enable the evaluation of genomic alterations on a genome-wide scale has significantly altered the field of genomic marker research in solid tumors. Researchers have moved away from the traditional model of identifying a particular genomic alteration and evaluating the association between this finding and a clinical outcome measure to a new approach involving the identification and measurement of multiple genomic markers simultaneously within clinical studies. This in turn has presented additional challenges in considering the use of genomic markers in oncology, such as clinical study design, reproducibility and interpretation and reporting of results. This Review will explore these challenges, focusing on microarray-based gene-expression profiling, and highlights some common failings in study design that have impacted on the use of putative genomic markers in the clinic. Despite these rapid technological advances there is still a paucity of genomic markers in routine clinical use at present. A rational and focused approach to the evaluation and validation of genomic markers is needed, whereby analytically validated markers are investigated in clinical studies that are adequately powered and have pre-defined patient populations and study endpoints. Furthermore, novel adaptive clinical trial designs, incorporating putative genomic markers into prospective clinical trials, will enable the evaluation of these markers in a rigorous and timely fashion. Such approaches have the potential to facilitate the implementation of such markers into routine clinical practice and consequently enable the rational and tailored use of cancer therapies for individual patients. © 2010 Macmillan Publishers Limited. All rights reserved.

HUMAN GENOME VARIATIONS AND EVOLUTION WITH A FOCUS ON THE ANALYSIS OF TRANSPOSABLE ELEMENTS

Genome sequence varies in numerous ways among individuals although the gross architecture is fixed for all humans. Retrotransposons create one of the most abundant structural variants in the human genome and are divided in many families, with certain members in some families, e.g., L1, Alu, SVA, and HERV-K, remaining active for transposition. Along with other types of genomic variants, retrotransponson-derived variants contribute to the whole spectrum of genome variants in humans. With the advancement of sequencing techniques, many human genomes are being sequenced at the individual level, fueling the comparative research on these variants among individuals. In this thesis, the evolution and functional impact of structural variations is examined primarily focusing on retrotransposons in the context of human evolution. The thesis comprises of three different studies on the topics that are presented in three data chapters. First, the recent evolution of all human specific AluYb members, representing the second most active subfamily of Alus, was tracked to identify their source/master copy using a novel approach. All human-specific AluYb elements from the reference genome were extracted, aligned with one another to construct clusters of similar copies and each cluster was analyzed to generate the evolutionary relationship between the members of the cluster. The approach resulted in identification of one major driver copy of all human specific Yb8 and the source copy of the Yb9 lineage. Three new subfamilies within the AluYb family – Yb8a1, Yb10 and Yb11 were also identified, with Yb11 being the youngest and most polymorphic. Second, an attempt to construct a relation between transposable elements (TEs) and tandem repeats (TRs) was made at a genome-wide scale for the first time. Upon sequence comparison, positional cross-checking and other relevant analyses, it was observed that over 20% of all TRs are derived from TEs. This result established the first connection between these two types of repetitive elements, and extends our appreciation for the impact of TEs on genomes. Furthermore, only 6% of these TE-derived TRs follow the already postulated initiation and expansion mechanisms, suggesting that the others are likely to follow a yet-unidentified mechanism. Third, by taking a combination of multiple computational approaches involving all types of genetic variations published so far including transposable elements, the first whole genome sequence of the most recent common ancestor of all modern human populations that diverged into different populations around 125,000-100,000 years ago was constructed. The study shows that the current reference genome sequence is 8.89 million base pairs larger than our common ancestor’s genome, contributed by a whole spectrum of genetic mechanisms. The use of this ancestral reference genome to facilitate the analysis of personal genomes was demonstrated using an example genome and more insightful recent evolutionary analyses involving the Neanderthal genome. The three data chapters presented in this thesis conclude that the tandem repeats and transposable elements are not two entirely distinctly isolated elements as over 20% TRs are actually derived from TEs. Certain subfamilies of TEs themselves are still evolving with the generation of newer subfamilies. The evolutionary analyses of all TEs along with other genomic variants helped to construct the genome sequence of the most recent common ancestor to all modern human populations which provides a better alternative to human reference genome and can be a useful resource for the study of personal genomics, population genetics, human and primate evolution.

Study of histone H3 lysine 56 deacetylation in saccharomyces cerevisiae

Chez la levure Saccharomyces cerevisiae, l'acétylation de l'histone H3 sur la lysine 56 (H3K56ac) est présente sur les histones néo-synthétisées déposées derrière les fourches de réplication et est essentielle pour préserver la viabilité cellulaire en réponse au dommage à l'ADN. La désacétylation d'H3K56 sur l'ensemble du génome catalysée par Hst3 et Hst4 et a lieu en phase G2 ou M. H3K56ac est une lame à double tranchant. L'absence d'H3K56ac rend les cellules sensibles aux dommages à l'ADN. En revanche, un excès d'acétylation d'H3K56 dans un mutant hst3Δ hst4Δ a des conséquences encore plus sévères tels que la thermo-sensibilité, l'hypersensibilité aux agents génotoxiques, l'instabilité génomique ainsi qu'une courte durée de vie réplicative. Les désacétylases Hst3 et Hst4 sont étroitement régulées au cours du cycle cellulaire afin de permettre à l'H3K56ac d'exercer son rôle en réponse aux dommages à l'ADN tout en évitant les conséquences néfastes de l'hyperacétylation d'H3K56. Dans cette thèse, nous avons identifié la machinerie moléculaire responsable de la dégradation de Hst3. De plus, nous avons exploré les raisons pour lesquelles l'absence de désacétylation donne lieu aux phénotypes du mutant hst3Δ hst4Δ. Au chapitre 2, nous démontrons que la dégradation d'Hst3 peut être complétée avant l'anaphase. Ceci suggère que la désacétylation de H3K56 a lieu durant une courte fenêtre du cycle cellulaire se situant entre la complétion de la phase S et la métaphase. De plus, nous avons identifié deux sites de phosphorylation d'Hst3 par la kinase cycline-dépendante 1 (Cdk1) et démontré que ces évènements de phosphorylation conduisent à la dégradation d'Hst3 in vivo. Nous avons aussi démontré que l'ubiquityltransférase Cdc34 et l'ubiquitine ligase SCFCdc4 sont requises pour la dégradation d'Hst3. Finalement, nous avons montré que la phosphorylation d'Hst3 par la kinase mitotique Clb2-Cdk1 peut directement entraîner l'ubiquitylation d'Hst3 par SCFCdc4 in vitro. Au chapitre 3, nous avons étudié les mécanismes moléculaires sous-jacents à la sensibilité extrême du mutant hst3Δ hst4Δ aux agents qui endommagent l'ADN. Nous avons établi qu'en raison de la présence anormale d'H3K56ac devant les fourches de réplication, le mutant hst3Δ hst4Δ exhibe une forte perte de viabilité lorsqu'exposé au méthyl méthanesulfonate (MMS) durant un seul passage à travers la phase S. Nous avons aussi découvert que, malgré le fait que le point de contrôle de réponse aux dommages à l'ADN est activé normalement dans le mutant hst3Δ hst4Δ, ce mutant est incapable de compléter la réplication de l'ADN et d'inactiver le point de contrôle pour une longue période de temps après exposition transitoire au MMS. L'ensemble de nos résultats suggère que les lésions à l'ADN induites par le MMS dans le mutant hst3Δ hst4Δ causent une forte perte de viabilité parce que ce mutant est incapable de compléter la réplication de l'ADN après une exposition transitoire au MMS. Dans la deuxième section du chapitre 3, nous avons employé une approche génétique afin d'identifier de nouveaux mécanismes de suppression de deux phénotypes prononcés du mutant hst3Δ hst4Δ. Nous avons découvert que la délétion de plusieurs gènes impliqués dans la formation de frontières entre l'hétérochromatine et de l'euchromatine atténue les phénotypes du mutant hst3Δ hst4Δ sans réduire l'hyperacétylation d'H3K56. Nos résultats indiquent aussi que l'abondante acétylation de l'histone H4 sur la lysine 16 (H4K16ac) est néfaste au mutant hst3Δ hst4Δ. Ce résultat suggère un lien génétique intriguant entre l'acétylation d'H3K56 et celle d'H4K16. L'existence de ce lien était jusqu'à présent inconnu. Nous avons identifié un groupe de suppresseurs spontanés où H3K56ac est indétectable, mais la majorité de nos suppresseurs ne montrent aucune réduction flagrante d'H3K56ac ou d'H4 K16ac par rapport aux niveaux observés dans le mutant hst3Δ hst4Δ. Une étude plus approfondie de ce groupe de suppresseurs est susceptible de mener à la découverte de nouveaux mécanismes génétiques ou épigénétiques permettant d'éviter les conséquences catastrophiques de l'hyperacétylation d'H3K56 chez le mutant hst3Δ hst4Δ. En résumé, cette thèse identifie la machinerie moléculaire responsable de la dégradation d'Hst3 (une désacétylase d'H3K56) durant une fenêtre de temps situées entre la fin de la phase S et la métaphase. Nos résultats permettent aussi d'expliquer pourquoi la dégradation d'Hst3 précède le début de la phase S durant laquelle l'acétylation d'H3K56 s'accumule derrière les fourches de réplication afin d'exercer son rôle de mécanisme de défense contre le dommage à l'ADN. De plus, nous avons identifié plusieurs suppresseurs qui permettent de contourner le rôle important d'Hst3 et Hst4 en réponse au dommage à l'ADN. Plusieurs suppresseurs révèlent un lien génétique inattendu entre deux formes abondantes d'acétylation des histones chez Saccharomyces cerevisiae, soit H3K56ac et H4K16ac.

Genomic architecture of sickle cell disease clinical variation in children from West Africa : a case-control study design

Contexte : L’anémie falciforme ou drépanocytose est un problème de santé important, particulièrement pour les patients d’origine africaine. La variation phénotypique de l’anémie falciforme est problématique pour le suivi et le traitement des patients. L’architecture génomique responsable de cette variabilité est peu connue. Principe : Mieux saisir la contribution génétique de la variation clinique de cette maladie facilitera l’identification des patients à risque de développer des phénotypes sévères, ainsi que l’adaptation des soins. Objectifs : L’objectif général de cette thèse est de combler les lacunes relatives aux connaissances sur l’épidémiologie génomique de l’anémie falciforme à l’aide d’une cohorte issue au Bénin. Les objectifs spécifiques sont les suivants : 1) caractériser les profils d’expressions génomiques associés à la sévérité de l’anémie falciforme ; 2) identifier des biomarqueurs de la sévérité de l’anémie falciforme ; 3) identifier la régulation génétique des variations transcriptionelles ; 4) identifier des interactions statistiques entre le génotype et le niveau de sévérité associé à l’expression ; 5) identifier des cibles de médicaments pour améliorer l’état des patients atteints d’anémie falciforme. Méthode : Une étude cas-témoins de 250 patients et 61 frères et soeurs non-atteints a été menée au Centre de Prise en charge Médical Intégré du Nourrisson et de la Femme Enceinte atteints de Drépanocytose, au Bénin entre février et décembre 2010. Résultats : Notre analyse a montré que des profils d’expressions sont associés avec la sévérité de l’anémie falciforme. Ces profils sont enrichis de génes des voies biologiques qui contribuent à la progression de la maladie : l’activation plaquettaire, les lymphocytes B, le stress, l’inflammation et la prolifération cellulaire. Des biomarqueurs transcriptionnels ont permis de distinguer les patients ayant des niveaux de sévérité clinique différents. La régulation génétique de la variation de l’expression des gènes a été démontrée et des interactions ont été identifiées. Sur la base de ces résultats génétiques, des cibles de médicaments sont proposées. Conclusion: Ce travail de thèse permet de mieux comprendre l’impact de la génomique sur la sévérité de l’anémie falciforme et ouvre des perspectives de développement de traitements ciblés pour améliorer les soins offerts aux patients.

Functional studies of the deubiquitinating enzymes USP19, USP4 and UCH-L1

El marcaje de proteínas con ubiquitina, conocido como ubiquitinación, cumple diferentes funciones que incluyen la regulación de varios procesos celulares, tales como: la degradación de proteínas por medio del proteosoma, la reparación del ADN, la señalización mediada por receptores de membrana, y la endocitosis, entre otras (1). Las moléculas de ubiquitina pueden ser removidas de sus sustratos gracias a la acción de un gran grupo de proteasas, llamadas enzimas deubiquitinizantes (DUBs) (2). Las DUBs son esenciales para la manutención de la homeostasis de la ubiquitina y para la regulación del estado de ubiquitinación de diferentes sustratos. El gran número y la diversidad de DUBs descritas refleja tanto su especificidad como su utilización para regular un amplio espectro de sustratos y vías celulares. Aunque muchas DUBs han sido estudiadas a profundidad, actualmente se desconocen los sustratos y las funciones biológicas de la mayoría de ellas. En este trabajo se investigaron las funciones de las DUBs: USP19, USP4 y UCH-L1. Utilizando varias técnicas de biología molecular y celular se encontró que: i) USP19 es regulada por las ubiquitin ligasas SIAH1 y SIAH2 ii) USP19 es importante para regular HIF-1α, un factor de transcripción clave en la respuesta celular a hipoxia, iii) USP4 interactúa con el proteosoma, iv) La quimera mCherry-UCH-L1 reproduce parcialmente los fenotipos que nuestro grupo ha descrito previamente al usar otros constructos de la misma enzima, y v) UCH-L1 promueve la internalización de la bacteria Yersinia pseudotuberculosis.

Búsqueda de selección en el polimorfismo 677C>T (c.665C>T) del gen de la metilentetrahidrofolato reductasa (MTHFR) en una población Colombiana

Se realizó un estudio genético – poblacional en dos grupos etarios de población colombiana con la finalidad de evaluar las diferencias genéticas relacionadas con el polimorfismo MTHFR 677CT en busca de eventos genéticos que soporten la persistencia de este polimorfismo en la especie humana debido que este ha sido asociado con múltiples enfermedades. De esta manera se genotipificaron los individuos, se analizaron los genotipos, frecuencias alélicas y se realizaron diferentes pruebas genéticas-poblacionales. Contrario a lo observado en poblaciones Colombianas revisadas se identificó la ausencia del Equilibrio Hardy-Weinberg en el grupo de los niños y estructuras poblacionales entre los adultos lo que sugiere diferentes historias demográficas y culturales entre estos dos grupos poblacionales al tiempo, lo que soporta la hipótesis de un evento de selección sobre el polimorfismo en nuestra población. De igual manera nuestros datos fueron analizados junto con estudios previos a nivel nacional y mundial lo cual sustenta que el posible evento selectivo es debido a que el aporte de ácido fólico se ha incrementado durante las últimas dos décadas como consecuencia de las campañas de fortificación de las harinas y suplementación a las embarazadas con ácido fólico, por lo tanto aquí se propone un modelo de selección que se ajusta a los datos encontrados en este trabajo se establece una relación entre los patrones nutricionales de la especie humana a través de la historia que explica las diferencias en frecuencias de este polimorfismo a nivel espacial y temporal.