Transcriptional profiling in Saccharomyces cerevisiae relevant for predicting alachlor mechanisms of toxicity

Alachlor has been a commonly applied herbicide and is a substance of ecotoxicological concern. The present study aims to identify molecular biomarkers in the eukaryotic model Saccharomyces cerevisiae that can be used to predict potential cytotoxic effects of alachlor, while providing new mechanistic clues with possible relevance for experimentally less accessible eukaryotes. It focuses on genome-wide expression profiling in a yeast population in response to two exposure scenarios exerting effects from slight to moderate magnitude at phenotypic level. In particular, 100 and 264 genes, respectively, were found as differentially expressed on a 2-h exposure of yeast cells to the lowest observed effect concentration (110 mg/L) and the 20% inhibitory concentration (200 mg/L) of alachlor, in comparison with cells not exposed to the herbicide. The datasets of alachlor-responsive genes showed functional enrichment in diverse metabolic, transmembrane transport, cell defense, and detoxification categories. In general, the modifications in transcript levels of selected candidate biomarkers, assessed by quantitative reverse transcriptase polymerase chain reaction, confirmed the microarray data and varied consistently with the growth inhibitory effects of alachlor. Approximately 16% of the proteins encoded by alachlor-differentially expressed genes were found to share significant homology with proteins from ecologically relevant eukaryotic species. The biological relevance of these results is discussed in relation to new insights into the potential adverse effects of alachlor in health of organisms from ecosystems, particularly in worst-case situations such as accidental spills or careless storage, usage, and disposal.

FROM GENES TO GENOMES: LOCAL ADAPTATION AND ADAPTIVE POTENTIAL IN TWO ARCTIC SEABIRDS

By investigating the mechanisms underlying the evolution and the maintenance of local adaptations we can help predict how species will adapt to future environmental change. In this thesis I investigate local adaptation and adaptive potential in thick-billed and common murres (Uria lomvia and U. aalge), two arctic seabirds of international conservation concern. Thanks to the recent development of new genomic methods, I address three major themes that are relevant for both the development of evolutionary theory and conservation: 1) the role of gene flow in the origin and maintenance of adaptation; 2) levels and distribution of standing genetic variation, and their contribution to adaptive potential; and 3) the genomic mechanisms maintaining an adaptive dimorphism within a single interbreeding population. First, I review the literature on genomics of local adaptation with gene flow and find that adaptation can be maintained despite gene flow, that gene flow itself can promote adaptation, and that genetic architecture is important in the origin and maintenance of local adaptations. Second, I genotype genome-wide markers and toll-like receptor genes (TLRs) to investigate local adaptation and adaptive potential in thick-billed murres. Thick-billed murres do not show signatures of local adaptation to their breeding grounds, but outlier loci group birds according to their non-breeding distributions, suggesting that selection and/or demographic connectivity in the winter may explain patterns of differentiation in this species. Genetic variation at TLRs does not decrease with increasing latitude as predicted, but tests of selection and measures of genetic diversity suggest differences in local selective regimes at most genes. Thick-billed murres show high levels of standing genetic variation and their adaptive potential will mostly depend on rate and magnitude of environmental change. Finally, I improve and annotate the assembly of the highly heterozygous genome of the thick-billed murre. Using this assembly as a reference, I perform whole genome analyses to investigate the genomic basis of an adaptive dimorphism in Atlantic common murres. I show for the first time that a 60 kb complex copy number variant in a non-coding region maintains differences in plumage and cold adaptation despite high gene flow.

Genomics to elucidate the molecular basis of calcific aortic valve disease

Le rétrécissement valvulaire aortique (RVA) est causé par une calcification et une fibrose progressive de la valve aortique. Le risque de développer la maladie augmente avec l’âge. À cause de l’augmentation de l’espérance de vie, le RVA est devenu un problème de santé publique. Le RVA est fatal en absence de traitement médical. Actuellement, la chirurgie est le seul traitement pour le stade sévère de la maladie, mais près de 50% des individus avec RVA n’y sont pas éligibles, principalement due à la présence de comorbidités. Plusieurs processus biologiques ont été associés à la maladie, mais les voies moléculaires spécifiques et les gènes impliqués dans le développement et la progression du RVA ne sont pas connus. Il est donc urgent de découvrir les gènes de susceptibilité pour le RVA afin d’identifier les personnes à risque ainsi que les biomarqueurs et les cibles thérapeutiques pouvant mener au développement de médicaments pour inverser ou limiter la progression de la maladie. L’objectif de cette thèse de doctorat était d’identifier la base moléculaire du RVA. Des approches modernes en génomique, incluant l’étude de gènes candidats et le criblage génomique par association (GWAS), ont été réalisées à l’aide de collections d’ADN provenant d’un grand nombre de patients bien caractérisés pour le RVA. Des études complémentaires en transciptomique ont comparé le profil d’expression global des gènes entre des valves calcifiées et non-calcifiées à l’aide de biopuces à ADN et de séquençage de l’ARN. Une première étude a identifié des variations dans le gène NOTCH1 et suggère pour la première fois la présence d’un polymorphisme commun dans ce gène conférant une susceptibilité au RVA. La deuxième étude a combiné par méta-analyse deux GWAS de patients provenant de la ville de Québec et Paris (France) aux données transcriptomiques. Cette étude de génomique intégrative a confirmé le rôle de RUNX2 dans le RVA et a permis l’identification d’un nouveau gène de susceptibilité, CACNA1C. Les troisième et quatrième études sur l’expression des gènes ont permis de mieux comprendre les bases moléculaires de la calcification des valves aortiques bicuspides et ainsi d’identifier de nouvelles cibles thérapeutiques pour le RVA. Les données générées par ce projet sont la base de futures découvertes importantes qui permettront d’améliorer les options de traitement et la qualité de vie des patients atteints du RVA.

Cartographie des cassures bicaténaires du remodelage chromatinien du spermatide et développement des outils techniques associés.

Résumé : La phase haploïde de la spermatogenèse (spermiogenèse) est caractérisée par une modification importante de la structure de la chromatine et un changement de la topologie de l’ADN du spermatide. Les mécanismes par lesquels ce changement se produit ainsi que les protéines impliquées ne sont pas encore complètement élucidés. Mes travaux ont permis d’établir la présence de cassures bicaténaires transitoires pendant ce remodelage par l’essai des comètes et l’électrophorèse en champ pulsé. En procédant à des immunofluorescences sur coupes de tissus et en utilisant un extrait nucléaire hautement actif, la présence de topoisomérases ainsi que de marqueurs de systèmes de réparation a été confirmée. Les protéines de réparation identifiées font partie de systèmes sujets à l’erreur, donc cette refonte structurale de la chromatine pourrait être génétiquement instable et expliquer le biais paternel observé pour les mutations de novo dans de récentes études impliquant des criblages à haut débit. Une technique permettant l’immunocapture spécifique des cassures bicaténaires a été développée et appliquée sur des spermatides murins représentant différentes étapes de différenciation. Les résultats de séquençage à haut débit ont montré que les cassures bicaténaires (hotspots) de la spermiogenèse se produisent en majorité dans l’ADN intergénique, notamment dans les séquences LINE1, l’ADN satellite et les répétions simples. Les hotspots contiennent aussi des motifs de liaisons des protéines des familles FOX et PRDM, dont les fonctions sont entre autres de lier et remodeler localement la chromatine condensée. Aussi, le motif de liaison de la protéine BRCA1 se trouve enrichi dans les hotspots de cassures bicaténaires. Celle-ci agit entre autres dans la réparation de l’ADN par jonction terminale non-homologue (NHEJ) et dans la réparation des adduits ADN-topoisomérase. De façon remarquable, le motif de reconnaissance de la protéine SPO11, impliquée dans la formation des cassures méiotiques, a été enrichi dans les hotspots, ce qui suggère que la machinerie méiotique serait aussi utilisée pendant la spermiogenèse pour la formation des cassures. Enfin, bien que les hotspots se localisent plutôt dans les séquences intergéniques, les gènes ciblés sont impliqués dans le développement du cerveau et des neurones. Ces résultats sont en accord avec l’origine majoritairement paternelle observée des mutations de novo associées aux troubles du spectre de l’autisme et de la schizophrénie et leur augmentation avec l’âge du père. Puisque les processus du remodelage de la chromatine des spermatides sont conservés dans l’évolution, ces résultats suggèrent que le remodelage de la chromatine de la spermiogenèse représente un mécanisme additionnel contribuant à la formation de mutations de novo, expliquant le biais paternel observé pour certains types de mutations.

Dissecting the Heat-Stress Tolerance QTLome in Durum Wheat

Heat stress negatively affects wheat performance during its entire cycle, particularly during the reproductive stage. In view of the climate change and the prediction of a continued increase in temperature in the new future, it is urgent to concentrate efforts to discover novel genetic sources able to improve the resilience of wheat to heat stress. In this direction, this study addressed two different experiments in durum wheat to identify novel QTLs suitable to be applied in marker-assisted selection for heat tolerance. Chlorophyll fluorescence (ChlF) is a valuable indicator of plant response to environmental changes allowing a detailed assessment of PSII activity in view of its non-invasive measurement and high-throughput phenotyping. In the first study (Chapter 2), the Light-Induced Fluorescence Transient (LIFT) method was used to access ChlF data to map QTLs for ChlF-related traits during the vegetative growth stage in durum wheat under heat stress condition. Our results provide evidence that LIFT consistently measures ChlF at the level of high-throughput phenotyping combined with high accuracy which is required for Genome-Wide Association Study (GWAS) aimed at identifying genomic regions affecting PSII activity. The 50 QTLs identified for ChlF-related traits under heat stress mostly clustered into five chromosomes hotspots unrelated to phenology, a feature that makes these QTLs a valuable asset for marker-assisted breeding programs across different latitudes. In the second study (Chapter 3), a set of 183 accessions suitable for GWAS, was exposed to optimal and high temperature during two crop seasons under field conditions. Important agronomic traits were evaluated in order to identify valuable QTLs for GY and its components. The GWAS analysis identified several QTLs in the single years as well as in the joint analysis. From the total QTLs identified, 13 QTL clusters can be highlighted to be affecting heat tolerance across different years and/or different traits.

Mapping new non-genetic dependencies in malignant pleural mesothelioma

Malignant Pleural Mesothelioma (MPM) is a very aggressive cancer whose incidence is growing worldwide. MPM escapes the classical models of carcinogenesis and lacks a distinctive genetic fingerprint, keeping obscure the molecular events that lead to tumorigenesis. This severely impacts on the limited therapeutic options and on the lack of specific biomarkers, concurring to make MPM one of the deadliest cancers. Here we combined a functional genome-wide loss of function CRISPR/Cas9 screening with patients’ transcriptomic and clinical data, to identify genes essential for MPM progression. Besides, we explored the role of non-coding RNAs to MPM progression by analysing gene expression profiles and clinical data from the MESO-TCGA dataset. We identified TRIM28 and the lncRNA LINC00941 as new vulnerabilities of MPM, associated with disease aggressiveness and bad outcome of patients. TRIM28 is a multi-domain protein involved in many processes, including transcription regulation. We showed that TRIM28 silencing impairs MPM cells’ growth and clonogenicity by blocking cells in mitosis. RNA-seq profiling showed that TRIM28 loss abolished the expression of major mitotic players. Our data suggest that TRIM28 is part of the B-MYB/FOXM1-MuvB complex that specifically drives the activation of mitotic genes, keeping the time of mitosis. In parallel, we found LINC00941 as strongly associated with reduced survival probability in MPM patients. LINC00941 KD profoundly reduced MPM cells’ growth, migration and invasion. This is accompanied by changes in morphology, cytoskeleton organization and cell-cell adhesion properties. RNA-seq profiling showed that LINC00941 KD impacts crucial functions of MPM, including HIF1α signalling. Collectively these data provided new insights into MPM biology and demonstrated that the integration of functional screening with patients’ clinical data is a powerful tool to highlight new non-genetic cancer dependencies that associate to a bad outcome in vivo, paving the way to new MPM-oriented targeted strategies and prognostic tools to improve patients risk-based stratification.

An Update On The Pharmacogenetics Of Treating Hypertension.

Hypertension is a leading cause of cardiovascular mortality, but only one third of patients achieve blood pressure goals despite antihypertensive therapy. Genetic polymorphisms may partially account for the interindividual variability and abnormal response to antihypertensive drugs. Candidate gene and genome-wide approaches have identified common genetic variants associated with response to antihypertensive drugs. However, there is no currently available pharmacogenetic test to guide hypertension treatment in clinical practice. In this review, we aimed to summarize the recent findings on pharmacogenetics of the most commonly used antihypertensive drugs in clinical practice, including diuretics, angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers, beta-blockers and calcium channel blockers. Notably, only a small percentage of the genetic variability on response to antihypertensive drugs has been explained, and the vast majority of the genetic variants associated with antihypertensives efficacy and toxicity remains to be identified. Despite some genetic variants with evidence of association with the variable response related to these most commonly used antihypertensive drug classes, further replication is needed to confirm these associations in different populations. Further studies on epigenetics and regulatory pathways involved in the responsiveness to antihypertensive drugs might provide a deeper understanding of the physiology of hypertension, which may favor the identification of new targets for hypertension treatment and genetic predictors of antihypertensive response.Journal of Human Hypertension advance online publication, 28 August 2014; doi:10.1038/jhh.2014.76.

Genetic Variants of Diabetes Risk and Incident Cardiovascular Events in Chronic Coronary Artery Disease

Objective: To determine whether information from genetic risk variants for diabetes is associated with cardiovascular events incidence. Methods: From the about 30 known genes associated with diabetes, we genotyped single-nucleotide polymorphisms at the 10 loci most associated with type-2 diabetes in 425 subjects from the MASS-II Study, a randomized study in patients with multi-vessel coronary artery disease. The combined genetic information was evaluated by number of risk alleles for diabetes. Performance of genetic models relative to major cardiovascular events incidence was analyzed through Kaplan-Meier curve comparison and Cox Hazard Models and the discriminatory ability of models was assessed for cardiovascular events by calculating the area under the ROC curve. Results: Genetic information was able to predict 5-year incidence of major cardiovascular events and overall-mortality in non-diabetic individuals, even after adjustment for potential confounders including fasting glycemia. Non-diabetic individuals with high genetic risk had a similar incidence of events then diabetic individuals (cumulative hazard of 33.0 versus 35.1% of diabetic subjects). The addition of combined genetic information to clinical predictors significantly improved the AUC for cardiovascular events incidence (AUC = 0.641 versus 0.610). Conclusions: Combined information of genetic variants for diabetes risk is associated to major cardiovascular events incidence, including overall mortality, in non-diabetic individuals with coronary artery disease.

Heteroduplex formation and S1 digestion for mapping alternative splicing sites

The identification of alternatively spliced transcripts has contributed to a better comprehension of developmental mechanisms, tissue-specific physiological processes and human diseases. Polymerase chain reaction amplification of alternatively spliced variants commonly leads to the formation of heteroduplexes as a result of base pairing involving exons common between the two variants. S1 nuclease cleaves single-stranded loops of heteroduplexes and also nicks the opposite DNA strand. In order to establish a strategy for mapping alternative splice-prone sites in the whole transcriptome, we developed a method combining the formation of heteroduplexes between 2 distinct splicing variants and S1 nuclease digestion. For 20 consensuses identified here using this methodology, 5 revealed a conserved splice site after inspection of the cDNA alignment against the human genome (exact splice sites). For 8 other consensuses, conserved splice sites were mapped at 2 to 30 bp from the border, called proximal splice sites; for the other 7 consensuses, conserved splice sites were mapped at 40 to 800 bp, called distal splice sites. These latter cases showed a nonspecific activity of S1 nuclease in digesting double-strand DNA. From the 20 consensuses identified here, 5 were selected for reverse transcription-polymerase chain reaction validation, confirming the splice sites. These data showed the potential of the strategy in mapping splice sites. However, the lack of specificity of the S1 nuclease enzyme is a significant obstacle that impedes the use of this strategy in large-scale studies.

Familial partial epilepsy with variable foci: A new partial epilepsy syndrome with suggestion of linkage to chromosome 2

Familial partial epilepsy with variable foci (FPEVF) joins the recently recognized group of inherited partial epilepsies. We describe an Australian family with 10 individuals with partial seizures over four generations. Detailed electroclinical studies were performed on all affected and 17 clinically unaffected family members. The striking finding was that the clinical features of the seizures and interictal electroencephalographic foci differed among family members and included frontal, temporal, occipital, and centroparietal seizures. Mean age of seizure onset was 13 years (range, 0.75-43 years). Two individuals without seizures had epileptiform abnormalities on electroencephalographic studies. Penetrance of seizures was 62%. A genome-wide search failed to demonstrate definitive linkage, but a suggestion of linkage was found on chromosome 2q with a LOD score of 2.74 at recombination fraction of zero with the marker D2S133. FPEVF differs from the other inherited partial epilepsies where partial seizures in different family members are clinically similar. The inherited nature of this new syndrome may be overlooked because of relatively low penetrance and because of the variability in age at onset and electroclinical features between affected family members.

De novo copy number variants identify new genes and loci in isolated sporadic tetralogy of Fallot

Tetralogy of Fallot (TOF), the most common severe congenital heart malformation, occurs sporadically, without other anomaly, and from unknown cause in 70% of cases. Through a genome-wide survey of 114 subjects with TOF and their unaffected parents, we identified 11 de novo copy number variants (CNVs) that were absent or extremely rare (<0.1%) in 2,265 controls. We then examined a second, independent TOF cohort (n = 398) for additional CNVs at these loci. We identified CNVs at chromosome 1q21.1 in 1% (5/512, P = 0.0002, OR = 22.3) of nonsyndromic sporadic TOF cases. We also identified recurrent CNVs at 3p25.1, 7p21.3 and 22q11.2. CNVs in a single subject with TOF occurred at six loci, two that encode known (NOTCH1, JAG1) disease-associated genes. Our findings predict that at least 10% (4.5-15.5%, 95% confidence interval) of sporadic nonsyndromic TOF cases result from de novo CNVs and suggest that mutations within these loci might be etiologic in other cases of TOF.

A transposon toolkit for gene transfer and mutagenesis in protozoan parasites

Protozoan parasites affect millions of people around the world. Treatment and control of these diseases are complicated partly due to the intricate biology of these organisms. The interactions of species of Plasmodium, Leishmania and trypanosomes with their hosts are mediated by an unusual control of gene expression that is not fully understood. The availability of the genome sequence of these protozoa sets the stage for using more comprehensive, genome-wide strategies to study gene function. Transposons are effective tools for the systematic introduction of genetic alterations and different transposition systems have been adapted to study gene function in these human pathogens. A mariner transposon toolkit for use in vivo or in vitro in Leishmania parasites has been developed and can be used in a variety of applications. These modified mariner elements not only permit the inactivation of genes, but also mediate the rescue of translational gene fusions, bringing a major contribution to the investigation of Leishmania gene function. The piggyBac and Tn5 transposons have also been shown to mobilize across Plasmodium spp. genomes circumventing the current limitations in the genetic manipulation of these organisms.

Familial varieties of primary aldosteronism

1. Improved approaches to screening and diagnosis have revealed primary aldosteronism (PAL) to be much more common than previously thought, with most patients normokalaemic. The spectrum of this disorder has been further broadened by the study of familial varieties. 2. Familial hyperaldosteronism type I (FH-I) is a glucocorticoid-remediable form of PAL caused by the inheritance of an adrenocorticotrophic hormone (ACTH)- regulated, hybrid CYP11B1/CYP11B2 gene. Diagnosis has been greatly facilitated by the advent of genetic testing. The severity of hypertension varies widely in FH-I, even among members of the same family, and has demonstrated relationships with gender, degree of biochemical disturbance and hybrid gene crossover point position. Hormone day curve studies show that the hybrid gene dominates over wild-type CYP11B2 in terms of aldosterone regulation. This may be due, in part, to a defect in wild-type CYP11B2-induced aldosterone production. Control of hypertension in FH-I requires only partial suppression of ACTH and much smaller glucocorticoid doses than previously recommended. 3. Familial hyperaldosteronism type II (FH-II) is not glucocorticoid remediable and is not associated with the hybrid gene mutation. Familial hyperaldosteronism type II is clinically, biochemically and morphologically indistinguishable from apparently non-familial PAL. Linkage studies in one informative family did not show segregation of FH-II with the CYP11B2, AT1 or MEN1 genes, but a genome-wide search has revealed linkage with a locus in chromosome 7. As has already occurred in FH-I, elucidation of causative mutations is likely to facilitate earlier detection of PAL.

Characterization of mouse profilaggrin: Evidence for nuclear engulfment and translocation of the profilaggrin B-domain during epidermal differentiation

Filaggrin is a keratin filament associated protein that is expressed in granular layer keratinocytes and derived by sequential proteolysis from a polyprotein precursor termed profilaggrin. Depending on the species, each profilaggrin molecule contains between 10 and 20 filaggrin subunits organized as tandem repeats with a calcium-binding domain at the N-terminal end. We now report the characterization of the complete mouse gene. The structural organization of the mouse gene is identical to the human profilaggrin gene and consists of three exons with a 4 kb intron within the 5' noncoding region and a 1.7 kb intron separating the sequences encoding the calcium-binding EF-hand motifs. A processed pseudogene was found embedded within the second intron. The third and largest exon encodes the second EF-hand, a basic domain (designated the B-domain) followed by 12 filaggrin repeats and a unique C-terminal tail domain. A polyclonal anti-body raised against the conceptually translated sequence of the B-domain specifically stained keratohyalin granules and colocalized with a filaggrin antibody in granular layer cells. In upper granular layer cells, B-domain containing keratohyalin granules were in close apposition to the nucleus and, in some cells, appeared to be completely engulfed by the nucleus. In transition layer cells, B-domain staining was evident in the nucleus whereas filaggrin staining remained cytoplasmic. Nuclear staining of the B-domain was also observed in primary mouse keratinocytes induced to differentiate. This study has also revealed significant sequence homology between the mouse and human promoter sequences and in the calcium-binding domain but the remainder of the protein-coding region shows substantial divergence.