The role of plasmids in carbapenem resistant E. coli in Alameda County, California.

BACKGROUND: Antimicrobial resistant infections continue to be a leading global public health crisis. Mobile genetic elements, such as plasmids, have been shown to play a major role in the dissemination of antimicrobial resistance (AMR) genes. Despite its ongoing threat to human health, surveillance of AMR in the United States is often limited to phenotypic resistance. Genomic analyses are important to better understand the underlying resistance mechanisms, assess risk, and implement appropriate prevention strategies. This study aimed to investigate the extent of plasmid mediated antimicrobial resistance that can be inferred from short read sequences of carbapenem resistant E. coli (CR-Ec) in Alameda County, California. E. coli isolates from healthcare locations in Alameda County were sequenced using an Illumina MiSeq and assembled with Unicycler. Genomes were categorized according to predefined multilocus sequence typing (MLST) and core genome multilocus sequence typing (cgMLST) schemes. Resistance genes were identified and corresponding contigs were predicted to be plasmid-borne or chromosome-borne using two bioinformatic tools (MOB-suite and mlplasmids). RESULTS: Among 82 of CR-Ec identified between 2017 and 2019, twenty-five sequence types (STs) were detected. ST131 was the most prominent (n = 17) followed closely by ST405 (n = 12). blaCTX-M were the most common ESBL genes and just over half (18/30) of these genes were predicted to be plasmid-borne by both MOB-suite and mlplasmids. Three genetically related groups of E. coli isolates were identified with cgMLST. One of the groups contained an isolate with a chromosome-borne blaCTX-M-15 gene and an isolate with a plasmid-borne blaCTX-M-15 gene. CONCLUSIONS: This study provides insights into the dominant clonal groups driving carbapenem resistant E. coli infections in Alameda County, CA, USA clinical sites and highlights the relevance of whole-genome sequencing in routine local genomic surveillance. The finding of multi-drug resistant plasmids harboring high-risk resistance genes is of concern as it indicates a risk of dissemination to previously susceptible clonal groups, potentially complicating clinical and public health intervention.

Imipenem-Relebactam Susceptibility and Genotypic Characteristics of Carbapenem-Resistant Enterobacterales Identified during Population-Based Surveillance.

Laboratories submit all carbapenem-resistant Enterobacter, Escherichia coli, and Klebsiella species to the Alameda County Public Health Department (ACPHD). ACPHD evaluated 75 isolates submitted during 9 months for susceptibility to imipenem-relebactam (I-R) and, using whole-genome sequencing, identified ß-lactamase genes. Of 60 (80%) isolates susceptible to I-R, 8 (13%) had detectable carbapenemase genes, including 4 KPC, two NDM, and two OXA-48-like; we described the relationship between the presence of ß-lactamase resistance genes and susceptibility to I-R.

Antibody detection by agglutination-PCR (ADAP) enables early diagnosis of HIV infection by oral fluid analysis.

Oral fluid (OF) is a highly effective substrate for population-based HIV screening efforts, as it is noninfectious and significantly easier to collect than blood. However, anti-HIV antibodies are found at far lower concentrations in OF compared with blood, leading to poor sensitivity and a longer period of time from infection to detection threshold. Thus, despite its inherent advantages in sample collection, OF is not widely used for population screening. Here we report the development of an HIV OF assay based on Antibody Detection by Agglutination-PCR (ADAP) technology. This assay is 1,000-10,000 times more analytically sensitive than clinical enzyme-linked immunoassays (EIAs), displaying both 100% clinical sensitivity and 100% specificity for detecting HIV antibodies within OF samples. We show that the enhanced analytical sensitivity enables this assay to correctly identify HIV-infected individuals otherwise missed by current OF assays. We envision that the attributes of this improved HIV OF assay can increase testing rates of at-risk individuals while enabling diagnosis and treatment at an earlier time point.

Simultaneous Evaluation of Diagnostic Assays for Pharyngeal and Rectal Neisseria gonorrhoeae and Chlamydia trachomatis Using a Master Protocol.

BACKGROUND: Pharyngeal and rectal Neisseria gonorrhoeae and Chlamydia trachomatis play important roles in infection and antibacterial resistance transmission, but no US Food and Drug Administration (FDA)-cleared assays for detection at these sites existed prior to this study. The objective was to estimate performance of assays to detect those infections in pharyngeal and rectal specimens to support regulatory submission. METHODS: We performed a cross-sectional, single-visit study of adults seeking sexually transmitted infection testing at 9 clinics in 7 states. We collected pharyngeal and rectal swabs from participants. The primary outcome was positive and negative percent agreement for detection of N. gonorrhoeae and C. trachomatis for 3 investigational assays compared to a composite reference. Secondary outcomes included positivity as well as positive and negative predictive values and likelihood ratios. Subgroup analyses included outcomes by symptom status and sex. RESULTS: A total of 2598 participants (79% male) underwent testing. We observed N. gonorrhoeae positivity of 8.1% in the pharynx and 7.9% in the rectum and C. trachomatis positivity of 2.0% in the pharynx and 8.7% in the rectum. Positive percent agreement ranged from 84.8% to 96.5% for different anatomic site infection combinations, whereas negative percent agreement was 98.8% to 99.6%. CONCLUSIONS: This study utilized a Master Protocol to generate diagnostic performance data for multiple assays from different manufacturers in a single study population, which ultimately supported first-in-class FDA clearance for extragenital assays. We observed very good positive percent agreement when compared to a composite reference method for the detection of both pharyngeal and rectal N. gonorrhoeae and C. trachomatis. CLINICAL TRIALS REGISTRATION: NCT02870101.

Assessing the Use of PCR To Screen for Shedding of Salmonella enterica in Infected Humans.

Recovery from enteric bacterial illness often includes a phase of organismal shedding over a period of days to months. The monitoring of this process through laboratory testing forms the foundation of public health action to prevent further transmission. Regulations in most jurisdictions in the United States exclude individuals who continue to shed certain organisms from sensitive occupations and situations, such as food handling, providing direct patient care, or attending day care. The burden that this creates for recovering patients and their families/coworkers is great, so any effort to provide efficiency to the testing process would be of significant benefit. We sought to assess the ability of PCR for the detection of Salmonella enterica shedding and to compare that ability to culture-based testing. PCR would be faster than culture and would allow results to be generated more quickly. Herein, we show data that indicate that, while PCR and culture testing agree in the majority of cases, there are incidents of discordance between the two tests, whereupon PCR shows positive results when culture indicates lack of detectable viable organisms. Using culture-based testing as the standard, the negative predictive value of PCR was found to be 100%, while the positive predictive value was 79%. The nature of this discordance is briefly investigated. We found that it is possible that PCR may not only detect nonviable organisms in stool but also viable organisms that remain undetectable by standard culture methods.

Implementation of a Rapid Genotypic Assay to Promote Targeted Ciprofloxacin Therapy of Neisseria gonorrhoeae in a Large Health System.

Multidrug-resistant Neisseria gonorrhoeae is a top threat to public health. In November 2015, UCLA Health introduced a rapid gyrase A (gyrA) genotypic assay for prediction of Neisseria gonorrhoeae susceptibility to ciprofloxacin. We found a significant reduction in ceftriaxone use with a concomitant increase in targeted therapy.

A Case of Persistent and Possibly Treatment Resistant Pharyngeal Gonorrhea.

An HIV-negative man with pharyngeal gonorrhea had a positive test-of-cure (nucleic acid amplification test) result 7 days after treatment with ceftriaxone/azithromycin. Neisseria gonorrhoeae Multi-Antigen Sequencing Type 1407 and mosaic pen A (XXXIV) gene were identified in the test-of-cure specimen, and culture was negative. Retreatment with ceftriaxone 500 mg intramuscularly plus azithromycin 2 g orally yielded a negative test-of-cure result.

Screening Yield of HIV Antigen/Antibody Combination and Pooled HIV RNA Testing for Acute HIV Infection in a High-Prevalence Population.

IMPORTANCE: Although acute HIV infection contributes disproportionately to onward HIV transmission, HIV testing has not routinely included screening for acute HIV infection. OBJECTIVE: To evaluate the performance of an HIV antigen/antibody (Ag/Ab) combination assay to detect acute HIV infection compared with pooled HIV RNA testing. DESIGN, SETTING, AND PARTICIPANTS: Multisite, prospective, within-individual comparison study conducted between September 2011 and October 2013 in 7 sexually transmitted infection clinics and 5 community-based programs in New York, California, and North Carolina. Participants were 12 years or older and seeking HIV testing, without known HIV infection. EXPOSURES: All participants with a negative rapid HIV test result were screened for acute HIV infection with an HIV Ag/Ab combination assay (index test) and pooled human immunodeficiency virus 1 (HIV-1) RNA testing. HIV RNA testing was the reference standard, with positive reference standard result defined as detectable HIV-1 RNA on an individual RNA test. MAIN OUTCOMES AND MEASURES: Number and proportion with acute HIV infections detected. RESULTS: Among 86,836 participants with complete test results (median age, 29 years; 75.0% men; 51.8% men who have sex with men), established HIV infection was diagnosed in 1158 participants (1.33%) and acute HIV infection was diagnosed in 168 participants (0.19%). Acute HIV infection was detected in 134 participants with HIV Ag/Ab combination testing (0.15% [95% CI, 0.13%-0.18%]; sensitivity, 79.8% [95% CI, 72.9%-85.6%]; specificity, 99.9% [95% CI, 99.9%-99.9%]; positive predictive value, 59.0% [95% CI, 52.3%-65.5%]) and in 164 participants with pooled HIV RNA testing (0.19% [95% CI, 0.16%-0.22%]; sensitivity, 97.6% [95% CI, 94.0%-99.4%]; specificity, 100% [95% CI, 100%-100%]; positive predictive value, 96.5% [95% CI, 92.5%-98.7%]; sensitivity comparison, P < .001). Overall HIV Ag/Ab combination testing detected 82% of acute HIV infections detectable by pooled HIV RNA testing. Compared with rapid HIV testing alone, HIV Ag/Ab combination testing increased the relative HIV diagnostic yield (both established and acute HIV infections) by 10.4% (95% CI, 8.8%-12.2%) and pooled HIV RNA testing increased the relative HIV diagnostic yield by 12.4% (95% CI, 10.7%-14.3%). CONCLUSIONS AND RELEVANCE: In a high-prevalence population, HIV screening using an HIV Ag/Ab combination assay following a negative rapid test detected 82% of acute HIV infections detectable by pooled HIV RNA testing, with a positive predictive value of 59%. Further research is needed to evaluate this strategy in lower-prevalence populations and in persons using preexposure prophylaxis for HIV prevention.

Stemming the tide of drug-resistant Neisseria gonorrhoeae: the need for an individualized approach to treatment.

Drug-resistant Neisseria gonorrhoeae poses a significant public health challenge. In recent years, gonococci resistant to first- and second-line antibiotics have spread worldwide and new strains have developed that are increasingly resistant to third-generation cephalosporins, which are currently our last line of available treatments. Given the timeline required to develop new drugs or an effective vaccine for N. gonorrhoeae, a top priority is to use the drugs that are available as effectively as possible. Currently, clinical management of gonorrhoea is based upon treatment guidelines informed by international gonococcal antimicrobial susceptibility surveillance programmes. This approach, although currently the most practical, is subject to a number of limitations since surveillance data inherently provide population-level information. As a result, basing treatment guidelines on these data can result in the prescription of more aggressive or broader treatment than is needed by individual patients and hence inadvertently contribute to the development and spread of resistance to important drugs. Clearly, methods are needed that provide patient-specific drug susceptibility information in a time frame that would allow clinicians to prescribe individualized treatment regimens for gonorrhoea. Fortunately, in recent years, there have been a number of advances in the development of rapid methods for characterizing both the genotype and the drug resistance phenotype of N. gonorrhoeae strains. Here, we review these advances and propose additional studies that would help facilitate a transition towards an individualized treatment approach for gonorrhoea.

Fluoroquinolone resistance in Neisseria gonorrhoeae after cessation of ciprofloxacin usage in San Francisco: using molecular typing to investigate strain turnover.

BACKGROUND: Ciprofloxacin resistance (CipR) among gonococcal strains in San Francisco (SF) increased between 2001 and 2006 and decreased between 2007 and 2009. Molecular typing of isolates obtained from 2005 to 2009 was performed to elucidate changes in CipR prevalence. METHODS: A total of 2526 samples were collected at the SF City Clinic between 2001 and 2009. Minimum inhibitory concentrations to ciprofloxacin were obtained by agar dilution. Prevalences of CipR strains were determined, with corresponding confidence intervals (CIs). Between 2005 and 2009, 460 isolates were selected for molecular typing using Neisseria gonorrhoeae multiantigen sequence typing. RESULTS: Between 2001 and 2006, the prevalence of CipR increased from 3.4% (95% CI, 1.3%-5.4%) to 44% (95% CI, 39%-50%). However, in 2007 prevalence began to decrease, reaching 9.6% (95% CI, 6.0%-13%) by 2009. Of the 203 strain types identified between 2005 and 2009, 126 genogroups of closely related strain types were formed (varying by ≤1% at both target loci). Levels of CipR within the data set correlate with the prevalence of 3 major genogroups (G): G437, G1407, and G3112. CONCLUSIONS: Molecular typing reveals that CipR within the tested population is maintained by strain turnover between resistant genogroups. Despite early recommendation in 2002 to stop ciprofloxacin use in California, CipR in SF increased through 2006. The subsequent decrease in CipR corresponds with the 2007 national recommendation to cease ciprofloxacin treatment of gonorrhea, which suggests that national recommendations are potentially more effective at reducing CipR than regional recommendations in areas with high strain turnover.

Impact of GeneXpert MTB/RIF on patients and tuberculosis programs in a low-burden setting. a hypothetical trial.

RATIONALE: Guidelines recommend routine nucleic-acid amplification testing in patients with presumed tuberculosis (TB), but these tests have not been widely adopted. GeneXpert MTB/RIF (Xpert), a novel, semiautomated TB nucleic-acid amplification test, has renewed interest in this technology, but data from low-burden countries are limited. OBJECTIVES: We sought to estimate Xpert's potential clinical and public health impact on empiric treatment, contact investigation, and housing in patients undergoing TB evaluation. METHODS: We performed a prospective, cross-sectional study with 2-month follow-up comparing Xpert with standard strategies for evaluating outpatients for active pulmonary TB at the San Francisco Department of Public Health TB Clinic between May 2010 and June 2011. We calculated the diagnostic accuracy of standard algorithms for initial empiric TB treatment, contact investigation, and housing in reference to three Mycobacterium tuberculosis sputum cultures, as compared with that of a single sputum Xpert test. We estimated the incremental diagnostic value of Xpert, and the hypothetical reductions in unnecessary treatment, contact investigation, and housing if Xpert were adopted to guide management decisions. MEASUREMENTS AND MAIN RESULTS: A total of 156 patients underwent Xpert testing. Fifty-nine (38%) received empiric TB treatment. Thirteen (8%) had culture-positive TB. Xpert-guided management would have hypothetically decreased overtreatment by 94%, eliminating a median of 44 overtreatment days (interquartile range, 43-47) per patient and 2,169 total overtreatment days (95% confidence interval, 1,938-2,400) annually, without reducing early detection of TB patients. We projected similar benefits for contact investigation and housing. CONCLUSIONS: Xpert could greatly reduce the frequency and impact of unnecessary empiric treatment, contact investigation, and housing, providing substantial patient and programmatic benefits if used in management decisions.

Neisseria gonorrhoeae strain with reduced susceptibilities to extended-spectrum cephalosporins.

The spread of Neisseria gonorrhoeae strains with reduced susceptibility to extended-spectrum cephalosporins is an increasing public health threat. Using Etest and multiantigen sequence typing, we detected sequence type 1407, which is associated with reduced susceptibilities to extended-spectrum cephalosporins, in 4 major populated regions in California, USA, in 2012.

Lymphogranuloma venereum causing a persistent genital ulcer.

Lymphogranuloma venereum (LGV) is a sexually transmitted cause of inguinal lymphadenopathy and proctocolitis. We report a patient with a persistent genital ulcer due to LGV (serovar L2b), an unusual presentation among US men who have sex with men. Lymphogranuloma venereum should be considered when evaluating persistent genital ulcers, and LGV-specific testing should be sought.

Characteristics of males infected with common Neisseria gonorrhoeae sequence types in the Gonococcal Isolate Surveillance Project, San Francisco, California, 2009.

We analyzed 265 urethral Neisseria gonorrhoeae specimens collected from symptomatic males at San Francisco's municipal sexually transmitted disease clinic, a participant in the Gonococcal Isolate Surveillance Project, during 2009. We used N. gonorrhoeae multiantigen sequence typing to describe characteristics of patients infected with common sequence type families. Specimens were classified into 6 homology-based families and 1 additional family of all other identified strains. Strain family results were combined with results of culture-based antibiotic sensitivity minimum inhibitory concentration, sociodemographic and behavioral risk data collected at the clinic, and presence or absence of the mosaic penicillin-binding protein 2 (penA) allele. Characteristics of patients were compared across strain families through the use of χ(2) statistics. Among men who have sex with men, strain distribution differed by those reporting receptive oral sex as their only urethral exposure (P = 0.04), by number of sex partners (P = 0.03), and by race/ethnicity (P < 0.001); there were no differences by age or human immunodeficiency virus status. Also, among men who have sex with men, strain family distributions differed for culture specimens with reduced susceptibility to a range of antibiotics, as well as with presence of the mosaic penA allele (all P < 0.001). The combination of molecular, phenotypic, and epidemiologic data on N. gonorrhoeae infection could help develop a more complete epidemiology of gonorrhea in the United States.

Neisseria gonorrhoeae and extended-spectrum cephalosporins in California: surveillance and molecular detection of mosaic penA.

BACKGROUND: The spread of Neisseria gonorrhoeae strains with mosaic penA alleles and reduced susceptibility to extended-spectrum cephalosporins is a major public health problem. While much work has been performed internationally, little is known about the genetics or molecular epidemiology of N. gonorrhoeae isolates with reduced susceptibility to extended-spectrum cephalosporins in the United States. The majority of N. gonorrhoeae infections are diagnosed without a live culture. Molecular tools capable of detecting markers of extended-spectrum cephalosporin resistance are needed. METHODS: Urethral N. gonorrhoeae isolates were collected from 684 men at public health clinics in California in 2011. Minimum inhibitory concentrations (MICs) to ceftriaxone, cefixime, cefpodoxime and azithromycin were determined by Etest and categorized according to the U.S. Centers for Disease Control 2010 alert value breakpoints. 684 isolates were screened for mosaic penA alleles using real-time PCR (RTPCR) and 59 reactive isolates were subjected to DNA sequencing of their penA alleles and Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST). To increase the specificity of the screening RTPCR in detecting isolates with alert value extended-spectrum cephalosporin MICs, the primers were modified to selectively amplify the mosaic XXXIV penA allele. RESULTS: Three mosaic penA alleles were detected including two previously described alleles (XXXIV, XXXVIII) and one novel allele (LA-A). Of the 29 isolates with an alert value extended-spectrum cephalosporin MIC, all possessed the mosaic XXXIV penA allele and 18 were sequence type 1407, an internationally successful strain associated with multi-drug resistance. The modified RTPCR detected the mosaic XXXIV penA allele in urethral isolates and urine specimens and displayed no amplification of the other penA alleles detected in this study. CONCLUSION: N. gonorrhoeae isolates with mosaic penA alleles and reduced susceptibility to extended-spectrum cephalosporins are currently circulating in California. Isolates with the same NG-MAST ST, penA allele and extended-spectrum cephalosporin MICs have caused treatment failures elsewhere. The RTPCR assay presented here may be useful for the detection of N. gonorrheoae isolates and clinical specimens with reduced extended-spectrum cephalosporin MICs in settings where antimicrobial susceptibility testing is unavailable. In an era of increasing antimicrobial resistance and decreasing culture capacity, molecular assays capable of detecting extended-spectrum cephalosporin of resistance are essential to public health.

Analysis of SARS-CoV-2 variants from patient specimens in Nevada from October 2020 to August 2021.

In early 2020, the emergence and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the human population quickly developed into a global pandemic. SARS-CoV-2 is the etiological agent of coronavirus disease 2019 (COVID-19) which has a broad range of respiratory illnesses. As the virus circulates, it acquires nucleotide changes. These mutations are potentially due to the inherent differences in the selection pressures within the human population compared to the original zoonotic reservoir of SARS-CoV-2 and formerly naïve humans. The acquired mutations will most likely be neutral, but some may have implications for viral transmission, disease severity, and resistance to therapies or vaccines. This is a follow-up study from our early report (Hartley et al. J Genet Genomics. 01202021;48(1):40-51) which detected a rare variant (nsp12, RdRp P323F) circulating within Nevada in mid 2020 at high frequency. The primary goals of the current study were to determine the phylogenetic relationship of the SARS-CoV-2 genomes within Nevada and to determine if there are any unusual variants within Nevada compared to the current database of SARS-CoV-2 sequences. Whole genome sequencing and analysis of SARS-CoV-2 from 425 positively identified nasopharyngeal/nasal swab specimens were performed from October 2020 to August 2021 to determine any variants that could result in potential escape from current therapeutics. Our analysis focused on nucleotide mutations that generated amino acid variations in the viral Spike (S) protein, Receptor binding domain (RBD), and the RNA-dependent RNA-polymerase (RdRp) complex. The data indicate that SARS-CoV-2 sequences from Nevada did not contain any unusual variants that had not been previously reported. Additionally, we did not detect the previously identified the RdRp P323F variant in any of the samples. This suggests that the rare variant we detected before was only able to circulate because of the stay-at-home orders and semi-isolation experience during the early months of the pandemic. IMPORTANCE: SARS-COV-2 continues to circulate in the human population. In this study, SARS-CoV-2 positive nasopharyngeal/nasal swab samples were used for whole genome sequencing to determine the phylogenetic relationship of SARS-CoV-2 sequences within Nevada from October 2020 to August 2021. The resulting data is being added to a continually growing database of SARS-CoV-2 sequences that will be important for understanding the transmission and evolution of the virus as it spreads around the globe.

GAMBIT (Genomic Approximation Method for Bacterial Identification and Tracking): A methodology to rapidly leverage whole genome sequencing of bacterial isolates for clinical identification.

Whole genome sequencing (WGS) of clinical bacterial isolates has the potential to transform the fields of diagnostics and public health. To realize this potential, bioinformatic software that reports identification results needs to be developed that meets the quality standards of a diagnostic test. We developed GAMBIT (Genomic Approximation Method for Bacterial Identification and Tracking) using k-mer based strategies for identification of bacteria based on WGS reads. GAMBIT incorporates this algorithm with a highly curated searchable database of 48,224 genomes. Herein, we describe validation of the scoring methodology, parameter robustness, establishment of confidence thresholds and the curation of the reference database. We assessed GAMBIT by way of validation studies when it was deployed as a laboratory-developed test in two public health laboratories. This method greatly reduces or eliminates false identifications which are often detrimental in a clinical setting.

Rapid Lineage Assignment of Severe Acute Respiratory Syndrome Coronavirus 2 Cases through Automated Library Preparation, Sequencing, and Bioinformatic Analysis.

The coronavirus disease 2019 (COVID-19) pandemic has provided a stage to illustrate that there is considerable value in obtaining rapid, whole-genome-based information about pathogens. This article describes the utility of a commercially available, automated severe acute respiratory syndrome associated coronavirus 2 (SARS-CoV-2) library preparation, genome sequencing, and a bioinformatics analysis pipeline to provide rapid, near-real-time SARS-CoV-2 variant description. This study evaluated the turnaround time, accuracy, and other quality-related parameters obtained from commercially available automated sequencing instrumentation, from analysis of continuous clinical samples obtained from January 1, 2021, to October 6, 2021. This analysis included a base-by-base assessment of sequencing accuracy at every position in the SARS-CoV-2 chromosome using two commercially available methods. Mean turnaround time, from the receipt of a specimen for SARS-CoV-2 testing to the availability of the results, with lineage assignment, was <3 days. Accuracy of sequencing by one method was 100%, although certain sites on the genome were found repeatedly to have been sequenced with varying degrees of read error rate.

The use of whole-genome sequencing and development of bioinformatics to monitor overlapping outbreaks of Candida auris in southern Nevada.

A Candida auris outbreak has been ongoing in Southern Nevada since August 2021. In this manuscript we describe the sequencing of over 200 C. auris isolates from patients at several facilities. Genetically distinct subgroups of C. auris were detected from Clade I (3 distinct lineages) and III (1 lineage). Open-source bioinformatic tools were developed and implemented to aid in the epidemiological investigation. The work herein compares three methods for C. auris whole genome analysis: Nullarbor, MycoSNP and a new pipeline TheiaEuk. We also describe a novel analysis method focused on elucidating phylogenetic linkages between isolates within an ongoing outbreak. Moreover, this study places the ongoing outbreaks in a global context utilizing existing sequences provided worldwide. Lastly, we describe how the generated results were communicated to the epidemiologists and infection control to generate public health interventions.

TheiaEuk: a species-agnostic bioinformatics workflow for fungal genomic characterization.

Introduction: The clinical incidence of antimicrobial-resistant fungal infections has dramatically increased in recent years. Certain fungal pathogens colonize various body cavities, leading to life-threatening bloodstream infections. However, the identification and characterization of fungal isolates in laboratories remain a significant diagnostic challenge in medicine and public health. Whole-genome sequencing provides an unbiased and uniform identification pipeline for fungal pathogens but most bioinformatic analysis pipelines focus on prokaryotic species. To this end, TheiaEuk_Illumina_PE_PHB (TheiaEuk) was designed to focus on genomic analysis specialized to fungal pathogens. Methods: TheiaEuk was designed using containerized components and written in the workflow description language (WDL) to facilitate deployment on the cloud-based open bioinformatics platform Terra. This species-agnostic workflow enables the analysis of fungal genomes without requiring coding, thereby reducing the entry barrier for laboratory scientists. To demonstrate the usefulness of this pipeline, an ongoing outbreak of C. auris in southern Nevada was investigated. We performed whole-genome sequence analysis of 752 new C. auris isolates from this outbreak. Furthermore, TheiaEuk was utilized to observe the accumulation of mutations in the FKS1 gene over the course of the outbreak, highlighting the utility of TheiaEuk as a monitor of emerging public health threats when combined with whole-genome sequencing surveillance of fungal pathogens. Results: A primary result of this work is a curated fungal database containing 5,667 unique genomes representing 245 species. TheiaEuk also incorporates taxon-specific submodules for specific species, including clade-typing for Candida auris (C. auris). In addition, for several fungal species, it performs dynamic reference genome selection and variant calling, reporting mutations found in genes currently associated with antifungal resistance (FKS1, ERG11, FUR1). Using genome assemblies from the ATCC Mycology collection, the taxonomic identification module used by TheiaEuk correctly assigned genomes to the species level in 126/135 (93.3%) instances and to the genus level in 131/135 (97%) of instances, and provided zero false calls. Application of TheiaEuk to actual specimens obtained in the course of work at a local public health laboratory resulted in 13/15 (86.7%) correct calls at the species level, with 2/15 called at the genus level. It made zero incorrect calls. TheiaEuk accurately assessed clade type of Candida auris in 297/302 (98.3%) of instances. Discussion: TheiaEuk demonstrated effectiveness in identifying fungal species from whole genome sequence. It further showed accuracy in both clade-typing of C. auris and in the identification of mutations known to associate with drug resistance in that organism.