Molecular genetics of rare hereditary diseases via genome-wide and integrated approaches

ABSTRACTThe Online Mendelian Inheritance in Man database (OMIM) reports about 3000 Mendelian diseases of known causal gene and about 2000 that remain to be mapped. These cases are often difficult to solve because of the rareness of the disease, the structure of the family (too big or too small) or the heterogeneity of the phenotype. The goal of this thesis is to explore the current genetic tools, before the advent of ultra high throughput sequencing, and integrate them in the attempt to map the genes behind the four studied cases. In this framework we have studied a small family with a recessive disease, a modifier gene for the penetrance of a dominant mutation, a large extended family with a cardiac phenotype and clinical and/or allelic heterogeneity and we have molecularly analyzed a balanced chromosomal translocation.RESUMELa base de données des maladies à transmission mendélienne, Online Mendelian Inheritance in Man (OMIM), contient environ 3000 affections à caractère mendélien pour lesquelles le gène responsable est connu et environ 2000 qui restent à élucider.Les cas restant à résoudre sont souvent difficiles soit par le caractère intrinsèquement rare de ces maladies soit à cause de difficultés structurelles (famille trop petite ou trop étendue) ou hétérogénéité du phénotype ou génétique. Cette thèse s'inscrit avant l'arrivée des nouveaux outils de séquençage à haut débit. Son but est d'explorer les outils génétiques actuels, et de les intégrer pour trouver les gènes impliqués dans quatre cas représentant chacun une situation génétique différente : nous avons étudié une famille de quatre individus avec une transmission récessive, recherché un gène modificateur de la pénétrance de mutations dominantes, étudié une famille étendue présentant un phénotype cardiaque cliniquement et/ou allèliquement hétérogène et nous avons fait l'analyse moléculaire d'une translocation chromosomique balancée.

Genome-Wide Analysis of Salicylate and Dibenzofuran Metabolism in Sphingomonas Wittichii RW1.

Sphingomonas wittichii RW1 is a bacterium isolated for its ability to degrade the xenobiotic compounds dibenzodioxin and dibenzofuran (DBF). A number of genes involved in DBF degradation have been previously characterized, such as the dxn cluster, dbfB, and the electron transfer components fdx1, fdx3, and redA2. Here we use a combination of whole genome transcriptome analysis and transposon library screening to characterize RW1 catabolic and other genes implicated in the reaction to or degradation of DBF. To detect differentially expressed genes upon exposure to DBF, we applied three different growth exposure experiments, using either short DBF exposures to actively growing cells or growing them with DBF as sole carbon and energy source. Genome-wide gene expression was examined using a custom-made microarray. In addition, proportional abundance determination of transposon insertions in RW1 libraries grown on salicylate or DBF by ultra-high throughput sequencing was used to infer genes whose interruption caused a fitness loss for growth on DBF. Expression patterns showed that batch and chemostat growth conditions, and short or long exposure of cells to DBF produced very different responses. Numerous other uncharacterized catabolic gene clusters putatively involved in aromatic compound metabolism increased expression in response to DBF. In addition, only very few transposon insertions completely abolished growth on DBF. Some of those (e.g., in dxnA1) were expected, whereas others (in a gene cluster for phenylacetate degradation) were not. Both transcriptomic data and transposon screening suggest operation of multiple redundant and parallel aromatic pathways, depending on DBF exposure. In addition, increased expression of other non-catabolic genes suggests that during initial exposure, S. wittichii RW1 perceives DBF as a stressor, whereas after longer exposure, the compound is recognized as a carbon source and metabolized using several pathways in parallel.

Peroxisome proliferator activated receptor BETA/DELTA as a central player in the skin response to ultra-violets radiation

Previous studies in the lab of Dr. Liliane Michalik, have shown thai the nuclear hormone receptor Peroxisome Proliferator Activated Receptor beta/delta (PPARß/ö) is an important regulator of skin homeostasis, being involved in the regulation of keratinocyte differentiation, inflammation, apoptosis, arid mouse skin wound healing. Studies of PPARß/ö knock out mice have suggested a possible role for this receptor in cancer. However, contradictory observations of the role for PPARß/ö on tumor growth have been published, depending on cellular contexts and biological models. Given the controversial role of PPARß/ö in skin carcinoma development, the main aim of this PhD work has been to further explore the implication of PPARß/ö in skin response to UV and skin tumor growth. This PhD dissertation is divided in four chapters. The first chapter describes the core part of the project, where I explored the changes in miRNA expression in the skin upon chronic UV irradiation of PPARß/ö wild type and knock-out mice. This analysis shed light on a miRNA- PPARß/ö signature and also predicted thai miR-21-3p (previously named miR-21*) is a key regulator of the PPARß/ö-dependent UV response in the pre-lesiona! skin. Using mice acutely UV-irradiated, ! further demonstrated that miR-21-3p is indirectly regulated by PPARß/ö through activation of Transforming Growth Factor (TGFß)-1 under UV exposure. I also show that miR-21-3p is deregulated in human cutaneous squamous celi carcinoma. In cultured keratinocytes, application of a miR-21 -3p mimic oligonucleotide sequence leads to the regulation of lipid metabolism-related pathway. In the second chapter, I demonstrate that the usage of an mRNA/miRNA combined bioinformatics analysis leads to the discovery of important pathways involved in the PPARß/ö-miRNA response of the skin to chronic UV irradiation, indeed, I validated angiogenesis and lipid metabolism as important functions regulated by PPARß/ö in this context. In the third chapter, we demonstrate that PPARß/5 knockout mice have decreased cutaneous squamous cell carcinomas incidence compared to wild type mice and that PPARß/5 directly activates the cSrc kinase gene. In the last chapter, we review novel insights into PPAR functions in keratinocytes and liver, with emphasis on PPARß/ö but also on PPARa. In summary, this PhD study shows that i) PPARß/5 is able to regulate biological function through regulation of miRNAs, and specifically through miR-21-3p, the passenger miRNA of the oncomiR miR-21, and that ii) the PPARß/5-dependent skin response to UV involves the regulation of angiogenesis and lipid metabolism. Furthermore, the bioinformatics study highlights the relevance of performing integrated mRNA and miRNA genome-wide studies in order to better screen mRNAs and/or miRNAs of interest in the biological context of diseases. - Des études préalables dans le laboratoire du Dr. Liliane Michalik ont démontré que le récepteur nucléaire PPARß/5 est un régulateur important de l'homéostasie de la peau, étant impliqué dans la régulation de la différenciation des keratinocytes, dans l'inflammation, dans l'apoptose et dans la cicatrisation de la peau chez !a souris. L'étude de souris knock-out pour le gène PPARß/5, ont suggérées un rôle possible de ce récepteur dans le cancer. Cependant, des observations opposées ont été publiées suggérant un rôle pro- ou anti- cancer selon le tissue impliqué et le type- cellulaire. En considérant cette controverse autour du rôle de PPARß/5 dans le développement des cancers de la peau, le but principal de mon projet de recherche aura été d'approfondir l'exploration du rôle de PPARß/5 dans la réponse de la peau aux UVs et dans le développement du cancer. Cette dissertation de thèse est divisée en quatre parties. Une première partie, représentant le coeur de mon travail de recherche, décrit la découverte de l'implication des microRNAs (rniRNAs) dans la réponse aux UVs de PPARß/ö et plus spécifiquement l'implication du miRNA miR- 21 -3p (précédemment nommé miR-21*). En étudiant un modèle de souris irradiées de manière aigüe aux UVs, nous montrons que ia régulation de miR-21-3p est PPARß/ö-däpenaante et que cette régulation à lieu par l'intermédiaire du facteur de transcription TGFß-1. Dans des cultures de keratinocytes Humains, la transfecticn d'une séquence oligonucléotidique similaire à celle de miR-21-3p (mimic), montre l'implication de rniR-21-3p dans des fonctions importantes pour le développement des cancers telles que le métabolisme des lipides. Dans un second chapitre, nous montrons que l'usage d'une méthode bioinformatique combinant l'expression des ARN messagers et des miRNAs permet de mettre en évidence des fonctions biologiques importantes lors de ia réponse de PPARß/ö à l'irradiation chronique. L'angiogenèse, le stress oxydatif et le métabolisme des lipides font partie de ces fonctions régulées par PPARß/5 dans la peau irradiée aux UVs. Nous mettons également en évidence la régulation du gène LpcatS par PPARß/5 dans la peau irradiée aux UV ainsi que dans des keratinocytes humains suggérant un rôle pour PPARß/5 dans le remodelage des lipides membranaires. Dans une troisième partie, nous établissons un lien entre la régulation de l'oncogène Src et l'activation de PPARß/5 dans les carcinomes spinocellulaires de la peau. Finalement dans un quatrième chapitre, nous faisons une revue des dernières recherches portées sur le rôle de PPARß/5 et de PPARa dans le foie et ia peau. En résumé ce projet de thèse représente un avancement pour la recherche sur rimplication de PPARß/5 dans la réponse aux UVs de la peau. Pour la première fois, un lien est établi entre ce facteur de transcription et la régulation de microRNAs dans le cadre du carcinome spinocellulare. Jusqu'alors resté dans l'ombre de rniR-21-5p, miR-21-3p est en fait fortement augmenté à la fois dans un modèle de souris d'irradiation aux UVs ainsi que dans ie carcinome spinocellulare chez i'humain. De nouvelles fonctions biologiques pour PPARß/5 ont été également mises en évidence dans ce travail, comme la régulation de l'angiogenèse ou du métabolisme des lipides dans Sa peau. De plus cette dissertation valorise l'intérêt d'une association entre le travail de laboratoire et celui de la bioinformatique.

Fast analysis of doping agents in urine by ultra-high-pressure liquid chromatography-quadrupole time-of-flight mass spectrometry I. Screening analysis.

The general strategy to perform anti-doping analyses of urine samples starts with the screening for a wide range of compounds. This step should be fast, generic and able to detect any sample that may contain a prohibited substance while avoiding false negatives and reducing false positive results. The experiments presented in this work were based on ultra-high-pressure liquid chromatography coupled to hybrid quadrupole time-of-flight mass spectrometry. Thanks to the high sensitivity of the method, urine samples could be diluted 2-fold prior to injection. One hundred and three forbidden substances from various classes (such as stimulants, diuretics, narcotics, anti-estrogens) were analysed on a C(18) reversed-phase column in two gradients of 9min (including two 3min equilibration periods) for positive and negative electrospray ionisation and detected in the MS full scan mode. The automatic identification of analytes was based on retention time and mass accuracy, with an automated tool for peak picking. The method was validated according to the International Standard for Laboratories described in the World Anti-Doping Code and was selective enough to comply with the World Anti-Doping Agency recommendations. In addition, the matrix effect on MS response was measured on all investigated analytes spiked in urine samples. The limits of detection ranged from 1 to 500ng/mL, allowing the identification of all tested compounds in urine. When a sample was reported positive during the screening, a fast additional pre-confirmatory step was performed to reduce the number of confirmatory analyses.

Genome-wide association study of metabolic traits reveals novel gene-metabolite-disease links.

Metabolic traits are molecular phenotypes that can drive clinical phenotypes and may predict disease progression. Here, we report results from a metabolome- and genome-wide association study on (1)H-NMR urine metabolic profiles. The study was conducted within an untargeted approach, employing a novel method for compound identification. From our discovery cohort of 835 Caucasian individuals who participated in the CoLaus study, we identified 139 suggestively significant (P<5×10(-8)) and independent associations between single nucleotide polymorphisms (SNP) and metabolome features. Fifty-six of these associations replicated in the TasteSensomics cohort, comprising 601 individuals from São Paulo of vastly diverse ethnic background. They correspond to eleven gene-metabolite associations, six of which had been previously identified in the urine metabolome and three in the serum metabolome. Our key novel findings are the associations of two SNPs with NMR spectral signatures pointing to fucose (rs492602, P = 6.9×10(-44)) and lysine (rs8101881, P = 1.2×10(-33)), respectively. Fine-mapping of the first locus pinpointed the FUT2 gene, which encodes a fucosyltransferase enzyme and has previously been associated with Crohn's disease. This implicates fucose as a potential prognostic disease marker, for which there is already published evidence from a mouse model. The second SNP lies within the SLC7A9 gene, rare mutations of which have been linked to severe kidney damage. The replication of previous associations and our new discoveries demonstrate the potential of untargeted metabolomics GWAS to robustly identify molecular disease markers.

Genetic variation in IL28B is associated with chronic hepatitis C and treatment failure: a genome-wide association study.

BACKGROUND & AIMS: Hepatitis C virus (HCV) induces chronic infection in 50% to 80% of infected persons; approximately 50% of these do not respond to therapy. We performed a genome-wide association study to screen for host genetic determinants of HCV persistence and response to therapy. METHODS: The analysis included 1362 individuals: 1015 with chronic hepatitis C and 347 who spontaneously cleared the virus (448 were coinfected with human immunodeficiency virus [HIV]). Responses to pegylated interferon alfa and ribavirin were assessed in 465 individuals. Associations between more than 500,000 single nucleotide polymorphisms (SNPs) and outcomes were assessed by multivariate logistic regression. RESULTS: Chronic hepatitis C was associated with SNPs in the IL28B locus, which encodes the antiviral cytokine interferon lambda. The rs8099917 minor allele was associated with progression to chronic HCV infection (odds ratio [OR], 2.31; 95% confidence interval [CI], 1.74-3.06; P = 6.07 x 10(-9)). The association was observed in HCV mono-infected (OR, 2.49; 95% CI, 1.64-3.79; P = 1.96 x 10(-5)) and HCV/HIV coinfected individuals (OR, 2.16; 95% CI, 1.47-3.18; P = 8.24 x 10(-5)). rs8099917 was also associated with failure to respond to therapy (OR, 5.19; 95% CI, 2.90-9.30; P = 3.11 x 10(-8)), with the strongest effects in patients with HCV genotype 1 or 4. This risk allele was identified in 24% of individuals with spontaneous HCV clearance, 32% of chronically infected patients who responded to therapy, and 58% who did not respond (P = 3.2 x 10(-10)). Resequencing of IL28B identified distinct haplotypes that were associated with the clinical phenotype. CONCLUSIONS: The association of the IL28B locus with natural and treatment-associated control of HCV indicates the importance of innate immunity and interferon lambda in the pathogenesis of HCV infection.

Ultra high performance supercritical fluid chromatography coupled with tandem mass spectrometry for screening of doping agents. II: Analysis of biological samples.

The potential and applicability of UHPSFC-MS/MS for anti-doping screening in urine samples were tested for the first time. For this purpose, a group of 110 doping agents with diverse physicochemical properties was analyzed using two separation techniques, namely UHPLC-MS/MS and UHPSFC-MS/MS in both ESI+ and ESI- modes. The two approaches were compared in terms of selectivity, sensitivity, linearity and matrix effects. As expected, very diverse retentions and selectivities were obtained in UHPLC and UHPSFC, proving a good complementarity of these analytical strategies. In both conditions, acceptable peak shapes and MS detection capabilities were obtained within 7min analysis time, enabling the application of these two methods for screening purposes. Method sensitivity was found comparable for 46% of tested compounds, while higher sensitivity was observed for 21% of tested compounds in UHPLC-MS/MS and for 32% in UHPSFC-MS/MS. The latter demonstrated a lower susceptibility to matrix effects, which were mostly observed as signal suppression. In the case of UHPLC-MS/MS, more serious matrix effects were observed, leading typically to signal enhancement and the matrix effect was also concentration dependent, i.e., more significant matrix effects occurred at the lowest concentrations.

Quantification of clenbuterol at trace level in human urine by ultra-high pressure liquid chromatography-tandem mass spectrometry.

Clenbuterol is a β2 agonist agent with anabolic properties given by the increase in the muscular mass in parallel to the decrease of the body fat. For this reason, the use of clenbuterol is forbidden by the World Anti-Doping Agency (WADA) in the practice of sport. This compound is of particular interest for anti-doping authorities and WADA-accredited laboratories due to the recent reporting of risk of unintentional doping following the eating of meat contaminated with traces of clenbuterol in some countries. In this work, the development and the validation of an ultra-high pressure liquid chromatography coupled to electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) method for the quantification of clenbuterol in human urine is described. The analyte was extracted from urine samples by liquid-liquid extraction (LLE) in basic conditions using tert butyl-methyl ether (TBME) and analyzed by UHPLC-MS/MS with a linear gradient of acetonitrile in 9min only. The simple and rapid method presented here was validated in compliance with authority guidelines and showed a limit of quantification at 5pg/mL and a linearity range from 5pg/mL to 300pg/mL. Good trueness (85.8-105%), repeatability (5.7-10.6% RSD) and intermediate precision (5.9-14.9% RSD) results were obtained. The method was then applied to real samples from eighteen volunteers collecting urines after single oral doses administration (1, 5 and 10μg) of clenbuterol-enriched yogurts.

Quality control and conduct of genome-wide association meta-analyses.

Rigorous organization and quality control (QC) are necessary to facilitate successful genome-wide association meta-analyses (GWAMAs) of statistics aggregated across multiple genome-wide association studies. This protocol provides guidelines for (i) organizational aspects of GWAMAs, and for (ii) QC at the study file level, the meta-level across studies and the meta-analysis output level. Real-world examples highlight issues experienced and solutions developed by the GIANT Consortium that has conducted meta-analyses including data from 125 studies comprising more than 330,000 individuals. We provide a general protocol for conducting GWAMAs and carrying out QC to minimize errors and to guarantee maximum use of the data. We also include details for the use of a powerful and flexible software package called EasyQC. Precise timings will be greatly influenced by consortium size. For consortia of comparable size to the GIANT Consortium, this protocol takes a minimum of about 10 months to complete.

Genome-wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America.

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome-wide next-generation sequencing of DNA pools ('pool-seq') from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome-wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well-known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.

Genome-wide linkage in a highly consanguineous pedigree reveals two novel loci on chromosome 7 for non-syndromic familial Premature Ovarian Failure.

BACKGROUND: The human condition known as Premature Ovarian Failure (POF) is characterized by loss of ovarian function before the age of 40. A majority of POF cases are sporadic, but 10-15% are familial, suggesting a genetic origin of the disease. Although several causal mutations have been identified, the etiology of POF is still unknown for about 90% of the patients.¦METHODOLOGY/PRINCIPAL FINDINGS: We report a genome-wide linkage and homozygosity analysis in one large consanguineous Middle-Eastern POF-affected family presenting an autosomal recessive pattern of inheritance. We identified two regions with a LOD(max) of 3.26 on chromosome 7p21.1-15.3 and 7q21.3-22.2, which are supported as candidate regions by homozygosity mapping. Sequencing of the coding exons and known regulatory sequences of three candidate genes (DLX5, DLX6 and DSS1) included within the largest region did not reveal any causal mutations.¦CONCLUSIONS/SIGNIFICANCE: We detect two novel POF-associated loci on human chromosome 7, opening the way to the identification of new genes involved in the control of ovarian development and function.

Permeability estimates from poro-elastic dispersion analyses of sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

Multiplex ultra-performance liquid chromatography-tandem mass spectrometry method for simultaneous quantification in human plasma of fluconazole, itraconazole, hydroxyitraconazole, posaconazole, voriconazole, voriconazole-N-oxide, anidulafungin, and caspofungin.

Therapeutic drug monitoring (TDM) may contribute to optimizing the efficacy and safety of antifungal therapy because of the large variability in drug pharmacokinetics. Rapid, sensitive, and selective laboratory methods are needed for efficient TDM. Quantification of several antifungals in a single analytical run may best fulfill these requirements. We therefore developed a multiplex ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) method requiring 100 μl of plasma for simultaneous quantification within 7 min of fluconazole, itraconazole, hydroxyitraconazole, posaconazole, voriconazole, voriconazole-N-oxide, caspofungin, and anidulafungin. Protein precipitation with acetonitrile was used in a single extraction procedure for eight analytes. After reverse-phase chromatographic separation, antifungals were quantified by electrospray ionization-triple-quadrupole mass spectrometry by selected reaction monitoring detection using the positive mode. Deuterated isotopic compounds of azole antifungals were used as internal standards. The method was validated based on FDA recommendations, including assessment of extraction yields, matrix effect variability (<9.2%), and analytical recovery (80.1 to 107%). The method is sensitive (lower limits of azole quantification, 0.01 to 0.1 μg/ml; those of echinocandin quantification, 0.06 to 0.1 μg/ml), accurate (intra- and interassay biases of -9.9 to +5% and -4.0 to +8.8%, respectively), and precise (intra- and interassay coefficients of variation of 1.2 to 11.1% and 1.2 to 8.9%, respectively) over clinical concentration ranges (upper limits of quantification, 5 to 50 μg/ml). Thus, we developed a simple, rapid, and robust multiplex UPLC-MS/MS assay for simultaneous quantification of plasma concentrations of six antifungals and two metabolites. This offers, by optimized and cost-effective lab resource utilization, an efficient tool for daily routine TDM aimed at maximizing the real-time efficacy and safety of different recommended single-drug antifungal regimens and combination salvage therapies, as well as a tool for clinical research.

A mega-analysis of genome-wide association studies for major depressive disorder.

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 × 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 × 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.