Using source-associated mobile genetic elements to identify zoonotic extraintestinal E. coli infections.

A one-health perspective may provide new and actionable information about Escherichia coli transmission. E. coli colonizes a broad range of vertebrates, including humans and food-production animals, and is a leading cause of bladder, kidney, and bloodstream infections in humans. Substantial evidence supports foodborne transmission of pathogenic E. coli strains from food animals to humans. However, the relative contribution of foodborne zoonotic E. coli (FZEC) to the human extraintestinal disease burden and the distinguishing characteristics of such strains remain undefined. Using a comparative genomic analysis of a large collection of contemporaneous, geographically-matched clinical and meat-source E. coli isolates (n = 3111), we identified 17 source-associated mobile genetic elements - predominantly plasmids and bacteriophages - and integrated them into a novel Bayesian latent class model to predict the origins of clinical E. coli isolates. We estimated that approximately 8 % of human extraintestinal E. coli infections (mostly urinary tract infections) in our study population were caused by FZEC. FZEC strains were equally likely to cause symptomatic disease as non-FZEC strains. Two FZEC lineages, ST131-H22 and ST58, appeared to have particularly high virulence potential. Our findings imply that FZEC strains collectively cause more urinary tract infections than does any single non-E. coli uropathogenic species (e.g., Klebsiella pneumoniae). Our novel approach can be applied in other settings to identify the highest-risk FZEC strains, determine their sources, and inform new one-health strategies to decrease the heavy public health burden imposed by extraintestinal E. coli infections.

Factor Xa cleaves SARS-CoV-2 spike protein to block viral entry and infection.

Serine proteases (SP), including furin, trypsin, and TMPRSS2 cleave the SARS-CoV-2 spike (S) protein, enabling the virus to enter cells. Here, we show that factor (F) Xa, an SP involved in blood coagulation, is upregulated in COVID-19 patients. In contrast to other SPs, FXa exerts antiviral activity. Mechanistically, FXa cleaves S protein, preventing its binding to ACE2, and thus blocking viral entry and infection. However, FXa is less effective against variants carrying the D614G mutation common in all pandemic variants. The anticoagulant rivaroxaban, a direct FXa inhibitor, inhibits FXa-mediated S protein cleavage and facilitates viral entry, whereas the indirect FXa inhibitor fondaparinux does not. In the lethal SARS-CoV-2 K18-hACE2 model, FXa prolongs survival yet its combination with rivaroxaban but not fondaparinux abrogates that protection. These results identify both a previously unknown function for FXa and an associated antiviral host defense mechanism against SARS-CoV-2 and suggest caution in considering direct FXa inhibitors for preventing or treating thrombotic complications in COVID-19 patients.

Multi-faceted metagenomic analysis of spacecraft associated surfaces reveal planetary protection relevant microbial composition.

The National Aeronautics and Space Administration (NASA) has been monitoring the microbial burden of spacecraft since the 1970's Viking missions. Originally culture-based and then focused 16S sequencing techniques were used, but we have now applied whole metagenomic sequencing to a variety of cleanroom samples at the Jet Propulsion Lab (JPL), including the Spacecraft Assembly Facility (SAF) with the goals of taxonomic identification and for functional assignment. Our samples included facility pre-filters, cleanroom vacuum debris, and surface wipes. The taxonomic composition was carried out by three different analysis tools to contrast marker, k-mer, and true alignment approaches. Hierarchical clustering analysis of the data separated vacuum particles from other SAF DNA samples. Vacuum particle samples were the most diverse while DNA samples from the ISO (International Standards Organization) compliant facilities and the SAF were the least diverse; all three were dominated by Proteobacteria. Wipe samples had higher diversity and were predominated by Actinobacteria, including human commensals Cutibacterium acnes and Corynebacterium spp. Taxa identified by the three methods were not identical, supporting the use of multiple methods for metagenome characterization. Likewise, functional annotation was performed using multiple methods. Vacuum particles and SAF samples contained strong signals of the tricarboxylic acid cycle and of amino acid biosynthesis, suggesting that many of the identified microorganisms have the ability to grow in nutrient-limited environments. In total, 18 samples generated high quality metagenome assembled genomes (MAG), which were dominated by Moraxella osloensis or Malassezia restricta. One M. osloensis MAG was assembled into a single circular scaffold and gene annotated. This study includes a rigorous quantitative determination of microbial loads and a qualitative dissection of microbial composition. Assembly of multiple specimens led to greater confidence for the identification of particular species and their predicted functional roles.

Genomic characterization of Francisella tularensis and other diverse Francisella species from complex samples.

Francisella tularensis, the bacterium that causes the zoonosis tularemia, and its genetic near neighbor species, can be difficult or impossible to cultivate from complex samples. Thus, there is a lack of genomic information for these species that has, among other things, limited the development of robust detection assays for F. tularensis that are both specific and sensitive. The objective of this study was to develop and validate approaches to capture, enrich, sequence, and analyze Francisella DNA present in DNA extracts generated from complex samples. RNA capture probes were designed based upon the known pan genome of F. tularensis and other diverse species in the family Francisellaceae. Probes that targeted genomic regions also present in non-Francisellaceae species were excluded, and probes specific to particular Francisella species or phylogenetic clades were identified. The capture-enrichment system was then applied to diverse, complex DNA extracts containing low-level Francisella DNA, including human clinical tularemia samples, environmental samples (i.e., animal tissue and air filters), and whole ticks/tick cell lines, which was followed by sequencing of the enriched samples. Analysis of the resulting data facilitated rigorous and unambiguous confirmation of the detection of F. tularensis or other Francisella species in complex samples, identification of mixtures of different Francisella species in the same sample, analysis of gene content (e.g., known virulence and antimicrobial resistance loci), and high-resolution whole genome-based genotyping. The benefits of this capture-enrichment system include: even very low target DNA can be amplified; it is culture-independent, reducing exposure for research and/or clinical personnel and allowing genomic information to be obtained from samples that do not yield isolates; and the resulting comprehensive data not only provide robust means to confirm the presence of a target species in a sample, but also can provide data useful for source attribution, which is important from a genomic epidemiology perspective.

Stenoparib, an inhibitor of cellular poly (ADP-ribose) polymerases (PARPs), blocks in vitro replication of SARS-CoV-2 variants.

We recently published a preliminary assessment of the activity of a poly (ADP-ribose) polymerase (PARP) inhibitor, stenoparib, also known as 2X-121, which inhibits viral replication by affecting pathways of the host. Here we show that stenoparib effectively inhibits a SARS-CoV-2 wild type (BavPat1/2020) strain and four additional variant strains; alpha (B.1.1.7), beta (B.1.351), delta (B.1.617.2) and gamma (P.1) in vitro, with 50% effective concentration (EC50) estimates of 4.1 µM, 8.5 µM, 24.1 µM, 8.2 µM and 13.6 µM, respectively. A separate experiment focusing on a combination of 10 µM stenoparib and 0.5 µM remdesivir, an antiviral drug, resulted in over 80% inhibition of the alpha variant, which is substantially greater than the effect achieved with either drug alone, suggesting at least additive effects from combining the different mechanisms of activity of stenoparib and remdesivir.

COVID-19 vaccination elicits an evolving, cross-reactive antibody response to epitopes conserved with endemic coronavirus spike proteins.

The COVID-19 pandemic has triggered the first widespread vaccination campaign against a coronavirus. Many vaccinated subjects are previously naive to SARS-CoV-2; however, almost all have previously encountered other coronaviruses (CoVs), and the role of this immunity in shaping the vaccine response remains uncharacterized. Here, we use longitudinal samples and highly multiplexed serology to identify mRNA-1273 vaccine-induced antibody responses against a range of CoV Spike epitopes, in both phylogenetically conserved and non-conserved regions. Whereas reactivity to SARS-CoV-2 epitopes shows a delayed but progressive increase following vaccination, we observe distinct kinetics for the endemic CoV homologs at conserved sites in Spike S2: these become detectable sooner and decay at later time points. Using homolog-specific antibody depletion and alanine-substitution experiments, we show that these distinct trajectories reflect an evolving cross-reactive response that can distinguish rare, polymorphic residues within these epitopes. Our results reveal mechanisms for the formation of antibodies with broad reactivity against CoVs.

Off-the-shelf CAR natural killer cells secreting IL-15 target spike in treating COVID-19.

Engineered natural killer (NK) cells represent a promising option for immune therapy option due to their immediate availability in allogeneic settings. Severe acute diseases, such as COVID-19, require targeted and immediate intervention. Here we show engineering of NK cells to express (1) soluble interleukin-15 (sIL15) for enhancing their survival and (2) a chimeric antigen receptor (CAR) consisting of an extracellular domain of ACE2, targeting the spike protein of SARS-CoV-2. These CAR NK cells (mACE2-CAR_sIL15 NK cells) bind to VSV-SARS-CoV-2 chimeric viral particles as well as the recombinant SARS-CoV-2 spike protein subunit S1 leading to enhanced NK cell production of TNF-α and IFN-γ and increased in vitro and in vivo cytotoxicity against cells expressing the spike protein. Administration of mACE2-CAR_sIL15 NK cells maintains body weight, reduces viral load, and prolongs survival of transgenic mice expressing human ACE2 upon infection with live SARS-CoV-2. These experiments, and the capacity of mACE2-CAR_sIL15 NK cells to retain their activity following cryopreservation, demonstrate their potential as an allogeneic off-the-shelf therapy for COVID-19 patients who are faced with limited treatment options.

COVID-19 vaccination recruits and matures cross-reactive antibodies to conserved epitopes in endemic coronavirus Spike proteins.

The COVID-19 pandemic has triggered the first widespread vaccination campaign against a coronavirus. Most vaccinated subjects are naïve to SARS-CoV-2, however almost all have previously encountered other coronaviruses (CoVs) and the role of this immunity in shaping the vaccine response remains uncharacterized. Here we use longitudinal samples and highly-multiplexed serology to identify mRNA-1273 vaccine-induced antibody responses against a range of CoV Spike epitopes and in both phylogenetically conserved and non-conserved regions. Whereas reactivity to SARS-CoV-2 epitopes showed a delayed but progressive increase following vaccination, we observed distinct kinetics for the endemic CoV homologs at two conserved sites in Spike S2: these became detectable sooner, and decayed at later timepoints. Using homolog-specific depletion and alanine-substitution experiments, we show that these distinctly-evolving specificities result from cross-reactive antibodies as they mature against rare, polymorphic residues within these epitopes. Our results reveal mechanisms for the formation of antibodies with broad reactivity against CoVs.

The K18-Human ACE2 Transgenic Mouse Model Recapitulates Non-severe and Severe COVID-19 in Response to an Infectious Dose of the SARS-CoV-2 Virus.

A comprehensive analysis and characterization of a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection model that mimics non-severe and severe coronavirus disease 2019 (COVID-19) in humans is warranted for understating the virus and developing preventive and therapeutic agents. Here, we characterized the K18-hACE2 mouse model expressing human (h)ACE2 in mice, controlled by the human keratin 18 (K18) promoter, in the epithelia, including airway epithelial cells where SARS-CoV-2 infections typically start. We found that intranasal inoculation with higher viral doses (2 × 103 and 2 × 104 PFU) of SARS-CoV-2 caused lethality of all mice and severe damage of various organs, including lung, liver, and kidney, while lower doses (2 × 101 and 2 × 102 PFU) led to less severe tissue damage and some mice recovered from the infection. In this hACE2 mouse model, SARS-CoV-2 infection damaged multiple tissues, with a dose-dependent effect in most tissues. Similar damage was observed in postmortem samples from COVID-19 patients. Finally, the mice that recovered from infection with a low dose of virus survived rechallenge with a high dose of virus. Compared to other existing models, the K18-hACE2 model seems to be the most sensitive COVID-19 model reported to date. Our work expands the information available about this model to include analysis of multiple infectious doses and various tissues with comparison to human postmortem samples from COVID-19 patients. In conclusion, the K18-hACE2 mouse model recapitulates both severe and non-severe COVID-19 in humans being dose-dependent and can provide insight into disease progression and the efficacy of therapeutics for preventing or treating COVID-19. IMPORTANCE The pandemic of coronavirus disease 2019 (COVID-19) has reached nearly 240 million cases, caused nearly 5 million deaths worldwide as of October 2021, and has raised an urgent need for the development of novel drugs and therapeutics to prevent the spread and pathogenesis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To achieve this goal, an animal model that recapitulates the features of human COVID-19 disease progress and pathogenesis is greatly needed. In this study, we have comprehensively characterized a mouse model of SARS-CoV-2 infection using K18-hACE2 transgenic mice. We infected the mice with low and high doses of SARS-CoV-2 to study the pathogenesis and survival in response to different infection patterns. Moreover, we compared the pathogenesis of the K18-hACE2 transgenic mice with that of the COVID-19 patients to show that this model could be a useful tool for the development of antiviral drugs and therapeutics.

FXa cleaves the SARS-CoV-2 spike protein and blocks cell entry to protect against infection with inferior effects in B.1.1.7 variant.

The ongoing coronavirus disease 2019 (COVID-19) pandemic is caused by infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Human natural defense mechanisms against SARS-CoV-2 are largely unknown. Serine proteases (SPs) including furin and TMPRSS2 cleave SARS-CoV-2 spike protein, facilitating viral entry. Here, we show that FXa, a SP for blood coagulation, is upregulated in COVID-19 patients compared to non-COVID-19 donors and exerts anti-viral activity. Mechanistically, FXa cleaves the SARS-CoV-2 spike protein, which prevents its binding to ACE2, and thus blocks viral entry. Furthermore, the variant B.1.1.7 with several mutations is dramatically resistant to the anti-viral effect of FXa compared to wild-type SARA-CoV-2 in vivo and in vitro. The anti-coagulant rivaroxaban directly inhibits FXa and facilitates viral entry, whereas the indirect inhibitor fondaparinux does not. In a lethal humanized hACE2 mouse model of SARS-CoV-2, FXa prolonged survival while combination with rivaroxaban but not fondaparinux abrogated this protection. These preclinical results identify a previously unknown SP function and associated anti-viral host defense mechanism and suggest caution in considering direct inhibitors for prevention or treatment of thrombotic complications in COVID-19 patients.

Dual roles of a novel oncolytic viral vector-based SARS-CoV-2 vaccine: preventing COVID-19 and treating tumor progression.

The ongoing coronavirus disease 2019 (COVID-19) pandemic is caused by infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Cancer patients are usually immunocompromised and thus are particularly susceptible to SARS-CoV-2 infection resulting in COVID-19. Although many vaccines against COVID-19 are being preclinically or clinically tested or approved, none have yet been specifically developed for cancer patients or reported as having potential dual functions to prevent COVID-19 and treat cancer. Here, we confirmed that COVID-19 patients with cancer have low levels of antibodies against the spike (S) protein, a viral surface protein mediating the entry of SARS-CoV-2 into host cells, compared with COVID-19 patients without cancer. We developed an oncolytic herpes simplex virus-1 vector-based vaccine named oncolytic virus (OV)-spike. OV-spike induced abundant anti-S protein neutralization antibodies in both tumor-free and tumor-bearing mice, which inhibit infection of VSV-SARS-CoV-2 and wild-type (WT) live SARS-CoV-2 as well as the B.1.1.7 variant in vitro. In the tumor-bearing mice, OV-spike also inhibited tumor growth, leading to better survival in multiple preclinical tumor models than the untreated control. Furthermore, OV-spike induced anti-tumor immune response and SARS-CoV-2-specific T cell response without causing serious adverse events. Thus, OV-spike is a promising vaccine candidate for both preventing COVID-19 and enhancing the anti-tumor response. ONE SENTENCE SUMMARY: A herpes oncolytic viral vector-based vaccine is a promising vaccine with dual roles in preventing COVID-19 and treating tumor progression.

St. Louis Encephalitis Virus in the Southwestern United States: A Phylogeographic Case for a Multi-Variant Introduction Event.

Since the reemergence of St. Louis Encephalitis (SLE) Virus (SLEV) in the Southwest United States, identified during the 2015 outbreak in Arizona, SLEV has been seasonally detected within Culex spp. populations throughout the Southwest United States. Previous work revealed the 2015 outbreak was caused by an importation of SLEV genotype III, which had only been detected previously in Argentina. However, little is known about when the importation occurred or the transmission and genetic dynamics since its arrival into the Southwest. In this study, we sought to determine whether the annual detection of SLEV in the Southwest is due to enzootic cycling or new importations. To address this question, we analyzed 174 SLEV genomes (142 sequenced as part of this study) using Bayesian phylogenetic analyses to estimate the date of arrival into the American Southwest and characterize the underlying population structure of SLEV. Phylogenetic clustering showed that SLEV variants circulating in Maricopa and Riverside counties form two distinct populations with little evidence of inter-county transmission since the onset of the outbreak. Alternatively, it appears that in 2019, Yuma and Clark counties experienced annual importations of SLEV that originated in Riverside and Maricopa counties. Finally, the earliest representatives of SLEV genotype III in the Southwest form a polytomy that includes both California and Arizona samples. We propose that the initial outbreak most likely resulted from the importation of a population of SLEV genotype III variants, perhaps in multiple birds, possibly multiple species, migrating north in 2013, rather than a single variant introduced by one bird.

Roles of bacteriophages, plasmids and CRISPR immunity in microbial community dynamics revealed using time-series integrated meta-omics.

Viruses and plasmids (invasive mobile genetic elements (iMGEs)) have important roles in shaping microbial communities, but their dynamic interactions with CRISPR-based immunity remain unresolved. We analysed generation-resolved iMGE-host dynamics spanning one and a half years in a microbial consortium from a biological wastewater treatment plant using integrated meta-omics. We identified 31 bacterial metagenome-assembled genomes encoding complete CRISPR-Cas systems and their corresponding iMGEs. CRISPR-targeted plasmids outnumbered their bacteriophage counterparts by at least fivefold, highlighting the importance of CRISPR-mediated defence against plasmids. Linear modelling of our time-series data revealed that the variation in plasmid abundance over time explained more of the observed community dynamics than phages. Community-scale CRISPR-based plasmid-host and phage-host interaction networks revealed an increase in CRISPR-mediated interactions coinciding with a decrease in the dominant 'Candidatus Microthrix parvicella' population. Protospacers were enriched in sequences targeting genes involved in the transmission of iMGEs. Understanding the factors shaping the fitness of specific populations is necessary to devise control strategies for undesirable species and to predict or explain community-wide phenotypes.

Integration of time-series meta-omics data reveals how microbial ecosystems respond to disturbance.

The development of reliable, mixed-culture biotechnological processes hinges on understanding how microbial ecosystems respond to disturbances. Here we reveal extensive phenotypic plasticity and niche complementarity in oleaginous microbial populations from a biological wastewater treatment plant. We perform meta-omics analyses (metagenomics, metatranscriptomics, metaproteomics and metabolomics) on in situ samples over 14 months at weekly intervals. Based on 1,364 de novo metagenome-assembled genomes, we uncover four distinct fundamental niche types. Throughout the time-series, we observe a major, transient shift in community structure, coinciding with substrate availability changes. Functional omics data reveals extensive variation in gene expression and substrate usage amongst community members. Ex situ bioreactor experiments confirm that responses occur within five hours of a pulse disturbance, demonstrating rapid adaptation by specific populations. Our results show that community resistance and resilience are a function of phenotypic plasticity and niche complementarity, and set the foundation for future ecological engineering efforts.

Persistence of Burkholderia thailandensis E264 in lung tissue after a single binge alcohol episode.

BACKGROUND: Binge drinking, an increasingly common form of alcohol use disorder, is associated with substantial morbidity and mortality; yet, its effects on the immune system's ability to defend against infectious agents are poorly understood. Burkholderia pseudomallei, the causative agent of melioidosis can occur in healthy humans, yet binge alcohol intoxication is increasingly being recognized as a major risk factor. Although our previous studies demonstrated that binge alcohol exposure increased B. pseudomallei near-neighbor virulence in vivo and increased paracellular diffusion and intracellular invasion, no experimental studies have examined the extent to which bacterial and alcohol dosage play a role in disease progression. In addition, the temporal effects of a single binge alcohol dose prior to infection has not been examined in vivo. PRINCIPAL FINDINGS: In this study, we used B. thailandensis E264 a close genetic relative of B. pseudomallei, as useful BSL-2 model system. Eight-week-old female C57BL/6 mice were utilized in three distinct animal models to address the effects of 1) bacterial dosage, 2) alcohol dosage, and 3) the temporal effects, of a single binge alcohol episode. Alcohol was administered comparable to human binge drinking (≤ 4.4 g/kg) or PBS intraperitoneally before a non-lethal intranasal infection. Bacterial colonization of lung and spleen was increased in mice administered alcohol even after bacterial dose was decreased 10-fold. Lung and not spleen tissue were colonized even after alcohol dosage was decreased 20 times below the U.S legal limit. Temporally, a single binge alcohol episode affected lung bacterial colonization for more than 24 h after alcohol was no longer detected in the blood. Pulmonary and splenic cytokine expression (TNF-α, GM-CSF) remained suppressed, while IL-12/p40 increased in mice administered alcohol 6 or 24 h prior to infection. Increased lung and not intestinal bacterial invasion was observed in human and murine non-phagocytic epithelial cells exposed to 0.2% v/v alcohol in vitro. CONCLUSIONS: Our results indicate that the effects of a single binge alcohol episode are tissue specific. A single binge alcohol intoxication event increases bacterial colonization in mouse lung tissue even after very low BACs and decreases the dose required to colonize the lungs with less virulent B. thailandensis. Additionally, the temporal effects of binge alcohol alters lung and spleen cytokine expression for at least 24 h after alcohol is detected in the blood. Delayed recovery in lung and not spleen tissue may provide a means for B. pseudomallei and near-neighbors to successfully colonize lung tissue through increased intracellular invasion of non-phagocytic cells in patients with hazardous alcohol intake.

Phylogenetic Analysis of Francisella tularensis Group A.II Isolates from 5 Patients with Tularemia, Arizona, USA, 2015-2017.

We examined 5 tularemia cases in Arizona, USA, during 2015-2017. All were caused by Francisella tularensis group A.II. Genetically similar isolates were found across large spatial and temporal distances, suggesting that group A.II strains are dispersed across long distances by wind and exhibit low replication rates in the environment.

Phylogenetic analysis of West Nile Virus in Maricopa County, Arizona: Evidence for dynamic behavior of strains in two major lineages in the American Southwest.

West Nile Virus (WNV) has been detected annually in Maricopa County, Arizona, since 2003. With this in mind, we sought to determine if contemporary strains are endemic to the county or are annually imported. As part of this effort, we developed a new protocol for tiled amplicon sequencing of WNV to efficiently attain greater than 99% coverage of 14 WNV genomes collected directly from positive mosquito pools distributed throughout Maricopa County between 2014 and 2017. Bayesian phylogenetic analyses revealed that contemporary genomes fall within two major lineages; NA/WN02 and SW/WN03. We found that all of the Arizona strains possessed an amino acid substitution known to be under positive selection, which has arisen independently at least four times in Arizona. The SW/WN03 strains exhibited transient behavior, with at least 10 separate introductions into Arizona when considering both historical and contemporary strains. However, NA/WN02 strains are geographically differentiated and appear to be endemic in Arizona, with two clades that have been circulating for four and seven years. This establishment in Maricopa County provides the first evidence of local overwintering by a WNV strain over the course of several years in Arizona. Within a national context, the placement of eleven contemporary Arizona strains in the NA/WN02 lineage indicates while WNV first entered the northeastern United States in 1999, the most ancestral extant strains of WNV are now circulating in the American southwest.

Phylogenetics of a Fungal Invasion: Origins and Widespread Dispersal of White-Nose Syndrome.

Globalization has facilitated the worldwide movement and introduction of pathogens, but epizoological reconstructions of these invasions are often hindered by limited sampling and insufficient genetic resolution among isolates. Pseudogymnoascus destructans, a fungal pathogen causing the epizootic of white-nose syndrome in North American bats, has exhibited few genetic polymorphisms in previous studies, presenting challenges for both epizoological tracking of the spread of this fungus and for determining its evolutionary history. We used single nucleotide polymorphisms (SNPs) from whole-genome sequencing and microsatellites to construct high-resolution phylogenies of P. destructans Shallow genetic diversity and the lack of geographic structuring among North American isolates support a recent introduction followed by expansion via clonal reproduction across the epizootic zone. Moreover, the genetic relationships of isolates within North America suggest widespread mixing and long-distance movement of the fungus. Genetic diversity among isolates of P. destructans from Europe was substantially higher than in those from North America. However, genetic distance between the North American isolates and any given European isolate was similar to the distance between the individual European isolates. In contrast, the isolates we examined from Asia were highly divergent from both European and North American isolates. Although the definitive source for introduction of the North American population has not been conclusively identified, our data support the origin of the North American invasion by P. destructans from Europe rather than Asia.IMPORTANCE This phylogenetic study of the bat white-nose syndrome agent, P. destructans, uses genomics to elucidate evolutionary relationships among populations of the fungal pathogen to understand the epizoology of this biological invasion. We analyze hypervariable and abundant genetic characters (microsatellites and genomic SNPs, respectively) to reveal previously uncharacterized diversity among populations of the pathogen from North America and Eurasia. We present new evidence supporting recent introduction of the fungus to North America from a diverse Eurasian population, with limited increase in genetic variation in North America since that introduction.

Low cost, low tech SNP genotyping tools for resource-limited areas: Plague in Madagascar as a model.

BACKGROUND: Genetic analysis of pathogenic organisms is a useful tool for linking human cases together and/or to potential environmental sources. The resulting data can also provide information on evolutionary patterns within a targeted species and phenotypic traits. However, the instruments often used to generate genotyping data, such as single nucleotide polymorphisms (SNPs), can be expensive and sometimes require advanced technologies to implement. This places many genotyping tools out of reach for laboratories that do not specialize in genetic studies and/or lack the requisite financial and technological resources. To address this issue, we developed a low cost and low tech genotyping system, termed agarose-MAMA, which combines traditional PCR and agarose gel electrophoresis to target phylogenetically informative SNPs. METHODOLOGY/PRINCIPAL FINDINGS: To demonstrate the utility of this approach for generating genotype data in a resource-constrained area (Madagascar), we designed an agarose-MAMA system targeting previously characterized SNPs within Yersinia pestis, the causative agent of plague. We then used this system to genetically type pathogenic strains of Y. pestis in a Malagasy laboratory not specialized in genetic studies, the Institut Pasteur de Madagascar (IPM). We conducted rigorous assay performance validations to assess potential variation introduced by differing research facilities, reagents, and personnel and found no difference in SNP genotyping results. These agarose-MAMA PCR assays are currently employed as an investigative tool at IPM, providing Malagasy researchers a means to improve the value of their plague epidemiological investigations by linking outbreaks to potential sources through genetic characterization of isolates and to improve understanding of disease ecology that may contribute to a long-term control effort. CONCLUSIONS: The success of our study demonstrates that the SNP-based genotyping capacity of laboratories in developing countries can be expanded with manageable financial cost for resource constraint laboratories. This is a practical formula that reduces resource-driven limitations to genetic research and promises to advance global collective knowledge of infectious diseases emanating from resource limited regions of the world.

Protein abundances can distinguish between naturally-occurring and laboratory strains of Yersinia pestis, the causative agent of plague.

The rapid pace of bacterial evolution enables organisms to adapt to the laboratory environment with repeated passage and thus diverge from naturally-occurring environmental ("wild") strains. Distinguishing wild and laboratory strains is clearly important for biodefense and bioforensics; however, DNA sequence data alone has thus far not provided a clear signature, perhaps due to lack of understanding of how diverse genome changes lead to convergent phenotypes, difficulty in detecting certain types of mutations, or perhaps because some adaptive modifications are epigenetic. Monitoring protein abundance, a molecular measure of phenotype, can overcome some of these difficulties. We have assembled a collection of Yersinia pestis proteomics datasets from our own published and unpublished work, and from a proteomics data archive, and demonstrated that protein abundance data can clearly distinguish laboratory-adapted from wild. We developed a lasso logistic regression classifier that uses binary (presence/absence) or quantitative protein abundance measures to predict whether a sample is laboratory-adapted or wild that proved to be ~98% accurate, as judged by replicated 10-fold cross-validation. Protein features selected by the classifier accord well with our previous study of laboratory adaptation in Y. pestis. The input data was derived from a variety of unrelated experiments and contained significant confounding variables. We show that the classifier is robust with respect to these variables. The methodology is able to discover signatures for laboratory facility and culture medium that are largely independent of the signature of laboratory adaptation. Going beyond our previous laboratory evolution study, this work suggests that proteomic differences between laboratory-adapted and wild Y. pestis are general, potentially pointing to a process that could apply to other species as well. Additionally, we show that proteomics datasets (even archived data collected for different purposes) contain the information necessary to distinguish wild and laboratory samples. This work has clear applications in biomarker detection as well as biodefense.