Stachybotrys chartarum is frequently isolated from damp building materials or improperly stored animal forage. Human and animal exposure to the secondary metabolites of this mold is linked to severe health effects. The mutually exclusive production of either satratoxins or atranones defines the chemotypes A and S. Based upon the genes (satratoxin cluster, SC1-3, sat or atranone cluster, AC1, atr) that are supposed to be essential for satratoxin and atranone production, S. chartarum can furthermore be divided into three genotypes: the S-type possessing all sat- but no atr-genes, the A-type lacking the sat- but harboring all atr-genes, and the H-type having only certain sat- and all atr-genes.
1. Introduction
The family of
Stachybotriaceae has been described as genetically hyperdiverse and heterogeneous concerning the different secondary metabolite classes they produce (e.g., atranones, phenylspirodrimanes, and trichothecenes)
[1,2][1][2]. According to the MycoBank database
[3[3][4],
4], the genus
Stachybotrys currently consists of 126 species. Up to now, only two of them (
S. chartarum and
S. chlorohalonata) are frequently isolated and have been associated with illnesses affecting humans and animals
[1,5,6][1][5][6].
S. chartarum occurs ubiquitously and is frequently found on dead plants (e.g., straw) and other cellulosic (e.g., culinary herbs) and building materials (e.g., gypsum and wallpaper)
[2,7,8,9][2][7][8][9].
S. chartarum is either subdivided into two distinct chemotypes (A and S) based on their production of either atranones (chemotype A) or macrocyclic trichothecenes (MT; chemotype S)
[10] or into the genotypes A, S, and H according to the presence or absence of genes that are presumed to encode the relevant enzymes for the biosynthesis of these mycotoxins (
atr1-14 and
sat1-21)
[11]. The highly cytotoxic chemotype S was found to produce MT, whereas cultures of the low-cytotoxic chemotype A produce atranones
[10,12,13,14][10][12][13][14]. MT represent the most cell-toxic trichothecenes currently known and include roridin, verrucarin, and satratoxins
[15].
By binding irreversibly to the 60S ribosomal subunit, satratoxins inhibit protein biosynthesis and can induce apoptosis in neuronal cell lines
[16,17,18,19][16][17][18][19]. The genotype S has been implicated in several types of diseases
[20,21,22,23,24,25,26][20][21][22][23][24][25][26]: in animals, stachybotryotoxicosis can occur after oral uptake of mycotoxins, especially in horses and less frequently in cattle and sheep
[27,28,29][27][28][29]. Humans, especially infants, are primarily at risk after exposure to toxins in water-damaged buildings
[20,30][20][30]. Exposure to MT may cause pulmonary hemorrhage in infants or symptoms related to the sick building syndrome complex
[20,31,32][20][31][32]. Data on the biological activity of atranones are still scarce, but the available information suggest that they are less cytotoxic than MT
[5]. Rand et al.
[33] showed that relatively high concentrations of atranone A and C (2.0–20 μg/animal) can cause an inflammatory response in mice, but clinical case reports for humans have not been published yet. Semeiks et al.
[34] described the satratoxin cluster 1–3 (SC 1–3) harboring the
sat genes 1 to 21 to be essential for satratoxin production and the atranone core cluster (CAC or AC1) that contains the
atr genes (
atr1–14) and may suffice for the production of atranones. The AC2 was described as a second atranone-specific gene cluster, but Semeiks et al.
[34] already noted that three of the six genes are also conserved in genotype S strains. The function of this gene cluster is still unclear.
Based upon recent genetic information, three genotypes of the fungus have been introduced: genotype A produces atranones but no satratoxins and possesses
atr but no
sat genes; genotype H (hybrid) produces no satratoxins and has a truncated cluster of
sat genes and all
atr genes; and genotype S produces satratoxins and has the complete set of
sat genes but no
atr genes
[5]. Remarkably, no
S. chartarum strains have been isolated to date that lack the
atr and the
sat genes. The three genotypes have common but also distinct features, and their evolutionary relationship is unknown. Currently, genomes of two satratoxin-producing strains (IBT 40293 and 7711) but only one genotype A strain (IBT 40288) are available in the NCBI (National Center for Biotechnology Information) database (ASM102136v1 - GCA_001021365.1 - strain 51-11, S40293v1 - GCA_000732565.1 - strain IBT 40293, S7711v1 - GCA_000730325.1 - strain IBT 7711). The fourth
S. chartarum genome that is available belongs to strain 51-11, for which no conclusive information exists so far concerning its chemotype. Up to now, no genome is available for a genotype H strain. Any analysis of the
S. chartarum genomes is hampered by the fact that they are not fully assembled. Only sets of genomic scaffolds are deposited in the database and gene annotations are available only for three of the four genomes (IBT 7711, IBT 40293, and IBT 40288).
2. Current Insights
Mycotoxins are produced by distinct and elaborated biosynthetic pathways that have often been elucidated in only one species
[40][35]. The class of trichothecenes is particular variable in its chemical structure and has therefore been divided, based on functional groups, in the types A to D
[18]. MTs are also known as type D trichothecenes; they are complex molecules with particularly high toxicity compared to the simpler trichothecenes of the types A to C
[18].
Certain MTs are also produced by different fungi other than
S. chartarum, but within the species, they are exclusively produced by genotype S strains. Atranones belong to a different mycotoxin class and are generated by genotype A strains of
S. chartarum as well as other
Stachybotrys species, e.g.,
S. chlorohalonata. Currently,
S. chartarum is not included in the MIBiG (Minimum Information about Biosynthetic Gene cluster)
[41][36] since a functional validation for the SC and AC gene clusters is still pending. MTs are the putative end products of the trichothecene biosynthesis pathway that is still elusive, but several possible mechanisms have been proposed in the literature
[18,34,42][18][34][37].
Biosynthetic gene clusters are highly variable, rapidly evolving, and often found on accessory chromosomes or at the ends of chromosomes
[43,44,45][38][39][40]. The available genomes of
S. chartarum are not fully annotated and exist as sets of genomic scaffolds. Hence,
wthe researche
rs know the precise position of an AC or SC within the respective scaffold but not on the corresponding chromosome. Most of these clusters are apparently not residing at the very end of a chromosome. This is less clear for AC1, for which only a short downstream sequence is present in the database.
The genome of
S. chartarum 51-11 comprises the most continuous sequences and thereby provides a framework to search for connections between individual scaffolds in the other
S. chartarum genomes. Horizontal gene transfer may be involved in the evolution of the different
S. chartarum genotypes, and this can result in a divergent G+C content
[46][41].
WThe researche
rs therefore analyzed the G+C content of the gene clusters and the flanking regions of the respective scaffolds. All gene clusters have a similar G+C content as the rest of the genome, suggesting that they originate from
Stachybotrys or a fungus with similar G+C content.
The genotype S is defined by the presence of the SC1-3 and the absence of the AC1
[34], criteria which are met by IBT 40293 and IBT 7711. No annotation is available for the genome sequence of strain 51-11. In this
res
tudy, weearch, the researchers have identified an SC1 and SC2 but no AC1 or SC3; hence, 51-11 does not fully match the criteria of a genotype S strain, and
wethe researchers therefore refer to it as an atypical S-type strain. The analysis of the flanking regions of the SC1-3 clusters indicates that their chromosomal positions are identical or at least very similar in the three S-type strains, which corroborates preliminary findings
[11]. From these data, it appears unlikely that the SCs are very mobile genetic elements in the
S. chartarum genome.
The situation is different for the AC. The AC1 cluster is present in the A-type strain IBT 40288 and H-type but not in S-type strains
[11,34][11][34]. AC2 is not an A-type-specific region since it is present in all A- and S-type genomes analyzed in this
res
tudyearch. According to
ourthe data, it is unlikely that the AC2-encoded proteins are directly involved in the production of atranones. No information is currently available on the presence and position of AC2 in H-type strains. In the A- and S-type genomes, AC2 resides in different positions, indicating that these clusters were acquired in several distinct events or that the respective regions underwent substantial genomic rearrangements.
A central question arising from the different genotypes of
S. chartarum is their evolutionary relationship.
WThe researche
rs have obtained several hints that shed light on the relations of the different genotypes and encouraged
uthe researchers to develop a model that is consistent with
ourthe findings.
For unknown reasons, all
S. chartarum strains
wthe researche
rs know of possess either an AC1 and/or specific satratoxin gene clusters. This is true for the A- and S-type but also for H-type strains
[11]. In a study comprising 105
S. chartarum strains, none possessed or lacked both the AC1 and the SC2
[47][42]. This strongly suggests that AC1 and SC2 are mutually exclusive. It is remarkable in this context that
the weresearchers have no evidence that
S. chartarum strains exist in which an intact SC2 coexists with an incomplete AC1 or vice versa. This finding is striking and suggests that acquisition of either of these gene clusters always entails a complete excision of the other. Although the precise pathway for the biosynthesis of MTs or atranones is unknown,
wthe researche
rs know that farnesylpyrophosphate is the starting compound for both processes
[14,18][14][18]. It appears that each strain has to choose whether it converts farnesylpyrophosphate into the direction of atranones or satratoxins. This is also supported by the genome of strain 51-11. Unfortunately, this strain is currently unavailable for culturing and MT analysis. The published data on the toxicity of 51-11 provide only hints but no concrete information on its MT production
[39][43]. Thus, it is unclear whether strain 51-11 is a typical chemotype S strain. It would be interesting to know whether this strain produces and releases satratoxins even though the supposed MFS transporter (encoded in the SC3) is missing. Strain 51-11 differs from the other S-type strains since it lacks not only the SC3 cluster but also most of the
sat1 gene of the SC1. This apparent loss of genetic information indicates that 55-11 developed from a typical S-type strain (
Figure 71).
Figure 71. Current model of the evolution of the S. chartarum genotypes S, A, and H. The different gene clusters are depicted as boxes. The size and the position of these boxes do not represent the correct size and position of the respective cluster in the genome. Satratoxin gene clusters are shown in red and atranone gene clusters in blue.
If AC1 and SC2 are mutually exclusive, it is conceivable that an ancestor harbored one of these clusters, but the question remains: which one was first? The genotype H may provide clues for a better understanding of the evolutionary processes that shaped the current genotypes since it combines elements present in A- or S-type strains. According to the hypothesis of Semeiks et al.
[34], the SC2 developed from a duplication of the CTC, which suggests that the S-type emerged from an A- or H-type ancestor. The fact that
wthe
researchers found truncated
atr1 sequences in all S-type strains strengthens this hypothesis. It is, in principle, possible that the duplication of the CTC took place in an ancestral strain that contained only the CTC cluster, but such a strain has not been identified yet. Following the duplication of the CTC, a functional divergence may have resulted in the emergence of the SC2 and the parallel loss of the AC1. In this context, it is remarkable that BlastP searches identified a gene cluster in
Monosporascus species, which encodes proteins that are homologous to specific SC2-encoded proteins but also trichodiensynthases that are homologous to proteins, which are encoded in the CTC of
S. chartarum (Table S2). Since
Monosporascus is an important spoilage agent in melons
[48][44] and was recently presumed to produce MT (roridin E)
[49][45], this finding may also be relevant with respect to potential adverse health effects to humans.
OuThe researcher
s' analyses also showed that the region located immediately downstream of the SC1 is particularly interesting. It is striking that the SC1-down-1 gene is well conserved in all S-type strains but truncated in the A-type strain and that the lost part of the gene would be most proximal to the SC1.
WeThe researchers assume that the SC1-down-1 gene was destroyed by an imprecise excision of the SC1, which implies that the A-type developed from a strain harboring the SC1, hence an H- or S-type strain.
The A-type strain IBT 40288 possesses a complete AC1 cluster, and
outhe researcher
s' previously published PCR data suggest the same for all H-type strains analyzed so far
[11]. Further data strongly indicate that H-type strains are genetically heterogeneous and often lack one or more of the SC1 and/or SC3 genes
[11,47][11][42]. This loss of genetic information can be construed as an erosive process, in which H-type strains derive from a strain harboring intact SC1 and SC3 and developed towards strains that lack these genetic elements. To integrate this into
the ourresearchers' model,
wethe researchers postulate the existence of an ancestral H-type variant with intact and functional SC1 and SC3 clusters that is represented by strain CBS 324.65.
Based on the currently available genomic information (summarized in
Figure 12,
Figure 53 and
Figure 64),
wthe researche
rs propose the model that is depicted in
Figure 71; it suggests that the H-type is the ancestor rather than the successor of the A-type. The H-type comprises a heterogeneous group of strains that have in common the presence of the SC1 and SC3 and the absence of the SC2 but differ with respect to the genetic erosion of the SC1 and SC3.
OurThe model takes four findings into account: (1) the presence of truncated
atr1 sequences in the S-type strains (
Figure S3), (2) the presence of incomplete SC3 clusters in H-type strains, (3) the truncated SC1-down-1 gene in the A-type strain, and (4) the loss of the SC3 and the
sat1 gene in the atypical S-type strain 51-11.
OurThe model implies that H-type strains represent a crucial element for a better understanding of the evolution of the
S. chartarum genotypes.
Figure 2. Satratoxin clusters 1–3 (SC1-3, (A–C)) and the core trichothecene cluster (CTC) of S. chartarum strains 51-11, IBT 40293, IBT 7711 (genotype S strains), and IBT 40288 (genotype A strain). Arrows depict genes; the different gene clusters are boxed. A dotted line indicates that a sequence is not existent in the respective genome. The genes in square brackets (Panel C) are annotated to a different scaffold than the other genes of these strains, but the PCR and sequencing results demonstrate for IBT 40293 that these regions are part of a continuous sequence. (Panel A) also shows a region of three homologous genes of the Purpureocillium lilacinum isolate PLBJ-1 for comparison. The brown arrows in (Panel C) indicate a set of four genes that are present in all four strains.
Figure 3. The atranone cluster 1 (AC1) region annotated for four S. chartarum strains. Arrows indicate genes. The boxed genes represent the AC1. Genes in square brackets are not part of the scaffold harboring AC1 and/or its upstream region but reside instead in other scaffolds. (IBT 7711: upstream = scaffolds 1183 and (1150)), IBT 40293: upstream = scaffold 1452; downstream = (scaffold 1474)), IBT 40288: scaffold 1 and 51-11: upstream = scaffold 47; downstream are five genes of scaffold (127) and two genes of scaffold (11)). An identical coloration indicates affiliation to a common scaffold in the respective strain.
Figure 4. The atranone cluster 2 (AC2) region annotated for S. chartarum strains IBT 7711, IBT 40293, IBT 40288, and 51-11. Arrows indicate genes. The boxed genes represent the AC2 and encode the hypothetical proteins: S40288_09035, S40288_11459, S40288_09036, S40288_09037, S40288_09038, and S40288_11460 (from left to right). Genes in square brackets are not part of the scaffold harboring the AC2 but reside in other scaffolds. An identical coloration indicates affiliation to a common scaffold in the respective strain that is distinct from that harboring the AC2.
Another essential question that arises from the current
res
tudyearch is how the different gene clusters collaborate to facilitate the production of mycotoxins, e.g., satratoxins and atranones. Usually, gene clusters contain so-called backbone biosynthesis genes for the production of individual mycotoxins. Genes that are engaged in the same secondary pathway commonly form gene clusters that are stable and rarely divided
[50,51][46][47]. In
S. chartarum, three clusters have been implicated in the biosynthesis of satratoxins despite the fact that they are spread in the genome. This situation is similar for the gene clusters that enable the synthesis of the T-2 toxin in
Fusarium sporotrichioides [52[48][49],
53], which is chemically closely related to the satratoxins. A physical separation of gene clusters suggests that they have different functions but does not exclude synergistic interactions. All three SC harbor one putative transcription factor (SC1: Sat9, SC2: Sat15, and SC3: Sat20). This pattern is often found for secondary metabolite clusters and allows a separate regulation for each set of genes. Apart from the biosynthesis, mycotoxins have to be released into the environment. Major Facilitator Superfamily (MFS) transporters have been implicated in this process
[54][50]. It is remarkable that the only MFS protein encoded by the
sat genes is Sat21, which is encoded in the SC3. Semeiks et al.
[34] therefore assumed that the SC3 plays an important role in the transport of trichothecenes, a notion that is further supported by the fact that the SC3 and the CTC reside next to each other. The fact that the atranone-producing strains either lack a SC3 cluster (A-type) or often lack the
sat21 gene (certain H-type strains)
[47][42], indicating that atranones are most likely exported via a mechanism that does not involve Sat21. How the proteins that are encoded by the SC clusters interact and cooperate is still an open question. Whether the mycotoxin biosynthesis takes place in the cytoplasm or whether certain steps are concentrated in certain intracellular compartments or organelles is another. Sat10 is a further noticeable SC-encoded protein. It comprises 16 ankyrin repeat domains and is therefore supposed to be a scaffolding protein that interacts with other Sat proteins to coordinate their activities. The
sat10 gene is the last gene of the SC1 cluster. Using BlastN,
wthe researche
rs found that the two
S. chartarum genes downstream of the SC1 are well conserved in the fungus
Purpureocillium lilacinum. Strikingly, a neighboring gene shows high homology to
sat10 and also encodes a protein with 16 ankyrin repeat proteins. The fact that these three genes are well conserved both with respect to their sequence and their orientation strongly suggests that they form a genetic element that was transferred from
S. chartarum to
P. lilacinum or vice versa. Such a horizontal gene transfer has been described for other fungal genes, e.g., for those involved in the sterigmatocystin synthesis in
Aspergillus species
[55,56][51][52]. Similar mechanism may have also enabled the transfer of MT-related genes to or from
Myrothecium species. The fact that the plant
Baccharis mesapotamica is able to produce roridins and verrucarins even suggests the possibility of an inter-kingdom transfer of genetic information
[57,58][53][54].
Three distinct mechanisms can drive the evolution of functionally diverse gene clusters: de novo assembly, functional divergence, and horizontal transfer
[43][38]. De novo assembly is the most time-consuming and therefore difficult step; it is required to generate an initial gene cluster that can then develop further.
OurThe BlastP searches identified proteins with substantial homology to Sat or ATR proteins in a variety of fungal species, and in most cases, the domains present in the Sat or Atr proteins are also well conserved. Remarkably, only few species harbored several of these genes in their genome.
S. elegans is one example, and its genome contains three genes encoding proteins that are homologous to Sat1, Sat2, and Sat3, but strikingly, these genes are not clustered. Hence, a pool of fungal genes exist that is scattered over many species and likely provided the raw material for the assembly of gene clusters, such as those comprising the
sat and
atr genes. It is likely that the assembly of such clusters occurs gradually starting with small clusters that are successively enlarged by additional genes. The identification of the gene cluster in
Monosporascus spp. and
Xylaria grammica indicate that gene assemblies that are related to SC2 and SC3 exist in other fungi. A closer analysis and comparison of these clusters may provide further important insights. Functional divergence could explain how the SC2 emerged from a duplication of the CTC. Alternatively, the CTC and SC2 may derive from a common ancestor, and the
Monosporascus gene cluster in fact appears to be a chimera of CTC and SC2. The presence of similar gene clusters in
S. chartarum and other fungi strongly suggest that horizontal gene transfer played an important role in the evolution of the
sat and
atr gene clusters.