Search

Whole-genome sequencing has advanced our understanding of the population structure of the pathogenic species complex Cryptococcus gattii, which has allowed for the phylogenomic specification of previously described major molecular type groupings and novel lineages. Recently, isolates collected in Mexico in the 1960s were determined to be genetically distant from other known molecular types and were classified as VGVI. We sequenced four clinical isolates and one veterinary isolate collected in the southwestern U.S. and Argentina during 2012-2021. Phylogenomic analysis groups these genomes with those of the Mexican VGVI isolates, expanding VGVI into a clade and establishing this molecular type as a clinically important population. These findings also potentially expand the known Cryptococcus ecological range with a previously unrecognized endemic area.

The lack of preparedness for detecting the highly infectious SARS-CoV-2 pathogen, the pathogen responsible for the COVID-19 disease, has caused enormous harm to public health and the economy. It took ∼60 days for the first reverse transcription quantitative polymerase chain reaction (RT-qPCR) tests for SARS-CoV-2 infection developed by the United States Centers for Disease Control (CDC) to be made publicly available. It then took >270 days to deploy 800,000 of these tests at a time when the estimated actual testing needs required over 6 million tests per day. Testing was therefore limited to individuals with symptoms or in close contact with confirmed positive cases. Testing strategies deployed on a population scale at ‘Day Zero’ i.e., at the time of the first reported case, would be of significant value. Next Generation Sequencing (NGS) has such Day Zero capabilities with the potential for broad and large-scale testing. However, it has limited detection sensitivity for low copy numbers of pathogens which may be present. Here we demonstrate that by using CRISPR-Cas9 to remove abundant sequences that do not contribute to pathogen detection, NGS detection sensitivity of COVID-19 is comparable to RT-qPCR. In addition, we show that this assay can be used for variant strain typing, co-infection detection, and individual human host response assessment, all in a single workflow using existing open-source analysis pipelines. This NGS workflow is pathogen agnostic, and therefore has the potential to transform how both large-scale pandemic response and focused clinical infectious disease testing are pursued in the future.SIGNIFICANCE STATEMENTThe lack of preparedness for detecting infectious pathogens has had a devastating effect on the global economy and society. Thus, a ‘Day Zero’ testing strategy, that can be deployed at the first reported case and expanded to population scale, is required. Next generation sequencing enables Day Zero capabilities but is inadequate for detecting low levels of pathogen due to abundant sequences of little biological interest. By applying the CRISPR-Cas system to remove these sequences in vitro, we show sensitivity of pathogen detection equivalent to RT-qPCR. The workflow is pathogen agnostic, and enables detection of strain types, co-infections and human host response with a single workflow and open-source analysis tools. These results highlight the potential to transform future large-scale pandemic response.

The Onchocerca lupi nematode infects dogs, cats, and humans, but whether it can be spread by coyotes has been unknown. We conducted surveillance for O. lupi nematode infection in coyotes in the southwestern United States. We identified multiple coyote populations in Arizona and New Mexico as probable reservoirs for this species.

Rhizopus spp. fungi are ubiquitous in the environment and a rare but substantial cause of infection in immunosuppressed persons and surgery patients. During 2005–2017, an abnormally high number of Rhizopus infections in surgery patients, with no apparent epidemiologic links, were reported in Argentina. To determine the likelihood of a common source of the cluster, we performed whole-genome sequencing on samples collected during 2006–2014. Most isolates were separated by >60 single-nucleotide polymorphisms, and we found no evidence for recombination or nonneutral mutation accumulation; these findings do not support common source or patient-to-patient transmission. Assembled genomes of most isolates were ≈25 Mbp, and multiple isolates had substantially larger assembled genomes (43–51 Mbp), indicative of infections with strain types that underwent genome expansion. Whole-genome sequencing has become an essential tool for studying epidemiology of fungal infections. Less discriminatory techniques may miss true relationships, possibly resulting in inappropriate attribution of point source.

The full geographic range of coccidioidomycosis is unknown, although it is most likely expanding with environmental change. We report an apparently autochthonous coccidioidomycosis patient from Spokane, Washington, USA, a location to which Coccidioides spp. are not known to be endemic.

Arizona's COVID-19 and Pets Program is a prospective surveillance study being conducted to characterize how SARS-CoV-2 impacts companion animals living in households with SARS-CoV-2-positive individuals. Among the enrolled pets, we identified a SARS-CoV-2-infected cat and dog from the same household; both animals were asymptomatic but had close contact with the symptomatic and SARS-CoV-2-positive owner. Whole genome sequencing of animal and owner specimens revealed identical viral genomes of the B.1.575 lineage, suggesting zoonotic transmission of SARS-CoV-2 from human to at least one pet. This is the first report of the B.1.575 lineage in companion animals. Genetically linking SARS-CoV-2 between people and animals, and tracking changes in SARS-CoV-2 genomes is essential to detect any cross-species SARS-CoV-2 transmission that may lead to more transmissible or severe variants that can affect humans. Surveillance studies, including genomic analyses of owner and pet specimens, are needed to further our understanding of how SARS-CoV-2 impacts companion animals., Highlights • Two companion animals from the same household tested positive for SARS-CoV-2. • Sequencing of the cat, dog, and pet owner samples revealed identical viral genomes. • Phylogenetic analysis shows that the viral genomes fall into a locally circulating clade of B.1.575. • This is the first known report of the B.1.575 lineage in companion animals. • Evidence of viral transmission from a pet owner to at least one pet in the same household was identified through genomic sequencing.

BackgroundEchinocandin resistance represents a great concern, as these drugs are recommended as first-line therapy for invasive candidiasis. Echinocandin resistance is conferred by mutations in FKS genes. Nevertheless, pathways are crucial for enabling tolerance, evolution, and maintenance of resistance. Therefore, understanding the biological processes and proteins involved in the response to caspofungin may provide clues indicating new therapeutic targets.ObjectivesWe determined the resistance mechanism and assessed the proteome response to caspofungin exposure. We then evaluated the phenotypic impact of calcineurin inhibition by FK506 and cephalosporine A (CsA) on caspofungin-resistant Candida glabrata isolates.MethodsTwenty-five genes associated with caspofungin resistance were analysed by NGS, followed by studies of the quantitative proteomic response to caspofungin exposure. Then, susceptibility testing of caspofungin in presence of FK506 and CsA was performed. The effects of calcineurin inhibitor/caspofungin combinations on heat stress (40°C), oxidative stress (0.2 and 0.4 mM menadione) and on biofilm formation (polyurethane catheter) were analysed. Finally, a Galleria mellonella model using blastospores (1 × 109 cfu/mL) was developed to evaluate the impact of the combinations on larval survival.ResultsF659-del was found in the FKS2 gene of resistant strains. Proteomics data showed some up-regulated proteins are involved in cell-wall biosynthesis, response to stress and pathogenesis, some of them being members of calmodulin–calcineurin pathway. Therefore, the impact of calmodulin inhibition was explored. Calmodulin inhibition restored caspofungin susceptibility, decreased capacity to respond to stress conditions, and reduced biofilm formation and in vivo pathogenicity.ConclusionsOur findings confirm that calmodulin-calcineurin-Crz1 could provide a relevant target in life-threatening invasive candidiasis.

Background Molecular assays are important tools for pathogen detection but need to be periodically re-evaluated with the discovery of additional genetic diversity that may cause assays to exclude target taxa or include non-target taxa. A single well-developed assay can find broad application across research, clinical, and industrial settings. Pathogen prevalence within a population is estimated using such assays and accurate results are critical for formulating effective public health policies and guiding future research. A variety of assays for the detection of Staphylococcus aureus are currently available. The utility of commercial assays for research is limited, given proprietary signatures and lack of transparent validation. Results In silico testing of existing peer-reviewed assays show that most suffer from a lack of sensitivity and specificity. We found no assays that were specifically designed and validated for quantitative use. Here we present a qPCR assay, SaQuant, for the detection and quantification of S. aureus as might be collected on sampling swabs. Sensitivity and specificity of the assay was 95.6 and 99.9 %, respectively, with a limit of detection of between 3 and 5 genome equivalents and a limit of quantification of 8.27 genome equivalents. The presence of DNA from non-target species likely to be found in a swab sample, did not impact qualitative or quantitative abilities of the assay. Conclusions This assay has the potential to serve as a valuable tool for the accurate detection and quantification of S. aureus collected from human body sites in order to better understand the dynamics of prevalence and transmission in community settings.

Individuals living in congregate settings, including those in group homes, have been disproportionately impacted by COVID-19 and may be at increased risk of exposure or infection due to underlying illness. In mid-May 2020, local public health officials responded to an outbreak of COVID-19 among staff and residents associated with a multi-residential group home that provides care for adults with intellectual and developmental disabilities. Samples were collected at 16 of the homes. In four of the homes all the residents tested positive, and in the remaining 12 houses where samples were collected, all residents tested negative. Of the 152 individuals tested, 15/58 (25.9%) residents and 27/94 (28.7%) staff were positive for SARS-CoV-2, including eight hospitalizations and four deaths. Phylogenetic analysis of genomes from this outbreak in the context of genomes from Northern Arizona shows that very few mutations separate the samples from this outbreak. A potential transmission network was developed to illustrate person-place epidemiologic linkages and further demonstrates the dynamic connections between staff and residents with respect to each group home location. Epidemiologic and genomic evidence correlate, and suggest that asymptomatic infected staff likely introduced and spread COVID-19 in this setting. Implementation of public health prevention measures alongside rapid genomic analysis can help guide policy development and guide management efforts to prevent and mitigate future outbreaks.

Antimicrobial resistance (AMR) in the nosocomial pathogen, Acinetobacter baumannii, is becoming a serious public health threat. While some mechanisms of AMR have been reported, understanding novel mechanisms of resistance is critical for identifying emerging resistance. One of the first steps in identifying novel AMR mechanisms is performing genotype/phenotype association studies. However, performing genotype/phenotype association studies is complicated by the plastic nature of the A. baumannii pan-genome. In this study, we compared the antibiograms of 12 antimicrobials associated with multiple drug families for 84 A. baumannii isolates, many isolated in Arizona, USA. in silico screening of these genomes for known AMR mechanisms failed to identify clear correlations for most drugs. We then performed a genome wide association study (GWAS) looking for associations between all possible 21-mers; this approach generally failed to identify mechanisms that explained the resistance phenotype. In order to decrease the genomic noise associated with population stratification, we compared four phylogenetically-related pairs of isolates with differing susceptibility profiles. RNA-Sequencing (RNA-Seq) was performed on paired isolates and differentially expressed genes were identified. In these isolate pairs, we identified four different potential mechanisms, highlighting the difficulty of broad AMR surveillance in this species. To verify and validate differential expression, amplicon sequencing was performed. These results suggest that a diagnostic platform based on gene expression rather than genomics alone may be beneficial in certain surveillance efforts. The implementation of such advanced diagnostics coupled with increased AMR surveillance will potentially improve A. baumannii infection treatment and patient outcomes.

Pathogen genomic data are increasingly used to characterize global and local transmission patterns of important human pathogens and to inform public health interventions. Yet, there is no current consensus on how to measure genomic variation. To test the effect of the variant-identification approach on transmission inferences for Mycobacterium tuberculosis, we conducted an experiment in which five genomic epidemiology groups applied variant-identification pipelines to the same outbreak sequence data. We compared the variants identified by each group in addition to transmission and phylogenetic inferences made with each variant set. To measure the performance of commonly used variant-identification tools, we simulated an outbreak. We compared the performance of three mapping algorithms, five variant callers and two variant filters in recovering true outbreak variants. Finally, we investigated the effect of applying increasingly stringent filters on transmission inferences and phylogenies. We found that variant-calling approaches used by different groups do not recover consistent sets of variants, which can lead to conflicting transmission inferences. Further, performance in recovering true variation varied widely across approaches. While no single variant-identification approach outperforms others in both recovering true genome-wide and outbreak-level variation, variant-identification algorithms calibrated upon real sequence data or that incorporate local reassembly outperform others in recovering true pairwise differences between isolates. The choice of variant filters contributed to extensive differences across pipelines, and applying increasingly stringent filters rapidly eroded the accuracy of transmission inferences and quality of phylogenies reconstructed from outbreak variation. Commonly used approaches to identify M. tuberculosis genomic variation have variable performance, particularly when predicting potential transmission links from pairwise genetic distances. Phylogenetic reconstruction may be improved by less stringent variant filtering. Approaches that improve variant identification in repetitive, hypervariable regions, such as long-read assemblies, may improve transmission inference.

In December of 2019, a novel coronavirus, SARS-CoV-2, emerged in the city of Wuhan, China causing severe morbidity and mortality. Since then, the virus has swept across the globe causing millions of confirmed infections and hundreds of thousands of deaths. To better understand the nature of the pandemic and the introduction and spread of the virus in Arizona, we sequenced viral genomes from clinical samples tested at the TGen North Clinical Laboratory, provided to us by the Arizona Department of Health Services, and at Arizona State University and the University of Arizona, collected as part of community surveillance projects. Phylogenetic analysis of 79 genomes we generated from across Arizona revealed a minimum of 9 distinct introductions throughout February and March. We show that >80% of our sequences descend from clades that were initially circulating widely in Europe but have since dominated the outbreak in the United States. In addition, we show that the first reported case of community transmission in Arizona descended from the Washington state outbreak that was discovered in late February. Notably, none of the observed transmission clusters are epidemiologically linked to the original travel-related cases in the state, suggesting successful early isolation and quarantine. Finally, we use molecular clock analyses to demonstrate a lack of identifiable, widespread cryptic transmission in Arizona prior to the middle of February 2020.

Altered intestinal microbiota is associated with systemic and intestinal diseases, such as inflammatory bowel disease (IBD). Dysbiotic microbiota with enhanced proinflammatory capacity is characterized by depletion of anaerobic commensals, increased proportion of facultatively anaerobic bacteria, as well as reduced diversity and stability. In this study, we developed a high-throughput in vitro screening assay to isolate intestinal commensal bacteria with anti-inflammatory capacity from a healthy fecal microbiota transplantation donor. Freshly isolated gut bacteria were screened for their capacity to attenuate Escherichia coli lipopolysaccharide (LPS)-induced interleukin 8 (IL-8) release from HT-29 cells. The screen yielded a number of Bacteroides and Parabacteroides isolates, which were identified as P. distasonis, B. caccae, B. intestinalis, B. uniformis, B. fragilis, B. vulgatus and B. ovatus using whole genome sequencing. We observed that a cell-cell contact with the epithelium was not necessary to alleviate in vitro inflammation as spent culture media from the isolates were also effective and the anti-inflammatory action did not correlate with the enterocyte adherence capacity of the isolates. The anti-inflammatory isolates also exerted enterocyte monolayer reinforcing action and lacked essential genes to synthetize hexa-acylated, proinflammatory lipid A, part of LPS. Yet, the anti-inflammatory effector molecules remain to be identified. The Bacteroides strains isolated and characterized in this study have potential to be used as so-called next-generation probiotics.

Water reclamation provides a valuable resource for meeting nonpotable water demands. However, little is known about the potential for wastewater reuse to disseminate antibiotic resistance genes (ARGs). Here, samples were collected seasonally in 2014–2015 from four U.S. utilities’ reclaimed and potable water distribution systems before treatment, after treatment, and at five points of use (POU). Shotgun metagenomic sequencing was used to profile the resistome (i.e., full contingent of ARGs) of a subset (n = 38) of samples. Four ARGs (qnrA, blaTEM, vanA, sul1) were quantified by quantitative polymerase chain reaction. Bacterial community composition (via 16S rRNA gene amplicon sequencing), horizontal gene transfer (via quantification of intI1 integrase and plasmid genes), and selection pressure (via detection of metals and antibiotics) were investigated as potential factors governing the presence of ARGs. Certain ARGs were elevated in all (sul1; p ≤ 0.0011) or some (blaTEM, qnrA; p ≤ 0.0145) reclaimed POU s...

The 2020 pandemic has transformed the world and elicited thousands of studies to better understand the SARS-CoV-2 virus. Viral load has been a common measure to monitor treatment therapies and associate viral dynamics with patient outcomes; however, methods associated with viral load have varied across studies. These variations have the potential to sacrifice the accuracy of findings as they often do not account for inter-assay variation or variation across samples. In a retrospective study of nasopharyngeal samples, we found a significant amount of variation within the DNA and RNA targets; for example, across time within a single patient, there was an average of a 32-fold change. Further, we explore the impacts of host normalization on 94 clinical samples using the TGen Quantitative SARS-CoV-2 assay, finding that without host normalization samples with the same viral concentration can have up to 100-fold variation in the viral load.

In December of 2019, a novel coronavirus, SARS-CoV-2, emerged in the city of Wuhan, China, causing severe morbidity and mortality. Since then, the virus has swept across the globe, causing millions of confirmed infections and hundreds of thousands of deaths. To better understand the nature of the pandemic and the introduction and spread of the virus in Arizona, we sequenced viral genomes from clinical samples tested at the TGen North Clinical Laboratory, the Arizona Department of Health Services, and those collected as part of community surveillance projects at Arizona State University and the University of Arizona. Phylogenetic analysis of 84 genomes from across Arizona revealed a minimum of 11 distinct introductions inferred to have occurred during February and March. We show that >80% of our sequences descend from strains that were initially circulating widely in Europe but have since dominated the outbreak in the United States. In addition, we show that the first reported case of community transmission in Arizona descended from the Washington state outbreak that was discovered in late February. Notably, none of the observed transmission clusters are epidemiologically linked to the original travel-related case in the state, suggesting successful early isolation and quarantine. Finally, we use molecular clock analyses to demonstrate a lack of identifiable, widespread cryptic transmission in Arizona prior to the middle of February 2020., mBio, 11 (5), ISSN:2150-7511, ISSN:2161-2129

Pathogen genomic data are increasingly used to characterize global and local transmission patterns of important human pathogens and to inform public health interventions. Yet there is no current consensus on how to measure genomic variation. We investigated the effects of variant identification approaches on transmission inferences for M. tuberculosis by comparing variants identified by five different groups in the same sequence data from a clonal outbreak. We then measured the performance of commonly used variant calling approaches in recovering variation in a simulated tuberculosis outbreak and tested the effect of applying increasingly stringent filters on transmission inferences and phylogenies. We found that variant calling approaches used by different groups do not recover consistent sets of variants, often leading to conflicting transmission inferences. Further, performance in recovering true outbreak variation varied widely across approaches. Finally, stringent filters rapidly eroded the accuracy of transmission inferences and quality of phylogenies reconstructed from outbreak variation. We conclude that measurements of genetic distance and phylogenetic structure are dependent on variant calling approach. Variant calling algorithms trained upon true sequence data outperform other approaches and enable inclusion of repetitive regions typically excluded from genomic epidemiology studies, maximizing the information gleaned from outbreak genomes.

There has been a considerable upsurge of extensive, treatment recalcitrant, dermatophytosis presenting as tinea corporis and tinea cruris in India since the past few years. Genome analysis of Trichophyton species causing severe superficial dermatophytosis in North India confirmed a unique clade related to the T. mentagrophytes/interdigitale complex, seeming to belong to an early diverging clade of the complex. The Indian Trichophyton species genomes were highly related showing only up to 42 SNPs between any two isolates confirming their clonal origin. Other genetic approaches such as ITS sequencing and multigene phylogeny used in this study were contradictory or inconclusive to show the differentiation of these isolates from T. mentagrophytes/T. interdigitale. Remarkably, high rates of resistance to all three commonly used oral antifungals, i.e., 36% for terbinafine (MICs 4 to ≥32 mg/L), 39.5% for fluconazole (MIC range 32 to ≥64 mg/L) and griseofulvin (Geometric mean MIC ≥ 4 mg/L) were observed. Two important amino acid substitutions (Leu393Phe or Phe397Leu) leading to a terbinafine resistant phenotype were found in the squalene epoxidase protein of all tested terbinafine resistant isolates. All 20 examined genomes presented a high mobility group (HMG) domain transcription factor gene corresponding to mating type (+). Of these, three isolates also showed positivity for both alpha-box and HMG in the genome which might indicate hybridization or an incomplete sexual cycle. Therefore, we highlight the potential of this organism to rapidly spread alleles that might be driving antifungal resistance among its population. This new population of Trichophyton with high rates of in vitro antifungal resistance seems to be driving an ongoing outbreak of dermatophytosis in India. Our study highlights difficulties in identifying isolates from the Trichophyton mentagrophytes/interdigitale clade of the genus using currently available molecular tools. High resistance rates of terbinafine warrant further clinical studies to assess its utility in the treatment of dermatophytosis caused by this clonal strain.

Coccidioidomycosis is a debilitating fungal disease caused by inhalation of arthroconidia. We developed a novel approach for detection of airborne Coccidioides and used it to investigate the distribution of arthroconidia across the Phoenix, Arizona, metropolitan area. Air filters were collected daily from 21 stationary air-sampling units across the area: the first set collected before, during and after a large dust storm on August 25, 2015, and the second over the 45-day period September 25–November 8, 2016. Analysis of DNA extracted from the filters demonstrated that the day of the dust storm was not associated with increase of Coccidioides in air samples, although evidence of the low-level polymerase chain reaction (PCR) inhibition was observed in DNA extracted from samples collected on the day of the dust storm. Testing over 45 days identified uneven geographic distribution suggesting Coccidioides hot spots. In 2016, highest daily concentration of arthroconidia was observed between September 25–October 20, and only sporadic low levels were detected after that. These results provide evidence of seasonality and uneven spatial distribution of Coccidioides in the air. Our results demonstrate that routine air monitoring for arthroconidia is possible and provides an important tool for Coccidioides surveillance, which can address important questions about environmental exposure and human infection.

Cryptococcus gattii molecular type VGII is one of the causes of cryptococcosis, a severe fungal disease that is acquired from the environment and affects many host species (including humans and koalas). In Australia, disease caused by C. gattii VGII is largely confined to western and central northern parts of the country, with sporadic cases reported in eastern Australia. We investigated an unusual case cluster of cryptococcosis, caused predominantly by C. gattii VGII, in a group of captive koalas in eastern Australia. This research identified that the movements of koalas between wildlife parks, including an initial transfer of a koala from Western Australia, introduced and subsequently spread C. gattii VGII in this captive environment. The spread of this pathogen by koalas could also impact other species, and these findings are significant in the implications they have for the management of koala transfers and captive environments., Cryptococcus gattii molecular type VGII is one of the etiologic agents of cryptococcosis, a systemic mycosis affecting a wide range of host species. Koalas (Phascolarctos cinereus) exhibit a comparatively high prevalence of cryptococcosis (clinical and subclinical) and nasal colonization, particularly in captivity. In Australia, disease associated with C. gattii VGII is typically confined to Western Australia and the Northern Territory (with sporadic cases reported in eastern Australia), occupying an enigmatic ecologic niche. A cluster of cryptococcosis in captive koalas in eastern Australia (five confirmed cases, a further two suspected), caused predominantly by C. gattii VGII, was investigated by surveying for subclinical disease, culturing koala nasal swabs and environmental samples, and genotyping cryptococcal isolates. URA5 restriction fragment length polymorphism analysis, multilocus sequence typing (MLST), and whole-genome sequencing (WGS) provided supportive evidence that the transfer of koalas from Western Australia and subsequently between several facilities in Queensland spread VGII into uncontaminated environments and environments in which C. gattii VGI was endemic. MLST identified VGII isolates as predominantly sequence type 7, while WGS further confirmed a limited genomic diversity and revealed a basal relationship with isolates from Western Australia. We hypothesize that this represents a founder effect following the introduction of a koala from Western Australia. Our findings suggest a possible competitive advantage for C. gattii VGII over VGI in the context of this captive koala environment. The ability of koalas to seed C. gattii VGII into new environments has implications for the management of captive populations and movements of koalas between zoos. IMPORTANCE Cryptococcus gattii molecular type VGII is one of the causes of cryptococcosis, a severe fungal disease that is acquired from the environment and affects many host species (including humans and koalas). In Australia, disease caused by C. gattii VGII is largely confined to western and central northern parts of the country, with sporadic cases reported in eastern Australia. We investigated an unusual case cluster of cryptococcosis, caused predominantly by C. gattii VGII, in a group of captive koalas in eastern Australia. This research identified that the movements of koalas between wildlife parks, including an initial transfer of a koala from Western Australia, introduced and subsequently spread C. gattii VGII in this captive environment. The spread of this pathogen by koalas could also impact other species, and these findings are significant in the implications they have for the management of koala transfers and captive environments.

Enteroviruses frequently result in respiratory and gastrointestinal illness; however, multiple subtypes, including poliovirus, can cause severe neurologic disease. Recent biennial increases (i.e., 2014, 2016, and 2018) in cases of non-polio acute flaccid paralysis have led to speculations that other enteroviruses, specifically enterovirus D68 (EV-D68), are emerging to fill the niche that was left from poliovirus eradication. A cluster of 11 suspect cases of pediatric acute flaccid myelitis (AFM) was identified in 2016 in Phoenix, AZ. Multiple genomic analyses identified the presence of EV-D68 in the majority of clinical AFM cases. Beyond limited detection of herpesvirus, no other likely etiologies were found in the cluster. These findings strengthen the likelihood that EV-D68 is a cause of AFM and show that the rapid molecular assays developed for this study are useful for investigations of AFM and EV-D68., Enteroviruses are a common cause of respiratory and gastrointestinal illness, and multiple subtypes, including poliovirus, can cause neurologic disease. In recent years, enterovirus D68 (EV-D68) has been associated with serious neurologic illnesses, including acute flaccid myelitis (AFM), frequently preceded by respiratory disease. A cluster of 11 suspect cases of pediatric AFM was identified in September 2016 in Phoenix, AZ. To determine if these cases were associated with EV-D68, we performed multiple genomic analyses of nasopharyngeal (NP) swabs and cerebrospinal fluid (CSF) material from the patients, including real-time PCR and amplicon sequencing targeting the EV-D68 VP1 gene and unbiased microbiome and metagenomic sequencing. Four of the 11 patients were classified as confirmed cases of AFM, and an additional case was classified as probable AFM. Real-time PCR and amplicon sequencing detected EV-D68 virus RNA in the three AFM patients from which NP swabs were collected, as well as in a fourth patient diagnosed with acute disseminated encephalomyelitis, a disease that commonly follows bacterial or viral infections, including enterovirus. No other obvious etiological causes for AFM were identified by 16S or RNA and DNA metagenomic sequencing in these cases, strengthening the likelihood that EV-D68 is an etiological factor. Herpes simplex viral DNA was detected in the CSF of the fourth case of AFM and in one additional suspect case from the cluster. Multiple genomic techniques, such as those described here, can be used to diagnose patients with suspected EV-D68 respiratory illness, to aid in AFM diagnosis, and for future EV-D68 surveillance and epidemiology.

Health care-acquired infections (HAIs) kill tens of thousands of people each year and add significantly to health care costs. Multidrug-resistant and epidemic strains are a large proportion of HAI agents, and multidrug-resistant strains of Klebsiella pneumoniae , a leading HAI agent, have caused an urgent public health crisis. In the health care environment, patient colonization by K. pneumoniae precedes infection, and transmission via colonization leads to outbreaks. Periodic patient screening for K. pneumoniae colonization has the potential to curb the number of HAIs. In this report, we describe the design and validation of KlebSeq, a highly informative screening tool that detects Klebsiella species and identifies clinically important strains and characteristics by using highly multiplexed amplicon sequencing without a live-culturing step. We demonstrate the utility of this tool on several complex specimen types, including urine, wound swabs and tissue, and several types of respiratory and fecal specimens, showing K. pneumoniae species and clonal group identification and antimicrobial resistance and virulence profiling, including capsule typing. Use of this amplicon sequencing tool to screen patients for Klebsiella carriage could inform health care staff of the risk of infection and outbreak potential. KlebSeq also serves as a model for next-generation molecular tools for public health and health care, as expansion of this tool can be used for several other HAI agents or applications.

The hyper-virulent emm59 genotype of invasive group A Streptococcus was identified in northern Arizona in 2015. Eighteen isolates belonging to a genomic cluster grouped most closely with recently identified isolates in New Mexico. The continued transmission of emm59 in the southwestern United States poses a public health concern.

Genomic epidemiology has proven successful for real-time and retrospective monitoring of small and large-scale outbreaks. Here, we report two genomic sequencing and analysis strategies for rapid-turnaround or high-throughput processing of metagenomic samples. The rapid-turnaround method was designed to provide a quick phylogenetic snapshot of samples at the heart of active outbreaks, and has a total turnaround time of st Century.

Need for global water security has spurred growing interest in wastewater reuse to offset demand for municipal water. While reclaimed (i.e., nonpotable) microbial water quality regulations target fecal indicator bacteria, opportunistic pathogens (OPs), which are subject to regrowth in distribution systems and spread via aerosol inhalation and other noningestion routes, may be more relevant. This study compares the occurrences of five OP gene markers ( Acanthamoeba spp., Legionella spp., Mycobacterium spp., Naegleria fowleri, Pseudomonas aeruginosa) in reclaimed versus potable water distribution systems and characterizes factors potentially contributing to their regrowth. Samples were collected over four sampling events at the point of compliance for water exiting treatment plants and at five points of use at four U.S. utilities bearing both reclaimed and potable water distribution systems. Reclaimed water systems harbored unique water chemistry (e.g., elevated nutrients), microbial community composition, and OP occurrence patterns compared to potable systems examined here and reported in the literature. Legionella spp. genes, Mycobacterium spp. genes, and total bacteria, represented by 16S rRNA genes, were more abundant in reclaimed than potable water distribution system samples ( p ≤ 0.0001). This work suggests that further consideration should be given to managing reclaimed water distribution systems with respect to nonpotable exposures to OPs.

Staphylococcus aureus is a major human pathogen worldwide in both community and health care settings. Surveillance for S. aureus strains is important to our understanding of their spread and to informing infection prevention and control. Confusion surrounding the strain nomenclature of one of the most prevalent lineages of S. aureus, clonal complex 8 (CC8), and the imprecision of current tools for typing S. aureus make surveillance and source tracing difficult and sometimes misleading. In this study, we clarify the CC8 strain designations and propose a new typing scheme for CC8 isolates that is rapid and easy to use. This typing scheme is based on relatively stable genomic markers, and we demonstrate its superiority over traditional typing techniques. This scheme has the potential to greatly improve epidemiological investigations of S. aureus., Strains of Staphylococcus aureus in clonal complex 8 (CC8), including USA300, USA500, and the Iberian clone, are prevalent pathogens in the United States, both inside and outside health care settings. Methods for typing CC8 strains are becoming obsolete as the strains evolve and diversify, and whole-genome sequencing has shown that some strain types fall into multiple sublineages within CC8. In this study, we attempt to clarify the strain nomenclature of CC8, classifying the major strain types based on whole-genome sequence phylogenetics using both methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) genomes. We show that isolates of the Archaic and Iberian clones from decades ago make up the most basal clade of the main CC8 lineages and that at least one successful lineage of CC8, made up mostly of MSSA, diverged before the other well-known strain types USA500 and USA300. We also show that the USA500 type includes two clades separated by the previously described “Canadian epidemic MRSA” strain CMRSA9, that one clade containing USA500 also contains the USA300 clade, and that the USA300-0114 strain type is not a monophyletic group. Additionally, we present a rapid, simple CC8 strain-typing scheme using real-time PCR assays that target single nucleotide polymorphisms (SNPs) derived from our CC8 phylogeny and show the significant benefit of using more stable genomic markers based on evolutionary lineages over traditional S. aureus typing techniques. This more accurate and accessible S. aureus typing system may improve surveillance and better inform the epidemiology of this very important pathogen. IMPORTANCE Staphylococcus aureus is a major human pathogen worldwide in both community and health care settings. Surveillance for S. aureus strains is important to our understanding of their spread and to informing infection prevention and control. Confusion surrounding the strain nomenclature of one of the most prevalent lineages of S. aureus, clonal complex 8 (CC8), and the imprecision of current tools for typing S. aureus make surveillance and source tracing difficult and sometimes misleading. In this study, we clarify the CC8 strain designations and propose a new typing scheme for CC8 isolates that is rapid and easy to use. This typing scheme is based on relatively stable genomic markers, and we demonstrate its superiority over traditional typing techniques. This scheme has the potential to greatly improve epidemiological investigations of S. aureus.

Coccidioides immitis and Coccidioides posadasii are soil fungi endemic to desert regions of the southwestern United States, and the causative agents of valley fever, or coccidioidomycosis. Studies have shown that the distribution of Coccidioides in soils is sporadic and cannot be explained by soil characteristics alone, suggesting that biotic and other abiotic factors should be examined. However, tools to reliably and robustly screen the large number of soils needed to investigate these potential associations have not been available. Thus, we developed a real-time polymerase chain reaction (PCR) assay for testing environmental samples by modifying CocciDx, an assay validated for testing clinical specimens to facilitate coccidioidomycosis diagnosis. For this study, we collected soil samples from previously established locations of C. posadasii in Arizona and new locations in fall 2013 and spring 2014, and screened the extracted DNA with the new assay known as CocciEnv. To verify the presence of Coccidioides in soil using an alternate method, we employed next generation amplicon sequencing targeting the ITS2 region. Results show our modified assay, CocciEnv, is a rapid and robust method for detecting Coccidioides DNA in complex environmental samples. The ability to test a large number of soils for the presence of Coccidioides is a much-needed tool in the understanding of the ecology of the organism and epidemiology of the disease and will greatly improve our understanding of this human pathogen.

Background In a patient with positive serum serology for coccidioidomycosis, the differential diagnosis of concurrent pleural effusions can be challenging. We, therefore, sought to clarify the performance characteristics of biochemical, serologic, and nucleic-acid-based testing in an attempt to avoid invasive procedures. The utility of adenosine deaminase (ADA), coccidioidal serology, and polymerase chain reaction (PCR) in the evaluation of pleuropulmonary coccidioidomycosis has not been previously reported. Methods Forty consecutive patients evaluated for pleuropulmonary coccidioidomycosis were included. Demographic data, pleural fluid values, culture results, and clinical diagnoses were obtained from patient chart review. ADA testing was performed by ARUP Laboratories, coccidioidal serologic testing was performed by the University of California-Davis coccidioidomycosis serology laboratory, and PCR testing was performed by the Translational Genomics Research Institute using a previously published methodology. Results Fifteen patients were diagnosed with pleuropulmonary coccidioidomycosis by European Organization for the Research and Treatment of Cancer/Mycoses Study Group criteria. Pleural fluid ADA concentrations were Conclusion Contrary to prior speculation, ADA levels in pleuropulmonary coccidioidomycosis were not elevated in this study. The sensitivity and specificity of coccidioidal serologic testing in nonserum samples remained high, but the clinical usefulness of PCR testing in pleural fluid was disappointing and was comparable to pleural fluid culture.

BackgroundPrevention of nosocomial transmission of infections is a central responsibility in the healthcare environment, and accurate identification of transmission events presents the first challenge. Phylogenetic analysis based on whole genome sequencing provides a high-resolution approach for accurately relating isolates to one another, allowing precise identification or exclusion of transmission events and sources for nearly all cases. We sequenced 24 methicillin-resistant Staphylococcus aureus (MRSA) genomes to retrospectively investigate a suspected point source of three surgical site infections (SSIs) that occurred over a one-year period. The source of transmission was believed to be a surgical team member colonized with MRSA, involved in all surgeries preceding the SSI cases, who was subsequently decolonized. Genetic relatedness among isolates was determined using whole genome single nucleotide polymorphism (SNP) data.ResultsWhole genome SNP typing (WGST) revealed 283 informative SNPs between the surgical team member’s isolate and the closest SSI isolate. The second isolate was 286 and the third was thousands of SNPs different, indicating the nasal carriage strain from the surgical team member was not the source of the SSIs. Given the mutation rates estimated for S. aureus, none of the SSI isolates share a common ancestor within the past 14 years, further discounting any common point source for these infections. The decolonization procedures and resources spent on the point source infection control could have been prevented if WGST was performed at the time of the suspected transmission, instead of retrospectively.ConclusionsWhole genome sequence analysis is an ideal method to exclude isolates involved in transmission events and nosocomial outbreaks, and coupling this method with epidemiological data can determine if a transmission event occurred. These methods promise to direct infection control resources more appropriately.

Healthcare-acquired infections (HAIs) kill tens of thousands of people each year and add significantly to healthcare costs. Multidrug resistant and epidemic strains are a large proportion of HAI agents, and multidrug resistant strains ofKlebsiella pneumoniae, a leading HAI agent, have become an urgent public health crisis. In the healthcare environment, patient colonization ofK. pneumoniaeprecedes infection, and transmission via colonization leads to outbreaks. Periodic patient screening forK. pneumoniaecolonization has cost-effective and life-saving potential. In this study, we describe the design and validation of KlebSeq, a highly informative screening tool that detectsKlebsiellaspecies and identifies clinically important strains and characteristics using highly multiplexed amplicon sequencing without a live culturing step. We demonstrate the utility of this tool on several complex specimen types including urine, wound swabs and tissue, several types of respiratory, and fecal, showingK. pneumoniaespecies and clonal group identification and antimicrobial resistance and virulence profiling, including capsule typing. Use of this amplicon sequencing tool can be used to screen patients forK. pneumoniaecarriage to assess risk of infection and outbreak potential, and the expansion of this tool can be used for several other HAI agents or applications.

Melioidosis is an emerging infectious disease caused by the soil bacterium Burkholderia pseudomallei . In diagnostic and forensic settings, molecular detection assays need not only high sensitivity with low limits of detection but also high specificity. In a direct comparison of published and newly developed TaqMan PCR assays, we found the TTS1- orf2 assay to be superior in detecting B. pseudomallei directly from clinical specimens. The YLF/BTFC multiplex assay (targeting the Yersinia -like fimbrial/ Burkholderia thailandensis -like flagellum and chemotaxis region) also showed high diagnostic sensitivity and provides additional information on possible geographic origin.

Francisella tularensis contains several highly pathogenic subspecies, including Francisella tularensis subsp. holarctica , whose distribution is circumpolar in the northern hemisphere. The phylogeography of these subspecies and their subclades was examined using whole-genome single nucleotide polymorphism (SNP) analysis, high-density microarray SNP genotyping, and real-time-PCR-based canonical SNP (canSNP) assays. Almost 30,000 SNPs were identified among 13 whole genomes for phylogenetic analysis. We selected 1,655 SNPs to genotype 95 isolates on a high-density microarray platform. Finally, 23 clade- and subclade-specific canSNPs were identified and used to genotype 496 isolates to establish global geographic genetic patterns. We confirm previous findings concerning the four subspecies and two Francisella tularensis subsp. tularensis subpopulations and identify additional structure within these groups. We identify 11 subclades within F. tularensis subsp. holarctica , including a new, genetically distinct subclade that appears intermediate between Japanese F. tularensis subsp. holarctica isolates and the common F. tularensis subsp. holarctica isolates associated with the radiation event (the B radiation) wherein this subspecies spread throughout the northern hemisphere. Phylogenetic analyses suggest a North American origin for this B-radiation clade and multiple dispersal events between North America and Eurasia. These findings indicate a complex transmission history for F. tularensis subsp. holarctica .

Optimising DNA extraction from clinical samples for Burkholderia pseudomallei Type III secretion system real-time PCR in suspected melioidosis patients confirmed that urine and sputum are useful diagnostic samples. Direct testing on blood remains problematic; testing DNA extracted from plasma was superior to DNA from whole blood or buffy coat.

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 bla KPC 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.

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.

The Staphylococcus aureus clonal lineage CC45 is a predominant colonizer of healthy individuals in northern Europe and constitutes a highly basal cluster of the S. aureus population. Here, we report the complete genome sequence of S. aureus strain CA-347 (NRS648), a representative of the methicillin-resistant USA600 clone predominantly found in the United States.

In a patient with positive serum serology for coccidioidomycosis, the differential diagnosis of concurrent pleural effusions can be challenging. We, therefore, sought to clarify the performance characteristics of biochemical, serologic, and nucleic-acid-based testing in an attempt to avoid invasive procedures. The utility of adenosine deaminase (ADA), coccidioidal serology, and polymerase chain reaction (PCR) in the evaluation of pleuropulmonary coccidioidomycosis has not been previously reported.Forty consecutive patients evaluated for pleuropulmonary coccidioidomycosis were included. Demographic data, pleural fluid values, culture results, and clinical diagnoses were obtained from patient chart review. ADA testing was performed by ARUP Laboratories, coccidioidal serologic testing was performed by the University of California-Davis coccidioidomycosis serology laboratory, and PCR testing was performed by the Translational Genomics Research Institute using a previously published methodology.Fifteen patients were diagnosed with pleuropulmonary coccidioidomycosis by European Organization for the Research and Treatment of Cancer/Mycoses Study Group criteria. Pleural fluid ADA concentrations were40 IU/L in all patients (range,1.0-28.6 IU/L; median, 4.7). The sensitivity and specificity of coccidioidal serologic testing was 100% in this study. The specificity of PCR testing was high (100%), although the overall sensitivity remained low, and was comparable to the experience of others in the clinical use of PCR for coccidioidal diagnostics.Contrary to prior speculation, ADA levels in pleuropulmonary coccidioidomycosis were not elevated in this study. The sensitivity and specificity of coccidioidal serologic testing in nonserum samples remained high, but the clinical usefulness of PCR testing in pleural fluid was disappointing and was comparable to pleural fluid culture.

Melioidosis is caused by Burkholderia pseudomallei, a Gram-negative bacillus, primarily found in soils in Southeast Asia and northern Australia. A recent case of melioidosis in non-endemic Arizona was determined to be the result of locally acquired infection, as the patient had no travel history to endemic regions and no previous history of disease. Diagnosis of the case was confirmed through multiple microbiologic and molecular techniques. To enhance the epidemiological analysis, we conducted several molecular genotyping procedures, including multi-locus sequence typing, SNP-profiling, and whole genome sequence typing. Each technique has different molecular epidemiologic advantages, all of which provided evidence that the infecting strain was most similar to those found in Southeast Asia, possibly originating in, or around, Malaysia. Advancements in new typing technologies provide genotyping resolution not previously available to public health investigators, allowing for more accurate source identification., Author Summary Melioidosis is a bacterial disease caused by percutaneous inoculation, aspiration or ingestion of the soil bacteria Burkholderia pseudomallei. Melioidosis is primarily found in Southeast Asia and Northern Australia, and, to a lesser degree, nearby regions. A recent case of melioidosis in Southwestern United States (Southern Arizona) prompted a detailed epidemiological and molecular investigation to discover the source of infection. The authors describe the use of multiple genomic analysis tools to aid in this investigation. The results of these analyses uniformly identified Southeast Asia as the source of the strain that infected the patient, however the epidemiological investigation had determined the patient had no international travel or known exposures to Southeast Asian products. New cutting edge technologies, such as next generation sequencing, are quickly being adapted into epidemiologic investigations, particularly for cases and outbreaks of unknown origin, although older, mature technologies with larger existing databases will still be needed for appropriate comparative analyses.

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.