Physicians frequently face challenges in predicting which bacterial subpopulations are likely to cause severe infections. A more accurate prediction of virulence would improve diagnostics and limit the extent of antibiotic resistance. Nowadays, bacterial pathogens can be typed with high accuracy with advanced genotyping technologies. However, effective translation of bacterial genotyping data into assessments of clinical risk remains largely unexplored. The discovery of unknown virulence genes is another key determinant of successful prediction of infectious disease outcomes. The trial-and-error method for virulence gene discovery is time-consuming and resource-intensive. Selecting candidate genes with higher precision can thus reduce the number of futile trials. Several in silico candidate gene prioritisation (CGP) methods have been proposed to aid the search for genes responsible for inherited diseases in human. It remains uninvestigated as to how the CGP concept can assist with virulence gene discovery in bacterial pathogens. The main contribution of this thesis is to demonstrate the value of translational bioinformatics methods to address challenges in virulence prediction and virulence gene discovery. This thesis studied an important perinatal bacterial pathogen, group B streptococcus (GBS), the leading cause of neonatal sepsis and meningitis in developed countries. While several antibiotic prophylactic programs have successfully reduced the number of early-onset neonatal diseases (infections that occur within 7 days of life), the prevalence of late-onset infections (infections that occur between 7??30 days of life) remained constant. In addition, the widespread use of intrapartum prophylactic antibiotics may introduce undue risk of penicillin allergy and may trigger the development of antibiotic-resistant microorganisms. To minimising such potential harm, a more targeted approach of antibiotic use is required. Distinguish virulent GBS strains from colonising counterparts thus lays the cornerstone of achieving the goal of tailored therapy. There are three aims of this thesis: 1. Prediction of virulence by analysis of bacterial genotype data: To identify markers that may be associated with GBS virulence, statistical analysis was performed on GBS genotype data consisting of 780 invasive and 132 colonising S. agalactiae isolates. From a panel of 18 molecular markers studied, only alp3 gene (which encodes a surface protein antigen commonly associated with serotype V) showed an increased association with invasive diseases (OR=2.93, p=0.0003, Fisher??s exact test). Molecular serotype II (OR=10.0, p=0.0007) was found to have a significant association with early-onset neonatal disease when compared with late-onset diseases. To investigate whether clinical outcomes can be predicted by the panel of genotype markers, logistic regression and machine learning algorithms were applied to distinguish invasive isolates from colonising isolates. Nevertheless, the predictive analysis only yielded weak predictive power (area under ROC curve, AUC: 0.56??0.71, stratified 10-fold cross-validation). It was concluded that a definitive predictive relationship between the molecular markers and clinical outcomes may be lacking, and more discriminative markers of GBS virulence are needed to be investigated. 2. Development of two computational CGP methods to assist with functional discovery of prokaryotic genes: Two in silico CGP methods were developed based on comparative genomics: statistical CGP exploits the differences in gene frequency against phenotypic groups, while inductive CGP applies supervised machine learning to identify genes with similar occurrence patterns across a range of bacterial genomes. Three rediscovery experiments were carried out to evaluate the CGP methods: a) Rediscovery of peptidoglycan genes was attempted with 417 published bacterial genome sequences. Both CGP methods achieved their best AUC >0.911 in Escherichia coli K-12 and >0.978 Streptococcus agalactiae 2603 (SA-2603) genomes, with an average improvement in precision of >3.2-fold and a maximum of >27-fold using statistical CGP. A median AUC of >0.95 could still be achieved with as few as 10 genome examples in each group in the rediscovery of the peptidoglycan metabolism genes. b) A maximum of 109-fold improvement in precision was achieved in the rediscovery of anaerobic fermentation genes. c) In the rediscovery experiment with genes of 31 metabolic pathways in SA-2603, 14 pathways achieved an AUC >0.9 and 28 pathways achieved AUC >0.8 with the best inductive CGP algorithms. The results from the rediscovery experiments demonstrated that the two CGP methods can assist with the study of functionally uncategorised genomic regions and the discovery of bacterial gene-function relationships. 3. Application of the CGP methods to discover GBS virulence genes: Both statistical and inductive CGP were applied to assist with the discovery of unknown GBS virulence factors. Among a list of hypothetical protein genes, several highly-ranked genes were plausibly involved in molecular mechanisms in GBS pathogenesis, including several genes encoding family 8 glycosyltransferase, family 1 and family 2 glycosyltransferase, multiple adhesins, streptococcal neuraminidase, staphylokinase, and other factors that may have roles in contributing to GBS virulence. Such genes may be candidates for further biological validation. In addition, the co-occurrence of these genes with currently known virulence factors suggested that the virulence mechanisms of GBS in causing perinatal diseases are multifactorial. The procedure demonstrated in this prioritisation task should assist with the discovery of virulence genes in other pathogenic bacteria.

Neisseria meningitidis (Nm) is one of the major causes of bacterial meningitis in the African meningitis belt (AMB). This organism is part of the genus Neisseria, which includes ten human restricted species, mostly harmless commensals of the nasopharynx; however, Nm is capable of causing invasive meningococcal disease. The transition from carriage to pathogenic state remains perplexing, and strict virulece factors have yet to be identified. It has been hypothesised that non-pathogenic Neisseria (NPN) carried asymptomatically in the oroopharynx could play a role in modulating carriage of Nm, and therefore, its likelihood of invasion. In chapter 3, the diversity of the genus was characterised within a collection of 46 034 nasopharyngeal samples obtained across the AMB: five different species were identified, with Nm and NPNs displaying inversely related risk factors fo carriage. Chapter 5 presents the whole genome sequence (WGS) analysis of 107 Neisseria isolates unclassified by other methods. This higher genetic resolution, complemented with the use of a novel speciation approach, revealed seven novel Neisseria species, mostly collected in African countries. The invasive potential may also be due to the presence of particular genetic factors in the meningococcal genome. Chapter 4 presents the WGS comparison of 23 carried and invasive serogroup A Nm collected in Chad during the 2011 meningitis epidemic. Isolates from both phenotypic groups were found to be part of the same bacterial populations; however, discrete clusters were identified, associated with distinct age groups. These results indicate that genomic analyses are essential to appropriately study Neisseria diversity, and that lower resolution methods have greatly underestimated the diversity of the genus in Africa. The identification of Nm clusters associated with certain niches and of the differences in carriage risk factors suggests that variation in the environment, including the presence of NPNs, may be key in modulating carriage of Nm.

Yes / The majority of Mycobacterium tuberculosis isolates resistant to isoniazid harbour a mutation in katG. Since these mutations cause a wide range of minimum inhibitory concentrations (MICs), largely below the serum level reached with higher dosing (15 mg/L upon 15–20 mg/kg), the drug might still remain partly active in presence of a katG mutation. We therefore investigated which genetic mutations predict the level of phenotypic isoniazid resistance in clinical M. tuberculosis isolates. To this end, the association between known and unknown isoniazid resistance-conferring mutations in whole genome sequences, and the isoniazid MICs of 176 isolates was examined. We found mostly moderate-level resistance characterized by a mode of 6.4 mg/L for the very common katG Ser315Thr mutation, and always very high MICs (≥19.2 mg/L) for the combination of katG Ser315Thr and inhA c-15t. Contrary to common belief, isolates harbouring inhA c-15t alone, partly also showed moderate-level resistance, particularly when combined with inhA Ser94Ala. No overt association between low-confidence or unknown mutations, except in katG, and isoniazid resistance (level) was found. Except for the rare katG deletion, line probe assay is thus not sufficiently accurate to predict the level of isoniazid resistance for a single mutation in katG or inhA. / European Research Council (Starting Grant INTERRUPTB 311725 to CM, LR and BdJ), The Damien Foundation

Chloromethane is the most abundant halocarbon in the environment, and responsible for substantial ozone destruction in the stratosphere. Sources and sinks of chloromethane are still poorly constrained. Although synthesized and used industrially, chloromethane is mainly produced naturally, with major emissions from vegetation and especially the phyllosphere, i.e. the aerial parts of plants. Some phyllosphere epiphytes are methylotrophic bacteria which can use single carbon compounds such as methanol and chloromethane as the sole source of carbon and energy for growth. Most chloromethane-degrading strains isolated so far utilize the cmu pathway for growth with chloromethane which was characterized by the team. The main objective of this work was to investigate whether epiphytes may act as filters for plant emissions of chloromethane, by using a laboratory bipartite system consisting of the model plant Arabidopsis thaliana, known to produce chloromethane mainly by way of the HOL1 gene, and the reference chloromethane-degrading bacterial strain Methylobacterium extorquens CM4, possessing the cmu pathway and of known genome sequence. Three A. thaliana Col-0 variants with different levels of expression of HOL1, i.e. the wild-type strain, its homozygous HOL1 knockout mutant hol1 and an HOL1-OX HOL1 overexpressor, were selected using PCR and qRT-PCR. Chloromethane-degrading strains were isolated from the A. thaliana phyllosphere, and shown to contain the cmu pathway. A plasmid-based bacterial bioreporter for chloromethane was constructed which exploits the promoter region of the conserved chloromethane dehalogenase gene cmuA of strain CM4. It yields rapid, highly sensitive, specific and methyl halide concentration-dependent fluorescence. Application of the bioreporter to the three A. thaliana variants differing in expression of HOL1 investigated in this work suggested that they indeed synthesize different levels of chloromethane. Analysis by qPCR and qRT-PCR of metagenomic DNA from the leaf surface of these variants showed that the relative proportion and expression of cmuA in this environment paralleled HOL1 gene expression. Taken together, the results obtained indicate that even minor amounts of chloromethane produced by A. thaliana in the face of large emissions of methanol may provide a selective advantage for chloromethane-degrading methylotrophic bacteria in the phyllosphere environment. This suggests that chloromethane-degrading epiphytes may indeed act as filters for emissions of chloromethane from plants. Further experiments are envisaged to further assess the adaptation mechanisms of chloromethane-degrading bacteria in the phyllosphere, building upon the comparative genomic analysis of chloromethane-degrading strains which was also performed in this work, and on the preliminary investigations using high-throughput sequencing that were initiated.

Le chlorométhane est le plus abondant des composés organo-halogénés dans l’atmosphère et il est impliqué dans la destruction de l’ozone dans la stratosphère. Les sources et les puits de chlorométhane restent mal évalués. Bien que synthétisé et utilisé de manière industrielle, il est principalement produit naturellement, avec comme sources majeures les émissions provenant des végétaux et plus particulièrement de la phyllosphère, qui correspond aux parties aériennes des plantes. Certaines bactéries épiphytes de la phyllosphère sont des méthylotrophes capables d’utiliser des composés organiques sans liaison carbone-carbone comme le méthanol et le chlorométhane comme unique source de carbone et d’énergie pour leur croissance. La plupart des bactéries chlorométhane-dégradantes isolées jusqu’à présent utilisent une voie métabolique pour leur croissance sur chlorométhane appelée voie cmu (pour chloromethane utilisation), caractérisée par l’équipe. L’objectif principal de cette thèse a été de déterminer si des bactéries de la phyllosphère peuvent jouer le rôle de filtre pour l’émission de chlorométhane par les plantes. Dans ce but, un modèle de laboratoire a été mis en place, constitué de la plante Arabidopsis thaliana connue pour produire du chlorométhane par une réaction impliquant le gène HOL1, et la bactérie Methylobacterium extorquens CM4, souche de référence pour l’étude du métabolisme de dégradation du chlorométhane, qui possède la voie cmu et dont le génome complet a été séquencé et analysé. Des variants d’A. thaliana avec différents niveaux d’expression du gène HOL1 (le type sauvage, le mutant homozygote « knock-out » hol1 et un variant HOL1-OX avec surexpression) ont été sélectionnés par PCR et qPCR. Des souches bactériennes chlorométhane-dégradantes ont été isolées à partir de la phyllosphère d’A. thaliana, dont il a été montré qu’elles possèdent la voie cmu. Un bio-rapporteur bactérien pour le chlorométhane a été construit à l’aide d’un plasmide exploitant la région promotrice du gène conservé de la déshalogénase (cmuA) de la souche M. extorquens CM4. Il présente une réponse fluorescente rapide, sensible, et spécifique aux méthyl-halogénés de manière concentration-dépendante. L’application du bio-rapporteur aux trois variants d’A. thaliana étudiés suggère des niveaux d’émissions de chlorométhane différents. L’analyse, par qPCR et qRT-PCR, de l’ADN métagénomique extrait de la surface des feuilles a montré une corrélation entre la proportion relative de bactéries portant le gène cmuA et l’exprimant dans cet environnement, et l’expression du gène HOL1. Ces résultats indiquent qu’une production de chlorométhane, même très modeste par rapport aux fortes émissions de méthanol par A. thaliana, confère un avantage sélectif pour les bactéries épiphytes chlorométhane-dégradantes. Ces dernières pourraient ainsi bien jouer un rôle de filtre pour les émissions de chlorométhane de la phyllosphère vers l’atmosphère. En perspective, de nouvelles expériences complémentaires, basées sur l’analyse par génomique comparative des souches chlorométhane-dégradantes également effectuée dans le cadre du projet et sur une analyse par séquençage à haut-débit initiée dans ce travail, sont proposées pour améliorer la compréhension des mécanismes d’adaptation des bactéries chlorométhane-dégradantes dans la phyllosphère. / Chloromethane is the most abundant halocarbon in the environment, and responsible for substantial ozone destruction in the stratosphere. Sources and sinks of chloromethane are still poorly constrained. Although synthesized and used industrially, chloromethane is mainly produced naturally, with major emissions from vegetation and especially the phyllosphere, i.e. the aerial parts of plants. Some phyllosphere epiphytes are methylotrophic bacteria which can use single carbon compounds such as methanol and chloromethane as the sole source of carbon and energy for growth. Most chloromethane-degrading strains isolated so far utilize the cmu pathway for growth with chloromethane which was characterized by the team. The main objective of this work was to investigate whether epiphytes may act as filters for plant emissions of chloromethane, by using a laboratory bipartite system consisting of the model plant Arabidopsis thaliana, known to produce chloromethane mainly by way of the HOL1 gene, and the reference chloromethane-degrading bacterial strain Methylobacterium extorquens CM4, possessing the cmu pathway and of known genome sequence. Three A. thaliana Col-0 variants with different levels of expression of HOL1, i.e. the wild-type strain, its homozygous HOL1 knockout mutant hol1 and an HOL1-OX HOL1 overexpressor, were selected using PCR and qRT-PCR. Chloromethane-degrading strains were isolated from the A. thaliana phyllosphere, and shown to contain the cmu pathway. A plasmid-based bacterial bioreporter for chloromethane was constructed which exploits the promoter region of the conserved chloromethane dehalogenase gene cmuA of strain CM4. It yields rapid, highly sensitive, specific and methyl halide concentration-dependent fluorescence. Application of the bioreporter to the three A. thaliana variants differing in expression of HOL1 investigated in this work suggested that they indeed synthesize different levels of chloromethane. Analysis by qPCR and qRT-PCR of metagenomic DNA from the leaf surface of these variants showed that the relative proportion and expression of cmuA in this environment paralleled HOL1 gene expression. Taken together, the results obtained indicate that even minor amounts of chloromethane produced by A. thaliana in the face of large emissions of methanol may provide a selective advantage for chloromethane-degrading methylotrophic bacteria in the phyllosphere environment. This suggests that chloromethane-degrading epiphytes may indeed act as filters for emissions of chloromethane from plants. Further experiments are envisaged to further assess the adaptation mechanisms of chloromethane-degrading bacteria in the phyllosphere, building upon the comparative genomic analysis of chloromethane-degrading strains which was also performed in this work, and on the preliminary investigations using high-throughput sequencing that were initiated.

Les Streptomyces sont des bactéries de la rhizosphère qui contribuent à la fertilité des sols (recyclage de la matière organique), et à la croissance et la santé des plantes. Elles possèdent parmi les plus grands génomes bactériens (12 Mb) et présentent une variabilité génétique importante. Cette variabilité connue au niveau interspécifique n’a jamais été abordée à l’échelle de la population, c’est-à-dire entre individus sympatriques appartenant à la même espèce (souches sœurs) au sein de la même niche écologique. L’objectif de ce travail est de rechercher cette diversité dans les populations de l’écosystème sol forestier, d’approcher sa dynamique et son rôle fonctionnel. Après séquençage et comparaison des génomes complets, nous avons observé une grande diversité génomique en termes de taille, de présence/absence d’éléments extrachromosomiques, mais également en terme de présence/absence de gènes le long du chromosome. Un grand nombre d’événements d’insertions et délétions (indels) comprenant de 1 à 241 gènes différencient les individus de la population. Au vu des liens phylogénétiques étroits entre les individus, l’ancêtre commun de la population est récent, aussi la diversité génomique résulterait d’un flux massif et rapide de gènes. La forte prévalence d’éléments conjugatifs intégrés dans la population suggère que la conjugaison est le moteur prépondérant de cette diversité génomique. La production différentielle de métabolites spécialisés (antibiotiques) a également été utilisée pour estimer l’impact de la diversité génétique sur le fonctionnement de la population. Nous avons pu montrer que cette production était liée à des gènes spécifiques de souches et qu’elle pouvait constituer un bien commun pour la population. Nous proposons que l’évolution rapide du génome participe au maintien des mécanismes de cohésion sociale chez ces bactéries du sol. / Streptomyces are rhizospheric bacteria that contribute to soil fertility (recycling of organic matter), plant growth and health. They have among the largest bacterial genomes (12 Mb) with a high genetic variability. The genome variability, observed at the interspecific level has never been addressed within a population, i.e. between sympatric individuals belonging to the same species (Conspecific strains) within the same ecological niche. The objective of this work was to investigate this diversity in the forest soil ecosystem, to estimate its dynamics and its potential functional roles. After sequencing and comparison of the complete genomes, we observed a wide genomic diversity in terms of size, presence/absence of extrachromosomal elements, but also in terms of presence/absence of genes along the chromosome. A large number of insertion and deletion events (indels) from 1 to 241 genes differentiate individuals in the population. Given the close phylogenetic relationship of these strains, the common ancestor of the population is recent, hence the genomic diversity would result from a massive and rapid gene flux. The high prevalence of integrative and conjugative elements in the population suggests that conjugation could act as a driving force of this diversity. Differential production of specialized metabolites (antibiotics) was also used to estimate the impact of genetic diversity on population’s ecology. We were able to show that this production was linked to strain specific genes and that it may constitute a « public good » for the population. We propose that the rapid evolution of the genome contributes to the maintenance of social cohesion mechanisms within these soil bacteria.

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.