Large Clusters of Invasive emm49 Group A Streptococcus Identified Within Arizona Health Care Facilities Through Statewide Genomic Surveillance System, 2019-2021.

A statewide genomic surveillance system for invasive Group A Streptococcus was implemented in Arizona in June 2019, resulting in 1046 isolates being submitted for genomic analysis to characterize emm types and identify transmission clusters. Eleven of the 32 identified distinct emm types comprised >80% of samples, with 29.7% of all isolates being typed as emm49 (and its genetic derivative emm151). Phylogenetic analysis initially identified an emm49 genomic cluster of 4 isolates that rapidly expanded over subsequent months (June 2019 to February 2020). Public health investigations identified epidemiologic links with 3 different long-term care facilities, resulting in specific interventions. Unbiased genomic surveillance allowed for identification and response to clusters that would have otherwise remained undetected.

Drug resistance and epidemiological success of modern Mycobacterium tuberculosis lineages in western India.

BACKGROUND: Transmission is contributing to the slow decline of tuberculosis (TB) incidence globally. Drivers of TB transmission in India, the country estimated to carry a quarter of the World's burden, are not well studied. We conducted a genomic epidemiology study to compare epidemiological success, host factors and drug resistance (DR) among the four major Mycobacterium tuberculosis (Mtb) lineages (L1-4) circulating in Pune, India. METHODS: We performed whole-genome sequencing (WGS) of Mtb sputum culture-positive isolates from participants in two prospective cohort studies and predicted genotypic susceptibility using a validated random forest model. We used maximum likelihood estimation to build phylogenies. We compared lineage specific phylogenetic and time-scaled metrics to assess epidemiological success. RESULTS: Of the 642 isolates that underwent WGS, 612 met sequence quality criteria. Most isolates belonged to L3 (44.6%). The majority (61.1%) of multidrug-resistant isolates belonged to L2 (P < 0.001). In molecular dating, L2 demonstrated a higher rate and more recent resistance acquisition. We measured higher clustering, and time-scaled haplotypic density (THD) for L4 and L2 compared to L3 and/or L1 suggesting higher epidemiological success. L4 demonstrated higher THD and clustering (OR 5.1 (95% CI 2.3-12.3) in multivariate models controlling for host factors and DR. CONCLUSION: L2 shows a higher frequency of DR and both L2 and L4 demonstrate evidence of higher epidemiological success than L3 or L1 in the study setting. Our findings highlight the need for contact tracing around TB cases, and heightened surveillance of TB DR in India.

Sequencing by binding rivals SMOR error-corrected sequencing by synthesis technology for accurate detection and quantification of minor (< 0.1%) subpopulation variants.

BACKGROUND: Detecting very minor (< 1%) subpopulations using next-generation sequencing is a critical need for multiple applications, including the detection of drug resistant pathogens and somatic variant detection in oncology. A recently available sequencing approach termed 'sequencing by binding (SBB)' claims to have higher base calling accuracy data "out of the box." This paper evaluates the utility of using SBB for the detection of ultra-rare drug resistant subpopulations in Mycobacterium tuberculosis (Mtb) using a targeted amplicon assay and compares the performance of SBB to single molecule overlapping reads (SMOR) error corrected sequencing by synthesis (SBS) data. RESULTS: SBS displayed an elevated error rate when compared to SMOR error-corrected SBS and SBB techniques. SMOR error-corrected SBS and SBB technologies performed similarly within the linear range studies and error rate studies. CONCLUSIONS: With lower sequencing error rates within SBB sequencing, this technique looks promising for both targeted and unbiased whole genome sequencing, leading to the identification of minor (< 1%) subpopulations without the need for error correction methods.

DETECTION OF SARS-COV-2 IN A SQUIRREL MONKEY (SAIMIRI SCIUREUS): A ONE HEALTH INVESTIGATION AND RESPONSE.

Through collaborative efforts, One Health partners have responded to outbreaks of COVID-19 among animals, including those in human care at zoos. Zoos have been faced with numerous challenges, including the susceptibility of many mammalian species, and therefore the need to heighten biosecurity measures rapidly. Robust One Health collaborations already exist in Arizona to address endemic and emerging zoonoses, but these have rarely included zoos. The pandemic shed light on this, and Arizona subsequently expanded its SARS-CoV-2 surveillance efforts to include zoo animals. Testing and epidemiologic support was provided to expedite the detection of and response to zoonotic SARS-CoV-2 infection in zoo animals, as well as to understand possible transmission events. Resulting from this program, SARS-CoV-2 was detected from a rectal swab collected from an 8-yr-old squirrel monkey (Saimiri sciureus) from a zoo in Southern Arizona. The animal had rapidly become ill with nonrespiratory symptoms and died in July 2022. Genomic sequencing from the swab revealed mutations consistent with the Omicron (BA.2) lineage. An epidemiologic investigation identified an animal caretaker in close proximity to the affected squirrel monkey who tested positive for COVID-19 the same day the squirrel monkey died. Critical One Health partners provided support to the zoo through engagement of local, state, and federal agencies. Necropsy and pathologic evaluation showed significant necrotizing colitis; the overall clinical and histopathological findings did not implicate SARS-CoV-2 infection alone as a causal or contributing factor in the squirrel monkey's illness and death. This report documents the first identification of SARS-CoV-2 in a squirrel monkey and highlights a successful and timely One Health investigation conducted through multisectoral collaboration.

Genomic Epidemiology Linking Nonendemic Coccidioidomycosis to Travel.

Coccidioidomycosis is a fungal infection endemic to hot, arid regions of the western United States, northern Mexico, and parts of Central and South America. Sporadic cases outside these regions are likely travel-associated; alternatively, an infection could be acquired in as-yet unidentified newly endemic locales. A previous study of cases in nonendemic regions with patient self-reported travel history suggested that infections were acquired during travel to endemic regions. We sequenced 19 Coccidioides isolates from patients with known travel histories from that earlier investigation and performed phylogenetic analysis to identify the locations of potential source populations. Our results show that those isolates were phylogenetically linked to Coccidioides subpopulations naturally occurring in 1 of the reported travel locales, confirming that these cases were likely acquired during travel to endemic regions. Our findings demonstrate that genomic analysis is a useful tool for investigating travel-related coccidioidomycosis.

Phylogeography and transmission of M. tuberculosis in Moldova: A prospective genomic analysis.

BACKGROUND: The incidence of multidrug-resistant tuberculosis (MDR-TB) remains critically high in countries of the former Soviet Union, where >20% of new cases and >50% of previously treated cases have resistance to rifampin and isoniazid. Transmission of resistant strains, as opposed to resistance selected through inadequate treatment of drug-susceptible tuberculosis (TB), is the main driver of incident MDR-TB in these countries. METHODS AND FINDINGS: We conducted a prospective, genomic analysis of all culture-positive TB cases diagnosed in 2018 and 2019 in the Republic of Moldova. We used phylogenetic methods to identify putative transmission clusters; spatial and demographic data were analyzed to further describe local transmission of Mycobacterium tuberculosis. Of 2,236 participants, 779 (36%) had MDR-TB, of whom 386 (50%) had never been treated previously for TB. Moreover, 92% of multidrug-resistant M. tuberculosis strains belonged to putative transmission clusters. Phylogenetic reconstruction identified 3 large clades that were comprised nearly uniformly of MDR-TB: 2 of these clades were of Beijing lineage, and 1 of Ural lineage, and each had additional distinct clade-specific second-line drug resistance mutations and geographic distributions. Spatial and temporal proximity between pairs of cases within a cluster was associated with greater genomic similarity. Our study lasted for only 2 years, a relatively short duration compared with the natural history of TB, and, thus, the ability to infer the full extent of transmission is limited. CONCLUSIONS: The MDR-TB epidemic in Moldova is associated with the local transmission of multiple M. tuberculosis strains, including distinct clades of highly drug-resistant M. tuberculosis with varying geographic distributions and drug resistance profiles. This study demonstrates the role of comprehensive genomic surveillance for understanding the transmission of M. tuberculosis and highlights the urgency of interventions to interrupt transmission of highly drug-resistant M. tuberculosis.

Comparative Performance of Genomic Methods for the Detection of Pyrazinamide Resistance and Heteroresistance in Mycobacterium tuberculosis.

Pyrazinamide is an important component of both drug-susceptible and drug-resistant tuberculosis treatment regimens. Although approximately 50% of rifampin-resistant isolates are also resistant to pyrazinamide, pyrazinamide susceptibility testing is not routinely performed due to the challenging nature of the assay. We investigated the diagnostic accuracy of genotypic and phenotypic methods and explored the occurrence of pyrazinamide heteroresistance. We assessed pyrazinamide susceptibility among 358 individuals enrolled in the South African EXIT-RIF cohort using Sanger and targeted deep sequencing (TDS) of the pncA gene, whole-genome sequencing (WGS), and phenotypic drug susceptibility testing. We calculated the diagnostic accuracy of the different methods and investigated the prevalence and clinical impact of pncA heteroresistance. True pyrazinamide susceptibility status was assigned to each isolate using the Köser classification and expert rules. We observed 100% agreement across genotypic methods for detection of pncA fixed mutations; only TDS confidently identified three isolates (0.8%) with minor variants. For the 355 (99.2%) isolates that could be assigned true pyrazinamide status with confidence, phenotypic DST had a sensitivity of 96.5% (95% confidence interval [CI], 93.8 to 99.3%) and specificity of 100% (95% CI, 100 to 100%), both Sanger sequencing and WGS had a sensitivity of 97.1% (95% CI, 94.6 to 99.6%) and specificity of 97.8% (95% CI, 95.7 to 99.9%), and TDS had sensitivity of 98.8% (95% CI, 97.2 to 100%) and specificity of 97.8% (95% CI, 95.7 to 99.9%). We demonstrate high sensitivity and specificity for pyrazinamide susceptibility testing among all assessed genotypic methods. The prevalence of pyrazinamide heteroresistance in Mycobacterium tuberculosis isolates was lower than that identified for other first-line drugs.

Genotypic, proteomic, and phenotypic approaches to decipher the response to caspofungin and calcineurin inhibitors in clinical isolates of echinocandin-resistant Candida glabrata.

BACKGROUND: Echinocandin 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. OBJECTIVES: We 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. METHODS: Twenty-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. RESULTS: F659-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. CONCLUSIONS: Our findings confirm that calmodulin-calcineurin-Crz1 could provide a relevant target in life-threatening invasive candidiasis.

Translating RNA sequencing into clinical diagnostics: opportunities and challenges.

With the emergence of RNA sequencing (RNA-seq) technologies, RNA-based biomolecules hold expanded promise for their diagnostic, prognostic and therapeutic applicability in various diseases, including cancers and infectious diseases. Detection of gene fusions and differential expression of known disease-causing transcripts by RNA-seq represent some of the most immediate opportunities. However, it is the diversity of RNA species detected through RNA-seq that holds new promise for the multi-faceted clinical applicability of RNA-based measures, including the potential of extracellular RNAs as non-invasive diagnostic indicators of disease. Ongoing efforts towards the establishment of benchmark standards, assay optimization for clinical conditions and demonstration of assay reproducibility are required to expand the clinical utility of RNA-seq.

Suspected Locally Acquired Coccidioidomycosis in Human, Spokane, Washington, USA.

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.

Coyotes as Reservoirs for Onchocerca lupi, United States, 2015-2018.

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 microsporus Infections Associated with Surgical Procedures, Argentina, 2006-2014.

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.

Molecular Evaluation of Fluoroquinolone Resistance in Serial Mycobacterium tuberculosis Isolates from Individuals Diagnosed with Multidrug-Resistant Tuberculosis.

Fluoroquinolones (FQ) are crucial components of multidrug-resistant tuberculosis (MDR TB) treatment. Differing levels of resistance are associated with specific mutations within the quinolone-resistance-determining region (QRDR) of gyrA We sequenced the QRDR from serial isolates of MDR TB patients in the Preserving Effective TB Treatment Study (PETTS) with baseline FQ resistance (FQR) or acquired FQ resistance (FQACQR) using an Ion Torrent Personal Genome Machine (PGM) to a depth of 10,000× and reported single nucleotide polymorphisms in ≥1% of reads. FQR isolates harbored 15 distinct alleles with 1.3 (maximum = 6) on average per isolate. Eighteen alleles were identified in FQACQR isolates with an average of 1.6 (maximum = 9) per isolate. Isolates from 78% of FQACQR individuals had mutant alleles identified within 6 months of treatment initiation. Asp94Gly was the predominant allele in the initial FQ-resistant isolates followed by Ala90Val. Seventy-seven percent (36/47) of FQACQR group patients had isolates with FQ resistance alleles prior to changes to the FQ component of their treatment. Unlike the individuals treated initially with other FQs, none of the 21 individuals treated initially with levofloxacin developed genotypic or phenotypic FQ resistance, although country of residence was likely a contributing factor since 69% of these individuals were from a single country. Initial detection of phenotypic resistance and genotypic resistance occurred simultaneously for most; however, phenotypic resistance occurred earlier in isolates harboring mixtures of alleles of very low abundance (<1% of reads), whereas genotypic resistance often occurred earlier for alleles associated with low-level resistance. Understanding factors influencing acquisition and evolution of FQ resistance could reveal strategies for improved treatment success.

Mitigating a COVID-19 Outbreak Among Major League Baseball Players - United States, 2020.

Mass gatherings have been implicated in higher rates of transmission of SARS-CoV-2, the virus that causes coronavirus disease 2019 (COVID-19), and many sporting events have been restricted or canceled to limit disease spread (1). Based on current CDC COVID-19 mitigation recommendations related to events and gatherings (2), Major League Baseball (MLB) developed new health and safety protocols before the July 24 start of the 2020 season. In addition, MLB made the decision that games would be played without spectators. Before a three-game series between teams A and B, the Philadelphia Department of Public Health was notified of a team A player with laboratory-confirmed COVID-19; the player was isolated as recommended (2). During the series and the week after, laboratory-confirmed COVID-19 was diagnosed among 19 additional team A players and staff members and one team B staff member. Throughout their potentially infectious periods, some asymptomatic team A players and coaches, who subsequently received positive SARS-CoV-2 test results, engaged in on-field play with teams B and C. No on-field team B or team C players or staff members subsequently received a clinical diagnosis of COVID-19. Certain MLB health and safety protocols, which include frequent diagnostic testing for rapid case identification, isolation of persons with positive test results, quarantine for close contacts, mask wearing, and social distancing, might have limited COVID-19 transmission between teams.

The detection of Coccidioides from ambient air in Phoenix, Arizona: Evidence of uneven distribution and seasonality.

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.

Coccidioidal Meningitis in New York Traced to Texas by Fungal Genomic Analysis.

A child developed hydrocephalus. Sixteen months later, it was discovered to be a complication of coccidioidal meningitis. The infection's source was uncertain until genomic analysis of the fungal isolate identified its origin to be a visit to Beeville, Texas. Improved national reporting of cases of coccidioidomycosis might reduce diagnostic delays.

Utility of Whole-Genome Sequencing to Ascertain Locally Acquired Cases of Coccidioidomycosis, Washington, USA.

Coccidioidomycosis is an emerging fungal infection in Washington, USA, and the epidemiology of the disease in this state is poorly understood. We used whole-genome sequencing to differentiate locally acquired cases in Washington on the basis of the previously identified phylogeographic population structure of Coccidioides spp. Clinical isolates from coccidioidomycosis cases involving possible Washington soil exposure were included. Of 17 human infections with epidemiologic evidence of possible local acquisition, 4 were likely locally acquired infections and 13 were likely acquired outside Washington. Isolates from locally acquired cases clustered within the previously established Washington clade of C. immitis. Genetic differences among these strains suggest multiple environmental reservoirs of C. immitis in the state.

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

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

Zoonotic Source Attribution of Salmonella enterica Serotype Typhimurium Using Genomic Surveillance Data, United States.

Increasingly, routine surveillance and monitoring of foodborne pathogens using whole-genome sequencing is creating opportunities to study foodborne illness epidemiology beyond routine outbreak investigations and case-control studies. Using a global phylogeny of Salmonella enterica serotype Typhimurium, we found that major livestock sources of the pathogen in the United States can be predicted through whole-genome sequencing data. Relatively steady rates of sequence divergence in livestock lineages enabled the inference of their recent origins. Elevated accumulation of lineage-specific pseudogenes after divergence from generalist populations and possible metabolic acclimation in a representative swine isolate indicates possible emergence of host adaptation. We developed and retrospectively applied a machine learning Random Forest classifier for genomic source prediction of Salmonella Typhimurium that correctly attributed 7 of 8 major zoonotic outbreaks in the United States during 1998-2013. We further identified 50 key genetic features that were sufficient for robust livestock source prediction.

Whole-genome and targeted sequencing of drug-resistant Mycobacterium tuberculosis on the iSeq100 and MiSeq: A performance, ease-of-use, and cost evaluation.

BACKGROUND: Accurate, comprehensive, and timely detection of drug-resistant tuberculosis (TB) is essential to inform patient treatment and enable public health surveillance. This is crucial for effective control of TB globally. Whole-genome sequencing (WGS) and targeted next-generation sequencing (NGS) approaches have potential as rapid in vitro diagnostics (IVDs), but the complexity of workflows, interpretation of results, high costs, and vulnerability of instrumentation have been barriers to broad uptake outside of reference laboratories, especially in low- and middle-income countries. A new, solid-state, tabletop sequencing instrument, Illumina iSeq100, has the potential to decentralize NGS for individual patient care. METHODS AND FINDINGS: In this study, we evaluated WGS and targeted NGS for TB on both the new iSeq100 and the widely used MiSeq (both manufactured by Illumina) and compared sequencing performance, costs, and usability. We utilized DNA libraries produced from Mycobacterium tuberculosis clinical isolates for the evaluation. We conducted WGS on three strains and observed equivalent uniform genome coverage with both platforms and found the depth of coverage obtained was consistent with the expected data output. Utilizing the standardized, cloud-based ReSeqTB bioinformatics pipeline for variant analysis, we found the two platforms to have 94.0% (CI 93.1%-94.8%) agreement, in comparison to 97.6% (CI 97%-98.1%) agreement for the same libraries on two MiSeq instruments. For the targeted NGS approach, 46 M. tuberculosis-specific amplicon libraries had 99.6% (CI 98.0%-99.9%) agreement between the iSeq100 and MiSeq data sets in drug resistance-associated SNPs. The upfront capital costs are almost 5-fold lower for the iSeq100 ($19,900 USD) platform in comparison to the MiSeq ($99,000 USD); however, because of difference in the batching capabilities, the price per sample for WGS was higher on the iSeq100. For WGS of M. tuberculosis at the minimum depth of coverage of 30x, the cost per sample on the iSeq100 was $69.44 USD versus $28.21 USD on the MiSeq, assuming a 2 × 150 bp run on a v3 kit. In terms of ease of use, the sequencing workflow of iSeq100 has been optimized to only require 27 minutes total of hands-on time pre- and post-run, and the maintenance is simplified by a single-use cartridge-based fluidic system. As these are the first sequencing attempts on the iSeq100 for M. tuberculosis, the sequencing pool loading concentration still needs optimization, which will affect sequencing error and depth of coverage. Additionally, the costs are based on current equipment and reagent costs, which are subject to change. CONCLUSIONS: The iSeq100 instrument is capable of running existing TB WGS and targeted NGS library preparations with comparable accuracy to the MiSeq. The iSeq100 has reduced sequencing workflow hands-on time and is able to deliver sequencing results in <24 hours. Reduced capital and maintenance costs and lower-throughput capabilities also give the iSeq100 an advantage over MiSeq in settings of individualized care but not in high-throughput settings such as reference laboratories, where sample batching can be optimized to minimize cost at the expense of workflow complexity and time.