Volume 30, Number 7—July 2024
Research Letter
Genomic Epidemiology of Large Blastomycosis Outbreak, Ontario, Canada, 2021
Abstract
Using phylogenomic analysis, we provide genomic epidemiology analysis of a large blastomycosis outbreak in Ontario, Canada, caused by Blastomyces gilchristii. The outbreak occurred in a locale where blastomycosis is rarely diagnosed, signaling a possible shift in geographically associated incidence patterns. Results elucidated fungal population genetic structure, enhancing understanding of the outbreak.
North American blastomycosis is an infection most commonly caused by environmental dimorphic fungi Blastomyces dermatitidis and B. gilchristii. Infections range from asymptomatic to severe, typically presenting as respiratory illness, with possible systemic dissemination (1,2). The geographic range of B. dermatitidis and B. gilchristii fungi spans the eastern half of North America, including Ontario, Canada (2,3). Although overall incidence rates are low, isolated cases of blastomycosis are diagnosed regularly among populations in endemic areas (2,3), and clusters and outbreaks occur due to common environmental exposures (1,4–7).
We describe a large genomic epidemiology investigation of a blastomycosis outbreak in Constance Lake First Nation, a small community (population <2,000) in northeastern Ontario, Canada, in a locale where blastomycosis has rarely been encountered (2). We studied samples from 181 patients that were received by the Public Health Ontario Laboratory during November 2021–May 2022. By August 2022, we identified B. gilchristii fungus, by using multilocus sequence typing (8), in cultures from 40 persons linked to the outbreak (37 community residents and 3 persons [deemed travel A, B, and C] who visited the community during the possible exposure window). We observed that most positive cultures (35/40) were derived from specimens collected during a 7-week period—mid-November 2021 through December 2021—and most (39/40) were obtained from respiratory specimens (Table). Patients spanned all age groups; 55% were male and 45% female (Table).
We performed whole genome short-read sequencing (Illumina, https://www.illumina.com) on outbreak isolates and 21 other randomly selected B. gilchristii isolates (Appendix) cultured from specimens of patients from northern and eastern Ontario in July 2019–July 2022 (BioProject no. PRJNA890593). This collection of sequenced isolates contributes substantially to the number of B. gilchristii genomes available to advance research on this important pathogen. We conducted single-nucleotide variant (SNV) analysis to determine genetic diversity and relatedness between the 61 Ontario isolates plus 4 isolates from patients in Minnesota, USA (BioProject no. PRJNA786864), by using MycoSNP v0.22 (https://github.com/CDCgov/mycosnp-nf) with reference genome B. gilchristii SLH14081 (GenBank accession no. GCA_000003855.2) and conducted phylogenomic analysis (neighbor-joining method) using MEGA11 (https://www.megasoftware.net) (1).
Phylogenomic analysis of SNVs suggested that 39 of the 40 outbreak isolates were genetically highly similar, including isolates from 2 nonresidents who had traveled to the impacted community (travel A and B) (Figure). By contrast, the isolate from travel C was genetically dissimilar from the other outbreak isolates and likely represents a sporadic case acquired elsewhere (Figure). High genetic similarity between outbreak isolates and the brief 7-week timeframe for symptom onset for 90% of cases suggests a shared discrete exposure window. The overall analysis depicts multiple genetically distinct populations of B. gilchristii fungi separated by thousands of SNVs, correlated with different geographic regions.
To validate the number of SNV differences between isolates, we separately investigated 9 pairs of biologic replicates (isolates from separate specimens from the same person) and 3 technical replicates from 2 other randomly selected isolates (Appendix). We found that 8 of 9 biologic replicates and all technical replicates averaged 12 (range 5–20) SNVs between corresponding isolates. In contrast, the number of SNVs between 39 outbreak isolates (excluding travel C) averaged 380 (range 96–778); 1 pair of biologic replicates differed by 192 SNVs. By comparison, 5 other clusters of isolates, including 2 clusters from Minnesota patients, had <100 SNVs between isolates. Taken together with the timing of the outbreak, the SNV differences between isolates suggests a shared exposure to a genetically similar yet heterogeneous fungal population, rather than a genetically identical point source. We theorize that the pair of biologic replicates differing by 192 SNVs might represent infection by 2 different strains of the fungal population. Of note, the outbreak locale possesses several environmental niches considered suitable reservoirs for Blastomyces spp. fungi, namely waterway-adjacent and forested areas, with coniferous trees supplying acidic soil rich in decaying organic material (7).
During the past 40 years, reports have documented several outbreaks and case clusters of North American blastomycosis (1,4–7). Most ascribe exposure to a discrete time and environmental locale, although that presumption is challenging to confirm because of the variable incubation time (1–6 months) (6) and patient travel. Recently, genomic epidemiologic information gleaned from whole-genome bioinformatics analysis has enhanced the description of outbreaks and case clusters of fungal pathogens (9,10), including B. gilchristii (1), enabling a more fulsome understanding of disease acquisition. In this study, we have interpreted the genomic epidemiology of a large outbreak of B. gilchristii within the context of several contemporaneous, genetically unrelated isolates to describe pathogen population genetic structure. We consider such genomic epidemiologic analyses useful in identifying cases caused by genetically related strains and, when combined with other epidemiologic factors including patient travel, pinpointing the timing and location of exposure. Epidemiologic investigation combined with pathogen phylogenomic analyses enables improved understanding of blastomycosis outbreaks and disease dynamics. This detailed information might help elucidate ecologic variables (possibly altered by climate change or modified land-use patterns) that influence disease acquisition (3). We believe that increased awareness of pathogen range and risk can aid prompt future case diagnoses.
Dr. McTaggart is a research technician at Public Health Ontario, Toronto, Ontario, Canada. Her research interests include mycology, infectious disease epidemiology, and the application of genomic analysis to pathogens of public health significance.
Acknowledgment
We respectfully acknowledge the devastating impact this outbreak has had on the Constance Lake First Nation community and beyond. We are grateful for the permission granted by the community to share these findings to enhance our collective understanding of blastomycosis to help mitigate illness. We also acknowledge the support from Public Health Ontario staff, Indigenous Services Canada, and Porcupine Health Unit, especially Jo Ann Majerovich and Lianne Catton.
References
- Bagal UR, Ireland M, Gross A, Fischer J, Bentz M, Berkow EL, et al. Molecular epidemiology of Blastomyces gilchristii clusters, Minnesota, USA. Emerg Infect Dis. 2022;28:1924–6. DOIPubMedGoogle Scholar
- Brown EM, McTaggart LR, Dunn D, Pszczolko E, Tsui KG, Morris SK, et al. Epidemiology and geographic distribution of blastomycosis, histoplasmosis, and coccidioidomycosis, Ontario, Canada, 1990–2015. Emerg Infect Dis. 2018;24:1257–66. DOIPubMedGoogle Scholar
- Ashraf N, Kubat RC, Poplin V, Adenis AA, Denning DW, Wright L, et al. Re-drawing the maps for endemic mycoses. Mycopathologia. 2020;185:843–65. DOIPubMedGoogle Scholar
- Roy M, Benedict K, Deak E, Kirby MA, McNiel JT, Sickler CJ, et al. A large community outbreak of blastomycosis in Wisconsin with geographic and ethnic clustering. Clin Infect Dis. 2013;57:655–62. DOIPubMedGoogle Scholar
- Dwight PJ, Naus M, Sarsfield P, Limerick B. An outbreak of human blastomycosis: the epidemiology of blastomycosis in the Kenora catchment region of Ontario, Canada. Can Commun Dis Rep. 2000;26:82–91.PubMedGoogle Scholar
- Klein BS, Vergeront JM, Weeks RJ, Kumar UN, Mathai G, Varkey B, et al. Isolation of Blastomyces dermatitidis in soil associated with a large outbreak of blastomycosis in Wisconsin. N Engl J Med. 1986;314:529–34. DOIPubMedGoogle Scholar
- Pfister JR, Archer JR, Hersil S, Boers T, Reed KD, Meece JK, et al. Non-rural point source blastomycosis outbreak near a yard waste collection site. Clin Med Res. 2011;9:57–65. DOIPubMedGoogle Scholar
- Brown EM, McTaggart LR, Zhang SX, Low DE, Stevens DA, Richardson SE. Phylogenetic analysis reveals a cryptic species Blastomyces gilchristii, sp. nov. within the human pathogenic fungus Blastomyces dermatitidis. PLoS One. 2013;8:
e59237 . DOIPubMedGoogle Scholar - Chow NA, Gade L, Tsay SV, Forsberg K, Greenko JA, Southwick KL, et al.; US Candida auris Investigation Team. Multiple introductions and subsequent transmission of multidrug-resistant Candida auris in the USA: a molecular epidemiological survey. Lancet Infect Dis. 2018;18:1377–84. DOIPubMedGoogle Scholar
- Oltean HN, Springer M, Bowers JR, Barnes R, Reid G, Valentine M, et al. Suspected locally acquired coccidioidomycosis in human, Spokane, Washington, USA. Emerg Infect Dis. 2020;26:606–9. DOIPubMedGoogle Scholar
Figure
Table
Cite This ArticleOriginal Publication Date: June 12, 2024
Table of Contents – Volume 30, Number 7—July 2024
EID Search Options |
---|
|
|
|
Please use the form below to submit correspondence to the authors or contact them at the following address:
Julianne V. Kus, Public Health Ontario, 661 University Ave, Toronto, ON M5G 1M1, Canada
Top