Genome sequencing and analysis of the Tasmanian devil and its transmissible cancer.

The Tasmanian devil (Sarcophilus harrisii), the largest marsupial carnivore, is endangered due to a transmissible facial cancer spread by direct transfer of living cancer cells through biting. Here we describe the sequencing, assembly, and annotation of the Tasmanian devil genome and whole-genome sequences for two geographically distant subclones of the cancer. Genomic analysis suggests that the cancer first arose from a female Tasmanian devil and that the clone has subsequently genetically diverged during its spread across Tasmania. The devil cancer genome contains more than 17,000 somatic base substitution mutations and bears the imprint of a distinct mutational process. Genotyping of somatic mutations in 104 geographically and temporally distributed Tasmanian devil tumors reveals the pattern of evolution and spread of this parasitic clonal lineage, with evidence of a selective sweep in one geographical area and persistence of parallel lineages in other populations.

Predicting Long Pendant Edges in Model Phylogenies, with Applications to Biodiversity and Tree Inference.

In the simplest phylogenetic diversification model (the pure-birth Yule process), lineages split independently at a constant rate $\lambda$ for time $t$. The length of a randomly chosen edge (either interior or pendant) in the resulting tree has an expected value that rapidly converges to $\frac{1}{2\lambda}$ as $t$ grows and thus is essentially independent of $t$. However, the behavior of the length $L$ of the longest pendant edge reveals remarkably different behavior: $L$ converges to $t/2$ as the expected number of leaves grows. Extending this model to allow an extinction rate $\mu$ (where $\mu<\lambda$), we also establish a similar result for birth-death trees, except that $t/2$ is replaced by $t/2 \cdot (1-\mu/\lambda)$. This "complete" tree may contain subtrees that have died out before time $t$; for the "reduced tree" that just involves the leaves present at time $t$ and their direct ancestors, the longest pendant edge length $L$ again converges to $t/2$. Thus, there is likely to be at least one extant species whose associated pendant branch attaches to the tree approximately half-way back in time to the origin of the entire clade. We also briefly consider the length of the shortest edges. Our results are relevant to phylogenetic diversity indices in biodiversity conservation, and to quantifying the length of aligned sequences required to correctly infer a tree. We compare our theoretical results with simulations and with the branch lengths from a recent phylogenetic tree of all mammals. [Birth-death process; phylogenetic diversification models; phylogenetic diversity.].

RCandy: an R package for visualizing homologous recombinations in bacterial genomes.

SUMMARY: Homologous recombination is an important evolutionary process in bacteria and other prokaryotes, which increases genomic sequence diversity and can facilitate adaptation. Several methods and tools have been developed to detect genomic regions recently affected by recombination. Exploration and visualization of such recombination events can reveal valuable biological insights, but it remains challenging. Here, we present RCandy, a platform-independent R package for rapid, simple and flexible visualization of recombination events in bacterial genomes. AVAILABILITY AND IMPLEMENTATION: RCandy is an R package freely available for use under the MIT license. It is platform-independent and has been tested on Windows, Linux and MacOSX. The source code comes together with a detailed vignette available on GitHub at https://github.com/ChrispinChaguza/RCandy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Land-use Suitability is Not an Intrinsic Property of a Land Parcel.

Agricultural production has economic, environmental, social and cultural consequences beyond farm boundaries, but information about these impacts is not readily available to decision makers. This study applied the land use suitability concept by carrying out an assessment of a region that has the potential for intensification of agricultural production, but where eutrophication of river and estuary receiving environments due to nitrogen enrichment is a significant issue. The assessment evaluated three indicators for each farmable land parcel in the region: productive potential (the inherent productive and economic potential of the parcel), relative contribution (the potential for the parcel to contribute nitrogen to receiving environments compared to other land parcels), and pressure (the load of nitrogen delivered to receiving environments compared to the loads that ensure environmental objectives are achieved). The assessment indicated that land with high suitability for land-use intensification in Southland is limited because areas with high productive potential and low relative contribution rarely coincide with receiving environments with low pressure. Existing data, methods and models can be used to calculate the indicators under different choices for regional land-use intensity and receiving environment objectives. However, the spatial resolution and accuracy that is achievable may preclude using assessment outputs to make land use decisions at small spatial scales such as individual farms. The study highlighted that land use suitability is not an intrinsic property of a land parcel because it is dependent on choices about land use elsewhere in the landscape and the environmental objectives, and that land use suitability is inherently subjective because of decisions that concern how indicators are combined and weighted.

Fast and flexible bacterial genomic epidemiology with PopPUNK.

The routine use of genomics for disease surveillance provides the opportunity for high-resolution bacterial epidemiology. Current whole-genome clustering and multilocus typing approaches do not fully exploit core and accessory genomic variation, and they cannot both automatically identify, and subsequently expand, clusters of significantly similar isolates in large data sets spanning entire species. Here, we describe PopPUNK (Population Partitioning Using Nucleotide K -mers), a software implementing scalable and expandable annotation- and alignment-free methods for population analysis and clustering. Variable-length k-mer comparisons are used to distinguish isolates' divergence in shared sequence and gene content, which we demonstrate to be accurate over multiple orders of magnitude using data from both simulations and genomic collections representing 10 taxonomically widespread species. Connections between closely related isolates of the same strain are robustly identified, despite interspecies variation in the pairwise distance distributions that reflects species' diverse evolutionary patterns. PopPUNK can process 103-104 genomes in a single batch, with minimal memory use and runtimes up to 200-fold faster than existing model-based methods. Clusters of strains remain consistent as new batches of genomes are added, which is achieved without needing to reanalyze all genomes de novo. This facilitates real-time surveillance with consistent cluster naming between studies and allows for outbreak detection using hundreds of genomes in minutes. Interactive visualization and online publication is streamlined through the automatic output of results to multiple platforms. PopPUNK has been designed as a flexible platform that addresses important issues with currently used whole-genome clustering and typing methods, and has potential uses across bacterial genetics and public health research.

The MITRE trial protocol: a study to evaluate the microbiome as a biomarker of efficacy and toxicity in cancer patients receiving immune checkpoint inhibitor therapy.

BACKGROUND: The gut microbiome is implicated as a marker of response to immune checkpoint inhibitors (ICI) based on preclinical mouse models and preliminary observations in limited patient series. Furthermore, early studies suggest faecal microbial transfer may have therapeutic potential, converting ICI non-responders into responders. So far, identification of specific responsible bacterial taxa has been inconsistent, which limits future application. The MITRE study will explore and validate a microbiome signature in a larger scale prospective study across several different cancer types. METHODS: Melanoma, renal cancer and non-small cell lung cancer patients who are planned to receive standard immune checkpoint inhibitors are being recruited to the MITRE study. Longitudinal stool samples are collected prior to treatment, then at 6 weeks, 3, 6 and 12 months during treatment, or at disease progression/recurrence (whichever is sooner), as well as after a severe (≥grade 3 CTCAE v5.0) immune-related adverse event. Additionally, whole blood, plasma, buffy coat, RNA and peripheral blood mononuclear cells (PBMCs) is collected at similar time points and will be used for exploratory analyses. Archival tumour tissue, tumour biopsies at progression/relapse, as well as any biopsies from body organs collected after a severe toxicity are collected. The primary outcome measure is the ability of the microbiome signature to predict 1 year progression-free survival (PFS) in patients with advanced disease. Secondary outcomes include microbiome correlations with toxicity and other efficacy end-points. Biosamples will be used to explore immunological and genomic correlates. A sub-study will evaluate both COVID-19 antigen and antibody associations with the microbiome. DISCUSSION: There is an urgent need to identify biomarkers that are predictive of treatment response, resistance and toxicity to immunotherapy. The data generated from this study will both help inform patient selection for these drugs and provide information that may allow therapeutic manipulation of the microbiome to improve future patient outcomes. TRIAL REGISTRATION: NCT04107168 , ClinicalTrials.gov, registered 09/27/2019. Protocol V3.2 (16/04/2021).

The medial condylar wall is a reliable landmark to kinematically align the femoral component in medial UKA: an in-silico study.

PURPOSE: Kinematic alignment (KA) aligns the femoral implant perpendicular to the cylindrical axis in the frontal and axial plane. Identification of the kinematic axes when using the mini-invasive sub-quadricipital approach is challenging in unicompartmental knee arthroplasty (UKA). This study aims to assess if the orientation of condylar walls may be suitable for use as an anatomical landmark to kinematically align the femoral component in medial UKA. It was hypothesised that the medial wall of the medial condyle would prove to be a reliable anatomical landmark to set both the frontal and axial alignment of the femoral component in medial UKA. METHODS: 73 patients undergoing medial UKA had pre-operative CT imaging to generate 3D models. Those with osteophytes that impaired visualisation of the condylar walls were excluded. 28 patients were included in the study. The ideal KA was determined using the cylindrical axis in the frontal and axial plane. Simulations using the medial wall of the medial condyle (MWMC) and the lateral wall of the medial condyle (LWMC) were performed to set the frontal alignment. To set the axial alignment, the MWMC, LWMC, medial wall of the lateral condyle (MWLC), and medial diagonal line (MDL) anatomical landmarks were investigated. Differences between the ideal measured KA values and values obtained using landmarks were investigated. RESULTS: Use of the MWMC let to similar frontal alignment compared to the ideal KA (2.9° valgus vs 3.4° valgus, p = 0.371) with 46.4% (13/28) of measurements being [Formula: see text] 1.0° different from the ideal KA and only 1 simulation with greater than 4.0° difference. Use of the MWMC led to very similar axial alignments compared to the ideal KA (0.5° internal vs 0.0°, p = 0.960) with 75.0% (21/28) of measurements being [Formula: see text] 1.0o different from the ideal KA, and a maximum difference of 3.0°. Use of the MWLC and MDL was associated with significant statistical differences when compared to the ideal KA (p < 0.001 for both). CONCLUSIONS: The native orientation of the medial condylar wall seems to be a reliable anatomical landmark for aligning the femoral component in medial KA UKA in both the axial plane and frontal planes. Other assessed landmarks were shown to not be reliable. Clinical and radiographic assessments of the reliability of using the MWMC to set the frontal and axial orientation of the femoral component when performing a medial KA UKA are needed.

The Distribution and Spread of Susceptible and Resistant Neisseria gonorrhoeae Across Demographic Groups in a Major Metropolitan Center.

BACKGROUND: Genomic epidemiology studies of gonorrhea in the United States have primarily focused on national surveillance for antibiotic resistance, and patterns of local transmission between demographic groups of resistant and susceptible strains are unknown. METHODS: We analyzed a convenience sample of genome sequences, antibiotic susceptibility, and patient data from 897 gonococcal isolates cultured at the New York City (NYC) Public Health Laboratory from NYC Department of Health and Mental Hygiene (DOHMH) Sexual Health Clinic (SHC) patients, primarily in 2012-2013. We reconstructed the gonococcal phylogeny, defined transmission clusters using a 10 nonrecombinant single nucleotide polymorphism threshold, tested for clustering of demographic groups, and placed NYC isolates in a global phylogenetic context. RESULTS: The NYC gonococcal phylogeny reflected global diversity with isolates from 22/23 of the prevalent global lineages (96%). Isolates clustered on the phylogeny by patient sexual behavior (P < .001) and race/ethnicity (P < .001). Minimum inhibitory concentrations were higher across antibiotics in isolates from men who have sex with men compared to heterosexuals (P < .001) and white heterosexuals compared to black heterosexuals (P < .01). In our dataset, all large transmission clusters (≥10 samples) of N. gonorrhoeae were susceptible to ciprofloxacin, ceftriaxone, and azithromycin, and comprised isolates from patients across demographic groups. CONCLUSIONS: All large transmission clusters were susceptible to gonorrhea therapies, suggesting that resistance to empiric therapy was not a main driver of spread, even as risk for resistance varied across demographic groups. Further study of local transmission networks is needed to identify drivers of transmission.

Contribution of thermogenic mechanisms by male and female mice lacking pituitary adenylate cyclase-activating polypeptide in response to cold acclimation.

Pituitary adenylate cyclase-activating polypeptide (PACAP) is a neuropeptide critical to the regulation of the stress response, including having a role in energy homeostasis. Mice lacking PACAP are cold-sensitive and have impaired adrenergic-induced thermogenesis. Interestingly, Pacap null mice can survive cold housing if acclimated slowly, similar to observations in uncoupling protein 1 (UCP1)-deficient mice. We hypothesized that Pacap null mice use alternate thermogenic pathways to compensate for impaired adaptive thermogenesis when acclimated to cold. Observations of behavior and assessment of fiber type in skeletal muscles did not show evidence of prolonged burst shivering or changes in oxidative metabolism in male or female Pacap-/- mice during cold acclimation compared with Pacap+/+ mice. Despite previous work that has established impaired capacity for adaptive thermogenesis in Pacap null mice, adaptive thermogenesis can be induced in mice lacking PACAP to support survival with cold housing. Interestingly, sex-specific morphological and molecular differences in adipose tissue remodeling were observed in Pacap null mice compared with controls. Thus, sexual dimorphisms are highlighted in adipose tissue remodeling and thermogenesis with cold acclimation in the absence of PACAP.NEW & NOTEWORTHY This manuscript adds to the literature of endocrine regulation of adaptive thermogenesis and energy balance. It specifically describes the role of pituitary adenylate cyclase-activating polypeptide on the regulation of brown adipose tissue via the sympathetic nervous system with a focus on compensatory mechanisms of thermogenesis. We highlight sex-specific differences in energy metabolism.

Genome-wide epistasis and co-selection study using mutual information.

Covariance-based discovery of polymorphisms under co-selective pressure or epistasis has received considerable recent attention in population genomics. Both statistical modeling of the population level covariation of alleles across the chromosome and model-free testing of dependencies between pairs of polymorphisms have been shown to successfully uncover patterns of selection in bacterial populations. Here we introduce a model-free method, SpydrPick, whose computational efficiency enables analysis at the scale of pan-genomes of many bacteria. SpydrPick incorporates an efficient correction for population structure, which adjusts for the phylogenetic signal in the data without requiring an explicit phylogenetic tree. We also introduce a new type of visualization of the results similar to the Manhattan plots used in genome-wide association studies, which enables rapid exploration of the identified signals of co-evolution. Simulations demonstrate the usefulness of our method and give some insight to when this type of analysis is most likely to be successful. Application of the method to large population genomic datasets of two major human pathogens, Streptococcus pneumoniae and Neisseria meningitidis, revealed both previously identified and novel putative targets of co-selection related to virulence and antibiotic resistance, highlighting the potential of this approach to drive molecular discoveries, even in the absence of phenotypic data.

Genomic evolution of Neisseria gonorrhoeae since the preantibiotic era (1928-2013): antimicrobial use/misuse selects for resistance and drives evolution.

BACKGROUND: Multidrug-resistant Neisseria gonorrhoeae strains are prevalent, threatening gonorrhoea treatment globally, and understanding of emergence, evolution, and spread of antimicrobial resistance (AMR) in gonococci remains limited. We describe the genomic evolution of gonococci and their AMR, related to the introduction of antimicrobial therapies, examining isolates from 1928 (preantibiotic era) to 2013 in Denmark. This is, to our knowledge, the oldest gonococcal collection globally. METHODS: Lyophilised isolates were revived and examined using Etest (18 antimicrobials) and whole-genome sequencing (WGS). Quality-assured genome sequences were obtained for 191 viable and 40 non-viable isolates and analysed with multiple phylogenomic approaches. RESULTS: Gonococcal AMR, including an accumulation of multiple AMR determinants, started to emerge particularly in the 1950s-1970s. By the twenty-first century, resistance to most antimicrobials was common. Despite that some AMR determinants affect many physiological functions and fitness, AMR determinants were mainly selected by the use/misuse of gonorrhoea therapeutic antimicrobials. Most AMR developed in strains belonging to one multidrug-resistant (MDR) clade with close to three times higher genomic mutation rate. Modern N. gonorrhoeae was inferred to have emerged in the late-1500s and its genome became increasingly conserved over time. CONCLUSIONS: WGS of gonococci from 1928 to 2013 showed that no AMR determinants, except penB, were in detectable frequency before the introduction of gonorrhoea therapeutic antimicrobials. The modern gonococcus is substantially younger than previously hypothesized and has been evolving into a more clonal species, driven by the use/misuse of antimicrobials. The MDR gonococcal clade should be further investigated for early detection of strains with predispositions to develop and maintain MDR and for initiation of public health interventions.

Genomic and Phenotypic Variability in Neisseria gonorrhoeae Antimicrobial Susceptibility, England.

Antimicrobial resistance (AMR) in Neisseria gonorrhoeae is a global concern. Phylogenetic analyses resolve uncertainties regarding genetic relatedness of isolates with identical phenotypes and inform whether AMR is due to new mutations and clonal expansion or separate introductions by importation. We sequenced 1,277 isolates with associated epidemiologic and antimicrobial susceptibility data collected during 2013-2016 to investigate N. gonorrhoeae genomic variability in England. Comparing genetic markers and phenotypes for AMR, we identified 2 N. gonorrhoeae lineages with different antimicrobial susceptibility profiles and 3 clusters with elevated MICs for ceftriaxone, varying mutations in the penA allele, and different epidemiologic characteristics. Our results indicate N. gonorrhoeae with reduced antimicrobial susceptibility emerged independently and multiple times in different sexual networks in England, through new mutation or recombination events and by importation. Monitoring and control for AMR in N. gonorrhoeae should cover the entire population affected, rather than focusing on specific risk groups or locations.

Systematic longitudinal survey of invasive Escherichia coli in England demonstrates a stable population structure only transiently disturbed by the emergence of ST131.

Escherichia coli associated with urinary tract infections and bacteremia has been intensively investigated, including recent work focusing on the virulent, globally disseminated, multidrug-resistant lineage ST131. To contextualize ST131 within the broader E. coli population associated with disease, we used genomics to analyze a systematic 11-yr hospital-based survey of E. coli associated with bacteremia using isolates collected from across England by the British Society for Antimicrobial Chemotherapy and from the Cambridge University Hospitals NHS Foundation Trust. Population dynamics analysis of the most successful lineages identified the emergence of ST131 and ST69 and their establishment as two of the five most common lineages along with ST73, ST95, and ST12. The most frequently identified lineage was ST73. Compared to ST131, ST73 was susceptible to most antibiotics, indicating that multidrug resistance was not the dominant reason for prevalence of E. coli lineages in this population. Temporal phylogenetic analysis of the emergence of ST69 and ST131 identified differences in the dynamics of emergence and showed that expansion of ST131 in this population was not driven by sequential emergence of increasingly resistant subclades. We showed that over time, the E. coli population was only transiently disturbed by the introduction of new lineages before a new equilibrium was rapidly achieved. Together, these findings suggest that the frequency of E. coli lineages in invasive disease is driven by negative frequency-dependent selection occurring outside of the hospital, most probably in the commensal niche, and that drug resistance is not a primary determinant of success in this niche.

Comprehensive global genome dynamics of Chlamydia trachomatis show ancient diversification followed by contemporary mixing and recent lineage expansion.

Chlamydia trachomatis is the world's most prevalent bacterial sexually transmitted infection and leading infectious cause of blindness, yet it is one of the least understood human pathogens, in part due to the difficulties of in vitro culturing and the lack of available tools for genetic manipulation. Genome sequencing has reinvigorated this field, shedding light on the contemporary history of this pathogen. Here, we analyze 563 full genomes, 455 of which are novel, to show that the history of the species comprises two phases, and conclude that the currently circulating lineages are the result of evolution in different genomic ecotypes. Temporal analysis indicates these lineages have recently expanded in the space of thousands of years, rather than the millions of years as previously thought, a finding that dramatically changes our understanding of this pathogen's history. Finally, at a time when almost every pathogen is becoming increasingly resistant to antimicrobials, we show that there is no evidence of circulating genomic resistance in C. trachomatis.

Evaluation of parameters affecting performance and reliability of machine learning-based antibiotic susceptibility testing from whole genome sequencing data.

Prediction of antibiotic resistance phenotypes from whole genome sequencing data by machine learning methods has been proposed as a promising platform for the development of sequence-based diagnostics. However, there has been no systematic evaluation of factors that may influence performance of such models, how they might apply to and vary across clinical populations, and what the implications might be in the clinical setting. Here, we performed a meta-analysis of seven large Neisseria gonorrhoeae datasets, as well as Klebsiella pneumoniae and Acinetobacter baumannii datasets, with whole genome sequence data and antibiotic susceptibility phenotypes using set covering machine classification, random forest classification, and random forest regression models to predict resistance phenotypes from genotype. We demonstrate how model performance varies by drug, dataset, resistance metric, and species, reflecting the complexities of generating clinically relevant conclusions from machine learning-derived models. Our findings underscore the importance of incorporating relevant biological and epidemiological knowledge into model design and assessment and suggest that doing so can inform tailored modeling for individual drugs, pathogens, and clinical populations. We further suggest that continued comprehensive sampling and incorporation of up-to-date whole genome sequence data, resistance phenotypes, and treatment outcome data into model training will be crucial to the clinical utility and sustainability of machine learning-based molecular diagnostics.

Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis.

Modern strains of Mycobacterium tuberculosis from the Americas are closely related to those from Europe, supporting the assumption that human tuberculosis was introduced post-contact. This notion, however, is incompatible with archaeological evidence of pre-contact tuberculosis in the New World. Comparative genomics of modern isolates suggests that M. tuberculosis attained its worldwide distribution following human dispersals out of Africa during the Pleistocene epoch, although this has yet to be confirmed with ancient calibration points. Here we present three 1,000-year-old mycobacterial genomes from Peruvian human skeletons, revealing that a member of the M. tuberculosis complex caused human disease before contact. The ancient strains are distinct from known human-adapted forms and are most closely related to those adapted to seals and sea lions. Two independent dating approaches suggest a most recent common ancestor for the M. tuberculosis complex less than 6,000 years ago, which supports a Holocene dispersal of the disease. Our results implicate sea mammals as having played a role in transmitting the disease to humans across the ocean.

Bayesian inference of ancestral dates on bacterial phylogenetic trees.

The sequencing and comparative analysis of a collection of bacterial genomes from a single species or lineage of interest can lead to key insights into its evolution, ecology or epidemiology. The tool of choice for such a study is often to build a phylogenetic tree, and more specifically when possible a dated phylogeny, in which the dates of all common ancestors are estimated. Here, we propose a new Bayesian methodology to construct dated phylogenies which is specifically designed for bacterial genomics. Unlike previous Bayesian methods aimed at building dated phylogenies, we consider that the phylogenetic relationships between the genomes have been previously evaluated using a standard phylogenetic method, which makes our methodology much faster and scalable. This two-step approach also allows us to directly exploit existing phylogenetic methods that detect bacterial recombination, and therefore to account for the effect of recombination in the construction of a dated phylogeny. We analysed many simulated datasets in order to benchmark the performance of our approach in a wide range of situations. Furthermore, we present applications to three different real datasets from recent bacterial genomic studies. Our methodology is implemented in a R package called BactDating which is freely available for download at https://github.com/xavierdidelot/BactDating.

Correction: Dynamics and impact of homologous recombination on the evolution of Legionella pneumophila.

[This corrects the article DOI: 10.1371/journal.pgen.1006855.].

Dynamics and impact of homologous recombination on the evolution of Legionella pneumophila.

Legionella pneumophila is an environmental bacterium and the causative agent of Legionnaires' disease. Previous genomic studies have shown that recombination accounts for a high proportion (>96%) of diversity within several major disease-associated sequence types (STs) of L. pneumophila. This suggests that recombination represents a potentially important force shaping adaptation and virulence. Despite this, little is known about the biological effects of recombination in L. pneumophila, particularly with regards to homologous recombination (whereby genes are replaced with alternative allelic variants). Using newly available population genomic data, we have disentangled events arising from homologous and non-homologous recombination in six major disease-associated STs of L. pneumophila (subsp. pneumophila), and subsequently performed a detailed characterisation of the dynamics and impact of homologous recombination. We identified genomic "hotspots" of homologous recombination that include regions containing outer membrane proteins, the lipopolysaccharide (LPS) region and Dot/Icm effectors, which provide interesting clues to the selection pressures faced by L. pneumophila. Inference of the origin of the recombined regions showed that isolates have most frequently imported DNA from isolates belonging to their own clade, but also occasionally from other major clades of the same subspecies. This supports the hypothesis that the possibility for horizontal exchange of new adaptations between major clades of the subspecies may have been a critical factor in the recent emergence of several clinically important STs from diverse genomic backgrounds. However, acquisition of recombined regions from another subspecies, L. pneumophila subsp. fraseri, was rarely observed, suggesting the existence of a recombination barrier and/or the possibility of ongoing speciation between the two subspecies. Finally, we suggest that multi-fragment recombination may occur in L. pneumophila, whereby multiple non-contiguous segments that originate from the same molecule of donor DNA are imported into a recipient genome during a single episode of recombination.

Interacting networks of resistance, virulence and core machinery genes identified by genome-wide epistasis analysis.

Recent advances in the scale and diversity of population genomic datasets for bacteria now provide the potential for genome-wide patterns of co-evolution to be studied at the resolution of individual bases. Here we describe a new statistical method, genomeDCA, which uses recent advances in computational structural biology to identify the polymorphic loci under the strongest co-evolutionary pressures. We apply genomeDCA to two large population data sets representing the major human pathogens Streptococcus pneumoniae (pneumococcus) and Streptococcus pyogenes (group A Streptococcus). For pneumococcus we identified 5,199 putative epistatic interactions between 1,936 sites. Over three-quarters of the links were between sites within the pbp2x, pbp1a and pbp2b genes, the sequences of which are critical in determining non-susceptibility to beta-lactam antibiotics. A network-based analysis found these genes were also coupled to that encoding dihydrofolate reductase, changes to which underlie trimethoprim resistance. Distinct from these antibiotic resistance genes, a large network component of 384 protein coding sequences encompassed many genes critical in basic cellular functions, while another distinct component included genes associated with virulence. The group A Streptococcus (GAS) data set population represents a clonal population with relatively little genetic variation and a high level of linkage disequilibrium across the genome. Despite this, we were able to pinpoint two RNA pseudouridine synthases, which were each strongly linked to a separate set of loci across the chromosome, representing biologically plausible targets of co-selection. The population genomic analysis method applied here identifies statistically significantly co-evolving locus pairs, potentially arising from fitness selection interdependence reflecting underlying protein-protein interactions, or genes whose product activities contribute to the same phenotype. This discovery approach greatly enhances the future potential of epistasis analysis for systems biology, and can complement genome-wide association studies as a means of formulating hypotheses for targeted experimental work.