High-throughput sequencing technologies are improving in quality, capacity, and costs, providing versatile applications in DNA and RNA research. For small genomes or fraction of larger genomes, DNA samples can be mixed and loaded together on the same sequencing track. This so-called multiplexing approach relies on a specific DNA tag, index, or barcode that is attached to the sequencing or amplification primer and hence accompanies every read. After sequencing, each sample read is identified on the basis of the respective barcode sequence. Alterations of DNA barcodes during synthesis, primer ligation, DNA amplification, or sequencing may lead to incorrect sample identification unless the error is revealed and corrected. This can be accomplished by implementing error correcting algorithms and codes. This barcoding strategy increases the total number of correctly identified samples, thus improving overall sequencing efficiency. Two popular sets of error-correcting codes are Hamming codes and codes based on the Levenshtein distance. Levenshtein-based codes operate only on words of known length. Since a DNA sequence with an embedded barcode is essentially one continuous long word, application of the classical Levenshtein algorithm is problematic. In this thesis we demonstrate the decreased error correction capability of Levenshtein-based codes in a DNA context and suggest an adaptation of Levenshtein-based codes that is proven of efficiently correcting nucleotide errors in DNA sequences. In our adaptation, we take any DNA context into account and impose more strict rules for the selection of barcode sets. In simulations we show the superior error correction capability of the new method compared to traditional Levenshtein and Hamming based codes in the presence of multiple errors. We present an adaptation of Levenshtein-based codes to DNA contexts capable of guaranteed correction of a pre-defined number of insertion, deletion, and substitution mutations. Our improved method is additionally capable of correcting on average more random mutations than traditional Levenshtein-based or Hamming codes. As part of this work we prepared software for the flexible generation of DNA codes based on our new approach. To adapt codes to specific experimental conditions, the user can customize sequence filtering, the number of correctable mutations and barcode length for highest performance. However, not every platform is susceptible to a large number of both indel and substitution errors. The Illumina “Sequencing by Synthesis” platform shows a very large number of substitution errors as well as a very specific shift of the read that results in inserted and deleted bases at the 5’-end and the 3’-end (which we call phaseshifts). We argue in this scenario that the application of Sequence-Levenshtein-based codes is not efficient because it aims for a category of errors that barely occurs on this platform, which reduces the code size needlessly. As a solution, we propose the “Phaseshift distance” that exclusively supports the correction of substitutions and phaseshifts. Additionally, we enable the correction of arbitrary combinations of substitution and phaseshift errors. Thus, we address the lopsided number of substitutions compared to phaseshifts on the Illumina platform. To compare codes based on the Phaseshift distance to Hamming Codes as well as codes based on the Sequence-Levenshtein distance, we simulated an experimental scenario based on the error pattern we identified on the Illumina platform. Furthermore, we generated a large number of different sets of DNA barcodes using the Phaseshift distance and compared codes of different lengths and error correction capabilities. We found that codes based on the Phaseshift distance can correct a number of errors comparable to codes based on the Sequence-Levenshtein distance while offering the number of DNA barcodes comparable to Hamming codes. Thus, codes based on the Phaseshift distance show a higher efficiency in the targeted scenario. In some cases (e.g., with PacBio SMRT in Continuous Long Read mode), the position of the barcode and DNA context is not well defined. Many reads start inside the genomic insert so that adjacent primers might be missed. The matter is further complicated by coincidental similarities between barcode sequences and reference DNA. Therefore, a robust strategy is required in order to detect barcoded reads and avoid a large number of false positives or negatives. For mass inference problems such as this one, false discovery rate (FDR) methods are powerful and balanced solutions. Since existing FDR methods cannot be applied to this particular problem, we present an adapted FDR method that is suitable for the detection of barcoded reads as well as suggest possible improvements.

Piscine francisellosis in an infectious emerging bacterial disease that affects several marine and fresh water fish species worldwide, including farmed salmon, wild and farmed cod, farmed tilapia and several ornamental species, for which no commercial treatment or vaccine exists. During 2011 and the first semester of 2012, chronic episodes of moderate to high levels of mortality with nonspecific clinical signs, and widespread multifocal white nodules as the most consistent gross pathological lesion were experienced by farmed tilapia fingerlings at two different locations in Northern Europe. In this study such outbreaks of granulomatous disease were diagnosed as francisellosis with a genus-specific PCR, and 10 new isolates of the bacterium including the one named STIR-GUS-F2f7, were recovered on a new selective “cysteine blood-tilapia” agar and cysteine heart agar with bovine haemoglobin. Ultrastructural observations of the pathogen in Nile tilapia (O. niloticus) tissues suggested the secretion of outer membrane vesicles (OMVs) by the bacterial cells during infection in these fish. This represented the first documented report of isolation of pathogenic Francisella strains from tilapia in Europe. The phenotypic characterisation indicated that isolates recovered were able to metabolise dextrin, N-acetyl-D glucosamine, D-fructose, α-D-glucose, D-mannose, methyl pyruvate, acetic acid, α-keto butyric acid, L-alaninamide, L-alanine, L-alanylglycine, L-asparagine, L-glutamic acid, L-proline, L-serine, L-threonine, inosine, uridine, glycerol, D L-α-glycerol phosphate, glucose-1-phosphate and glucose-6-phosphate. The predominant structural fatty acids of the isolates were 24:1 (20.3%), 18:1n-9 (16.9%), 24:0 (13.1%) 14:0 (10.9%), 22:0 (7.8%), 16:0 (7.6%) and 18:0 (5.5%). Anti-microbial resistance analyses indicated that STIR-GUS-F2f7 was susceptible to neomycin, novobiocin, amikacin, ciprofloxacin, imipenem, gatifloxacin, meropenem, tobramycin, nitrofurantoin, and levofloxacin using the quantitative broth micro-dilution method, while the qualitative disc diffusion method indicated susceptibility to enrofloxacin, kanamycin, gentamicin, tetracycline, oxytetracycline, florfenicol, oxolinic acid and streptomycin. The use of the following housekeeping genes: mdh, dnaA, mutS, 16SrRNA-ITS-23SrRNA, prfB putA rpoA, rpoB and tpiA indicated 100% similarity with other isolates belonging to the subspecies F. noatunensis orientalis (Fno). Koch’s postulates were successfully fulfilled by establishing an intraperitoneal injection (IP) challenge model with STIR-GUS-F2f7 in Nile tilapia. Moreover, the challenge model was used to investigate the susceptibility of 3 genetic groups of tilapia to STIR-GUS-F2f7. The lowest amount of bacteria required to cause mortality was 12 CFU/ml and this was seen as early as only 24 hours post infection in the red Nile tilapia and in the wild type after 26 days, no mortalities were seen in the species O. mossambicus with this dose. The mortality in red O. niloticus was significantly higher than that of the other two tilapia groups when 12 and 120 CFU/fish were injected. It was also observed that when a dose of 1200 CFU/ml was used, the mortality in O. niloticus wild type was significantly lower than that of the other two tilapia groups and no differences were seen among the 3 groups when the highest dose (1.2 x105 CFU/fish) was used. The median lethal dose (LD50) of O. niloticus wild type was the most stable during the experiment (values around 104 CFU/ml) and the highest of the three groups after day 25 post infection. At the end of the experiment (day 45) the LD50 was 30 CFU/ml in the red Nile tilapia, 2.3x104 CFU/ml for the wild type and 3.3x102 CFU/ml for O. mossambicus. This pattern, where the LD50 of the red tilapia was lower than that of the other two groups, was observed during the whole experiment. The outcomes of these experiments suggested that the red Nile tilapia family appeared to be the most susceptible while the wild type Nile tilapia family the most resistant. The complete genome of STIR-GUS-F2f7 was sequenced using next generation sequencing (NGS) Illumina Hi-Seq platform™, and the annotation of the assembled genome predicted 1970 protein coding sequences and 63 non-coding rRNA sequences distributed in 328 sub-systems. The taxonomy of the species Francisella noatunensis was revised using genomic-derived parameters form STIR-GUS-F2f7 and other strains in combination with a polyphasic approach that included ecologic, chemotaxonomic and phenotypic analyses. The results indicated that STIR-GUS-F2f7 and all the other strains from warm water fish represent a new bacterial species for which the name Francisella orientalis was assigned. Moreover the description of F. noatunensis was emended and the creation of a new subspecies within this taxon i.e. Francisella noatunensis subsp. chilense was proposed. The results of this study led to the development of a highly efficacious vaccine to protect tilapia against francisellosis.

Le cancer épithélial de l’ovaire (CEO) est le cancer gynécologique le plus létal. Plus de 70% des patientes diagnostiquées avec une tumeur de stade avancé rechutent suite aux traitements chimiothérapeutiques de première ligne, la survie à cinq ans étant ainsi très faible. Afin de mieux comprendre l’évolution de la maladie, nous avons recherché de nouveaux gènes, responsables de l’initiation et de la progression du CEO. Précédemment, des lignées cellulaires ont été dérivées à partir de la tumeur primaire et récurrente et/ou d’ascites de trois patientes. Le séquençage de l’ARN de ces lignées par la technologie de séquençage de nouvelle génération (TSNG) nous a permis d’identifier des mutations ponctuelles qui pourraient nous indiquer des gènes dérégulés dans le CEO. La TSNG est un bon outil qui permet d’identifier et de cribler à grande échelle des mutations. Nous avons sélectionné PLEC1, SCRIB, NCOR2, SEMA6C, IKBKB, GLCE et ITGAE comme gènes candidats présentant des mutations dans nos lignées et ayant une relation fonctionnelle avérée avec le cancer. Étant donné que la TSNG est une technique à taux de fiabilité limité, nous avons validé ces mutations par séquençage Sanger. Ensuite, nous avons étudié l’effet de ces mutations sur la structure protéique et l’expression de PLEC1, de SCRIB et de SEMA6C. Seules certaines mutations dans les gènes PLEC1, SCRIB et SEMA6C ont pu être confirmées. PLEC1 et SCRIB sont deux protéines d’échafaudage dont la mutation, rapportée dans plusieurs cancers, pourrait induire des changements de leurs conformations et affecter leurs interactions et leurs fonctions. Les conséquences de ces mutations sur la tumorigenèse de l’ovaire devront être étudiées. / Epithelial ovarian cancer (EOC) is the most lethal gynecological cancer. Over 70% of the patients diagnosed with advanced stage of cancer relapse following first-line chemotherapy treatments; consequently the five-year survival is very low. To better understand the evolution of the disease, our aim was to identify new genes responsible for the initiation and progression of EOC. Previously, cell lines derived from solid tumors or ascites were developed from the primary and recurrent tumor or ascites of three patients. RNA sequencing of these cell lines by next-generation sequencing technology (NGST) allowed us to identify mutations that might point to genes whose deregulation is important in EOC. Mutations were detected in PLEC1, SCRIB, NCOR2, SEMA6C, IKBKB, GLCE and ITGAE. We selected these genes for further studies as they have previously been identified as being associated with cancer. First, we validated these mutations by Sanger sequencing in order to determine the concordance with NGST data. Secondly, we studied the impact of the validated mutations on protein structure and gene expression. Only certain mutations in PLEC1, SCRIB and SEMA6C were confirmed. Of interest, PLEC1 and SCRIB are two scaffold proteins, where mutations have been reported in several cancers and, possibly leading to changes in their conformation and thereby affecting their interactions and functions. The consequences of these mutations on ovarian tumorigenesis remain to be determined.

Chez les plantes, le génome plastidique est continuellement exposé à divers stress mutagènes, tels l’oxydation des bases et le blocage des fourches de réplication. Étonnamment, malgré ces menaces, le génome du plastide est reconnu pour être très stable, sa stabilité dépassant même celle du génome nucléaire. Néanmoins, les mécanismes de réparation de l’ADN et du maintien de la stabilité du génome plastidique sont encore peu connus. Afin de mieux comprendre ces processus, nous avons développé une approche, basée sur l’emploi de la ciprofloxacine, qui nous permet d’induire des bris d’ADN double-brins (DSBs) spécifiquement dans le génome des organelles. En criblant, à l’aide de ce composé, une collection de mutants d’Arabidopsis thaliana déficients pour des protéines du nucléoïde du plastide, nous avons identifié 16 gènes vraisemblablement impliqués dans le maintien de la stabilité génomique de cette organelle. Parmi ces gènes, ceux de la famille Whirly jouent un rôle primordial dans la protection du génome plastidique face aux réarrangements dépendants de séquences de microhomologie. Deux autres familles de gènes codant pour des protéines plastidiques, soit celle des polymérases de types-I et celle des recombinases, semblent davantage impliquées dans les mécanismes conservateurs de réparation des DSBs. Les relations épistatiques entre ces gènes et ceux des Whirly ont permis de définir les bases moléculaires des mécanismes de la réparation dépendante de microhomologies (MHMR) dans le plastide. Nous proposons également que ce type de mécanismes servirait en quelque sorte de roue de secours pour les mécanismes conservateurs de réparation. Finalement, un criblage non-biaisé, utilisant une collection de plus de 50,000 lignées mutantes d’Arabidopsis, a été réalisé. Ce criblage a permis d’établir un lien entre la stabilité génomique et le métabolisme des espèces réactives oxygénées (ROS). En effet, la plupart des gènes identifiés lors de ce criblage sont impliqués dans la photosynthèse et la détoxification des ROS. Globalement, notre étude a permis d’élargir notre compréhension des mécanismes du maintien de la stabilité génomique dans le plastide et de mieux comprendre l’importance de ces processus. / The plant plastidial genome is constantly threatened by many mutagenic stresses, such as base oxidation and replication fork stalling. Despite these threats, the plastid genome has long been known to be more stable than the nuclear genome, suggesting that alterations of its structure would have dramatic consequences on plant fitness. At the moment, little is known about the genes and the pathways allowing such conservation of the organelle genome sequences. To gain insight into these mechanisms, we developed an assay which uses ciprofloxacin, a gyrase inhibitor, to generate DNA double-strand breaks (DSBs) exclusively in plant organelles. By screening mutants deficient for proteins composing the plastid nucleoid on ciprofloxacin, we were able to identify 16 candidate genes, most likely involved in the repair of DSBs in plastid. Among these genes, those of the Whirly family of single-stranded DNA binding proteins are shown to be key factors in protecting the genome from error-prone microhomology mediated repair (MHMR). Two other family of proteins, the plastid type-I polymerases and the plastid recombinases, seem to be involved in the conservative repair pathways. The evaluation of the epistatic relationship between those two genes and the Whirly genes led us to define the molecular basis of MHMR and to propose that they might act as a backup system for conservative repair pathways. Finally, a non-biased screen, using 50,000 different insertion lines, allowed the identification of numerous genes that were already associated with ROS homeostasis, suggesting a link between DNA repair and ROS imbalance. Globally, our study shed light on the mechanisms that allow the maintenance of plastid genome, while explaining the importance of such conservation of the plastid genome.

Le but de ce projet était de développer des méthodes d'assemblage de novo dans le but d'assembler de petits génomes, principalement bactériens, à partir de données de séquençage de nouvelle-génération. Éventuellement, ces méthodes pourraient être appliquées à l'assemblage du génome de StachEndo, une Alpha-Protéobactérie inconnue endosymbiote de l'amibe Stachyamoeba lipophora. Suite à plusieurs analyses préliminaires, il fut observé que l’utilisation de lectures Illumina avec des assembleurs par graphe DeBruijn produisait les meilleurs résultats. Ces expériences ont également montré que les contigs produits à partir de différentes tailles de k-mères étaient complémentaires pour la finition des génomes. L’ajout de longues paires de lectures chevauchantes se montra essentiel pour la finition complète des grandes répétitions génomiques. Ces méthodes permirent d'assembler le génome de StachEndo (1,7 Mb). L'annotation de ce génome permis de montrer que StachEndo possède plusieurs caractéristiques inhabituelles chez les endosymbiotes. StachEndo constitue une espèce d'intérêt pour l'étude du développement endosymbiotique. / The goal of this project was to develop de novo genome assembly methods adapted to small genomes, especially bacterial, using next-generation sequencing data. Eventually, these methods could be used to assemble the genome of StachEndo, an unknown Alpha-Proteobacteria ensymbiont of the Stachyamoeba lipophora amoeba. Preliminary findings showed that the use of Illumina reads with DeBruijn graph assemblers yielded the best results. These experiments also showed that contigs produced with k-mers of various sizes were complementary in genome finishing assays. The addition of long-range paired-end reads proved necessary to fully close genomic assembly gaps. These methods made the assembly of StachEndo’s genome (1.7 Mb) possible. Through the annotation of StachEndo’s genes, several features that are unusal for endosymbionts were identified. StachEndo seems to be an interesting species for the study of endosymbiotic evolution.

L'hypothyroïdie congénitale par dysgénésie thyroïdienne (HCDT, ectopie dans plus de 80 %) a une prévalence de 1 cas sur 4000 naissances vivantes. L’HCDT est la conséquence d'une défaillance de la thyroïde embryonnaire à se différencier, à se maintenir ou à migrer vers sa localisation anatomique (partie antérieure du cou), qui aboutit à une absence totale de la thyroïde (athyréose) ou à une ectopie thyroïdienne (linguale ou sublinguale). Les HCDT sont principalement non-syndromiques (soit 98% des cas sont non-familiale), ont un taux de discordance de 92% chez les jumeaux monozygotes, et ont une prédominance féminine et ethnique (i.e., Caucasienne). La majorité des cas d’HCDT n’a pas de cause connue, mais est associée à un déficit sévère en hormones thyroïdiennes (hypothyroïdie). Des mutations germinales dans les facteurs de transcription liés à la thyroïde (NKX2.1, FOXE1, PAX8, NKX2.5) ont été identifiées dans seulement 3% des patients atteints d’HCDT sporadiques et l’analyse de liaisons exclue ces gènes dans les rares familles multiplex avec HCDT. Nous supposons que le manque de transmission familiale claire d’HCDT peut résulter de la nécessité d’au moins deux « hits » génétiques différents dans des gènes importants pour le développement thyroïdien. Pour répondre au mieux nos questions de recherche, nous avons utilisé deux approches différentes: 1) une approche gène candidat, FOXE1, seul gène impliqué dans l’ectopie dans le modèle murin et 2) une approche en utilisant les techniques de séquençage de nouvelle génération (NGS) afin de trouver des variants génétiques pouvant expliquer cette pathologie au sein d’une cohorte de patients avec HCDT. Pour la première approche, une étude cas-contrôles a été réalisée sur le promoteur de FOXE1. Il a récemment été découvert qu’une région du promoteur de FOXE1 est différentiellement méthylée au niveau de deux dinucléotides CpG consécutifs, définissant une zone cruciale de contrôle de l’expression de FOXE1. L’analyse d’association basée sur les haplotypes a révélé qu’un haplotype (Hap1: ACCCCCCdel1C) est associé avec le HCDT chez les Caucasiens (p = 5x10-03). Une réduction significative de l’activité luciférase est observée pour Hap1 (réduction de 68%, p<0.001) comparé au promoteur WT de FOXE1. Une réduction de 50% de l’expression de FOXE1 dans une lignée de cellules thyroïdienne humaine est suffisante pour réduire significativement la migration cellulaire (réduction de 55%, p<0.05). Un autre haplotype (Hap2: ACCCCCCC) est observé moins fréquemment chez les Afro-Américain comparés aux Caucasiens (p = 1.7x10-03) et Hap2 diminue l’activité luciférase (réduction de 26%, p<0.001). Deux haplotypes distincts sont trouvés fréquemment dans les contrôles Africains (Black-African descents). Le premier haplotype (Hap3: GTCCCAAC) est fréquent (30.2%) chez les contrôles Afro-Américains comparés aux contrôles Caucasiens (6.3%; p = 2.59 x 10-9) tandis que le second haplotype (Hap4: GTCCGCAC) est trouvé exclusivement chez les contrôles Afro-Américains (9.4%) et est absent chez les contrôles Caucasiens (P = 2.59 x 10-6). Pour la deuxième approche, le séquençage de l’exome de l’ADN leucocytaire entre les jumeaux MZ discordants n’a révélé aucune différence. D'où l'intérêt du projet de séquençage de l’ADN et l’ARN de thyroïdes ectopiques et orthotopiques dans lesquelles de l'expression monoallélique aléatoire dans a été observée, ce qui pourrait expliquer comment une mutation monoallélique peut avoir des conséquences pathogéniques. Finalement, le séquençage de l’exome d’une cohorte de 36 cas atteints d’HCDT a permis d’identifier de nouveaux variants probablement pathogéniques dans les gènes récurrents RYR3, SSPO, IKBKE et TNXB. Ces quatre gènes sont impliqués dans l’adhésion focale (jouant un rôle dans la migration cellulaire), suggérant un rôle direct dans les défauts de migration de la thyroïde. Les essais de migration montrent une forte diminution (au moins 60% à 5h) de la migration des cellules thyroïdiennes infectées par shRNA comparés au shCtrl dans 2 de ces gènes. Des zebrafish KO (-/- et +/-) pour ces nouveaux gènes seront réalisés afin d’évaluer leur impact sur l’embryologie de la thyroïde. / Congenital hypothyroidism by thyroid dysgenesis (CHTD) is a common disorder with prevalence at birth of 1 in 4000 live births. CHTD is the consequence of a failure of embryonic thyroid to differentiate or to migrate to its anatomical location (front of the neck), which leads to a total lack of thyroid (athyreosis) or an ectopic thyroid (lingual or sublingual). The most common category is ectopic thyroid diagnosis (up 85%). Most cases of CHTD have no known cause, but are associated with severe deficiency of thyroid hormones (hypothyroidism). The clinical diagnosis of hypothyroidism is usually possible only when permanent brain damage is already present. On the other hand, biochemical screening on the second day of life allows initiating replacement therapy from the second week of life, pre-empting severe intellectual deficit associated with the congenital hypothyroidism. Even with early treatment (an average of 9 days), loss of IQ, which is not exclusively due to the severity of hypothyroidism, can still be observed. Molecular markers may identify patients at risk for intellectual deficit (by e.g., genes involved in neuronal migration and the thyroid during development). These patients might benefit from early intervention to stimulate their neurocognitive development. Cases of CHTD are mainly non-syndromic and sporadic (in 98% of cases, there is no other affected in the family), have a discordant rate of 92% in monozygotic twins, and a female and ethnic (Caucasian) dominance. Germline mutations in thyroid-related transcription factors have been identified in only 3% of patients with sporadic CHTD, and linkage analysis has excluded these genes in rare multiplex families with CHTD. In addition, non-penetrating mutations among close relatives (for Nkx2.5 gene) suggest that modifying genes as germline variants de novo copy number (CNV) and / or somatic mutations are associated with CHTD. To respond to this research questions, we used two different approaches: 1) a candidate gene approach studying FOXE1, the only gene involved in ectopic thyroid in the mouse model and, 2) an approach using next generation sequencing techniques (NGS) to find genetic variants that could explain this pathology using a cohort of mostly sporadic CHTD. Variants and genes discovered by these two different approaches have been validated and their functional impact on the thyroid gland was evaluated by several experiments.

Dead wood plays an important role in forest ecosystems in the context of C dynamics, nutrient cycling, forest regeneration and biodiversity. Decaying wood sustains biodiversity by providing habitats and energy for fungi, bacteria, invertebrates, and many other organisms. Dead wood is resistant to decomposition and its decay is driven mainly by filamentous fungi. Community structure of wood- inhabiting fungi changes during decomposition, but the relationship between substrate quality and decomposer community is still poorly understood. This work studied fungal community composition with respect to tree species, stage of decay, volume and physico-chemical properties (such as pH, carbon and nitrogen content) of dead wood. Fungi were identified using next generation sequencing approaches - 454-pyrosequencing and Illumina MiSeq sequencing. Tree species, volume of dead wood (branches x logs) and stage of decay were the main variables affecting fungal community composition. Higher enzyme activities and content of fungal biomass indicate faster colonization of small branches than tree trunks by fungi. Fungal community composition, wood chemical properties and enzyme activities changed during decomposition. Both content of nitrogen and fungal biomass increased during decomposition. Enzyme activites peaked...

This work is part of a project focused on the study of the variability of human cytomegalovirus (HCMV) among clinical isolates with the aim to map the geographical distribution of HCMV genotypes, reveal the relationships between genotypes and the severity of HCMV-associated diseases, and identify regions in the HCMV genome with a potential for use as diagnostic and therapeutic targets. Attention was paid to the development of the methodology for the preparation of the material for next-generation sequencing (NGS) from HCMV clinical isolates and evaluation of the obtained sequencing data. Blood and urine samples collected from hematopoietic stem cell transplantat recipients and congenitally infected children were analyzed. Samples suitable for NGS were sequenced by the Illumina platform and sequences were created by de novo assembly followed by mapping assembly. Urine samples in comparison to blood samples had higher yield of material for NGS. Of the samples positive for HCMV DNA (7 of 50) after amplification in the cell cultures, only one sample had high purity of the viral DNA (98%) while six samples had purity of less than 7%. The sample containing 98% of the viral DNA was fully sequenced and the sequence was compared to the sequences of other clinical isolates from Belgium in 11 polymorphic...

This thesis is a part of the project aimed at sequencing the genome and transcriptome of the bank vole (Clethrionomys glareolus). The role of natural selection in the evolution of mitochondrial DNA (mtDNA) has been subject to much discussion; while some studies did not provide evidence that selection affected the phylogeography of the studied species, other considered adaptive evolution important. The bank vole is the key model we use to study the adaptation to climate change. As with other species, the phylogeography of the bank vole has been based on the variation of a small part of mtDNA. The goal of the thesis was to sequence the entire mitochondrial genome for representatives of all main mtDNA lineages of the bank vole using the Sanger and Illumina technologies, and to assess the role of selection and adaptation in the evolution and phylogeography of this species. The adaptive evolution in mtDNA probably was not the main driving force during the postlacial colonization of Europe. However, signatures of adaptive evolution have been found - an amino acid change with possible functional consequences in one gene and an excess of radical changes in physical- chemical properties of amino acids in populations at the latitudinal (northern and southern) extremes of the bank vole distribution. Key...

Les systèmes toxine-antitoxine (TA) sont de petits modules génétiques largement présents dans les génomes bactériens. Ils codent pour une petite protéine toxique et une antitoxine. Ils sont classés en six types en fonction de la nature et du mode d'action de l'antitoxine. Ce travail a porté sur l'étude du type I, pour lequel l'antitoxine est un ARN antisens qui cible l'ARNm de la toxine afin de réprimer son expression. Au cours de cette thèse, nous avons étudié le système aapA3/IsoA3, codé sur le chromosome du pathogène gastrique humain Helicobacter pylori. À ce jour, la plupart des systèmes TA ont été étudiés à l'aide de systèmes d'expression artificiels, qui ne permettent pas de caractériser la régulation transcriptionnelle ou post-transcriptionnelle. En utilisant la létalité induite par l’expression chromosomique de la toxine obtenue en absence d’antitoxine, nous avons développé une sélection génétique de mutants suppresseurs révélés par séquençage haut-débit. Cette approche, appelée FASTBAC-Seq, nous a permis de cartographier une myriade de déterminants de toxicité localisés dans les régions codantes et non codantes du gène de la toxine AapA3. En particulier, certaines de ces mutations ont révélé l'existence de tige-boucles ARN transitoires qui agissent de manière co-transcriptionnelle pour empêcher l'initiation de la traduction pendant la synthèse de l'ARNm codant pour la toxine. Ces structures ARN métastables fonctionnelles sont nécessaires pour découpler les processus de transcription et de traduction et permettent la présence de ces gènes toxiques sur le chromosome bactérien. Bien que les ARNm non traduits deviennent rapidement instables, nos travaux ont également révélé l'existence de deux tige-boucles protectrices situées aux deux extrémités de l'ARNm. Ces structures secondaires empêchent des activités exonucléolytiques agissant en 5' et 3'. Dans l’ensemble, notre travail met en évidence les conséquences de la forte pression de sélection pour limiter l'expression des toxines sous laquelle évoluent les systèmes TA. Cela nous a permis de mieux comprendre l’influence du repliement secondaire des ARNm, non seulement lors de la régulation posttranscriptionnelle, mais aussi co-transcriptionnelle de l’expression de cette famille particulière de gènes. Ces caractéristiques de régulation basées sur l'ARN peuvent être exploitées à l'avenir pour des applications biotechnologiques (p. ex., production accrue de protéines par stabilisation d'ARNm) ou biomédicales (p.ex., développement de stratégies antimicrobiennes alternatives pour l'activation de la synthèse de toxines). / Toxin-antitoxin (TA) systems are small genetic modules widely present in bacterial genomes. They usually code for a small toxic protein and its cognate antitoxin and can be classified into six types depending on the nature and mode of action of the antitoxin. This work focuses on the study of type I, for which the antitoxin is an antisense RNA that targets the toxin mRNA to inhibit its expression. We characterized the aapA3/IsoA3 system, encoded on the chromosome of the human gastric pathogen Helicobacter pylori. To date, most TAs have been studied using artificial expression systems, which do not allow the characterization of transcriptional or post-transcriptional regulation. Taking advantage of the lethality induced by the toxin chromosomal expression in the absence of antitoxin, we developed a high-throughput genetic selection of suppressor mutations revealed by Next-Generation Sequencing. This approach, named FASTBAC-Seq, allowed us to map a myriad of toxicity determinants located in both, coding and noncoding regions, of the aapA3 toxic gene. More precisely, some suppressor mutations revealed the existence of transient RNA hairpins that act co-transcriptionally to prevent translation initiation while the toxinencoding mRNA is being made. Such functional RNA metastable structures are essential to uncouple the transcription and translation processes and allow the presence of these toxic genes on bacterial chromosomes. Although untranslated mRNAs become rapidly unstable, our work also revealed the presence of two protective stem-loops located at both mRNA ends that prevent from both, 5’ and 3’ exonucleolytic activity. Altogether, our work evidenced the consequences of the strong selection pressure to silence toxin expression under which the TAs evolve, and highlighted the key role of mRNA folding in the co- and post-transcriptional regulation of this family of genes. These RNA-based regulatory mechanisms may be exploited in the future for biotechnological (e.g., increased protein production through mRNA stabilization) or biomedical (e.g., development of alternative antimicrobial strategies aiming at the activation of toxin synthesis) applications.