# The genus Madurella: Molecular identification and epidemiology in Sudan

Elhadi A. Ahmed, Bakri Y. M. Nour, Adam D. Abakar, Samirah Hamid, Ahmed A. Mohamadani, Mohamed Daffalla, Mogahid Mahmoud, Hisham N. Altayb, Marie Desnos-Ollivier, Sybren de Hoog, et al.

► To cite this version:

Elhadi A. Ahmed, Bakri Y. M. Nour, Adam D. Abakar, Samirah Hamid, Ahmed A. Mohamadani, et al.. The genus Madurella: Molecular identification and epidemiology in Sudan. PLoS Neglected Tropical Diseases, 2020, 14 (7), pp.e0008420. 10.1371/journal.pntd.0008420. pasteur-03243542

HAL Id: pasteur-03243542

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

Submitted on 31 May 2021

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>
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# The genus *Madurella*: Molecular identification and epidemiology in Sudan

Elhadi A. Ahmed<sup>1*</sup>, Bakri Y. M. Nour<sup>2</sup>, Adam D. Abakar<sup>3</sup>, Samirah Hamid<sup>2</sup>, Ahmed A. Mohamadani<sup>4</sup>, Mohamed Daffalla<sup>5</sup>, Mogahid Mahmoud<sup>5</sup>, Hisham N. Altayb<sup>6</sup>, Marie Desnos-Ollivier<sup>7</sup>, Sybren de Hoog<sup>8,9</sup>, Sarah A. Ahmed<sup>8,9,10</sup>

1 Department of Medical Microbiology, Faculty of Medical Laboratory Sciences, University of Gezira, Sudan,

2 Blue Nile National Institute for Communicable Diseases (BNNICD), University of Gezira, Sudan,

3 Department of Medical Parasitology, Faculty of Medical Laboratory Sciences, University of Gezira, Sudan,

4 Department of Pathology, Faculty of Medicine, University of Gezira, Sudan, 5 Department of Surgery, Faculty of Medicine, University of Gezira, Sudan, 6 Biochemistry Department, Faculty of Sciences, King Abdulaziz University, Saudi Arabia, 7 Institut Pasteur, Molecular Mycology Unit, National Reference Center for Invasive Mycoses and Antifungals, Paris, France, 8 Foundation Atlas of Clinical Fungi, Hilversum, The Netherlands, 9 Center of Expertise in Mycology of Radboud University Medical Center / Canisius Wilhelmina Hospital, Nijmegen, The Netherlands, 10 Faculty of Medical Laboratory Sciences, University of Khartoum, Khartoum, Sudan

* hadilabone@yahoo.com

## Abstract

Eumycetoma (mycotic mycetoma) is the fungal form of mycetoma, a subcutaneous infection occurring in individuals living in endemic areas of the disease. The Sudan is hyperendemic for mycetoma, with the highest incidence being reported from Gezira State, Central Sudan. The present study was conducted at the Gezira Mycetoma Center and aimed to determine the cause of black-grain eumycetoma in the state and describe its epidemiology. Black-grain specimens were collected during the surgical operation and direct detection of the causative agent was performed using *M. mycetomatis* species-specific PCR and ITS PCR followed by sequencing. Black-grain was reported from 93.3% of all confirmed mycetoma cases (n = 111/119), with a prevalence in young males. Of the 91 samples subjected to direct PCR, 90.1% (n = 82) gave positive results. The predominant species (88.2%) was *Madurella mycetomatis*. One sample was identified as *M. fahalii*, one as *M. tropicana*, and one matched the phytopathogenic species *Sphaerulina rhododendricola*. The highest endemic zones were Southern Gezira (76.6%) and Northern Sinnar (23.4%). The study confirmed that direct molecular detection on grains provides rapid and specific diagnosis of agents of eumycetoma.

## Author summary

Eumycetoma is a neglected fungal disease endemic in Africa, India, and Latin America. Black-grain eumycetoma is the most common type in Africa and is mainly caused by *Madurella* spp. The Sudan, and in particular Gezira State, central Sudan is hyperendemic for black-grain eumycetoma. Patients with this type of mycetoma are treated with surgery in combination with antifungal therapy. In this study, we collected surgical biopsies from

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

**OPEN ACCESS**

**Citation:** Ahmed EA, Nour BYM, Abakar AD, Hamid S, Mohamadani AA, Daffalla M, et al. (2020) The genus *Madurella*: Molecular identification and epidemiology in Sudan. PLoS Negl Trop Dis 14(7): e0008420. https://doi.org/10.1371/journal.pntd.0008420

**Editor:** Abdallah M. Samy, Faculty of Science, Ain Shams University (ASU), EGYPT

**Received:** February 12, 2020

**Accepted:** May 23, 2020

**Published:** July 30, 2020

**Copyright:** © 2020 Ahmed et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** The authors received no specific funding for this work.

**Competing interests:** The authors have declared that no competing interests exist.
patients attending Gezira Mycetoma Center to directly identify the etiology of black-grain eumycetoma in this state. We also studied the epidemiology of the disease based on the demography of the patients' population. Our result showed that the highest endemic regions were Southern Gezira (76.6%) and Northern Sinnar (23.4%). By applying direct PCR and sequencing we confirmed that the most common etiology of the disease is *Madurella mycetomatis* (88.2%). In addition, we found one case of *M. fahalii* and the first Sudanese case of *M. tropicana* and *Sphaerulina rhododendricola*.

**Introduction**

Black-grain mycetoma is caused exclusively by fungi, and to date around 23 species have been reported as being responsible for this infection [1, 2, 3]. Agents of black grain mycetoma are affiliated to three orders of Ascomycota, viz. *Pleosporales*, *Sordariales* and *Chaetothyriales* [1]. Although eight species of *Chaetothyriales*, belonging to *Exophiala* and *Cladophialophora* have been incriminated in mycetoma, members of this order are less frequent. The largest diversity of black-grain agents is found in *Pleosporales*: ten species belonging to eight genera have been reported [1, 4, 5, 6]. In contrast, *Sordariales* harbor only a small number of species able to cause mycetoma, but among these are the predominant agents worldwide in the genus *Madurella* [7, 8].

Irrespective of phylogenetic positions, all agents of black-grain mycetoma produce similar clinical features, i.e. progressive tumefaction, sinuses, and black grains in discharge or embedded in tissue [9]. For this reason, it is difficult to differentiate the 23 species histologically [1,10]. After culturing of the grains, members of *Chaetothyriales* produce dark olivaceous black colonies and microscopically yeast cells or more or less characteristic conidia [4]. Species of *Pleosporales* produce dark grey or olive green velvety or woolly colonies, but clinical strains may remain sterile or produce ascomata or conidiomata with delay [6, 11]. Members of *Sordariales* produce brownish to dark yellow or grey colonies which are mostly cerebriform and glabrous or with very short aerial mycelium. *Sordariales* species usually remain sterile, only chlamydospores or large swollen hyphae being observed [4]. Culture differences can be used for approximate differentiation of agents of the three orders, while identification down to the species level requires molecular techniques [5, 10, 12].

The importance of species identification in black-grain eumycetoma infection has been indicated in several studies [10, 11, 12, 13]. The pathogens differ in susceptibility to antifungal agents, as well as in their epidemiology [12, 14]. For example, in *Pleosporales*, *Medicopsis romeroi* showed low susceptibility to several antifungal agents including amphotericin B, fluconazole, itraconazole and caspofungin [14, 15, 16]. In addition, differences are noted between the closely interrelated species of *Madurella* in their susceptibility to the main drug used for treating eumycetoma, itraconazole. While *M. mycetomatis*, the predominant species worldwide, is inhibited by low concentration of itraconazole, *M. fahalii* has high Minimum Inhibitory Concentrations (MICs) to this drug [12, 17]. Of note, the epidemiology of the two species may be overlapping, as both species are found in Sudan [12, 18].

Sudan is a hyperendemic region for mycetoma and it is considered as the disease homeland [2, 19]. The incidence of newly diagnosed cases per year in the Mycetoma Research Centre in Khartoum is 355 [18]. Of the 6,792 patients reported with myceoma in this center, 2,476 patients (37%) were from Gezira State, Central Sudan [18]. Therefore, in 2012, the Gezira Mycetoma Center was established to provide medical care for mycetoma patients in that region. During the period January 2013 to December 2016, a total of 1,654 patients were seen
at the Gezira Mycetoma Center; of them, 584 underwent surgery, including amputation in 71 patients (12.2%). Diagnoses were made by clinical evaluation, imaging, fine-needle aspiration cytology and histology [20].

Although culture and molecular identification of the causative pathogen is considered the gold standard for diagnosis of mycetoma, culturing is rarely performed even in highly specialized centers [10, 5]. One of the reasons is that mycetoma agents require several weeks to grow, and risks of absence of growth or contamination are very high [10, 21]. Recently, culture-free methods have been introduced for direct detection of agents of eumycetoma in clinical materials. These include recombinase polymerase amplification (RPA), loop-mediated isothermal amplification (LAMP), and direct PCR amplification followed by sequencing [22]. Despite the high accuracy of these methods, none of them has been installed for routine diagnostics in endemic regions.

In the present study, we applied direct PCR and sequencing for diagnosis of black-grain eumycetoma cases seen at the Gezira Mycetoma Center. Since *M. mycetomatis* is known as the main agent of mycetoma in Sudan, we also applied species-specific PCR using primers described by Ahmed *et al.* [25]. Furthermore, we studied the epidemiology of the disease in the Gezira State based on stored information on the patient population.

## Materials and methods

### Ethical considerations

The ethical approval for this study was obtained from Faculty of Medical Laboratory Sciences, University of Gezira and from The Ministry of Health, Sudan. Since the majority of the patients cannot read or write, verbal consents were obtained before collecting the data and specimens. The consent was approved by the ethical review board. In case of children the consent was taken from parents. Data was processed and analyzed anonymously.

### Study area and population

The study was conducted at the Gezira Mycetoma Center (GMC, Sudan), which is a specialized center for mycetoma diagnosis and treatment. Patients admitted to the center were mainly from Gezira and neighboring States (Sinnar State), but patients from other parts of Sudan were also seen. Gezira State is located in the central-eastern of Sudan, about 116 km from Khartoum, and it is bordered to the South by Sinnar State. Gezira State comprises eight localities with estimated population of 4.9 million, while Sinnar State comprises seven localities with approximately 1.5 million inhabitants. All patients clinically suspected to have mycetoma and underwent surgical intervention between January 2015 and January 2017, were included in the study. Mycetoma was defined by the presence of grains in a discharge or biopsy. In the absence of visible grains, histological sections were used to confirm mycetoma. Data on demography, and clinical information of the patients were collected.

### Sample collection and processing

Biopsies were collected from the patients and the color of the grains was recorded. Of note, the surgical procedure was performed as a part of treatment and diagnostic purposes. Samples with black grains were further processed for specific identification of the etiologic agent. For 17 patients, samples were inadequate and thus were used for diagnostic purposes only. Ninety-four black-grain specimens were collected from excised mycetoma biopsies during surgery and transferred into sterile physiological saline solution. Grains were washed immediately for
several times to remove tissue remnants, blood and inflammatory exudates. The grains were
stored in the saline solution at -20 °C until used.

Genomic DNA extraction

Approximately, 0.5 to 2 g of black-grain material was transferred to a sterile mortar, ground to fine pieces with pestle, and scrapped into a sterile 2 mL Eppendorf tube (Fig 1). CTAB (cetyl trimethyl ammonium bromide) /chloroform-isoamyl alcohol extraction protocol for fungi was followed [24]. In the procedure; 490 µL freshly prepared CTAB buffer and 10 µL of proteinase K (iNtron, Seongnam, Korea) were added to the black-grain material. Mixture was homogenized by vortexing for few minutes and then incubated for one hour at 60°C. After incubation, tubes were shaken to ensure mycelial disruption. An amount of 500 µL of chloroform: isoamyl-alcohol (CDH) (24: 1) were added to the tubes followed by shaking again for 1-2 minutes

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Sample collection and homogenization (A) Mycetoma lesion, (B) Deep excised biopsy, (C) Black-grains specimen, (D) Manual homogenization of the black-grains.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0008420.g001
to form an emulsion. The tubes were spun at 14,000 rpm in a microfuge for 10 minutes. The upper aqueous layer was collected in a new sterile Eppendorf tube containing 0.55 volume ice-cold iso-propanol (stored at -20°C) and the tubes were spun again. Finally, the DNA pellets were washed with 70% ethanol, air-dried and re-suspended in 100 µL TE buffer (10 mM Tris, 1 mM EDTA, pH 8.0). The DNA was stored frozen at -20°C. DNA concentration (µg/mL) and purity were determined using NanoDrop spectrophotometer (Bibby Scientific, UK) at 260 nm wavelength and (OD) ratio at wavelength of 260/280. DNA quality was also assessed in agarose electrophoresis.

## PCR and sequencing

Extracted genomic DNA was amplified using universal fungal primers ITS4/ITS5 and *Madurella mycetomatis* species-specific primers, 26.1A/28.3A (Eurofins Genomics, Ebersberg, Germany) [23, 25]. The reactions were performed using iNtRON's Maxime PCR PreMix Kit with reaction volume of 20 µL that contains in addition to Maxime PCR PreMix: 3 µL genomic DNA, 0.8 µL reverse primer, 0.8 µL forward primer (10 pmol each) and 15.4 µL distilled water. PCR were performed using GeneAmp PCR system 9700 (Applied Biosystems) with the following conditions: primary denaturation (95 for 2 min), then 35 cycles of alternating denaturation (94°C for 30 s), annealing (58°C for 30 s), and extension (72°C for 30 s) with final extension for 2 min. The resulted amplicons were visualized in 1.5% agarose gel. For *M. mycetomatis* specific PCR, a positive control (*M. mycetomatis* isolate previously identified by sequencing of ITS) that showed a band size of approximately 420 bp was used. For ITS PCR, 71 randomly selected amplicons were chosen for sequencing using primers ITS4 and ITS5. Standard Sanger sequencing was performed by Macrogen (Seoul, South Korea).

## Sequences alignment and strains identification

Sequence data were assembled using BioEdit 7.2.5 [26] and the identity was accomplished using BLAST in GenBank. To confirm the identity of the causative pathogens, we assessed the phylogeny of our sequences to sequences from three reference collections, i. e. the CBS-KNAW reference collection of the Westerdijk Fungal Biodiversity Institute, The Netherlands; the Institute Pasteur Collection of Fungi (UMIP), France; and the National Reference Center for Invasive Mycoses and Antifungals (NRCMA), France. The multiple-sequence alignments were generated using MUSCLE [27]. Phylogenetic trees were constructed with MEGA v7.0.18 using maximum likelihood approach with Kimura 2-parameter model [28].

## Statistics

Descriptive analysis was performed to determine the frequency of patients sex, ages, localities and infected sites as well as the frequency of identified *Madurella* species.

## Results

### Demographic and clinical data

In total, 138 patients were suspected to have mycetoma and underwent surgical intervention during the study period. The number of confirmed mycetoma cases was 119 (86.2%); of them 86 were male and 28 female, the ratio being 2.6:1. The most affected age group was from 11 to 30 years, while the overall range was from 1 to 71 years old. Affected body sites included: foot 70% (n = 83), hand 16% (n = 19), leg 4% (n = 5), knee 3% (n = 4), back 2.5%, (n = 3), bottom 2.5% (n = 3) and perianal 2% (n = 2) (Fig 2). Patients presented with different stages of the disease ranging from early and small lesions to massive lesions and bone involvement (Fig 3).
<figure><img src="image_5.png" /><figcaption>Fig 2. Frequency of sex, age groups and infected body sites among studied mycetoma patients.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0008420.g002

Bone involvement, mycetoma relapse and positive family history were recorded in 21% (n = 25), 12% (n = 14) and 12% (n = 14) of cases, respectively.

Black-grain mycetoma constituted 93.3% (n = 111/119) of the cases, 76.6% (n = 85/111) were from Gezira State and distributed among the eight localities. Highest number of cases was from South El Gezira: 30.6% (n = 26/85) and from Um Elqura: 27.1% (n = 23/85), while El Managil showed the lowest rate with 2.4% (n = 2/85). Black-grain mycetoma cases from Sinnar State were 23.4% (n = 26/111) and came from two localities: East Sinnar 76.9% (n = 20/26) and Sinnar 23.1% (n = 6/26) (Fig 4).

<figure><img src="image_6.png" /><figcaption>Fig 3. Clinical presentation of mycetoma cases. (A) Advance cases with multiple sinuses, (B) Recurrent mycetoma, (C) Bone involvement.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0008420.g003
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Frequency of black-grain mycetoma in Gezira and Sinnar States. Upper left: Gezira State localities, Lower right: Sinnar State localities.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0008420.g004

## Molecular identification

In the present study, a total of 94 black-grain specimens were processed for direct fungal DNA extraction following the CTAB protocol. Of them, 3 samples were lost during processing while 91 samples were maintained. The purity of the extracted DNA was evaluated and 75.8%
<figure><img src="image_8.png" /><figcaption>**Fig 5.** Gel representation of direct ITS and *Madurella mycetomatis*-specific PCR. (A) PCR of the ITS region using ITS4/ITS5 primers. A double band was noted in sample 4; (B) PCR of the ITS region. Lanes: 1 DNA ladder, 2 *M. mycetomatis*, 3 *M. tropicana*, 4 *M. fahalii*, 5 *Sphaerulina rhododendricola*, 6 and 7 are negative; (C) PCR using *Madurella mycetomatis*-Specific primers. Lanes: (2–6) *M. mycetomatis*, 7 *M. tropicana*, 8 *M. fahalii*.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0008420.g005

(n = 69) of the samples had 260/280 ratio between 1.6–2.0; for 15 (16.5%) samples the ratio was <1.6, while for 7 (7.7%) samples the ratio was >2. Time required for extraction including the homogenization step was approximately 1.5 hrs.

Of the 91 DNA specimens, 90.1% (n = 82) were successfully amplified with the ITS primers and yielded amplicons of approximately 600 bp. Negative specimen (n = 9) couldn't be amplified even after using other universal primer for ribosomal DNA amplification, namely V9D and LS266 [25]. Furthermore, one sample showed double DNA bands (Fig 5). By using *Madurella mycetomatis* species-specific primers, 79 specimens (86.8%) were successfully amplified and showed positive results with band sizes of approximately 420 bp (Fig 5). Only 3 specimens were positive with ITS PCR and negative with *Madurella mycetomatis* species-specific primers.

ITS sequencing was performed for 68 randomly selected samples, in addition to the three samples which showed positive ITS PCR but negative *M. mycetomatis* species-specific PCR. Sequences were deposited in GenBank database under the accession numbers shown in Fig 6.
<figure><img src="image_9.png" /></figure>
**Fig 6. Maximum likelihood tree of ITS sequences obtained from black-grain samples.** Blue square, samples from Gezira; red square, samples from Sinnar. Type strains are shown in bold.

https://doi.org/10.1371/journal.pntd.0008420.g006

Of the 71 samples, 67 sequences were obtained and analyzed, while for 4 samples no electro-
pherogram was obtained. With the BLAST search using ITS sequences, *M. mycetomatis* was
identified in 64 samples which is in 100% agreement with *M. mycetomatis* species-specific
PCR. Of the three *M. mycetomatis* species-specific PCR negative samples, one was identified as
*M. tropicana* and one as *M. fahalii*. One sample didn't match any of the known mycetoma
causative species, instead, it showed 99.37% similarity to the plant pathogen *Sphaerulina rho-
dodendricola* (ITS accession of the type strain, NR_137839).

An ITS alignment was constructed for 80 sequences with total length of 544 bp including
gaps. The species of *Madurella* were clearly separated into 4 clades each representing a single
species (Fig 6). The largest clade contained samples that showed identity to the type strain of
*M. mycetomatis*. No intraspecific genetic differences were observed in the ITS of *M. mycetoma-
tis* samples of the patients from different areas in Gezira State. In agreement with BLAST
results, we found one sample of *M. fahalii* and this represents the second record of *M. fahalii*
from Sudan. Besides the type strain of *M. fahalii*, other strains originated from Mali, British
Columbia, and Canada. We found the first Sudanese *M. tropicana* case in El Managil region of
Gezira State; other records of *M. tropicana* were from Indonesia, and Curacao. The sample
that matched *Sphaerulina rhododendricola*, *Capnodiales* was at a large distance from *Madur-
ella*, *Sordariales*.

Discussion

This study dealt with black-grain mycetoma patients operated at the Gezira Mycetoma Center (GMC) in the years 2015 to 2017. In the study group, male patients represented 72% of the total number and the most affected age group was between 11–30 years. In Sudan, a similar observation was made by Fahal *et al*. [18] in which 76% of the mycetoma patients were found to be male and 64% were under the age of 30 years. As mycetoma is an implantation mycosis, infection of males at this age group was previously explained by an increased risk of trauma during field work, and contact with soil and thorns in agricultural settings [29]. However, since both sexes might have a similar chance of contact with the causative pathogen, another hypothesis is, that males may express a higher susceptibility to the infection [30]. Although not proven for mycetoma, the endocrinological and immunological differences between male and female was found to be determinant for susceptibility to many infectious diseases including some fungal infections [30, 31, 32]. In the present study, we found that the most frequently infected body sites were foot 70% (n = 83) and hand 16% (n = 19). In addition, we encountered two cases of uncommon, perianal infected sites; one of these patients was a horse carriage driver, while the other was a laborer. The enigmatic route of infection emphasizes the need of taking more protective measures for workers in endemic regions [7].

Despite the large diversity of agents of black-grain eumycetoma, the present study showed that *M. mycetomatis* is responsible for 96% of the cases of eumycetoma in the Sudanese states Gezira and Sinnar. Its prevalence is comparable to observed frequencies of the species in the Mycetoma Reference Center in Khartoum, where *M. mycetomatis* was found to cause 70% of the infections; 4,754 out of 6,792 patients being infected by this fungus [18]. There is overlap in the data, as most of the cases treated in the Khartoum center (37%) were from Gezira State. In Sudan, Gezira, Sinnar, White Nile, Kordofan, and Khartoum States were the most affected regions according to data from Mycetoma Reference Center [18]. The high endemicity in these regions was explained by the ecological predilection of the causative agent in these
<figure><img src="image_10.png" /><figcaption>**Fig 7.** Sudan map, mycetoma endemic areas. A) Gezira State (Red), Sinnar State (Purple). B) Thorns and animal dung in an affected village in Gezira State.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0008420.g007

semiarid regions [33]. Although the exact niches of *Madurella* species remains as yet unknown, ecological niche modeling revealed a correlation between *Acacia* trees and mycetoma occurrence in Sudan [33]. However, testing *Acacia* thorns from the same endemic region for the presence of *Madurella* DNA yielded negative results [34]. The role of the thorny bushes of *Acacia* might be implantation of *Madurella* into the human body, rather than being a habitat for this fungus. The phylogenetic position of *Madurella* in *Chaetomiaceae* (order *Sordariales*) has generated another hypothesis, i.e. that *Madurella* might be a coprophilous fungus, with animal dung as a possible habitat [1,7]. This hypothesis was also supported by the fact that in villages endemic for mycetoma, animals are kept in close vicinity to the houses and in sheds surrounded by thorny bushes (Fig 7) [35]. The origins of the cases reported in the present study seemed to be in line with this hypothesis, as the highest number of the cases originated from South El Gezira (23.4%) and Um Elqura (20%), which are rural regions. Greater Wad Madani is located between these provinces but only 11 cases (10%) were recovered from there. This is probably due to urbanization, as Wad Madani is the most developed, and is the economic center of the Gezira State. Of the seven locations in Sinnar, only two harbored cases of eumycetoma, with East Sinnar showing highest number of cases (18%). Interestingly, 77.5% of overall cases (n = 86/111) were aggregated in a zone of five adjacent localities of the Gezira and Sinnar States (Fig 4). These five localities have similar environmental conditions, and there is extensive traffic between the states.

Until recently, all black-grain eumycetomata caused by *Madurella* were attributed to *M. mycetomatis*, as this was the only accepted species in the genus beside *M. grisea* which now has been transferred to *Trematosphaeria* as *T. grisea* [6, 8, 36]. However, molecular studies have revealed larger diversity in the genus with three additional species, viz. *M. fahalii*, *M. pseudo-mycetomatis*, and *M. tropicana* [12, 37]. *Madurella mycetomatis* remains the prevalent species causing mycetoma worldwide, while the others are rare [2]. Of *M. fahalii*, a single infection has been reported from Sudan [12], but we discovered several more cases caused by this species. The one from our study concerned a patient from South El Gezira, while other cases were observed in Mali, British Columbia, and Canada (Fig 6). Furthermore, two strains of *M. fahalii*
have been reported in patients from Saudi Arabia [38]. Given the origin of the patients infected
with *M. fahalii*, it seems that the epidemiology of the species indeed overlap with *M. mycetomatis*. One of our patients was infected with *M. tropicana*, a species which has as yet not been
reported from Sudan [12]. Thus, all *Madurella* species except *M. pseudomycetomatis* are pres-
ent in Sudan causing eumycetoma. Given the differences in antifungal susceptibility, identifi-
cation at the species level is recommended, which can be achieved by molecular methods [12].
By applying such methods, we have identified a new agent of black grain eumycetoma in
Sudan, *Sphaerulina rhododendricola*. Although agents of eumycetoma are distributed into at
least eight orders in the Ascomycota, there are no reports for species in the order *Capnodiales*.
*Sphaerulina rhododendricola* is the first Capnodialean fungus to be reported [39]. Further-
more, there might be other species from this order, but the use of methods with insufficient
specificity for diagnosis has limited the detection of such species.

Our study applies molecular, culture-free methods for the diagnosis of black-grain eumycetoma using a large sample size from hyperendemic regions in Sudan. K et al. [40], used direct PCR and sequencing for diagnosis of a culture-negative black-grain mycetoma case in India. One of the challenges of direct molecular detection is DNA extraction, as the grains may be recalcitrant to homogenization. Ahmed et al. used metal beads in a tissue lyser and a DNA extraction kit [22]. For less equipped laboratories in developing countries, we recommend cheapest extraction method with CTAB and homogenization by grinding [24]. To achieve a high DNA concentration of sufficient quality, grains should carefully be washed to remove dead tissue and inflammatory exudates. With ITS4 and ITS5, 90.1% (n = 82/91) of our samples were successfully amplified. Absence of amplification might be due to loss of DNA during purification, or to presence of PCR inhibitors. Although double bands were observed with ITS PCR in one of the samples, this sample proved to be positive with *M. mycetomatis* specific PCR. Of note, all our samples were from excised mycetoma biopsies during surgery, rather than from discharge through fistels, minimizing contamination. Nevertheless, evaluation of the same detection method using discharged grains is needed to verify whether the approach can be simplified abandoning invasive procedures.

In summary, we confirmed that the conclusive diagnosis of eumycetoma can be achieved
by direct detection of the causative pathogen from the grains using PCR and sequencing. The
method can be applied for routine diagnosis in reference centers such as the Gezira Mycetoma
Center. The identification of uncommon agents of eumycetoma in this study, such as *M. fahali*, *M. tropicana*, and *S. rhododendricola* highlighted the importance of species identifica-
tion in clinical settings. These species might differ in their susceptibility to antifungal agents
and their epidemiology thus specific identification is required for both patient's management
and to properly reveal the epidemiology of the species [12].

Author Contributions

**Conceptualization:** Elhadi A. Ahmed, Bakri Y. M. Nour.

**Data curation:** Elhadi A. Ahmed, Marie Desnos-Ollivier, Sybren de Hoog, Sarah A. Ahmed.

**Formal analysis:** Adam D. Abakar, Hisham N. Altayb, Sarah A. Ahmed.

**Investigation:** Elhadi A. Ahmed.

**Methodology:** Adam D. Abakar, Mogahid Mahmoud.

**Resources:** Mohamed Daffalla, Mogahid Mahmoud.

**Software:** Sarah A. Ahmed.

**Supervision:** Bakri Y. M. Nour, Adam D. Abakar.
**Validation:** Mohamed Daffalla.

**Visualization:** Hisham N. Altayb.

**Writing – original draft:** Elhadi A. Ahmed, Sarah A. Ahmed.

**Writing – review & editing:** Elhadi A. Ahmed, Marie Desnos-Ollivier, Sybren de Hoog.

## References

1. Ahmed S. New insights into a disfiguring fungal disease, eumycetoma. Ph.D. Thesis, University of Amsterdam. 2016. Available from: https://hdl.handle.net/11245/1.509518.

2. van de Sande WW. Global burden of human mycetoma: a systematic review and meta-analysis. PLoS Negl Trop Dis. 2013; 7: e2550. https://doi.org/10.1371/journal.pntd.0002550 PMID: 24244780

3. Desnos-Ollivier M, Bretagne S, Dromer F, Lortholary O, Dannaoui E. Molecular identification of black-grain mycetoma agents. J Clin Microbiol. 2006; 44: 3517–3523. https://doi.org/10.1128/JCM.00862-06 PMID: 17021076

4. de Hoog GS, Guarro J, Gené J, Ahmed S, Al-Hatmi AMS, Figueras MJ, Vitale RG. Atlas of clinical fungi. 4th web-ed. Westerdijk Institute, Utrecht and Universitat Rovira i Virgili, Reus: 2019.

5. van de Sande W, Fahal A, Ahmed SA, Serrano JA, Bonifaz A, Zijlstra E, et al. Closing the mycetoma knowledge gap. Med Mycol. 2018; 56: 153–164. https://doi.org/10.1093/mmy/myx061 PMID: 28992217

6. Ahmed SA, van de Sande WW, Stevens DA, Fahal A, van Diepeningen AD, Menken SB, et al. Revision of agents of black-grain eumycetoma in the order Pleosporales. Persoonia. 2014; 33: 141–154. https://doi.org/10.3767/003158514X684744 PMID: 25737597

7. de Hoog GS, Ahmed SA, Najafzadeh MJ, Sutton DA, Keisari MS, Fahal AH, et al. Phylogenetic findings suggest possible new habitat and routes of infection of human eumycetoma. PLoS Negl Trop Dis. 2013; 7: e2229. https://doi.org/10.1371/journal.pntd.0002229 PMID: 23696914

8. de Hoog GS, Adelmann D, Ahmed AO, van Belkum A. Phylogeny and typification of Madurella mycetomatis, with a comparison of other agents of eumycetoma. Mycoses. 2004; 47: 121–130. https://doi.org/10.1111/j.1439-0507.2004.00964.x PMID: 15078428

9. Fahal AS, Suliman SH. Clinical presentation of mycetoma. Sudan Med J. 1994; 32: 46–66.

10. van de Sande WW, Fahal AH, Goodfellow M, Mahgoub el S, Welsh O, Zijlstra EE. Merits and pitfalls of currently used diagnostic tools in mycetoma. PLoS Negl Trop Dis. 2014; 8: e2918. https://doi.org/10.1371/journal.pntd.0002918 PMID: 24992636

11. Borman AM, Desnos-Ollivier M, Campbell CK, Bridge PD, Dannaoui E, Johnson EM. Novel taxa associated with human fungal black-grain mycetomas: Emarellia grisea gen. nov., sp. nov., and Emarellia paragrisea sp. nov. J Clin Microbiol. 2016; 54: 1738–1745. https://doi.org/10.1128/JCM.00477-16 PMID: 27076666

12. de Hoog GS, van Diepeningen AD, el Mahgoub S, van de Sande WW. New species of Madurella, causative agents of black-grain mycetoma. J Clin Microbiol. 2012; 50: 988–994. https://doi.org/10.1128/JCM.05477-11 PMID: 22205798

13. Ahmed SA, Gonzalez GM, Tirado-Sanchez A, Moreno-Lopez LM, de Hoog S, Bonifaz A. Nigrograna mackinnonii, not Trematosphaeria grisea (syn., Madurella grisea), is the main agent of black grain eumycetoma in Latin America. J Clin Microbiol. 2018; 56. pii: e01723-17.

14. Ahmed SA, de Hoog GS, Stevens DA, Fahal AH, van de Sande WW. In vitro antifungal susceptibility of coelomycete agents of black grain eumycetoma to eight antifungals. Med Mycol. 2015; 53: 295–301. https://doi.org/10.1093/mmy/myu098 PMID: 25631481

15. Cerar D, Malallah YM, Howard SJ, Bowyer P, Denning DW. Isolation, identification and susceptibility of Pyrenochaeta romeroi in a case of eumycetoma of the foot in the UK. Int J Antimicrob Agents. 2009; 34: 617–618. https://doi.org/10.1016/j.ijantimicag.2009.08.004 PMID: 19783408

16. Los-Arcos I, Royuela M, Martin-Gomez MT, Alastruey-Izquierdo A, Sellares J, Perelló M, et al. Phaeo-hyphomycosis caused by Medicopsis romeroi in solid organ transplant recipients: Report of two cases and comprehensive review of the literature. Transpl Infect Dis. 2019; 21: e13072. https://doi.org/10.1111/tid.13072 PMID: 30865352

17. van de Sande WW, Luijendijk A, Ahmed AO, Bakker-Woudenberg IA, van Belkum A. Testing of the in vitro susceptibilities of Madurella mycetomatis to six antifungal agents by using the Sensititre system in comparison with a viability-based 2,3-bis(2-methoxy-4-nitro-5-sulfophenyl)-5-[(phenylamino)carbonyl]-2H-tetrazolium hydroxide (XTT) assay and a modified NCCLS method. Antimicrob Agents Chemother. 2005; 49: 1364–1368. https://doi.org/10.1128/AAC.49.4.1364-1368.2005 PMID: 15793113
18. Fahal A, Mahgoub E, El Hassan AM, Abdel-Rahman ME. Mycetoma in the Sudan: An update from the Mycetoma Research Centre, University of Khartoum, Sudan. PLoS Negl Trop Dis. 2015; 9: e0003679. https://doi.org/10.1371/journal.pntd.0003679 PMID: 25816316

19. van Belkum A, Fahal A, van de Sande WW. Mycetoma caused by Madurella mycetomatis: a completely neglected medico-social dilemma. Adv Exp Med Biol. 2013; 764: 179–189. https://doi.org/10.1007/978-1-4614-4726-9_15 PMID: 23654067

20. Gismalla MDA, Ahmed GMA, Mohamed Ali MM, Taha SM, Mohamed TA, Ahmed AE, et al. Surgical management of eumycetoma: experience from Gezira Mycetoma Center, Sudan. Trop Med Health. 2019; 47: 6. https://doi.org/10.1186/s41182-018-0129-2 PMID: 30675125

21. Hazra N, Gupta RM, Verma R, et al. Black grain eumycetoma revisited: enhanced fungal isolation by use of a novel technique. J Bacteriol Mycol Open Access. 2017; 4:159–163.

22. Ahmed SA, van de Sande WW, Desnos-Ollivier M, Fahal AH, Mhmoud NA, de Hoog GS. Application of isothermal amplification techniques for identification of Madurella mycetomatis, the prevalent agent of human mycetoma. J Clin Microbiol. 2015; 53: 3280–3285. https://doi.org/10.1128/JCM.01544-15 PMID: 26246484

23. Ahmed AO, Mukhtar MM, Kools-Sijmons M, Fahal AH, de Hoog S, Gerrits van den Ende B, et al. Development of a species-specific PCR-restriction fragment length polymorphism analysis procedure for identification of Madurella mycetomatis. J Clin Microbiol. 1999; 37: 3175–3178. https://doi.org/10.1128/JCM.37.10.3175-3178.1999 PMID: 10488173

24. Moller EM, Bahnweg G, Sandermann H, Geiger HH. A simple and efficient protocol for isolation of high-molecular-weight DNA from filamentous fungi, fruit bodies, and infected-plant tissues. Nucleic Acids Res. 1992; 20: 6115–6116. https://doi.org/10.1093/nar/20.22.6115 PMID: 1461751

25. White T J, Bruns T, Lee S, Taylor J W. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In Innis M. A., Gelfand D. H., Sninsky J. J., and White T. J., editors. PCR protocols: a guide to methods and applications. New York Academic Press; 1990. pp. 315–322.

26. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Window 95/98/NT. Nucleic Acids Symp Ser. 1999; 41: 95–98.

27. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32: 1792–1797. https://doi.org/10.1093/nar/gkh340 PMID: 15034147

28. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016; 33: 1870–1874. https://doi.org/10.1093/molbev/msw054 PMID: 27004904

29. Queiroz-Telles F, Fahal AH, Falci DR, Caceres DH, Chiller T, Pasqualotto AC. Neglected endemic mycoses. Lancet Infect Dis. 2017; 17: E367–E377. https://doi.org/10.1016/S1473-3099(17)30306-7 PMID: 28774696

30. van de Sande WW, Fahal A., Tavakol M, van Belkum A. Polymorphisms in catechol-O-methyltransferase and cytochrome p450 subfamily 19 genes predispose towards Madurella mycetomatis-induced mycetoma susceptibility. Med Mycol. 2010; 48: 959–968. https://doi.org/10.3109/13693781003636680 PMID: 20184498

31. Giefing-Kroll C, Berger P, Lepperdinger G, Grubeck-Loebenstein B. How sex and age affect immune responses, susceptibility to infections, and response to vaccination. Aging Cell. 2015; 14: 309–321. https://doi.org/10.1111/acel.12326 PMID: 25720438

32. Guess TE, Rosen JA, McClelland EE. An overview of sex bias in C. neoformans infections. J Fungi (Basel). 2018; 4. pii: E49.

33. Samy AM, van de Sande WW, Fahal AH, Peterson AT. Mapping the potential risk of mycetoma infection in Sudan and South Sudan using ecological niche modeling. PLoS Negl Trop Dis. 2014; 8: e3250. https://doi.org/10.1371/journal.pntd.0003250 PMID: 25330098

34. Ahmed A, Adelmann D, Fahal A, Verbrugh H, van Belkum A, de Hoog S. Environmental occurrence of Madurella mycetomatis, the major agent of human eumycetoma in Sudan. J Clin Microbiol. 2002; 40: 1031–1036. https://doi.org/10.1128/jcm.40.3.1031-1036.2002 PMID: 11880433

35. Kulneff CA. A comparative study of urban rural dairy management systems in Sudan. Ph.D. Thesis, Swedish University of Agricultural Science. 2006. Available from: http://epsilon.slu.se/10092933.pdf.

36. Ahmed AO, van Leeuwen W, Fahal A, van de Sande W, Verbrugh H, van Belkum A. Mycetoma caused by Madurella mycetomatis: a neglected infectious burden. Lancet Infect Dis. 2004; 4: 566–574. https://doi.org/10.1016/S1473-3099(04)01131-4 PMID: 15336224

37. Yan J, Deng J, Zhou CJ, Zhong BY, Hao F. Phenotypic and molecular characterization of Madurella pseudomycetomatis sp. nov., a novel opportunistic fungus possibly causing black-grain mycetoma. J Clin Microbiol. 2010; 48: 251–257. https://doi.org/10.1128/JCM.00018-09 PMID: 19923486
38. Fraser M, Borman AM, Johnson EM. Rapid and robust identification of the agents of black-grain myce-toma by Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry. J Clin Micro-biol. 2017; 55: 2521–2528. https://doi.org/10.1128/JCM.00417-17 PMID: 28592546

39. Rathnayaka AR and Jayawardena RS. Checklist of order Capnodiales in Thailand. Plant Pthol. & Quar. 2019; 9:166–184.

40. K S, Das S, Pandhi D, Rai G, Ansari MA, Gupta C, et al. Challenges in culture-negative cases of Madur-ella mycetomatis: A case report re-accentuating PCR as an essential diagnostic tool. J Mycol Med. 2017; 27:577–581. https://doi.org/10.1016/j.mycmed.2017.09.004 PMID: 29102309