Bacterial microevolution and the Pangenome

The comparison of multiple genome sequences sampled from a bacterial population reveals considerable diversity in both the core and the accessory parts of the pangenome. This diversity can be analysed in terms of microevolutionary events that took place since the genomes shared a common ancestor, especially deletion, duplication, and recombination. We review the basic modelling ingredients used implicitly or explicitly when performing such a pangenome analysis. In particular, we describe a basic neutral phylogenetic framework of bacterial pangenome microevolution, which is not incompatible with evaluating the role of natural selection. We survey the different ways in which pangenome data is summarised in order to be included in microevolutionary models, as well as the main methodological approaches that have been proposed to reconstruct pangenome microevolutionary history

Identification of hidden population structure in time-scaled phylogenies

Abstract Population structure influences genealogical patterns, however data pertaining to how populations are structured are often unavailable or not directly observable. Inference of population structure is highly important in molecular epidemiology where pathogen phylogenetics is increasingly used to infer transmission patterns and detect outbreaks. Discrepancies between observed and idealised genealogies, such as those generated by the coalescent process, can be quantified, and where significant differences occur, may reveal the action of natural selection, host population structure, or other demographic and epidemiological heterogeneities. We have developed a fast non-parametric statistical test for detection of cryptic population structure in time-scaled phylogenetic trees. The test is based on contrasting estimated phylogenies with the theoretically expected phylodynamic ordering of common ancestors in two clades within a coalescent framework. These statistical tests have also motivated the development of algorithms which can be used to quickly screen a phylogenetic tree for clades which are likely to share a distinct demographic or epidemiological history. Epidemiological applications include identification of outbreaks in vulnerable host populations or rapid expansion of genotypes with a fitness advantage. To demonstrate the utility of these methods for outbreak detection, we applied the new methods to large phylogenies reconstructed from thousands of HIV-1 partial pol sequences. This revealed the presence of clades which had grown rapidly in the recent past, and was significantly concentrated in young men, suggesting recent and rapid transmission in that group. Furthermore, to demonstrate the utility of these methods for the study of antimicrobial resistance, we applied the new methods to a large phylogeny reconstructed from whole genome Neisseria gonorrhoeae sequences. We find that population structure detected using these methods closely overlaps with the appearance and expansion of mutations conferring antimicrobial resistance

The population structure of Pseudomonas aeruginosa is characterized by genetic isolation of exoU+ and exoS+ lineages

The diversification of microbial populations may be driven by many factors including adaptation to distinct ecological niches and barriers to recombination. We examined the population structure of the bacterial pathogen Pseudomonas aeruginosa by analyzing whole-genome sequences of 739 isolates from diverse sources. We confirmed that the population structure of P. aeruginosa consists of two major groups (referred to as Groups A and B) and at least two minor groups (Groups C1 and C2). Evidence for frequent intra-group but limited inter-group recombination in the core genome was observed, consistent with sexual isolation of the groups. Likewise, accessory genome analysis demonstrated more gene flow within Groups A and B than between these groups, and a few accessory genomic elements were nearly specific to one or the other group. In particular, the exoS gene was highly over-represented in Group A compared to Group B isolates (99.4% vs. 1.1%) and the exoU gene was highly over-represented in Group B compared to Group A isolates (95.2% vs. 1.8%). The exoS and exoU genes encode effector proteins secreted by the P. aeruginosa type III secretion system. Together these results suggest that the major P. aeruginosa groups defined in part by the exoS and exoU genes are divergent from each other, and that these groups are genetically isolated and may be ecologically distinct. Although both groups were globally distributed and caused human infections, certain groups predominated in some clinical contexts

Sustainable management of scab control through the integration of apple resistant cultivars in a low-fungicide input system

Evaluation of the sustainability of disease control strategies through experimental field studies is poorly documented. Plant genetic resistance to pathogens offers an interesting alternative to the use of pesticides, but pathogen populations are able to adapt, thus frequently resulting in the breakdown of the resistance. Partial resistance is considered to provide more durable resistance than major genes. However, partial resistance does not confer complete protection and its efficiency can also decrease. Developing appropriate strategies which integrate resistant cultivars into crop systems is therefore needed to increase the efficiency and durability of the resistance, whatever the kind of resistance. The aim of this study was to evaluate the relevance of the association of control methods in terms of increasing the efficiency and durability of two kinds of resistances: (i) partial resistance in the apple cultivar Reine des Reinettes and (ii) major resistance (Rvi6) in the apple cultivar Ariane, when planted in a region where the climatic conditions are very favourable to the disease. It was found that the removal of leaf litter in autumn together with spraying of fungicides in the case of moderate or high risks of scab infection resulted in a sustainable control of scab on Reine des Reinettes over a five-year period and delayed the breakdown of the major resistance Rvi6 of Ariane by virulent isolates

Genomic signatures of pre-resistance in Mycobacterium tuberculosis

Recent advances in bacterial whole-genome sequencing have resulted in a comprehensive catalog of antibiotic resistance genomic signatures in Mycobacterium tuberculosis. With a view to pre-empt the emergence of resistance, we hypothesized that pre-existing polymorphisms in susceptible genotypes (pre-resistance mutations) could increase the risk of becoming resistant in the future. We sequenced whole genomes from 3135 isolates sampled over a 17-year period. After reconstructing ancestral genomes on time-calibrated phylogenetic trees, we developed and applied a genome-wide survival analysis to determine the hazard of resistance acquisition. We demonstrate that M. tuberculosis lineage 2 has a higher risk of acquiring resistance than lineage 4, and estimate a higher hazard of rifampicin resistance evolution following isoniazid mono-resistance. Furthermore, we describe loci and genomic polymorphisms associated with a higher risk of resistance acquisition. Identifying markers of future antibiotic resistance could enable targeted therapy to prevent resistance emergence in M. tuberculosis and other pathogens

Direct characterization of circulating DNA in blood plasma using μLAS technology

Circulating cell-free DNA (cfDNA) is a powerful cancer biomarker for establishing targeted therapies or monitoring patients' treatment. However, current cfDNA characterization is severely limited by its low concentration, requiring the extensive use of amplification techniques. Here we report that the μLAS technology allows us to quantitatively characterize the size distribution of purified cfDNA in a few minutes, even when its concentration is as low as 1 pg/μL. Moreover, we show that DNA profiles can be directly measured in blood plasma with a minimal conditioning process to speed up considerably speed up the cfDNA analytical chain

Rapid detection of mobilized colistin resistance using a nucleic acid based lab-on-a-chip diagnostic system

The increasing prevalence of antimicrobial resistance is a serious threat to global public health. One of the most concerning trends is the rapid spread of Carbapenemase-Producing Organisms (CPO), where colistin has become the last-resort antibiotic treatment. The emergence of colistin resistance, including the spread of mobilized colistin resistance (mcr) genes, raises the possibility of untreatable bacterial infections and motivates the development of improved diagnostics for the detection of colistin-resistant organisms. This work demonstrates a rapid response for detecting the most recently reported mcr gene, mcr−9, using a portable and affordable lab-on-a-chip (LoC) platform, offering a promising alternative to conventional laboratory-based instruments such as real-time PCR (qPCR). The platform combines semiconductor technology, for non-optical real-time DNA sensing, with a smartphone application for data acquisition, visualization and cloud connectivity. This technology is enabled by using loop-mediated isothermal amplification (LAMP) as the chemistry for targeted DNA detection, by virtue of its high sensitivity, specificity, yield, and manageable temperature requirements. Here, we have developed the first LAMP assay for mcr−9 - showing high sensitivity (down to 100 genomic copies/reaction) and high specificity (no cross-reactivity with other mcr variants). This assay is demonstrated through supporting a hospital investigation where we analyzed nucleic acids extracted from 128 carbapenemase-producing bacteria isolated from clinical and screening samples and found that 41 carried mcr−9 (validated using whole genome sequencing). Average positive detection times were 6.58 ± 0.42 min when performing the experiments on a conventional qPCR instrument (n = 41). For validating the translation of the LAMP assay onto a LoC platform, a subset of the samples were tested (n = 20), showing average detection times of 6.83 ± 0.92 min for positive isolates (n = 14). All experiments detected mcr−9 in under 10 min, and both platforms showed no statistically significant difference (p-value > 0.05). When sample preparation and throughput capabilities are integrated within this LoC platform, the adoption of this technology for the rapid detection and surveillance of antimicrobial resistance genes will decrease the turnaround time for DNA detection and resistotyping, improving diagnostic capabilities, patient outcomes, and the management of infectious diseases

The bounded coalescent model : conditioning a genealogy on a minimum root date

The coalescent model represents how individuals sampled from a population may have originated from a last common ancestor. The bounded coalescent model is obtained by conditioning the coalescent model such that the last common ancestor must have existed after a certain date. This conditioned model arises in a variety of applications, such as speciation, horizontal gene transfer or transmission analysis, and yet the bounded coalescent model has not been previously analysed in detail. Here we describe a new algorithm to simulate from this model directly, without resorting to rejection sampling. We show that this direct simulation algorithm is more computationally efficient than the rejection sampling approach. We also show how to calculate the probability of the last common ancestor occurring after a given date, which is required to compute the probability density of realisations under the bounded coalescent model. Our results are applicable in both the isochronous (when all samples have the same date) and heterochronous (where samples can have different dates) settings. We explore the effect of setting a bound on the date of the last common ancestor, and show that it affects a number of properties of the resulting phylogenies. All our methods are implemented in a new R package called BoundedCoalescent which is freely available online