Leonardos Mageiros

A reference pan-genome approach to comparative bacterial genomics: identification of novel epidemiological markers in pathogenic Campylobacter.

The increasing availability of hundreds of whole bacterial genomes provides opportunities for enhanced understanding of the genes and alleles responsible for clinically important phenotypes and how they evolved. However, it is a significant challenge to develop easy-to-use and scalable methods for characterizing these large and complex data and relating it to disease epidemiology. Existing approaches typically focus on either homologous sequence variation in genes that are shared by all isolates, or non-homologous sequence variation - focusing on genes that are differentially present in the population. Here we present a comparative genomics approach that simultaneously approximates core and accessory genome variation in pathogen populations and apply it to pathogenic species in the genus Campylobacter. A total of 7 published Campylobacter jejuni and Campylobacter coli genomes were selected to represent diversity across these species, and a list of all loci that were present at least once was compiled. After filtering duplicates a 7-isolate reference pan-genome, of 3,933 loci, was defined. A core genome of 1,035 genes was ubiquitous in the sample accounting for 59% of the genes in each isolate (average genome size of 1.68 Mb). The accessory genome contained 2,792 genes. A Campylobacter population sample of 192 genomes was screened for the presence of reference pan-genome loci with gene presence defined as a BLAST match of ≥70% identity over ≥50% of the locus length - aligned using MUSCLE on a gene-by-gene basis. A total of 21 genes were present only in C. coli and 27 only in C. jejuni, providing information about functional differences associated with species and novel epidemiological markers for population genomic analyses. Homologs of these genes were found in several of the genomes used to define the pan-genome and, therefore, would not have been identified using a single reference strain approach.

Ecological Overlap and Horizontal Gene Transfer in Staphylococcus aureus and Staphylococcus epidermidis

The opportunistic pathogens Staphylococcus aureus and Staphylococcus epidermidis represent major causes of severe nosocomial infection, and are associated with high levels of mortality and morbidity worldwide. These species are both common commensals on the human skin and in the nasal pharynx, but are genetically distinct, differing at 24% average nucleotide divergence in 1,478 core genes. To better understand the genome dynamics of these ecologically similar staphylococcal species, we carried out a comparative analysis of 324 S. aureus and S. epidermidis genomes, including 83 novel S. epidermidis sequences. A reference pan-genome approach and whole genome multilocus-sequence typing revealed that around half of the genome was shared between the species. Based on a BratNextGen analysis, homologous recombination was found to have impacted on 40% of the core genes in S. epidermidis, but on only 24% of the core genes in S. aureus. Homologous recombination between the species is rare, with a maximum of nine gene alleles shared between any two S. epidermidis and S. aureus isolates. In contrast, there was considerable interspecies admixture of mobile elements, in particular genes associated with the SaPIn1 pathogenicity island, metal detoxification, and the methicillin-resistance island SCCmec. Our data and analysis provide a context for considering the nature of recombinational boundaries between S. aureus and S. epidermidis and, the selective forces that influence realized recombination between these species.

Enhanced biofilm formation and multi-host transmission evolve from divergent genetic backgrounds in Campylobacter jejuni

Summary Multicellular biofilms are an ancient bacterial adaptation that offers a protective environment for survival in hostile habitats. In microaerophilic organisms such as C ampylobacter, biofilms play a key role in transmission to humans as the bacteria are exposed to atmospheric oxygen concentrations when leaving the reservoir host gut. Genetic determinants of biofilm formation differ between species, but little is known about how strains of the same species achieve the biofilm phenotype with different genetic backgrounds. Our approach combines genome‐wide association studies with traditional microbiology techniques to investigate the genetic basis of biofilm formation in 102 C ampylobacter jejuni isolates. We quantified biofilm formation among the isolates and identified hotspots of genetic variation in homologous sequences that correspond to variation in biofilm phenotypes. Thirteen genes demonstrated a statistically robust association including those involved in adhesion, motility, glycosylation, capsule production and oxidative stress. The genes associated with biofilm formation were different in the host generalist ST‐21 and ST‐45 clonal complexes, which are frequently isolated from multiple host species and clinical samples. This suggests the evolution of enhanced biofilm from different genetic backgrounds and a possible role in colonization of multiple hosts and transmission to humans.

Genome-wide association of functional traits linked with Campylobacter jejuni survival from farm to fork

Campylobacter jejuni is a major cause of bacterial gastroenteritis worldwide, primarily associated with the consumption of contaminated poultry. C. jejuni lineages vary in host range and prevalence in human infection, suggesting differences in survival throughout the poultry processing chain. From 7343 MLST-characterised isolates, we sequenced 600 C. jejuni and C. coli isolates from various stages of poultry processing and clinical cases. A genome-wide association study (GWAS) in C. jejuni ST-21 and ST-45 complexes identified genetic elements over-represented in clinical isolates that increased in frequency throughout the poultry processing chain. Disease-associated SNPs were distinct in these complexes, sometimes organised in haplotype blocks. The function of genes containing associated elements was investigated, demonstrating roles for cj1377c in formate metabolism, nuoK in aerobic survival and oxidative respiration, and cj1368-70 in nucleotide salvage. This work demonstrates the utility of GWAS for investigating transmission in natural zoonotic pathogen populations and provides evidence that major C. jejuni lineages have distinct genotypes associated with survival, within the host specific niche, from farm to fork.

Population-genomic insights into emergence, crop adaptation and dissemination of Pseudomonas syringae pathogens.

Many bacterial pathogens are well characterized but, in some cases, little is known about the populations from which they emerged. This limits understanding of the molecular mechanisms underlying disease. The crop pathogen Pseudomonas syringae sensu lato has been widely isolated from the environment, including wild plants and components of the water cycle, and causes disease in several economically important crops. Here, we compared genome sequences of 45 P. syringae crop pathogen outbreak strains with 69 closely related environmental isolates. Phylogenetic reconstruction revealed that crop pathogens emerged many times independently from environmental populations. Unexpectedly, differences in gene content between environmental populations and outbreak strains were minimal with most virulence genes present in both. However, a genome-wide association study identified a small number of genes, including the type III effector genes hopQ1 and hopD1, to be associated with crop pathogens, but not with environmental populations, suggesting that this small group of genes may play an important role in crop disease emergence. Intriguingly, genome-wide analysis of homologous recombination revealed that the locus Psyr 0346, predicted to encode a protein that confers antibiotic resistance, has been frequently exchanged among lineages and thus may contribute to pathogen fitness. Finally, we found that isolates from diseased crops and from components of the water cycle, collected during the same crop disease epidemic, form a single population. This provides the strongest evidence yet that precipitation and irrigation water are an overlooked inoculum source for disease epidemics caused by P. syringae.

Genome-wide identification of host-segregating epidemiological markers for source attribution in Campylobacter jejuni

ABSTRACT Campylobacter is among the most common worldwide causes of bacterial gastroenteritis. This organism is part of the commensal microbiota of numerous host species, including livestock, and these animals constitute potential sources of human infection. Molecular typing approaches, especially multilocus sequence typing (MLST), have been used to attribute the source of human campylobacteriosis by quantifying the relative abundance of alleles at seven MLST loci among isolates from animal reservoirs and human infection, implicating chicken as a major infection source. The increasing availability of bacterial genomes provides data on allelic variation at loci across the genome, providing the potential to improve the discriminatory power of data for source attribution. Here we present a source attribution approach based on the identification of novel epidemiological markers among a reference pan-genome list of 1,810 genes identified by gene-by-gene comparison of 884 genomes of Campylobacter jejuni isolates from animal reservoirs, the environment, and clinical cases. Fifteen loci involved in metabolic activities, protein modification, signal transduction, and stress response or coding for hypothetical proteins were selected as host-segregating markers and used to attribute the source of 42 French and 281 United Kingdom clinical C. jejuni isolates. Consistent with previous studies of British campylobacteriosis, analyses performed using STRUCTURE software attributed 56.8% of British clinical cases to chicken, emphasizing the importance of this host reservoir as an infection source in the United Kingdom. However, among French clinical isolates, approximately equal proportions of isolates were attributed to chicken and ruminant reservoirs, suggesting possible differences in the relative importance of animal host reservoirs and indicating a benefit for further national-scale attribution modeling to account for differences in production, behavior, and food consumption. IMPORTANCE Accurately quantifying the relative contribution of different host reservoirs to human Campylobacter infection is an ongoing challenge. This study, based on the development of a novel source attribution approach, provides the first results of source attribution in Campylobacter jejuni in France. A systematic analysis using gene-by-gene comparison of 884 genomes of C. jejuni isolates, with a pan-genome list of genes, identified 15 novel epidemiological markers for source attribution. The different proportions of French and United Kingdom clinical isolates attributed to each host reservoir illustrate a potential role for local/national variations in C. jejuni transmission dynamics.

Biofilm Morphotypes and Population Structure among Staphylococcus epidermidis from Commensal and Clinical Samples

Bacterial species comprise related genotypes that can display divergent phenotypes with important clinical implications. Staphylococcus epidermidis is a common cause of nosocomial infections and, critical to its pathogenesis, is its ability to adhere and form biofilms on surfaces, thereby moderating the effect of the host’s immune response and antibiotics. Commensal S. epidermidis populations are thought to differ from those associated with disease in factors involved in adhesion and biofilm accumulation. We quantified the differences in biofilm formation in 98 S. epidermidis isolates from various sources, and investigated population structure based on ribosomal multilocus typing (rMLST) and the presence/absence of genes involved in adhesion and biofilm formation. All isolates were able to adhere and form biofilms in in vitro growth assays and confocal microscopy allowed classification into 5 biofilm morphotypes based on their thickness, biovolume and roughness. Phylogenetic reconstruction grouped isolates into three separate clades, with the isolates in the main disease associated clade displaying diversity in morphotype. Of the biofilm morphology characteristics, only biofilm thickness had a significant association with clade distribution. The distribution of some known adhesion-associated genes (aap and sesE) among isolates showed a significant association with the species clonal frame. These data challenge the assumption that biofilm-associated genes, such as those on the ica operon, are genetic markers for less invasive S. epidermidis isolates, and suggest that phenotypic characteristics, such as adhesion and biofilm formation, are not fixed by clonal descent but are influenced by the presence of various genes that are mobile among lineages.

Recombination-Mediated Host Adaptation by Avian Staphylococcus aureus

Staphylococcus aureus are globally disseminated among farmed chickens causing skeletal muscle infections, dermatitis, and septicaemia. The emergence of poultry-associated lineages has involved zoonotic transmission from humans to chickens but questions remain about the specific adaptations that promote proliferation of chicken pathogens. We characterized genetic variation in a population of genome-sequenced S. aureus isolates of poultry and human origin. Genealogical analysis identified a dominant poultry-associated sequence cluster within the CC5 clonal complex. Poultry and human CC5 isolates were significantly distinct from each other and more recombination events were detected in the poultry isolates. We identified 44 recombination events in 33 genes along the branch extending to the poultry-specific CC5 cluster, and 47 genes were found more often in CC5 poultry isolates compared with those from humans. Many of these gene sequences were common in chicken isolates from other clonal complexes suggesting horizontal gene transfer among poultry associated lineages. Consistent with functional predictions for putative poultry-associated genes, poultry isolates showed enhanced growth at 42 °C and greater erythrocyte lysis on chicken blood agar in comparison with human isolates. By combining phenotype information with evolutionary analyses of staphylococcal genomes, we provide evidence of adaptation, following a human-to-poultry host transition. This has important implications for the emergence and dissemination of new pathogenic clones associated with modern agriculture.

Comparative Genomics Study of Staphylococcus epidermidis Isolates from Orthopedic-Device-Related Infections Correlated with Patient Outcome

ABSTRACT Staphylococcus epidermidis has emerged as an important opportunistic pathogen causing orthopedic-device-related infections (ODRI). This study investigated the association of genome variation and phenotypic features of the infecting S. epidermidis isolate with the clinical outcome for the infected patient. S. epidermidis isolates were collected from 104 patients with ODRI. Their clinical outcomes were evaluated, after an average of 26 months, as either “cured” or “not cured.” The isolates were tested for antibiotic susceptibility and biofilm formation. Whole-genome sequencing was performed on all isolates, and genomic variation was related to features associated with “cured” and “not cured.” Strong biofilm formation and aminoglycoside resistance were associated with a “not-cured” outcome (P = 0.031 and P < 0.001, respectively). Based on gene-by-gene analysis, some accessory genes were more prevalent in isolates from the “not-cured” group. These included the biofilm-associated bhp gene, the antiseptic resistance qacA gene, the cassette chromosome recombinase-encoding genes ccrA and ccrB, and the IS256-like transposase gene. This study identifies biofilm formation and antibiotic resistance as associated with poor outcome in S. epidermidis ODRI. Whole-genome sequencing identified specific genes associated with a “not-cured” outcome that should be validated in future studies. (The study has been registered at ClinicalTrials.gov with identifier NCT02640937.)

Lineage-specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co‐evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large‐scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co‐evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.