Volume 31, Number 8—August 2025
CME ACTIVITY - Research
Scheffersomyces spartinae Fungemia among Pediatric Patients, Pakistan, 2020–2024
Introduction
In support of improving patient care, this activity has been planned and implemented by Medscape, LLC and Emerging Infectious Diseases. Medscape, LLC is jointly accredited with commendation by the Accreditation Council for Continuing Medical Education (ACCME), the Accreditation Council for Pharmacy Education (ACPE), and the American Nurses Credentialing Center (ANCC), to provide continuing education for the healthcare team.
Medscape, LLC designates this Journal-based CME activity for a maximum of 1.00 AMA PRA Category 1 Credit(s)™. Physicians should claim only the credit commensurate with the extent of their participation in the activity.
Successful completion of this CME activity, which includes participation in the evaluation component, enables the participant to earn up to 1.0 MOC points in the American Board of Internal Medicine's (ABIM) Maintenance of Certification (MOC) program. Participants will earn MOC points equivalent to the amount of CME credits claimed for the activity. It is the CME activity provider's responsibility to submit participant completion information to ACCME for the purpose of granting ABIM MOC credit.
All other clinicians completing this activity will be issued a certificate of participation. To participate in this journal CME activity: (1) review the learning objectives and author disclosures; (2) study the education content; (3) take the post-test with a 75% minimum passing score and complete the evaluation at https://www.medscape.org/qna/processor/74834?showStandAlone=true&src=prt_jcme_eid_mscpedu; and (4) view/print certificate.
NOTE: It is the policy of Medscape Education to avoid the mention of brand names or specific manufacturers in accredited educational activities. However, trade and manufacturer names in this activity are provided in an effort to provide clarity. The use of brand or manufacturer names should not be viewed as an endorsement by Medscape of any specific product or manufacturer.
Upon completion of this activity, participants will be able to:
Assess characteristics of S. spartinae
Distinguish the median age of patients with S. spartinae infection in the current study
Analyze the antifungal sensitivity profile of S. spartinae
Evaluate the genetic signature of S. spartinae isolates in the current study
Amy J. Guinn, BA, MA, Technical Writer/Editor, Emerging Infectious Diseases. Disclosure: Amy J. Guinn, BA, MA, has no relevant financial relationships.
Charles P. Vega, MD, Health Sciences Clinical Professor of Family Medicine, University of California, Irvine School of Medicine, Irvine, California. Disclosure: Charles P. Vega, MD, has the following relevant financial relationships: served as consultant or advisor for Boehringer Ingelheim Pharmaceuticals, Inc.; Exact Sciences Corporation.
Kauser Jabeen, MBBS, MSc; Joveria Farooqi, MBBS, MSc; Lacy M. Simons, BS; Judd F. Hultquist, PhD; Ramon Lorenzo-Redondo, PhD ; Charlesnika T. Evans, PhD, MPH; Erica M. Hartmann, PhD; Mohammad Hanif, MBBS; Zahra Hasan, PhD; Syed Faisal Mahmood, MBBS; Javaria Ashraf, PhD; Hassan Ghayas, PhD; Sadaf Zaka, MPhil; Noureen Saeed, MSc; Sayed Ali Raza Shah Bukhari, PhD; Mehreen Arshad, MBBS; Larry K. Kociolek, MD, MSCI; Sameer J. Patel, MD; Rumina Hasan, MBBS, PhD; Egon A. Ozer, MD, PhD.
Abstract
Prevalence of emerging fungal infections is increasing, particularly among immunocompromised persons, children, and older persons. We report 108 cases of Scheffersomyces spartinae infection in pediatric patients from Karachi and other cities in Pakistan, of which 107 were identified from blood cultures. Cultures were initially misidentified as Clavispora lusitaniae by a biochemical assay before speciation as S. spartinae by whole-genome sequencing. All isolates were from children <12 years of age, and >69% were from children <1 month of age. Isolates were genetically distinct across regions of Pakistan; however, genetic diversity was low in isolates from patients in Karachi and nearby Nawabshah and had median differences of just 9 pairwise nucleotide variants. This study demonstrates S. spartinae is a potentially emerging pathogen in neonates and young infants in Pakistan. The findings highlight the limitations of phenotypic identification for detecting emerging fungal infections and underscore the value of molecular identification approaches.
Globally, fungi continue to emerge as serious threats to human health, marked by an increase in incidence of invasive fungal infections and deaths, the growing burden of cutaneous and allergic fungal diseases, and the rapid emergence of antifungal resistance in fungi pathogenic to humans (1). The number of fungi associated with invasive disease is rising as new pathogenic fungi emerge and modern sequencing and analysis approaches drive insights into fungal taxonomy (2,3). Emergence of new pathogenic species and increased infections caused by previously rare or region-specific pathogens might be attributed in various cases to climate change (4–7), more frequent or severe natural disasters (8,9), expanding host populations (10–12), increased use of antifungal agents in clinical or industrial practice, or other unknown epidemiologic drivers.
Considering the immense diversity of fungi and the emergence of new pathogens, ongoing epidemiologic and molecular surveillance of fungal infections is essential. Current phenotypic and biochemical methods do not reliably identify rare or emerging fungal pathogens, and misclassification can occur (13–15). Molecular approaches, such as internal transcribed spacer (ITS) sequencing, and genomic approaches, such as whole-genome sequencing (WGS), can provide greater specificity for distinguishing new or uncommon fungal infections; however, those methods might not be available to clinical microbiology laboratories in some resource-limited settings (16).
In late 2022, the Aga Khan University Hospital (AKUH) clinical mycology laboratory in Karachi, Pakistan, noticed increased cases of yeast, biochemically identified as Clavispora lusitaniae, cultured from the blood of infants admitted to neonatal intensive care units. Most isolates were subsequently identified by WGS as Scheffersomyces spartinae, an environmental yeast with high tolerance to temperature, pH, and salinity changes (17,18). S. spartinae has been isolated from brackish waters, marshes, and oceans and can survive temperatures of 5°C–35°C, pH of 3–9, and salinity of 0–60 g/L NaCl (19,20). We conducted genomic analyses of available isolates to identify epidemiologic characteristics of S. spartinae infections in Pakistan.
Study Setting and Duration
The study was performed at AKUH during 2020–2024. The AKUH clinical microbiology laboratory received specimens from 315 satellite collection units in 99 cities and towns across all 4 provinces of Pakistan. Specimens were sent at the primary physician’s request; the laboratory had no control over specimen ordering.
Isolate Collection
The AKUH clinical microbiology laboratory followed a standardized approach for processing fungal cultures and selecting specimens for archiving. The laboratory received and cultured blood samples in a BacT/ALERT continuous monitoring system (bioMérieux) for 5 days. If samples were flagged positive, results of Gram stain from the broth determined the choice of solid media for subculture. Specimens with yeast on Gram stain were inoculated on BD BBL sheep blood agar, BD BBL Saboraud’s dextrose agar, and BD BBL CHROMagar chromogenic agar plates (all from Becton, Dickinson and Company) and incubated at 35°C–37°C. In addition, the laboratory staff struck yeast on BiGGY agar and Corn Meal Agar with Tween 80 (all Oxoid) for morphological identification, then performed biochemical identification by using the VITEK 2 YST ID card system (bioMérieux).
Antifungal susceptibility was determined by using the Sensititre YeastONE system (Thermo Fisher Scientific). Laboratory staff regularly reviewed sample processing records and selected isolates resistant to fluconazole, echinocandins, or amphotericin B and uncommon or rare yeast isolates for banking. Selected isolates were recultured on potato dextrose agar and incubated at 37°C for 48 hours. Staff transferred colonies from pure cultures to vials of 50% glycerol phosphate broth and banked isolates at −80°C.
Retrospective Data Review
We reviewed the laboratory database to identify all cases of C. lusitaniae in blood cultures from 2015–2024. We retrieved patient age, sex, location of residence, and contact details from the clinical laboratory’s integrated laboratory information system. For most cases, the contact details were for patients or their guardians and not for the treating physicians. We obtained informed consent from the patients or their guardians by telephone in the presence of a witness before obtaining medical histories. We recorded information in a database on REDCap. The AKUH ethical review committee provided ethical review and approval (study nos. 2022-6798-20372 and 2019-0438-2659).
Whole-Genome Sequencing
We inoculated banked isolates in sterile brain heart infusion broth for overnight incubation at 37°C and subcultured on Sabouraud’s dextrose agar for DNA extraction. We performed cell lysis from colonies by using glass beads in lysis buffer and proteinase K and 2 cycles of vortexing at 3,000 rpm for 1 minute each. We then incubated cells at 56°C for 15 minutes. We used the QIAquick DNA Minikit (QIAGEN) to extract genomic DNA from the lysate using the default protocol, eluted in AE buffer, and either stored at −20°C for <2 weeks before sequencing or stored at −80°C if sequencing was not performed within 2 weeks. Sequencing was performed at Aga Khan University (AKU) or Northwestern University (NU; Chicago, Illinois, USA). We used the Speed Vac vacuum concentrator (Thermo Fisher Scientific) to dry extracted DNA specimens before shipping to NU at ambient temperature. NU performed WGS by using the NovaSeq X or MiSeq platforms (Illumina) after library preparation with the plexWell Library Prep Kit (SeqWell). AKU sequenced a subset of the specimens on the Illumina MiniSeq platform after library preparation using the Nextera XT Kit (Illumina).
We trimmed Illumina read sequences by using fastp version 0.23.2 and performed de novo assembly by using SPAdes version 3.9.1. We removed contigs <200 bp or that had <5× mean read coverage. We used BUSCO and the saccharomyces_odb10 database (21) to assess completeness of genome assembly. We performed species identification by extracting ITS sequences from draft assemblies by using in silico PCR and previously described primer sequences ITS5 and ITS4 (22). We compared resulting sequences against the National Center for Biotechnology Information (NCBI) Fungal ITS RNA sequence database (BioProject accession no. PRJNA177353) using BLAST (23). We generated ITS sequence alignments by using MAFFT (24), and generated neighbor-joining trees by using RapidNJ (25).
Long-Read Sequencing
We randomly selected 1 isolate, 128-CS, for long-read sequencing on the MinION sequencer at AKU by using a Flongle Flow Cell R10.4.1 after genomic DNA library preparation using the Ligation Sequencing Kit V14 (all Oxford Nanopore). We used the 400-bp super accuracy model for basecalling, then performed genome assembly by using Trycycler version 0.5.0 (26) as follows. We filtered raw nanopore reads by using Filtlong version 0.2.1 to remove reads of <1,000 bases and the 5% of reads with the lowest quality. We used Trycycler to subsample reads into 12 subsets, then separately reassembled read subsets by using Raven version 1.5.3 (27). We used Trycycler to perform clustering, reconciliation, circularization, multiple sequence alignment, and consensus sequence generation from the 12 assemblies. We polished the resulting consensus assembly with nanopore long reads and medaka version 1.4.3 using the r1041_e82_400bps_sup_v4.1.0 model. We aligned Illumina reads to the assembly by using BWA version 0.7.15 (28) and corrected assembly errors by using Polypolish version 0.5.0 (29) and the POLCA module of MaSuRCA version 4.0.9 (30,31). We deposited sequences in the NCBI BioProject database (accession no. PRJNA1164417) (Appendix Table 1).
Single-Nucleotide Variant Identification and Phylogenetic Analysis
We aligned Illumina reads from each isolate to the genome sequence of strain 128-CS by using BWA version 0.7.15 (28). We used bcftools version 1.9 and a haploid model to identify single-nucleotide variants (SNVs) and skipped bases with quality <25 or alignment quality <30. We further filtered SNVs by using the bcftools_filter software, as previously described (32), to remove variants with SNV quality scores of <200, read consensus <75%, read depth <5, read numbers of <1 in each direction, or locations within repetitive regions as defined by BLAST alignment of the reference genome sequence against itself. We generated a maximum-likelihood phylogenetic tree from the core genome alignment by using IQ-TREE version 2.2.0 and the ModelFinder function to estimate the best fit nucleotide substitution model by means of Bayesian information criterion (33,34). We used the ultrafast booststrap method (35) with 1,000 replicates to assess tree topology. We used phylo-treetime 0.11.4 to estimate time-scaled phylogeny by using the previously inferred maximum-likelihood phylogeny and incorporating specimen sampling dates using the covariation and stochastic-resolve options, an autocorrelated molecular clock with the relax 1.0 0.5 option, and then rerooted the tree by using the least-squares method (36). We visualized and annotated the tree in R version 4.2.2 (The R Project for Statistical Computing) by using the ggtree version 3.8.0 and ggtreeExtra version 1.8.1 packages (37,38).
During 2015–2023, the AKUH clinical microbiology laboratory identified 432 presumed Clavispora (formerly Candida) lusitaniae fungemia cases. Presumed C. lusitaniae cases detected from blood cultures increased from 10–20 cases/year during 2015–2019 to 30 cases in 2020 and 56 cases in 2021. The number of cases detected from cultures from blood specimens further increased to 108 in 2022 and to 150 in 2023 (Figure 1, panel A), representing <0.11% of the total annual blood culture specimens processed by the laboratory (Appendix Table 2). Of the 432 C. lusitaniae cases, 415 patients (96.1%) were <13 years (4,745 days) of age and 365 (84.5%) were <1 year (365 days) of age: 326 (75.4%) <3 months of age, 274 (63.4%) <1 month of age, and 68 (15.7%) <7 days of age. Reflective of the primary referral base for the AKUH laboratory, most (n = 369; 85.4%) C. lusitaniae bloodstream cases were from Karachi, and 241 (55.8%) samples were from 1 tertiary care public sector pediatric hospital.
Using the laboratory practices for isolate banking described herein, 136 presumed C. lusitaniae isolates had been banked: 21 of 31 isolates cultured in 2020, 0 of 56 isolates cultured in 2021, 29 of 108 isolates cultured in 2022, 80 of 150 isolates cultured in 2023, and 6 of 15 isolates cultured in January of 2024. To investigate those cases, we performed WGS on banked isolates, 121 of which yielded high-quality assemblies. We excluded 15 isolates because total contig size exceeded 12.3 million bases (n = 8), total contig number exceeded 800 (n = 11), genome coverage was <95% as determined by BUSCO (n = 9), or a combination of those factors. BUSCO analyses showed equivalent rates of assembly completeness, fragmentation, and missingness between the isolates sequenced at either institution. Querying the internal ITS sequences from the whole-genome assemblies against the NCBI ITS database showed that just 13 (11%) of the 121 high-quality genome sequences were true C. lusitaniae, whereas ITS sequences from the other 108 (89%) sequences had 99.3%–100% sequence identity with the ITS sequence of Scheffersomyces spartinae CBS 6059 (NCBI accession no. NR_111290.1) (Figure 1, panel B; Appendix Table 1).
The median age of S. spartinae patients was 19.5 days (range 2 days–12 years); 68.5% of cases were among infants <4 weeks of age (Table 1). For most cases, specimens were obtained in Karachi, reflective of catchment, but >8% of cases were identified in patients from other regions. Most (69.4%) cases were identified at a single hospital (hospital H) in Karachi specializing in neonatal and pediatric care. Blood cultures from that hospital constituted 1.4% of the total blood cultures processed by the AKUH clinical microbiology laboratory since 2020, and 2.9% of blood cultures from that hospital were C. lusitaniae positive by biochemical assay (Appendix Table 2).
We cultured S. spartinae from specimens received from 14 different institutions but could not determine the origin of 11% of specimens (Table 1). Of patients whose families or guardians could be contacted and consented to provide clinical information (92/108), all were hospitalized at the time of specimen collection. Treating physicians provided admitting diagnoses and brief clinical histories for 63 (58.3%) patients, 43 (68.2%) of whom were admitted for sepsis. MICs for amphotericin, azole, and echinocandin antifungal drugs were low among the 108 S. spartinae isolates tested, and MIC90 values were <1 µg/mL for all agents in those classes (Table 2).
Isolate 128-CS, which we selected for long-read sequencing, came from a sample collected in Karachi in December 2022. Assembly from long and short reads of 128-CS produced 8 linear chromosomes ranging in length from 1.19 to 2.05 Mbp and 1 circular mitochondrial genome, for a total genome size of 12.27 Mbp (Appendix Table 4). Comparison of the ITS sequence of isolate 128-CS to the NCBI ITS database reconfirmed S. spartinae (Figure 2).
Phylogenetic analysis of all 108 S. spartinae sequences aligned to 128-CS revealed that most sequences belonged to a clade consisting of nearly all specimens obtained from Karachi and 2 specimens from Nawabshah (Figure 3). That large clade (clade A) contained 2 subclades, clades A-1 and A-2 (Figure 4). We also identified 2 smaller clades: clade B, consisting of 1 isolate each from Karachi and Hafizabad and all 5 isolates from Lahore, and clade C, consisting of 1 genetically distant sequence isolated from a patient in Multan (Figure 3). Although the Multan isolate differed by >213,000 SNVs from 128-CS, the ITS regions of the 2 isolates were 100% identical. All isolates in clade B were collected from infants <1 month of age; the clade C isolate was collected from a 6-month-old child (Figure 3, panel A).
Because most isolates cultured during 2020–2024 were closely related, we examined the characteristics of those 100 specimens. Within that clade, the median pairwise genetic distance across the 12.2 million base alignment was just 9 (range 0–22) SNVs. The time-scaled maximum-likelihood phylogenetic tree showed that, despite the small genetic distances between sequences, subclades A-1 and A-2 are separated with >80% bootstrap support (Figure 4). Isolates from 2020 were all in clade A-2 and clustered more closely with each other than with most later isolates; we excluded 1 isolate, 2020-498, which did not follow the fitted molecular clock. All 5 isolates from hospital B were in the same subclade; otherwise, we noted no apparent association between sampling location and either isolation date or genetic similarity.
Isolates from patients <1 month of age predominated in subclade A-2, comprising 93.6% of that clade compared with 38.5% of clade A-1 (p<0.001 by Fisher exact test) (Table 3). No association between patient age and genetic similarity was noted among isolates in either subclade. Only isolates from hospital H were represented in both subclades and were equally distributed (Table 3). Those results suggest that genetic distances are small among S. spartinae in that group; that isolates were diverse across time, geography, and patient ages; and that closely related but independent lineages were cocirculating within the same region and institution, potentially indicating separate sources of infection.
We describe S. spartinae as a cause of fungemia in pediatric patients in Pakistan. S. spartinae has been previously described as an environmental organism capable of surviving extreme conditions, but isolation from clinical specimens represents emergence as a human pathogen. The genus Scheffersomyces was proposed by Kurtzman and Suzuki on the basis of D1-D2 large subunit rRNA and small subunit rRNA sequencing; species assigned to the Scheffersomyces genus were originally included in the genus Pichia (39). Most Scheffersomyces species ferment xylose and are used in industrial applications (39–41). The assignment of S. spartinae to the genus was originally considered uncertain because of weak bootstrap support in the D1–D2 small subunit tree and because S. spartinae does not ferment xylose. According to our literature review, S. spartinae has not previously been reported as a cause of human infections nor cultured in healthcare facilities. We speculate S. spartinae might be emerging in humans because of introduction from environmental sources into healthcare environments. However, given the broad distribution of healthcare facilities from which we cultured S. spartinae in Pakistan, we cannot rule out nonhealthcare community sources of infection.
The reasons for the emergence of some fungi as agents of human infection have been multifactorial. Fungi often survive extreme conditions, adapt to selection pressures, and develop enhanced thermotolerance, virulence, and antifungal drug resistance. Emergence of Candida auris in human infections has been attributed to global warming (6,7); some researchers postulate that C. auris inhabits an environmental niche and has only recently become a human pathogen (4). Low virulence fungi such as Saccharomyces cerevisiae, Saprochaeta clavate, and Rhodotorula spp. can contaminate food, medical products, and healthcare environments and cause infections in susceptible hosts (10–12). Alterations in the geographic range of Cryptococcus deuterogattii from tropical regions to temperate climates have led to human infections in nonendemic areas, attributed to human activity, thermal adaptation, and climate change (5).
Another reason for increased emergence of fungal infections is environmental disruption caused by natural disasters, which can lead to a wider distribution of fungi (9). An example is increased human Coccidioides immitis infections after earthquake-related landslides and sandstorms in previously low-prevalence areas (8). Whether similar or other factors contributed to the emergence of S. spartinae as a human pathogen is unclear.
Cases of C. lusitaniae fungemia identified by the AKUH laboratory began increasing in 2020, and a further substantial rise in cases began in 2022. Most C. lusitaniae patients were children <1 year of age who were admitted to hospitals in the city of Karachi. During 2015–2019, rare Candida species were reported as causes of invasive infections in neonatal and pediatric populations in Pakistan (42). Isolation of C. lusitaniae was infrequent, however, so rising case numbers from 2022 onwards were flagged by the laboratory, which more frequently banked cultured isolates. Biochemical identification of Candida by the VITEK 2 YST ID card system inaccurately identified the organism, and most of the isolates were ultimately identified as S. spartinae by whole-genome and ITS sequencing. The high percentage of isolates subsequently identified as S. spartinae suggests that at least some of the isolates identified as C. lusitaniae from before 2020 also might have been S. spartinae, but we cannot verify that hypothesis because earlier specimens were not banked.
Species misidentification is concerning and implies that phenotypic identification systems cannot reliably identify emerging pathogens. Misidentification underscores the need for using molecular approaches, such as ITS sequencing, for species identification. Decreasing costs and increasing availability of WGS promise improved species identification as well as strain typing and genomic epidemiology, but the technology remains difficult to access in many low-resource areas.
Although evidence implies that cases of S. spartinae infection represent the emergence of a new pathogen, the possibility of a pseudo-outbreak (i.e., clustering of positive cultures because of contamination of clinical specimens incorrectly attributed to infections) must also be considered. Compared with bacterial pseudo-outbreaks, fungal pseudo-outbreaks have rarely been reported. Documented pseudo-outbreaks involving Candida guilliermondii were related to inadequate sterile technique during blood culture collection (43), contamination of heparin solution used to flush needles before blood draws (44), and an anaerobic holding chamber in the clinical microbiology laboratory (45). Pseudo-outbreaks reported with other Candida species include C. versatilis contamination in olive oil used for culture media supplementation (46) and C. parapsilosis contamination of a salt solution used for grinding tissues (47).
Unlike prior pseudo-outbreak reports that involved temporally clustered samples from single institutions, we report isolates collected over the course of 4 years from >14 healthcare institutions across Pakistan. Furthermore, we noted considerable demographic similarity between patients from which the positive specimens were collected: >68% were <1 month of age, and >80% were <1 year of age (Table 1). Although all the isolates were cultured and identified by a single clinical microbiology laboratory at AKUH, all processes and equipment used for processing clinical specimens in the laboratory were independent of the age of the patient from which the specimen was collected. If contamination in the laboratory explained most or all observed cases, we would expect detection in cultures from patients across all age ranges and a larger number of specimens. Even when C. lusitaniae cases increased during 2022 and 2023, those cases constituted only 0.57% of 16,505 blood cultures processed from children <1 year of age in 2022 and only 0.60% of 18,051 cultures from children in that age group in 2023 (Appendix Table 2).
The genetic diversity of specimens across and within geographic locations (Figures 3, 4) reduces the likelihood of a single infection or contamination source. Although genetic diversity between isolates within a geographic area of sampling was low, that diversity was not inconsistent with the degree of genetic variation among geographically clustered infections caused by other yeast species, such as C. auris (48,49). Until we learn more about the genomic variability and evolution of S. spartinae, the degree of expected genetic difference among epidemiologically linked or unlinked specimens will remain unclear; however, cocirculation of at least 2 lineages in Karachi over the course of >4 years suggests multiple introductions into the region.
One limitation of this study is that, before 2022, AKUH prioritized its finite laboratory resources for investigating and banking invasive antimicrobial-resistant isolates and isolates from known pathogens of epidemiologic concern. As a result, isolates identified as C. lusitaniae were not consistently banked, and few isolates from before 2022 were available for species verification and genomic analysis. Another limitation is that detailed patient clinical and demographic information was either unavailable or difficult to obtain because of variability in record-keeping and data availability across the many institutions that provided specimens to AKUH. Many treating physicians could not be contacted because in Pakistan laboratory tests are paid out of pocket by the patients or their guardians, who often provide their own contact details rather than their physicians’ contact information at sample collection. That practice greatly hinders epidemiologic investigations and transmission reconstructions to determine common exposures or shared risk factors among affected persons, to establish patients’ clinical manifestations, and to assess responses to treatments or outcomes. Efforts are ongoing to establish cooperative agreements and study protocols for the most affected institutions in Pakistan for standardized data collection and sharing for S. spartinae fungemia cases or other unusual yeasts.
In summary, our findings underscore the value of epidemiologic monitoring for identifying infection clusters and of genomic and molecular surveillance for identifying rare and emerging pathogens. Future studies will be directed at characterizing mediators of pathogenicity and virulence factors that could contribute to S. spartinae emergence as a human pathogen, as well as exploring potential environmental reservoirs or other sources of infection. Nonetheless, this study demonstrates that building capacity for specimen identification and banking, WGS, and bioinformatic analysis in low- and middle-income countries like Pakistan is imperative for early detection and study of emerging infectious disease threats.
Author contributions: conceptualization: K.J., J.F., R.H., E.A.O.; methodology: K.J., J.F., L.M.S., R.L.R., J.A., H.G., S.A.R.S.B., R.H., E.A.O.; software: R.L.R., E.A.O.; validation: K.J., J.F., L.M.S., N.S., E.A.O.; formal analysis: K.J., J.F., R.L.R., E.A.O.; investigation: K.J., J.F., L.M.S., J.F.H., R.L.R., C.T.E., E.M.H., M.H., S.F.M., J.A., H.G., S.Z., N.S., S.A.R.S.B., M.A., L.K.K., S.J.P., R.H., E.A.O.; resources: K.J., J.F., J.F.H., M.H., Z.H., R.H., E.A.O.; data curation: K.J., L.M.S., N.S., R.H., E.A.O.; writing original draft: K.J., E.A.O.; writing, review, and editing: K.J., J.F., L.M.S., J.F.H., R.L.R., C.T.E., E.M.H., Z.H., S.F.M., N.S., S.A.R.S.B., M.A., L.K.K., S.J.P., R.H., E.A.O.; visualization: K.J., E.A.O.; supervision: K.J., J.F., J.F.H., Z.H., L.K.K., R.H., E.A.O.; project administration: L.M.S., J.F.H., Z.H., R.H., E.A.O.; funding acquisition: J.F.H., R.H., E.A.O.
Acknowledgments
We thank Safia Bibi for her contribution.
This study was supported by the US Centers for Disease Control and Prevention (award no. NU3HCK000007) and by the US National Institutes of Health (award no. S10OD032243). Support was also provided by the Antimicrobial Resistance Working Group in the Northwestern University Buffett Institute for Global Affairs. This work was supported in part through the computational resources and staff contributions provided by Genomics Compute Cluster, which is jointly supported by the Feinberg School of Medicine, the Center for Genetic Medicine, and Feinberg’s Department of Biochemistry and Molecular Genetics, the Office of the Provost, the Office for Research, and Northwestern Information Technology. The Genomics Compute Cluster is part of Quest, Northwestern University’s high-performance computing facility, with the purpose to advance research in genomics. The funding sources had no role in the study design, data collection, analysis, interpretation, or writing of the report.
References
- Fisher MC, Gurr SJ, Cuomo CA, Blehert DS, Jin H, Stukenbrock EH, et al. Threats posed by the fungal kingdom to humans, wildlife, and agriculture. MBio. 2020;11:e00449–20. DOIPubMedGoogle Scholar
- Wiederhold NP. Emerging fungal infections: new species, new names, and antifungal resistance. Clin Chem. 2021;68:83–90. DOIPubMedGoogle Scholar
- World Health Organization. WHO fungal priority pathogens list to guide research, development and public health action. Geneva: The Organization; 2022.
- Casadevall A, Kontoyiannis DP, Robert V. On the emergence of Candida auris: climate change, azoles, swamps, and birds. MBio. 2019;10:e01397–19. DOIPubMedGoogle Scholar
- Bartlett KH, Kidd SE, Kronstad JW. The emergence of Cryptococcus gattii in British Columbia and the Pacific Northwest. Curr Infect Dis Rep. 2008;10:58–65. DOIPubMedGoogle Scholar
- Chakrabarti A, Sood P. On the emergence, spread and resistance of Candida auris: host, pathogen and environmental tipping points. J Med Microbiol. 2021;70:
001318 . DOIPubMedGoogle Scholar - Du H, Bing J, Hu T, Ennis CL, Nobile CJ, Huang G. Candida auris: Epidemiology, biology, antifungal resistance, and virulence. PLoS Pathog. 2020;16:
e1008921 . DOIPubMedGoogle Scholar - Schneider E, Hajjeh RA, Spiegel RA, Jibson RW, Harp EL, Marshall GA, et al. A coccidioidomycosis outbreak following the Northridge, Calif, earthquake. JAMA. 1997;277:904–8. DOIPubMedGoogle Scholar
- Seidel D, Wurster S, Jenks JD, Sati H, Gangneux JP, Egger M, et al. Impact of climate change and natural disasters on fungal infections. Lancet Microbe. 2024;5:e594–605. DOIPubMedGoogle Scholar
- Vaux S, Criscuolo A, Desnos-Ollivier M, Diancourt L, Tarnaud C, Vandenbogaert M, et al.; Geotrichum Investigation Group. Multicenter outbreak of infections by Saprochaete clavata, an unrecognized opportunistic fungal pathogen. MBio. 2014;5:e02309–14. DOIPubMedGoogle Scholar
- Wirth F, Goldani LZ. Epidemiology of Rhodotorula: an emerging pathogen. Interdiscip Perspect Infect Dis. 2012;2012:
465717 . DOIPubMedGoogle Scholar - Muñoz P, Bouza E, Cuenca-Estrella M, Eiros JM, Pérez MJ, Sánchez-Somolinos M, et al. Saccharomyces cerevisiae fungemia: an emerging infectious disease. Clin Infect Dis. 2005;40:1625–34. DOIPubMedGoogle Scholar
- Wickes BL, Wiederhold NP. Molecular diagnostics in medical mycology. Nat Commun. 2018;9:5135. DOIPubMedGoogle Scholar
- Lücking R, Aime MC, Robbertse B, Miller AN, Ariyawansa HA, Aoki T, et al. Unambiguous identification of fungi: where do we stand and how accurate and precise is fungal DNA barcoding? IMA Fungus. 2020;11:14. DOIPubMedGoogle Scholar
- Steenwyk JL, Balamurugan C, Raja HA, Gonçalves C, Li N, Martin F, et al. Phylogenomics reveals extensive misidentification of fungal strains from the genus Aspergillus. Microbiol Spectr. 2024;12:
e0398023 . DOIPubMedGoogle Scholar - Chen SC, Perfect J, Colombo AL, Cornely OA, Groll AH, Seidel D, et al. Global guideline for the diagnosis and management of rare yeast infections: an initiative of the ECMM in cooperation with ISHAM and ASM. Lancet Infect Dis. 2021;21:e375–86. DOIPubMedGoogle Scholar
- Breyer E, Espada-Hinojosa S, Reitbauer M, Karunarathna SC, Baltar F. Physiological properties of three pelagic fungi isolated from the Atlantic Ocean. J Fungi (Basel). 2023;9:439. DOIPubMedGoogle Scholar
- Urano N, Shirao A, Naito Y, Okai M, Ishida M, Takashio M. Molecular phylogeny and phenotypic characterization of yeasts with a broad range of pH tolerance isolated from natural aquatic environments. Adv Microbiol. 2019;9:56–73. DOIGoogle Scholar
- Hagler AN, Mendonça-Hagler LC, Pagnocca FC. Yeasts in aquatic ecotone habitats. In: Buzzini P, Lachance M-A, Yurkov A, editors. Yeasts in natural ecosystems: diversity. Cham: Springer International Publishing; 2017. p. 63–85.
- Frigon MD, Liu D. Effect of high salinity on yeast activated sludge reactor operation. Water Sci Technol. 2016;74:2124–34. DOIPubMedGoogle Scholar
- Manni M, Berkeley MR, Seppey M, Simão FA, Zdobnov EM. BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol Biol Evol. 2021;38:4647–54. DOIPubMedGoogle Scholar
- White T, Bruns T, Lee S, Taylor J, Innis M, Gelfand D, et al. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, editors. PCR protocols: a guide to methods and applications, volume 18. San Diego: Academic Press; 1990. p. 315–22.
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10. DOIPubMedGoogle Scholar
- Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80. DOIPubMedGoogle Scholar
- Simonsen M, Mailund T, Pedersen CNS. Rapid neighbour-joining. In: Cradall KA, Lagergren J, editors. WABI ’08: Proceedings of the 8th international workshop on algorithms in bioinformatics. Berlin: Springer Berlin Heidelberg; 2008. p. 113–22.
- Wick RR, Judd LM, Cerdeira LT, Hawkey J, Méric G, Vezina B, et al. Trycycler: consensus long-read assemblies for bacterial genomes. Genome Biol. 2021;22:266. DOIPubMedGoogle Scholar
- Vaser R, Šikić M. Time- and memory-efficient genome assembly with Raven. Nat Comput Sci. 2021;1:332–6. DOIPubMedGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60. DOIPubMedGoogle Scholar
- Wick RR, Holt KE. Polypolish: Short-read polishing of long-read bacterial genome assemblies. PLOS Comput Biol. 2022;18:
e1009802 . DOIPubMedGoogle Scholar - Zimin AV, Salzberg SL. The genome polishing tool POLCA makes fast and accurate corrections in genome assemblies. PLOS Comput Biol. 2020;16:
e1007981 . DOIPubMedGoogle Scholar - Zimin AV, Marçais G, Puiu D, Roberts M, Salzberg SL, Yorke JA. The MaSuRCA genome assembler. Bioinformatics. 2013;29:2669–77. DOIPubMedGoogle Scholar
- Bulman ZP, Krapp F, Pincus NB, Wenzler E, Murphy KR, Qi C, et al. Genomic features associated with the degree of phenotypic resistance to carbapenems in carbapenem-resistant Klebsiella pneumoniae. mSystems. 2021;6:
e0019421 . DOIPubMedGoogle Scholar - Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9. DOIPubMedGoogle Scholar
- Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37:1530–4. DOIPubMedGoogle Scholar
- Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35:518–22. DOIPubMedGoogle Scholar
- Sagulenko P, Puller V, Neher RA. TreeTime: Maximum-likelihood phylodynamic analysis. Virus Evol. 2018;4:
vex042 . DOIPubMedGoogle Scholar - Yu G. Using ggtree to visualize data on tree-like structures. Curr Protoc Bioinformatics. 2020;69:
e96 . DOIPubMedGoogle Scholar - Xu S, Dai Z, Guo P, Fu X, Liu S, Zhou L, et al. ggtreeExtra: compact visualization of richly annotated phylogenetic data. Mol Biol Evol. 2021;38:4039–42. DOIPubMedGoogle Scholar
- Kurtzman CP, Suzuki M. Phylogenetic analysis of ascomycete yeasts that form coenzyme Q-9 and the proposal of the new genera Babjeviella, Meyerozyma, Millerozyma, Priceomyces, and Scheffersomyces. Mycoscience. 2010;51:2–14. DOIGoogle Scholar
- Kurtzman CP, Suzuki M. Chapter 65: Scheffersomyces. In: Kurtzman CP, Fell JW, Boekhout T, editors. The yeasts, fifth edition. London: Elsevier; 2011. p. 773–7.
- Tan L, He M, Song L, Fu X, Shi S. Aerobic decolorization, degradation and detoxification of azo dyes by a newly isolated salt-tolerant yeast Scheffersomyces spartinae TLHS-SF1. Bioresour Technol. 2016;203:287–94. DOIPubMedGoogle Scholar
- Rattani S, Memon S, Jabeen K, Farooqi J. Spectrum of invasive yeast infections in neonates, children and adults in Pakistan over five years: 2015–2019. Infect Dis J Pakistan. 2021;30:14–8.
- Medeiros EA, Lott TJ, Colombo AL, Godoy P, Coutinho AP, Braga MS, et al. Evidence for a pseudo-outbreak of Candida guilliermondii fungemia in a university hospital in Brazil. J Clin Microbiol. 2007;45:942–7. DOIPubMedGoogle Scholar
- Yagupsky P, Dagan R, Chipman M, Goldschmied-Reouven A, Zmora E, Karplus M. Pseudooutbreak of Candida guilliermondii fungemia in a neonatal intensive care unit. Pediatr Infect Dis J. 1991;10:928–32.PubMedGoogle Scholar
- Kirby JE, Branch-Elliman W, LaSalvia MT, Longhi L, MacKechnie M, Urman G, et al. Investigation of a Candida guilliermondii pseudo-outbreak reveals a novel source of laboratory contamination. J Clin Microbiol. 2017;55:1080–9. DOIPubMedGoogle Scholar
- Brandt ME, Benjamin LE, Steinkraus GE. Pseudooutbreak of Candida versatilitis fungemia in a microbiology laboratory. Diagn Microbiol Infect Dis. 2003;46:73–5. DOIPubMedGoogle Scholar
- Deresinski SC, Clemons KV, Kemper CA, Roesch K, Walton B, Stevens DA. Genotypic analysis of pseudoepidemic due to contamination of Hanks’ balanced salt solution with Candida parapsilosis. J Clin Microbiol. 1995;33:2224–6. DOIPubMedGoogle Scholar
- Gorzalski A, Ambrosio FJ III, Massic L, Scribner MR, Siao DD, Hua C, et al. The use of whole-genome sequencing and development of bioinformatics to monitor overlapping outbreaks of Candida auris in southern Nevada. Front Public Health. 2023;11:
1198189 . DOIPubMedGoogle Scholar - Roberts SC, Zembower TR, Ozer EA, Qi C. Genetic evaluation of nosocomial Candida auris transmission. J Clin Microbiol. 2021;59:e02252–20. DOIPubMedGoogle Scholar
Figures
Tables
Follow Up
To obtain credit, you should first read the journal article. After reading the article, you should be able to answer the following, related, multiple-choice questions. To complete the questions (with a minimum 75% passing score) and earn continuing medical education (CME) credit, please go to https://www.medscape.org/journal/eid. Credit cannot be obtained for tests completed on paper, although you may use the worksheet below to keep a record of your answers.
You must be a registered user on http://www.medscape.org. If you are not registered on http://www.medscape.org, please click on the “Register” link on the right hand side of the website.
Only one answer is correct for each question. Once you successfully answer all post-test questions, you will be able to view and/or print your certificate. For questions regarding this activity, contact the accredited provider,
Which of the following statements regarding Scheffersomyces spartinae is the most accurate?
It is found in dry soil
It generally will die when temperatures drop below 10 °C
It only thrives in environments with pH below 7
It has not been considered an infectious agent among humans
What was the approximate median age of patients with positive cultures for S. spartinae in the current study?
19 days
4 months
6 years
16 years
S. spartinae demonstrated resistance to which of the following antifungal classes in the current study?
Amphotericin
Azoles
Echinocandins
None of the above
Which of the following statements regarding genetic studies of S. spartinae in the current study is the most accurate?
Only 1 clade was identified
There was a large degree of genetic variability among isolates
Isolates were diverse across time and geographic location
There was no co-circulation of independent lineages in particular locations
Original Publication Date: July 24, 2025
Related Links
Table of Contents – Volume 31, Number 8—August 2025
EID Search Options |
---|
|
|
|
Please use the form below to submit correspondence to the authors or contact them at the following address:
Egon A. Ozer, Northwestern University Feinberg School of Medicine, 303 E Superior St, Ste 900, Chicago, IL 60611, USA
Top