Search

It is not possible to systematically screen the environment for rabies virus (RABV) using current approaches. We sought to determine under what conditions RABV is detectable from feces and other accessible samples from infected wildlife to broaden the number of biological samples that could be used to test for RABV. We employed a recently-developed quantitative RT-PCR assay called the "LN34 panlyssavirus real-time RT-PCR assay", which is highly sensitive and specific for all variants of RABV. We harvested and tested brain tissue, fecal, and/or mouth swab samples from 25 confirmed RABV positive bats of six species. To determine if rabies RNA lasts in feces sufficiently long post-defecation to use it as a surveillance tool, we tested fecal samples from 10 bats at the time of sample collection and after 24 hours of exposure to ambient conditions, with an additional test on six bats out to 72 hours. To assess whether we could pool fecal pellets and still detect a positive, we generated dilutions of known positives at 1:1, 1:10, 1:50, and 1:200. For six individuals for which matched brain, mouth swab, and fecal samples were tested, results were positive for 100%, 67%, and 67%, respectively. For the first time test to 24 hours, 63% of feces that were positive at time 0 were still positive after 24 hours, and 50% of samples at 72 hours were positive across all three replicates. Pooling tests revealed that fecal positives were detected at 1:10 dilution, but not at 1:50 or 1:200. Our preliminary results suggest that fecal samples hold promise for a rapid and non-invasive environmental screening system.

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.

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.

Staphylococcus aureus is an important clinical pathogen worldwide and understanding this organism's phylogeny and, in particular, the role of recombination, is important both to understand the overall spread of virulent lineages and to characterize outbreaks. To further elucidate the phylogeny of S. aureus, 35 diverse strains were sequenced using whole genome sequencing. In addition, 29 publicly available whole genome sequences were included to create a single nucleotide polymorphism (SNP)-based phylogenetic tree encompassing 11 distinct lineages. All strains of a particular sequence type fell into the same clade with clear groupings of the major clonal complexes of CC8, CC5, CC30, CC45 and CC1. Using a novel analysis method, we plotted the homoplasy density and SNP density across the whole genome and found evidence of recombination throughout the entire chromosome, but when we examined individual clonal lineages we found very little recombination. However, when we analyzed three branches of multiple lineages, we saw intermediate and differing levels of recombination between them. These data demonstrate that in S. aureus, recombination occurs across major lineages that subsequently expand in a clonal manner. Estimated mutation rates for the CC8 and CC5 lineages were different from each other. While the CC8 lineage rate was similar to previous studies, the CC5 lineage was 100-fold greater. Fifty known virulence genes were screened in all genomes in silico to determine their distribution across major clades. Thirty-three genes were present variably across clades, most of which were not constrained by ancestry, indicating horizontal gene transfer or gene loss.

Undetected and untreated, low-levels of drug resistant (DR) subpopulations in clinical Mycobacterium tuberculosis (Mtb) infections may lead to development of DR-tuberculosis, potentially resulting in treatment failure. Current phenotypic DR susceptibility testing has a theoretical potential for 1% sensitivity, is not quantitative, and requires several weeks to complete. The use of "single molecule-overlapping reads" (SMOR) analysis with next generation DNA sequencing for determination of ultra-rare target alleles in complex mixtures provides increased sensitivity over standard DNA sequencing. Ligation free amplicon sequencing with SMOR analysis enables the detection of resistant allele subpopulations at ≥0.1% of the total Mtb population in near real-time analysis. We describe the method using standardized mixtures of DNA from resistant and susceptible Mtb isolates and the assay's performance for detecting ultra-rare DR subpopulations in DNA extracted directly from clinical sputum samples. SMOR analysis enables rapid near real-time detection and tracking of previously undetectable DR sub-populations in clinical samples allowing for the evaluation of the clinical relevance of low-level DR subpopulations. This will provide insights into interventions aimed at suppressing minor DR subpopulations before they become clinically significant.

Multidrug-resistant Klebsiella pneumoniae producing the KPC carbapenemase have rapidly spread throughout the world, causing severe healthcare-associated infections with limited antimicrobial treatment options. Dissemination of KPC-producing K. pneumoniae is largely attributed to expansion of a single dominant strain, ST258. In this study, we explore phylogenetic relationships and evolution within ST258 and its clonal group, CG258, using whole genome sequence analysis of 167 isolates from 20 countries collected over 17 years. Our results show a common ST258 ancestor emerged from its diverse parental clonal group around 1995 and likely acquired blaKPC prior to dissemination. Over the past two decades, ST258 has remained highly clonal despite diversity in accessory elements and divergence in the capsule polysaccharide synthesis locus. Apart from the large recombination event that gave rise to ST258, few mutations set it apart from its clonal group. However, one mutation occurs in a global transcription regulator. Characterization of outer membrane protein sequences revealed a profile in ST258 that includes a truncated OmpK35 and modified OmpK37. Our work illuminates potential genomic contributors to the pathogenic success of ST258, helps us better understand the global dissemination of this strain, and identifies genetic markers unique to ST258.

With the decreasing cost of next-generation sequencing, deep sequencing of clinical samples provides unique opportunities to understand host-associated microbial communities. Among the primary challenges of clinical metagenomic sequencing is the rapid filtering of human reads to survey for pathogens with high specificity and sensitivity. Metagenomes are inherently variable due to different microbes in the samples and their relative abundance, the size and architecture of genomes, and factors such as target DNA amounts in tissue samples (i.e. human DNA versus pathogen DNA concentration). This variation in metagenomes typically manifests in sequencing datasets as low pathogen abundance, a high number of host reads, and the presence of close relatives and complex microbial communities. In addition to these challenges posed by the composition of metagenomes, high numbers of reads generated from high-throughput deep sequencing pose immense computational challenges. Accurate identification of pathogens is confounded by individual reads mapping to multiple different reference genomes due to gene similarity in different taxa present in the community or close relatives in the reference database. Available global and local sequence aligners also vary in sensitivity, specificity, and speed of detection. The efficiency of detection of pathogens in clinical samples is largely dependent on the desired taxonomic resolution of the organisms. We have developed an efficient strategy that identifies "all against all" relationships between sequencing reads and reference genomes. Our approach allows for scaling to large reference databases and then genome reconstruction by aggregating global and local alignments, thus allowing genetic characterization of pathogens at higher taxonomic resolution. These results were consistent with strain level SNP genotyping and bacterial identification from laboratory culture.

Acinetobacter baumannii is an emergent and global nosocomial pathogen. In addition to A. baumannii, other Acinetobacter species, especially those in the Acinetobacter calcoaceticus-baumannii (Acb) complex, have also been associated with serious human infection. Although mechanisms of attachment, persistence on abiotic surfaces, and pathogenesis in A. baumannii have been identified, the genetic mechanisms that explain the emergence of A. baumannii as the most widespread and virulent Acinetobacter species are not fully understood. Recent whole genome sequencing has provided insight into the phylogenetic structure of the genus Acinetobacter. However, a global comparison of genomic features between Acinetobacter spp. has not been described in the literature. In this study, 136 Acinetobacter genomes, including 67 sequenced in this study, were compared to identify the acquisition and loss of genes in the expansion of the Acinetobacter genus. A whole genome phylogeny confirmed that A. baumannii is a monophyletic clade and that the larger Acb complex is also a well-supported monophyletic group. The whole genome phylogeny provided the framework for a global genomic comparison based on a blast score ratio (BSR) analysis. The BSR analysis demonstrated that specific genes have been both lost and acquired in the evolution of A. baumannii. In addition, several genes associated with A. baumannii pathogenesis were found to be more conserved in the Acb complex, and especially in A. baumannii, than in other Acinetobacter genomes; until recently, a global analysis of the distribution and conservation of virulence factors across the genus was not possible. The results demonstrate that the acquisition of specific virulence factors has likely contributed to the widespread persistence and virulence of A. baumannii. The identification of novel features associated with transcriptional regulation and acquired by clades in the Acb complex presents targets for better understanding the evolution of pathogenesis and virulence in the expansion of the genus.

Case reports of Apophysomyces spp. in immunocompetent hosts have been a result of traumatic deep implantation of Apophysomyces spp. spore-contaminated soil or debris. On May 22, 2011 a tornado occurred in Joplin, MO, leaving 13 tornado victims with Apophysomyces trapeziformis infections as a result of lacerations from airborne material. We used whole genome sequence typing (WGST) for high-resolution phylogenetic SNP analysis of 17 outbreak Apophysomyces isolates and five additional temporally and spatially diverse Apophysomyces control isolates (three A. trapeziformis and two A. variabilis isolates). Whole genome SNP phylogenetic analysis revealed three clusters of genotypically related or identical A. trapeziformis isolates and multiple distinct isolates among the Joplin group; this indicated multiple genotypes from a single or multiple sources. Though no linkage between genotype and location of exposure was observed, WGST analysis determined that the Joplin isolates were more closely related to each other than to the control isolates, suggesting local population structure. Additionally, species delineation based on WGST demonstrated the need to reassess currently accepted taxonomic classifications of phylogenetic species within the genus Apophysomyces.

Burkholderia pseudomallei is a gram-negative bacterium that causes the serious human disease, melioidosis. There is no vaccine against melioidosis and it can be fatal if not treated with a specific antibiotic regimen, which typically includes the third-generation cephalosporin, ceftazidime (CAZ). There have been several resistance mechanisms described for B. pseudomallei, of which the best described are amino acid changes that alter substrate specificity in the highly conserved class A β-lactamase, PenA. In the current study, we sequenced penA from isolates sequentially derived from two melioidosis patients with wild-type (1.5 µg/mL) and, subsequently, resistant (16 or ≥256 µg/mL) CAZ phenotypes. We identified two single-nucleotide polymorphisms (SNPs) that directly increased CAZ hydrolysis. One SNP caused an amino acid substitution (C69Y) near the active site of PenA, whereas a second novel SNP was found within the penA promoter region. In both instances, the CAZ resistance phenotype corresponded directly with the SNP genotype. Interestingly, these SNPs appeared after infection and under selection from CAZ chemotherapy. Through heterologous cloning and expression, and subsequent allelic exchange in the native bacterium, we confirmed the role of penA in generating both low-level and high-level CAZ resistance in these clinical isolates. Similar to previous studies, the amino acid substitution altered substrate specificity to other β-lactams, suggesting a potential fitness cost associated with this mutation, a finding that could be exploited to improve therapeutic outcomes in patients harboring CAZ resistant B. pseudomallei. Our study is the first to functionally characterize CAZ resistance in clinical isolates of B. pseudomallei and to provide proven and clinically relevant signatures for monitoring the development of antibiotic resistance in this important pathogen.

A recent emergence of Cryptococcus gattii in the Pacific Northwest involves strains that fall into three primarily clonal molecular subtypes: VGIIa, VGIIb and VGIIc. Multilocus sequence typing (MLST) and variable number tandem repeat analysis appear to identify little diversity within these molecular subtypes. Given the apparent expansion of these subtypes into new geographic areas and their ability to cause disease in immunocompetent individuals, differentiation of isolates belonging to these subtypes could be very important from a public health perspective. We used whole genome sequence typing (WGST) to perform fine-scale phylogenetic analysis on 20 C. gattii isolates, 18 of which are from the VGII molecular type largely responsible for the Pacific Northwest emergence. Analysis both including and excluding (289,586 SNPs and 56,845 SNPs, respectively) molecular types VGI and VGIII isolates resulted in phylogenetic reconstructions consistent, for the most part, with MLST analysis but with far greater resolution among isolates. The WGST analysis presented here resulted in identification of over 100 SNPs among eight VGIIc isolates as well as unique genotypes for each of the VGIIa, VGIIb and VGIIc isolates. Similar levels of genetic diversity were found within each of the molecular subtype isolates, despite the fact that the VGIIb clade is thought to have emerged much earlier. The analysis presented here is the first multi-genome WGST study to focus on the C. gattii molecular subtypes involved in the Pacific Northwest emergence and describes the tools that will further our understanding of this emerging pathogen.

Widespread resistance among circulating influenza A strains to at least one of the anti-influenza drugs is a major public health concern. A triple combination antiviral drug (TCAD) regimen comprised of amantadine, oseltamivir, and ribavirin has been shown to have synergistic and broad spectrum activity against influenza A strains, including drug resistant strains. Here, we used mathematical modeling along with three different experimental approaches to understand the effects of single agents, double combinations, and the TCAD regimen on resistance in influenza in vitro, including: 1) serial passage at constant drug concentrations, 2) serial passage at escalating drug concentrations, and 3) evaluation of the contribution of each component of the TCAD regimen to the suppression of resistance. Consistent with the modeling which demonstrated that three drugs were required to suppress the emergence of resistance in influenza A, treatment with the TCAD regimen resulted in the sustained suppression of drug resistant viruses, whereas treatment with amantadine alone or the amantadine-oseltamivir double combination led to the rapid selection of resistant variants which comprised ∼100% of the population. Furthermore, the TCAD regimen imposed a high genetic barrier to resistance, requiring multiple mutations in order to escape the effects of all the drugs in the regimen. Finally, we demonstrate that each drug in the TCAD regimen made a significant contribution to the suppression of virus breakthrough and resistance at clinically achievable concentrations. Taken together, these data demonstrate that the TCAD regimen was superior to double combinations and single agents at suppressing resistance, and that three drugs at a minimum were required to impede the selection of drug resistant variants in influenza A virus. The use of mathematical modeling with multiple experimental designs and molecular readouts to evaluate and optimize combination drug regimens for the suppression of resistance may be broadly applicable to other infectious diseases.

A real-time PCR assay, BurkDiff, was designed to target a unique conserved region in the B. pseudomallei and B. mallei genomes containing a SNP that differentiates the two species. Sensitivity and specificity were assessed by screening BurkDiff across 469 isolates of B. pseudomallei, 49 isolates of B. mallei, and 390 isolates of clinically relevant non-target species. Concordance of results with traditional speciation methods and no cross-reactivity to non-target species show BurkDiff is a robust, highly validated assay for the detection and differentiation of B. pseudomallei and B. mallei.