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.

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.

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.

Candida auris Discovery through Community Wastewater Surveillance during Healthcare Outbreak, Nevada, USA, 2022.

Candida auris transmission is steadily increasing across the United States. We report culture-based detection of C. auris in wastewater and the epidemiologic link between isolated strains and southern Nevada, USA, hospitals within the sampled sewershed. Our results illustrate the potential of wastewater surveillance for containing C. auris.

Comparison of a Prototype SARS-CoV-2 Lateral Flow IMMUNOASSAY with the BinaxNOWTM COVID-19 Antigen CARD.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the virus responsible for the COVID-19 pandemic. From the onset of the pandemic, rapid antigen tests have quickly proved themselves to be an accurate and accessible diagnostic platform. The initial (and still most commonly used antigen tests) for COVID-19 diagnosis were constructed using monoclonal antibodies (mAbs) specific to severe acute respiratory syndrome coronavirus (SARS-CoV) nucleocapsid protein (NP). These mAbs are able to bind SARS-CoV-2 NP due to high homology between the two viruses. However, since first being identified in 2019, SARS-CoV-2 has continuously mutated, and a multitude of variants have appeared. These mutations have an elevated risk of leading to possible diagnostic escape when using tests produced with SARS-CoV-derived mAbs. Here, we established a library of 18 mAbs specific to SARS-CoV-2 NP and used two of these mAbs (1CV7 and 1CV14) to generate a prototype antigen-detection lateral flow immunoassay (LFI). A side-by-side analysis of the 1CV7/1CV14 LFI and the commercially available BinaxNOWTM COVID-19 Antigen CARD was performed. Results indicated the 1CV7/1CV14 LFI outperformed the BinaxNOWTM test in the detection of BA.2, BA.2.12.1, and BA.5 Omicron sub-variants when testing remnant RT-PCR positive patient nasopharyngeal swabs diluted in viral transport media.

Detecting SARS-CoV-2 variants in wastewater and their correlation with circulating variants in the communities.

Detection of SARS-CoV-2 viral load in wastewater has been highly informative in estimating the approximate number of infected individuals in the surrounding communities. Recent developments in wastewater monitoring to determine community prevalence of COVID-19 further extends into identifying SARS-CoV-2 variants, including those being monitored for having enhanced transmissibility. We sequenced genomic RNA derived from wastewater to determine the variants of coronaviruses circulating in the communities. Wastewater samples were collected from Truckee Meadows Water Reclamation Facility (TMWRF) from November 2020 to June 2021. SARS-CoV-2 variants resulting from wastewater were compared with the variants detected in infected individuals' clinical specimens (nasal/nasopharyngeal swabs) during the same period and found conclusively in agreement. Therefore, wastewater monitoring for SARS-CoV-2 variants in the community is a feasible strategy as a complementary tool to clinical specimen testing in the latter's absence.

Rapid repeat infection of SARS-CoV-2 by two highly distinct delta-lineage viruses.

An instance of sequential infection of an individual with, firstly, the Delta variant and secondly a Delta-sub-lineage has been identified. The individual was found positive for the AY.26 lineage 22 days after being found positive for the Delta [B.1.617.2] variant. The viruses associated with the cases showed dramatic genomic difference, including 31 changes that resulted in deletions or amino acid substitutions. Seven of these differences were observed in the Spike protein. The patient in question was between 30 and 35 years old and had no underlying health conditions. Though singular, this case illustrates the possibility that infection with the Delta variant may not itself be fully protective against a population of SARS-CoV-2 variants that are becoming increasingly diverse.

Rapid Diagnostic Testing for Response to the Monkeypox Outbreak - Laboratory Response Network, United States, May 17-June 30, 2022.

As part of public health preparedness for infectious disease threats, CDC collaborates with other U.S. public health officials to ensure that the Laboratory Response Network (LRN) has diagnostic tools to detect Orthopoxviruses, the genus that includes Variola virus, the causative agent of smallpox. LRN is a network of state and local public health, federal, U.S. Department of Defense (DOD), veterinary, food, and environmental testing laboratories. CDC developed, and the Food and Drug Administration (FDA) granted 510(k) clearance* for the Non-variola Orthopoxvirus Real-time PCR Primer and Probe Set (non-variola Orthopoxvirus [NVO] assay), a polymerase chain reaction (PCR) diagnostic test to detect NVO. On May 17, 2022, CDC was contacted by the Massachusetts Department of Public Health (DPH) regarding a suspected case of monkeypox, a disease caused by the Orthopoxvirus Monkeypox virus. Specimens were collected and tested by the Massachusetts DPH public health laboratory with LRN testing capability using the NVO assay. Nationwide, 68 LRN laboratories had capacity to test approximately 8,000 NVO tests per week during June. During May 17-June 30, LRN laboratories tested 2,009 specimens from suspected monkeypox cases. Among those, 730 (36.3%) specimens from 395 patients were positive for NVO. NVO-positive specimens from 159 persons were confirmed by CDC to be monkeypox; final characterization is pending for 236. Prompt identification of persons with infection allowed rapid response to the outbreak, including isolation and treatment of patients, administration of vaccines, and other public health action. To further facilitate access to testing and increase convenience for providers and patients by using existing provider-laboratory relationships, CDC and LRN are supporting five large commercial laboratories with a national footprint (Aegis Science, LabCorp, Mayo Clinic Laboratories, Quest Diagnostics, and Sonic Healthcare) to establish NVO testing capacity of 10,000 specimens per week per laboratory. On July 6, 2022, the first commercial laboratory began accepting specimens for NVO testing based on clinician orders.

Use of wastewater surveillance for early detection of Alpha and Epsilon SARS-CoV-2 variants of concern and estimation of overall COVID-19 infection burden.

A decline in diagnostic testing for SARS-CoV-2 is expected to delay the tracking of COVID-19 variants of concern and interest in the United States. We hypothesize that wastewater surveillance programs provide an effective alternative for detecting emerging variants and assessing COVID-19 incidence, particularly when clinical surveillance is limited. Here, we analyzed SARS-CoV-2 RNA in wastewater from eight locations across Southern Nevada between March 2020 and April 2021. Trends in SARS-CoV-2 RNA concentrations (ranging from 4.3 log10 gc/L to 8.7 log10 gc/L) matched trends in confirmed COVID-19 incidence, but wastewater surveillance also highlighted several limitations with the clinical data. Amplicon-based whole genome sequencing (WGS) of 86 wastewater samples identified the B.1.1.7 (Alpha) and B.1.429 (Epsilon) lineages in December 2020, but clinical sequencing failed to identify the variants until January 2021, thereby demonstrating that 'pooled' wastewater samples can sometimes expedite variant detection. Also, by calibrating fecal shedding (11.4 log10 gc/infection) and wastewater surveillance data to reported seroprevalence, we estimate that ~38% of individuals in Southern Nevada had been infected by SARS-CoV-2 as of April 2021, which is significantly higher than the 10% of individuals confirmed through clinical testing. Sewershed-specific ascertainment ratios (i.e., X-fold infection undercounts) ranged from 1.0 to 7.7, potentially due to demographic differences. Our data underscore the growing application of wastewater surveillance in not only the identification and quantification of infectious agents, but also the detection of variants of concern that may be missed when diagnostic testing is limited or unavailable.

Detecting SARS-CoV-2 Variants in Wastewater and Their Correlation With Circulating Variants in the Communities.

Detection of SARS-CoV-2 viral load in wastewater has been highly informative in estimating the approximate number of infected individuals in the surrounding communities. Recent developments in wastewater monitoring to determine community prevalence of COVID-19 further extends into identifying SARS-CoV-2 variants, including those being monitored for having enhanced transmissibility. We sequenced genomic RNA derived from wastewater to determine the variants of coronaviruses circulating in the communities. Wastewater samples were collected from Truckee Meadows Water Reclamation Facility (TMWRF) from November 2021 to June 2021 were analyzed for SARS-CoV-2 variants and were compared with the variants detected in the clinical specimens (nasal/nasopharyngeal swabs) of infected individuals during the same period. The comparison was found to be conclusively in agreement. Therefore, wastewater monitoring for SARS-CoV-2 variants in the community is a feasible strategy both as a complementary tool to clinical specimen testing and in the latter's absence.

Genomic surveillance of Nevada patients revealed prevalence of unique SARS-CoV-2 variants bearing mutations in the RdRp gene.

Patients with signs of COVID-19 were tested through diagnostic RT-PCR for SARS-CoV-2 using RNA extracted from the nasopharyngeal/nasal swabs. To determine the variants of SARS-CoV-2 circulating in the state of Nevada, specimens from 200 COVID-19 patients were sequenced through our robust sequencing platform, which enabled sequencing of SARS-CoV-2 from specimens with even very low viral loads, without the need of culture-based amplification. High genome coverage allowed the identification of single and multi-nucleotide variants in SARS-CoV-2 in the community and their phylogenetic relationships with other variants present during the same period of the outbreak. We report the occurrence of a novel mutation at 323aa (314aa of orf1b) of nsp12 (RNA-dependent RNA polymerase) changed to phenylalanine (F) from proline (P), in the first reported isolate of SARS-CoV-2, Wuhan-Hu-1. This 323F variant was present at a very high frequency in Northern Nevada. Structural modeling determined this mutation in the interface domain, which is important for the association of accessory proteins required for the polymerase. In conclusion, we report the introduction of specific SARS-CoV-2 variants at very high frequency in distinct geographic locations, which is important for understanding the evolution and circulation of SARS-CoV-2 variants of public health importance, while it circulates in humans.

Genomic evidence for reinfection with SARS-CoV-2: a case study.

BACKGROUND: The degree of protective immunity conferred by infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is currently unknown. As such, the possibility of reinfection with SARS-CoV-2 is not well understood. We describe an investigation of two instances of SARS-CoV-2 infection in the same individual. METHODS: A 25-year-old man who was a resident of Washoe County in the US state of Nevada presented to health authorities on two occasions with symptoms of viral infection, once at a community testing event in April, 2020, and a second time to primary care then hospital at the end of May and beginning of June, 2020. Nasopharyngeal swabs were obtained from the patient at each presentation and twice during follow-up. Nucleic acid amplification testing was done to confirm SARS-CoV-2 infection. We did next-generation sequencing of SARS-CoV-2 extracted from nasopharyngeal swabs. Sequence data were assessed by two different bioinformatic methodologies. A short tandem repeat marker was used for fragment analysis to confirm that samples from both infections came from the same individual. FINDINGS: The patient had two positive tests for SARS-CoV-2, the first on April 18, 2020, and the second on June 5, 2020, separated by two negative tests done during follow-up in May, 2020. Genomic analysis of SARS-CoV-2 showed genetically significant differences between each variant associated with each instance of infection. The second infection was symptomatically more severe than the first. INTERPRETATION: Genetic discordance of the two SARS-CoV-2 specimens was greater than could be accounted for by short-term in vivo evolution. These findings suggest that the patient was infected by SARS-CoV-2 on two separate occasions by a genetically distinct virus. Thus, previous exposure to SARS-CoV-2 might not guarantee total immunity in all cases. All individuals, whether previously diagnosed with COVID-19 or not, should take identical precautions to avoid infection with SARS-CoV-2. The implications of reinfections could be relevant for vaccine development and application. FUNDING: Nevada IDEA Network of Biomedical Research, and the National Institute of General Medical Sciences (National Institutes of Health).

Characteristics of viral specimens collected from asymptomatic and fatal cases of COVID-19.

We sought to determine the characteristics of viral specimens associated with fatal cases, asymptomatic cases and non-fatal symptomatic cases of COVID-19. This included the analysis of 1264 specimens found reactive for at least two SARS-CoV-2 specific loci from people screened for infection in Northern Nevada in March-May of 2020. Of these, 30 were specimens from fatal cases, while 23 were from positive, asymptomatic cases. We assessed the relative amounts of SARS-CoV-2 RNA from sample swabs by real-time PCR and use of the threshold crossing value (Ct). Moreover, we compared the amount of human RNase P found on the same swabs. A considerably higher viral load was found to be associated with swabs from cases involving fatality and the difference was found to be strongly statistically significant. Noting this difference, we sought to assess whether any genetic correlation could be found in association with virus from fatal cases using whole genome sequencing. While no common genetic elements were discerned, one branch of epidemiologically linked fatal cases did have two point mutations, which no other of 156 sequenced cases from northern Nevada had. The mutations caused amino acid changes in the 3'-5' exonuclease protein, and the product of the gene, orf8.

Genomic surveillance of Nevada patients revealed prevalence of unique SARS-CoV-2 variants bearing mutations in the RdRp gene.

Patients with signs of COVID-19 were tested with CDC approved diagnostic RT-PCR for SARS-CoV-2 using RNA extracted from nasopharyngeal/nasal swabs. In order to determine the variants of SARS-CoV-2 circulating in the state of Nevada, 200 patient specimens from COVID-19 patients were sequenced through our robust protocol for sequencing SARS-CoV-2 genomes. Our protocol enabled sequencing of SARS-CoV-2 genome directly from the specimens, with even very low viral loads, without the need of culture-based amplification. This allowed the identification of specific nucleotide variants including those coding for D614G and clades defining mutations. These sequences were further analyzed for determining SARS-CoV-2 variants circulating in the state of Nevada and their phylogenetic relationships with other variants present in the united states and the world during the same period of the outbreak. Our study reports the occurrence of a novel variant in the nsp12 (RNA dependent RNA Polymerase) protein at residue 323 (314aa of orf1b) to Phenylalanine (F) from Proline (P), present in the original isolate of SARS-CoV-2 (Wuhan-Hu-1). This 323F variant is found at a very high frequency (46% of the tested specimen) in Northern Nevada. Functional significance of this unique and highly prevalent variant of SARS-CoV-2 with RdRp mutation is currently under investigation but structural modeling showed this 323aa residue in the interface domain of RdRp, which is required for association with accessory proteins. In conclusion, we report the introduction of specific SARS-CoV-2 variants at a very high frequency within a distinct geographic location, which is important for clinical and public health perspectives in understanding the evolution of SARS-CoV-2 while in circulation.

High-Throughput Transcription-mediated amplification on the Hologic Panther is a highly sensitive method of detection for SARS-CoV-2.

BACKGROUND: As the demand for laboratory testing for SARS-CoV-2 increases, additional varieties of testing methodologies are being considered. While real time polymerase chain reaction (RT-PCR) has performed as the main method for virus detection, other methods are becoming available, including transcription mediated amplification (TMA). The Hologic Aptima SARS-CoV-2 Assay utilizes TMA as a target amplification mechanism, and it has only recently received Emergency Use Authorization (EUA) by the Food and Drug Administration (FDA). OBJECTIVES: We sought to compare the sensitivity and specificity of the Aptima SARS-CoV-2 Assay to RTPCR as a means of SARS-CoV-2 detection in a diagnostic setting. STUDY DESIGN: We performed a limit-of-detection study (LoD) to assess the analytical sensitivity of TMA and RT-PCR. This preceded a comparison of the methods using previously evaluated clinical specimens (nasopharyngeal swabs) using 116 human specimens tested by both methodologies. Specimens included sixty-one (61) specimens found reactive by real-time PCR, fifty-one (51) found non-reactive, and four (4) deemed inconclusive. RESULTS: The Aptima SARS-CoV-2 Assay showed a markedly higher analytical sensitivity than RT-PCR by LoD study. Evaluation of clinical specimens resulted in fewer inconclusive results by the SARS-CoV-2 assay, leading to potentially higher clinical sensitivity. CONCLUSIONS: Higher analytical sensitivity may explain TMA's ability to ascertain for the presence of SARS-CoV-2 genome in human specimens deemed inconclusive by real-time PCR. TMA provides an effective, highly sensitive means of detection of SARS-CoV-2 in nasopharyngeal specimens.

Transmission of novel Klebsiella pneumoniae carbapenemase-producing Escherichia coli sequence type 1193 among residents and caregivers in a community-based, residential care setting-Nevada, 2018.

We describe transmission of Klebsiella pneumoniae carbapenemase-producing Escherichia coli sequence type (ST) 1193 in a group home. E. coli ST1193 is an emerging multidrug-resistant clone not previously shown to carry carbapenemases in the United States. Our investigation illustrates the potential of residential group homes to amplify rare combinations of pathogens and resistance mechanisms.

Efficient genome editing by CRISPR-Mb3Cas12a in mice.

As an alternative and complementary approach to Cas9-based genome editing, Cas12a has not been widely used in mammalian cells largely due to its strict requirement for the TTTV protospacer adjacent motif (PAM) sequence. Here, we report that Mb3Cas12a (Moraxella bovoculi AAX11_00205) can efficiently edit the mouse genome based on the TTV PAM sequence with minimal numbers of large on-target deletions or insertions. When TTTV PAM sequence-targeting CRISPR (cr)RNAs of 23 nt spacers are used, >70% of the founders obtained are edited. Moreover, the use of Mb3Cas12a tagged to monomeric streptavidin (mSA) in conjunction with biotinylated DNA donor template leads to high knock-in efficiency in two-cell mouse embryos, with 40% of founders obtained containing the desired knock-in sequences.