Wallemia Johan-Olsen (Wallemiales, Wallemiomycetes) is a genus of cosmopolitan xerophilic fungi that can be found in a wide variety of environments that are characterized by low a_w (Samson et al., 2004; Zalar et al., 2005b). Its phylogenetic position was unclear until recently, and previously it has been placed in various positions in the Basidiomycota phylogenetic tree, from the root of basidiomycetes (Zalar et al., 2005b), to incertae sedis (Hibbett et al., 2007), to being a sister group of the Agaricomycotina and Ustilaginomycotina (Matheny et al., 2006). Genome sequencing has shown that it is indeed a sister group of the Agaricomycotina (Padamsee et al., 2012; Zajc et al., 2013). The Wallemiomycetes split from the Agaricomycotina ancestors an estimated 250 million years ago (Zajc et al., 2013). Initially the genus contained only one species, but it was later segregated into three species based on differences in conidial size, xerotolerance, and sequence data: W. ichthyophaga, Wallemia sebi and Wallemia muriae (Zalar et al., 2005b). To date, only around 20 strains of W. ichthyophaga have been isolated from hypersaline waters of solar salterns, bitterns (i.e., magnesium-rich residual solutions in salt production from sea water) and salted meat (Zalar et al., 2005b). In addition to phylogenetic differences, W. ichthyophaga is also distinguished from the other two representatives of this genus by its characteristic morphology and halophilic physiology (Zalar et al., 2005b; Kralj Kunčič et al., 2010).

Although xerotolerance is rare in the Basidiomycota, all three Wallemia spp. are among the most xerophilic fungal taxa known to date (Zalar et al., 2005b). However, while W. sebi and W. muriae strongly prefer high concentrations of non-ionic solutes over those of NaCl (Kralj Kunčič et al., 2013), the opposite is true for W. ichthyophaga (Zalar et al., 2005b). For growth, W. ichthyophaga requires at least 1.5 M NaCl, or some other osmolyte at an equivalent a_w. Such a narrow ecological amplitude of salt concentrations is common for specialized archaeal halophiles, but it is exceptional in the fungal kingdom. Hence, W. ichthyophaga is a rare fungal example of an obligate extremophilic specialist (Gostinčar et al., 2010), and it is considered to be the most halophilic fungus known to date. Although it even thrives in saturated NaCl solution, its in vitro growth optimum is between 2.6 M and 3.5 M NaCl, which is the highest described among fungi (Zajc et al., 2014). It also tolerates high concentrations of salts other than NaCl; e.g., MgCl₂ (our unpublished data).

In H. werneckii and W. ichthyophaga, several components of the HOG pathway have been identified and characterized (Lenassi and Plemenitas, 2007; Fettich et al., 2011; Konte and Plemenitaš, 2013). Since their sequenced genomes became available, the presence of some novel HOG components has been confirmed through homology searches. Altogether, there are many similarities between the HOG pathways in H. werneckii, W. ichthyophaga, and S. cerevisiae, although there are also some important differences that might explain the different halotolerant/halophilic characters of H. werneckii and W. ichthyophaga.

The presence of homologs of Sho1, Ste20, and Ste11 has been confirmed for the genomes of H. werneckii and W. ichthyophaga, with two copies of each component in H. werneckii (Figure 1; Lenassi et al., 2013; Zajc et al., 2013; our unpublished data). Two isoforms of the H. werneckii putative osmosensor protein HwSho1A and HwSho1B fully complement the function of the homologous S. cerevisiae Sho1 protein and they can activate the HOG pathway under osmotic stress in S. cerevisiae. Structurally, when compared to other fungal Sho1 homologs, they contain a conserved SH3 domain and a divergent Ste11-binding motif (Fettich et al., 2011). On the other hand, the SH3 domain of W. ichthyophaga (Wi)Sho1 is functional when it is attached to the N-terminal part of S. cerevisiae Sho1, although the whole sequence of WiSho1 does not appear to function correctly in S. cerevisiae (our unpublished data). Regardless of the complementation of the HwSho1 protein and the WiSho1 SH3 domain, data from recent preliminary investigations addressing the role of the SHO1 branch pathway in osmo-adaptation do not support the involvement of this branch in the signal transfer downstream to heterologously expressed HwPbs2 and WiPbs2 in S. cerevisiae. WiSte11 also failed to complement ScSte11 in S. cerevisiae ste11Δssk2ssk22Δ cells, further supporting this hypothesis (our unpublished data).

The other branch of the HOG pathway, which was named as SLN1 after the transmembrane Sln1 hybrid histidine kinase, transmits its signals via a Sln1–Ypd1–Ssk1 phosphorelay. Sln1 kinase is inactive under hyperosmolar conditions, where Ssk1 is dephosphorylated and therefore binds to the autoinhibitory region of Ssk2 and Ssk22, which triggers their autophosphorylation (Saito and Posas, 2012). The H. werneckii histidine kinases HwHhk7A and HwHhk7B have been identified and characterized in more detail (Figure 1; Lenassi and Plemenitas, 2007). HwHhk7A and HwHhk7B lack the transmembrane domain, but otherwise they have a typical eukaryotic hybrid histidine-kinase-domain composition. Their transcription in H. werneckii depends on the extracellular salt concentration, and when they are expressed in S. cerevisiae, they increase its osmotolerance (Lenassi and Plemenitas, 2007). Both the H. werneckii and W. ichthyophaga genomes also contain homologs of the group III histidine kinases. This group of cytosolic histidine kinases can act as osmosensors through their HAMP domain repeats (Meena et al., 2010). On the other hand, the membrane-spanning Sln1-like histidine kinase is present in the genome of H. werneckii, again in two copies, but there is no evidence for it in the genome of W. ichthyophaga (Figure 1), which suggests that the group III histidine kinases are probably involved in osmosensing in this halophilic fungus. Of the other proteins involved in the SLN1 branch, two forms of each of the S. cerevisiae homologs have been found in the genome of H. werneckii, HwYpd1A/B, HwSsk1A/B, and two homologs of the MAPKK kinase Ssk2, HwSsk2A/B (Table 2) (Lenassi et al., 2013; our unpublished data). Also W. ichthyophaga has the proteins WiYpd1, WiSsk1, and WiSsk2 (Figure 1, Table 2; Zajc et al., 2013; our unpublished data).

The signals from both branches of the HOG pathway in S. cerevisiae converge at the MAPKK Pbs2 scaffold, which transmits the signals further to Hog1 (Saito and Posas, 2012). Two gene copies of the MAPKK HwPbs2 have been identified in H. werneckii and one in W. ichthyophaga, WiPbs2 (Figure 1; Lenassi et al., 2013; Zajc et al., 2013; our unpublished data). However, preliminary data show that the kinases HwPbs2 and WiPbs2 do not interact with the S. cerevisiae Sho1 protein (our unpublished data). This suggests that the SHO1 branch is not involved in HOG pathway activation in H. werneckii and W. ichthyophaga.

We had previously identified and characterized only one isoform of the final MAPK HwHog1 (Turk and Plemenitaš, 2002; Lenassi et al., 2007); however, the H. werneckii whole genome sequence revealed another copy of the HwHOG1 gene, the characterization of which is currently in progress. We have also identified two Hog1-like kinase paralogs in W. ichthyophaga, although all of the other HOG pathway components are represented by only single gene copies (Figure 1; Konte and Plemenitaš, 2013).

HwHog1A, HwHog1B, WiHog1A, and WiHog1B are all considerably shorter than ScHog1, although they contain the conserved domains and motifs that are characteristic of the MAPKs, such as the ATP-binding region, Asp in the active site, a TGY phosphorylation motif, a common docking domain, and a Pbs2-binding domain. While HwHog1A, HwHog1B, and WiHog1B are fully functional kinases in the S. cerevisiae hog1Δ background (Lenassi et al., 2007; Konte and Plemenitaš, 2013; our unpublished data), WiHog1A can only partly restore the osmotolerance of the hog1Δ strain. Lower phosphorylation levels, lower GPD1 induction, and greater cross-talk with the mating pathway indicate that WiHog1A cannot interact optimally with the protein partners in S. cerevisiae. WiHog1B, on the other hand, is a fully functional kinase in the S. cerevisiae hog1Δ background. Moreover, WiHog1B even improves the salt tolerance of S. cerevisiae (Konte and Plemenitaš, 2013). We have also demonstrated that in contrast to S. cerevisiae, where the levels of HOG1 mRNA remain unchanged when the cells are exposed to osmotic shock (Brewster et al., 1993), the transcript levels of HOG1-like kinases in H. werneckii and W. ichthyophaga are salt-dependent (Lenassi et al., 2007; Konte and Plemenitaš, 2013).

Upon hyperosmotic shock in S. cerevisiae, Hog1 is rapidly phosphorylated and it translocates into the nucleus. After the cell adapts to the higher osmolarity, Hog1 is dephosphorylated by phosphatases in a negative-feedback manner (Hohmann et al., 2007). Phosphorylation patterns in the extremely halotolerant H. werneckii and the obligate halophile W. ichthyophaga appear to be more complex. In H. werneckii, we have observed a phosphorylation mechanism that is similar to that of S. cerevisiae, although the HwHog1 kinase is noticeably phosphorylated only when the H. werneckii cells were exposed to ≥3 M NaCl (Turk and Plemenitaš, 2002). While in S. cerevisiae constitutive Hog1 phosphorylation is lethal (Maeda et al., 1994), in W. ichthyophaga this is not the case. Even more interestingly, W. ichthyophaga has a completely “opposite” phosphorylation pattern to that of S. cerevisiae: WiHog1 kinase is dephosphorylated after hypo-osmotic or hyperosmotic shock in W. ichthyophaga, and it is constitutively phosphorylated under optimal osmotic conditions (3.4 M NaCl). These data indicate an important role for the phosphatases in the regulation of the HOG pathway in W. ichthyophaga (Konte and Plemenitaš, 2013). This model has already been reported for Cryptococcus neoformans, where some serotypes show inverted W. ichthyophaga-like phosphorylation patterns (Bahn et al., 2007).

When activated, HwHog1 is translocated into the nucleus, where it associates with the chromatin of osmoresponsive genes and induces or represses their expression (Vaupotič and Plemenitaš, 2007). A transcriptional response to hyperosmolar stress of 95 differentially expressed genes has been reported for the comparison of moderately (3 M NaCl) and extremely (4.5 M) osmolar environments. Data from the ChIP method show that 36 of these genes physically interact with HwHog1 in long-term adaptation to extreme environments (Vaupotič and Plemenitaš, 2007). In 17 out of these 36 genes, simultaneous co-localization of RNA polymerase II was seen. More than half of differentially expressed genes are related to general metabolism and energy production, and the other osmoresponsive genes are involved in the biogenesis of mitochondria, protein biosynthesis, protein quality control, transport facilitation, the cell cycle, and the cell wall (Vaupotič and Plemenitaš, 2007). Thirteen of these 95 genes could not be classified. Certain osmoresponsive genes controlled by MAPK HwHog1 have been studied in greater detail. Genes that code for the P-type ATPases HwEna1 and HwEna2 are the S. cerevisiae ENA1 homologs, and therefore they are believed to be involved in the maintenance of a low intracellular K⁺/Na⁺ ratio. Their transcription is salt regulated (Gorjan and Plemenitaš, 2006). Two homologs of the S. cerevisiae key enzyme in glycerol biosynthesis, the glycerol-3-phosphate dehydrogenase Gpd1, have been characterized in H. werneckii. These both show similar transcription profiles in response to different salt concentrations (Lenassi et al., 2011). We have also demonstrated that the MAPK WiHog1 can up-regulate the transcription of GPD1 in S. cerevisiae. The regulation of other osmoresponsive genes that are potential targets of WiHog1 remains to be defined.

The genome of H. werneckii was recently sequenced and it has been deposited at DDBJ/EMBL/GenBank under the accession number AIJO00000000 (Lenassi et al., 2013). The genome statistics are summarized in Table 1. The genome of H. werneckii has a size of 51.6 Mb, which is relatively large. In species belonging to the same order as H. werneckii (Capnodiales), the genome sizes are very variable, as they range from 21.88 to 74.12 Mb. The larger genome sizes are mostly due to a substantial amount of repetitive sequences. However, in H. werneckii, despite its large genome size, the proportion of repetitive sequences is only 1.02%. On the other hand, it contains 23,333 predicted genes, which is twice as many as the average number of predicted genes in other related fungi (approx. 11,955 genes) (Ohm et al., 2012). This large number of genes can be attributed to a relatively recent whole genome duplication, which resulted in two nearly identical copies of almost every protein of H. werneckii (Lenassi et al., 2013). This discovery is in line with our previous studies of several individual genes from H. werneckii that were present in two copies (Gorjan and Plemenitaš, 2006; Lenassi and Plemenitas, 2007; Fettich et al., 2011). In most cases, the expression of both of the gene copies is salt dependent, although their expression profiles differ (Lenassi and Plemenitas, 2007). It may well be that as a consequence of this whole genome duplication, H. werneckii can benefit from the potential advantages of large genetic redundancy, even though it is formally in a haploid stage (i.e., it is not a diploid that has resulted from the mating of two strains with opposite mating types, which would regain the haploid stage with meiosis before the next mating event; see below).

The H. werneckii genome sequence has offered the opportunity to gain insight into the genetic information on the mating type(s) and on the mating strategy. To date, no sexual cycle has been described for H. werneckii. Using M. graminicola proteins that contain the alpha1 domain (Mat1-1-1) and the HMG domain (Mat1-1-2), we identified the putative HwMAT1-1-1A and HwMAT1-1-1B genes (Figure 1, Table 2), both of which are translated into 358 amino-acid proteins that contain the alpha1 domain and have an overall amino-acid sequence identity of 87.5% (Lenassi et al., 2013). Importantly, no homologs of the HMG-domain-containing Mat1-1-2 protein were found in H. werneckii, which indicates that this species is heterothallic, and that if it can still undergo sexual reproduction, this requires a strain that codes for the opposite mating type (in the case of the sequenced strain, this would be a strain with a Mat1-1-2 homolog; Lenassi et al., 2013).

Eukaryotic microorganisms have developed numerous plasma-membrane transport systems to maintain their appropriate alkali cation levels, and in particular, to eliminate any surplus of toxic Na⁺ ions. The alkali-cation transport systems in S. cerevisiae and in non-conventional yeast have recently been reviewed (Arino et al., 2010; Ramos et al., 2011).

In S. cerevisiae, the plasma-membrane transporters Trk1 and Trk2 for K⁺ uptake, the Tok1 K⁺ channel, the Pho98 inorganic phosphate (P_i)-Na⁺ symporter, the Ena Na⁺-ATPases, and the Nha1 Na⁺/H⁺ antiporter have all been well characterized (Arino et al., 2010). Together with these, non-specific protein transporters (e.g., Pm3, Qdr2) have been described to be involved in K⁺/Na⁺ fluxes across the plasma membrane (Arino et al., 2010). Trk transporters for K⁺ uptake, Nha antiporters, Ena ATPases, and Tok1 channels have also been identified in non-conventional yeast, together with the Hak K⁺/H⁺ symporters and the rare K⁺/Na⁺-uptake ATPase Acu (Ramos et al., 2011).

Physiological studies have shown that H. werneckii maintains very low intracellular K⁺ and Na⁺ levels (Kogej et al., 2005), even when it grows in the presence of 4.5 M NaCl, which suggested that it can effectively extrude Na⁺ ions and also prevent their influx. Analysis of H. werneckii genome has revealed considerable expansion of families of genes that encode plasma-membrane metal cation transporters, as presented schematically in Figure 1 and summarized in Table 2 (Lenassi et al., 2013). We identified eight homologs of the Trk1 and Trk2 K⁺ channels, with each containing the conserved TrkH domain that is typical for cation transport proteins. In general, they show low homology to the Trk1 protein, but the amino-acid sequence identity increases in the TrkH domain. We also identified four homologs of the Tok1 K⁺ channels, each of which contains two conserved transmembrane helices that are typical of this ion-channel family. Again, the homology to the Tok1 protein is low, but the identity is high in the transmembrane helices. The presence of eight homologs of the Nha1 Na⁺/K⁺, H⁺ antiporters was demonstrated, each of which contains a transmembrane region at the N-terminal, which is conserved through the Na⁺/K⁺, H⁺ exchanger family, and only two of them additionally contain the C-terminal cytoplasmic region. Extensive expansion has also been observed for the Pho89 homologs in H. werneckii, as we identified six homologs of the Pho89 Na⁺, P_i symporter, with each homolog containing at least one PHO4 domain. In contrast with the abundant transporter families mentioned, only four homologs of three S. cerevisiae Ena Na⁺ P-type ATPases have been identified in the H. werneckii genome (Lenassi et al., 2013). Previously, we identified and characterized two Ena-like P-ATPases (Gorjan and Plemenitaš, 2006). Analysis of the Ena Na⁺ P-type ATPases identified in the genome of H. werneckii reveals that each homolog contains all four of the conserved domains found in the S. cerevisiae Ena proteins. Considering their multiplication, it appears that Nha transporters are more important than Ena. On the other hand, based on our previous data that demonstrated that HwENA genes are highly induced at alkaline pH (Gorjan and Plemenitaš, 2006), we speculate that they have complementary functions: Ena ATPases are more important at high pH, where the Nha antiporters cannot function correctly.

As well as the important role of plasma-membrane transport systems in ion homeostasis, in the cytosol, K⁺ homeostasis and Na⁺ detoxification are also connected to cation transport across the organelle membranes (Arino et al., 2010). In S. cerevisiae, endosomal Nhx1 (Nass and Rao, 1999) and Kha1 from the Golgi apparatus (Maresova and Sychrova, 2005) are Na⁺/H⁺ exchangers, similar to Nha1 at the plasma membrane (Prior et al., 1996). The vacuolar Vnx1 (Cagnac et al., 2007) and the mitochondrial Mdm38 and Mrs7 (Nowikovsky et al., 2004; Zotova et al., 2010) have similar Na⁺/K⁺, H⁺ exchanger functions, but different structures.

We found that homologs of Nhx1 and Kha1 are duplicated in the H. werneckii genome, all of which contain the domains that are typical for the Na⁺/H⁺ exchanger family. The same has been observed for the Kha1 homologs. We also identified two homologs of transporters with high homology to the Mrs7 and Mdm38 transporters from S. cerevisiae. Of the intracellular cation transporters, only the homologs of the vacuolar Vnx1 are enriched in H. werneckii in comparison to S. cerevisiae. We identified eight homologs of the Vnx1 Na⁺/K⁺, H⁺ antiporter, but the homology of the HwVnx proteins compared to Vnx1 is low (Lenassi et al., 2013).

The activities of many transporters are closely connected to the proton gradients across the membranes, which are generated by the Pma1 P-type ATPase at the plasma membrane (Serrano et al., 1986; Ambesi et al., 2000) and the V-type ATPase at the vacuolar membrane (Graham et al., 2000). Different P-type ATPases use ATP hydrolysis as a source of energy for the transport of ions through the membrane, and they are structurally similar (Kuhlbrandt, 2004).

The enrichment of the transporters responsible for supplying the energy for the cation transporters in H. werneckii supports the importance of the complex cation transporter system for combating high environmental Na⁺. We identified four homologs of Pma1 in H. werneckii, with each homolog containing three conserved domains that are also found in the S. cerevisiae Pma1 and Pma2 proteins. The importance of all four of the H. werneckii Pma homologs for cation homeostasis is also supported by expression analysis, as the expression profiles of the PMA1 and PMA2 homologs show different levels of response to saline conditions (Lenassi et al., 2013). Comparisons of the expression profiles of the PMA genes in H. werneckii with those described in S. cerevisiae have shown that in S. cerevisiae, PMA1 is not induced by salt stress (Yale and Bohnert, 2001), while in H. werneckii, both PMA1 and PMA2 have salt-regulated transcription.

The yeast vacuolar ATPase does not only have a crucial role in the acidification of the vacuolar lumen, but it is also important for the correct functioning of other organelles (Arino et al., 2010). In the H. werneckii genome, we found homologs of all of the subunits of the S. cerevisiae V-ATPase complex. The H. werneckii vacuolar subunits in general share a lot of similarity with the S. cerevisiae subunits, which is not surprising, as their structures and function have been highly conserved through evolution (Graham et al., 2000). S. cerevisiae vacuolar ATPases are localized at different cellular locations; however, it remains to be determined where they are specifically localized in H. werneckii. In contrast to the transcription of the HwPMAs, which is salt regulated, no such trends have been seen for the expression of the VMA homologs under different salinities.

The genome of W. ichthyophaga has been deposited as a Whole Genome Shotgun project at DDBJ/EMBL/GenBank under the accession number APLC00000000 (Zajc et al., 2013). The genome of W. ichthyophaga is 9.6 Mb in size, and the sequence currently consists of 101 contigs and 82 scaffolds (Table 1). Most basidiomycetous haploid genomes are more than twice this size (and in some cases, larger by 40-fold or more; Gregory et al., 2007). The closely related species W. sebi also has a slightly larger genome (9.8 Mb) (Padamsee et al., 2012). Of the species investigated thus far, only the dandruff- and seborrhoeic-dermatitis-causing Malassezia globosa has a smaller genome (9.0 Mb; Gregory et al., 2007). The compactness of the genome of W. ichthyophaga is reflected in its low level of repetitive sequences (1.67%), and high density of genes (514 genes/Mb scaffold) (Zajc et al., 2013). This is only slightly lower than W. sebi (538 genes/Mb), but more than in M. globosa (476 genes/Mb). This means that the coding DNA sequences in W. ichthyophaga cover almost three quarters of the genome. The GC content in W. ichthyophaga is 45.35%, while in W. sebi this is even lower, at 40.01% (Table 1). The absolute number of predicted proteins in W. ichthyophaga (4884; Zajc et al., 2013) is also unusually small for a basidiomycete (where more than 10,000 proteins are not uncommon), and is in the range observed for Escherichia coli (Lukjancenko et al., 2010). For comparison, the W. sebi genome codes for 5284 proteins, while M. globosa contains 4285 proteins. Interestingly, the reduction in genome size and gene number is not accompanied by a reduction in intron number, such as has been reported for some other fungi with small genomes (Kelkar and Ochman, 2012).

It has not been possible to assign functions for an unproportionally large number of the proteins that are found in W. ichthyophaga but not in W. sebi. With searches through the Pfam database, three quarters of these proteins could not be classified into any of the protein families (Zajc et al., 2013). Among those that could be identified, there were several proteins related to DNA processing and DNA damage.

The sequencing of the transcriptomes of W. ichthyophaga grown in 10% and 30% (w/v) NaCl has revealed that 13.1% of the genes are differentially expressed under these conditions (Zajc et al., 2013). Of these, two thirds are more expressed at lower salinity. Alternative splicing, which has been identified as intron retention, was detected for 15.0% of the genes, and in more than half of the cases (51.6%), alternative splicing was detected only at one of the two tested salinities (Zajc et al., 2013).

The strategy of osmo-adaptation of both the extremely halotolerant H. werneckii and the halophilic W. ichthyophaga is the accumulation of a mixture of polyols that act as compatible solutes (Figure 1). The main osmotically regulated polyol of both H. werneckii and W. ichthyophaga is glycerol, the levels of which are increased with increasing salinity and decreased after hypo-osmotic shock (Kogej et al., 2007; Zajc et al., 2014). In addition to glycerol, we have reported erythritol, arabitol, and mannitol in H. werneckii (Plemenitaš et al., 2008), and smaller amounts of arabitol, and traces of mannitol in W. ichthyophaga (Zajc et al., 2014). The genes for the enzymes known to be involved in compatible solute management are found in the genome of W. ichthyophaga (Table 2). These are present in several copies, with the exception of the glycerol-3-phosphatase Gpp. W. ichthyophaga contains a homolog of GPD1, WiGPD1, the expression of which is salt-induced (Lenassi et al., 2011). A second homolog was also found by searching the genome. When compared with the homologs of H. werneckii, the expression level of WiGPD1 is lower, and the response to hyperosmotic shock is slower (Lenassi et al., 2011). Expression of WiGPD1 in S. cerevisiae boosted the osmotolerance of the gpd1 and gpd1gpd2 mutants. As was reported for homologs from H. werneckii, WiGPD1 lacks the N-terminal peroxisomal targeting (PTS2) sequence (Lenassi et al., 2011), which is important for peroxisome localization of WiGpd1 (Jung et al., 2010). This might mean that WiGpd1 remains in the cytosol, which would be an advantage when living in extremely saline environments, as it is this fraction that is important for the synthesis of the compatible solutes (Lenassi et al., 2011). During hyperosmotic shock, S. cerevisiae counteracts glycerol leakage by its active re-import using the protein Stl1, a glycerol/H⁺ symporter in the plasma membrane (Ferreira et al., 2005). Similarly, the aquaglyceroporin channel Fps1 remains closed (while it opens during hypo-osmotic shock, to facilitate expulsion of excess glycerol) (Luyten et al., 1995). In W. ichthyophaga, four homologs of Stl1 have been found (Figure 1, Table 2), as well as three aquaglyceroporin-related proteins.

In the basidiomycete Agaricus bisporus, the solute D-mannitol is synthesized from fructose via a reduction step that is catalysed by two NADP-dependent mannitol dehydrogenases (Stoop and Mooibroek, 1998). W. ichthyophaga also contains two homologs of D-arabinitol-2-dehydrogenases, which are used in other fungi for the production of arabitol from an intermediate of the pentose phosphate pathway, D-ribulose-5-phosphate.

In line with the strategy of compatible solutes, the intracellular levels of K⁺ and Na⁺ in W. ichthyophaga remain low at constant salinities (not above 30 nmol/mg dry biomass) (Zajc et al., 2013), even when compared to H. werneckii (not above 180 nmol/mg dry biomass) (Kogej et al., 2007). However, when under hyperosmotic shock, the levels of both cations increase significantly in W. ichthyophaga, indicating its poor capability to adjust to changing environments. The ratio between these cations decreases with increasing salinity, due to the rising levels of Na⁺ and the lowering of K⁺. However, the intracellular K⁺/Na⁺ ratio is higher across the whole salinity range in W. ichthyophaga compared to H. werneckii (Figure 1) and some other halotolerant fungi (e.g., Aureobasidium pullulans, Debaryomyces hansenii). In addition, the K⁺/Na⁺ decrease over the salinity range is less steep in W. ichthyophaga compared to H. werneckii. The growth performance of W. ichthyophaga is greatest when the Na⁺ content exceeds that of K⁺ (Zajc et al., 2014). This indicates that these intracellular concentrations of Na⁺ ions are not toxic to the cells.

Data from the genome show that there are only a low number of cation transporters, except for the enriched protein family of P-type ATPases (Figure 1, Table 2). Also, expression of the cation transporters is low and independent of salt, with only three minor exceptions (described below). This is probably associated with the life of W. ichthyophaga at relatively constant (although extremely high) salinities. Nevertheless, in its genome we observed a significant enrichment of the cation-transporting ATPases family. The identified proteins of this family are three H⁺ and two Na⁺ P-type ATPases (all of which are assumed to be located at the plasma membrane), two Ca²⁺ P-type ATPases (vacuolar Pmc1 and Pmr1 from the Golgi apparatus), and a putative transporter of unknown specificity (Zajc et al., 2013). The W. ichthyophaga genome encodes three putative Pma proton pumps, while W. sebi contains only two (Zajc et al., 2013).

In environments with high concentrations of Na⁺ salts, the cell must prevent the intracellular accumulation of the highly toxic Na⁺, without lowering the levels of K⁺. This is achieved by a variety of other secondary active transporters. W. ichthyophaga contains homologs of most known transporters from S. cerevisiae (Arino et al., 2010) and unconventional yeast (Ramos et al., 2011), as either those located on intracellular membranes (Kha1, Mrs7/Mdm37, Nhx1, Pmc1, Pmr1, Vnx1, Vma1) or at the plasma membrane (Ena, Nha1, Pho89, Pma, Trk1).

Judging by the genes that encode alkali metal cation transporters, the extremely halotolerant ascomycete H. werneckii and W. ichthyophaga use different salt-combating strategies. As described above, in H. werneckii, the numbers of most of the plasma-membrane alkali cation transporters are substantially increased. In W. ichthyophaga this is not the case (Figure 1, Table 2). W. ichthyophaga contains only one Trk homolog (inward K⁺ transporter, with eight copies in H. werneckii) and no Tok homologs (outward K⁺ channel, with four copies in H. werneckii). Similarly, W. ichthyophaga has only two Nha homologs (Na⁺/K⁺ proton antiporters) and one Pho89 (Na⁺/P_i symporter), while H. werneckii contains eight and six, respectively (Lenassi et al., 2013; Zajc et al., 2013).

Active import of K⁺ might contribute to ion homeostasis in hypersaline environments, and this would complement the action of passive K⁺ channels. The known active transporters are K⁺-H⁺symporters (Hak symporters) and K⁺(Na⁺)-ATPase (Acu, alkali cation uptake transporters) (Benito et al., 2004; Ramos et al., 2011). While H. werneckii has no homologs of either of these transporter types (Lenassi et al., 2013), W. ichthyophaga contains two possible homologs of the otherwise rare Acu ATPases, of which only one contains a P-type ATPase domain (Zajc et al., 2013).

Several different membrane alkali metal transporters (mainly cation/H⁺ antiporters) are also located on the organelle membranes: the Golgi apparatus (Kha1), mitochondria (Mdm38 or Mrs7), endosomes (Nhx1), and vacuole (Vnx1) (Arino et al., 2010). In W. ichthyophaga, there are single-copy genes of all of these proteins (with the exception of the duplicated Kha1).

W. ichthyophaga can live at extremely high salinity; however, the above-described findings indicate that the salinity remains relatively constant and thus no rapid responses are needed. Under constant conditions, there is no need for a quick release of surplus K⁺, and thus the absence of the Tok outward K⁺ channel might not be detrimental. On the other hand, W. ichthyophaga contains three homologs of aquaglyceroporins (Zajc et al., 2013), which fulfil the need for rapid expulsion of the accumulated compatible solute glycerol (Luyten et al., 1995). Two explanations are possible: (1) W. ichthyophaga might deal with hypo-osmotic shock-related K⁺ expulsion in other ways than in other fungi; and (2) the aquaglyceroporin channels might serve some other functions than expulsion of glycerol during the shock.

In cells grown at 10 and 30% (w/v) NaCl, the large majority of the genes that encode metal-cation transporters are not differentially expressed. This is not as expected, since these transporters are believed to have crucial roles in adaptation to salt. The only exceptions are a homolog of the Pho89 Na⁺/P_i symporter, which shows elevated expression at high salinity, and a putative P-type Na⁺ ATPase and a possible Acu K⁺ importer, with the expression of both of these latter higher at low salinity. This is in stark contrast, for example, with other halotolerant fungi (e.g., D. hansenii, H. werneckii), where even at different constant salinities, differential expression of P-type H⁺ and Na⁺ ATPases has been observed (Almagro et al., 2001; Gorjan and Plemenitaš, 2006; Lenassi et al., 2013). Furthermore, in W. ichthyophaga the expression of transporter coding genes is relatively low: of the total of 4884 genes, neither the Na⁺-exporting P-type ATPases nor the two Na⁺/H⁺ antiporters are among the 2000 most-expressed genes at high salinity. The possible post-transcriptional control of all of these genes remains to be investigated.

Continuous removal and/or compartmentalization of Na⁺ at a constant high salinity is extremely demanding energetically. The low numbers of transporters in W. ichthyophaga might reflect the relatively low adaptive potential of this species to changes in salt concentrations and/or its specialization with other mechanisms that are energetically more efficient. The apparent transcriptional non-responsiveness of transporters to salt can lead to similar conclusions. This would mean that the halophilic strategy of W. ichthyophaga, which is a unique example of a narrowly specialized fungal halophile (Gostinčar et al., 2010), is substantially different from that of H. werneckii, which contains a collection of K⁺ channels and can adapt to a wide salinity range (Zajc et al., 2013).

The analysis of the W. ichthyophaga genome has revealed a significant expansion of seven protein families and contraction of 19. The most interesting of the expanded families are the hydrophobins, which are proteins that potentially have a role in the particular morphological adaptations of W. ichthyophaga (Figure 1; Zajc et al., 2013).

W. ichthyophaga has a characteristic morphology that can be seen in many stress-tolerant species: compact multicellular clumps that are similar to sarcinae (Zalar et al., 2005b). This morphology has been observed in, and it is believed to enhance survival in, high-stress environments (Wollenzien et al., 1995; Palkova and Vachova, 2006; Gostinčar et al., 2011). These cells have an abundant cover of extracellular polysaccharides (Kralj Kunčič et al., 2010) that can serve as protectants during desiccation (Selbmann et al., 2005), and possibly also when exposed to high concentrations of salt (Zajc et al., 2013). Additionally, at high salinity the morphology of W. ichthyophaga cells changes strikingly. The size of the meristematic cell clumps increases four-fold, and the cell walls become three-fold thicker, which also substantially decreases the intracellular volume (Kralj Kunčič et al., 2010).

The hydrophobins are proteins in the cell wall of filamentous fungi. These small (≤20 kDa) and amphipathic molecules (Linder et al., 2005) are involved in a range of processes of cellular growth and development (Wosten, 2001). Possibly as a reflection of the many different roles they have, the hydrophobin genes are often present in multiple different copies. From an estimated 15 hydrophobin genes in the last common ancestor of W. ichthyophaga and W. sebi, these genes are enriched to 26 in W. ichthyophaga, while the number has fallen to 12 in W. sebi. This is the most significant protein family expansion in the genome of W. ichthopyhaga (Zajc et al., 2013).

The hydrophobins that have been identified in W. ichthyophaga and W. sebi contain the characteristic hydrophobin pattern of the conserved spacing of eight cysteine residues (Zajc et al., 2013) that form four disulphide bridges (Hektor and Scholtmeijer, 2005; Linder et al., 2005). These hydrophobins in both W. ichthyophaga and W. sebi contain a high proportion of acidic amino acids compared to homologs from other fungi (Zajc et al., 2013). This is similar to the archaeal halophilic proteins (Madern et al., 2000) and it might represent an adaptation to salt exposure. If acidic amino acids are exposed on a protein surface, they can bind salt and water, and thus help to avoid salt-induced changes in conformation that would lead to loss of activity (Siglioccolo et al., 2011). This might also be the case for the hydrophobins of W. ichthyophaga, as these are among the few of its proteins that are actually exposed to high concentrations of NaCl and are not protected in the intracellular compatible-solute-rich environment. Interestingly, half of the genes encoding hydrophobins are differentially expressed during growth of W. ichthyophaga at different salinities, although their responses are not the same: at high salinity some of these hydrophobins have higher expression, and some have lower expression (Zajc et al., 2013).

The hydrophobins can self-assemble into amphipathic monolayers on hydrophobic–hydrophilic interfaces, which is crucial for their diverse functions. With their help, the cell can attach to hydrophobic surfaces, break through a water-air interface, and avoid water-logging without impeding the exchange of gasses. Hydrophobins can also strengthen the cell wall and make it more rigid, and they also impact on the movement of solutes (Wosten, 2001; Bayry et al., 2012). It is not difficult to imagine that these functions will be beneficial to cells exposed to high salinity (Zajc et al., 2013). Under these conditions, the leakage of compatible solutes and intrusion of toxic salt ions are among the greatest challenges to the cell. Increased strength and rigidity of the cell wall will help the cell to survive the structural stress that can be imposed by changes in the environmental osmolarity.

The hydrophobins might also have a role in the characteristic sarcina-like morphology of W. ichthyophaga, which is similar to Fusarium verticillioides, where the hydrophobins trigger microconidial chain formation (Fuchs et al., 2004). This might additionally be involved in the aggregation of W. ichthyophaga cells into compact cell clusters, as is characteristic of W. ichthyophaga, and which is probably a form of salt-stress response (Gostinčar et al., 2010). This characteristic sarcina-like morphology is considered to be among the main haloadaptations of W. ichthyophaga. Indeed, both the changes in the cell wall and the formation of multicellular structures have previously been suggested to be among the main adaptations of W. ichthyophaga to hypersaline environments (Kralj Kunčič et al., 2010).

Traditional mycological approaches have not lead to any descriptions of mating behaviors for either H. werenckii or the Wallemia spp.. In both cases, no descriptions of fruiting bodies or reports of teleomorphs can be found in the literature. Basidomycetous fungi have a tetrapolar MAT locus. Two additional unlinked loci encode the homeodomain-containing transcription factors and pheromone/pheromone receptors. In some cases, other mating architectures have evolved through the expansion and fusion of these MAT loci (Lee et al., 2010). The genome of W. sebi contains a single mating-type locus and lacks only a few meiosis-specific genes, which suggested that it can undergo sexual reproduction (Padamsee et al., 2012). This is not the case for W. ichthyophaga (Zajc et al., 2013). Searches of the W. ichthyophaga genome for proteins that are similar to the gene products involved in mating in other basidiomycetes has resulted in the identification of only a few proteins, and even these are very dissimilar to proteins from other fungi. No discernible mating locus has been identified, and only three of eight meiosis-specific genes (as listed in Malik et al., 2008) have been found (Zajc et al., 2013).

Asexuality in fungi is not unusual. Asexual reproduction saves energy, as it eliminates the need to produce gametes and attractants, and this might be useful in extreme environments where careful management of energy is of key importance. In habitats that do not undergo major changes over longer time scales, this would also prevent the drowning of specific adaptations of local populations in the larger gene pool of the species (Gostinčar et al., 2010; Sun and Heitman, 2011). Therefore, in the case of these extremophilic fungal species, an asexual lifestyle might have evolutionary advantages (Gostinčar et al., 2010). This appears to be the case for W. ichthyophaga, and this might also be true for H. werneckii. While a putative mating type locus has been found in the genome of H. werneckii, no sexual stage of this species has been described to date. As mating genes can have roles outside of mating (Srikantha et al., 2012), their presence in itself is not confirmation of sexual reproduction in a species. The potential role of sexual reproduction in H. werneckii therefore needs to be further investigated.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.