Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Molecular epidemiology of human sporotrichosis in Venezuela reveals high frequency of Sporothrix globosa

  • Emma Camacho1,
  • Isabel León-Navarro1,
  • Sabrina Rodríguez-Brito1,
  • Mireya Mendoza2 and
  • Gustavo A Niño-Vega1Email author
BMC Infectious Diseases201515:94

https://doi.org/10.1186/s12879-015-0839-6

Received: 18 November 2014

Accepted: 16 February 2015

Published: 25 February 2015

Abstract

Background

Sporotrichosis is a cutaneous and subcutaneous fungal disease of humans and other mammals, known to be caused by the Sporothrix schenckii species complex, which comprises four species of clinical importance: S. brasiliensis, S. globosa, S. luriei, and S. schenckii sensu stricto. Of them, S. globosa and S. schenckii s. str. show global distribution and differences in global frequency as causal agents of the disease. In the Americas, only three species are present: S. schenckii s. str., S. brasiliensis (so far, only reported in Brazil), and S. globosa. In Venezuela, since the first case of sporotrichosis reported in 1935, S. schenckii have been considered its unique etiological agent. In the present work, the presence of more than one species in the country was evaluated.

Methods

By phenotypic key features and molecular phylogeny analyses, we re-examined 30 isolates from diverse Venezuelan regions belonging to the fungi collection of Instituto de Biomedicina, Caracas, Venezuela, and national reference center for skin diseases. All isolates were collected between 1973 and 2013, and maintained in distilled water.

Results

Sporotrichosis in Venezuela is mainly caused by S. schenckii s. str. (70%). However, a significant proportion (30%) of sporotrichosis cases in the country can be attributable to S. globosa. A correlation between intraspecific genotypes and clinical presentation is proposed.

Conclusions

Our data suggest that sporotrichosis various clinical forms might be related to genetic diversity of isolates, and possibly, to diverse virulence profiles previously reported in the S. schenckii species complex. Sporothrix globosa was found to be the causative agent of 30% of sporotrichosis for the Venezuelan cases re-examined, the highest frequency of this species so far reported in the Americas. The high genetic variability presented by S. schenckii s. str. indicates that species distinction based on phenotypic key features could be a challenging and uncertain task; molecular identification should be always employed.

Keywords

Sporothrix schenckii species complex Sporothrix globosa Sporotrichosis Clinical forms Phenotypic identification Phylogenetic analysis American continent Venezuela

Background

Over the last few years, the Sporothrix schenckii species complex (formerly known as S. schenckii) has received special interest due to the increasing number of infections caused by it worldwide, particularly in immunocompromised patients [1], and epidemic outbreaks in cats in Brazil [2]. It is responsible for sporotrichosis, a chronic granulomatous subcutaneous mycosis of humans and other mammals. Primary infection has been reported as caused by traumatic inoculation of environmental material carrying saprophytic hyphae and conidia of the etiologic agents [2,3]. Once into the warm-blooded host, these dimorphic fungi are capable to convert into pathogenic yeasts, causing infections that range from fixed cutaneous localized lesions to severe disseminated sporotrichosis [3-5]. However, zoonotic transmission is also possible, by scratches or bites from asymptomatic or infected animals, with cats being the main vectors through which the disease is transmitted to humans or other animals [6-8]. The last makes direct inoculation of pathogenic yeast a plausible source of infection. Although the disease has a global occurrence, endemic areas are mostly located in tropical and subtropical countries [6,9]; furthermore, Sporothrix infections may take epidemic proportions [7,8,10] and its distinct etiological agents differ in virulence profiles [11-13], antifungal susceptibility [14,15], and geographic distribution [16].

Molecular phylogenetic analyses have led to description of at least four cryptic species of clinical relevance within the S. schenckii species complex, comprising S. brasiliensis, S. schenckii sensu stricto (s. str), S. globosa and S. luriei [2,8,17-19]. Scattered reports of Sporothrix spp. outside these four clades, causing clinical cases have been published, such as S. mexicana [16,20] and S. pallida [21], or by close relatives in the genus Ophiostoma, i.e. O. piceae [22] and O. stenoceras [23]. However, these species appear to lack human-pathogenic potential, thus they are proposed to be accidentally pathogenic [17]. Besides molecular phylogeny, morphological and physiological key features have been proposed for species recognition within the Sporothrix species complex [18,19].

In Venezuela, sporotrichosis is the second most common subcutaneous mycosis after chromomycosis. It was first described in 1935, nevertheless the exact prevalence in the country is unknown [24]. According to conventional mycological procedures, as well as epidemiological data, all isolates related to sporotrichosis have been identified as Sporothrix schenckii.

The aim of the present study was to determine if, besides S. schenckii, other Sporothrix species of the complex might be present in the country. We re-examined 30 isolates, by phenotypic and molecular methods, discovering that S. globosa was the etiological agent in one third of the cases.

Methods

Fungal isolates and strains used

Thirty isolates (29 clinical and 1 environmental) from different Venezuelan regions (Coastal range; n = 22, The Andes; n = 7; The Plains; n = 2) were examined in this study (Table 1). All isolates were previously identified as S. schenckii by means of morphology (macro and microscopic studies), ability of isolates to reverse to yeast-like cells at 37°C, and serological tests to patients from where they were isolated, all performed in the Mycology Laboratory at Instituto de Biomedicina, Caracas, Venezuela, a national reference center for skin diseases. Isolates were taken as part of standard patient care, and no ethical approval was required for their use. They have been kept as part of the laboratory fungal collection over a period of 40 years (from 1973 to 2013) and were selected for inclusion in this study based on their geographic distribution. They were preserved in distilled water (Castellani’s method), and recovered by growth on Sabouraud dextrose agar (SDA) complemented with 150 μg ml−1chloramphenicol at room temperature during 7 days. As reference, strains S. schenckii s. str. ATCC-MYA 4820, and S. brasiliensis ATCC-MYA 4823 (provided by Dr. L. Bezerra, Universidade do Estado do Rio de Janeiro, Brazil), as well as S. pallida CBS 302.73T, S. globosa FMR 9023, and S. mexicana FMR 9108 kindly provided by Dr. J. Cano (Reus, Spain) were used.
Table 1

Isolates used in this study

Isolate

Species (by molecular genotyping)

Source

Origin

GenBank Accession N°

CAL

ITS

C0014

S. schenckii

Clinical, Human (lymphocutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478909

1KJ999877

C0329

S. globosa

Clinical, Human (fixed cutaneous)

The Andes (Táchira, Venezuela)

1KF478892

1KJ999878

C0596

S. schenckii

Clinical, Human (fixed cutaneous)

Coastal range (Miranda, Venezuela

1KF478896

1KJ999879

C1163

S. schenckii

Clinical, Human (lymphocutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478910

1KJ999880

C2606

S. schenckii

NA

Coastal range (Distrito Capital, Venezuela)

1KF478887

1KJ999881

C2656

S. schenckii

Clinical, Human (fixed cutaneous)

The Andes (Trujillo, Venezuela)

1KF478897

1KJ999882

C2745

S. schenckii

Clinical, Human (fixed cutaneous)

The Plains (Monagas, Venezuela)

1KF478891

1KJ999883

C2888

S. schenckii

Clinical, Human (lymphocutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478886

1KJ999884

C3037

S. schenckii

Clinical, Human (disseminated)

Coastal range (Distrito Capital, Venezuela)

1KF478888

1KJ999885

C3527

S. schenckii

NA

Coastal range (Distrito Capital, Venezuela)

1KF478898

1KJ999886

C3673

S. globosa

Clinical, Human (fixed cutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478892

1KJ999887

C3705

S. schenckii

Clinical, Human (lymphocutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478884

1KJ999888

C4346

S. schenckii

Clinical, Human (lymphocutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478889

1KJ999889

C5000

S. schenckii

Clinical, Human (lymphocutaneous)

Coastal range (Aragua, Venezuela)

1KF478885

1KJ999890

C5859

S. globosa

Clinical, Human (fixed cutaneous)

Coastal range (Vargas, Venezuela)

1KF478893

1KJ999891

C7012

S. schenckii

Clinical, Human (lymphocutaneous)

The Andes (Táchira, Venezuela)

1KF478899

1KJ999892

C8213

O. stenoceras

Clinical, Human (fixed cutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478913

1KJ999893

C8287

S. schenckii

Clinical, Human (fixed cutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478900

1KJ999894

C8704

S. schenckii

NA

Coastal range (Distrito Capital, Venezuela)

1KF478890

1KJ999895

C8775

S. globosa

Clinical, Human (fixed cutaneous)

Coastal range (Aragua, Venezuela)

1KF478901

1KJ999896

C8857

S. globosa

Clinical, Human (fixed cutaneous)

The Andes (Trujillo, Venezuela)

1KF478902

1KJ999897

C8888

S. schenckii

Clinical, Human (lymphocutaneous)

The Andes (Trujillo, Venezuela)

1KF478903

1KJ999898

C8962

S. globosa

Clinical, Human (fixed cutaneous)

The Plains (Cojedes, Venezuela)

1KF478904

1KJ999899

C8981

S. globosa

Clinical, Human (fixed cutaneous)

The Andes (Mérida, Venezuela)

1KF478905

1KJ999900

C9254

S. globosa

Clinical, Human (fixed cutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478911

1KJ999901

C9300

S. schenckii

Clinical, Human (lymphocutaneous)

The Andes (Trujillo, Venezuela)

1KF478906

1KJ999902

C9526

S. schenckii

NA

Coastal range (Distrito Capital, Venezuela)

1KF478907

1KJ999903

C9862

S. schenckii

Clinical, Human (fixed cutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478912

1KJ999904

C9887

S. globosa

Clinical, Human (fixed cutaneous)

Coastal range (Distrito Capital, Venezuela)

1KF478908

1KJ999905

A0001

S. schenckii

Environmental, Soil

Coastal range (Aragua, Venezuela)

1KF478894

1KJ999876

ATCC MYA-4820

S. schenckii

Clinical, Human

Brazil

2 JF313361

3 JQ070111

ATCC MYA-4823

S. brasiliensis

Feline sporothricosis

Brazil

2 JF313351

3 JQ070114

CBS 302.73 T

S. pallida

Environmental

United Kingdom

4 AM398396

1 KJ999906

FMR 9023

S. globosa

Clinical, Human

Japan

4 AM398393

1 KJ999907

FMR 9108

S. mexicana

Environmental

Mexico

4 AM398393

5 FN549906

Previously genotyped strains used as reference, in bold.

1This study; 2Teixeira, M. and Lopes-Bezerra, L.M., direct submission; 3Suh, S.-O. and Zhou, J.J, direct submission; 4ATCC Mycology Authentication Project.

5Madrid et al. [33]. NA, data non-available.

Phenotypic characterization

Isolates were phenotypically identified according to Marimon et al. [18]. Macroscopic features of colonies were studied by culturing isolates on potato dextrose agar (PDA, HIMEDIA, India) plates, incubated at room temperature (growth control), 30°C and 37°C in dark. Petri dishes (90 × 15 mm, with 8 mm PDA) were centrally inoculated with portions of the colonies of the fungi, approximately 1 mm in diameter, and incubated upside down. After 14 and 21 days, colony diameters (in millimeters) were measured in duplicate for growth temperatures of 30°C and 37°C, and the mean of the recorded measured diameters taken for final analysis. Microscopic features of conidia were found to be the same for reference strains and a sample of Venezuelan isolates when lactrimel agar (LA) was used instead of the corn meal agar (CMA) used by Marimon et al. [18]. So for the rest of the study, microscopic features of conidia were determined from slide cultures made on lactrimel agar (10 g honey, 20 g wheat flour, 200 ml milk, 9.5 g agar, 800 ml distilled water and 150 μg ml−1 chloramphenicol) after 12–15 days of incubation in the dark in a humid chamber. Coverslips were mounted in lactophenol cotton blue (0.05% (w/v)), 2% phenol, 2% lactic acid and 4% glycerin), and examined under a light microscope (Leica DM2000) at 40X and 100X. Carbohydrate assimilation tests for sucrose and raffinose were conducted in 96-well microplates, using 150-μl of a working solution composed of yeast nitrogen base (YNB) medium (Sigma-Aldrich, USA) with corresponding test sugar, or plain (only base medium), as described by Marimon et al. [18]. Briefly, the isolates were grown on potato dextrose broth for 7 days at 23°C, or on Brain-Heart Infusion agar (BHI, OXOID, England), at 30°C for isolates that did not sporulated with the previous conditions. Conidia were collected by filtering 10 to 20 ml liquid culture through sterile gauze and quantified in a spectrophotometer at 600 nm. Each suspension was adjusted to an optical density that ranged from 0.21 to 0.29, corresponding to a final inoculum in the microplate of 2×105 to 2×106 CFU/ml. Each well was inoculated with 50 μl inoculum and 150 μl YNB medium with corresponding carbohydrate, at a final concentration of 3.75%. YNB medium without sugar was used as negative control, and with glucose as positive control. Microplates were read at 5 and 10 days of incubation at 25°C, experiments were conducted by duplicate. The viability of the conidia was verified by plating 20 μl of a 1:500 dilution of each inoculum onto YPD (0.5% yeast, 0.5% peptone, 1.5% dextrose) plates. In case of discordant results, experiments were repeated at least two additional times on different days.

Molecular characterization

Total genomic DNA was isolated from 3-day-old mycelial cultures grown on YPD broth, following the protocol described in [25]. Partial amplification of the calmodulin (CAL) locus was obtained by using degenerated primers CL1 and CL2A [26], amplifying a region corresponding to exons 3 through 5 of the CAL gene. ITS regions were amplified from genomic DNA with universal primers ITS5 and ITS4 as described by Zhou et al. [17]. Following PCR, amplicons were purified with the CONCERT™ Rapid PCR purification system (Life Technologies, USA) and sent for sequencing on both strands to Macrogen Inc (Seoul, Korea), with the same primers used for PCR. Assembly of the sequences was done with the Contig program within the Vector NTI suite (Vector NTI, InforMax, Inc, USA). Homology searches were performed on GenBank database using BLAST [27]. The newly reported sequences generated in this study were deposited in the GenBank/EMBL/DDBJ database under accession numbers listed in Table 1.

Phylogenetic analysis

Sporothrix-calmudulin-related and Sporothrix-ITS-related sequences, previously reported to be obtained with the same two sets of primers here used, were retrieved from Genbank and included in the analyses (Figures 1 and 2) [16,18]. As an outgroup, the saprophytic fungus Grosmannia serpens (Ophiostomataceae) [28] was included for the CAL analysis. The multiple nucleotide sequence alignment was performed using the ClustalW algorithm implemented in MEGA5.2 software [29]. Evolutionary analyses were also conducted in MEGA5.2 as described by Rodrigues et al. [16].
Figure 1

Phylogenetic relationships ofSporothrix schenckiicomplex isolates, inferred from CAL sequences. The phylogenetic tree was built with MEGA 5.1, by the Maximum Likelihood method based on the Kimura 2-parameter model. Thirty sequences from this study and sixty retrieved from previous studies were used for analysis. Number close to branches represents Bootstraps support values based on 1000 bootstrap replications. GenBank accessions numbers are indicated next to strain code in branch labels.

Figure 2

Phylogenetic relationships ofSporothrix schenckiicomplex isolates inferred from ITS sequences. The phylogenetic tree was built with MEGA 5.1, by the Maximum Likelihood method based on the Kimura 2-parameter model. Thirty sequences from this study and forty four from previous reports were used for analysis. Number close to branches represents Bootstraps support values based on 1000 bootstrap replications. GenBank accessions numbers are indicated next to strain code.

Genetic variation analysis

To evaluate genetic variation, the CAL locus was used as nuclear marker. Nucleotide (π) and haplotype (Hd) diversities [30] were estimated using DnaSP software, version 5.10.01 [31]. Sites containing gaps and missing data were not considered in the analysis.

Results

A total of 30 previously identified S. schenckii Venezuelan isolates (29 clinical and 1 environmental) (Table 1) were re-examined on the basis of phenotypic key features (i.e. growth at various temperatures, macroscopic and microscopic features, and carbohydrate assimilation, according to Marimon et al. [18]), as well as molecular phylogeny inferred from CAL and ITS sequences. Based on the result obtained for reference strain S. globosa FMR 9023 (Table 2), restricted growth at 37°C was set at ≤ 6.4 mm colony diameter (maximum diameter reached in our study by this reference strain, previously identified by Marimon et al. [18]). Only 8 out of 30 isolates (26.67%) were clearly identified when the summary of key features for Sporothrix species differentiation, proposed by Marimon et al. [18], was strictly followed (Table 2).
Table 2

Isolates characterized by phenotypic and molecular methods included in the study

Isolate

Mean colony diameter (mm) ± SD on PDA during 21 days

Assimilation test result

Morphology of dark sessile conidia

Identification on physiological tets

Identification on CAL and ITS sequences

30°C

37°C

Sucrose

Raffinose

C0014

45.5 ± 0.7

12.0 ± 0.0

+

-

Triangular

NC

S. schenckii

C0329

46.5 ± 2.1

6.5 ± 0.7

+

-

Subglobose

NC

S. globosa

C0596

53.5 ± 0.7

23.0 ± 0.0

+

-

Triangular

NC

S. schenckii

C1163

52.5 ± 2.1

11.5 ± 0.7

+

+

Triangular

NC

S. schenckii

C2606

56.0 ± 0.0

11.5 ± 0.7

+

-

Triangular

NC

S. schenckii

C2656

52.5 ± 0.7

13.5 ± 0.7

+

-

Triangular

NC

S. schenckii

C2745

45.5 ± 0.7

10.0 ± 0.0

+

-

Triangular

NC

S. schenckii

C2888

53.5 ± 0.7

12.5 ± 0.7

+

-

Triangular

NC

S. schenckii

C3037

57.0 ± 1.4

15.0 ± 0.0

+

+

Triangular

NC

S. schenckii

C3527

54.0 ± 1.4

17.0 ± 1.4

+

+

Triangular

NC

S. schenckii

C3673

41.5 ± 0.7

5.5 ± 0.7

+

-

Globose

S. globosa

S. globosa

C3705

50.0 ± 1.4

17.0 ± 1.4

+

+

Triangular

S. schenckii

S. schenckii

C4346

52.0 ± 1.4

15.0 ± 0.0

+

+

Triangular

NC

S. schenckii

C5000

55.5 ± 0.7

16.5 ± 0.7

+

+

Triangular

NC

S. schenckii

C5859

38.0 ± 0.0

4.0 ± 0.0

+

-

Globose

S. globosa

S. globosa

C7012

58.0 ± 0.0

22.0 ± 0.0

+

+

Triangular

NC

S. schenckii

†C8213

25.5 ± 4.9

12.7 ± 3.1

+

-

NP

NC

O. stenoceras

†C8287

30.0 ± 0.0

3.0 ± 0.0

+

-

Triangular

S. globosa

S. schenckii

C8704

44.5 ± 0.7

7.0 ± 1.4

+

-

Triangular

NC

S. schenckii

†C8775

23.0 ± 0.0

3.0 ± 2.8

-

-

Subglobose

NC

S. globosa

C8857

33.5 ± 0.7

3.0 ± 2.8

+

-

Subglobose

S. globosa

S. globosa

C8888

52.5 ± 0.7

14.5 ± 0.7

+

-

Triangular

NC

S. schenckii

C8962

41.0 ± 0.0

5.5 ± 0.7

+

-

Subglobose

S. globosa

S. globosa

C8981

38.0 ± 0.0

3.0 ± 0.0

+

-

Globose

S. globosa

S. globosa

C9254

44.5 ± 0.7

8.0 ± 0.0

+

-

Globose

NC

S. globosa

C9300

44.5 ± 0.7

13.5 ± 0.7

+

-

Triangular

NC

S. schenckii

C9526

42.0 ± 0.0

13.5 ± 0.7

+

+

Triangular

S. schenckii

S. schenckii

C9862

52.5 ± 2.1

12.0 ± 1.4

+

-

Triangular

NC

S. schenckii

C9887

37.5 ± 0.7

3.5 ± 0.7

+

-

Globose

S. globosa

S. globosa

A0001

58.5 ± 0.7

15.5 ± 0.7

+

-

Triangular

NC

S. schenckii

ATCC MYA-4820

41.0 ±  0.0

13.0  ±  1.4

+

+

Triangular

S. schenckii

S. schenckii

ATCC MYA-4823

43.5 ±  0.7

13.0  ±  1.4

-

-

Globose

S. brasiliensis

S. brasiliensis

CBS 302.73 T

76.5 ±  0.7

11.0  ±  0.0

+

-

NP

S. pallida

S. pallida

FMR 9023

33.5  ±  0.7

5.0  ±  1.4

+

-

Subglobose

S. globosa

S. globosa

FMR 9108

73.0  ±  1.4

16.5  ±  0.7

+

+

Ellipsoidal

S. mexicana

S. mexicana

Reference strains in bold.

NC, non-conclusive; NP, not present; †Isolates showing atypical phenotypic profiles.

One isolate, S. schenckii C8287, was misidentified by following the phenotypic features solely, requiring molecular genotyping for proper identification; it was labeled as S. globosa due to its restrictive growth at 37°C, and positive/negative assimilation profile for sucrose and raffinose. The remaining 21 isolates (70.00%) showed inconclusive results, requiring molecular genotyping for proper identification (Table 2). The most ambiguous results (colony growth over 50 mm diameter at 30°C and/or negative assimilation of raffinose) were shown by S. schenckii s. str. isolates, for which just 2 out of 20 (10%) were properly identified without the need of molecular phylogeny analysis. However, it is important to notice that when the microscopic features of sessile conidia (pigmentation and morphology) were taken solely, correlation between identification by these features and molecular identification was found (Table 2).

Genotyping of the isolates was performed by using the amplified fragment from the CAL locus. The complete alignment included 90 sequences, 30 generated in this study and 60 retrieved from previous studies, respectively [16,18]. The aligned CAL sequences were 595 bp long, including 302 invariable characters, 239 variable parsimony-informative sites (40.17%) and 43 singletons. Twenty nine (29) of the 30 Venezuelan isolates clustered within the clinical clades (as defined by previous phylogenetic analyses of Sporothrix spp. based on CAL [6,16,32]). Of those 29 isolates, 20 (>60%) grouped as S. schenckii s. str., while 9 (30%) grouped as S. globosa (Figure 1). One isolate (C8213), clustered within the environmental clades as O. stenoceras (Figure 1). To corroborate these results, a second phylogenetic analysis was performed, using ITS sequences generated by the universal primers ITS4 and ITS5. Seventy four (74) nucleotide sequences, 30 generated in this study and 44 previously reported [17] were aligned. The aligned ITS sequences were 506 bp long, including 407 invariable characters, 85 variable parsimony-informative sites (16.79%) and 12 singletons. The Venezuelan isolates were distributed in complete agreement with the CAL-generated phylogeny, with the Venezuelan isolates of S. schenckii s. str. and S. globosa (pathogenic species), and O. stenoceras (accidentally pathogenic species) showing strong, statistically supported separation (99% bootstrap) between clinical and environmental clades (Figure 2).

Assessment of the Venezuelan Sporothrix population genetic diversity was explored by using the DnaSP software [31] and CAL as nuclear marker. The haplotype analysis of calmodulin sequences (n = 65) divided the isolates into 17 Hap groups. A total of 15 and 5 different types were detected for S. schenckii s. str. and S. globosa, respectively. The majority of haplotypes (Hd = 0.85) belonged to S. schenckii s. str., demonstrating a highly diverse group (π = 0.020). Meanwhile, S. globosa presented low genetic variation (Hd = 0.59; π = 0.005).

According to epidemiological data of clinical cases for the Venezuelan isolates, most of them were geographically restricted to the coastal range region (68.97%). An analysis of the distribution by gender showed a higher occurrence in male patients (64%) compared to women (36%). The most predominant clinical form in this study was fixed cutaneous (60%), followed by lymphocutaneous (36%) and one case of disseminated sporotrichosis (4%).

An association between clinical forms and phylogenetic species was found (Figure 3). All lymphocutaneous cases (100%) were caused by S. schenckii s. str., while fixed cutaneous cases were attributable to both, S. schenckii s. str. (33.3%) and S. globosa (60.0%), plus a single case by O. stenoceras. The only case with disseminated manifestations was related to S. schenckii s. str.
Figure 3

Frequency of clinical forms of sporotrichosis vs. phylogenetic species. DS, disseminated; LC, lymphocutaneous; FC, fixed cutaneous. Information corresponds to clinical cases with available epidemiological data (Table 1).

Discussion

Even though sporotrichosis is a cosmopolitan subcutaneous mycosis, the species grouped in the clinical clades have shown a geographical distribution [17]. S. schenckii s. str. presents a frequency over 50% in America, while S. globosa has prevailed in Asia and Europe (56% and 28%, respectively), followed by the Americas (11%) [17]. In the present study, thirty fungal isolates related to sporotrichosis in Venezuela, and previously identified by morphological studies as S. schenckii, were characterized by phenotypic key features and molecular phylogeny. The results show that in the country, at least two species are circulating: S. schenckii s. str. and S. globosa. Although S. schenckii s. str. is the main circulating species, with over 60% of evaluated cases, S. globosa could be considered the second most frequent species causing sporotrichosis in Venezuela. The incidence of S. globosa, with 30% of cases evaluated, is the highest reported so far in the Americas for this species, and is similar to its reported incidence in Europe [17]. An interesting finding is the fact that all S. globosa isolates here characterized, were related to fixed cutaneous sporotrichosis, which could be related to the low virulence reported for this species [33,34]. On the other hand, 100% of the lymphocutaneous cases, and the single disseminated case were all related to S. schenckii s. str., as well as 33% of fixed cutaneous sporotrichosis.

In a recent study [13], the virulence profiles of eight S. schenckii s. str. were evaluated, and different degrees of virulence were found for this species, from high virulence to non-virulence. High genetic variation have been reported for S. schenckii s. str. [15,33,35]. Here, we also found high genetic variation for the Venezuelan S. schenckii s. str. isolates evaluated (Hd = 0.85; π = 0.020), while low genetic variation was obtained for the S. globosa isolates (Hd = 0.59; π = 0.005). Summing up the preceding results, we could hypothesize that a relationship might be present between the high genetic variation, different clinical presentations, and multiple virulence levels found for S. schenckii s. str., while the low genetic variation of S. globosa isolates might be related to a single clinical presentation (fixed cutaneous sporotrichosis), and low virulence profiles related to this species [11,13]. However, in Northeast China, clinical cases of sporotrichosis (from fixed cutaneous to disseminated forms), have been attributable exclusively to S. globosa [36]. To confirm or discard our hypothesis, further studies are required, evaluating whether correlations can be found between virulence profiles, species, patient’s immunological conditions, and clinical presentation of isolate’s source.

Earlier characterization reports of Sporothrix spp. comparing the phenotypic key features for species identification proposed by Marimon et al. [18], and molecular phylogeny methodologies, have shown few disagreements for species identification [13,16,32,37]. In the present study, a high proportion (70%) of Venezuelan isolates related to sporotrichosis, were not identified according to their phenotype key features (as proposed by Marimon et al. [18]). We discarded a technical problem, since all reference strains were accurately identified by this method (Table 2). So, at least for the Venezuelan isolates here studied, the phenotypic characterization was not useful for species identification.

In our molecular phylogenetic analyses based on CAL and ITS sequences, a strong separation between clinical and environmental clades was evident (Figures 1 and 2), in correspondence with previous reports [6,16-18,33]. Isolate C8213, identified as O. stenoceras was isolated from an immunocompetent patient, supporting the idea that mammal-pathogenicity of Ophiostomales, outside the S. schenckii spp. complex, although highly exceptional [17], still can occur. To the best of our knowledge, this is the first report of O. stenoceras misidentified as causative agent of sporotrichosis in the Americas. Further studies into the mechanisms of opportunistic fungal pathogenesis would help to understand why a harmless fungus, as O. stenoceras, could become an infective agent, even for immunocompetent patients. This becomes even more relevant, when we consider the increasing number of reports of Sporothrix infection in immunocompromised patients, mainly the HIV-infected population [1].

Even when selection of the isolates here studied was made trying to cover a broad geographical area of the country, most isolates originated from clinical cases of the coastal region, which can be explained, besides it being a possible ecological niche for the species, to the fact that although this area covers only 3% of the total surface of the country, it is also the most heavily populated region. Further work with larger data on molecular epidemiology of sporotrichosis in Venezuela is required to verify whether correlation between clinical forms and phylogenetic species exists, and gain better knowledge about S. schenckii species complex spatial distribution in the country.

Conclusions

The present study shows that sporotrichosis in Venezuela is caused at least by two species: S. schenckii s. str. and S. globosa, the latter representing the highest incidence reported in the Americas (30%). Our results firmly demonstrate that phenotypic key features for species distinction should be used with extreme caution and molecular methods should be always employed. Sporotrichosis clinical forms could be related to genetic variation of S. schenckii spp. complex and their virulence profiles. We strongly support application of ITS as a fungal barcoding gene for species distinction of all clinical and environmental clades of Sporothrix, and accidental infective agents within the Ophiostomales.

Declarations

Acknowledgments

We are grateful to Drs. L. Bezerra (Universidade do Estado do Rio de Janeiro, Brazil), and J. Cano (Tarragona, Spain) for providing references isolates. This work was partially supported by Research Project N° 112 to GAN-V from Instituto Venezolano de Investigaciones Cientificas (IVIC, Venezuela) and grant PEI2012000827 to EC from Fondo Nacional de Ciencia y Tecnologia (FONACIT, Venezuela).

Authors’ Affiliations

(1)
Centro de Microbiología y Biología Celular, Instituto Venezolano de Investigaciones Científicas
(2)
Laboratorio de Micología, Instituto de Biomedicina

References

  1. Lopez-Romero E, Reyes-Montes Mdel R, Perez-Torres A, Ruiz-Baca E, Villagomez-Castro JC, Mora-Montes HM, et al. Sporothrix schenckii complex and sporotrichosis, an emerging health problem. Future Microbiol. 2011;6(1):85–102.View ArticlePubMedGoogle Scholar
  2. Rodrigues AM, de Melo Teixeira M, de Hoog GS, Schubach TM, Pereira SA, Fernandes GF, et al. Phylogenetic analysis reveals a high prevalence of Sporothrix brasiliensis in feline sporotrichosis outbreaks. PLoS Negl Trop Dis. 2013;7(6):e2281.View ArticlePubMedPubMed CentralGoogle Scholar
  3. Mackinnon JE, Conti-Diaz IA, Gezuele E, Civila E, da Luz S. Isolation of Sporothrix schenckii from nature and considerations on its pathogenicity and ecology. Sabouraudia. 1969;7(1):38–45.View ArticlePubMedGoogle Scholar
  4. Kauffman CA. Sporotrichosis. Clin Infect Dis. 1999;29(2):231–6.View ArticlePubMedGoogle Scholar
  5. Lopes-Bezerra LM, Schubach A, Costa RO. Sporothrix schenckii and sporotrichosis. An Acad Bras Cienc. 2006;78(2):293–308.View ArticlePubMedGoogle Scholar
  6. Romeo O, Scordino F, Criseo G. New insight into molecular phylogeny and epidemiology of Sporothrix schenckii species complex based on calmodulin-encoding gene analysis of Italian isolates. Mycopathologia. 2011;172(3):179–86.View ArticlePubMedGoogle Scholar
  7. Dixon DM, Salkin IF, Duncan RA, Hurd NJ, Haines JH, Kemna ME, et al. Isolation and characterization of Sporothrix schenckii from clinical and environmental sources associated with the largest U.S. epidemic of sporotrichosis. J Clin Microbiol. 1991;29(6):1106–13.PubMedPubMed CentralGoogle Scholar
  8. Marimon R, Gené J, Cano J, Trilles L, Dos Santos Lazera M, Guarro J. Molecular phylogeny of Sporothrix schenckii. J Clin Microbiol. 2006;44(9):3251–6.View ArticlePubMedPubMed CentralGoogle Scholar
  9. Barros MB, de Almeida Paes R, Schubach AO. Sporothrix schenckii and Sporotrichosis. Clin Microbiol Rev. 2011;24(4):633–54.View ArticlePubMedPubMed CentralGoogle Scholar
  10. Verma S, Verma GK, Singh G, Kanga A, Shanker V, Singh D, et al. Sporotrichosis in sub-himalayan India. PLoS Negl Trop Dis. 2012;6(6):e1673.View ArticlePubMedPubMed CentralGoogle Scholar
  11. Arrillaga-Moncrieff I, Capilla J, Mayayo E, Marimon R, Mariné M, Gené J, et al. Different virulence levels of the species of Sporothrix in a murine model. Clin Microbiol Infect. 2009;15(7):651–5.View ArticlePubMedGoogle Scholar
  12. Fernández-Silva F, Capilla J, Mayayo E, Guarro J. Virulence of Sporothrix luriei in a murine model of disseminated infection. Mycopathologia. 2012;173(4):245–9.View ArticlePubMedGoogle Scholar
  13. Fernandes GF, dos Santos PO, Rodrigues AM, Sasaki AA, Burger E, de Camargo ZP. Characterization of virulence profile, protein secretion and immunogenicity of different Sporothrix schenckii sensu stricto isolates compared with S. globosa and S. brasiliensis species. Virulence. 2013;4(3):241–9.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Marimon R, Serena C, Gené J, Cano J, Guarro J. In vitro antifungal susceptibilities of five species of Sporothrix. Antimicrob Agents Chemother. 2008;52(2):732–4.View ArticlePubMedGoogle Scholar
  15. Rodrigues AM, de Hoog GS, de Cassia Pires D, Brihante RS, Sidrim JJ, Gadelha MF, et al. Genetic diversity and antifungal susceptibility profiles in causative agents of sporotrichosis. BMC Infect Dis. 2014;14:219.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Rodrigues AM, de Hoog S, de Camargo ZP. Emergence of pathogenicity in the Sporothrix schenckii complex. Med Mycol. 2013;51(4):405–12.View ArticlePubMedGoogle Scholar
  17. Zhou X, Rodrigues A, Feng P, de Hoog GS. Global ITS diversity in the Sporothrix schenckii complex. Fungal Divers. 2014;66(1):153–65.Google Scholar
  18. Marimon R, Cano J, Gené J, Sutton DA, Kawasaki M, Guarro J. Sporothrix brasiliensis, S. globosa, and S. mexicana, three new Sporothrix species of clinical interest. J Clin Microbiol. 2007;45(10):3198–206.View ArticlePubMedPubMed CentralGoogle Scholar
  19. Marimon R, Gené J, Cano J, Guarro J. Sporothrix luriei: a rare fungus from clinical origin. Med Mycol. 2008;46(6):621–5.View ArticlePubMedGoogle Scholar
  20. Dias NM, Oliveira MM, Santos C, Zancope-Oliveira RM, Lima N. Sporotrichosis caused by Sporothrix Mexicana, Portugal. Emerg Infect Dis. 2011;17(10):1975–6.View ArticlePubMedPubMed CentralGoogle Scholar
  21. Morrison AS, Lockhart SR, Bromley JG, Kim JY, Burd EM. An environmental Sporothrix as a cause of corneal ulcer. Medical Mycol Case Reports. 2013;2:88–90.View ArticleGoogle Scholar
  22. Bommer M, Hutter ML, Stilgenbauer S, de Hoog GS, de Beer ZW, Wellinghausen N. Fatal Ophiostoma piceae infection in a patient with acute lymphoblastic leukaemia. J Med Microbiol. 2009;58(Pt 3):381–5.View ArticlePubMedGoogle Scholar
  23. Mariat F. Adaptation of Ceratocystis to a parasitic life in animals–aquisition of a pathogenicity comparable to Sporothrix schenckii. Sabouraudia. 1971;9(3):191–205.View ArticlePubMedGoogle Scholar
  24. Mata-Essayag S, Delgado A, Colella MT, Landaeta-Nezer ME, Rosello A, Perez de Salazar C, et al. Epidemiology of sporotrichosis in Venezuela. Int J Dermatol. 2013;52(8):974–80.View ArticlePubMedGoogle Scholar
  25. Calcagno AM, Nino-Vega G, San-Blas F, San-Blas G. Geographic discrimination of Paracoccidioides brasiliensis strains by randomly amplified polymorphic DNA analysis. J Clin Microbiol. 1998;36(6):1733–6.PubMedPubMed CentralGoogle Scholar
  26. O’Donnell K. Molecular phylogeny of the Nectria haematococca-Fusarium solani species complex. Mycologia. 2000;92:919–38.View ArticleGoogle Scholar
  27. Zhang J, Madden TL. PowerBLAST: a new network BLAST application for interactive or automated sequence analysis and annotation. Genome Res. 1997;7(6):649–56.PubMedPubMed CentralGoogle Scholar
  28. Duong TA, de Beer ZW, Wingfield BD, Wingfield MJ. Phylogeny and taxonomy of species in the Grosmannia serpens complex. Mycologia. 2012;104(3):715–32.View ArticlePubMedGoogle Scholar
  29. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
  30. Nei M. Molecular evolutionary genetics. New York: Columbia University Press; 1987.Google Scholar
  31. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.View ArticlePubMedGoogle Scholar
  32. Oliveira MM, Almeida-Paes R, Muniz MM, Gutierrez-Galhardo MC, Zancope-Oliveira RM. Phenotypic and molecular identification of Sporothrix isolates from an epidemic area of sporotrichosis in Brazil. Mycopathologia. 2011;172(4):257–67.View ArticlePubMedGoogle Scholar
  33. Rodrigues AM, de Hoog G, Zhang Y, de Camargo ZP. Emerging sporotrichosis is driven by clonal and recombinant Sporothrix species. Emerg Microbes Infect. 2014;3:e32.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Madrid H, Cano J, Gené J, Bonifaz A, Toriello C, Guarro J. Sporothrix globosa, a pathogenic fungus with widespread geographical distribution. Rev Iberoam Micol. 2009;26(3):218–22.View ArticlePubMedGoogle Scholar
  35. Sasaki AA, Fernandes GF, Rodrigues AM, Lima FM, Marini MM, Dos SFL, et al. Chromosomal polymorphism in the Sporothrix schenckii complex. PLoS One. 2014;9(1):e86819.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Yu X, Wan Z, Zhang Z, Li F, Li R, Liu X. Phenotypic and molecular identification of Sporothrix isolates of clinical origin in Northeast China. Mycopathologia. 2013;176(1–2):67–74.View ArticlePubMedPubMed CentralGoogle Scholar
  37. Ottonelli Stopiglia CD, Magagnin CM, Castrillon MR, Mendes SD, Heidrich D, Valente P, et al. Antifungal susceptibilities and identification of species of the Sporothrix schenckii complex isolated in Brazil. Med Mycol. 2014;52(1):56–64.PubMedGoogle Scholar

Copyright

© Camacho et al.; licensee BioMed Central. 2015

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement