# Outbreak of Invasive Wound Mucormycosis in a Burn Unit Due to Multiple Strains of Mucor circinelloides f. circinelloides Resolved by Whole-Genome Sequencing

Dea Garcia-Hermoso, Alexis Criscuolo, Soo Chan Lee, Matthieu Legrand, Marc Chaouat, Blandine Denis, Matthieu Lafaurie, Martine Rouveau, Charles Soler, Jean-Vivien Schaal, et al.

► To cite this version:

Dea Garcia-Hermoso, Alexis Criscuolo, Soo Chan Lee, Matthieu Legrand, Marc Chaouat, et al. Outbreak of Invasive Wound Mucormycosis in a Burn Unit Due to Multiple Strains of Mucor circinelloides f. circinelloides Resolved by Whole-Genome Sequencing. mBio, 2018, 9 (2), pp.e00573-18. 10.1128/mBio.00573-18. pasteur-01854232

HAL Id: pasteur-01854232

https://pasteur.hal.science/pasteur-01854232

Submitted on 6 Aug 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Outbreak of Invasive Wound Mucormycosis in a Burn Unit Due to Multiple Strains of *Mucor circinelloides* f. *circinelloides* Resolved by Whole-Genome Sequencing

Dea Garcia-Hermoso,ª,ᵇ Alexis Criscuolo,ᶜ Soo Chan Lee,ᵈ Matthieu Legrand,ᵉ,ᵉ Marc Chaouat,ᵉ,ᵉ Blandine Denis,ʰ Matthieu Lafaurie,ʰ Martine Rouveau,<fnref n="1" /> Charles Soler,<fnref n="1" /> Jean-Vivien Schaal,⁵ Maurice Mimoun,ᵉ,ᵉ Alexandre Mebazaa,ᵉ,ᵉ Joseph Heitman,<fnref n="1" /> Françoise Dromer,ª,ᵇ Sylvain Brisse,⁸ Stéphane Bretagne,ª,ᵇ,ᵉ,ⁿ Alexandre Alanioª,ᵇ,ᵉ,ⁿ

**ABSTRACT** Mucorales are ubiquitous environmental molds responsible for mucor-mycosis in diabetic, immunocompromised, and severely burned patients. Small outbreaks of invasive wound mucormycosis (IWM) have already been reported in burn units without extensive microbiological investigations. We faced an outbreak of IWM in our center and investigated the clinical isolates with whole-genome sequencing (WGS) analysis. We analyzed *M. circinelloides* isolates from patients in our burn unit (BU1, Hôpital Saint-Louis, Paris, France) together with nonoutbreak isolates from Burn Unit 2 (BU2, Paris area) and from France over a 2-year period (2013 to 2015). A total of 21 isolates, including 14 isolates from six BU1 patients, were analyzed by whole-genome sequencing (WGS). Phylogenetic classification based on *de novo* assembly and assembly free approaches showed that the clinical isolates clustered in four highly divergent clades. Clade 1 contained at least one of the strains from the six epidemiologically linked BU1 patients. The clinical isolates were specific to each patient. Two patients were infected with more than two strains from different clades, suggesting that an environmental reservoir of clonally unrelated isolates was the source of contamination. Only two patients from BU1 shared one strain, which could correspond to direct transmission or contamination with the same environmental source. In conclusion, WGS of several isolates per patients coupled with precise epidemiological data revealed a complex situation combining potential cross-transmission

**Received** 13 March 2018 **Accepted** 29 March 2018 **Published** 24 April 2018

**Citation** Garcia-Hermoso D, Criscuolo A, Lee SC, Legrand M, Chaouat M, Denis B, Lafaurie M, Rouveau M, Soler C, Schaal J-V, Mimoun M, Mebazaa A, Heitman J, Dromer F, Brisse S, Bretagne S, Alanio A. 2018. Outbreak of invasive wound mucormycosis in a burn unit due to multiple strains of *Mucor circinelloides* f. *circinelloides* resolved by whole-genome sequencing. *mBio* 9:e00573-18. https://doi.org/10.1128/mBio.00573-18.

**Editor** Anuradha Chowdhary, Vallabhbhai Patel Chest Institute

**Copyright** © 2018 Garcia-Hermoso et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Address correspondence to Alexandre Alanio, alexandre.alanio@pasteur.fr.

D.G.-H. and A.C. contributed equally to this article.

This article is a direct contribution from a Fellow of the American Academy of Microbiology. Solicited external reviewers: Shawn Lockhart, Centers for Disease Control and Prevention; Kerstin Voigt, Leibniz Institute for Natural Product and Infection Biology.

<footnote>ªInstitut Pasteur, National Reference Center for Invasive Mycoses and Antifungals (NRCMA), Molecular Mycology Unit, Paris, France</footnote>

<footnote>ᵇCNRS UMR2000, Paris, France</footnote>

<footnote>ᶜInstitut Pasteur-Bioinformatics and Biostatistics Hub-C3BI, USR 3756 IP CNRS-Paris, France</footnote>

<footnote>ᵈSouth Texas Center for Emerging Infectious Diseases (STCEID), Department of Biology, the University of Texas at San Antonio, San Antonio, Texas, USA</footnote>

<footnote>ᵉUMR 942, INSERM, Service d'Anesthésie-Réanimation, Assistance Publique-Hôpitaux de Paris (AP-HP), Groupe Hospitalier Saint-Louis-Lariboisière-Fernand-Widal, Paris, France</footnote>

<footnote>ᵉV Université Paris Diderot, Sorbonne Paris Cité, Paris, France</footnote>

<footnote>ᵉService de Chirurgie Plastique, Reconstructrice, Esthétique et Traitement Chirurgical des Brûlés, Assistance Publique-Hôpitaux de Paris (AP-HP), Groupe Hospitalier Saint-Louis-Lariboisière-Fernand-Widal, Paris, France</footnote>

<footnote>ʰService de Maladies Infectieuses et Tropicales, Assistance Publique-Hôpitaux de Paris (AP-HP), Groupe Hospitalier Saint-Louis-Lariboisière-Fernand-Widal, Paris, France</footnote>

<footnote n="1">¹Laboratoire de Microbiologie, Hôpital Saint-Louis, Assistance Publique-Hôpitaux de Paris (AP-HP), Groupe Hospitalier Saint-Louis-Lariboisière-Fernand-Widal, Paris, France</footnote>

<footnote>JService de Biologie Médicale, Hôpital d'Instruction des Armées, Clamart, France</footnote>

<footnote>kCentre de Traitement des Brûlés, Hôpital d'Instruction des Armées, Clamart, France</footnote>

<footnote n="1">¹Molecular Genetics and Microbiology, Duke University Medical Center, Durham, North Carolina, USA</footnote>

<footnote>mInstitut Pasteur, Biodiversity and Epidemiology of Bacterial Pathogens, Paris, France</footnote>

<footnote>nLaboratoire de Parasitologie-Mycologie, Hôpital Saint-Louis, Assistance Publique-Hôpitaux de Paris (AP-HP), Groupe Hospitalier Saint-Louis-Lariboisière-Fernand-Widal, Paris, France</footnote>
between patients and multiple contaminations with a heterogeneous pool of strains
from a cryptic environmental reservoir.

**IMPORTANCE** Invasive wound mucormycosis (IWM) is a severe infection due to environmental molds belonging to the order Mucorales. Severely burned patients are particularly at risk for IWM. Here, we used whole-genome sequencing (WGS) analysis to resolve an outbreak of IWM due to *Mucor circinelloides* that occurred in our hospital (BU1). We sequenced 21 clinical isolates, including 14 from BU1 and 7 unrelated isolates, and compared them to the reference genome (1006PhL). This analysis revealed that the outbreak was mainly due to multiple strains that seemed patient specific, suggesting that the patients were more likely infected from a pool of diverse strains from the environment rather than from direct transmission among them. This study revealed the complexity of a Mucorales outbreak in the settings of IWM in burn patients, which has been highlighted based on WGS combined with careful sampling.

**KEYWORDS** *Mucor circinelloides*, Mucorales, mucormycosis, strains, burn, clade, genome, genotype, mixed infection, outbreak, whole-genome sequencing, wound

Mucormycosis is a rare and life-threatening infection caused by Mucorales belonging to the subphylum *Mucoromycotina* (1). These molds are ubiquitously distributed in the environment and mostly disseminated through airborne spores, which can be considered infective propagules responsible mainly for respiratory (lung and sinuses), wound, and skin infections (2).

Patients at risk for mucormycosis are immunocompromised (hematological malignancies, hematopoietic stem cell transplantation, solid organ transplant, and steroid therapy) or have diabetes mellitus, deferoxamine treatment, trauma, or severe burns (2-4). Among skin-related infections, contaminated materials (Elastoplast bandages, tape, tongue depressors, ostomy bags, and linens) have been implicated as the causes of local or disseminated infections in patients with various underlying diseases (5). Specifically, in burn patients, invasive wound mucormycosis (IWM) has been reported in both small series (6) and epidemiological surveys (7-10).

Over the past 10 years, outbreaks of mucormycosis have been increasingly reported in various environments. In humans, nosocomial outbreak cases occurring before 2008 were reviewed by Antoniadou in 2009 (11). Antoniadou found 12 reported outbreaks and two pseudoepidemics of cases since 1977. Mucormycosis outbreaks have been reported in the United States, United Kingdom, and Europe (11). Since 2008, outbreaks or clustered cases have been reported after the tornadoes in Joplin, MO (13 patients) (12, 13), in an intensive care unit (ICU) in France (3 patients) (14), in adults (6 patients) (15) or infants (5 patients) (16) exposed to contaminated linens in the United States, in infants in Egypt (5 patients) (17), and in patients undergoing arthroscopy in Argentina (40 patients) (18). More specifically, in a Belgian burn unit, Christiaens et al. described an outbreak of *Lichtheimia corymbifera* associated with nonsterile Elastoplast bandage contamination in seven burn patients, including five with infection and two with colonization (19). In this study, the authors did not have evidence for the genotypic relatedness of the strains between patients and material strains. In addition, a large outbreak due to yogurt contamination in the United States responsible for digestive symptoms (nausea, cramps, vomiting, and diarrhea) in about 300 individuals has been described recently (20). In this study, phylogenetic analysis and whole-genome sequence (WGS) analysis yielded new information on the genetic structure of *Mucor circinelloides*. *Mucor circinelloides* is a single species consisting of four different *formae* (*f. circinelloides*, *f. griseocyanus*, *f. janssenii*, and *f. lusitanicus*), with *forma circinelloides* the most commonly involved in human mucormycosis (20, 21).

Between 2013 and 2015, we faced an outbreak of proven IWM due to *M. circinelloides* f. *circinelloides* in a burn unit (BU) in Hôpital Saint-Louis (SLS), Paris, France, involving six patients raising the hypothesis of a common source of contamination. The outbreak was not suspected until *M. circinelloides* was recovered from wounds of
<figure><img src="image_2.png" /><figcaption>TABLE 1 Isolates sequenced in this study</figcaption></figure>

<table><thead><tr><th>Hospital and ward<sup>a</sup></th><th>Isolate ID</th><th>Date of recovery (day/mo/yr)</th><th>Day from index isolate</th><th>Sample type</th><th>Patient ID</th><th>Site</th><th>EORTC/MSG classification</th><th>Outcome</th><th>Cluster no.</th><th>Strain no.<sup>b</sup></th></tr></thead><tbody><tr><td rowspan="15">SLS BU1</td><td>P01_617_BU1_SLS</td><td>21/03/13</td><td>0</td><td>Wound</td><td>P01</td><td>Wound</td><td>No</td><td>Alive</td><td>C3</td><td>S12</td></tr><tr><td>P02_783_BU1_SLS</td><td>18/10/13</td><td>211</td><td>Bone</td><td>P02</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S5</td></tr><tr><td>P03_592_BU1_SLS</td><td>20/08/14</td><td>517</td><td>Wound (shoulder)</td><td>P03</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S1</td></tr><tr><td>P03_594_BU1_SLS</td><td>2/9/14</td><td>530</td><td>Wound (hand)</td><td>P03</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S1</td></tr><tr><td>P04_559_BU1_SLS</td><td>29/08/14</td><td>526</td><td>Wound (hand)</td><td>P04</td><td>Wound</td><td>Proven</td><td>Death</td><td>C2</td><td>S9</td></tr><tr><td>P04_601_BU1_SLS</td><td>2/9/14</td><td>530</td><td>Wound (shoulder)</td><td>P04</td><td>Wound</td><td>Proven</td><td>Death</td><td>C2</td><td>S9</td></tr><tr><td>P04_602_BU1_SLS</td><td>4/9/14</td><td>532</td><td>Wound (forearm)</td><td>P04</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S7</td></tr><tr><td>P04_603_BU1_SLS</td><td>8/9/14</td><td>536</td><td>Wound (leg)</td><td>P04</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S7</td></tr><tr><td>P05_599_BU1_SLS</td><td>5/9/14</td><td>533</td><td>Wound (arm 1)</td><td>P05</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S1</td></tr><tr><td>P05_598_BU1_SLS</td><td>5/9/14</td><td>533</td><td>Wound (thigh)</td><td>P05</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S8</td></tr><tr><td>P05_600_BU1_SLS</td><td>5/9/14</td><td>533</td><td>Wound (arm 2)</td><td>P05</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S8</td></tr><tr><td>P05_622_BU1_SLS</td><td>16/09/14</td><td>544</td><td>Wound (forearm)</td><td>P05</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S3</td></tr><tr><td>P06_023_BU1_SLS</td><td>8/1/15</td><td>658</td><td>Wound (leg)</td><td>P06</td><td>Wound</td><td>Proven</td><td>Alive</td><td>C1</td><td>S6</td></tr><tr><td>P06_032_BU1_SLS</td><td>12/01/15</td><td>662</td><td>Wound (knee)</td><td>P06</td><td>Wound</td><td>Proven</td><td>Alive</td><td>C1</td><td>S6</td></tr><tr><td>Pneumology</td><td>P07_621_SLS</td><td>29/04/13</td><td>39</td><td>Sputum</td><td>P07</td><td>Lung</td><td>No</td><td></td><td>C4</td><td>S14</td></tr><tr><td>Hematology<sup>c</sup></td><td>E01_615_SLS</td><td>12/03/14</td><td>356</td><td>Surface swab</td><td></td><td></td><td></td><td></td><td>C2</td><td>S10</td></tr><tr><td rowspan="4">PER BU2</td><td>P08_701_BU2_PER</td><td>3/2/13</td><td>-46</td><td>Wound</td><td>P08</td><td>Wound</td><td>Proven</td><td>Death</td><td>C3</td><td>S13</td></tr><tr><td>P09_704_BU2_PER</td><td>17/02/13</td><td>-32</td><td>Wound</td><td>P09</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S2</td></tr><tr><td>P10_703_BU2_PER</td><td>27/02/13</td><td>-22</td><td>Wound</td><td>P010</td><td>Wound</td><td>Proven</td><td>Death</td><td>C1</td><td>S2</td></tr><tr><td>P11_702_BU2_PER</td><td>25/07/14</td><td>491</td><td>Wound</td><td>P011</td><td>Wound</td><td>Proven</td><td>Death</td><td>C2</td><td>S11</td></tr><tr><td>STR Kidney transplant</td><td>P12_579_STR</td><td>23/07/14</td><td>489</td><td>Kidney biopsy</td><td>P012</td><td>Kidney</td><td>Proven</td><td>Death</td><td>C1</td><td>S4</td></tr></tbody></table>

patient P04. Over the same period of time, 4 additional cases that occurred in a burn unit (Burn Unit 2 [BU2]) of another hospital in a Paris suburb (Hôpital d'Instruction des Armées, Clamart, France [PER]) were reported to the National Reference Center for Invasive Mycoses and Antifungals (NRCMA).

Our aim was to clarify the origin/source of infection in both BU1 and BU2. In the absence of genotyping markers for this organism, WGS analysis was performed on 21 isolates (14 from the outbreaks and 7 unrelated) to investigate the links between clinical isolates, understand the epidemiology of the outbreak, and identify and eliminate the potential source of the infections.

## RESULTS

**Clinical and microbiological investigations.** Three patients (P03, P04, and P05) developed proven IWM due to *M. circinelloides* f. *circinelloides* within 18 days after admission in BU1 (between 18 August and 5 September 2014) and subsequently died from these infections (Table 1). The outbreak was suspected when a positive culture was observed in P04 (11 days after the first positive sample in P03). Sequential samples from the wounds were prospectively obtained starting with P03. A few months later (114 days), another patient (P06) developed proven IWM, and *M. circinelloides* f. *circinelloides* was also involved (Fig. 1). We retrospectively noticed isolates from the same species had already been identified in 2013 from two patients in BU1 (P01 with no infection and P02 with proven IWM). A total of seven patients were exposed to (P01 and P07) and/or infected with (P02 to P06) this species in BU1, with P01 as the putative index case. Infection control measures were implemented locally to avoid potential nosocomial transmission to other patients of the unit. More than 30 environmental samples were cultured. All were negative. DNA amplified with the *Mucor/Rhizopus* PCR test (22) was detected only in the Bair Hugger filters that were used during the hospitalization of P03, P04, and P05.

<footnote><sup>a</sup>SLS, Hôpital Saint-Louis, Paris, France; PER, Hôpital d'Instruction des Armées, Clamart (a Paris suburb); STR, Strasbourg, France; BU1 and -2, Burn Units 1 and 2.</footnote>

<footnote><sup>b</sup>Based on strain relatedness after WGS analysis.</footnote>

<footnote><sup>c</sup>Environmental isolate.</footnote>
<figure><img src="image_3.png" /><figcaption>**FIG 1** Epidemiological map of 13 patients whose isolates were selected in our study. Isolates from several geographic areas in France (from dark to light gray) have been studied: Burn Unit 1 (BU1) of hospital 1, Paris France; wards of SLS, Paris, France; BU2 of PER in the Paris area and eastern France (one isolate from Strasbourg [STR]). Isolates were prospectively collected for patients P03 to P06. Analyzed isolates are depicted in dark open circles. The index case from BU1 (solid arrow) was thought to be P01, and the outbreak was recognized after P04 got infected (dashed arrow).</figcaption></figure>

To prevent transmission to other patients, we needed to investigate whether these strains were clonal and needed unrelated isolates from other geographic areas. The additional cases corresponded to an outbreak in BU2 (PER) involving four patients with IWM, as well as one case of proven invasive mucormycosis with kidney invasion in a transplant recipient in Strasbourg (STR) about whom the NRCMA was notified.

Overall, the 12 patients (21 clinical isolates) included 10 cases identified in two burn units (Table 1; Fig. 1).

**Phylogenetic analyses of three loci.** Internal transcribed spacer (ITS), D1/D2, and RPB1 sequences were analyzed as separate (data not shown) and combined data sets. The topologies of the multilocus data set by 3 methods (neighbor-joining [NJ], maximum likelihood [ML], and Bayesian inference) were comparable between the individual trees of the three genes analyzed. Four clades (Fig. 2), respectively, denoted C1 (14 isolates, including 11 from BU1), C2 (4 isolates, including 2 from BU1), C3 (2 strains in addition to the reference strain, including 1 from BU1), and C4 (1 isolate), were identified from the analysis of the combined data set, which yielded a significant support (≥95% bootstrap for NJ and ML; 1.0 for Bayesian inference). Isolates recovered from BU1 were distributed in three clades (C2, C3, and C4). All patients from BU1 had at least one isolate included in C1.

**Whole-genome analysis.** To better resolve the diversity of the strains within the four clades and because no further genotyping methods existed for this organism, whole-genome sequencing was performed. Because the biology and the genetics of this organism are poorly understood, we first checked the reproducibility of the sequencing process and the stability of the genome, to be able to define genetically identical strains.

**Establishing the genetic threshold to determine genetically identical strains.** For the three strains isolated from single-spore colonies (i.e., P05_600_BU1_SLS, P04_603_BU1_SLS, and P03_594_BU1_SLS), two complementary approaches were used to estimate the genetic distance between the parent and single-spore colony: estimates of both the evolutionary distance (defined as the proportion of nucleotide differences:
<figure><img src="image_4.png" /><figcaption>**FIG 2** PhyML tree constructed by the maximum likelihood method. The tree was inferred from concatenated 3-locus data set (ITS, 28S, and RPB1). Bootstrap support values from PhyML greater than 70% (left) and a Bayesian posterior probability of >0.80 (right) are shown at the nodes. The 21 clinical isolates are grouped in clades C1 to C4. The scale bar indicates 0.001 nucleotide substitution per character.</figcaption></figure>

i.e., 0.00035, 0.00044, and 0.00029, respectively) and the number of SNP differences (i.e., 4,731, 5,296, and 4,187, respectively). This information gave us the expected genetic proximity measures between pairs of genomes arising from identical strains and independently sequenced isolates. The largest of the three pairs of measures (i.e., 0.00044 mismatch per nucleotide and 5,296 SNPs) was therefore selected as a cutoff below which two compared isolates were defined as belonging to the same strain.

**Experimental investigation of the potential genetic drift of *M. circinelloides* f. *circinelloides*.** Experimental investigation of the 1006PhL genome upon iterative sub-culturing on agar (n = 3) and three passages in mice (n = 3) revealed no acquisition of SNPs during this process, suggesting that the genome of *M. circinelloides* f. *circinelloides* was stable upon iterative passages.

**Whole-genome phylogenetic classification of the 21 clinical isolates.** A phylogenetic classification of the whole genome of the 21 isolates was then performed
<figure><img src="image_5.png" /><figcaption>**FIG 3** Minimum-evolution phylogenetic tree of the whole genome of 21 clinical isolates and the 1006PhL reference strain. The 21 isolates are grouped in different strains clustered in clades C1 to C4. Clinical isolates belonging to the same strain are highlighted in light gray. The scale bar indicates 0.004 nucleotide substitution per character.</figcaption></figure>

(Fig. 3). This phylogenetic tree allows classification of the genomes in four main clades corresponding exactly to the same clades (C1 to C4) described for the analysis of three loci (Fig. 2). Clades C2, C3, and C4 contained isolates that are clearly distinct from those inside clade C1: e.g., the average estimated evolutionary distances and SNP numbers between isolates from C1 and those inside C2, C3, and C4 are 0.0187, 0.0385, and 0.0390 and 291,467, 702,061, and 715,595, respectively, whereas among C1 isolates, the average pairwise distance and SNP number are 0.00146 and 19,588, respectively (see Table S1 in the supplemental material).

Furthermore, increased resolution of WGS allowed robust identification of strains (as defined above and in Materials and Methods) and understanding of which clinical isolates belong to which strain. As a result, the 21 investigated isolates and the reference 1006PhL could be partitioned into 14 distinct strains (S1 to S14 [Fig. 3; Table 1]).
**Outbreak dynamics.** In BU1, the isolate of the potential index case P01 (P01_617_BU1_SLS, S12) was different from the isolates subsequently recovered in BU1. The isolate from P02, corresponding to the specific strain S5 clustered in C1 with isolates from P03, P04, P05, and P06. P03 and P06 were also infected with two isolates from S1 (594 and 592) and S6 (032 and 023), respectively. P04 was infected over 10 days with two strains, S9 and S7, belonging to C2 (isolates 601 and 559) and C1 (isolates 602 and 603), respectively. P05 was infected over 11 days with isolates belonging to three strains from C1, S8 (isolates 598 and 600), S3 (isolate 622), and S1 (isolate 599). Indeed, P04 and P05 had mixed infections during the course of their disease, suggesting initial contamination with a mixture of strains, the latter strain also recovered in P03. This suggests cross contamination or common infection in P3 and P5. This has also been observed in BU2, where two patients (P09 and P10) shared strain S2 (isolates 703 and 704), which clustered in C1. Two patients (P11 and P08) were also infected with strains that belonged to C2 and C3, respectively. The environmental (isolate 615) and the colonization (isolate 621) isolates from our hospital clustered in C2 and C4, respectively. The patient from eastern France clustered in C1 (S4), but with a specific strain different from the other C1 strains, as expected for a geographically unrelated infection.

## DISCUSSION

Because isolation of Mucorales is rare in the hospital, the observation of the same species in two independent samples or patients has long been considered a sufficient criterion to suspect and assess transmission or common contamination. Here, we investigated further outbreak-related and -unrelated isolates based on WGS analysis. To our knowledge, this is the first time WGS analysis has been used to resolve an outbreak of invasive Mucorales infection in the context of nosocomial acquisition.

Because WGS was applied for the first time in this setting, we first sought to evaluate reproducibility of the sequencing process and intraculture variation/stability during *in vivo* passage. As the genome of three selected strains was sequenced and assembled twice, we were able to compare contig sets belonging *a priori* to the same strain. This method led to the definition of a pairwise distance cutoff. Therefore, if two genome sequences belonging to different isolates have a pairwise distance below this cutoff, it was inferred *a posteriori* that the two isolates corresponded to the same strain. This cutoff should vary as a function of the organism, the method of sequencing, the bioinformatics pipeline, and the pathophysiology of the disease, suggesting that such data should be obtained each time an investigation of an outbreak due to rare organisms is undertaken. Of note, using pairwise evolutionary distances could be considered as a fast but accurate alternative to the well-known average nucleotide identity (ANI) approach to compare genomes (23, 24) because both were shown to be linearly correlated (25).

Contrary to the initial hypothesis of a single-strain transmission in BU1, we observed that all of the patients from the BU1 outbreak (P01 to P06) were infected by different strains. Surprisingly, our data revealed that each strain was patient specific in BU1, except for S1, suggesting that the outbreak in BU1 was due to multiple strains present in and acquired from a local environmental "reservoir" containing clonally unrelated isolates. This hypothesis is reinforced by two patients (P04 and P05) with IWM coinfected by more than 2 genetically distinct strains. Another hypothesis is that the patients could have been exposed to specific strains or a mixture of strains before arriving in BU1. However, a delay between admission and the first positive culture was 16 days (median), making the hypothesis that exposure occurred in the environment of BU1 more likely. In the settings of severe burns, where invasive fungi do sporulate on the wounds, transmission by air and transmission by the hands of health care workers to other patients are both possible. A major point to emphasize is that only two patients from BU1 (P03 and P05) shared the same strain, S1. This feature has also been observed in BU2, suggesting that transmission between patients is possible. However, this does not rule out the hypothesis of a contamination by the same strain from the
environment. Thus, while we cannot entirely exclude that mating and also recombination could occur on the surface of the wound, zygospores of *Mucor circinelloides* have never been observed in wound tissue and are not known to germinate under any conditions tested thus far *in vitro*. Moreover, the conditions of the host environment (elevated temperature and light) are not optimal for sexual reproduction, and thus it seems unlikely that mating and recombination would be occurring in the setting of infection.

Environmental investigation of the outbreak in BU1 failed to identify the source of infection using culture of multiple samples from the environment, as well as by PCR. DNA amplified with the *Mucor/Rhizopus* PCR (22) was detected only in the Bair Hugger filters that were used during the hospitalization of IWM patients P03, P04, and P05. Despite the negative result of this investigation, it is likely that the contamination came from a local source because this has already been described with linens or Elastoplast in burn units (5, 15, 16).

Our findings about the genetic structure of *M. circinelloides* f. *circinelloides* are reminiscent of the WGS investigations of the *Apophysomyces* sp. outbreak in Joplin (13) or of the *Saprochaete clavata* outbreak in France (26), which revealed that several genetic groups can be responsible for infections over the same period of time. In our case, in a given restricted area (BU1), we identified a large diversity of isolates responsible for IWM and were not able to find isolates belonging to unique strains recovered in different places. At the other end of this spectrum, for *Apophysomyces trapeziformis* several genetically identical isolates were recovered in different places at a distance of several miles (13). In the case of *Exserohilum rostratum*, all outbreak isolates have closely related genomes, suggesting that a unique strain was responsible for the outbreak (27). These differences could be explained by a completely different pathophysiology of the disease and mode of fungal transmission between these fungal organisms.

This outbreak is illustrative of the importance of the sampling strategy. Repeated sampling is paramount, as is avoiding the assumption that one patient should harbor only one strain. Indeed, in our case two patients were infected by mixtures of strains concomitantly, as already described for cryptococcosis (28). Mixed infections and the impact of genetic heterogeneity of isolates from single patients on the interpretation of transmission chains are an emerging theme in molecular epidemiology in the genomic era (29). Our investigation fully supports the view that multisampling is critical to decipher transmission patterns, especially in the context of outbreaks with such long time frames. Here, multiple sampling was only performed prospectively for four patients (P03, P04, P05, and P06), with one isolate stored and studied for all of the remaining patients. The proportion of mixed infection may thus have been higher than detected here if all of the isolates from all of the patients had been investigated. Based on our experience, to investigate a cluster of wound mucormycosis, we recommend sampling and analyzing several wound lesions and several isolates per patient and analyzing isolates recovered at different times. This will have the potential to find the maximal diversity of the strains if present and have a better view of the outbreak dynamics. In outbreaks involving rare organisms for which the biology and the genetics are mostly unknown, we recommend testing the stability of the genome under different conditions, including *in vivo* models to determine the threshold above which to strains can be delineated.

To prevent potential transmission between patients, all patients with a positive culture were isolated and had a dedicated nurse. To avoid persistence of spores in the environment, after a patient's death or discharge, rooms should be cleaned by using vaporization of 6% H₂O₂. In addition, our management of this outbreak led to the implementation of a twice-weekly screening of serum samples from patients hospitalized in BU1 using Mucorales PCR and by immediate prescription of antifungal treatment when the PCR test was positive (22), in addition to isolation of all culture-positive and/or PCR-positive patients and dedicated nurses to prevent the risk of transmission to other patients. So far, with 2 years of hindsight, no additional cases of *M. circinelloides* f. *circinelloides* IWM have been observed. This illustrates how prevention through hygiene
together with early diagnosis of mucormycosis could improve patient management and
help control dramatic outbreaks in hospital settings among populations at risk.

## MATERIALS AND METHODS

**Isolates and patients.** In order to study the genetic relatedness between the BU1 clinical isolates, the following additional isolates identified as *M. circinelloides* from other sources were selected (Table 1). (i) The first group of isolates were sequential isolates prospectively collected from the skin lesions of 4 of the patients, as previously described (22). (ii) The next isolates were collected in Hôpital Saint-Louis (SLS) but in another ward (one from the environment and another one colonizing a patient). (iii) The next isolates (n = 4) were recovered in another outbreak that involved 4 patients hospitalized in another burn unit (BU2) in the Hopital d'Instruction des Armées, Clamart (PER), located in the suburb of Paris, over the same period of time. These isolates have been sent to the French National Reference Center for Invasive Mycoses and Antifungals (NRCMA). (iv) Finally, one isolate was obtained from invasive mucormycosis recovered from the kidney biopsy specimen of a kidney transplant recipient in Strasbourg (STR) in the eastern part of France.

All positive slants were subcultured once on Sabouraud dextrose agar with gentamicin and chlor-
amphenicol (Bio-Rad, Marnes-la-Coquette, France) at 30°C using the bulk and never single colonies.

Overall, 21 isolates were selected for further analysis: 15 clinical isolates from BU1 (SLS) (12 bulk cultures, including 1 from colonization), 1 environmental isolate from SLS, 4 clinical isolates from BU2 (PER), and 1 clinical isolate from STR.

The sequence of the *Mucor circinelloides* 1006PhL strain was used as the reference genome (30).

Five of the six patients had proven IWM according to a modified version of the EORTC/MSG criteria (22, 31), and one was considered colonized and was not treated.

**Environmental investigation in BU1.** Extensive environmental sampling was performed, and mycological contamination was investigated by culture methods on Sabouraud agar (Bio-Rad, Marnes-la-Coquette, France) and 2% malt extract incubated at 30 and 37°C. Overall, 30 specimens from nonsterile material, air, surfaces, and aeration machineries, the technical room, Bair Hugger machines, and dedicated nonsterile material were tested. The Bair Hugger machine allows active warming of the patient. It pulses warm air through a plastic pipe into a blanket applied on the patient.

**Polyphasic identification of isolates.** The 21 isolates were sent to the NRCMA, where the purity was verified and identification to the species-level performed using phenotypic and molecular identification on the bulk culture. In detail, microscopic examination was performed on 5- to 7-day-old cultures growth on 2% malt agar at 30°C. Amplification and sequencing of the ITS1-5.8S-ITS2 region and the D1/D2 region of the large subunit (LSU) ribosomal DNA (rDNA) were performed as described previously (32). The amplification of the RPB1 gene (RNA polymerase II largest subunit) was made with primers RPB1Ac and RPB1Cr (33). The PCR products were then sequenced, and the consensus sequences were obtained as previously described (32). Sequences were subjected to pairwise alignments against curated fungal reference databases available at the online MycoBank database (http://www.mycobank.org/).

**Sequence alignment and phylogenetic analysis.** Multilocus sequence alignments of partial sequences of 3 different loci (ITS, 28S, and RPB1) were performed (32, 34).

The sequences were aligned using MAFFT v.7.308 with default settings. Data from each gene were analyzed separately and combined as a concatenated 3-locus data set. For the multilocus data set, neighbor-joining (NJ) phylogenetic trees were constructed using MEGA6 software (35) with the Tamura 3-parameter substitution model and 1,000 bootstrap replicates. The program PhyML v3.0.1 (36) was used to infer maximum likelihood (ML) phylogeny using the TN93 substitution model and 1,000 bootstrap repetitions. Bayesian analysis with default prior of MrBayes v.3.2 (37) was conducted to determine posterior probabilities. Two analyses were done by running 10⁶ generations in four chains, sampling every 100 generations.

**Whole-genome sequencing and assembly.** Because genotyping methods for *M. circinelloides* f. *circinelloides* were lacking, WGS was performed to compare isolates. Libraries were constructed using a Nextera XT DNA sample preparation kit (Illumina) and sequenced with an Illumina NextSeq 500 sequencing system with a 2× 150-nucleotide paired-end protocol. Two lanes from these tagged libraries resulted in ~12.6 million read pairs per strain on average. Statistics of assembly for all the isolates are summarized in Table 2.

All sequenced reads were clipped and trimmed with AlienTrimmer v.0.4.0 (38), corrected with Musket v.1.1 (39), merged (if needed) with FLASH v.1.2.11 (40), and subjected to a digital normalization procedure with Khmer v.2.0 (41). For each sample, processed reads were finally assembled with SPAdes v.3.10 (42).

**Whole-genome analysis.** For each pair of assembled genomes, an evolutionary distance (proportion of aligned nucleotide differences) was estimated with Mash v.1.0.2 (sketch size = 100,000) (25). For each sample, sequenced reads were also aligned with minimap2 v.2.9 (https://github.com/lh3/minimap2) against the assembled contigs of the reference strain 1006PhL in order to estimate the number of SNPs between each pair of isolates. Both pairwise evolutionary distance and SNP number difference estimates are available in Table S1 in the supplemental material. Finally, the matrix of evolutionary distances was used to infer a minimum-evolution phylogenetic tree with FastME v.2.1.5 (43, 44).

**Reproducibility of the sequencing process.** In a specific experiment dedicated to determine the reproducibility of all parts of the sequencing process (from extraction to sequence analysis), single-spore isolation was performed for 3 isolates recovered from patients P03, P04, and P05 in BU1 (P05_600hti_SLS, P04_603hti_SLS, and P03_594hti_SLS [Table 1]). Spore suspensions were prepared in sterile water with 0.05% Tween 80 and serially diluted, and spores were enumerated to
<figure><img src="image_6.png" /><figcaption>TABLE 2 Statistics of assembly of the strains sequenced in this study</figcaption></figure>

<table><thead><tr><th>Isolate ID</th><th>No. of contigs</th><th>No. of bases</th><th>N<sub>50</sub></th><th>Coverage depth</th></tr></thead><tbody><tr><td>P01_6171_SLS</td><td>8,462</td><td>34,880,780</td><td>13,803</td><td>67x</td></tr><tr><td>P02_783_111_SLS</td><td>4,642</td><td>34,558,450</td><td>25,365</td><td>150x</td></tr><tr><td>P03_592_111_SLS</td><td>5,295</td><td>34,215,243</td><td>22,844</td><td>80x</td></tr><tr><td>P03_594_111_SLS</td><td>4,637</td><td>34,551,072</td><td>26,027</td><td>162x</td></tr><tr><td>P04_559_111_SLS</td><td>5,585</td><td>34,509,945</td><td>21,800</td><td>103x</td></tr><tr><td>P04_601_111_SLS</td><td>4,795</td><td>34,726,648</td><td>24,551</td><td>150x</td></tr><tr><td>P04_602_111_SLS</td><td>5,041</td><td>34,155,516</td><td>24,450</td><td>85x</td></tr><tr><td>P04_603_111_SLS</td><td>5,479</td><td>34,215,113</td><td>23,173</td><td>76x</td></tr><tr><td>P05_598_111_SLS</td><td>5,257</td><td>34,139,893</td><td>23,651</td><td>70x</td></tr><tr><td>P05_599_111_SLS</td><td>5,247</td><td>34,181,601</td><td>23,820</td><td>94x</td></tr><tr><td>P05_600_111_SLS</td><td>5,562</td><td>34,319,442</td><td>21,377</td><td>117x</td></tr><tr><td>P05_622_111_SLS</td><td>5,760</td><td>35,087,509</td><td>19,734</td><td>55x</td></tr><tr><td>P06_023_111_SLS</td><td>4,648</td><td>34,540,329</td><td>25,831</td><td>174x</td></tr><tr><td>P06_032_111_SLS</td><td>4,950</td><td>34,519,384</td><td>22,633</td><td>130x</td></tr><tr><td>P07_621_SLS</td><td>7,249</td><td>34,741,059</td><td>15,678</td><td>63x</td></tr><tr><td>E01_615_SLS</td><td>5,446</td><td>34,322,336</td><td>18,029</td><td>68x</td></tr><tr><td>P08_701_112_PER</td><td>6,530</td><td>35,386,009</td><td>19,860</td><td>202x</td></tr><tr><td>P09_704_112_PER</td><td>5,090</td><td>34,266,372</td><td>19,215</td><td>124x</td></tr><tr><td>P10_703_112_PER</td><td>5,081</td><td>34,495,044</td><td>20,978</td><td>162x</td></tr><tr><td>P11_702_112_PER</td><td>4,100</td><td>33,816,248</td><td>27,537</td><td>133x</td></tr><tr><td>P12_579_STR</td><td>4,757</td><td>34,508,712</td><td>24,438</td><td>168x</td></tr><tr><td>1006PHL</td><td>4,316</td><td>34,938,514</td><td>43,860</td><td>148x</td></tr><tr><td>P03_594_111_SLS_CLONE</td><td>5,649</td><td>34,391,090</td><td>17,573</td><td>61x</td></tr><tr><td>P04_603_111_SLS_CLONE</td><td>5,849</td><td>34,227,008</td><td>19,061</td><td>74x</td></tr><tr><td>P05_600_111_SLS_CLONE</td><td>5,684</td><td>34,478,201</td><td>18,116</td><td>73x</td></tr><tr><td>SEC_INF1</td><td>4,162</td><td>34,642,463</td><td>42,044</td><td>77x</td></tr><tr><td>SEC_INF3</td><td>3,914</td><td>34,582,644</td><td>43,611</td><td>99x</td></tr><tr><td>SEC_VEG1</td><td>4,293</td><td>34,619,149</td><td>39,025</td><td>52x</td></tr><tr><td>SEC_VEG2</td><td>3,986</td><td>34,622,087</td><td>44,282</td><td>109x</td></tr><tr><td>SEC_VEG3</td><td>4,104</td><td>34,662,622</td><td>43,314</td><td>73x</td></tr></tbody></table>

approximately plate five spores per potato dextrose agar (PDA) plate. After 12 h of incubation at 30°C, plates were inspected under a microscope at ×10 magnification, and germinating spores were then transferred to new PDA plates using a stereomicroscope. One colony was thus selected from each isolate for additional sequencing.

**Genome stability experiments.** Genome stability during vegetative growth and host infections was analyzed with the 1006PhL strain. In brief, the strain was grown on PDA for 1 day. The colony was streaked to isolate a single colony, which was then transferred onto a new PDA plate. After 4 days of incubation at 30°C under the light, the spores were collected and subjected to the same procedures up to three passages (VEG1, VEG2, and VEG3 isolates).

For infection passages, spores were suspended in sterile phosphate-buffered saline (PBS). Male 8-week-old BALB/c mice were infected with 10⁶ spores in 200 µl of sterile PBS via tail vein injection. At day 3 postinoculation, the mice were sacrificed, and their brains were collected and placed onto PDA after brain homogenization. After 1 day of incubation, the fungal colony emerging from the brain was streaked to isolate single colonies. One colony was transferred onto PDA and incubated at 30°C under light for 4 days. Spores were then used as the inoculum for the next infection passage, up to three passages (INF1, INF2, and INF3 isolates). Genomic DNAs of these isolates were used to construct libraries with 180-base fragments and 2- to 3-kb jumps, which were sequenced by an Illumina HiSeq 2000 platform at the University of North Carolina High-Throughput Sequencing Facility (HTSF). Data were analyzed using the same pipelines as that implemented for the clinical isolates.

**Ethics statement.** The murine infection experiment was conducted at the Duke University Medical Center in full compliance with all of the guidelines of the Duke University Medical Center Institutional Animal Care and Use Committee (IACUC) and in full compliance with the United States Animal Welfare Act (Public Law 98-198). The Duke University Medical Center IACUC approved all of the vertebrate animal studies under protocol no. A061-12-03. The studies were conducted in the Division of Laboratory Animal Resources (DLAR) facilities, which are accredited by the Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC).

**Data availability.** All FASTQ files were made available on GenBank under accession no. ERS2032288 to ERS2032317 in BioProject PRJEB23686.

## SUPPLEMENTAL MATERIAL

Supplemental material for this article may be found at https://doi.org/10.1128/mBio .00573-18.

TABLE S1, XLSX file, 0.1 MB.
<figure><img src="image_7.png" /></figure>

ACKNOWLEDGMENTS

We thank Vincent Enouf from the P2M (Plateforme de Microbiology Mutualisée), Institut Pasteur, Paris, for supervising and organizing the sequencing of the *M. circinelloides* f. *circinelloides* isolates. We also thank Damien Hoinard (French National Reference Center for Invasive Mycoses and Antifungals) for technical assistance.

This study has been supported by Institut Pasteur and Santé Publique France. J.H. is supported by NIH/NIAID R37 MERIT award AI39115-20, NIH/NIAID R01 award AI050113-13, and Astellas Award CRES-15G01.

REFERENCES

1. Hibbett DS, Binder M, Bischoff JF, Blackwell M, Cannon PF, Eriksson OE, Huhndorf S, James T, Kirk PM, Lücking R, Thorsten Lumbsch H, Lutzoni F, Matheny PB, Mclaughlin DJ, Powell MJ, Redhead S, Schoch CL, Spatafora JW, Stalpers JA, Vilgalys R, Aime MC, Aptroot A, Bauer R, Begerow D, Benny GL, Castlebury LA, Crous PW, Dai YC, Gams W, Geiser DM, Griffith GW, Gueidan C, Hawksworth DL, Hestmark G, Hosaka K, Humber RA, Hyde KD, Ironside JE, Köljalg U, Kurtzman CP, Larsson KH, Lichtwardt R, Longcore J, Miadlikowska J, Miller A, Moncalvo JM, Mozley-Standridge S, Oberwinkler F, Parmasto E, Reeb V, Rogers JD, Roux C, Ryvarden L, Sampaio JP, Schüssler A, Sugiyama J, Thorn RG, Tibell L, Untereiner WA, Walker C, Wang Z, Weir A, Weiss M, White MM, Winka K, Yao YJ, Zhang N. 2007. A higher-level phylogenetic classification of the Fungi. Mycol Res 111:509–547. https://doi.org/10.1016/j.mycres.2007.03.004.

2. Roden MM, Zaoutis TE, Buchanan WL, Knudsen TA, Sarkisova TA, Schaufele RL, Sein M, Sein T, Chiou CC, Chu JH, Kontoyiannis DP, Walsh TJ. 2005. Epidemiology and outcome of zygomycosis: a review of 929 reported cases. Clin Infect Dis 41:634–653. https://doi.org/10.1086/432579.

3. Petrikkos G, Skiada A, Lortholary O, Roilides E, Walsh TJ, Kontoyiannis DP. 2012. Epidemiology and clinical manifestations of mucormycosis. Clin Infect Dis 54:523–534. https://doi.org/10.1093/cid/cir866.

4. Lanternier F, Dannaoui E, Morizot G, Elie C, Garcia-Hermoso D, Huerre M, Bitar D, Dromer F, Lortholary O, The French Mycosis Study Group. 2012. A global analysis of mucormycosis in France: the RetroZygo Study (2005–2007). Clin Infect Dis 54:S35–S43. https://doi.org/10.1093/cid/cir880.

5. Rammaert B, Lanternier F, Zahar JR, Dannaoui E, Bougnoux ME, Lecuit M, Lortholary O. 2012. Healthcare-associated mucormycosis. Clin Infect Dis 54:544–554. https://doi.org/10.1093/cid/cir867.

6. Rabin ER, Lundberg GD, Mitchell ET. 1961. Mucormycosis in severely burned patients. Report of two cases with extensive destruction of the face and nasal cavity. N Engl J Med 264:1286–1289. https://doi.org/10.1056/NEJM196106222642504.

7. Church D, Elsayed S, Reid O, Winston B, Lindsay R. 2006. Burn wound infections. Clin Microbiol Rev 19:403–434. https://doi.org/10.1128/CMR.19.2.403-434.2006.

8. Schaal JV, Leclerc T, Pasquier P, Bargues L. 2012. Epidemiology of fungal infection in burns: therapeutic implications. Burns 38:942–944. https://doi.org/10.1016/j.burns.2012.02.028.

9. Skiada A, Pagano L, Groll A, Zimmerli S, Dupont M, Lagrou K, Lass-Flori C, Bouza E, Klimko N, Gaustad P, Richardson M, Hamal P, Akova M, Meis JF, Rodriguez-Tudela JL, Roilides E, Mitrousia-Ziouva A, Petrikkos G, European Confederation of Medical Mycology Working Group on Zygomycosis. 2011. Zygomycosis in Europe: analysis of 230 cases accrued by the registry of the European Confederation of Medical Mycology (ECMM) Working Group on Zygomycosis between 2005 and 2007. Clin Microbiol Infect 17:1859–1867. https://doi.org/10.1111/j.1469-0691.2010.03456.x.

10. Becker WK, Cioffi WG, McManus AT, Kim SH, McManus WF, Mason AD, Pruitt BA. 1991. Fungal burn wound infection. A 10-year experience. Arch Surg 126:44–48. https://doi.org/10.1001/archsurg.1991.01410250048008.

11. Antoniadou A. 2009. Outbreaks of zygomycosis in hospitals. Clin Microbiol Infect 15:55–59. https://doi.org/10.1111/j.1469-0691.2009.02982.x.

12. Neblett Fanfair R, Benedict K, Bos J, Bennett SD, Lo YC, Adebanjo T, Etienne K, Deak E, Derado G, Shieh WJ, Drew C, Zaki S, Sugerman D, Gade L, Thompson EH, Sutton DA, Engelthaler DM, Schupp JM, Brandt ME, Harris JR, Lockhart SR, Turabelidze G, Park BJ. 2012. Necrotizing cutaneous mucormycosis after a tornado in Joplin, Missouri, in 2011. N Engl J Med 367:2214–2225. https://doi.org/10.1056/NEJMoa1204781.

13. Etienne KA, Gillece J, Hilsabeck R, Schupp JM, Colman R, Lockhart SR, Gade L, Thompson EH, Sutton DA, Neblett-Fanfair R, Park BJ, Turabelidze

G, Keim P, Brandt ME, Deak E, Engelthaler DM. 2012. Whole genome sequence typing to investigate the *Apophysomyces* outbreak following a tornado in Joplin, Missouri, 2011. PLoS One 7:e49989. https://doi.org/10.1371/journal.pone.0049989.

14. Poirier P, Pourrisson C, Gibold L, Chalus E, Guelon D, Descamp S, Traore O, Cambon M, Aumeran C. 2013. Three cases of cutaneous mucormycosis with Lichtheimia spp. (ex Absidia/Mycocladus) in ICU. Possible cross-transmission in an intensive care unit between 2 cases. J Mycol Med 23:265–269. https://doi.org/10.1016/j.mycmed.2013.09.002.

15. Cheng VCC, Chen JHK, Wong SCY, Leung SSM, So SYC, Lung DC, Lee WM, Trendell-Smith NJ, Chan WM, Ng D, To L, Lie AKW, Yuen KY. 2016. Hospital outbreak of pulmonary and cutaneous zygomycosis due to contaminated linen items from substandard laundry. Clin Infect Dis 62:714–721. https://doi.org/10.1093/cid/civ1006.

16. Duffy J, Harris J, Gade L, Sehulster L, Newhouse E, O'Connell H, Noble-Wang J, Rao C, Balajee SA, Chiller T. 2014. Mucormycosis outbreak associated with hospital linens. Pediatr Infect Dis J 33:472–476. https://doi.org/10.1097/INF.0000000000000261.

17. El-Mahallawy HA, Khed R, Taha H, Shalaby L, Mostafa A. 2016. Investigation and management of a Rhizomucor outbreak in a pediatric cancer hospital in Egypt. Pediatr Blood Cancer 63:171–173. https://doi.org/10.1002/pbc.25673.

18. Chaves MS, Franco D, Nanni JC, Basaldúa ML, Boleas M, Aphalo G, Feltes Silva P, Elgadban MC, Emery F, Gamarra S, Garcia-Effron G, Bantar C. 2016. Control of an outbreak of postoperative bone mucormycosis: an intervention study of contiguous cohorts. Am J Infect Control 44: 1715–1717. https://doi.org/10.1016/j.ajic.2016.06.017.

19. Christiaens G, Hayette MP, Jacquemin D, Melin P, Mutsers J, De Mol P. 2005. An outbreak of Absidia corymbifera infection associated with bandage contamination in a burns unit. J Hosp Infect 61:88. https://doi.org/10.1016/j.jhin.2004.12.011.

20. Lee SC, Billmyre RB, Li A, Carson S, Sykes SM, Huh EY, Mieczkowski P, Ko DC, Cuomo CA, Heitman J. 2014. Analysis of a food-borne fungal pathogen outbreak: virulence and genome of a Mucor circinelloides isolate from yogurt. mBio 5:e01390-14. https://doi.org/10.1128/mBio.01390-14.

21. Walther G, Pawłowska J, Alastruey-Izquierdo A, Wrzosek M, Rodriguez-Tudela JL, Dolatabadi S, Chakrabarti A, de Hoog GS. 2013. DNA barcoding in Mucorales: an inventory of biodiversity. Persoonia 30:11–47. https://doi.org/10.3767/003158513X665070.

22. Legrand M, Gits-Muselli M, Boutin L, Garcia-Hermoso D, Maurel V, Soussi S, Benyamina M, Ferry A, Chaussard M, Hamane S, Denis B, Touratier S, Guigue N, Fréalle E, Jeanne M, Shaal JV, Soler C, Mimoun M, Chouat M, Lafaurie M, Mebazaa A, Bretagne S, Alanio A. 2016. Detection of circulating Mucorales DNA in critically ill burn patients: preliminary report of a screening strategy for early diagnosis and treatment. Clin Infect Dis 63:1312–1317. https://doi.org/10.1093/cid/ciw563.

23. Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. 2007. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol 57:81–91. https://doi.org/10.1099/ijs.0.64483-0.

24. Wibberg D, Rupp O, Blom J, Jelonek L, Kröber M, Verwaaijen B, Goesmann A, Albaum S, Grosch R, Pühler A, Schlüter A. 2015. Development of a Rhizoctonia solani AG1-IB specific gene model enables comparative genome analyses between phytopathogenic R. solani AG1-IA, AG1-IB, AG3 and AG8 isolates. PLoS One 10:e0144769. https://doi.org/10.1371/journal.pone.0144769.

25. Ondov BD, Treangen TJ, Melsted P, Mallonee AB, Bergman NH, Koren S, Phillippy AM. 2016. Mash: fast genome and metagenome distance esti-
mation using MinHash. Genome Biol 17:132. https://doi.org/10.1186/
s13059-016-0997-x.

26. Vaux S, Criscuolo A, Desnos-Ollivier M, Diancourt L, Tarnaud C, Vanden-
bogaert M, Brisse S, Coignard B, Dromer F, Geotrichum Investigation
Group. 2014. Multicenter outbreak of infections by *Saprochaete clavata*,
an unrecognized opportunistic fungal pathogen. mBio 5:e02309-14.
https://doi.org/10.1128/mBio.02309-14.

27. Litvintseva AP, Hurst S, Gade L, Frace MA, Hilsabeck R, Schupp JM, Gillece JD, Roe C, Smith D, Keim P, Lockhart SR, Changayil S, Weil MR, MacCannell DR, Brandt ME, Engelthaler DM. 2014. Whole-genome analysis of *Exserohilum rostratum* from an outbreak of fungal meningitis and other infections. J Clin Microbiol 52:3216-3222. https://doi.org/10.1128/JCM.00936-14.

28. Desnos-Ollivier M, Patel S, Spaulding AR, Charlier C, Garcia-Hermoso D, Nielsen K, Dromer F. 2010. Mixed infections and in vivo evolution in the human fungal pathogen *Cryptococcus neoformans*. mBio 1:e00091-10. https://doi.org/10.1128/mBio.00091-10.

29. Didelot X, Fraser C, Gardy J, Colijn C. 2017. Genomic infectious disease epidemiology in partially sampled and ongoing outbreaks. Mol Biol Evol 34:997-1007. https://doi.org/10.1093/molbev/msw275.

30. Findley K, Oh J, Yang J, Conlan S, Deming C, Meyer JA, Schoenfeld D, Nomicos E, Park M, NIH Intramural Sequencing Center Comparative Sequencing Program, Kong HH, Segre JA. 2013. Topographic diversity of fungal and bacterial communities in human skin. Nature 498:367-370. https://doi.org/10.1038/nature12171.

31. de Pauw B, Walsh TJ, Donnelly JP, Stevens DA, Edwards JE, Calandra T, Pappas PG, Maertens J, Lortholary O, Kauffman CA, Denning DW, Patterson TF, Maschmeyer G, Bille J, Dismukes WE, Herbrecht R, Hope WW, Kibbler CC, Kullberg BJ, Marr KA, Muñoz P, Odds FC, Perfect JR, Restrepo A, Ruhnke M, Segal BH, Sobel JD, Sorrell TC, Viscoli C, Wingard JR, Zaoutis T, Bennett JE, European Organization for Research and Treatment of Cancer/Invasive Fungal Infections Cooperative Group, National Institute of Allergy and Infectious Diseases Mycoses Study Group (EORTC/MSG) Consensus Group. 2008. Revised definitions of invasive fungal disease from the European Organization for Research and Treatment of Cancer/Invasive Fungal Infections Cooperative Group and the National Institute of Allergy and Infectious Diseases Mycoses Study Group (EORTC/MSG) Consensus Group. Clin Infect Dis 46:1813-1821. https://doi.org/10.1086/588660.

32. Garcia-Hermoso D, Hoinard D, Gantier JC, Grenouillet F, Dromer F, Dannaoui E. 2009. Molecular and phenotypic evaluation of Lichtheimia corymbifera (formerly Absidia corymbifera) complex isolates associated with human mucormycosis: rehabilitation of L. ramosa. J Clin Microbiol 47:3862-3870. https://doi.org/10.1128/JCM.02094-08.

33. Stiller JW, Hall BD. 1997. The origin of red algae: implications for plastid evolution. Proc Natl Acad Sci U S A 94:4520-4525. https://doi.org/10.1073/pnas.94.9.4520.

34. Alvarez E, Sutton DA, Cano J, Fothergill AW, Stchigel A, Rinaldi MG, Guarro J. 2009. Spectrum of zygomycete species identified in clinically significant specimens in the United States. J Clin Microbiol 47: 1650-1656. https://doi.org/10.1128/JCM.00036-09.

35. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. 2013. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol 30: 2725-2729. https://doi.org/10.1093/molbev/mst197.

36. Guindon S, Gascuel O. 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696-704. https://doi.org/10.1080/10635150390235520.

37. Ronquist F. 2004. Bayesian inference of character evolution. Trends Ecol Evol 19:475-481. https://doi.org/10.1016/j.tree.2004.07.002.

38. Criscuolo A, Brisse S. 2013. AlienTrimmer: a tool to quickly and accurately trim off multiple short contaminant sequences from high-throughput sequencing reads. Genomics 102:500-506. https://doi.org/10.1016/j.ygeno.2013.07.011.

39. Liu Y, Schröder J, Schmidt B. 2013. Musket: a multistage k-mer spectrum-based error corrector for Illumina sequence data. Bioinformatics 29: 308-315. https://doi.org/10.1093/bioinformatics/bts690.

40. Magoč T, Salzberg SL. 2011. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27:2957-2963. https://doi.org/10.1093/bioinformatics/btr507.

41. Crusoe MR, Alameldin HF, Awad S, Boucher E, Caldwell A, Cartwright R, Charbonneau A, Constantinides B, Edvenson G, Fay S, Fenton J, Fenzl T, Fish J, Garcia-Gutierrez L, Garland P, Gluck J, González I, Guermond S, Guo J, Gupta A, Herr JR, Howe A, Hyer A, Härpfer A, Irber L, Kidd R, Lin D, Lippi J, Mansour T, McA'Nulty P, McDonald E, Mizzi J, Murray KD, Nahum JR, Nanlohy K, Nederbragt AJ, Ortiz-Zuazaga H, Ory J, Pell J, Pepe-Ranney C, Russ ZN, Schwarz E, Scott C, Seaman J, Sievert S, Simpson J, Skennerton CT, Spencer J, Srinivasan R, Standage D, Stapleton JA, Steinman SR, Stein J, Taylor B, Trimble W, Wiencko HL, Wright M, Wyss B, Zhang Q, Zyme E, Brown CT. 2015. The Khmer software package: enabling efficient nucleotide sequence analysis. F1000Res 4:900. https://doi.org/10.12688/f1000research.6924.1.

42. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455-477. https://doi.org/10.1089/cmb.2012.0021.

43. Desper R, Gascuel O. 2002. Fast and accurate phylogeny reconstruction algorithms based on the minimum-evolution principle. J Comput Biol 9:687-705. https://doi.org/10.1089/106652702761034136.

44. Lefort V, Desper R, Gascuel O. 2015. FastME 2.0: a comprehensive, accurate, and fast distance-based phylogeny inference program. Mol Biol Evol 32:2798-2800. https://doi.org/10.1093/molbev/msv150.