Genome sequencing of four Aureobasidium pullulans varieties: biotechnological potential, stress tolerance, and description of new species
© Gostinčar et al.; licensee BioMed Central Ltd. 2014
Received: 5 March 2014
Accepted: 20 June 2014
Published: 1 July 2014
Aureobasidium pullulans is a black-yeast-like fungus used for production of the polysaccharide pullulan and the antimycotic aureobasidin A, and as a biocontrol agent in agriculture. It can cause opportunistic human infections, and it inhabits various extreme environments. To promote the understanding of these traits, we performed de-novo genome sequencing of the four varieties of A. pullulans.
The 25.43-29.62 Mb genomes of these four varieties of A. pullulans encode between 10266 and 11866 predicted proteins. Their genomes encode most of the enzyme families involved in degradation of plant material and many sugar transporters, and they have genes possibly associated with degradation of plastic and aromatic compounds. Proteins believed to be involved in the synthesis of pullulan and siderophores, but not of aureobasidin A, are predicted. Putative stress-tolerance genes include several aquaporins and aquaglyceroporins, large numbers of alkali-metal cation transporters, genes for the synthesis of compatible solutes and melanin, all of the components of the high-osmolarity glycerol pathway, and bacteriorhodopsin-like proteins. All of these genomes contain a homothallic mating-type locus.
The differences between these four varieties of A. pullulans are large enough to justify their redefinition as separate species: A. pullulans, A. melanogenum, A. subglaciale and A. namibiae. The redundancy observed in several gene families can be linked to the nutritional versatility of these species and their particular stress tolerance. The availability of the genome sequences of the four Aureobasidium species should improve their biotechnological exploitation and promote our understanding of their stress-tolerance mechanisms, diverse lifestyles, and pathogenic potential.
KeywordsAureobasidium pullulans Dothideomycetes Genome Stress Haloadaptation Halotolerance Polyextremotolerant New species Opportunistic mycosis Pullulan
Aureobasidium pullulans (de Bary) G. Arnaud is a polyextremotolerant black yeast of considerable biotechnological importance, and it has exceptional stress tolerance and increasing medical relevance . A. pullulans is well known for its production of pullulan, a neutral polysaccharide of repeating maltotriose units, which has numerous applications in medicine, pharmacy, the food industry, and other fields [2, 3]. A. pullulans also produces a β-glucan that shows high reactivity to human IgG antibodies  and possible beneficial immunomodulatory effects . A. pullulans has an unusually large spectrum of extracellular enzymatic activities [6, 7]. Several of these are of biotechnological interest, and include: amylases, cellulases, lipases, proteases, xylanases, β-fructofuranosidases, maltosyltransferases, mannanases, and laccases (for review, see ). At least some of these appear to have interesting traits that are different from their homologues in other species [8, 9].
A strain of A. pullulans is used for the production of a cyclic peptide that has specific antifungal activity: aureobasidin A . Due to its strong antagonistic activity against other microorganisms, A. pullulans is used as a biocontrol agent in agriculture . Additionally, a recent study reported that some strains of A. pullulans can produce an antibacterial compound, exophilin A, as well as high yields of liamocins, and heavy oils with previously unknown acylated mannitol structures, which have possible industrial applications as surfactants .
The occurrence of A. pullulans is widespread in tropical, temperate and polar areas. It is frequently found in association with diverse plants, such as in the phyllosphere [13, 14], as an epiphyte or endophyte, on stored barley grain , and in coconut water . A. pullulans has also been found in numerous other habitats, some of which are particularly unusual, such as coastal hypersaline water [17, 18], glacial ice [19, 20], other polar environments [21, 22], polluted water , refrigerated, frozen, salt-preserved and dried foods [24, 25], various indoor habitats (e.g., bathroom surfaces , house dust , dishwashers , tap water ), the surface of human skin , aviation fuel tanks , and the surface of synthetic polymers  and of degrading polyurethane and PVC plastics . A. pullulans has been reported to cause a variety of localised infections in humans, and even systemic infections, although very rarely (for review, see [34, 35]). It has been suggested that this infection potential is at least partially supported by the production of extracellular enzymes , and by the pronounced stress tolerance of A. pullulans.
Aureobasidium pullulans has evolved an exceptional tolerance for a broad range of ecological conditions. It is considered a polyextremotolerant organism [1, 37], and it can survive hypersaline , acidic, basic [38, 39], cold and oligotrophic  conditions. Some of the A. pullulans adaptations to stress (and especially to elevated salt concentrations) are associated with rigorous management of intracellular concentrations of alkali-metal cations , synthesis of compatible solutes  and of mycosporines , and adaptation at the level of the membrane-lipid composition [44, 45]. A. pullulans can also mitigate environmental stress by rapid dimorphic switching from small colourless yeast cells to thick-walled, heavily melanised, meristematic forms .
Various loci have been sequenced in the past to infer the taxonomy and phylogeny of the taxa in Aureobasidium, such as internal transcribed spacer (ITS) rDNA, intergenic spacer 1, translation elongation factor-1α, β-tubulin, and RNA polymerase II [19, 47–49]. Based on a multilocus analysis of a worldwide selection of A. pullulans-like isolates, it was confirmed that the earlier described variety A. pullulans var. melanogenum Hermanides-Nijhov is distinct from A. pullulans (de Bary) G. Arnaud 1918. In addition, two new varieties have been described: A. pullulans var. subglaciale Zalar, de Hoog & Gunde-Cimerman, and A. pullulans var. namibiae Zalar, de Hoog & Gunde-Cimerman. A. pullulans var. aubasidani Yurlova (Yurlova and de Hoog), which had been previously described due to its structurally unique polysaccharide, aubasidan, was synonymised with A. pullulans var. pullulans.
In 2008, the infraspecies classification of the varieties was retained in the redefinition of the species to prevent confusion linked to the use of new or additional epithets in applied fields . Nevertheless, substantial phenotypic differences in growth-temperature range, melanisation, and tolerance to NaCl had already been observed among the varieties. A. pullulans var. pullulans is found mostly in mildly osmotic environments, it is frequently associated with plants, and it can tolerate up to 17% NaCl (w/v), the highest salt concentration of all four of the A. pullulans varieties . A. pullulans var. melanogenum has been isolated mainly from oligotrophic, aqueous environments, and it grows at 37°C, while the other varieties cannot. On the other hand, A. pullulans var. subglaciale is unique for its psychrotolerant nature (it grows at 4°C) and its occurrence in glacial habitats in Svalbard (Norway). Finally, A. pullulans var. namibiae was named based on a single isolate from Namib Desert marble . Polymorphism is evident within A. pullulans, even within the varieties (Figure 1), and the species has been proposed as a model for the investigation of fungal phenotypic plasticity, because of the unique colony morphologies, the frequent changes in appearance and growth-temperature regimes, and the different nutrient requirements and morphological responses to light among its strains and varieties .
The genome sequence of one strain of each of the above-described varieties of A. pullulans was obtained: A. pullulans var. pullulans, A. pullulans var. melanogenum, A. pullulans var. subglaciale (type strain), and A. pullulans var. namibiae (type strain). The results of the genome analyses performed are presented and discussed in the following sections.
Results and discussion
The sequencing of the genomes of the four varieties of A. pullulans was aimed to uncover the genetic basis of their many interesting and useful traits. This included three aspects in particular: (i) their biotechnological potential (e.g., genes for extracellular enzymes, and for production of pullulan, aureobasidin A, and other compounds of interest); (ii) their exceptional stress tolerance (mainly focusing on genes associated with the synthesis of compatible solutes and water management, transport of alkali-metal ions and water, the high osmolarity glycerol [HOG] pathway, and melanin synthesis); and (iii) evidence in favour of describing these varieties as separate species.
Genome properties: sequencing and assembly of the genomes
Genome assembly and annotation statistics for the four Aureobasidium pullulans varieties
Value per A. pullulans variety*
Assembly length (Mbp)
Contig length total (Mbp)
Number of contigs
Contig L50 (kbp)
Number of scaffolds
Scaffold L50 (Mbp)
Number of scaffold gaps
Scaffolds gaps length (bp)
Percentage of scaffolds in gaps (%)
Number of repeat-covered regions
Length of repeat-covered regions (bp)
Percentage assembly covered by repeats (%)
GC content (%)
Number of genes
Protein length (amino acids, median)
Exon length (bp, median)
Gene length (bp, median)
Transcript length (bp, median)
Intron length (bp, median)
Number of genes with intron
Percentage of genes with an intron
Introns per gapped gene (median)
Intergenic distance (bp, median)
GC content of exons (%)
Genes with KEGG annotation [n, (%)]
Genes with KOG annotation [n, (%)]
Genes with Swissprot hit [n, (%)]
Genes with Pfam domain [n, (%)]
Genes with transmembrane domain [n, (%)]
Aureobasidium pullulans var. pullulans contains more repetitive sequences than the other A. pullulans varieties. However, the level of the repetitive sequences is low across all four of these varieties (0.78%-1.45%), and cannot explain the differences in their genome sizes. A. pullulans var. pullulans is the most halotolerant of the four investigated varieties , and it was isolated from hypersaline water. Therefore, it is interesting to note that investigations of Saccharomyces cerevisiae have revealed a tendency towards increased genome size as a response to high concentrations of salt [52, 53]. An increased genome size was also observed for the extremely halotolerant black yeast Hortaea werneckii, which has experienced recent whole-genome duplication . However, the increase in the genome size of A. pullulans var. pullulans is less extensive, as its genome is only 14.8% larger than the genome of its closest relative, A. pullulans var. subglaciale.
(i) Biotechnological potential of the Aureobasidium pullulansvarieties
Synthesis of pullulan
Pullulan is a linear α-D-glucan that is made of maltotriose units connected with α-1,6 linkages, while the glucose units within the maltotriose are connected with α-1,4-glycosidic bonds. This alternation of the bonds gives the molecule its distinctive physical properties: flexibility, solubility, adhesive ability, biodegradability, and the ability to form viscous solutions, and oxygen-impermeable and transparent fibres and films. Pullulan thus has broad application value . Despite its economic importance, relatively little is known about the biosynthesis of pullulan. Duan et al.  proposed a biosynthetic pathway in which the key enzymes for converting glucose units into pullulan were α-phosphoglucose mutase, uridine diphosphoglucose pyrophosphorylase, and glucosyltransferase.
There are single-copy genes for phosphoglucose mutase and uridine diphosphoglucose pyrophosphorylase in all four A. pullulans varieties. This is expected, as both of these enzymes catalyse important metabolic reactions, and they have been evolutionarily conserved. The enzyme phosphoglucose mutase catalyses a key step in hexose metabolism: the interconversion of glucose 1-phosphate and glucose 6-phosphate. The enzyme uridine diphosphoglucose pyrophosphorylase catalyses the reversible formation of uridine diphosphoglucose from glucose 1-phosphate and UTP. As described below, the predicted A. pullulans secretomes contain representatives of several Carbohydrate-Active enZYme (CAZy) database families , which include members with glucosyltransferase activities. Thus all four A. pullulans varieties contain all of the putative enzymes that were proposed to be involved in the biosynthesis of pullulan by Duan et al. .
Kang et al.  showed that disruption of the putative pullulan synthetase gene (pul) of A. pullulans [GenBank:AF470619] reduces its exopolysaccharide production to a pure β-glucan. Each of the predicted proteomes of the four A. pullulans varieties contains one similar protein (JGI Protein IDs: 349889, 3372, 53761, 64747, in A. pullulans var. pullulans, subglaciale, namibiae, and melanogenum, respectively; e-values between 10-22 and 10-27). All of these are predicted as being secreted. The low similarity scores can be attributed to the differences in annotation, as the pul gene prediction in GenBank contains an atypically long intron (almost 600 bp) and its second exon lies in a poorly conserved region (Additional file 1). When the annotation of the pul gene [GenBank:AF470619] was corrected to correspond to our annotation, it contained a one-nucleotide frameshift–insertion at nucleotide position 62 (Additional file 1, blue). After this was removed, the resulting protein-coding sequences and the predicted protein sequences were more than 70% identical. It is possible that the insertion in AF470619 was an artefact of the sequencing and that it led to the differences in the above-described gene annotation. If the pul gene is indeed involved in pullulan production, this observation is of significant importance.
Synthesis of aureobasidin A and siderophores
Strains of A. pullulans are used efficiently as biocontrol agents of post-harvest diseases . Not only numerous hydrolytic enzymes , but also various antimicrobial compounds, such as exophilin A , siderophores  and aureobasidin A , might have roles in the strong antagonistic effects of A. pullulans towards other species.
The antibiotic aureobasidin A is a cyclic nonadepsipeptide  that shows strong fungicidal activity, including against Candida species, Cryptococcus neoformans, and some Aspergillus species . Aureobasidin A has been shown to inhibit the phosphatidylinositol:ceramide phosphoinositol transferase that is involved in sphingolipid synthesis . The synthesis of aureobasidin A is catalysed by the 11659-amino-acid-long biosynthesis complex that is encoded by the intronless gene aba1 of the A. pullulans strain R106 .
The four sequenced A. pullulans varieties have not been tested previously for aureobasidin A production. However, their genomes do not include any homologues of the aba1 gene. This is not entirely surprising, as the production of aureobasidin A is not a universal trait of A. pullulans.
However, in a search for aureobasidin A synthase, a group of similar non-ribosomal peptide synthases was identified, although these proteins are substantially shorter than aureobasidin A synthase (1123–2920 amino acids) and are most similar to synthases for siderophores. A. pullulans var. melanogenum and A. pullulans var. namibiae have one putative synthase for siderophores, whereas A. pullulans var. subglaciale has two copies, and A. pullulans var. pullulans has three.
Siderophores are iron-chelating compounds with substantial biotechnological potential that are known to be produced by A. pullulans. Siderophores have been reported to act as antimicrobials . Their presence is also in line with the oligotrophic nature of A. pullulans, as they are beneficial under conditions of iron-depletion. Additionally, in other fungi, siderophores have roles in virulence, fungal–host interactions, and resistance to oxidative stress .
Secondary metabolite biosynthesis clusters
Fungi produce a multitude of low-molecular-mass compounds, known as secondary metabolites, which have roles in a range of cellular processes, such as transcription, development, and intercellular communication . Many of these compounds have important applications as antibiotics or immunosuppressants , and they are used in medicine or for plant protection . Genome mining investigations have indicated that the ability of fungi to produce secondary metabolites has been substantially underestimated, because many of the fungal secondary metabolite biosynthesis gene clusters are not expressed under standard cultivation conditions .
Secondary metabolite biosynthetic clusters for the four Aureobasidium pullulans varieties
Secondary metabolite biosynthetic cluster
Number per A. pullulans variety*
Type I PKS
Type III PKS
Total number of clusters
An unusually large spectrum of extracellular enzymatic activities has been described for A. pullulans, many of which are of considerable biotechnological interest [6–8]. Some of these, such as alkaline serine proteases, glucanases and chitinases, are also believed to have roles in the above-described antagonistic effects of A. pullulans against phytopathogenic fungi .
Predicted secreted proteins of the four Aureobasidium pullulans varieties
Predicted secreted proteins
Number per A. pullulans variety*
All four of these A. pullulans varieties also contain a large number of enzyme families, the members of which are known to be involved in the degradation of plant material, as designated by . These are discussed below.
The enzymes involved in the degradation of pectin are represented by three families of polysaccharide lyases (PL1, PL3, PL4), and the families of GH28 (polygalacturonases), GH78 (rhamnosidases), GH105 (rhamnogalacturonyl hydrolases) and the carbohydrate esterases CE8 (pectinesterases). The GH88 family (d-4,5-unsaturated β-glucuronyl hydrolases) is only present in A. pullulans var. namibiae, although it is not recognised as secreted. No representatives of the polysaccharide lyases from the PL9 and PL11 families were identified.
In all four of these A. pullulans varieties we identified five out of seven protein families with members that are involved in (among other activities) the degradation of cellulose (they contain the GH6, GH7, GH12, GH45, and AA9 families, but not the GH74 and GH94 families). A. pullulans var. namibiae contains all of the 15 families that contain hemicellulases (GH10, GH11, GH27, GH29, GH35, GH36, GH39, GH43, GH51, GH53, GH54, GH62, GH67, GH93, GH115), while the other three A. pullulans varieties lack the GH39 family (Additional file 3). This diversity is comparable to, and in some cases even larger than, other plant pathogenic fungi , including those from Dothideomycetes (Additional file 4) , and it is larger than in most of the saprophytic fungi, especially in the case of the hemicellulases (Additional file 4) .
Although the sequence-based families of carbohydrate-active enzymes frequently group together enzymes of varying substrate specificities, they are nevertheless considered to be good reporters of fungal lifestyles, especially when considering broad substrate categories, such as cellulose, hemicellulose or pectin . The richness of the enzyme families that are involved in the degradation of plant material is thus interesting in light of the epiphytic lifestyle of A. pullulans. When comparing individual varieties, proteins containing PFAM protein family domains (http://pfam.sanger.ac.uk/), which are characteristic for the GH3 (PF00933, PF01915), GH5 (PF00150), GH16 (PF00722), GH28 (PF01915), and GH43 (PF04616) families, are significantly enriched in A. pullulans var. pullulans compared to the other varieties (using the Computational Analysis of Gene Family Evolution [CAFE] software ), while the GH3 and GH5 families are depleted in A. pullulans var. melanogenum. This corresponds to the plant-related ecology of A. pullulans var. pullulans and to the largest phylogenetic distance between A. pullulans var. pullulans and A. pullulans var. melanogenum.
The numbers of genes and families of secreted proteases in the genomes of A. pullulans are also comparable to those found in the dothideomycete plant pathogens . All of the subfamilies of the MEROPS database of proteolytic enzymes , which are expected to efficiently digest proteins and/or to work in inhospitable environments of the extracellular matrix (A01, C13, G01, M35, M20, S08, S09, S10; ), are represented in these predicted secretomes. The only exception is C13 protease, which is recognised as being secreted only in A. pullulans var. melanogenum, whereas the homologues from the other varieties are not. The family of scytalidoglutamic peptidases (G01) contains even more representatives than in the analysed plant pathogens. Aspartic proteases (PF00026) and the A4 family (PF01828) are enriched in A. pullulans var. pullulans.
The number of lipases in the predicted secretomes is relatively low, although the number of cutinases, which are enzymes that are important for the degradation of the plant cuticle, is in the range observed in related plant pathogens  (Table 3). However, some of the putative cutinases are not designated as secreted with the criteria used in the present analysis, so the actual number of these enzymes might be even higher.
Tannases are another group of enzymes that are significantly enriched within the genome of A. pullulans var. pullulans. Tannases, or tannin acyl hydrolases (EC, 126.96.36.199), catalyse the hydrolysis of ester bonds in the hydrolysable tannins and gallic acid esters. These enzymes are used industrially as catalysts in the manufacture of gallic acid, and they also have potential use in beverage and food processing .
Given the previously recognised substantial biotechnological potential of the large spectrum of the extracellular enzymatic activities of this fungus [6–8], which is confirmed here by the predicted secretome, their application should be brought to a new level with the availability of these genomic data. The wide array of enzymatic activities should also be considered when investigating the pathogenicity potential of A. pullulans. As already noted by Chan et al. , lipases, phospholipases, proteases and β-lactamases are important virulence factors that might have roles in A. pullulans as an emerging opportunistic human pathogen.
Biodegradation of plastic and aromatic compounds
Black yeast, including A. pullulans, are known for their ability to degrade aromatic pollutants [74, 75] and plastic [76, 77]. In the A. pullulans AY4 strain, Chan et al.  reported the presence of genes coding for 2-monooxygenases and catechol dioxygenases (biodegradation of aromatic substances) and depolymerase (biodegradation of plastic). There are two families of 2-monooxygenases in these four A. pullulans varieties. The larger of these has six representatives in each of A. pullulans var. pullulans and A. pullulans var. subglaciale, five in A. pullulans var. melanogenum, and four in A. pullulans var. namibiae. The second family contains only one representative in A. pullulans var. pullulans and A. pullulans var. subglaciale. A family of proteins with catechol 1,2-dioxygenase activity is also present, with four representatives in A. pullulans var. melanogenum and five in each of the other three A. pullulans varieties. There are single-copy genes for poly-β-hydroxybutyrate depolymerase, an enzyme that can degrade polyhydroxyalkanoates , in all four of these A. pullulans varieties. The presence of these genes underlines the considerable potential of A. pullulans for use in bioremediation or biodegradation processes.
Major facilitator superfamily sugar transporters
The major facilitator superfamily (MFS) of sugar transporters is the largest group of secondary transmembrane carriers. This consists of 74 protein families, each of which is specialised for the transport of a certain type of substrate. They have different transport modes, from uniport, to solute:cation (H+ or Na+) symport and/or solute:proton or solute:solute antiport. Most of these have 400 to 600 amino acyl residues, with two six-transmembrane segment repeat units . According to the Transporter Classification (TC) database , two representatives in this group, the MFS family (TC no. 2.1) and the glycoside-pentoside-hexuronide:cation symporter family (TC no. 2.2), contain sugar-specific families of secondary active transporters. In Eukarya, only two families of MFS sugar transport substrates are known, of which the sugar porter family (TC no. 2.1.1) is the largest of all of the MFS families, with several hundred known members. The sugar porter family is diverse with respect to their substrate specificity, and they can work either as uniporters or as cation symporters. The other familiy of sugar transporters is the sialate:H1 symporter family (TC no. 2.1.12) .
Secondary active sugar transporters within the major facilitator superfamily (MFS) in the four Aureobasidium pullulans varieties
Transporter classification (TC) no.
Number per A. pullulans variety and yeast*
Major facilitator superfamily
Sugar porter family
Quinate: H+ symporter or
Myo-inositol: H+ symporter or
Lactose, galactose: H+ symporter or
Maltose: H+ symporter or
General α-glucoside: H+ symporter
Fructose: H+ symporter
Glycerol: H+ symporter
High affinity glucose transporter (probably H+ symporter)
Xylose: H+ symporter
High affinity glucose: H+ symporter (monosaccharides including xylose)
Low affinity glucose: H+ symporter
Galactose: H+ symporter
MFS monosaccharide transporter
MFS sugar transporter
Low-affinity glucose transporter Hxt1/3
High-affinity glucose transporter Hxt2
Hexose uniporter Hxt 10
Low affinity, constitutive, glucose (hexose; xylose) uniporter Hxt4
Hexose transporter with moderate affinity for glucose Hxt5
High affinity, glucose-repressible, glucose (hexose) uniporter Hxt6/7
Hexose transporter Hxt8
Hexose transporter Hxt9/11
Hexose transporter Hxt13/17
Hexose transporter Hxt14
Hexose transporter Hxt15/16
Galactose, glucose uniporter Gal2
Vacuolar protein sorting-associated protein 73 or
Probable metabolite transport protein YBR241C or
Probable metabolite transport protein YFL040W
The sialate: H + symporter family
Lactate/pyruvate: H+ symporter
The glycoside-pentoside-hexuronide: cation symporter family
Maltose/sucrose H+ symporter Sut1
The high numbers of sugar transporters correspond to the nutritional diversity of A. pullulans, while the differences between the four A. pullulans varieties might reflect their ecological preferences. A. pullulans var. pullulans is most frequently associated with plants, and it thus has access to various sugars of plant origin. On the other hand, A. pullulans var. melanogenum has significantly fewer sugar transporters, and it is associated more with freshwater habitats, where simple sugars are not as readily available. It is also separated from A. pullulans var. pullulans by the largest evolutionary distance of all of these A. pullulans varieties. The most notable difference between these varieties is in the number of the putative maltose or α-glucoside: H+ symporters (TC no. 2.A.1.1.10 or 2.A.1.1.11), as A. pullulans var. melanogenum has almost half the number of copies compared to the other varieties.
Components of the high-osmolarity glycerol pathway
The mitogen-activated protein kinases (MAPKs) are involved in many cellular processes, such as stress responses and regulation of differentiation and proliferation, and they are highly conserved in eukaryotes. This system consists of three kinases that phosphorylate one another in a signalling cascade. One of the best-characterised is the high-osmolarity glycerol (HOG) pathway, which is a branched MAPK signal-transduction system, the physiological role for which in fungi is primarily the mediation of cellular adaptation to increased osmolarity of the surrounding medium. The HOG pathway is also required for adaptation to other stress conditions, such as oxidative, heavy-metal, and hot or cold stress  and in the virulence of pathogenic fungi .
The HOG pathway in fungi consists of two major signalling modules. One is the MAPK module that comprises MAPK, MAPK kinase (MAPKK), and MAPKK kinase (MAPKKK). The other module is a two-component phospho-relay system that is composed of hybrid sensor kinases, a histidine-containing phosphotransfer protein, and response regulators; this system senses and relays environmental signals, and subsequently activates the HOG pathway [84, 85].
As expected, the main components  in the four A. pullulans varieties are very similar to those in other fungi. These genomes contain one homologue of each kinase in the MAPK module (Hog1 MAPK, Pbs2 MAPKK, and two MAPKKKs, Ssk2/22 and Ste11), one Ste50 kinase homologue, and one homologue of all of the components in two sensory branches Sln1 and Sho1 (Sln1, Ypd1, Ssk1, Sho1). Interestingly, A. pullulans var. pullulans and A. pullulans var. namibiae have two predicted homologues of Ste20 or Cla4 kinases (activators of Ste11 MAPKKK), while the other two varieties have only one. One of the duplicated homologues is missing the pleckstrin-homology domain (PF00169). This domain has roles in the recruitment of proteins to different cellular membranes, and thus enables the proteins to interact with other components of the signal-transduction pathways . Whether the observed difference has a role in the polyextremotolerance of A. pullulans remains a question for future studies.
Genes involved in the biosynthesis of compatible solutes
Stress tolerance is frequently associated with accumulation of small organic molecules that can act as protective compounds against different stress factors (for review, see ). In A. pullulans, exposure to high temperatures or high salt concentrations, or a combination thereof, leads to increased intracellular concentrations of the polyhydroxy compounds trehalose, mannitol and glycerol . The concentrations of trehalose in A. pullulans increase in heat-stressed cells, and in simultaneously heat-stressed and salt-stressed cells, but not in cells subjected to salt stress alone. Mannitol increases under all of the stress conditions that have been examined, while an increase in intracellular glycerol has only been detected in salt-stressed cells . Indeed, in response to a saline stress, glycerol is the most abundant compatible solute (Kogej, unpublished data). In other microbes, glycerol, trehalose and other polyols are also very important for life in extremely cold environments . This might also be the case in the psychrotolerant A. pullulans, however, the role of compatible solutes at low temperatures has not yet been studied in this species.
Glycerol is synthesised from dihydroxyacetone phosphate, a glycolytic intermediate, via two reaction steps that are catalysed by nicotinamide adenine dinucleotide (NAD)-dependent glycerol-3-phosphate dehydrogenase (Gpd) and glycerol-3-phosphatase (Gpp) . Both of the genes that encode these enzymes required for glycerol biosynthesis are present in the genomes of all four of these A. pullulans varieties, as single-copy genes. The yeast S. cerevisiae and the halophilic basidiomycete Wallemia ichthyophaga each have two copies of Gpd , while the extremely halotolerant H. werneckii has four, two of which can be attributed to a whole-genome duplication event . The predicted Gpd proteins from A. pullulans lack the N-terminal PTS2 sequence that is important for peroxisome localisation, a trait also noted in other halotolerant fungi . It was proposed that the constant cytosolic localisation of the Gpd1 homologues is advantageous for the organisms that live in extremely saline environments , as the osmoprotective role of Gpd is dependent on the cytosol and nuclear fractions of this protein .
The biosythesis of trehalose is a two-step process in which glucose 6-phosphate and UDP-glucose are first converted to α,α-trehalose 6-phosphate by trehalose-6-phosphate synthase (Tps), and then to trehalose by trehalose-6-phosphate phosphatase (Tpp) . In Ascomycota, mannitol is also synthesised in two steps: fructose 6-phosphate is first reduced to mannitol 1-phosphate by NAD-dependent mannitol-1-phosphate dehydrogenase (Mpd), and then dephosphorylated to mannitol, by mannitol-1-phosphate phosphatase (Mpp) . One copy of trehalose-phosphatase was found in each of the four investigated A. pullulans varieties, and two relatively dissimilar (~30% identical amino acids) trehalose synthases were identified. Two proteins in each variety have predicted mannitol dehydrogenase activities. One is similar to the homologues that are crucial for mannitol biosynthesis in other fungi, while the other is probably involved in the first step in the catabolism of mannitol . Interestingly, no homologues of known fungal mannitol-1-phosphate phosphatases were found in any of these four A. pullulans genomes.
In the genomes of all four of the A. pullulans varieties, we identified a putative D-xylose reductase (cluster 19, 48900/38638, 70467, 36428 and 36838 for A. pullulans var. pullulans, var. subglaciale, var. namibiae and var. melanogenum, respectively). D-xylose reductase converts xylose to xylitol, a sweetener compound. The production of xylitol is a commercially interesting process [96, 97], which has rarely been contemplated in A. pullulans.
Melanin biosynthesis genes
The cells of A. pullulans produce a black pigment that has long been known to be 1,8-dihydroxynaphthalene (DHN)-melanin . Melanin is a high-molecular-weight, dark brown or black pigment that can be produced by numerous fungi , and it has a known protective role under various stress conditions (for review, see ), and also under hypersaline conditions . Melanin is located in the fungal cell wall, either enmeshed within the structure of the cell wall, or as its outermost layer . Microbes predominantly produce melanin pigment via tyrosinases, laccases, catecholases, and the polyketide synthase pathway .
Dothideomycetes mostly produce DHN-melanin via the polyketide synthase pathway [104, 105]. DHN-melanin biosynthesis starts with polyketide synthase, with acetyl coenzyme A or malonyl coenzyme A as a precursor. The polyketide synthase produces 1,3,6,8-tetrahydroxynaphthalene, which is reduced by hydroxynaphthalene reductase to form scytalone. Dehydration of scytalone by scytalone dehydratase forms 1,3,8-trihydroxynaphthalene. In turn, 1,3,8-trihydroxynaphthalene reductase converts 1,3,8-trihydroxynaphthalene to vermelone, which is further dehydrated by scytalone dehydratase, to 1,8-DHN. Subsequent steps are believed to involve dimerisation of the 1,8-DHN molecules, followed by polymerisation that is catalysed by p-diphenol oxidase .
One melanin polyketide synthase gene was identified in each of these four A. pullulans varieties. This is in agreement with findings in the other DHN-melanin–producing fungi, which also have one melanin polyketide synthase [107–109]. Additionally, there are one of each of the trihydroxynaphthalene reductase-like gene and tetrahydroxynaphthalene reductase-like gene. Each of the two reductases from these four A. pullulans varieties falls within either a trihydroxynaphthalene-like or a tetrahydroxynaphthalene-like group, as has been observed in other melanised fungi. One of the traits not expected is the presence of two scytalone dehydratase-like genes in all of the investigated varieties except A. pullulans var. melanogenum. These duplicated genes form two distinct phylogenetic groups. Furthermore, our phylogenetic analysis of the scytalone dehydratases from these four A. pullulans varieties and the scytalone dehydratase sequences collected from the NCBI-NR protein database show the presence of two scytalone dehydratases in several other fungal species as well (Additional file 5). These duplicated genes form two distinct phylogenetic groups: the larger one (cluster I) contains scytalone dehydratases from known DHN-melanin–producing fungi, and proteins 67660, 42009, 345205 and 70470 (from A. pullulans var. melanogenum, var. subglaciale, var. pullulans, and var. namibiae, respectively). The second cluster has two subclusters: one (IIa) contains the three duplicate proteins from A. pullulans var. subglaciale, A. pullulans var. pullulans and A. pullulans var. namibiae (676013, 282919, 48527), together with both scytalone dehydratase proteins from Fusarium, Botryotinia and Metarhizium. The other subcluster (IIb) contains proteins that produce bluish-green pigment through the DHN-melanin pathway , and it has no homologues in the four A. pullulans varieties (Additional file 5). This is the first report of several scytalone dehydratase genes and proteins in a single fungal species. Thus, we have identified here genes that possibly encode essential components of the DHN-melanin biosynthesis pathway in all four of these A. pullulans varieties.
A large family of major intrisic proteins (MIPs), which are membrane-channel proteins that are selective for the transport of water (aquaporins) or water plus glycerol (aquaglyceroporins), has been found in diverse life forms . All of the aquaporins are transmembrane proteins with six transmembrane domains, and with their N-terminus and C-terminus in the cytosol. These are water channels that are 109-fold faster compared to ion channels and transporters. The presence of solute-permeable aquaporins represents a new concept in terms of absorption, not only because of the rate of the process, but also because of its osmotic implications .
Orthodox aquaporins mediate rapid and selective fluxes of water across biological membranes, and hence they have important roles in the osmoregulation of cells and organisms. Aquaglyceroporins, on the other hand, facilitate transmembrane transport of small uncharged molecules, like polyols, urea, arsenite, and many more, thereby having roles in nutrient uptake, osmoregulation, and probably other processes [113, 114]. Petterson et al.  divided the fungal aquaporins into four groups: one group of orthodox aquaporins, and three groups of aquaglyceroporins (Fps1-like, Yfl054-like, and aquaglyceroporins that do not fit into these groups). A later phylogenetic analysis of 229 major intrisic fungal proteins by Xu et al.  also classified the aquaporin-like genes into four clusters, delineated by functionally characterised major intrisic fungal proteins: orthodox aquaporins, aquaglyceroporins, facultative fungal aquaporins, and X intrinsic proteins.
Alkali-metal cation transporters
As a halotolerant species, A. pullulans has to maintain its intracellular cation homeostasis during changing and occasionally high environmental salinities. Under such conditions, the maintenance of a high and stable K+ content and the elimination of toxic Na+ ions is crucial for survival [117, 118]. However, it is not clear if and how the number of alkali-metal cation transporter genes is correlated with the halotolerance of an organism. The genome sequencing of the extremely halotolerant black yeast H. werneckii revealed an extensive enrichment of genes that encode transporters of alkali-metal cations . These appear to allow H. werneckii to maintain the very low internal Na+ concentrations that are observed even at high external NaCl concentrations . On the other hand, genome sequencing of the most halophilic fungus known, the basidomycete W. ichthyophaga, uncovered only modest numbers of ion-transporter genes, and additionally showed that their transcription is relatively low and non-responsive to different salt concentrations .
Numbers of specific types of alkali-metal cation transporters of the four Aureobasidium pullulans varieties
Substrate specificity/main function
Number per A. pullulans variety*
Antiporter with H+
K+ uptake (Hak)
K+ uptake (Acu)
Na+ (and Li+) efflux
Antiporter with H+
Na+ (and K+) efflux
Ca2+ and Mn2+ transport into Golgi apparatus
Na + (and K+) efflux
Na+, Ca2+/H+ exchange
V-type ATPase (4 subunits)
H+ in vacuole
Depleting cytosol of Ca2+ ions
S. cerevisiae contains two typical K+-transport systems: Trk1 and Trk2, which are high-affinity channels for K+ uptake , while the Tok1 membrane-depolarisation-activated channel is for K+ efflux . Ena1, Ena2 and Ena5 P-type ATPases [121–123] and the Nha1 antiporter, which is driven by a proton gradient, can also contribute to K+ efflux, although they are primarily known as exporters of Na+. In A. pullulans, both Trk and Tok channels are duplicated (Figure 4). A. pullulans var. pullulans contains another homologue of Trk that belongs to cluster 2, which is still present in H. werneckii and C. heterostrophus, but not in other fungi analysed. The same is true for A. pullulans var. subglaciale, although its homologue is uncharacteristically short (it contains only 137 amino acids, compared to the average of 808 for the other A. pullulans Trk proteins).
Unexpectedly, in addition to Trk, all of these four A. pullulans varieties contain two less-common K+ uptake systems (Table 5): the Hak H+/K+ symporters, which under some conditions appear to act as Na+/K+ symporters ; and the Acu P-type ATPases , which are encoded by not one, but two relatively divergent genes in each of these genomes (Figure 4). Both of these types of transporters are present in only a few other Dothideomycetes. According to the predicted phylogeny, it is parsimonious to assume that both of these K+ uptake systems were present in the ancestor of H. werneckii and were lost later in the evolution of this by-far-the-most halotolerant of all of the species analysed. Therefore, the primary purpose of the redundancy of these K+ uptake systems might not be adaptation to salt.
In other groups of transporters two duplications were observed in the above-mentioned Nha Na+/H+ antiporter, and one duplication in the Ena Na+ exporters (Figure 4). The smaller number of Ena genes compared to Nha (Table 5) has already been noted for H. werneckii. Two duplication events and the subsequent loss of genes (Figure 4) led to different numbers of Pho Na+/Pi symporters in these individual A. pullulans varieties (Table 5). In salty environments, the Na+ gradient possibly represents an alternative energy source for the influx of Pi.
The most halotolerant of all of the species that have been analysed, H. werneckii, contains the highest number of transporter genes (30 genes), although it should be noted that the number is doubled due to the recent whole-genome duplication of H. werneckii. The second largest number of transporter genes is present in C. heterostrophus (19 genes), which is not a halotolerant species (although, of note, it has not been studied in this respect), and is not found in hypersaline habitats. The high number of plasma-membrane K+ transporters of A. pullulans (seven or eight transporters across the four varieties) is in even starker contrast with the most halophilic fungus known, the basidiomycete W. ichthyophaga. The latter has only one Trk-encoding gene, and has no homologues of the Tok, Hak or Acu transporters .
When comparing only these four A. pullulans varieties, the most halotolerant A. pullulans var. pullulans has neither the highest number nor the largest diversity of transporters. Furthermore, although in S. cerevisiae the Ena pumps are the major determinant of salt tolerance (for review, see ), A. pullulans var. pullulans even has one Ena gene homologue less than the other A. pullulans varieties. Although the number of alkali-metal cation-transporter genes in A. pullulans is higher than in many related fungi, this difference is not particularly outstanding. It appears that the diversity and the total number of alkali-metal cation-transporter genes are (at best) only approximate predictors of halotolerance.
Only one gene for the plasma-membrane proton-exporting ATPase Pma is present in each of these four A. pullulans varieties. In contrast, H. werneckii has four copies of this gene, and their expression is salt responsive . Pma pumps are responsible for establishing the electrochemical gradient of protons across the plasma membrane, and through this, they supply the energy to the secondary active symporters and antiporters [125, 126], among which there are also those responsible for maintaining ion homeostasis under hypersaline conditions.
In addition to the Pma pumps, other transporters might be responsible for the building of the electrochemical gradient of protons. Type 1 microbial rhodopsins can function as proton pumps to establish the electrochemical gradient for ATP production. This process has already been demonstrated in another dothideomycetous fungus, Leptosphaeria maculans. An alternative way to build this proton-motive force would not only be useful in hypersaline conditions, as has already been proposed for the extremely halotolerant H. werneckii, but it might also be of importance in the light-exposed, oligotrophic biology of A. pullulans; e.g., to rapidly resume metabolism after cryptobiosis. The search for bacteriorhodopsins revealed at least two gene copies in the genomes of all four of the A. pullulans varieties (with three copies in A. pullulans var. melanogenum). The copies within these genomes differ in terms of their spliceosomal intron numbers and positions. According to our reconstructed phylogeny, these bacteriorhodopsin copies with one and two introns formed monophyletic groups (orthologues; see Additional file 9), and appear to represent old duplications before the split of the lineages within A. pullulans. However, A. pullulans var. melanogenum contains an additional third copy, suggesting either a recent duplication in this A. pullulans variety or a loss of a paralogue in the other A. pullulans varieties.
Although bacteriorhodopsins from A. pullulans have not yet been used for biotechnological purposes, they can be considered for the many different applications that this family of proteins has; e.g., mainly for optical appliances, but also for therapeutic/medical applications and research, including in optogenetics and bioelectronics [129, 130].
(iii) Description of new Aureobasidiumspecies
Elevation of these four A. pullulans varieties to the species level
A revised concept for the species A. pullulans is presented here. These four varieties of A. pullulans are elevated to the species level based primarily on the phylogenetic analysis of their whole genome sequences, as well as on their phylogenetically relevant genes.
Aureobasidium melanogenum (Hermanides-Nijhof) Zalar, Gostincar, Gunde-Cimerman, stat. nov. MycoBank MB 807698.
Basionym: Aureobasidium pullulans var. melanogenum Hermanides-Nijhof, Stud. Mycol. 15: 161, 1977. MycoBank MB352628.
Aureobasidium subglaciale (Zalar, de Hoog & Gunde-Cimerman) Zalar, Gostincar, Gunde-Cimerman, stat. nov. MycoBank MB 807700.
Basionym: Aureobasidium pullulans var. subglaciale Zalar, de Hoog & Gunde-Cimerman, Stud. Mycol. 61: 33, 2008. Mycobank MB512380.
Aureobasidium namibiae (Zalar, de Hoog & Gunde-Cimerman) Zalar, Gostincar, Gunde-Cimerman, stat. nov. MycoBank MB807701.
Basionym: Aureobasidium pullulans var. namibiae Zalar, de Hoog & Gunde-Cimerman, Stud. Mycol. 61: 34, 2008. Mycobank MB512381.
For the synonyms, see the previously published taxonomic study of Zalar et al. .
Taxonomic markers for distinguishing the new Aureobasidium species
Species/varieties within the genus Aureobasidium have most often been distinguished based on multilocus analyses of different regions of the rDNA gene clusters. The ITS region has been most commonly used [19, 47–49], although it is recognised as insufficient in certain genera . Based on the ITS rDNA, the distinction of Aureobasidium pullulans was possible, but the resolution based on other household genes was better . Therefore, the different household genes used earlier in the phylogenetic analyses of Aureobasidium and of other genera were compared here. The selected genes encode for the following proteins: actin (ACT), β-tubulin (BTUB), calmodulin (CAL), chytin synthase (CHS), NAD-dependent glycerol-3-phosphate dehydrogenase (GPD), mini-chromosome maintenance proteins essential for initiation of eukaryotic genome replication (MCM7), RNA polymerase 2, the largest subunit (RPB1), RNA polymerase 2, the second-largest subunit (RPB2), and translation elongation factor 1α (TEF1α). CAL and CHS are present in several different copies in these four Aureobasidium taxa and are thus less appropriate for use in phylogenetic analyses. Other markers were identified as single-copy genes in the genomes.
Household genes established as fungal phylogenetic markers
Number of exons/introns
Region used as marker/length/number of variable positions (vp) in Aureobasidium species
Previous use in Aureobasidium taxonomy/phylogeny (analysed region and reference)
Exon 5/722 bp/39 vp
Exon 2/1115 bp/272 vp
Exon 2/2285 bp/562 vp
Exon 3/exon 5/700 bp/103 vp
Exon 3/1139 bp/87 vp
Exon 3/4898 bp/1123 vp
Exon 2 to partial exon 3 (644 bp) 
Exon 1/3608/803 vp
The sequences of ACT, GPD and MCM7 have not yet been used in phylogenetic studies of Aureobasidium, although they represent appropriate candidate regions for such studies . The positions of the introns in these genes are conserved in all four of the investigated genomes. To evaluate the usefulness of the listed gene exons and introns as barcodes for species identification or as phylogenetic marker loci, further analyses based on a larger set of Aureobasidium strains are required.
Pairwise distances for the set of single-copy household genes in the four Aureobasidium species
Pairwise distances of genes
The genome sequences generated here also allowed the correct identification of two other ‘A. pullulans’ strains for which the genomes had been sequenced in the past and are publicly available. Strain ATCC 62921 (http://genome.fungalgenomics.ca/), which is of an unknown origin, but it has been used in numerous industrial applications, is confirmed as A. pullulans. The sequenced strain AY4, which was isolated from a human patient , is identified as A. melanogenum according to the above-presented species concept. The new species A. melanogenum also includes all other infection-causing ‘A. pullulans’ strains. Of the other commonly used ‘A. pullulans’ strains, strain NRRL Y-6220 is one of the commonly used strains for pullulan production, and it is related to A. melanogenum. On the other hand, the strain R106, which is used for aureobasidin A production, could not be reliably classified with the publicly available sequences, since in the phylogenetic analyses it was placed in a lineage separate from the here-described four species. The two Aureobasidium strains that are sold as commercial agricultural biocontrol products (DSM14940, DSM14941) remain in the species A. pullulans – for this purpose it is important to note that the species A. pullulans does not contain any potentially opportunistic pathogens, as these are placed in the newly described A. melanogenum.
Homothallic MAT loci and evidence for sexuality
A possible teleomorph of A. pullulans, Columnosphaeria fagi (H.J. Huds.) M.E. Barr, was suggested on the basis of the ITS sequence similarity , although this was not further confirmed either by a multilocus analysis or by a single-spore re-isolation of the fungus. Due to the lack of the recombination between the four previous A. pullulans varieties, it was concluded that A. pullulans is strictly clonal, although recombination within the four phylogenetic clusters of A. pullulans strains (at that time recognised as varieties) was not definitively ruled out . However, the redefinition of these varieties as the four Aureobasidium species (see above) implies that the lack of the recombination between the varieties might be simply a consequence of interspecific reproduction barriers and not a sign of asexuality. The genome sequences of the four reclassified Aureobasidium species reveal additional details about their possible mating strategies.
The sexual reproduction of ascomycetous fungi is orchestrated by the presence of different arrangements of mating-type genes that encode the key transcription factors at one or more of the MAT loci: the MAT1-1-1 gene encodes a protein with an alpha1 domain, and the MAT1-2-1 gene encodes a protein with a high mobility group (HMG) box .
There is an open reading frame of an unknown function between the MAT1-1-1 and MAT1-2-1 genes in three of the studied genomes, with the genome of A. subglaciale being the exception. In A. subglaciale, this open reading frame is still present as a pseudogene, but due to a substantial degeneration, its structure can be discerned only approximately. In Sordariomycetes and Eurotiomycetes, the MAT locus frequently lies between the homologues of the S. cerevisiae APN2 and SLA2 genes [54, 137, 139]. In the here investigated Aureobasidium species, this is the case for the APN2 homologues (Figure 6); while the homologues of SLA2 are located on the same scaffold, they are separated from the MAT genes by a substantial distance (70 to 500 kbp).
The single-copy homologues of other genes involved in mating were found outside of the mating loci, including G-protein-coupled receptor for alpha-factor (Ste2 in S. cerevisiae) and a-factor (Ste3), and the transcription factor that activates genes involved in mating or pseudohyphal/ invasive growth pathways (Ste12). Out of nine meiosis-specific genes listed by Malik et al. , five were identified (homologues of S. cerevisiae SPO11, MSH4, MSH5, MER3, REC8). Interestingly, the homologue of MER3, a DNA helicase that is involved in crossing over , is present in three copies in A. melanogenum, and in two copies in A. pullulans, A. subglaciale and A. namibiae.
The roles of absent genes might be provided by alternative proteins. For example, the species A. pullulans contains two homologues of Rad51, which is a protein involved in the repair of double-strand DNA breaks. This might replace Dmc1, which performs the same function, but only in meiosis. Furthermore, the absence of the DMC1 homologue in sexual fungi is not unusual .
According to their genetic design, all four of the Aureobasidium species can be considered as homothallic (self-compatible, or self-fertile), which implies that a single strain is capable of sexual reproduction even in an axenic culture. However, the existence of mating genes is not conclusive evidence for sexual recombination, as these genes might also have other roles, as has been reported in other species of fungi . As a teleomorph of the previously recognised A. pullulans has not been described yet, the ability of these four Aureobasidium species to undergo sexual reproduction needs to be confirmed experimentally.
Aureobasidium pullulans, as it was previously defined, was considered a species of exceptional adaptability. This was reflected in its great phenotypic plasticity, polyextremotolerance, and the ability to survive across a wide variety of habitats. A part of this diversity can be explained by the relatively wide definition of this species. However, with the accumulation of new knowledge, it had become clear that smaller, well-defined groups of strains can be recognised within A. pullulans. These were initially described as varieties , but this comparison of their genomic data shows that the differences between these varieties are large enough to justify their redefinition as four separate Aureobasidium species: A. pullulans, A. melanogenum, A. subglaciale and A. namibiae.
Even single strains of A. pullulans can show wide adaptability, and this is reflected in the great diversity of several groups of genes. This diversity can reach, or even surpass, the gene diversity of more specialised related species, even though almost all of them have larger genomes. The investigated genomes of these four Aureobasidium species contain an abundance of different families of extracellular enzymes, which are important for the degradation of the plant cell-wall material. Their diversity is comparable to the phylogenetically related plant pathogens. The differences in the abundance of the extracellular enzymes and also in the higher numbers of sugar transporters reflect the ecological preferences of the four studied Aureobasidium species, especially for the plant-associated lifestyle of A. pullulans. Furthermore, for example, all four of the Aureobasidium genomes presented here have genes for the alternative K+ importers Hak and Acu, which are not found even in the extremely halotolerant H. werneckii or in the halophilic W. ichthyophaga. Also, when considering the mating strategy, these four Aureobasidium species appear to be more flexible than the specialised H. werneckii and W. ichthyophaga: the four Aureobasidium genomes all contain a homothallic mating type locus, whereas H. werneckii is heterothallic, and W. ichthyophaga lacks the mating genes altogether.
The genome sequences of the described four Aureobasidium species are expected to facilitate the exploitation of the substantial biotechnological potential of these fungi. For most of the interesting genes, at least four different genome copies are now available for study, with potentially different traits of the proteins that they encode. Additionally, we propose the elimination of the decade-old annotation error in the putative pullulan synthase gene. Furthermore, the genes for the degradation of a variety of compounds from plant polysaccharides to plastics and aromatic compounds should lead to more efficient industrial use of these Aureobasidium species, as well as help in the selection or construction of the best strains for any given purpose.
The polyextremotolerant character of the investigated Aureobasidium species indicates that they have efficient stress-combating mechanisms, some of which have been described here. These mechanisms can potentially be used in the future to improve the stress tolerance of economically important microorganisms or plants. Last but not least, as the opportunistic human pathogens belong only to A. melanogenum, the phylogenetic redefinition of the previously recognised A. pullulans will facilitate the identification of the more problematic strains and help in limiting their use in applications in which they might represent potential health risks for workers or consumers.
Strains, growth conditions and microscopy
Aureobasidium pullulans var. pullulans (EXF-150, CBS 100280) was isolated from the hypersaline waters of the Sečovlje solar saltern (Slovenia). Aureobasidium pullulans var. melanogenum (EXF-3378, CBS 110374) was isolated from a public fountain in Bangkok (Thailand). Aureobasidium pullulans var. subglaciale (type strain; EXF-2481, CBS 123387) was isolated from the subglacial ice of the Kongsvegen glacier on Spitsbergen (Svalbard, Norway). Aureobasidium pullulans var. namibiae (type strain; EXF-3398, CBS 147.97) was isolated from dolomitic marble in the Namib Desert (Namibia). All of these strains are preserved in the Ex Culture Collection of the Department of Biology, Biotechnical Faculty, University of Ljubljana (Infrastructural Centre Mycosmo, MRIC UL) and at the Centraalbureau voor Schimmelcultures (CBS, Utrecht, The Netherlands).
These cultures were grown at 28°C in a rotary shaker (180 rpm) in the defined yeast nitrogen base medium (YNB) of 0.17% (w/v) yeast nitrogen base, 0.08% (w/v) complete supplement mixture (both Qbiogene), 0.5% (w/v) ammonium sulphate, and 2% (w/v) glucose, in deionised water. The pH was adjusted to 7.0 prior to autoclaving. For the isolation of RNA for transcriptome sequencing of A. pullulans var. pullulans and A. pullulans var. subglaciale, which was used for gene annotation, a variety of conditions was used: (1) YNB; (2) YNB with additional 10% (w/v) NaCl; (3) YNB with 2% (w/v) raffinose or (4) 31.3% (w/v) glycerol instead of glucose; (5) YNB with 49% (w/v) sorbitol; (6) yeast peptone glucose medium (1% [w/v] yeast extract, 2% [w/v] peptone, and 2% [w/v] glucose, in deionised water); (7) minimal medium (2% [w/v] glucose, macroelements [w/v]: 0.6% NaNO3, 0.15% KH2PO4, 0.05% MgSO4 × 7H2O, 0.05% KCl, microelements [w/v]: 0.01% EDTA, 0.0044% ZnSO4, 0.001% MnCl2 × 4H2O, 0.00032% CoCl2 × 6H2O, 0.00032% CuSO4 × 5H2O, 0.00022% (NH4)6Mo7O24 × 4H2O, 0.00147% CaCl2 × 2H2O, 0.001% FeSO4 × 7H2O, in deionised water); (8) minimal medium with 2% (w/v) cellulose instead of glucose; (9) YNB and an incubation temperature of 12°C; (10) YNB and incubation until mid-exponential phase followed by heat shock at 42°C for 45 min; (11) YNB and incubation until stationary phase; (12) YNB and incubation until mid-exponential phase, followed by addition of 0.32 mM H2O2 for 5 min; and (13) YNB and incubation until mid-exponential phase, followed by centrifugation and further resuspension either in YNB at pH 2 and a 4-h incubation (one generation), or (14) in Tris base at pH 10 for 30 min. Unless otherwise noted, the cells were harvested in the mid-exponential growth phase, using a 10-min centrifugation at 5000× g, and immediately frozen in liquid nitrogen.
For photography of the macromorphology, the four A. pullulans varieties were grown on malt extract agar (MEA) according to Blakeslee, containing 2% malt extract, 0.1% peptone, 2% glucose and 2% agar (all Merck), for 30 days at 24°C.
For microscopic images, slide cultures were prepared on MEA blocks in moist chambers, and they were incubated for 7 days at 25°C in the dark. The cover slips were carefully removed and mounted in 60% lactic acid for microscopic observation under the Olympus BX51 microscope, using an Olympus DP12 camera, and the DP software.
DNA and RNA isolation
Cells frozen in liquid nitrogen were homogenised using a pestle and a mortar. Their DNA was isolated using the phenol/ chloroform/ isoamyl alcohol method, modified for DNA isolation from filamentous fungi, as described previously . The quality and quantity of the DNA was evaluated spectrophotometrically with NanoDrop 2000 (Thermo Fisher Scientific, USA), and on standard 1% agarose gel electrophoresis, with the molecular weight also checked.
The RNA was isolated using TRI Reagent (Sigma-Aldrich, Germany), according to the manufacturer instructions. Possible DNA contamination was removed with DNAse I (Thermo Fisher Scientific – Fermentas, Lithuania), and the integrity and purity of the RNA was evaluated spectrophotometrically with NanoDrop 2000 (Thermo Fisher Scientific, USA) and by capillary electrophoresis (Agilent 2100 Bioanalyser; Agilent Technologies, USA).
Genome sequencing and assembly
The draft genome of A. pullulans var. pullulans was sequenced with Illumina HiSeq using two libraries, fragment, 270-bp-insert size (2× 150 bp reads) and CLIP, 4-kbp-insert size (2× 100 bp reads). Each fastq file was QC filtered for artifacts/ process contamination, and subsequently assembled together with AllPathsLG version R38445 , to produce a 172× coverage 29.6 Mbp assembly in 186 scaffolds (L50 = 1.2 Mbp) and 209 contigs (L50 = 779.8 kbp). The 37121 bp mitochondrion genome was assembled separately in a single contig with AllPathsLG.
The genomes of A. pullulans var. namibiae, A. pullulans var. melanogenum, and A. pullulans var. subglaciale were sequenced using a single fragment Illumina library with 270-bp-insert size (2× 150 bp reads). After filtering for artifacts/ process contamination, the sequence data were assembled with Velvet 1.2.03 . The resulting assembly was used to create a long mate-pair library with an insert of 3000 +/-300 bp, which was then assembled together with the original Illumina library with AllPathsLG release version R42328 . The assembly statistics are summarised in Table 1.
Transcriptome sequencing and assembly
mRNA was purified from total RNA using Absolutely mRNA™ purification kits (Stratagene) and Dynabeads® mRNA purification kits (Invitrogen), and chemically fragmented to 200 bp to 250 bp (Ambion). The mRNA was reverse transcribed with SuperScript II, using random hexamers. Second-strand cDNA was synthesised using a dNTP/dUTP mix (Thermo Scientific), E. coli DNA ligase, E. coli DNA polymerase I, and E coli RnaseH (Invitrogen). The fragmented cDNA was treated with end-pair, A-tailing, and adapter ligation, using TruSeq Sample Preparation kits (Illumina) and KAPA-Illumina library creation kits (KAPA biosystems). Second-strand cDNA was removed by AmpErase UNG (Applied Biosystems), to generate strandedness. qPCR was used to determine the concentrations of the unamplified libraries. The libraries were sequenced on Illumina Hiseq.
In all, 120 million paired-end 100-bp Illumina HiSeq 2000 reads of stranded RNA-seq data were used as input for the de-novo assembly of A. pullulans var. subglaciale EXF-2481. One-hundred and sixty-four million paired-end 150-bp Illumina HiSeq 2000 reads of stranded RNA-seq data were used as input for the de-novo assembly of A. pullulans var. pullulans EXF-150. The reads were assembled into consensus sequences using Rnnotator 2.4.14 , which takes short read sequences as input, and outputs assembled transcript contigs. The reads were trimmed and filtered for quality, low-complexity, adapter and duplicates, and assembled with Velvet 1.2.03 . The minimum contig length was set at 100 bp. The read depth minimum was set to 3 reads. Redundant contigs were removed using Vmatch 2.1.3, and contigs with significant overlap were further assembled using Minimus2, with a minimum overlap of 40 bp. Contig post-processing included splitting misassembled contigs, and contig extension and polishing using the strand information of the reads. Single base errors were corrected by aligning the reads back to each contig with BWA, to generate a consensus nucleotide sequence. Post-processed contigs were clustered into loci and putative transcript precursors were identified. The accuracy, completeness and contiguity of the assembled contigs were checked using Blat alignment  of the contigs to the available genomic reference. Here, 95.87% of the EXF-2481 contigs mapped to the genome, and 90.24% of the EXF-150 contigs mapped to the genome (threshold ≥95% of the contigs).
All four of these genome assemblies of A. pullulans were annotated using the Pipeline JGI annotation , which combines several gene prediction and annotation methods, and integrates the annotated genome into the web-based fungal resource MycoCosm , for comparative genomics.
Before gene prediction, assembly scaffolds were masked using RepeatMasker  and RepBase library , with the most frequent (>150 times) repeats recognised by RepeatScout . The following combination of gene predictors was run on the masked assembly: ab-initio Fgenesh  and GeneMark , trained for specific genomes; homology-based Fgenesh +  and Genewise , seeded by BLASTx alignments  against the NCBI-NR protein database; and transcriptome-based CombEST (Zhou et al., personal communication). In addition to protein-coding genes, tRNAs were predicted using tRNAscan-SE . All of the predicted proteins were functionally annotated using SignalP  for signal sequences, TMHMM  for transmembrane domains, InterProScan  for the integrated collection of functional and structured protein domains, and protein alignments to NCBI-NR, SwissProt (http://www.expasy.org/sprot/), KEGG  for metabolic pathways, and KOG  for eukaryotic clusters of orthologues. Interpro and SwissProt hits were also used to map the gene-ontology terms . For each genomic locus, the best representative gene model was selected based on a combination of protein similarity and EST support.
Shared, unique and duplicated genes
The numbers of shared and unique proteins of all four of the A. pullulans varieties were determined using all-against-all blastp (included in blast 2.2.25+) of their whole proteomes, with the e-value cut-off at 10-10. Additionally, multigene families were predicted with the MCL, to cluster the proteins based on sequence similarity, using blastp alignment scores between proteins as a similarity metric . The numbers of unique and shared clusters were retrieved from MycoCosm .
Mating loci, secreted proteins, membrane transporters
Mating loci were identified by using the JGI online blastp programme (e-value cut-off, 10-5). The S. cerevisiae Mat1-1-1 and Mat1-2-1 proteins were used as sequence queries. Additionally, other proteins characteristic of mating type loci were identified for each of the A. pullulans varieties, as described in . A blastp search was performed for pheromone precursors (α and α-factors), MAPK scaffold (Ste5), transcription factors (Ste12) and pheromone factor receptors (Ste2, Ste3) , and also for the representative set of homologues of meiotic genes from S. cerevisiae (i.e., a ‘meiosis detection toolkit’) .
The prediction of secreted proteins was performed in several steps. First SignalP  was used to search for signal peptides in all of the predicted proteins of each of the four A. pullulans varieties, with the default cut-off D-value of 0.43. Proteins that were predicted by TargetP  to localise in other cellular compartments, as well as those containing predicted transmembrane regions by TMHMM , were removed from the selection. This produced a list of proteins that was expanded with the help of the multigene families predicted with the Markov clustering algorithm (MCL), as described below. Clusters that contained at least half of the proteins predicted as secreted (above) were automatically included in the list of secreted proteins (2506 proteins). This covered the cases where homologues of one or two of the A. pullulans varieties were not included in the original list due to annotation errors, or for other reasons. Clusters that contained less than half of the predicted secreted proteins (2274 proteins) were further analysed with online predictors CELLO  and WegoLoc , the latter with fungal BaCelLo and Höglund databases, with an e-value threshold of 1 × e-10 and a multiplex threshold of 1. The proteins were considered as secreted if they were predicted as such by at least three of four methods (our method described above, and three online predictors). The functions of predicted secreted proteins were assigned manually with the help of the available data about each of these proteins: placement into multigene families, gene ontology terms, presence of PFAM domains, functions of similar proteins from the NCBI-NR protein database, and results of CAZymes Analysis Toolkit  and MEROPS .
The analysis of the major facilitator transporters included all of the A. pullulans predicted protein sequences containing the PFAM domain (PF00083), which is characteristic of sugar transport proteins. The proteins were classified into sugar-specific families based on the phylogenetic analyses (performed as described below), which included the related proteins from other species. These were identified with blastp against the NCBI-NR protein database , or were retrieved from the Transporter Classification Database (http://www.tcdb.org/superfamily.php, 1. 8. 2013), and the MFS family sugar transporters of S. cerevisiae S288C were retrieved from the TransportDB database (http://www.membranetransport.org/, 1. 8. 2013).
Genes for alkali-metal cation transporters from Mycosphaerella graminicola, S. cerevisiae and C. neoformans were recovered from public databases, as described in , while the genes from H. werneckii were published in . Homologues from other Dothideomycetes with sequenced genomes (the same as those used for phylogenetic analyses; see below) were recovered from the Genome Portal of the JGI, and were considered homologous to the analysed proteins from A. pullulans if they formed the same MCL cluster . A second set of trees (for comparison purposes; not shown) was generated by applying a maximum parsimony method as implemented in the Mega software, version 5.05 . The specificity of P-type pumps was inferred by phylogenetic analysis of all of the P-type cation transporters found from A. pullulans and their homologues with known function from other fungi, as described in .
Melanin polyketide synthase genes were identified by searching for a homologue of ALM1 from Alternaria alternata [GenBank:AB665444]. Homologues from Magnaporthe grisea [GenBank:L22309, GenBank:AY846877] were used as queries in the search for the trihydroxynaphthalene reductase-like gene and tetrahydroxynaphthalene reductase-like genes. Protein sequences of orthodox aquaporins Aqy1p and Aqy2p, and of aquaglyceroporins Fps1p and Yfl054p, from S. cerevisiae (sequences from the Saccharomyces genome database; http://www.yeastgenome.org/), were used as queries for the aquaporin genes. Other genes discussed were found with the help of blastp and MCL clustering data on the Genome Portal of the JGI , and analysed with various search strategies against the NCBI-NR protein database and the Saccharomyces Genome database.
For the analysis of phylogenetic relationships of the four A. pullulans varieties and related species, a super alignment of the selected fungal proteomes was constructed with the Hal pipeline , allowing for no missing data. In addition to all four of the A. pullulans varieties, several other proteomes available on the Genome Portal of the Department of Energy Joint Genome Institute  were included: Baudoinia compniacensis, C. fulvum, C. heterostrophus, Hysterium pulicare, M. graminicola, N. crassa, Sclerotinia sclerotiorum, and Tuber melanosporum. A. nidulans was used as an outgroup. Hal produced a conservative alignment of 812692 amino acids. Poorly aligned positions and positions with gaps were removed with Gblocks 0.91b . Stringent parameters were used: the maximum number of contiguous non-conserved positions was limited to 5 amino acids, and the minimum length of a block to 15 amino acids. This produced a 373159-bp-long alignment, which was used for the estimation of phylogeny.
The best protein evolution model was estimated with ProtTest 3.2.1 . The species tree was generated with the PhyML 3.1 software . Approximate Bayes (aBayes) branch supports were calculated. The analysis was run using the LG model of evolution. The ProtTest estimate of the α-parameter of the γ-distribution of six substitution rate categories (0.791), and the determined proportion of invariable sites (0.305) were used. The tree was then calibrated with the r8s software , by assigning the root of the tree to an arbitrary value of 1.
Gene phylogenies (alkali-metal cation transporters, major facilitator superfamily transporters, bacteriorhodopsins, components of HOG, aquaporins, scytalone dehydratases) were estimated by aligning the protein sequences with the MAFFT software in the ‘--auto’ mode . The tree was generated with the PhyML 3.1 software , as described above, by using the model of protein evolution, the α-parameter, and the proportion of invariable sites, as estimated by ProtTest 3.2.1 .
Genome alignment and genome distances
The whole genome alignments between all of the possible pairs of these four genomes were calculated with the promer algorithm, as implemented in Mummer 3.23, and plotted with the mummerplot utility . All of the parameters were the same as described by Hane et al. , except that instead of discarding scaffolds <500 kbp, only scaffolds <200 kbp were discarded.
The distances between these four A. pullulans varieties and between four Saccharomyces species (for comparison purposes) were estimated by using two methods. First, the K r distances were calculated with the genomediff software included in the GenomeTools library . This alignment-free calculation produces Jukes-Cantor corrected divergences between the pairs of genomes as output. As a second approach, the proteomes were aligned with the Hal pipeline . The conservative alignment was used for the estimation of the number of amino-acid substitutions per site, in Mega software version 5.05 . The Poisson correction model was used, and all of the positions containing gaps or missing data were removed prior to analysis, which resulted in a 698915-amino-acid-long alignment.
Analysis of protein family expansions and contractions was performed with the CAFE 3.0 software . Pfam domains of selected fungal proteomes were identified with a stand-alone Pfam scanner and a database downloaded on 30.1.2013 . Additionally, multigene families predicted with the MCL, as described above, were used as a second variant of grouping the predicted proteins. This was used to produce a Table with the numbers of proteins belonging to specific protein families in each of the analysed species. The Table was used as the input of CAFE, together with the chronogram constructed from the whole proteomes, as described above.
This produced a list of groups of proteins with significant expansions/contractions. This was manually checked, and the groups with potential roles in stress tolerance or biotechnological use were selected for further analyses.
The antiSMASH 2.0 software  was used for the identification of the secondary metabolite biosynthesis clusters in the genome sequences. antiSMASH is the first comprehensive pipeline that can identify biosynthetic loci that cover the whole range of known secondary metabolite compound classes (polyketides, non-ribosomal peptides, terpenes, aminoglycosides, aminocoumarins, indolocarbazoles, lantibiotics, bacteriocins, nucleosides, β-lactams, butyrolactones, siderophores, melanins, and others) . Nucleotide sequences of all four of the A. pullulans varieties were uploaded to the antiSMASH online site (http://antismash.secondarymetabolites.org/) in FASTA format, and the parameters were set to default, with the additional setting for the DNA of eukaryotic origin.
The genome assemblies and annotations can be interactively accessed through the JGI fungal genome portal MycoCosm  at http://genome.jgi-psf.org/Aurpu_var_pul1/Aurpu_var_pul1.home.html, http://genome.jgi-psf.org/Aurpu_var_sub1/Aurpu_var_sub1.home.html, http://genome.jgi-psf.org/Aurpu_var_nam1/Aurpu_var_nam1.home.html, and http://genome.jgi-psf.org/Aurpu_var_mel1/Aurpu_var_mel1.home.html, and have also been deposited with DDBJ/EMBL/GenBank under accession [GenBank:AYEO00000000, GenBank:AYYB00000000, GenBank:AYEM00000000, GenBank:AYEN00000000], for A. pullulans, A. subglaciale, A. namibiae, A. melanogenum, respectively.
Antibiotics and secondary metabolites analysis shell
Computational analysis of gene family evolution
High osmolarity glycerol
Major facilitator superfamily
Major intrinsic protein
Mitogen-activated protein kinase
National center for biotechnology information non-redundant
The work conducted by the U.S. Department of Energy Joint Genome Institute is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. The authors acknowledge financial support from the state budget through the Slovenian Research Agency (Infrastructural Centre Mycosmo, MRIC UL, Postdoctoral Project Z4-5531 to CG, and Young Researcher Grant to JZ). The study was also partly financed via operation “Centre of excellence for integrated approaches in chemistry and biology of proteins” number OP188.8.131.52.02.0005, financed by European Regional Development Fund (85% share of financing) and by the Slovenian Ministry of Higher Education, Science and Technology (15% share of financing). The authors would like to thank Chris Berrie for language editing assistance.
- Gostinčar C, Grube M, Gunde-Cimerman N: Evolution of fungal pathogens in domestic environments?. Fungal Biol. 2011, 115: 1008-1018.PubMedGoogle Scholar
- Leathers TD: Biotechnological production and applications of pullulan. Appl Microbiol Biotechnol. 2003, 62: 468-473.PubMedGoogle Scholar
- Cheng KC, Demirci A, Catchmark JM: Pullulan: biosynthesis, production, and applications. Appl Microbiol Biotechnol. 2011, 92: 29-44.PubMedGoogle Scholar
- Tada R, Tanioka A, Ishibashi K, Adachi Y, Tsubaki K, Ohno N: Involvement of branched units at position 6 in the reactivity of a unique variety of beta-D-glucan from Aureobasidium pullulans to antibodies in human sera. Biosci Biotechnol Biochem. 2009, 73: 908-911.PubMedGoogle Scholar
- Muramatsu D, Iwai A, Aoki S, Uchiyama H, Kawata K, Nakayama Y, Nikawa Y, Kusano K, Okabe M, Miyazaki T: β-glucan derived from Aureobasidium pullulans is effective for the prevention of influenza in mice. Plos One. 2012, 7:Google Scholar
- Molnarova J, Vadkertiova R, Stratilova E: Extracellular enzymatic activities and physiological profiles of yeasts colonizing fruit trees. J Basic Microbiol. 2013, 53:Google Scholar
- Buzzini P, Martini A: Extracellular enzymatic activity profiles in yeast and yeast-like strains isolated from tropical environments. J Appl Microbiol. 2002, 93: 1020-1025.PubMedGoogle Scholar
- Chi Z, Wang F, Yue L, Liu G, Zhang T: Bioproducts from Aureobasidium pullulans, a biotechnologically important yeast. Appl Microbiol Biotechnol. 2009, 82: 793-804.PubMedGoogle Scholar
- Rich JO, Leathers TD, Anderson AM, Bischoff KM, Manitchotpisit P: Laccases from Aureobasidium pullulans. Enzyme Microb Technol. 2013, 53: 33-37.PubMedGoogle Scholar
- Takesako K, Ikai K, Haruna F, Endo M, Shimanaka K, Sono E, Nakamura T, Kato I, Yamaguchi H: Aureobasidins, new antifungal antibiotics taxonomy, fermentation, isolation, and properties. J Antibiot (Tokyo). 1991, 44: 919-924.Google Scholar
- Sharma RR, Singh D, Singh R: Biological control of postharvest diseases of fruits and vegetables by microbial antagonists: A review. Biol Control. 2009, 50: 205-221.Google Scholar
- Price NPJ, Manitchotpisit P, Vermillion KE, Bowman MJ, Leathers TD: Structural characterization of novel extracellular liamocins (mannitol oils) produced by Aureobasidium pullulans strain NRRL 50380. Carbohydr Res. 2013, 370: 24-32.PubMedGoogle Scholar
- Andrews JH, Spear RN, Nordheim EV: Population biology of Aureobasidium pullulans on apple leaf surfaces. Can J Microbiol. 2002, 48: 500-513.PubMedGoogle Scholar
- Grube M, Schmid F, Berg G: Black fungi and associated bacterial communities in the phyllosphere of grapevine. Fungal Biol. 2011, 115: 978-986.PubMedGoogle Scholar
- Olstorpe M, Schnurer J, Passoth V: Microbial changes during storage of moist crimped cereal barley grain under Swedish farm conditions. Anim Feed Sci Technol. 2010, 156: 37-46.Google Scholar
- Maciel NOP, Pilo FB, Freitas LFD, Gomes FCO, Johann S, Nardi RMD, Lachance MA, Rosa CA: The diversity and antifungal susceptibility of the yeasts isolated from coconut water and reconstituted fruit juices in Brazil. Int J Food Microbiol. 2013, 160: 201-205.PubMedGoogle Scholar
- Gunde-Cimerman N, Zalar P, de Hoog S, Plemenitaš A: Hypersaline waters in salterns - natural ecological niches for halophilic black yeasts. FEMS Microbiol Ecol. 2000, 32: 235-240.Google Scholar
- Oren A, Gunde-Cimerman N: Fungal life in the Dead Sea. Biology of Marine Fungi. Edited by: Raghukumar C. 2012Google Scholar
- Zalar P, Gostinčar C, de Hoog GS, Uršič V, Sudhadham M, Gunde-Cimerman N: Redefinition of Aureobasidium pullulans and its varieties. Stud Mycol. 2008, 61: 21-38.PubMed CentralPubMedGoogle Scholar
- Branda E, Turchetti B, Diolaiuti G, Pecci M, Smiraglia C, Buzzini P: Yeast and yeast-like diversity in the southernmost glacier of Europe (Calderone Glacier, Apennines, Italy). FEMS Microbiol Ecol. 2010, 72: 354-369.PubMedGoogle Scholar
- Vaz ABM, Rosa LH, Vieira MLA, de Garcia V, Brandao LR, Teixeira LCRS, Moline M, Libkind D, van Broock M, Rosa CA: The diversity, extracellular enzymatic activities and photoprotective compounds of yeasts isolated in Antarctica. Braz J Microbiol. 2011, 42: 937-947.PubMed CentralPubMedGoogle Scholar
- D'Elia T, Veerapaneni R, Theraisnathan V, Rogers SO: Isolation of fungi from Lake Vostok accretion ice. Mycologia. 2009, 101: 751-763.PubMedGoogle Scholar
- Vadkertiova R, Slavikova E: Killer activity of yeasts isolated from the water environment. Can J Microbiol. 1995, 41: 759-766.PubMedGoogle Scholar
- Pitt JI, Hocking AD: Fungi and food spoilage. 1999, Gaithersburg, Maryland: Aspen Publishers, Inc., 2Google Scholar
- Nisiotou AA, Chorianopoulos N, Nychas GJE, Panagou EZ: Yeast heterogeneity during spontaneous fermentation of black Conservolea olives in different brine solutions. J Appl Microbiol. 2010, 108: 396-405.PubMedGoogle Scholar
- Lotrakul P, Deenarn P, Prasongsuk S, Punnapayak H: Isolation of Aureobasidium pullulans from bathroom surfaces and their antifungal activity against some Aspergilli. Afr J Microbiol Res. 2009, 3: 253-257.Google Scholar
- Kaarakainen P, Rintala H, Vepsalainen A, Hyvarinen A, Nevalainen A, Meklin T: Microbial content of house dust samples determined with qPCR. Sci Total Environ. 2009, 407: 4673-4680.PubMedGoogle Scholar
- Zalar P, Novak M, De Hoog GS, Gunde-Cimerman N: Dishwashers - A man-made ecological niche accommodating human opportunistic fungal pathogens. Fungal Biol. 2011, 115: 997-1007.PubMedGoogle Scholar
- Arvanitidou M, Kanellou K, Constantinides TC, Katsouyannopoulos V: The occurrence of fungi in hospital and community potable waters. Lett Apl Microbiol. 1999, 29: 81-84.Google Scholar
- Zhang ES, Tanaka T, Tajima M, Tsuboi R, Nishikawa A, Sugita T: Characterization of the skin fungal microbiota in patients with atopic dermatitis and in healthy subjects. Microbiol Immunol. 2011, 55: 625-632.PubMedGoogle Scholar
- Rauch ME, Graef HW, Rozenzhak SM, Jones SE, Bleckmann CA, Kruger RL, Naik RR, Stone MO: Characterization of microbial contamination in United States Air Force aviation fuel tanks. J Ind Microbiol Biotechnol. 2006, 33: 29-36.PubMedGoogle Scholar
- Cappitelli F, Sorlini C: Microorganisms attack synthetic polymers in items representing our cultural heritage. Appl Environ Microbiol. 2008, 74: 564-569.PubMed CentralPubMedGoogle Scholar
- Shah AA, Hasan F, Hameed A, Ahmed S: Biological degradation of plastics: a comprehensive review. Biotechnol Adv. 2008, 26: 246-265.PubMedGoogle Scholar
- Hawkes M, Rennie R, Sand C, Vaudry W: Aureobasidium pullulans infection: fungemia in an infant and a review of human cases. Diagn Microbiol Infect Dis. 2005, 51: 209-213.PubMedGoogle Scholar
- Chan GF, Puad MSA, Chin CF, Rashid NAA: Emergence of Aureobasidium pullulans as human fungal pathogen and molecular assay for future medical diagnosis. Folia Microbiol (Praha). 2011, 56: 459-467.Google Scholar
- Chan GF, Bamadhaj HM, Gan HM, Rashid NAA: Genome sequence of Aureobasidium pullulans AY4, an emerging opportunistic fungal pathogen with diverse biotechnological potential. Eukaryot Cell. 2012, 11: 1419-1420.PubMed CentralPubMedGoogle Scholar
- Gostinčar C, Grube M, de Hoog GS, Zalar P, Gunde-Cimerman N: Extremotolerance in fungi: evolution on the edge. FEMS Microbiol Ecol. 2010, 71: 2-11.PubMedGoogle Scholar
- Ranta HM: Effect of simulated acid rain on quantity of epiphytic microfungi on Scots pine (Pinus sylvestris L.) needles. Environ Pollut. 1990, 67: 349-359.PubMedGoogle Scholar
- Shiomi N, Yasuda T, Inoue Y, Kusumoto N, Iwasaki S, Katsuda T, Katoh S: Characteristics of neutralization of acids by newly isolated fungal cells. J Biosci Bioeng. 2004, 97: 54-58.PubMedGoogle Scholar
- Onofri S: Antarctic Microfungi. Enigmatic microorganisms and life in extreme environments. Edited by: Seckbach J. 1999, Dordrecht; London: Kluwer Academic, 323-336. Cellular origin and life in extreme habitats; v.1,Google Scholar
- Kogej T, Ramos J, Plemenitas A, Gunde-Cimerman N: The halophilic fungus Hortaea werneckii and the halotolerant fungus Aureobasidium pullulans maintain low intracellular cation concentrations in hypersaline environments. Appl Environ Microbiol. 2005, 71: 6600-6605.PubMed CentralPubMedGoogle Scholar
- Managbanag JR, Torzilli AP: An analysis of trehalose, glycerol, and mannitol accumulation during heat and salt stress in a salt marsh isolate of Aureobasidium pullulans. Mycologia. 2002, 94: 384-391.PubMedGoogle Scholar
- Kogej T, Gostinčar C, Volkmann M, Gorbushina AA, Gunde-Cimerman N: Mycosporines in extremophilic fungi — novel complementary osmolytes?. Environ Chem. 2006, 3: 105-110.Google Scholar
- Turk M, Mejanelle L, Šentjurc M, Grimalt JO, Gunde-Cimerman N, Plemenitaš A: Salt-induced changes in lipid composition and membrane fluidity of halophilic yeast-like melanized fungi. Extremophiles. 2004, 8: 53-61.PubMedGoogle Scholar
- Gostinčar C, Turk M, Trbuha T, Vaupotič T, Plemenitaš A, Gunde-Cimerman N: Expression of fatty-acid-modifying enzymes in the halotolerant black yeast Aureobasidium pullulans (de Bary) G. Arnaud under salt stress. Stud Mycol. 2008, 61: 51-59.PubMed CentralPubMedGoogle Scholar
- Bermejo JM, Dominguez JB, Goni FM, Uruburu F: Influence of pH on the transition from yeast-like cells to chlamydospores in Aureobasidium pullulans. Antonie Van Leeuwenhoek. 1981, 47: 385-392.PubMedGoogle Scholar
- Crous PW, Summerell BA, Swart L, Denman S, Taylor JE, Bezuidenhout CM, Palm ME, Marincowitz S, Groenewald JZ: Fungal pathogens of Proteaceae. Persoonia. 2011, 27: 20-45.PubMed CentralPubMedGoogle Scholar
- Peterson SW, Manitchotpisit P, Leathers TD: Aureobasidium thailandense sp. nov. isolated from leaves and wooden surfaces. Int J Syst Evol Microbiol. 2013, 63: 790-795.PubMedGoogle Scholar
- Manitchotpisit P, Leathers TD, Peterson SW, Kurtzman CP, Li XL, Eveleigh DE, Lotrakul P, Prasongsuk S, Dunlap CA, Vermillion KE, Punnapayak H: Multilocus phylogenetic analyses, pullulan production and xylanase activity of tropical isolates of Aureobasidium pullulans. Mycol Res. 2009, 113: 1107-1120.PubMedGoogle Scholar
- Slepecky RA, Starmer WT: Phenotypic plasticity in fungi: a review with observations on Aureobasidium pullulans. Mycologia. 2009, 101: 823-832.PubMedGoogle Scholar
- Ohm RA, Feau N, Henrissat B, Schoch CL, Horwitz BA, Barry KW, Condon BJ, Copeland AC, Dhillon B, Glaser F, Hesse CN, Kosti I, Labutti K, Lindquist EA, Lucas S, Salamov AA, Bradshaw RE, Ciuffetti L, Hamelin RC, Kema GH, Lawrence C, Scott JA, Spatafora JW, Turgeon BG, de Wit PJ, Zhong S, Goodwin SB, Grigoriev IV: Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog. 2012, 8: e1003037-PubMed CentralPubMedGoogle Scholar
- Gerstein AC, Chun HJE, Grant A, Otto SP: Genomic convergence toward diploidy in Saccharomyces cerevisiae. PLoS Genet. 2006, 2: 1396-1401.Google Scholar
- Dhar R, Sagesser R, Weikert C, Yuan J, Wagner A: Adaptation of Saccharomyces cerevisiae to saline stress through laboratory evolution. J Evol Biol. 2011, 24: 1135-1153.PubMedGoogle Scholar
- Lenassi M, Gostinčar C, Jackman S, Turk M, Sadowski I, Nislow C, Jones S, Birol I, Gunde-Cimerman N, Plemenitaš A: Whole genome duplication and enrichment of metal cation transporters revealed by de novo genome sequencing of extremely halotolerant black yeast Hortaea werneckii. PLoS ONE. 2013, 8: e71328-PubMed CentralPubMedGoogle Scholar
- Duan XH, Chi ZM, Wang L, Wang XH: Influence of different sugars on pullulan production and activities of alpha-phosphoglucose mutase, UDPG-pyrophosphorylase and glucosyltransferase involved in pullulan synthesis in Aureobasidium pullulans Y68. Carbohydr Polym. 2008, 73: 587-593.PubMedGoogle Scholar
- Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B: The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014, 42: D490-D495.PubMed CentralPubMedGoogle Scholar
- Kang BK, Yang HJ, Choi NS, Ahn KH, Park CS, Yoon BD, Kim MS: Production of pure beta-glucan by Aureobasidium pullulans after pullulan synthetase gene disruption. Biotechnol Lett. 2010, 32: 137-142.PubMedGoogle Scholar
- Zhang DP, Spadaro D, Valente S, Garibaldi A, Gullino ML: Cloning, characterization, expression and antifungal activity of an alkaline serine protease of Aureobasidium pullulans PL5 involved in the biological control of postharvest pathogens. Int J Food Microbiol. 2012, 153: 453-464.PubMedGoogle Scholar
- Chi Z, Wang XX, Ma ZC, Buzdar MA, Chi ZM: The unique role of siderophore in marine-derived Aureobasidium pullulans HN6.2. Biometals. 2012, 25: 219-230.PubMedGoogle Scholar
- Takesako K, Kuroda H, Inoue T, Haruna F, Yoshikawa Y, Kato I, Uchida K, Hiratani T, Yamaguchi H: Biological properties of aureobasidin A, a cyclic depsipeptide antifungal antibiotic. J Antibiot (Tokyo). 1993, 46: 1414-1420.Google Scholar
- Slightom JL, Metzger BP, Luu HT, Elhammer AP: Cloning and molecular characterization of the gene encoding the aureobasidin A biosynthesis complex in Aureobasidium pullulans BP-1938. Gene. 2009, 431: 67-79.PubMedGoogle Scholar
- Wang WL, Chi ZM, Chi Z, Li J, Wang XH: Siderophore production by the marine-derived Aureobasidium pullulans and its antimicrobial activity. Bioresour Technol. 2009, 100: 2639-2641.PubMedGoogle Scholar
- Johnson L: Iron and siderophores in fungal-host interactions. Mycol Res. 2008, 112: 170-183.PubMedGoogle Scholar
- Brakhage AA: Regulation of fungal secondary metabolism. Nat Rev Microbiol. 2013, 11: 21-32.PubMedGoogle Scholar
- Scharf DH, Brakhage AA: Engineering fungal secondary metabolism: A roadmap to novel compounds. J Biotechnol. 2013, 163: 179-183.PubMedGoogle Scholar
- Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, Weber T, Takano E, Breitling R: antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011, 39: 339-346.Google Scholar
- Zhong XY, Tian YQ, Niu GQ, Tan HR: Assembly and features of secondary metabolite biosynthetic gene clusters in Streptomyces ansochromogenes. Sci China Life Sci. 2013, 56: 609-618.PubMedGoogle Scholar
- Knerr PJ, van der Donk WA: Discovery, biosynthesis, and engineering of lantipeptides. Annu Rev Biochem. 2012, 81: 479-505.PubMedGoogle Scholar
- Choi J, Park J, Kim D, Jung K, Kang S, Lee YH: Fungal Secretome Database: Integrated platform for annotation of fungal secretomes. BMC Genomics. 2010, 11:Google Scholar
- Zhao ZT, Liu HQ, Wang CF, Xu JR: Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genomics. 2013, 14:Google Scholar
- Han MV, Thomas GWC, Lugo-Martinez J, Hahn MW: Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol Biol Evol. 2013, 30: 1987-1997.PubMedGoogle Scholar
- Rawlings ND, Barrett AJ, Bateman A: MEROPS: the database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2012, 40: D343-D350.PubMed CentralPubMedGoogle Scholar
- Belmares R, Contreras-Esquivel JC, Rodriguez-Herrera R, Coronel AR, Aguilar CN: Microbial production of tannase: an enzyme with potential use in food industry. Lebensm Wiss Technol. 2004, 37: 857-864.Google Scholar
- Isola D, Selbmann L, de Hoog GS, Fenice M, Onofri S, Prenafeta-Boldu FX, Zucconi L: Isolation and screening of black fungi as degraders of volatile aromatic hydrocarbons. Mycopathologia. 2013, 175: 369-379.PubMedGoogle Scholar
- Prenafeta-Boldu FX, Guivernau M, Gallastegui G, Vinas M, de Hoog GS, Elias A: Fungal/bacterial interactions during the biodegradation of TEX hydrocarbons (toluene, ethylbenzene and p-xylene) in gas biofilters operated under xerophilic conditions. FEMS Microbiol Ecol. 2012, 80: 722-734.PubMedGoogle Scholar
- Webb JS, Nixon M, Eastwood IM, Greenhalgh M, Robson GD, Handley PS: Fungal colonization and biodeterioration of plasticized polyvinyl chloride. Appl Environ Microbiol. 2000, 66: 3194-3200.PubMed CentralPubMedGoogle Scholar
- Howard GT: Biodegradation of polyurethane: a review. Int Biodeterior Biodegradation. 2002, 49: 245-252.Google Scholar
- Shivakumar S: Poly-beta-hydroxybutyrate (PHB) Depolymerase from Fusarium solani Thom. J Chem. 2013Google Scholar
- Reddy VS, Shlykov MA, Castillo R, Sun EI, Saier MH: The major facilitator superfamily (MFS) revisited. FEBS J. 2012, 279: 2022-2035.PubMed CentralPubMedGoogle Scholar
- Saier MH, Tran CV, Barabote RD: TCDB: the Transporter Classification Database for membrane transport protein analyses and information. Nucleic Acids Res. 2006, 34: D181-D186.PubMed CentralPubMedGoogle Scholar
- Saier MH: Families of transmembrane sugar transport proteins. Mol Microbiol. 2000, 35: 699-710.PubMedGoogle Scholar
- Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, Rehman B, Elkins T, Engels R, Wang S, Nielsen CB, Butler J, Endrizzi M, Qui D, Ianakiev P, Bell-Pedersen D, Nelson MA, Werner-Washburne M, Selitrennikoff CP, Kinsey JA, Braun EL, Zelter A, Schulte U, Kothe GO, Jedd G, Mewes W, et al: The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003, 422: 859-868.PubMedGoogle Scholar
- Leandro MJ, Fonseca C, Goncalves P: Hexose and pentose transport in ascomycetous yeasts: an overview. FEMS Yeast Res. 2009, 9: 511-525.PubMedGoogle Scholar
- Hohmann S: Control of high osmolarity signalling in the yeast Saccharomyces cerevisiae. FEBS Lett. 2009, 583: 4025-4029.PubMedGoogle Scholar
- Bahn YS: Master and commander in fungal pathogens: the two-component system and the HOG signaling pathway. Eukaryot Cell. 2008, 7: 2017-2036.PubMed CentralPubMedGoogle Scholar
- Krantz M, Becit E, Hohmann S: Comparative genomics of the HOG-signalling system in fungi. Curr Genet. 2006, 49: 137-151.PubMedGoogle Scholar
- Musacchio A, Gibson T, Rice P, Thompson J, Saraste M: The PH Domain - a common piece in the structural patchwork of signaling proteins. Trends Biochem Sci. 1993, 18: 343-348.PubMedGoogle Scholar
- Gostinčar C, Muggia L, Grube M: Polyextremotolerant black fungi: oligotrophism, adaptive potential, and a link to lichen symbioses. Front Microbiol. 2012, 3: 390-PubMed CentralPubMedGoogle Scholar
- Robinson CH: Cold adaptation in Arctic and Antarctic fungi. New Phytol. 2001, 151: 341-353.Google Scholar
- Ansell R, Granath K, Hohmann S, Thevelein JM, Adler L: The two isoenzymes for yeast NAD(+)-dependent glycerol-3-phosphate dehydrogenase encoded by GPD1 and GPD2 have distinct roles in osmoadaptation and redox regulation. EMBO J. 1997, 16: 2179-2187.PubMed CentralPubMedGoogle Scholar
- Zajc J, Liu Y, Dai W, Yang Z, Hu J, Gostinčar C, Gunde-Cimerman N: Genome and transcriptome sequencing of the halophilic fungus Wallemia ichthyophaga: haloadaptations present and absent. BMC Genomics. 2013, 14: 617-PubMed CentralPubMedGoogle Scholar
- Lenassi M, Zajc J, Gostinčar C, Gorjan A, Gunde-Cimerman N, Plemenitaš A: Adaptation of the glycerol-3-phosphate dehydrogenase Gpd1 to high salinities in the extremely halotolerant Hortaea werneckii and halophilic Wallemia ichthyophaga. Fungal Biol. 2011, 115: 959-970.PubMedGoogle Scholar
- Jung S, Marelli M, Rachubinski RA, Goodlett DR, Aitchison JD: Dynamic changes in the subcellular distribution of Gpd1p in response to cell stress. J Biol Chem. 2010, 285: 6739-6749.PubMed CentralPubMedGoogle Scholar
- Francois J, Parrou JL: Reserve carbohydrates metabolism in the yeast Saccharomyces cerevisiae. FEMS Microbiol Rev. 2001, 25: 125-145.PubMedGoogle Scholar
- Ruijter GJG, Bax M, Patel H, Flitter SJ, van de Vondervoort PJI, de Vries RP, van Kuyk PA, Visser J: Mannitol is required for stress tolerance in Aspergillus niger conidiospores. Eukaryot Cell. 2003, 2: 690-698.PubMed CentralPubMedGoogle Scholar
- Silva SS, Felipe MGA, Mancilha IM: Factors that affect the biosynthesis of xylitol by xylose-fermenting yeasts - A review. Appl Biochem Biotechnol. 1998, 70–72: 331-339.PubMedGoogle Scholar
- Granstrom TB, Izumori K, Leisola M: A rare sugar xylitol. Part I: the biochemistry and biosynthesis of xylitol. Appl Microbiol Biotechnol. 2007, 74: 277-281.PubMedGoogle Scholar
- Sugai JK, Veiga LA: Purification and properties of the xylitol dehydrogenase from Pullularia pullulans. An Acad Bras Cienc. 1981, 53: 183-193.PubMedGoogle Scholar
- Siehr DJ: Melanin biosynthesis in Aureobasidium pullulans. J Coating Tech. 1981, 53: 23-25.Google Scholar
- Bell AA, Wheeler MH: Biosynthesis and functions of fungal melanins. Annu Rev Phytopathol. 1986, 24: 411-451.Google Scholar
- Kogej T, Stein M, Volkmann M, Gorbushina AA, Galinski EA, Gunde-Cimerman N: Osmotic adaptation of the halophilic fungus Hortaea werneckii: role of osmolytes and melanization. Microbiol. 2007, 153: 4261-4273.Google Scholar
- Butler MJ, Day AW: Fungal melanins: a review. Can J Microbiol. 1998, 44: 1115-1136.Google Scholar
- Taborda CP, da Silva MB, Nosanchuk JD, Travassos LR: Melanin as a virulence factor of Paracoccidioides brasiliensis and other dimorphic pathogenic fungi: a minireview. Mycopathologia. 2008, 165: 331-339.PubMed CentralPubMedGoogle Scholar
- Yurlova NA, de Hoog GS, Fedorova LG: The influence of ortho- and para-diphenoloxidase substrates on pigment formation in black yeast-like fungi. Stud Mycol. 2008, 39-49.Google Scholar
- Kogej T, Wheeler MH, Lanisnik Rizner T, Gunde-Cimerman N: Evidence for 1,8-dihydroxynaphthalene melanin in three halophilic black yeasts grown under saline and non-saline conditions. FEMS Microbiol Lett. 2004, 232: 203-209.PubMedGoogle Scholar
- Langfelder K, Streibel M, Jahn B, Haase G, Brakhage AA: Biosynthesis of fungal melanins and their importance for human pathogenic fungi. Fungal Genet Biol. 2003, 38: 143-158.PubMedGoogle Scholar
- Takano Y, Kubo Y, Kawamura C, Tsuge T, Furusawa I: The Alternaria alternata melanin biosynthesis gene restores appressorial melanization and penetration of cellulose membranes in the melanin-deficient albino mutant of Colletotrichum lagenarium. Fungal Genet Biol. 1997, 21: 131-140.Google Scholar
- Takano Y, Kubo Y, Shimizu K, Mise K, Okuno T, Furusawa I: Structural analysis of Pks1, a polyketide synthase gene involved in melanin biosynthesis in Colletotrichum lagenarium. Mol Gen Genet. 1995, 249: 162-167.PubMedGoogle Scholar
- Zhang A, Lu P, Dahl-Roshak AM, Paress PS, Kennedy S, Tkacz JS, An Z: Efficient disruption of a polyketide synthase gene (pks1) required for melanin synthesis through Agrobacterium-mediated transformation of Glarea lozoyeasis. Mol Genet Genomics. 2003, 268: 645-655.PubMedGoogle Scholar
- Tsai HF, Wheeler MH, Chang YC, Kwon-Chung KJ: A developmentally regulated gene cluster involved in conidial pigment biosynthesis in Aspergillus fumigatus. J Bacteriol. 1999, 181: 6469-6477.PubMed CentralPubMedGoogle Scholar
- Borgnia MJ, Agre P: Reconstitution and functional comparison of purified GlpF and AqpZ, the glycerol and water channels from Escherichia coli. Proc Natl Acad Sci USA. 2001, 98: 2888-2893.PubMed CentralPubMedGoogle Scholar
- Alleva K, Chara O, Amodeo G: Aquaporins: Another piece in the osmotic puzzle. FEBS Lett. 2012, 586: 2991-2999.PubMedGoogle Scholar
- Borgnia M, Nielsen S, Engel A, Agre P: Cellular and molecular biology of the aquaporin water channels. Annu Rev Biochem. 1999, 68: 425-458.PubMedGoogle Scholar
- Hohmann S, Bill RM, Kayingo G, Prior BA: Microbial MIP channels. Trends Microbiol. 2000, 8: 33-38.PubMedGoogle Scholar
- Pettersson N, Filipsson C, Becit E, Brive L, Hohmann S: Aquaporins in yeasts and filamentous fungi. Biol Cell. 2005, 97: 487-500.PubMedGoogle Scholar
- Xu H, Cooke JEK, Zwiazek JJ: Phylogenetic analysis of fungal aquaporins provides insight into their possible role in water transport of mycorrhizal associations. Botany-Botanique. 2013, 91: 495-504.Google Scholar
- Arino J, Ramos J, Sychrova H: Alkali metal cation transport and homeostasis in yeasts. Microbiol Mol Biol Rev. 2010, 74: 95-120.PubMed CentralPubMedGoogle Scholar
- Shabala S, Cuin TA: Potassium transport and plant salt tolerance. Physiol Plant. 2008, 133: 651-669.PubMedGoogle Scholar
- Ko CH, Gaber RF: Trk1 and Trk2 encode structurally related K+ transporters in Saccharomyces cerevisiae. Mol Cell Biol. 1991, 11: 4266-4273.PubMed CentralPubMedGoogle Scholar
- Ketchum KA, Joiner WJ, Sellers AJ, Kaczmarek LK, Goldstein SA: A new family of outwardly rectifying potassium channel proteins with two pore domains in tandem. Nature. 1995, 376: 690-695.PubMedGoogle Scholar
- Haro R, Garciadeblas B, Rodriguez-Navarro A: A novel P-type ATPase from yeast involved in sodium-transport. FEBS Lett. 1991, 291: 189-191.PubMedGoogle Scholar
- Garciadeblas B, Rubio F, Quintero FJ, Banuelos MA, Haro R, Rodriguez-Navarro A: Differential expression of two genes encoding isoforms of the ATPase involved in sodium efflux in Saccharomyces cerevisiae. Mol Gen Genet. 1993, 236: 363-368.PubMedGoogle Scholar
- Wieland J, Nitsche AM, Strayle J, Steiner H, Rudolph HK: The PMR2 gene-cluster encodes functionally distinct isoforms of a putative Na+ pump in the yeast plasma-membrane. EMBO J. 1995, 14: 3870-3882.PubMed CentralPubMedGoogle Scholar
- Ramos J, Arino J, Sychrova H: Alkali-metal-cation influx and efflux systems in nonconventional yeast species. FEMS Microbiol Lett. 2011, 317: 1-8.PubMedGoogle Scholar
- Serrano R, Kiellandbrandt MC, Fink GR: Yeast plasma-membrane ATPase is essential for growth and has homology with (Na++K+), K+- and Ca2+-ATPases. Nature. 1986, 319: 689-693.PubMedGoogle Scholar
- Ambesi A, Miranda M, Petrov VV, Slayman CW: Biogenesis and function of the yeast plasma-membrane H+-ATPase. J Exp Biol. 2000, 203: 155-160.PubMedGoogle Scholar
- Waschuk SA, Bezerra AG, Shi L, Brown LS: Leptosphaeria rhodopsin: Bacteriorhodopsin-like proton pump from a eukaryote. Proc Natl Acad Sci USA. 2005, 102: 6879-6883.PubMed CentralPubMedGoogle Scholar
- Gostinčar C, Lenassi M, Gunde-Cimerman N, Plemenitaš A: Fungal adaptation to extremely high salt concentrations. Adv Appl Microbiol. 2011, 77: 71-96.PubMedGoogle Scholar
- Wagner NL, Greco JA, Ranaghan MJ, Birge RR: Directed evolution of bacteriorhodopsin for applications in bioelectronics. J R Soc Interface. 2013, 10:Google Scholar
- Zhang F, Vierock J, Yizhar O, Fenno LE, Tsunoda S, Kianianmomeni A, Prigge M, Berndt A, Cushman J, Polle J, Magnuson J, Hegemann P, Deisseroth K: The microbial opsin family of optogenetic tools. Cell. 2011, 147: 1446-1457.PubMed CentralPubMedGoogle Scholar
- Schoch CL, Seifert KA, Huhndorf S, Robert V, Spouge JL, Levesque CA, Chen W, Bolchacova E, Voigt K, Crous PW, Miller AN, Wingfield MJ, Aime MC, An KD, Bai FY, Barreto RW, Begerow D, Bergeron MJ, Blackwell M, Boekhout T, Bogale M, Boonyuen N, Burgaz AR, Buyck B, Cai L, Cai Q, Cardinali G, Chaverri P, Coppins BJ, Crespo A, et al: Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc Natl Acad Sci USA. 2012, 109: 6241-6246.PubMed CentralPubMedGoogle Scholar
- Manitchotpisit P, Skory CD, Peterson SW, Price NPJ, Vermillion KE, Leathers TD: Poly(beta-L-malic acid) production by diverse phylogenetic clades of Aureobasidium pullulans. J Ind Microbiol Biotechnol. 2012, 39: 125-132.PubMedGoogle Scholar
- Spatafora JW, Sung GH, Johnson D, Hesse C, O'Rourke B, Serdani M, Spotts R, Lutzoni F, Hofstetter V, Miadlikowska J, Reeb V, Gueidan C, Fraker E, Lumbsch T, Lucking R, Schmitt I, Hosaka K, Aptroot A, Roux C, Miller AN, Geiser DM, Hafellner J, Hestmark G, Arnold AE, Budel B, Rauhut A, Hewitt D, Untereiner WA, Cole MS, Scheidegger C, et al: A five-gene phylogeny of Pezizomycotina. Mycologia. 2006, 98: 1018-1028.PubMedGoogle Scholar
- Liu YJ, Hall BD: Body plan evolution of ascomycetes, as inferred from an RNA polymerase II, phylogeny. Proc Natl Acad Sci USA. 2004, 101: 4507-4512.PubMed CentralPubMedGoogle Scholar
- Schmitt I, Crespo A, Divakar PK, Fankhauser JD, Herman-Sackett E, Kalb K, Nelsen MP, Nelson NA, Rivas-Plata E, Shimp AD, Widhelm T, Lumbsch HT: New primers for promising single-copy genes in fungal phylogenetics and systematics. Persoonia. 2009, 23: 35-40.PubMed CentralPubMedGoogle Scholar
- de Hoog GS, Guarro J, Gene J, Figueras MJ: Atlas of clinical fungi. 2000, Utrecht: ASM Press, 2Google Scholar
- Debuchy R, Turgeon BG: Mating-type structure, evolution, and function in Euascomycetes. The Mycota: A Comprehensive Treatise on Fungi as Experimental Systems for Basic and Applied Research: Growth, Differentiation and Sexuality, Volume 1. Edited by: Kües U, Fischer R. 2006, Berlin, Heidelberg: Springer, 293-323.Google Scholar
- Ni M, Feretzaki M, Sun S, Wang XY, Heitman J: Sex in Fungi. Annu Rev Genet. 2011, 45: 405-430.PubMed CentralPubMedGoogle Scholar
- Butler G, Kenny C, Fagan A, Kurischko C, Gaillardin C, Wolfe KH: Evolution of the MAT locus and its Ho endonuclease in yeast species. Proc Natl Acad Sci USA. 2004, 101: 1632-1637.PubMed CentralPubMedGoogle Scholar
- Malik SB, Pightling AW, Stefaniak LM, Schurko AM, Logsdon JM: An expanded inventory of conserved meiotic genes provides evidence for sex in Trichomonas vaginalis. Plos One. 2008, 3:Google Scholar
- Nakagawa T, Kolodner RD: Saccharomyces cerevisiae Mer3 is a DNA helicase involved in meiotic crossing over. Mol Cell Biol. 2002, 22: 3281-3291.PubMed CentralPubMedGoogle Scholar
- Gioti A, Nystedt B, Li WJ, Xu J, Andersson A, Averette AF, Munch K, Wang XY, Kappauf C, Kingsbury JM, Kraak B, Walker LA, Johansson HJ, Holm T, Lehtio J, Stajich JE, Mieczkowski P, Kahmann R, Kennell JC, Cardenas ME, Lundeberg J, Saunders CW, Boekhout T, Dawson TL, Munro CA, de Groot PWJ, Butler G, Heitman J, Scheynius A: Genomic insights into the atopic eczema-associated skin commensal yeast Malassezia sympodialis. Mbio. 2013, 4:Google Scholar
- Srikantha T, Daniels KJ, Pujol C, Sahni N, Yi S, Soll DR: Nonsex genes in the mating type locus of Candida albicans play roles in a/alpha biofilm formation, including impermeability and fluconazole resistance. PLoS Pathog. 2012, 8:Google Scholar
- Rozman D, Komel R: Isolation of genomic DNA from filamentous fungi with high glucan level. BioTechniques. 1994, 16: 382-383.PubMedGoogle Scholar
- Gnerre S, MacCallum I, Przybylski D, Ribeiro FJ, Burton JN, Walker BJ, Sharpe T, Hall G, Shea TP, Sykes S, Berlin AM, Aird D, Costello M, Daza R, Williams L, Nicol R, Gnirke A, Nusbaum C, Lander ES, Jaffe DB: High-quality draft assemblies of mammalian genomes from massively parallel sequence data. Proc Natl Acad Sci USA. 2011, 108: 1513-1518.PubMed CentralPubMedGoogle Scholar
- Zerbino DR, Birney E: Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829.PubMed CentralPubMedGoogle Scholar
- Martin J, Bruno VM, Fang ZD, Meng XD, Blow M, Zhang T, Sherlock G, Snyder M, Wang Z: Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads. BMC Genomics. 2010, 11:Google Scholar
- Kent WJ: BLAT - The BLAST-like alignment tool. Genome Res. 2002, 12: 656-664.PubMed CentralPubMedGoogle Scholar
- Grigoriev I, Martines D, Salamov A: Fungal genomic annotation. Applied mycology and biotechnology, Volume 6. Edited by: Aurora D, Berka R, Singh G. 2006, Amsterdam London: Elsevier, 123-142.Google Scholar
- Grigoriev IV, Nikitin R, Haridas S, Kuo A, Ohm R, Otillar R, Riley R, Salamov A, Zhao X, Korzeniewski F, Smirnova T, Nordberg H, Dubchak I, Shabalov I: MycoCosm portal: gearing up for 1000 fungal genomes. Nucleic Acids Res. 2014, 42: D699-D704.PubMed CentralPubMedGoogle Scholar
- RepeatMasker Open-3.0. http://www.repeatmasker.org,
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467.PubMedGoogle Scholar
- Price AL, Jones NC, Pevzner PA: De novo identification of repeat families in large genomes. Bioinformatics. 2005, 21: I351-I358.PubMedGoogle Scholar
- Salamov AA, Solovyev VV: Ab initio gene finding in Drosophila genomic DNA. Genome Res. 2000, 10: 516-522.PubMed CentralPubMedGoogle Scholar
- Ter-Hovhannisyan V, Lomsadze A, Chernoff YO, Borodovsky M: Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. Genome Res. 2008, 18: 1979-1990.PubMed CentralPubMedGoogle Scholar
- Birney E, Clamp M, Durbin R: GeneWise and Genomewise. Genome Res. 2004, 14: 988-995.PubMed CentralPubMedGoogle Scholar
- Altschul SF, Madden TL, Shaffer AA, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402.PubMed CentralPubMedGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.PubMed CentralPubMedGoogle Scholar
- Nielsen H, Engelbrecht J, Brunak S, von Heijne G: Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Eng. 1997, 10: 1-6.PubMedGoogle Scholar
- Melen K, Krogh A, von Heijne G: Reliability measures for membrane protein topology prediction algorithms. J Mol Biol. 2003, 327: 735-744.PubMedGoogle Scholar
- Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: Protein domains identifier. Nucleic Acids Res. 2005, 33: W116-W120.PubMed CentralPubMedGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36: D480-D484.PubMed CentralPubMedGoogle Scholar
- Koonin EV, Fedorova ND, Jackson JD, Jacobs AR, Krylov DM, Makarova KS, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Rogozin IB, Smirnov S, Sorokin AV, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA: A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes. Genome Biol. 2004, 5:Google Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G, Consortium GO: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29.PubMed CentralPubMedGoogle Scholar
- Enright AJ, Van Dongen S, Ouzounis CA: An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 2002, 30: 1575-1584.PubMed CentralPubMedGoogle Scholar
- Grigoriev IV, Nordberg H, Shabalov I, Aerts A, Cantor M, Goodstein D, Kuo A, Minovitsky S, Nikitin R, Ohm RA, Otillar R, Poliakov A, Ratnere I, Riley R, Smirnova T, Rokhsar D, Dubchak I: The genome portal of the Department of Energy Joint Genome Institute. Nucleic Acids Res. 2012, 40: 26-32.Google Scholar
- Merlini L, Dudin O, Martin SG: Mate and fuse: how yeast cells do it. Open Biology. 2013, 3:Google Scholar
- Petersen TN, Brunak S, von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8: 785-786.PubMedGoogle Scholar
- Emanuelsson O, Nielsen H, Brunak S, von Heijne G: Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J Mol Biol. 2000, 300: 1005-1016.PubMedGoogle Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer ELL: Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J Mol Biol. 2001, 305: 567-580.PubMedGoogle Scholar
- Yu CS, Lin CJ, Hwang JK: Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci. 2004, 13: 1402-1406.PubMed CentralPubMedGoogle Scholar
- Chi SM, Nam D: WegoLoc: accurate prediction of protein subcellular localization using weighted Gene Ontology terms. Bioinformatics. 2012, 28: 1028-1030.PubMedGoogle Scholar
- Park BH, Karpinets TV, Syed MH, MR L u, Uberbacher EC: CAZymes Analysis Toolkit (CAT): Web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database. Glycobiology. 2010, 20: 1574-1584.PubMedGoogle Scholar
- 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: 2731-2739.PubMed CentralPubMedGoogle Scholar
- Robbertse B, Yoder RJ, Boyd A, Reeves J, Spatafora JW: Hal: an automated pipeline for phylogenetic analyses of genomic data. PLoS Currents. 2011, 3: RRN1213-PubMed CentralPubMedGoogle Scholar
- de Wit PJGM, van der Burgt A, Okmen B, Stergiopoulos I, Abd-Elsalam KA, Aerts AL, Bahkali AH, Beenen HG, Chettri P, Cox MP, Datema E, de Vries RP, Dhillon B, Ganley AR, Griffiths SA, Guo YA, Hamelin RC, Henrissat B, Kabir MS, Jashni MK, Kema G, Klaubauf S, Lapidus A, Levasseur A, Lindquist E, Mehrabi R, Ohm RA, Owen TJ, Salamov A, Schwelm A: The genomes of the fungal plant pathogens Cladosporium fulvum and Dothistroma septosporum reveal adaptation to different hosts and lifestyles but also signatures of common ancestry. PLoS Genet. 2012, 8:Google Scholar
- Goodwin SB, Ben M'Barek S, Dhillon B, Wittenberg AHJ, Crane CF, Hane JK, Foster AJ, Van der Lee TAJ, Grimwood J, Aerts A, Antoniw J, Bailey A, Bluhm B, Bowler J, Bristow J, van der Burgt A, Canto-Canche B, Churchill ACL, Conde-Ferraez L, Cools HJ, Coutinho PM, Csukai M, Dehal P, De Wit P, Donzelli B, van de Geest HC, Van Ham RCHJ, Hammond-Kosack KE, Henrissat B, Kilian A, et al: Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis. PLoS Genet. 2011, 7:Google Scholar
- Amselem J, Cuomo CA, van Kan JAL, Viaud M, Benito EP, Couloux A, Coutinho PM, de Vries RP, Dyer PS, Fillinger S, Fournier E, Gout L, Hahn M, Kohn L, Lapalu N, Plummer KM, Pradier JM, Quevillon E, Sharon A, Simon A, ten Have A, Tudzynski B, Tudzynski P, Wincker P, Andrew M, Anthouard V, Beever RE, Beffa R, Benoit I, Bouzid O, et al: Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea. PLoS Genet. 2011, 7:Google Scholar
- Martin F, Kohler A, Murat C, Balestrini R, Coutinho PM, Jaillon O, Montanini B, Morin E, Noel B, Percudani R, Porcel B, Rubini A, Amicucci A, Amselem J, Anthouard V, Arcioni S, Artiguenave F, Aury JM, Ballario P, Bolchi A, Brenna A, Brun A, Buee M, Cantarel B, Chevalier G, Couloux A, Da Silva C, Denoeud F, Duplessis S, Ghignone S, et al: Périgord black truffle genome uncovers evolutionary origins and mechanisms of symbiosis. Nature. 2010, 464: 1033-1038.PubMedGoogle Scholar
- Arnaud MB, Chibucos MC, Costanzo MC, Crabtree J, Inglis DO, Lotia A, Orvis J, Shah P, Skrzypek MS, Binkley G, Miyasato SR, Wortman JR, Sherlock G: The Aspergillus Genome Database, a curated comparative genomics resource for gene, protein and sequence information for the Aspergillus research community. Nucleic Acids Res. 2010, 38: 420-427.Google Scholar
- Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007, 56: 564-577.PubMedGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D: ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011, 27: 1164-1165.PubMedGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321.PubMedGoogle Scholar
- Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19: 301-302.PubMedGoogle Scholar
- Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9: 286-298.PubMedGoogle Scholar
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Biol. 2004, 5:Google Scholar
- Hane JK, Rouxel T, Howlett BJ, Kema GHJ, Goodwin SB, Oliver RP: A novel mode of chromosomal evolution peculiar to filamentous Ascomycete fungi. Genome Biol. 2011, 12:Google Scholar
- Gremme G, Steinbiss S, Kurtz S: GenomeTools: A comprehensive software library for efficient processing of structured genome annotations. IEEE/ACM Trans Comput Biol Bioinform. 2013, 10: 645-656.PubMedGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer ELL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database. Nucleic Acids Res. 2012, 40: 290-301.Google Scholar
- Blin K, Medema MH, Kazempour D, Fischbach MA, Breitling R, Takano E, Weber T: antiSMASH 2.0-a versatile platform for genome mining of secondary metabolite producers. Nucleic Acids Res. 2013, 41: 204-212.Google Scholar
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.