Efficient inference of recombination hot regions in bacterial genomes

Epistasis, core-genome disharmony, and adaptation in recombining bacteria

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

Homologous Recombination Shapes the Architecture and Evolution of Bacterial Genomes

Homologous recombination is a key evolutionary force that varies considerably across bacterial species. However, how the landscape of homologous recombination varies across genes and within individual genomes has only been studied in a few species. Here, we used Approximate Bayesian Computation to estimate the recombination rate along the genomes of 145 bacterial species. Our results show that homologous recombination varies greatly along bacterial genomes and shapes many aspects of genome architecture and evolution. The genomic landscape of recombination presents several key signatures: rates are highest near the origin of replication in most species, patterns of recombination generally appear symmetrical in both replichores (i.e. replicational halves of circular chromosomes) and most species have genomic hotpots of recombination. Furthermore, many closely related species share conserved landscapes of recombination across orthologs indicating that recombination landscapes are conserved over significant evolutionary distances. We show evidence that recombination drives the evolution of GC-content through increasing the effectiveness of selection and not through biased gene conversion, thereby contributing to an ongoing debate. Finally, we demonstrate that the rate of recombination varies across gene function and that many hotspots of recombination are associated with adaptive and mobile regions often encoding genes involved in pathogenicity.

RecGraph: recombination-aware alignment of sequences to variation graphs

MOTIVATION

Bacterial genomes present more variability than human genomes, which requires important adjustments in computational tools that are developed for human data. In particular, bacteria exhibit a mosaic structure due to homologous recombinations, but this fact is not sufficiently captured by standard read mappers that align against linear reference genomes. The recent introduction of pangenomics provides some insights in that context, as a pangenome graph can represent the variability within a species. However, the concept of sequence-to-graph alignment that captures the presence of recombinations has not been previously investigated.

RESULTS

In this paper, we present the extension of the notion of sequence-to-graph alignment to a variation graph that incorporates a recombination, so that the latter are explicitly represented and evaluated in an alignment. Moreover, we present a dynamic programming approach for the special case where there is at most a recombination-we implement this case as RecGraph. From a modelling point of view, a recombination corresponds to identifying a new path of the variation graph, where the new arc is composed of two halves, each extracted from an original path, possibly joined by a new arc. Our experiments show that RecGraph accurately aligns simulated recombinant bacterial sequences that have at most a recombination, providing evidence for the presence of recombination events.

AVAILABILITY AND IMPLEMENTATION

Our implementation is open source and available at https://github.com/AlgoLab/RecGraph.

Dynamics of bacterial recombination in the human gut microbiome

Horizontal gene transfer (HGT) is a ubiquitous force in microbial evolution. Previous work has shown that the human gut is a hotspot for gene transfer between species, but the more subtle exchange of variation within species-also known as recombination-remains poorly characterized in this ecosystem. Here, we show that the genetic structure of the human gut microbiome provides an opportunity to measure recent recombination events from sequenced fecal samples, enabling quantitative comparisons across diverse commensal species that inhabit a common environment. By analyzing recent recombination events in the core genomes of 29 human gut bacteria, we observed widespread heterogeneities in the rates and lengths of transferred fragments, which are difficult to explain by existing models of ecological isolation or homology-dependent recombination rates. We also show that natural selection helps facilitate the spread of genetic variants across strain backgrounds, both within individual hosts and across the broader population. These results shed light on the dynamics of in situ recombination, which can strongly constrain the adaptability of gut microbial communities.

The Helicobacter pylori Genome Project: insights into H. pylori population structure from analysis of a worldwide collection of complete genomes

Helicobacter pylori, a dominant member of the gastric microbiota, shares co-evolutionary history with humans. This has led to the development of genetically distinct H. pylori subpopulations associated with the geographic origin of the host and with differential gastric disease risk. Here, we provide insights into H. pylori population structure as a part of the Helicobacter pylori Genome Project (HpGP), a multi-disciplinary initiative aimed at elucidating H. pylori pathogenesis and identifying new therapeutic targets. We collected 1011 well-characterized clinical strains from 50 countries and generated high-quality genome sequences. We analysed core genome diversity and population structure of the HpGP dataset and 255 worldwide reference genomes to outline the ancestral contribution to Eurasian, African, and American populations. We found evidence of substantial contribution of population hpNorthAsia and subpopulation hspUral in Northern European H. pylori. The genomes of H. pylori isolated from northern and southern Indigenous Americans differed in that bacteria isolated in northern Indigenous communities were more similar to North Asian H. pylori while the southern had higher relatedness to hpEastAsia. Notably, we also found a highly clonal yet geographically dispersed North American subpopulation, which is negative for the cag pathogenicity island, and present in 7% of sequenced US genomes. We expect the HpGP dataset and the corresponding strains to become a major asset for H. pylori genomics.

The Helicobacter pylori Genome Project: insights into H. pylori population structure from analysis of a worldwide collection of complete genomes

Helicobacter pylori, a dominant member of the gastric microbiota, shares co-evolutionary history with humans. This has led to the development of genetically distinct H. pylori subpopulations associated with the geographic origin of the host and with differential gastric disease risk. Here, we provide insights into H. pylori population structure as a part of the Helicobacter pylori Genome Project (HpGP), a multi-disciplinary initiative aimed at elucidating H. pylori pathogenesis and identifying new therapeutic targets. We collected 1011 well-characterized clinical strains from 50 countries and generated high-quality genome sequences. We analysed core genome diversity and population structure of the HpGP dataset and 255 worldwide reference genomes to outline the ancestral contribution to Eurasian, African, and American populations. We found evidence of substantial contribution of population hpNorthAsia and subpopulation hspUral in Northern European H. pylori. The genomes of H. pylori isolated from northern and southern Indigenous Americans differed in that bacteria isolated in northern Indigenous communities were more similar to North Asian H. pylori while the southern had higher relatedness to hpEastAsia. Notably, we also found a highly clonal yet geographically dispersed North American subpopulation, which is negative for the cag pathogenicity island, and present in 7% of sequenced US genomes. We expect the HpGP dataset and the corresponding strains to become a major asset for H. pylori genomics.

Systems-Based Approach for Optimization of Assembly-Free Bacterial MLST Mapping

Epidemiological surveillance of bacterial pathogens requires real-time data analysis with a fast turnaround, while aiming at generating two main outcomes: (1) species-level identification and (2) variant mapping at different levels of genotypic resolution for population-based tracking and surveillance, in addition to predicting traits such as antimicrobial resistance (AMR). Multi-locus sequence typing (MLST) aids this process by identifying sequence types (ST) based on seven ubiquitous genome-scattered loci. In this paper, we selected one assembly-dependent and one assembly-free method for ST mapping and applied them with the default settings and ST schemes they are distributed with, and systematically assessed their accuracy and scalability across a wide array of phylogenetically divergent Public Health-relevant bacterial pathogens with available MLST databases. Our data show that the optimal k-mer length for stringMLST is species-specific and that genome-intrinsic and -extrinsic features can affect the performance and accuracy of the program. Although suitable parameters could be identified for most organisms, there were instances where this program may not be directly deployable in its current format. Next, we integrated stringMLST into our freely available and scalable hierarchical-based population genomics platform, ProkEvo, and further demonstrated how the implementation facilitates automated, reproducible bacterial population analysis.

Zoonotic Transmission of Campylobacter jejuni to Caretakers From Sick Pen Calves Carrying a Mixed Population of Strains With and Without Guillain Barré Syndrome-Associated Lipooligosaccharide Loci

Campylobacter jejuni causes foodborne gastroenteritis and may trigger acute autoimmune sequelae including Guillain Barré Syndrome. Onset of neuromuscular paralysis is associated with exposure to C. jejuni lipooligosaccharide (LOS) classes A, B, C, D, and E that mimic and evoke antibodies against gangliosides on myelin and axons of peripheral nerves. Family members managing a Michigan dairy operation reported recurring C. jejuni gastroenteritis. Because dairy cattle are known to shed C. jejuni, we hypothesized that calves in the sick pen were the source of human infections. Fecal samples obtained from twenty-five calves, one dog, and one asymptomatic family member were cultured for Campylobacter. C. jejuni isolates were obtained from thirteen calves and the family member: C. coli from two calves, and C. hyointestinalis from two calves. Some calves had diarrhea; most were clinically normal. Typing of lipooligosaccharide biosynthetic loci showed that eight calf C. jejuni isolates fell into classes A, B, and C. Two calf isolates and the human isolate possessed LOS class E, associated mainly with enteric disease and rarely with Guillain Barré Syndrome. Multi-locus sequence typing, porA and flaA typing, and whole genome comparisons of the thirteen C. jejuni isolates indicated that the three LOS class E strains that included the human isolate were closely related, indicating zoonotic transmission. Whole-genome comparisons revealed that isolates differed in virulence gene content, particularly in loci encoding biosynthesis of surface structures. Family members experienced diarrheal illness repeatedly over 2 years, yet none experienced GBS despite exposure to calves carrying invasive C. jejuni with LOS known to elicit antiganglioside autoantibodies.

Horizontal gene transfer and adaptive evolution in bacteria

Horizontal gene transfer (HGT) is arguably the most conspicuous feature of bacterial evolution. Evidence for HGT is found in most bacterial genomes. Although HGT can considerably alter bacterial genomes, not all transfer events may be biologically significant and may instead represent the outcome of an incessant evolutionary process that only occasionally has a beneficial purpose. When adaptive transfers occur, HGT and positive selection may result in specific, detectable signatures in genomes, such as gene-specific sweeps or increased transfer rates for genes that are ecologically relevant. In this Review, we first discuss the various mechanisms whereby HGT occurs, how the genetic signatures shape patterns of genomic variation and the distinct bioinformatic algorithms developed to detect these patterns. We then discuss the evolutionary theory behind HGT and positive selection in bacteria, and discuss the approaches developed over the past decade to detect transferred DNA that may be involved in adaptation to new environments.

Natural Horizontal Gene Transfer of Antimicrobial Resistance Genes in Campylobacter spp. From Turkeys and Swine

Antibiotic-resistant Campylobacter constitutes a serious threat to public health. The clonal expansion of resistant strains and/or the horizontal spread of resistance genes to other strains and species can hinder the clinical effectiveness of antibiotics to treat severe campylobacteriosis. Still, gaps exist in our understanding of the risks of acquisition and spread of antibiotic resistance in Campylobacter. While the in vitro transfer of antimicrobial resistance genes between Campylobacter species via natural transformation has been extensively demonstrated, experimental studies have favored the use of naked DNA to obtain transformants. In this study, we used experimental designs closer to real-world conditions to evaluate the possible transfer of antimicrobial resistance genes between Campylobacter strains of the same or different species (Campylobacter coli or Campylobacter jejuni) and originating from different animal hosts (swine or turkeys). This was evaluated in vitro through co-culture experiments and in vivo with dual-strain inoculation of turkeys, followed by whole genome sequencing of parental and newly emerged strains. In vitro, we observed four independent horizontal gene transfer events leading to the acquisition of resistance to beta-lactams (blaOXA), aminoglycosides [aph(2′′)-If and rpsL] and tetracycline [tet(O)]. Observed events involved the displacement of resistance-associated genes by a mutated version, or the acquisition of genomic islands harboring a resistance determinant by homologous recombination; we did not detect the transfer of resistance-carrying plasmids even though they were present in some strains. In vivo, we recovered a newly emerged strain with dual-resistance pattern and identified the replacement of an existing non-functional tet(O) by a functional tet(O) in the recipient strain. Whole genome comparisons allowed characterization of the events involved in the horizontal spread of resistance genes between Campylobacter following in vitro co-culture and in vivo dual inoculation. Our study also highlights the potential for antimicrobial resistance transfer across Campylobacter species originating from turkeys and swine, which may have implications for farms hosting both species in close proximity.

Quantifying bacterial evolution in the wild: A birthday problem for Campylobacter lineages

Measuring molecular evolution in bacteria typically requires estimation of the rate at which nucleotide changes accumulate in strains sampled at different times that share a common ancestor. This approach has been useful for dating ecological and evolutionary events that coincide with the emergence of important lineages, such as outbreak strains and obligate human pathogens. However, in multi-host (niche) transmission scenarios, where the pathogen is essentially an opportunistic environmental organism, sampling is often sporadic and rarely reflects the overall population, particularly when concentrated on clinical isolates. This means that approaches that assume recent common ancestry are not applicable. Here we present a new approach to estimate the molecular clock rate in Campylobacter that draws on the popular probability conundrum known as the 'birthday problem'. Using large genomic datasets and comparative genomic approaches, we use isolate pairs that share recent common ancestry to estimate the rate of nucleotide change for the population. Identifying synonymous and non-synonymous nucleotide changes, both within and outside of recombined regions of the genome, we quantify clock-like diversification to estimate synonymous rates of nucleotide change for the common pathogenic bacteria Campylobacter coli (2.4 x 10-6 s/s/y) and Campylobacter jejuni (3.4 x 10-6 s/s/y). Finally, using estimated total rates of nucleotide change, we infer the number of effective lineages within the sample time frame-analogous to a shared birthday-and assess the rate of turnover of lineages in our sample set over short evolutionary timescales. This provides a generalizable approach to calibrating rates in populations of environmental bacteria and shows that multiple lineages are maintained, implying that large-scale clonal sweeps may take hundreds of years or more in these species.

Of Traits and Trees: Probabilistic Distances under Continuous Trait Models for Dissecting the Interplay among Phylogeny, Model, and Data

Stochastic models of character trait evolution have become a cornerstone of evolutionary biology in an array of contexts. While probabilistic models have been used extensively for statistical inference, they have largely been ignored for the purpose of measuring distances between phylogeny-aware models. Recent contributions to the problem of phylogenetic distance computation have highlighted the importance of explicitly considering evolutionary model parameters and their impacts on molecular sequence data when quantifying dissimilarity between trees. By comparing two phylogenies in terms of their induced probability distributions that are functions of many model parameters, these distances can be more informative than traditional approaches that rely strictly on differences in topology or branch lengths alone. Currently, however, these approaches are designed for comparing models of nucleotide substitution and gene tree distributions, and thus, are unable to address other classes of traits and associated models that may be of interest to evolutionary biologists. Here, we expand the principles of probabilistic phylogenetic distances to compute tree distances under models of continuous trait evolution along a phylogeny. By explicitly considering both the degree of relatedness among species and the evolutionary processes that collectively give rise to character traits, these distances provide a foundation for comparing models and their predictions, and for quantifying the impacts of assuming one phylogenetic background over another while studying the evolution of a particular trait. We demonstrate the properties of these approaches using theory, simulations, and several empirical data sets that highlight potential uses of probabilistic distances in many scenarios. We also introduce an open-source R package named PRDATR for easy application by the scientific community for computing phylogenetic distances under models of character trait evolution.[Brownian motion; comparative methods; phylogeny; quantitative traits.].

Gene pool transmission of multidrug resistance among Campylobacter from livestock, sewage and human disease

The use of antimicrobials in human and veterinary medicine has coincided with a rise in antimicrobial resistance (AMR) in the food-borne pathogens Campylobacter jejuni and Campylobacter coli. Faecal contamination from the main reservoir hosts (livestock, especially poultry) is the principal route of human infection but little is known about the spread of AMR among source and sink populations. In particular, questions remain about how Campylobacter resistomes interact between species and hosts, and the potential role of sewage as a conduit for the spread of AMR. Here, we investigate the genomic variation associated with AMR in 168 C. jejuni and 92 C. coli strains isolated from humans, livestock and urban effluents in Spain. AMR was tested in vitro and isolate genomes were sequenced and screened for putative AMR genes and alleles. Genes associated with resistance to multiple drug classes were observed in both species and were commonly present in multidrug-resistant genomic islands (GIs), often located on plasmids or mobile elements. In many cases, these loci had alleles that were shared among C. jejuni and C. coli consistent with horizontal transfer. Our results suggest that specific antibiotic resistance genes have spread among Campylobacter isolated from humans, animals and the environment.

The recombination-cold region as an epidemiological marker of recombinogenic opportunistic pathogen Mycobacterium avium

Background

The rapid identification of lineage remains a challenge in the genotyping of clinical isolates of recombinogenic pathogens. The chromosome of Mycobacterium avium subsp. hominissuis (MAH), an agent of Mycobacterium avium complex (MAC) lung disease, is often mosaic and is composed of chromosomal segments originating from different lineages. This makes it difficult to infer the MAH lineage in a simple experimental set-up. To overcome this difficulty, we sought to identify chromosomal marker genes containing lineage-specific alleles by genome data mining.

Results

We conducted genetic population structure analysis, phylogenetic analysis, and a survey of historical recombination using data from 125 global MAH isolates. Six MAH lineages (EA1, EA2, SC1, SC2, SC3, and SC4) were identified in the current dataset. One P-450 gene (locus_tag MAH_0788/MAV_0940) in the recombination-cold region was found to have multiple alleles that could discriminate five lineages. By combining the information about allele type from one additional gene, the six MAH lineages as well as other M. avium subspecies were distinguishable. A recombination-cold region of 116 kb contains an insertion hotspot and is flanked by a mammalian cell-entry protein operon where allelic variants have previously been reported to occur. Hence, we speculate that the acquisition of lineage- or strain-specific insertions has introduced homology breaks in the chromosome, thereby reducing the chance of interlineage recombination.

Conclusions

The allele types of the newly identified marker genes can be used to predict major lineages of M. avium. The single nucleotide polymorphism typing approach targeting multiallelic loci in recombination-cold regions will facilitate the epidemiological study of MAC, and may also be useful for equivalent studies of other nontuberculous mycobacteria potentially carrying mosaic genomes.

Patterns of inter- and intrasubspecific homologous recombination inform eco-evolutionary dynamics of Xylella fastidiosa

High rates of homologous recombination (HR) in the bacterial plant pathogen Xylella fastidiosa have been previously detected. This study aimed to determine the extent and explore the ecological significance of HR in the genomes of recombinants experimentally generated by natural transformation and wild-type isolates. Both sets of strains displayed widespread HR and similar average size of recombined fragments consisting of random events (2-10 kb) of inter- and intrasubspecific recombination. A significantly higher proportion and greater lengths (>10 kb, maximum 31.5 kb) of recombined fragments were observed in subsp. morus and in strains isolated in Europe from intercepted coffee plants shipped from the Americas. Such highly recombinant strains pose a serious risk of emergence of novel variants, as genetically distinct and formerly geographically isolated genotypes are brought in close proximity by global trade. Recently recombined regions in wild-type strains included genes involved in regulation and signaling, host colonization, nutrient acquisition, and host evasion, all fundamental traits for X. fastidiosa ecology. Identification of four recombinant loci shared between wild-type and experimentally generated recombinants suggests potential hotspots of recombination in this naturally competent pathogen. These findings provide insights into evolutionary forces possibly affecting the adaptive potential to colonize the host environments of X. fastidiosa.

Transposon Mutagenesis in Chlamydia trachomatis Identifies CT339 as a ComEC Homolog Important for DNA Uptake and Lateral Gene Transfer

Transposon mutagenesis is a widely applied and powerful genetic tool for the discovery of genes associated with selected phenotypes. Chlamydia trachomatis is a clinically significant, obligate intracellular bacterium for which many conventional genetic tools and capabilities have been developed only recently. This report describes the successful development and application of a Himar transposon mutagenesis system for generating single-insertion mutant clones of C. trachomatis This system was used to generate a pool of 105 transposon mutant clones that included insertions in genes encoding flavin adenine dinucleotide (FAD)-dependent monooxygenase (C. trachomatis 148 [ct148]), deubiquitinase (ct868), and competence-associated (ct339) proteins. A subset of Tn mutant clones was evaluated for growth differences under cell culture conditions, revealing that most phenocopied the parental strain; however, some strains displayed subtle and yet significant differences in infectious progeny production and inclusion sizes. Bacterial burden studies in mice also supported the idea that a FAD-dependent monooxygenase (ct148) and a deubiquitinase (ct868) were important for these infections. The ct339 gene encodes a hypothetical protein with limited sequence similarity to the DNA-uptake protein ComEC. A transposon insertion in ct339 rendered the mutant incapable of DNA acquisition during recombination experiments. This observation, along with in situ structural analysis, supports the idea that this protein is playing a role in the fundamental process of lateral gene transfer similar to that of ComEC. In all, the development of the Himar transposon system for Chlamydia provides an effective genetic tool for further discovery of genes that are important for basic biology and pathogenesis aspects.IMPORTANCE Chlamydia trachomatis infections have an immense impact on public health; however, understanding the basic biology and pathogenesis of this organism has been stalled by the limited repertoire of genetic tools. This report describes the successful adaptation of an important tool that has been lacking in Chlamydia studies: transposon mutagenesis. This advance enabled the generation of 105 insertional mutants, demonstrating that numerous gene products are not essential for in vitro growth. Mammalian infections using these mutants revealed that several gene products are important for infections in vivo Moreover, this tool enabled the investigation and discovery of a gene critical for lateral gene transfer; a process fundamental to the evolution of bacteria and likely for Chlamydia as well. The development of transposon mutagenesis for Chlamydia has broad impact for the field and for the discovery of genes associated with selected phenotypes, providing an additional avenue for the discovery of molecular mechanisms used for pathogenesis and for a more thorough understanding of this important pathogen.

Analysis of genetic recombination and the pan-genome of a highly recombinogenic bacteriophage species

Bacteriophages are the most prevalent biological entities impacting on the ecosystem and are characterized by their extensive diversity. However, there are two aspects of phages that have remained largely unexplored: genetic flux by recombination between phage populations and characterization of specific phages in terms of the pan-genome. Here, we examined the recombination and pan-genome in Helicobacter pylori prophages at both the genome and gene level. In the genome-level analysis, we applied, for the first time, chromosome painting and fineSTRUCTURE algorithms to a phage species, and showed novel trends in inter-population genetic flux. Notably, hpEastAsia is a phage population that imported a higher proportion of DNA fragments from other phages, whereas the hpSWEurope phages showed weaker signatures of inter-population recombination, suggesting genetic isolation. The gene-level analysis showed that, after parameter tuning of the prokaryote pan-genome analysis program, H. pylori phages have a pan-genome consisting of 75 genes and a soft-core genome of 10 genes, which includes genes involved in the lytic and lysogenic life cycles. Quantitative analysis of recombination events of the soft-core genes showed no substantial variation in the intensity of recombination across the genes, but rather equally frequent recombination among housekeeping genes that were previously reported to be less prone to recombination. The signature of frequent recombination appears to reflect the host–phage evolutionary arms race, either by contributing to escape from bacterial immunity or by protecting the host by producing defective phages.

DCA for genome-wide epistasis analysis: the statistical genetics perspective

Direct coupling analysis (DCA) is a now widely used method to leverage statistical information from many similar biological systems to draw meaningful conclusions on each system separately. DCA has been applied with great success to sequences of homologous proteins, and also more recently to whole-genome population-wide sequencing data. We here argue that the use of DCA on the genome scale is contingent on fundamental issues of population genetics. DCA can be expected to yield meaningful results when a population is in the quasi-linkage equilibrium (QLE) phase studied by Kimura and others, but not, for instance, in a phase of clonal competition. We discuss how the exponential (Potts model) distributions emerge in QLE, and compare couplings to correlations obtained in a study of about 3000 genomes of the human pathogen Streptococcus pneumoniae.

Prediction of Hot Regions in PPIs Based on Improved Local Community Structure Detecting

The hot regions in PPIs are some assembly regions which are composed of the tightly packed HotSpots. The discovery of hot regions helps to understand life activities and has very important value for biological applications. The identification of hot regions is the basis for protein design and cancer prevention. The existing algorithms of predicting hot regions often have some defects, such as low accuracy and unstability. This paper proposes a novel hot region prediction method based on diverse biological characteristics. First, feature evaluation is employed by using an impoved mRMR method. Then, SVM is adopted to create cassification model based on the features selected. In addition, a new clustering algorithm, namely LCSD (Local community structure detecting), is developed to detect and analyze the conformation of hot regions. In the clustering process, the link similarity of protein residues is introduced to handle the boundary nodes. This algorithm can effectively deal with the missing residue nodes and control the local community boundaries. The results indicate that the spatial structure of hot regions can be obtained more effectively, and that our method is more effective than previous methods for precise identification of hot regions.

Population genomics of bacterial host adaptation

Some bacteria can transfer to new host species, and this poses a risk to human health. Indeed, an estimated 60% of all human pathogens have originated from other animal species. Similarly, human-to-animal transitions are recognized as a major threat to sustainable livestock production, and emerging pathogens impose an increasing burden on crop yield and global food security. Recent advances in high-throughput sequencing technologies have enabled comparative genomic analyses of bacterial populations from multiple hosts. Such studies are providing new insights into the evolutionary processes that underpin the establishment of bacteria in new host niches. A better understanding of the genetic and mechanistic basis for bacterial host adaptation may reveal novel targets for controlling infection or inform the design of approaches to limit the emergence of new pathogens.

Mitonuclear genomics challenges the theory of clonality in Trypanosoma congolense: Reply to Tibayrenc and Ayala

We recently published the first genomic diversity study of Trypanosoma congolense, a major aetiological agent of Animal African Trypanosomiasis. We demonstrated striking levels of SNP and indel diversity in the Eastern province of Zambia as a consequence of hybridization between divergent trypanosome lineages. We concluded that these and earlier findings in T. congolense challenge the predominant clonal evolution (PCE) model. In a recent comment, Tibayrenc and Ayala claim that there are many features in T. congolense supporting their theory of clonality. While we can follow the reasoning of the authors, we also identify major limitations in their theory and interpretations that resulted in incorrect conclusions. First, we argue that each T. congolense subgroup should be analysed independently as they may represent different (sub)species rather than "near-clades". Second, the authors neglect major findings of two robust population genetic studies on Savannah T. congolense that provide clear evidence of frequent recombination. Third, we reveal additional events of introgressive hybridization in T. congolense by analysing the maxicircle coding region using next-generation sequencing analyses. At last, we pinpoint two important misinterpretations by the authors and show that there are no spatially and temporally widespread clones in T. congolense. We stand by our earlier conclusions that the clonal framework is unlikely to accurately model the population structure of T. congolense. Other theoretical frameworks such as Maynard Smith's epidemic model may better represent the complex ancestry seen in T. congolense, where clones delimited in space and time arise against a background of recombination.

Investigation of recombination-intense viral groups and their genes in the Earth's virome

Bacteriophages (phages), or bacterial viruses, are the most abundant and diverse biological entities that impact the global ecosystem. Recent advances in metagenomics have revealed their rampant abundance in the biosphere. A fundamental aspect of bacteriophages that remains unexplored in metagenomic data is the process of recombination as a driving force in evolution that occurs among different viruses within the same bacterial host. Here, we systematically examined signatures of recombination in every gene from 211 species-level viral groups in a recently obtained dataset of the Earth’s virome that contain corresponding information on the host bacterial species. Our study revealed that signatures of recombination are widespread (84%) among the diverse viral groups. We identified 25 recombination-intense viral groups, widely distributed across the viral taxonomy, and present in bacterial species living in the human oral cavity. We also revealed a significant inverse association between the recombination-intense viral groups and Type II restriction endonucleases, that could be effective in reducing recombination among phages in a cell. Furthermore, we identified recombination-intense genes that are significantly enriched for encoding phage morphogenesis proteins. Changes in the viral genomic sequence by recombination may be important to escape cleavage by the host bacterial immune systems.

Population Structure and Local Adaptation of MAC Lung Disease Agent Mycobacterium avium subsp. hominissuis

Mycobacterium avium subsp. hominissuis (MAH) is one of the most common nontuberculous mycobacterial species responsible for chronic lung disease in humans. Despite increasing worldwide incidence, little is known about the genetic mechanisms behind the population evolution of MAH. To elucidate the local adaptation mechanisms of MAH, we assessed genetic population structure, the mutual homologous recombination, and gene content for 36 global MAH isolates, including 12 Japanese isolates sequenced in the present study. We identified five major MAH lineages and found that extensive mutual homologous recombination occurs among them. Two lineages (MahEastAsia1 and MahEastAsia2) were predominant in the Japanese isolates. We identified alleles unique to these two East Asian lineages in the loci responsible for trehalose biosynthesis (treS and mak) and in one mammalian cell entry operon, which presumably originated from as yet undiscovered mycobacterial lineages. Several genes and alleles unique to East Asian strains were located in the fragments introduced via recombination between East Asian lineages, suggesting implication of recombination in local adaptation. These patterns of MAH genomes are consistent with the signature of distribution conjugative transfer, a mode of sexual reproduction reported for other mycobacterial species.

Coalescent framework for prokaryotes undergoing interspecific homologous recombination

Coalescent process for prokaryote species is theoretically considered. Prokaryotes undergo homologous recombination with individuals of the same species (intraspecific recombination) and with individuals of other species (interspecific recombination). This work particularly focuses on interspecific recombination because intraspecific recombination has been well incorporated in coalescent framework. We present a simulation framework for generating SNP (single-nucleotide polymorphism) patterns that allows external DNA integration into host genome from other species. Using this simulation tool, msPro, we observed that the joint processes of intra- and interspecific recombination generate complex SNP patterns. The direct effect of interspecific recombination includes increased polymorphism. Because interspecific recombination is very rare in nature, it generates regions with exceptionally high polymorphism. Following interspecific recombination, intraspecific recombination cuts the integrated external DNA into small fragments, generating a complex SNP pattern that appears as if external DNA was integrated multiple times. The insight gained from our work using the msPro simulator will be useful for understanding and evaluating the relative contributions of intra- and interspecific recombination events in generating complex SNP patters in prokaryotes.

A Robust Framework for Microbial Archaeology

Microbial archaeology is flourishing in the era of high-throughput sequencing, revealing the agents behind devastating historical plagues, identifying the cryptic movements of pathogens in prehistory, and reconstructing the ancestral microbiota of humans. Here, we introduce the fundamental concepts and theoretical framework of the discipline, then discuss applied methodologies for pathogen identification and microbiome characterization from archaeological samples. We give special attention to the process of identifying, validating, and authenticating ancient microbes using high-throughput DNA sequencing data. Finally, we outline standards and precautions to guide future research in the field.

Combining pseudo dinucleotide composition with the Z curve method to improve the accuracy of predicting DNA elements: a case study in recombination spots

Pseudo dinucleotide composition (PseDNC) and Z curve showed excellent performance in the classification issues of nucleotide sequences in bioinformatics. Inspired by the principle of Z curve theory, we improved PseDNC to give the phase-specific PseDNC (psPseDNC). In this study, we used the prediction of recombination spots as a case to illustrate the capability of psPseDNC and also PseDNC fused with Z curve theory based on a novel machine learning method named large margin distribution machine (LDM). We verified that combining the two widely used approaches could generate better performance compared to only using PseDNC with a support vector machine based (SVM-based) model. The best Mathew's correlation coefficient (MCC) achieved by our LDM-based model was 0.7037 through the rigorous jackknife test and improved by ∼6.6%, ∼3.2%, and ∼2.4% compared with three previous studies. Similarly, the accuracy was improved by 3.2% compared with our previous iRSpot-PseDNC web server through an independent data test. These results demonstrate that the joint use of PseDNC and Z curve enhances performance and can extract more information from a biological sequence. To facilitate research in this area, we constructed a user-friendly web server for predicting hot/cold spots, HcsPredictor, which can be freely accessed from . In summary, we provided a united algorithm by integrating Z curve with PseDNC. We hope this united algorithm could be extended to other classification issues in DNA elements.

A Locus Encoding Variable Defense Systems against Invading DNA Identified in Streptococcus suis

Streptococcus suis, an important zoonotic pathogen, is known to have an open pan-genome and to develop a competent state. In S. suis, limited genetic lineages are suggested to be associated with zoonosis. However, little is known about the evolution of diversified lineages and their respective phenotypic or ecological characteristics. In this study, we performed comparative genome analyses of S. suis, with a focus on the competence genes, mobile genetic elements, and genetic elements related to various defense systems against exogenous DNAs (defense elements) that are associated with gene gain/loss/exchange mediated by horizontal DNA movements and their restrictions. Our genome analyses revealed a conserved competence-inducing peptide type (pherotype) of the competence system and large-scale genome rearrangements in certain clusters based on the genome phylogeny of 58 S. suis strains. Moreover, the profiles of the defense elements were similar or identical to each other among the strains belonging to the same genomic clusters. Our findings suggest that these genetic characteristics of each cluster might exert specific effects on the phenotypic or ecological differences between the clusters. We also found certain loci that shift several types of defense elements in S. suis. Of note, one of these loci is a previously unrecognized variable region in bacteria, at which strains of distinct clusters code for different and various defense elements. This locus might represent a novel defense mechanism that has evolved through an arms race between bacteria and invading DNAs, mediated by mobile genetic elements and genetic competence.

Rapid evolution of distinct Helicobacter pylori subpopulations in the Americas

For the last 500 years, the Americas have been a melting pot both for genetically diverse humans and for the pathogenic and commensal organisms associated with them. One such organism is the stomach-dwelling bacterium Helicobacter pylori, which is highly prevalent in Latin America where it is a major current public health challenge because of its strong association with gastric cancer. By analyzing the genome sequence of H. pylori isolated in North, Central and South America, we found evidence for admixture between H. pylori of European and African origin throughout the Americas, without substantial input from pre-Columbian (hspAmerind) bacteria. In the US, strains of African and European origin have remained genetically distinct, while in Colombia and Nicaragua, bottlenecks and rampant genetic exchange amongst isolates have led to the formation of national gene pools. We found three outer membrane proteins with atypical levels of Asian ancestry in American strains, as well as alleles that were nearly fixed specifically in South American isolates, suggesting a role for the ethnic makeup of hosts in the colonization of incoming strains. Our results show that new H. pylori subpopulations can rapidly arise, spread and adapt during times of demographic flux, and suggest that differences in transmission ecology between high and low prevalence areas may substantially affect the composition of bacterial populations.

Helicobacter pylori is one of the best studied examples of an intimate association between bacteria and humans, due to its ability to colonize the stomach for decades and to transmit from generation to generation. A number of studies have sought to link diversity in H. pylori to human migrations but there are some discordant signals such as an “out of Africa” dispersal within the last few thousand years that has left a much stronger signal in bacterial genomes than in human ones. In order to understand how such discrepancies arise, we have investigated the evolution of H. pylori during the recent colonization of the Americas. We find that bacterial populations evolve quickly and can spread rapidly to people of different ethnicities. Distinct new bacterial subpopulations have formed in Colombia from a European source and in Nicaragua and the US from African sources. Genetic exchange between bacterial populations is rampant within Central and South America but is uncommon within North America, which may reflect differences in prevalence. Our results also suggest that adaptation of bacteria to particular human ethnic groups may be confined to a handful of genes involved in interaction with the immune system.

A new subclass of intrinsic aminoglycoside nucleotidyltransferases, ANT(3")-II, is horizontally transferred among Acinetobacter spp. by homologous recombination

The emergence and spread of antibiotic resistance among Acinetobacter spp. have been investigated extensively. Most studies focused on the multiple antibiotic resistance genes located on plasmids or genomic resistance islands. On the other hand, the mechanisms controlling intrinsic resistance are still not well understood. In this study, we identified the novel subclass of aminoglycoside nucleotidyltransferase ANT(3")-II in Acinetobacter spp., which comprised numerous variants distributed among three main clades. All members of this subclass can inactivate streptomycin and spectinomycin. The three ant(3")-II genes, encoding for the three ANT(3")-II clades, are widely distributed in the genus Acinetobacter and always located in the same conserved genomic region. According to their prevalence, these genes are intrinsic in Acinetobacter baumannii, Acinetobacter pittii, and Acinetobacter gyllenbergii. We also demonstrated that the ant(3")-II genes are located in a homologous recombination hotspot and were recurrently transferred among Acinetobacter species. In conclusion, our findings demonstrated a novel mechanism of natural resistance in Acinetobacter spp., identified a novel subclass of aminoglycoside nucleotidyltransferase and provided new insight into the evolutionary history of intrinsic resistance genes.

The level of interest in intrinsic resistance genes has increased recently, and one of reasons is that their mobilization could lead to emergence of resistant pathogens. Insertion sequences (ISs) or plasmids can capture intrinsic resistance genes and disseminate them in bacterial populations. In this study, we identified a novel subclass of aminoglycoside nucleotidyltransferases which are intrinsic in A. baumannii and other Acinetobacter species. The genes encoding the aminoglycoside nucleotidyltransferase were frequently horizontally transferred between different Acinetobacter species by homologous recombination. This work reports a novel mechanism of natural resistance in Acinetobacter and an overlooked pathway for the dissemination of resistance among species in this genus.

Reconstructing the Ancestral Relationships Between Bacterial Pathogen Genomes

Following recent developments in DNA sequencing technology, it is now possible to sequence hundreds of whole genomes from bacterial isolates at relatively low cost. Analyzing this growing wealth of genomic data in terms of ancestral relationships can reveal many interesting aspects of the evolution, ecology, and epidemiology of bacterial pathogens. However, reconstructing the ancestry of a sample of bacteria remains challenging, especially for the majority of species where recombination is frequent. Here, we review and describe the computational techniques currently available to infer ancestral relationships, including phylogenetic methods that either ignore or account for the effect of recombination, as well as model-based and model-free phylogeny-independent approaches.

gmos: Rapid Detection of Genome Mosaicism over Short Evolutionary Distances

Prokaryotic and viral genomes are often altered by recombination and horizontal gene transfer. The existing methods for detecting recombination are primarily aimed at viral genomes or sets of loci, since the expensive computation of underlying statistical models often hinders the comparison of complete prokaryotic genomes. As an alternative, alignment-free solutions are more efficient, but cannot map (align) a query to subject genomes. To address this problem, we have developed gmos (Genome MOsaic Structure), a new program that determines the mosaic structure of query genomes when compared to a set of closely related subject genomes. The program first computes local alignments between query and subject genomes and then reconstructs the query mosaic structure by choosing the best local alignment for each query region. To accomplish the analysis quickly, the program mostly relies on pairwise alignments and constructs multiple sequence alignments over short overlapping subject regions only when necessary. This fine-tuned implementation achieves an efficiency comparable to an alignment-free tool. The program performs well for simulated and real data sets of closely related genomes and can be used for fast recombination detection; for instance, when a new prokaryotic pathogen is discovered. As an example, gmos was used to detect genome mosaicism in a pathogenic Enterococcus faecium strain compared to seven closely related genomes. The analysis took less than two minutes on a single 2.1 GHz processor. The output is available in fasta format and can be visualized using an accessory program, gmosDraw (freely available with gmos).

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.

SNP synteny analysis of Staphylococcus aureus and Pseudomonas aeruginosa population genomics

Genomic sequence diversity of a bacterial species mainly results from the frequency distribution of single nucleotide polymorphisms (SNPs). Here we report on an SNP matrix-based binary algorithm to determine the intraclonal or interclonal genomic diversity by the number of shared sequential SNPs, the so-called SNP synteny or haplotype. All SNP positions and the frequency and length distribution of haplotypes are determined from pairwise alignment of completely sequenced genomes. This metric is invariant regarding the reference genome chosen. Information is obtained about the size of haplotypes, genomic gradients of recombination frequency, relatedness of strains and population composition of a taxon or clonal populations. The approach is illustrated with whole genome data sets of Staphylococcus aureus and Pseudomonas aeruginosa strains.

Inferences from tip-calibrated phylogenies: a review and a practical guide

Molecular dating of phylogenetic trees is a growing discipline using sequence data to co‐estimate the timing of evolutionary events and rates of molecular evolution. All molecular‐dating methods require converting genetic divergence between sequences into absolute time. Historically, this could only be achieved by associating externally derived dates obtained from fossil or biogeographical evidence to internal nodes of the tree. In some cases, notably for fast‐evolving genomes such as viruses and some bacteria, the time span over which samples were collected may cover a significant proportion of the time since they last shared a common ancestor. This situation allows phylogenetic trees to be calibrated by associating sampling dates directly to the sequences representing the tips (terminal nodes) of the tree. The increasing availability of genomic data from ancient DNA extends the applicability of such tip‐based calibration to a variety of taxa including humans, extinct megafauna and various microorganisms which typically have a scarce fossil record. The development of statistical models accounting for heterogeneity in different aspects of the evolutionary process while accommodating very large data sets (e.g. whole genomes) has allowed using tip‐dating methods to reach inferences on divergence times, substitution rates, past demography or the age of specific mutations on a variety of spatiotemporal scales. In this review, we summarize the current state of the art of tip dating, discuss some recent applications, highlight common pitfalls and provide a ‘how to’ guide to thoroughly perform such analyses.

Comparative Genomic Analyses of the Moraxella catarrhalis Serosensitive and Seroresistant Lineages Demonstrate Their Independent Evolution

The bacterial speciesMoraxella catarrhalishas been hypothesized as being composed of two distinct lineages (referred to as the seroresistant [SR] and serosensitive [SS]) with separate evolutionary histories based on several molecular typing methods, whereas 16S ribotyping has suggested an additional split within the SS lineage. Previously, we characterized whole-genome sequences of 12 SR-lineage isolates, which revealed a relatively small supragenome when compared with other opportunistic nasopharyngeal pathogens, suggestive of a relatively short evolutionary history. Here, we performed whole-genome sequencing on 18 strains from both ribotypes of the SS lineage, an additional SR strain, as well as four previously identified highly divergent strains based on multilocus sequence typing analyses. All 35 strains were subjected to a battery of comparative genomic analyses which clearly show that there are three lineages-the SR, SS, and the divergent. The SR and SS lineages are closely related, but distinct from each other based on three different methods of comparison: Allelic differences observed among core genes; possession of lineage-specific sets of core and distributed genes; and by an alignment of concatenated core sequences irrespective of gene annotation. All these methods show that the SS lineage has much longer interstrain branches than the SR lineage indicating that this lineage has likely been evolving either longer or faster than the SR lineage. There is evidence of extensive horizontal gene transfer (HGT) within both of these lineages, and to a lesser degree between them. In particular, we identified very high rates of HGT between these two lineages for ß-lactamase genes. The four divergent strains aresui generis, being much more distantly related to both the SR and SS groups than these other two groups are to each other. Based on average nucleotide identities, gene content, GC content, and genome size, this group could be considered as a separate taxonomic group. The SR and SS lineages, although distinct, clearly form a single species based on multiple criteria including a large common core genome, average nucleotide identity values, GC content, and genome size. Although neither of these lineages arose from within the other based on phylogenetic analyses, the question of how and when these lineages split and then subsequently reunited in the human nasopharynx is explored.

Identifying lineage effects when controlling for population structure improves power in bacterial association studies

Bacteria pose unique challenges for genome-wide association studies because of strong structuring into distinct strains and substantial linkage disequilibrium across the genome(1,2). Although methods developed for human studies can correct for strain structure(3,4), this risks considerable loss-of-power because genetic differences between strains often contribute substantial phenotypic variability(5). Here, we propose a new method that captures lineage-level associations even when locus-specific associations cannot be fine-mapped. We demonstrate its ability to detect genes and genetic variants underlying resistance to 17 antimicrobials in 3,144 isolates from four taxonomically diverse clonal and recombining bacteria: Mycobacterium tuberculosis, Staphylococcus aureus, Escherichia coli and Klebsiella pneumoniae. Strong selection, recombination and penetrance confer high power to recover known antimicrobial resistance mechanisms and reveal a candidate association between the outer membrane porin nmpC and cefazolin resistance in E. coli. Hence, our method pinpoints locus-specific effects where possible and boosts power by detecting lineage-level differences when fine-mapping is intractable.

Genome-wide survey of codons under diversifying selection in a highly recombining bacterial species, Helicobacter pylori

Selection has been a central issue in biology in eukaryotes as well as prokaryotes. Inference of selection in recombining bacterial species, compared with clonal ones, has been a challenge. It is not known how codons under diversifying selection are distributed along the chromosome or among functional categories or how frequently such codons are subject to mutual homologous recombination. Here, we explored these questions by analysing genes present in >90% among 29 genomes of Helicobacter pylori, one of the bacterial species with the highest mutation and recombination rates. By a method for recombining sequences, we identified codons under diversifying selection (dN/dS> 1), which were widely distributed and accounted for ∼0.2% of all the codons of the genome. The codons were enriched in genes of host interaction/cell surface and genome maintenance (DNA replication,recombination, repair, and restriction modification system). The encoded amino acid residues were sometimes found adjacent to critical catalytic/binding residues in protein structures.Furthermore, by estimating the intensity of homologous recombination at a single nucleotide level, we found that these codons appear to be more frequently subject to recombination.We expect that the present study provides a new approach to population genomics of selection in recombining prokaryotes.

The Landscape of Realized Homologous Recombination in Pathogenic Bacteria

Recombination enhances the adaptive potential of organisms by allowing genetic variants to be tested on multiple genomic backgrounds. Its distribution in the genome can provide insight into the evolutionary forces that underlie traits, such as the emergence of pathogenicity. Here, we examined landscapes of realized homologous recombination of 500 genomes from ten bacterial species and found all species have "hot" regions with elevated rates relative to the genome average. We examined the size, gene content, and chromosomal features associated with these regions and the correlations between closely related species. The recombination landscape is variable and evolves rapidly. For example in Salmonella, only short regions of around 1 kb in length are hot whereas in the closely related species Escherichia coli, some hot regions exceed 100 kb, spanning many genes. Only Streptococcus pyogenes shows evidence for the positive correlation between GC content and recombination that has been reported for several eukaryotes. Genes with function related to the cell surface/membrane are often found in recombination hot regions but E. coli is the only species where genes annotated as "virulence associated" are consistently hotter. There is also evidence that some genes with "housekeeping" functions tend to be overrepresented in cold regions. For example, ribosomal proteins showed low recombination in all of the species. Among specific genes, transferrin-binding proteins are recombination hot in all three of the species in which they were found, and are subject to interspecies recombination.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

Coupling the CRISPR/Cas9 System with Lambda Red Recombineering Enables Simplified Chromosomal Gene Replacement in Escherichia coli

To date, most genetic engineering approaches coupling the type II Streptococcus pyogenes clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system to lambda Red recombineering have involved minor single nucleotide mutations. Here we show that procedures for carrying out more complex chromosomal gene replacements in Escherichia coli can be substantially enhanced through implementation of CRISPR/Cas9 genome editing. We developed a three-plasmid approach that allows not only highly efficient recombination of short single-stranded oligonucleotides but also replacement of multigene chromosomal stretches of DNA with large PCR products. By systematically challenging the proposed system with respect to the magnitude of chromosomal deletion and size of DNA insertion, we demonstrated DNA deletions of up to 19.4 kb, encompassing 19 nonessential chromosomal genes, and insertion of up to 3 kb of heterologous DNA with recombination efficiencies permitting mutant detection by colony PCR screening. Since CRISPR/Cas9-coupled recombineering does not rely on the use of chromosome-encoded antibiotic resistance, or flippase recombination for antibiotic marker recycling, our approach is simpler, less labor-intensive, and allows efficient production of gene replacement mutants that are both markerless and "scar"-less.

Tackling Drug Resistant Infection Outbreaks of Global Pandemic Escherichia coli ST131 Using Evolutionary and Epidemiological Genomics

High-throughput molecular screening is required to investigate the origin and diffusion of antimicrobial resistance in pathogen outbreaks. The most frequent cause of human infection is Escherichia coli, which is dominated by sequence type 131 (ST131)-a set of rapidly radiating pandemic clones. The highly infectious clades of ST131 originated firstly by a mutation enhancing conjugation and adhesion. Secondly, single-nucleotide polymorphisms occurred enabling fluoroquinolone-resistance, which is near-fixed in all ST131. Thirdly, broader resistance through beta-lactamases has been gained and lost frequently, symptomatic of conflicting environmental selective effects. This flexible approach to gene exchange is worrying and supports the proposition that ST131 will develop an even wider range of plasmid and chromosomal elements promoting antimicrobial resistance. To stop ST131, deep genome sequencing is required to understand the origin, evolution and spread of antimicrobial resistance genes. Phylogenetic methods that decipher past events can predict future patterns of virulence and transmission based on genetic signatures of adaptation and gene exchange. Both the effect of partial antimicrobial exposure and cell dormancy caused by variation in gene expression may accelerate the development of resistance. High-throughput sequencing can decode measurable evolution of cell populations within patients associated with systems-wide changes in gene expression during treatments. A multi-faceted approach can enhance assessment of antimicrobial resistance in E. coli ST131 by examining transmission dynamics between hosts to achieve a goal of pre-empting resistance before it emerges by optimising antimicrobial treatment protocols.

Measurably evolving pathogens in the genomic era

Current sequencing technologies have created unprecedented opportunities for studying microbial populations. For pathogens with comparatively low per-site mutation rates, such as DNA viruses and bacteria, whole-genome sequencing can reveal the accumulation of novel genetic variation between population samples taken at different times. The concept of 'measurably evolving populations' and related analytical approaches have provided powerful insights for fast-evolving RNA viruses, but their application to other pathogens is still in its infancy. We argue that previous distinctions between slow- and fast-evolving pathogens become blurred once evolution is assessed at a genome-wide scale, and we highlight important analytical challenges to be overcome to infer pathogen population dynamics from genomic data.

Characterizing the genetic basis of bacterial phenotypes using genome-wide association studies: a new direction for bacteriology

Genome-wide association studies (GWASs) have become an increasingly important approach for eukaryotic geneticists, facilitating the identification of hundreds of genetic polymorphisms that are responsible for inherited diseases. Despite the relative simplicity of bacterial genomes, the application of GWASs to identify polymorphisms responsible for important bacterial phenotypes has only recently been made possible through advances in genome sequencing technologies. Bacterial GWASs are now about to come of age thanks to the availability of massive datasets, and because of the potential to bridge genomics and traditional genetic approaches that is provided by improving validation strategies. A small number of pioneering GWASs in bacteria have been published in the past 2 years, examining from 75 to more than 3,000 strains. The experimental designs have been diverse, taking advantage of different processes in bacteria for generating variation. Analysis of data from bacterial GWASs can, to some extent, be performed using software developed for eukaryotic systems, but there are important differences in genome evolution that must be considered. The greatest experimental advantage of bacterial GWASs is the potential to perform downstream validation of causality and dissection of mechanism. We review the recent advances and remaining challenges in this field and propose strategies to improve the validation of bacterial GWASs.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

A phylogeny-based sampling strategy and power calculator informs genome-wide associations study design for microbial pathogens

Whole genome sequencing is increasingly used to study phenotypic variation among infectious pathogens and to evaluate their relative transmissibility, virulence, and immunogenicity. To date, relatively little has been published on how and how many pathogen strains should be selected for studies associating phenotype and genotype. There are specific challenges when identifying genetic associations in bacteria which often comprise highly structured populations. Here we consider general methodological questions related to sampling and analysis focusing on clonal to moderately recombining pathogens. We propose that a matched sampling scheme constitutes an efficient study design, and provide a power calculator based on phylogenetic convergence. We demonstrate this approach by applying it to genomic datasets for two microbial pathogens: Mycobacterium tuberculosis and Campylobacter species.

Mobile elements drive recombination hotspots in the core genome of Staphylococcus aureus

Horizontal gene transfer is an important driver of bacterial evolution, but genetic exchange in the core genome of clonal species, including the major pathogen Staphylococcus aureus, is incompletely understood. Here we reveal widespread homologous recombination in S. aureus at the species level, in contrast to its near-complete absence between closely related strains. We discover a patchwork of hotspots and coldspots at fine scales falling against a backdrop of broad-scale trends in rate variation. Over megabases, homoplasy rates fluctuate 1.9-fold, peaking towards the origin-of-replication. Over kilobases, we find core recombination hotspots of up to 2.5-fold enrichment situated near fault lines in the genome associated with mobile elements. The strongest hotspots include regions flanking conjugative transposon ICE6013, the staphylococcal cassette chromosome (SCC) and genomic island νSaα. Mobile element-driven core genome transfer represents an opportunity for adaptation and challenges our understanding of the recombination landscape in predominantly clonal pathogens, with important implications for genotype–phenotype mapping.

Horizontal gene transfer occurs in most bacteria, yet it is unclear whether it happens in clonal species. Here, Everitt et al. show widespread within-species recombination, driven by mobile elements, in the genome of the pathogen Staphylococcus aureus, but no recombination between closely related strains.