Jeffrey T. Foster

Staphylococcus aureus CC398: Host Adaptation and Emergence of Methicillin Resistance in Livestock

ABSTRACT Since its discovery in the early 2000s, methicillin-resistant Staphylococcus aureus (MRSA) clonal complex 398 (CC398) has become a rapidly emerging cause of human infections, most often associated with livestock exposure. We applied whole-genome sequence typing to characterize a diverse collection of CC398 isolates (n = 89), including MRSA and methicillin-susceptible S. aureus (MSSA) from animals and humans spanning 19 countries and four continents. We identified 4,238 single nucleotide polymorphisms (SNPs) among the 89 core genomes. Minimal homoplasy (consistency index = 0.9591) was detected among parsimony-informative SNPs, allowing for the generation of a highly accurate phylogenetic reconstruction of the CC398 clonal lineage. Phylogenetic analyses revealed that MSSA from humans formed the most ancestral clades. The most derived lineages were composed predominantly of livestock-associated MRSA possessing three different staphylococcal cassette chromosome mec element (SCCmec) types (IV, V, and VII-like) including nine subtypes. The human-associated isolates from the basal clades carried phages encoding human innate immune modulators that were largely missing among the livestock-associated isolates. Our results strongly suggest that livestock-associated MRSA CC398 originated in humans as MSSA. The lineage appears to have undergone a rapid radiation in conjunction with the jump from humans to livestock, where it subsequently acquired tetracycline and methicillin resistance. Further analyses are required to estimate the number of independent genetic events leading to the methicillin-resistant sublineages, but the diversity of SCCmec subtypes is suggestive of strong and diverse antimicrobial selection associated with food animal production. IMPORTANCE Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production. Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production.

Multidrug-Resistant Staphylococcus aureus in US Meat and Poultry

We characterized the prevalence, antibiotic susceptibility profiles, and genotypes of Staphylococcus aureus among US meat and poultry samples (n = 136). S. aureus contaminated 47% of samples, and multidrug resistance was common among isolates (52%). S. aureus genotypes and resistance profiles differed significantly among sample types, suggesting food animal–specific contamination.

Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer.

BackgroundPhylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction.ResultsBayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallaces Line, a biogeographic boundary between Southeast Asia and Australia.ConclusionWe describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

Whole-Genome-Based Phylogeny and Divergence of the Genus Brucella

Brucellae are worldwide bacterial pathogens of livestock and wildlife, but phylogenetic reconstructions have been challenging due to limited genetic diversity. We assessed the taxonomic and evolutionary relationships of five Brucella species-Brucella abortus, B. melitensis, B. suis, B. canis, and B. ovis-using whole-genome comparisons. We developed a phylogeny using single nucleotide polymorphisms (SNPs) from 13 genomes and rooted the tree using the closely related soil bacterium and opportunistic human pathogen, Ochrobactrum anthropi. Whole-genome sequencing and a SNP-based approach provided the requisite level of genetic detail to resolve species in the highly conserved brucellae. Comparisons among the Brucella genomes revealed 20,154 orthologous SNPs that were shared in all genomes. Rooting with Ochrobactrum anthropi reveals that the B. ovis lineage is basal to the rest of the Brucella lineage. We found that B. suis is a highly divergent clade with extensive intraspecific genetic diversity. Furthermore, B. suis was determined to be paraphyletic in our analyses, only forming a monophyletic clade when the B. canis genome was included. Using a molecular clock with these data suggests that most Brucella species diverged from their common B. ovis ancestor in the past 86,000 to 296,000 years, which precedes the domestication of their livestock hosts. Detailed knowledge of the Brucella phylogeny will lead to an improved understanding of the ecology, evolutionary history, and host relationships for this genus and can be used for determining appropriate genotyping approaches for rapid detection and diagnostic assays for molecular epidemiological and clinical studies.

Phylogenetic understanding of clonal populations in an era of whole genome sequencing

Phylogenetic hypotheses using whole genome sequences have the potential for unprecedented accuracy, yet a failure to understand issues associated with discovery bias, character sampling, and strain sampling can lead to highly erroneous conclusions. For microbial pathogens, phylogenies derived from whole genome sequences are becoming more common, as large numbers of characters distributed across entire genomes can yield extremely accurate phylogenies, particularly for strictly clonal populations. The availability of whole genomes is increasing as new sequencing technologies reduce the cost and time required for genome sequencing. Until entire sample collections can be fully sequenced, harnessing the phylogenetic power from whole genome sequences in more than a small subset of fully sequenced strains requires the integration of whole genome and partial genome genotyping data. Such integration involves discovering evolutionarily stable polymorphic characters by whole genome comparisons, then determining allelic states across a wide panel of isolates using high-throughput genotyping technologies. Here, we demonstrate how such an approach using single nucleotide polymorphisms (SNPs) yields highly accurate, but biased phylogenetic reconstructions and how the accuracy of the resulting tree is compromised by incomplete taxon and character sampling. Despite recent phylogenetic work detailing the strengths and biases of integrating whole genome and partial genome genotype data, these issues are relatively new and remain poorly understood by many researchers. Here, we revisit these biases and provide strategies for maximizing phylogenetic accuracy. Although we write this review with bacterial pathogens in mind, these concepts apply to any clonally reproducing population or indeed to any evolutionarily stable marker that is inherited in a strictly clonal manner. Understanding the ways in which current and emerging technologies can be used to maximize phylogenetic knowledge is advantageous only with a complete understanding of the strengths and weaknesses of these methods.

Host and pathogen ecology drive the seasonal dynamics of a fungal disease, white-nose syndrome

Seasonal patterns in pathogen transmission can influence the impact of disease on populations and the speed of spatial spread. Increases in host contact rates or births drive seasonal epidemics in some systems, but other factors may occasionally override these influences. White-nose syndrome, caused by the emerging fungal pathogen Pseudogymnoascus destructans, is spreading across North America and threatens several bat species with extinction. We examined patterns and drivers of seasonal transmission of P. destructans by measuring infection prevalence and pathogen loads in six bat species at 30 sites across the eastern United States. Bats became transiently infected in autumn, and transmission spiked in early winter when bats began hibernating. Nearly all bats in six species became infected by late winter when infection intensity peaked. In summer, despite high contact rates and a birth pulse, most bats cleared infections and prevalence dropped to zero. These data suggest the dominant driver of seasonal transmission dynamics was a change in host physiology, specifically hibernation. Our study is the first, to the best of our knowledge, to describe the seasonality of transmission in this emerging wildlife disease. The timing of infection and fungal growth resulted in maximal population impacts, but only moderate rates of spatial spread.

Development and validation of Burkholderia pseudomallei-specific real-time PCR assays for clinical, environmental or forensic detection applications.

The bacterium Burkholderia pseudomallei causes melioidosis, a rare but serious illness that can be fatal if untreated or misdiagnosed. Species-specific PCR assays provide a technically simple method for differentiating B. pseudomallei from near-neighbor species. However, substantial genetic diversity and high levels of recombination within this species reduce the likelihood that molecular signatures will differentiate all B. pseudomallei from other Burkholderiaceae. Currently available molecular assays for B. pseudomallei detection lack rigorous validation across large in silico datasets and isolate collections to test for specificity, and none have been subjected to stringent quality control criteria (accuracy, precision, selectivity, limit of quantitation (LoQ), limit of detection (LoD), linearity, ruggedness and robustness) to determine their suitability for environmental, clinical or forensic investigations. In this study, we developed two novel B. pseudomallei specific assays, 122018 and 266152, using a dual-probe approach to differentiate B. pseudomallei from B. thailandensis, B. oklahomensis and B. thailandensis-like species; other species failed to amplify. Species specificity was validated across a large DNA panel (>2,300 samples) comprising Burkholderia spp. and non-Burkholderia bacterial and fungal species of clinical and environmental relevance. Comparison of assay specificity to two previously published B. pseudomallei-specific assays, BurkDiff and TTS1, demonstrated comparable performance of all assays, providing between 99.7 and 100% specificity against our isolate panel. Last, we subjected 122018 and 266152 to rigorous quality control analyses, thus providing quantitative limits of assay performance. Using B. pseudomallei as a model, our study provides a framework for comprehensive quantitative validation of molecular assays and provides additional, highly validated B. pseudomallei assays for the scientific research community.

Single nucleotide polymorphisms for assessing genetic diversity in castor bean ( Ricinus communis )

BackgroundCastor bean (Ricinus communis) is an agricultural crop and garden ornamental that is widely cultivated and has been introduced worldwide. Understanding population structure and the distribution of castor bean cultivars has been challenging because of limited genetic variability. We analyzed the population genetics of R. communis in a worldwide collection of plants from germplasm and from naturalized populations in Florida, U.S. To assess genetic diversity we conducted survey sequencing of the genomes of seven diverse cultivars and compared the data to a reference genome assembly of a widespread cultivar (Hale). We determined the population genetic structure of 676 samples using single nucleotide polymorphisms (SNPs) at 48 loci.ResultsBayesian clustering indicated five main groups worldwide and a repeated pattern of mixed genotypes in most countries. High levels of population differentiation occurred between most populations but this structure was not geographically based. Most molecular variance occurred within populations (74%) followed by 22% among populations, and 4% among continents. Samples from naturalized populations in Florida indicated significant population structuring consistent with local demes. There was significant population differentiation for 56 of 78 comparisons in Florida (pairwise population ϕPT values, p < 0.01).ConclusionLow levels of genetic diversity and mixing of genotypes have led to minimal geographic structuring of castor bean populations worldwide. Relatively few lineages occur and these are widely distributed. Our approach of determining population genetic structure using SNPs from genome-wide comparisons constitutes a framework for high-throughput analyses of genetic diversity in plants, particularly in species with limited genetic diversity.

Genetic structure and evolved malaria resistance in Hawaiian honeycreepers

Infectious diseases now threaten wildlife populations worldwide but population recovery following local extinction has rarely been observed. In such a case, do resistant individuals recolonize from a central remnant population, or do they spread from small, perhaps overlooked, populations of resistant individuals? Introduced avian malaria (Plasmodium relictum) has devastated low‐elevation populations of native birds in Hawaii, but at least one species (Hawaii amakihi, Hemignathus virens) that was greatly reduced at elevations below about 1000 m tolerates malaria and has initiated a remarkable and rapid recovery. We assessed mitochondrial and nuclear DNA markers from amakihi and two other Hawaiian honeycreepers, apapane (Himatione sanguinea) and iiwi (Vestiaria coccinea), at nine primary study sites from 2001 to 2003 to determine the source of re‐establishing birds. In addition, we obtained sequences from tissue from amakihi museum study skins (1898 and 1948–49) to assess temporal changes in allele distributions. We found that amakihi in lowland areas are, and have historically been, differentiated from birds at high elevations and had unique alleles retained through time; that is, their genetic signature was not a subset of the genetic variation at higher elevations. We suggest that high disease pressure rapidly selected for resistance to malaria at low elevation, leaving small pockets of resistant birds, and this resistance spread outward from the scattered remnant populations. Low‐elevation amakihi are currently isolated from higher elevations (> 1000 m) where disease emergence and transmission rates appear to vary seasonally and annually. In contrast to results from amakihi, no genetic differentiation between elevations was found in apapane and iiwi, indicating that slight variation in genetic or life‐history attributes can determine disease resistance and population recovery. Determining the conditions that allow for the development of resistance to disease is essential to understanding how species evolve resistance across a landscape of varying disease pressures.

Genomics in Conservation: Case Studies and Bridging the Gap between Data and Application

We agree with Shafer et al. [1] that there is a need for well-documented case studies of the application of genomics in conservation and management as well as increased communication between academics and natural resource managers. However, we challenge Shafer et al.’s [1] relatively pessimistic assertion that ‘conservation genomics is far from seeing regular application’. Here we illustrate by examples that conservation practitioners utilize more genomic research than is often apparent. In addition, we highlight the work of nonacademic laboratories [government and nongovernmental organizations (NGOs)], some of which are not always well represented in peer-reviewed literature.