Diversity and evolution of ABC proteins in mycorrhiza-forming fungi
© Kovalchuk et al. 2015
Received: 13 April 2015
Accepted: 28 October 2015
Published: 28 December 2015
Transporter proteins are predicted to have an important role in the mycorrhizal symbiosis, due to the fact that this type of an interaction between plants and fungi requires a continuous nutrient and signalling exchange. ABC transporters are one of the large groups of transporter proteins found both in plants and in fungi. The crucial role of plant ABC transporters in the formation of the mycorrhizal symbiosis has been demonstrated recently. Some of the fungal ABC transporter-encoding genes are also induced during the mycorrhiza formation. However, no experimental evidences of the direct involvement of fungal ABC transporters in this process are available so far. To facilitate the identification of fungal ABC proteins with a potential role in the establishment of the mycorrhizal symbiosis, we have performed an inventory of the ABC protein-encoding genes in the genomes of 25 species of mycorrhiza-forming fungi.
We have identified, manually annotated and curated more than 1300 gene models of putative ABC protein-encoding genes. Out of those, more than 1000 models are predicted to encode functional proteins, whereas about 300 models represent gene fragments or putative pseudogenes. We have also performed the phylogenetic analysis of the identified sequences. The sets of ABC proteins in the mycorrhiza-forming species were compared to the related saprotrophic or plant-pathogenic fungal species. Our results demonstrate the high diversity of ABC genes in the genomes of mycorrhiza-forming fungi. Via comparison of transcriptomics data from different species, we have identified candidate groups of ABC transporters that might have a role in the process of the mycorrhiza formation.
Results of our inventory will facilitate the identification of fungal transporters with a role in the mycorrhiza formation. We also provide the first data on ABC protein-coding genes for the phylum Glomeromycota and for orders Pezizales, Atheliales, Cantharellales and Sebacinales, contributing to the better knowledge of the diversity of this protein family within the fungal kingdom.
KeywordsABC transporters Fungi Mycorrhiza Membrane transport Protein family evolution
Biological functions of different groups of fungal ABC proteins
Biological function in fungi unknown, but Magnaporthe grisea transporter Abc4 is required for appresoria formation. Mammalian ABC-A transporters are involved in lipid transport and metabolism
Export of mating pheromones
Multidrug resistance; a specific group of ABCB transporters in ascomycetes (but not in Agaricomycotina) is associated with the biosynthesis cluster for the hydroxamate-type siderophores and might have a role in their metabolism.
Export of the mitochondria-synthesized precursors of iron-sulfur (Fe/S) clusters; localised to the inner membrane of mitochondria
Heavy-metal resistance; localized to the vacuolar membrane
Transporter of the inner mitochondrial membrane involved in the export of peptides released upon proteolysis by m-AAA protease
Vacuolar import of phosphatidylcholine; bile acid transport
Unknown; often associated with the secondary metabolism gene clusters
Heavy metal resistance; vacuolar uptake of glutathione conjugates; localized to vacuolar membrane
Drug resistance; heavy metal resistance
Import of long-chain fatty acids into peroxisomes
Import of long-chain fatty acids into peroxisomes
Ribosome biogenesis and assembly of translational pre-initiation complexes
Biogenesis of 40S and 60S ribosomal subunits
Regulation of mRNA translation; positive regulator of Gcn2p kinase activity
mRNA export factor
Translation elongation factor
Pleiotropic drug resistance; weak organic acid tolerace; export of toxic metabolites, sterol uptake, translocation of membrane phospholipids and quorum sensing
Unknown; are characterised by a presence of EGF domain
Part of the CCR4-NOT transcriptional regulatory complex
Numerous species of fungi are known to form mycorrhizae, symbiotic interactions between fungi and higher plants. It is estimated that more than 90 % of all plant species, including forest trees, wild grasses and many crops are participating in mycorrhizal interactions . Mycorrhiza is a rather broad term that describes structurally heterogeneous types of interactions between fungi and roots of higher plants. In case of ectomycorrhiza, fungal hyphae surround the root tips and can grow between epidermal cells, but never enter cell lumen. On contrary, endomycorrhiza-forming fungi penetrate inside living cells of the root epidermis and cortex . Endomycorrhizae can be further subdivided into ericoid, arbutoid, orchid and arbuscular mycorrhiza. Mycorrhiza-forming fungi are also very diverse taxonomically, as they represent different evolutionary lineages from three fungal phyla, Glomeromycota, Ascomycota and Basidiomycota. Phylogenetic analysis indicates that mycorrhizal symbionts have evolved independently in different lineages of fungi . Establishment of mycorrhizal interaction is a complex process that requires a continuous nutrient and signal exchange between the partners . Transporter proteins are expected to play a key role in this process. Details of those interactions are gradually being elucidated, but many of the involved molecular components remain unknown. A number of works on different plant models have shown that specific plant ABC transporters are required for the formation of arbuscular mycorrhiza, a widespread type of mutualistic interaction between plants and fungi of the phylum Glomeromycota [8–10], and, furthermore, it was demonstrated that the corresponding transporter acts as strigolactone exporter in petunia . It is not known yet whether any of fungal ABC transporters participate in the formation of mycorrhizal interactions.
Distribution of subfamilies of ABC proteins among analysed species of mycorrhiza-forming fungi
Rhizophagus irregularis DAOM 181602
Meliniomyces bicolor E
Meliniomyces variabilis F
Oidiodendron majus Zn
Cenococcum geophilum 1.58
Tuber melanosporum Mel28
Choiromyces venosus 120613-1
Terfezia boudieri S1
Amanita muscaria Koide
Cortinarius glaucopus AT 2004 276
Hebeloma cylindrosporum h7
Laccaria amethystina LaAM-08-1
Laccaria bicolor S238N-H82
Tricholoma matsutake 945
Paxillus involutus ATCC 200175
Paxillus rubicundus Ve08.2 h10
Pisolithus microcarpus 441
Pisolithus tinctorius Marx 270
Scleroderma citrinum Foug A
Suillus luteus UH-Slu-Lm8-n1
Piloderma croceum F 1598
Tulasnella calospora AL13/4D
Piriformospora indica DSM 11827
Sebacina vermifera MAFF 305830
List of ABC-protein encoding genes significantly up-regulated in mycorrhiza-forming mycelium (MYC) as compared to free-living mycelium (FLM)
List of ABC-protein encoding genes significantly down-regulated in mycorrhiza-forming mycelium (MYC) as compared to free-living mycelium (FLM)
We were interested in how the fungal ABC genes, and, in particular, ABC transporters were influenced by the adaptations of fungi to their symbiotic lifestyle. These adaptations have evolved independently in different lineages of mycorrhiza-forming fungi , and they show a remarkable degree of evolutionary convergence. If ABC transporters do play a role in the process of mycorrhiza formation, their recruitment also occurred independently, and we might expect to find different transporters in different lineages depending on the type of mycorrhizal symbiosis. We would also expect to see the lineage-specific differences when comparing ABC transporter repertoire between mycorrhiza-forming fungi and their non-symbiotic relatives. Therefore, when possible, we have compared the sets of ABC genes found in the genomes of mycorrhiza-forming fungi, with the situation in related saprotrophic or plant-pathogenic species, paying attention to both quantitative and qualitative differences in the composition of ABC genes.
We have carried out a comprehensive analysis of ABC protein-encoding genes in the 25 fungal genomes that were sequenced by the MGI . We have identified all predicted ABC genes and gene fragments, manually curated and annotated the gene models and classified them into the subfamilies according to the HUGO-proposed scheme . The identified genes were further subdivided into smaller groups defined in our previous analyses [14, 15]. We have also analysed the data on expression of ABC transporter-encoding genes provided by members of the MGI consortium in order to identify those genes that might play a role in the establishment of mycorrhizal symbiosis. We could not identify any group of ABC transporters that would be unique for mycorrhiza-forming fungi. However, functions of many of identified genes are unknown, and it is possible that some of them play a certain role in the establishment of mycorrhizal interactions. Our data provide an insight into the diversity of ABC proteins in the newly sequenced genomes of mycorrhiza-forming fungi. Currently, it is not clear whether any fungal ABC transporters are implicated in the mycorrhiza formation. Results of our surveyshould provide a support for further studies on the possible role of fungal ABC transporters in the establishment of mycorrhizal symbiosis. Additionally, they will contribute to the better understanding of the evolution and dynamics of this protein family in the fungal kingdom.
Below we will present and discuss the most important findings of our analysis separately for each taxonomic group.
Phylum Glomeromycota is one of the oldest fungal lineages that is highly unusual in many ways. Members of this group are obligate mutualistic symbionts forming a particular type of mycorrhizal interaction, known as arbuscular mycorrhiza (AM) [16–18]. We have identified 34 genes and 19 pseudogenes or gene fragments in the genome assembly of Rhizophagus irregularis, a species of Glomeromycota sequenced by the MGI consortium . The distribution of the identified genes among the recognised subfamilies is shown in Table 2 and Additional file 2. Whereas the total number of ABC genes in Rh. irregularis is comparable with the numbers in other species included in our analysis, we could not identify representatives of some of the groups that are commonly found in the genomes of asco- and basidiomycetes (Additional file 1). Absence of genes homologous to Saccharomyces cerevisiae MDL1, YEF3 and YDR061W is particularly remarkable as those groups are nearly ubiquitously present in the genomes of asco- and basidiomycetes analysed so far [14, 15], and YEF3 (encoding the translation elongation factor 3) has shown to be essential for vegetative growth in S. cerevisiae . We cannot exclude the possibility that the mentioned genes are present in the genome, but for some reason they are missing from the available genome assembly. Alternatively, some of them may have been lost by Rh. irregularis as a result of the adaptation to its obligate symbiotic lifestyle.
Some aspects of the organisation of ABC genes in Rh. irregularis show remarkable similarity with the zygomycete Rhizopus oryzae (Additional file 1) . The similarity in the genome organization of Rh. irregularis and species of Mucoromycotina has been noticed previously . One of the most prominent common features of Rhizophagus and Rhizopus is a presence of four ABC-D genes. Nearly all asco- and basidiomycetes have only two ABC-D genes. The products of those genes act as a heterodimeric peroxisomal long-chain fatty acid transporter. Remarkably, some of the species representing the ancient evolutionary lineages within fungal kingdom (e.g., Rhizopus, Rhizophagus and chytridiomycete Spizellomyces punctatus ) harbour additional genes of the subfamily ABC-D that in phylogenetic analysis are placed separately from the proteins of asco- and basidiomycetes (Fig. 3). The origin and function of those additional genes remains enigmatic.
Little is known about the role of Rhizophagus ABC transporters in the formation of arbuscular mycorrhiza. However, the transcript profiling has shown that one of the Rh. irregularis ABC-B transporters is arbuscle-specific, i.e. it is highly expressed in intraradical mycelium, but not during other stages of the Rhizophagus life cycle .
The unusual genetic organisation has substantially hampered the genome assembly of Rh. irregularis , and the available assembly is highly fragmented . Therefore, we cannot exclude that some of the identified genes represent the different alleles of the same gene, and the actual number of ABC genes might be somewhat lower. At the same time, the detected gene fragments might be extended to full-length genes with further improvement of the genome assembly. It is also possible that some of the genes are still missing from the assembly (in particular, we would expect the presence of the ScYEF3 homologue in Rh. irregularis genome).
Class Leotiomycetes/order Helotiales
The number of ABC genes in the three mycorrhiza-forming species of the order Helotiales is considerably higher than in two phytopathogenic species from the same order, Sclerotinia sclerotiorum and Botryotinia fuckeliana (anamorph Botrytis cinerea)  (Additional file 1). The higher number of genes representing the group ABCC5 is probably linked with the higher number of secondary metabolism gene clusters in those species. Representatives of all four groups over-represented in mycorrhiza-forming Helotiales were identified among the genes with the increased transcript abundance in mycorrhiza-forming mycelium (fold change > 2, FDR < 0.05) (Table 3; Additional file 6) . This observation indicates that at least some of them might have a role in the establishment of the mycorrhiza symbiosis.
Cenococcum geophilum is the only known ectomycorrhizal species within the class Dothideomycetes. We have detected 35 genes encoding predicted ABC protein together with 5 putative pseudogenes or gene fragments in the genome of C. geophilum (Table 2 and Additional file 1; Additional file 2). Contrary to the situation observed within the order Helotiales, the number of ABC genes in C. geophilum is somewhat lower than the corresponding numbers in the genomes of the phytopathogenic dothideomycetes Pyrenophora tritici-repentis and Stagonospora nodorum . Those differences are mainly due to the lower number of full-length ABC-B and ABC-G transporters in C. geophilum. Nevertheless, the sets of ABC genes in the genomes of all three representatives of Dothideomycetes show high degree of similarity and there are only a few noteworthy differences (Additional file 1).
Class Pezizomycetes/ order Pezizales
The order Pezizales is represented in our analysis by three species, Tuber melanosporum (family Tuberaceae), Choiromyces venosus (Tuberaceae) and Terfezia boudieri (Pezizaceae). All three species are characterised by rather low number of ABC protein-encoding genes (Table 2 and Additional file 1; Additional file 2), and T. boudieri has the lowest number of ABC genes (29 genes) among mycorrhiza-forming species included in our analysis and, as far as known, among all analysed filamentous ascomycetes . In particular, numbers of ABC-G, ABC-C and full-length ABC-B transporters in the analysed species of Pezizales are lower than in other filamentous ascomycetes (Additional file 1).
The existence of 27 ABC transporters in the genome of T. melanosporum has been reported previously . After careful examination and manual curation of the available gene models, we have identified five cases of misannotation, when a single ABC transporter-encoding gene was split into two separate non-overlapping models (Additional file 2). After their correction, we concluded that the genome of T. melanosporum has 24 genes for ABC transporters and 8 genes encoding ABC proteins without transmembrane domains.
The sets of ABC proteins in all three analysed species of Pezizales show a high level of similarity, and in most cases the corresponding homologues form well-supported groups in our phylogenetic analysis. The presence of the genes representing the group ABCC4 in all three species of Pezizales is noteworthy, as this type of transporters is quite uncommon in ascomycetes.
The transcript abundance of one of the T. melanosporum ABC-G transporters was increased in the ectomycorrhizal root tips compared to the free-living mycelium , indicating its potential role in the establishment of the mycorrhizal symbiosis.
Phylum Basidiomycota/class Agaricomycetes
We have analysed genome sequences of six mycorrhiza-forming species of Agaricales: Amanita muscaria, Cortinarius glaucopus, Hebeloma cylindrosporum, Laccaria amethystina, L. bicolor and Tricholoma matsutake. The data on L. bicolor ABC genes were taken from our previous analysis . All six species are engaged in ectomycorrhiza symbiosis with a wide range of forest trees, including both conifers and broad-leaved trees. The numbers of identified ABC genes in the genomes of mycorrhiza-forming Agaricales range from 33 genes in A. muscaria to 47 genes in C. glaucopus. All of the analysed species, except for H. cylindrosporum, are characterised by the presence of high numbers of gene fragments and putative pseudogenes related to the functional ABC protein-encoding genes, ranging from 15 in A. muscaria to 53 in L. amethystina (Additional file 1).
There is a considerable similarity between the sets of ABC genes in the genomes of the analysed mycorrhiza-forming species and of saprotrophic representatives of Agaricales, Coprinopsis cinerea  and Agaricus bisporus (our unpublished data), with the differences being mostly quantitative. However, both saprotrophic species have a considerably lower number of putative pseudogenes or gene fragments in their genomes. The set of ABC genes in A. muscaria is quite distinct, as its genome lacks some of the genes present in the remaining species (Additional file 1).
The distinctive features of C. glaucopus genome are the high number of the genes belonging to the groups ABCC4 and ABCC7, and the absence of ABC-D transporters and of the group ABCC6. The group ABCC4 is made up by proteins of unknown function nearly ubiquitously present in Agaricomycotina. There are 11 genes of this group in the genome of C. glaucopus, whereas their number in other analysed species of Agaricales ranges from 1 to 5. The group has undergone a remarkable amplification in Polyporales  with some species having up to 16 genes.
A common feature of the analysed species of Agaricales (including saprotrophs) is the high number of genes belonging to the group ABCB6 of half-size ABC-B transporters. To our current knowledge, the group is found only in the species of Agaricomycotina , and the biological role of the corresponding proteins remains obscure.
Analyses of the gene expression during mycorrhiza formation have been performed in two species of Agaricales, L. bicolor and H. cylindrosporum (Tables 3 and 4; Additional file 6) . The identification of a group ABCB6 half-size ABC-B transporter and ScYDR061W homologues in L. bicolor experiment is noteworthy, as biological functions of those two groups have not been elucidated in any of the fungal species yet.
We have analysed the genomes of seven mycorrhiza-forming species of Boletales representing four distinct evolutionary lineages within the order : Suillineae (Suillus luteus), Sclerodermatineae (Scleroderma citrinum, Pisolithus microcarpus and P. tinctorius), Paxillineae (Paxillus involutus and P. rubicundus), and Boletinae (Boletus edulis). We have also compared the data from these 7 species with the sets of ABC transporters identified in two wood-degrading species of Boletales (Coniophora puteana and Serpula lacrymans ) representing the basal lineages within the order . The numbers of ABC-protein encoding genes in the analysed mycorrhiza-forming species of Boletales range from 30 in B. edulis to 42 in P. tinctorius (Table 2 and Additional file 1; Additional file 2). Genomes of many species (in particular, P. involutus and P. microcarpus) also harbour a high number of gene fragments and putative pseudogenes derived from ABC transporter-encoding genes.
There are no pronounced differences in sets of ABC proteins between wood-degrading and mycorrhiza-forming species of Boletales. However, some of the ABC transporters identified in Boletales have noteworthy features. One of the interesting finding of our analysis was the observation that the predicted proteins of the group ABCB4 of half-size ABC-B transporters in B. edulis, Paxillus spp., Pisolithus spp. and S. citrinum differ considerably from their counterparts in other species of Agaricomycotina, including the remaining members of Boletales (S. luteus, C. puteana and S. lacrymans). Sequences from the six mentioned species were also clearly separated in the phylogenetic analysis, forming an isolated branch on the phylogenetic tree (Additional file 3). Remarkably, the sequences of their predicted transmembrane domains had a higher level of sequence similarity to the ABCB4 proteins from other fungi than the sequences of their nucleotide-binding domains. Most plausible explanation for this observation would be that the changes in the nucleotide sequence ABCB4 gene occurred in the common ancestor of Boletineae, Paxillineae and Sclerodermatineae after the separation of the lineage of Suillineae. It is not entirely clear what could be the cause of these changes, and whether the derived genes retained their biological function. The assignment of these proteins to the group ABCB4 is tentative, and further analysis will be required to clarify their origin and phylogenetic affinities.
Another distinct feature of Boletales is a presence of the gene encoding an unusual ABC-C transporter with a predicted AMP-binding N-terminal domain. This gene was identified before in S. lacrymans, C. puteana and polyporoid fungus Ceriporiopsis subvermispora . We could detect this gene in all the analysed species of Boletales except for S. luteus. The identified sequences were placed in a separate isolated group in our phylogenetic analysis. Remarkably, in all species the sequence coding for the predicted AMP-binding domain is present, further strengthening our suggestion that this domain indeed makes part of the transporter protein. It has been hypothesised that this unusual structure could be essential for tight coupling of the biosynthesis of yet unknown metabolite with its transport, and indeed, the gene is associated with an uncharacterised NRPS clusters in C. puteana and S. lacrymans . However, we could not detect similar associations of this gene with secondary metabolism gene clusters in mycorrhiza-forming species of Boletales.
Similarly to Agaricales, species of Boletales have a diverse set of genes of unknown function belonging to the Agaricomycotina-specific group ABCB6 of half-size ABC-B transporters. Transcriptomics experiments have identified several ABC genes of Boletales up-regulated in ectomycorrhiza-forming mycelium (Table 3; Additional file 6) .
The small order Atheliales is represented in our analysis by a single species, a broad host-range mycorrhizal symbiont Piloderma croceum, widespread in boreal and temperate forests. Atheliales are most closely related to Boletales and Agaricales . However, the number of ABC genes in Piloderma (48) exceeds the numbers found in the analysed species of those two orders. It also has a considerable number (15) of gene fragments and putative pseudogenes. The set of ABC transporters found in the genome of P. croceum is characterised in particular by the high diversity of ABC-G transporters, with some of the genes without counterparts in other analysed species of Agaricomycotina (Additional file 5). Other noteworthy features of P. croceum are a high number of full-length ABC-B transporters and a diverse set of half-size ABCB6 genes. Three P. croceum ABC genes were significantly induced during mycorrhiza formation: full-length ABC-B (ABCB2) transporter, half-size ABC-B transporter (ABCB6), and a half-size ABC-G transporter (ABCG6) (Table 3; Additional file 6) .
This order is one of the ancient lineages within Agaricomycotina . It is represented in our analysis by Tulasnella calospora, a widespread symbiont of green orchids with the world-wide distribution. In association with its host plants, T. calospora forms a distinct type of endomycorrhiza known as orchid-type mycorrhiza. The genome of T. calospora harbours 40 predicted ABC genes and 3 gene fragments likely derived from the intact ABC genes. The identified set of ABC genes is in many aspects “ordinary”, i.e. there are no prominent distinctive features. All of the groups of ABC proteins commonly found in Agaricomycotina are represented, and there are no “odd” genes. We could only mention here the relatively high number of full-length ABC-B genes, and the lower diversity of the group ABCB6 of half-size ABC-B transporters as compared with the species of Agaricales and Boletales.
Transcriptomics analysis has shown that 9 genes encoding ABC transporters and 2 genes encoding soluble ABC proteins were up-regulated in T. calospora during mycorrhiza formation (Table 3; Additional file 6) .
Together with the order Cantharellales, this is one of the basal lineages within Agaricomycotina. In our analysis, we have included the orchid symbiont Sebacina vermifera, which was compared with the root endophyte Piriformospora indica . The sets of ABC genes in the analysed species show considerable similarity (Additional file 1). The sequences of the half-size ABC-B transporters from the two species, which we provisionally assign to group ABCB4, are highly unusual and differ significantly from the corresponding sequences found in other Agaricomycotina, in particular in the amino acid sequence of their nucleotide-binding domain. Both proteins are also considerably shorter than ABCB4 proteins in other species of fungi. In our phylogenetic analysis the two sequences formed an isolated group of unclear relationship (Additional file 4). Their affinities with the group ABCB4 require additional studies.
Transcriptomics analysis has identified a number of S. vermifera ABC genes as being significantly up-regulated in mycorrhiza-forming mycelium (Table 3).
Studies on several model plant species have demonstrated that particular plant ABC transporters are required for the formation of the arbuscular mycorrhiza [8–10]. At the same time, the extensive exchange of nutrients between the both partners engaged in mycorrhiza formation probably requires the action of additional transporter proteins, not only on the host side, but also on the side of the colonising fungus. Currently, it remains unknown whether any of fungal ABC transporters contribute to the formation of any type of mycorrhizal symbiosis.
Our inventory provides an important insight into the diversity of ABC transporters and soluble ABC proteins in the genomes of mycorrhiza-forming fungi. However, the interpretation of the obtained results is considerably hampered by the scarcity of the available information on the biological functions of fungal ABC transporters. Results of our inventory indicate that there are no principal differences in sets of ABC proteins encoded by mycorrhiza-forming fungi, on one hand, and saprotrophic or phytopathogenic fungi, on other hand. We could not identify any ABC transporters that could be classified as “mycorrhiza-specific”. Moreover, the differences between members of various evolutionary lineages of fungi are much more pronounced than between the species occupying different ecological niches. Thus, the principal question, whether any of fungal ABC transporters play a role in the establishment and functioning of mycorrhizal symbiosis, remains open. However, the available data on gene expression in several fungal species indicate that some ABC transporter-encoding genes are transcriptionally up-regulated during the mycorrhiza formation [12, 17, 20, 24]. The performed transcriptomics studies have repeatedly identified members of some groups of ABC transporters as induced upon mycorrhiza formation, e.g. full-length ABC-B transporters of the group ABCB2, half-size ABC-B transporters of the group ABCB6, the transporters belonging to the group ABCG6 and the soluble ABC proteins of the group ABCI2. The group ABCB2 is widely distributed among fungi and nearly ubiquitously present in members of different lineages. Transporters of this group are also known as MDR transporters, as some of them contribute to the resistance to antifungal compounds. The formation of paralogous genes within this group is evident from our phylogenetic analysis (Fig. 2). It is tempting to speculate that some of those paralogous copies have acquired new functions as a result of an adaptation to the symbiotic life style.
The group ABCB6 is particularly interesting as it seems to be specific to the species of Agaricomycotina. No members of this group have been identified up till now from other fungal lineages. Members of this group have a typical structure of half-size ABC transporters, with N-terminal transmembrane domain and C-terminal nucleotide binding domain. They show a considerable degree of diversification within Agaricomycotina. In our phylogenetic analysis, four clusters with high bootstrap support could be recognised within the group. The highest diversity of its members is observed in the crown groups of Agaricomycotina, i.e. orders of Agaricales and Boletales, whereas the diversity is lower in basal lineages (like Sebacinales and Cantharellales). ABCB6 genes have been identified in all analysed Agaricomycotina regardless of their life style and ecological niche with the sole exception of opportunistic pathogen Cryptococcus neoformans. None of the members of this group has been characterised experimentally so far, but their wide distribution and a high diversity indicate that they might play an important role in basidiomycetes.
ABCG6 is another group encompassing transporters of unknown function. It is widely but unevenly distributed among fungal lineages, as it is completely missing, for example, from all analysed species of Boletales. The last group ABCI2 is made up by the homologues of the uncharacterised S. cerevisiae gene YDR061W. Members of this group are ubiquitously present in fungal genomes, but their biological role remains obscure.
Our analysis and the available transcriptomics data have pinpointed those of ABC transporters that might contribute to the mycorrhiza formation in different fungal lineages. Whereas some of the transporters have been repeatedly identified in several species, there is no uniform picture, and many of the genes have been up-regulated in just a single species. This is not particularly surprising, as mycorrhiza-forming fungi are very diverse phylogenetically, and the ability to form mycorrhiza-like symbiosis has apparently evolved independently in several evolutionary lineages of fungi . Furthermore, there are considerable structural and physiological differences between various types of mycorrhiza (ectomycorrhiza, ericoid, orchid and arbuscular mycorrhiza), and fungi engaged in each of those type of symbiotic interactions likely have evolved their own specific adaptations. Additional work will be required to clarify the differences in nutrient flow and exchange and in signalling in various types of mycorrhiza symbiosis and, in particular, contribution of different classes of transporter proteins to those processes.
We report here the results of our large-scale effort towards the complete inventory of ABC protein-encoding genes in the genomes of mycorrhiza-forming fungi. Results of our inventory show that the sets of ABC genes differ considerably in various phylogenetic lineages, both in total number of genes and in their distribution among the subfamilies. Genomes of Helotiales are particularly rich in ABC transporter-encoding genes, having almost twice as many corresponding genes as species of Pezizales or some Boletales. Therefore, it is difficult to draw a common pattern for all the analysed groups. The comparison between the mycorrhiza-forming species and their saprotrophic or parasitic relatives also revealed the considerable differences between various evolutionary lineages of fungi. Thus, mycorrhiza-forming species of Helotiales have higher numbers of ABC genes compared with the analysed plant pathogens of the same order, whereas in Dothideomycetes we have observed an opposite trend (however, here the compared species were not very closely related). In Agaricales and Boletales the differences between the saprotrophs and mycorrhiza-formers were even less pronounced. However, in many of the analysed species we have observed the formation of species- or lineage-specific groups of paralogous genes. Neofunctionalisation within such groups might be one of the forces contributing to the adaptation to the symbiotic life style.
Our results also contribute to the better knowledge of the evolution and diversity of fungal soluble ABC proteins and ABC transporters. We provide the first data on ABC protein-coding genes for the members of one fungal phylum (Glomeromycota) and four fungal orders (Pezizales, Atheliales, Cantharellales and Sebacinales). The information on ABC genes from the genome of the glomeromycete Rhizomyces irregularis is particularly important for the understanding of the evolution of this gene family in the ancient lineages of the fungal evolutionary tree before the divergence of lineages leading to modern Ascomycota and Basidiomycota.
Finally, the reported data should provide a contribution to prioritising target genes for the forthcoming functional characterisation, and in this way essentially advance the studies of fungal ABC transporters in general.
The genomic DNA sequences of the mycorrhiza-forming fungi (except for the previously sequenced species Laccaria bicolor , Tuber melanosporum  and Piriformospora indica ) were obtained at the Joint Genome Institute (JGI) of the US Department of Energy (DOE) in Walnut Creek (California), as part of the Mycorrhizal Genomic Initiative (MGI) project. The genomes were produced as described by . The results from gene prediction and annotation in the above genomes are available for searching in Mycocosm  at the JGI portal (http://jgi.doe.gov/fungi).
ABC genes identification
To identify gene loci encoding ABC proteins in the fungal genomes, multiple tblastn and blastp searches against selected genomes were performed at the website of the Fungal Genomics Program of the Department of Energy Joint Genome Institute (JGI) . Sequences of S. cerevisiae and Coprinopsis cinerea ABC proteins representing all known subfamilies were used as queries. All hits producing E-values below 10-4 were further analysed. Gene models were manually curated and, when necessary, the positions of N and C termini and exon-intron boundaries were adjusted.
Phylogenetic analysis was performed with the program package MEGA6 . Initial assignment of identified sequences to the subfamilies was done separately for every species in a following way. All proteins identified in a particular genome were aligned with known ABC proteins from S. cerevisiae, Aspergillus nidulans, C. cinerea and Ustilago maydis using MUSCLE algorithm integrated in MEGA6 package. Neighbor-joining phylogenetic trees were reconstituted based on the obtained alignments. All recognized subfamilies of ABC proteins could be separated on trees produced in such way, allowing the fast assignment of the analysed sequences to the particular subfamily. The detailed analysis was performed separately for each subfamily and, in the case of ABC-B transporters, also for full-length and half-size proteins. Multiple sequence alignments were constructed with the MUSCLE algorithm integrated in MEGA6 package using default settings. Alignments were quality trimmed with Gblocks 0.91 b . Both full-length and trimmed alignments were used to produce phylogenetic trees. Neighbor-joining trees were constructed using Jones-Taylor-Thornton substitution model with 500 bootstrap replications. Maximum-likelihood trees were obtained using Jones-Taylor-Thornton model  with 100 bootstrap replications. Although the topologies of the phylogenetic trees produced with the different algorithms showed some minor differences, the same major groups of ABC proteins were recognized in all reconstructions.
Availability of supporting data
The data sets supporting the results of this article are available in the TreeBase repository under the following link: http://purl.org/phylo/treebase/phylows/study/TB2:S17878.
Academy of Finland is gratefully acknowledged for research funding. The members of the Mycorrhizal Genomics Initiative consortium are acknowledged for making available their unpublished data on the expression of ABC protein-encoding genes.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Dassa E, Bouige P. The ABC of ABCS: a phylogenetic and functional classification of ABC systems in living organisms. Res Microbiol. 2001;152(3-4):211–29.View ArticlePubMedGoogle Scholar
- Lage H. ABC-transporters: implications on drug resistance from microorganisms to human cancers. Int J Antimicrob Agents. 2003;22(3):188–99.View ArticlePubMedGoogle Scholar
- Piddock LJ. Multidrug-resistance efflux pumps - not just for resistance. Nat Rev Microbiol. 2006;4(8):629–36.View ArticlePubMedGoogle Scholar
- Lubelski J, Konings WN, Driessen AJ. Distribution and physiology of ABC-type transporters contributing to multidrug resistance in bacteria. Microbiol Mol Biol Rev. 2007;71(3):463–76.PubMed CentralView ArticlePubMedGoogle Scholar
- Jungwirth H, Kuchler K. Yeast ABC transporters-- a tale of sex, stress, drugs and aging. FEBS Lett. 2006;580(4):1131–8.View ArticlePubMedGoogle Scholar
- Bonfante P, Genre A. Mechanisms underlying beneficial plant-fungus interactions in mycorrhizal symbiosis. Nat Commun. 2010;1:48.View ArticlePubMedGoogle Scholar
- Hibbett DS, Gilbert LB, Donoghue MJ. Evolutionary instability of ectomycorrhizal symbioses in basidiomycetes. Nature. 2000;407(6803):506–8.View ArticlePubMedGoogle Scholar
- Zhang Q, Blaylock LA, Harrison MJ. Two Medicago truncatula half-ABC transporters are essential for arbuscule development in arbuscular mycorrhizal symbiosis. Plant Cell. 2010;22(5):1483–97.PubMed CentralView ArticlePubMedGoogle Scholar
- Gutjahr C, Radovanovic D, Geoffroy J, Zhang Q, Siegler H, Chiapello M, et al. The half-size ABC transporters STR1 and STR2 are indispensable for mycorrhizal arbuscule formation in rice. Plant J. 2012;69(5):906–20.View ArticlePubMedGoogle Scholar
- Kretzschmar T, Kohlen W, Sasse J, Borghi L, Schlegel M, Bachelier JB, et al. A petunia ABC protein controls strigolactone-dependent symbiotic signalling and branching. Nature. 2012;483(7389):341–4.View ArticlePubMedGoogle Scholar
- Plett JM, Martin F. Blurred boundaries: lifestyle lessons from ectomycorrhizal fungal genomes. Trends Genet. 2011;27(1):14–22.View ArticlePubMedGoogle Scholar
- Kohler A, Kuo A, Nagy LG, Morin E, Barry KW, Buscot F, et al. Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. Nat Genet. 2015;47(4):410–5.View ArticlePubMedGoogle Scholar
- Dean M, Rzhetsky A, Allikmets R. The human ATP-binding cassette (ABC) transporter superfamily. Genome Res. 2001;11(7):1156–66.View ArticlePubMedGoogle Scholar
- Kovalchuk A, Driessen AJ. Phylogenetic analysis of fungal ABC transporters. BMC Genomics. 2010;11:177.PubMed CentralView ArticlePubMedGoogle Scholar
- Kovalchuk A, Lee YH, Asiegbu FO. Diversity and evolution of ABC proteins in basidiomycetes. Mycologia. 2013;105(6):1456–70.View ArticlePubMedGoogle Scholar
- Martin F, Gianinazzi-Pearson V, Hijri M, Lammers P, Requena N, Sanders IR, et al. The long hard road to a completed Glomus intraradices genome. New Phytol. 2008;180(4):747–50.View ArticlePubMedGoogle Scholar
- Tisserant E, Kohler A, Dozolme-Seddas P, Balestrini R, Benabdellah K, Colard A, et al. The transcriptome of the arbuscular mycorrhizal fungus Glomus intraradices (DAOM 197198) reveals functional tradeoffs in an obligate symbiont. New Phytol. 2012;193(3):755–69.View ArticlePubMedGoogle Scholar
- Tisserant E, Malbreil M, Kuo A, Kohler A, Symeonidi A, Balestrini R, et al. Genome of an arbuscular mycorrhizal fungus provides insight into the oldest plant symbiosis. Proc Natl Acad Sci U S A. 2013;110(50):20117–22.PubMed CentralView ArticlePubMedGoogle Scholar
- Sandbaken MG, Lupisella JA, Di Domenico B, Chakraburtty K. Protein synthesis in yeast. Structural and functional analysis of the gene encoding elongation factor 3. J Biol Chem. 1990;265(26):15838–44.PubMedGoogle Scholar
- Martin F, Kohler A, Murat C, Balestrini R, Coutinho PM, Jaillon O, et al. Perigord black truffle genome uncovers evolutionary origins and mechanisms of symbiosis. Nature. 2010;464(7291):1033–8.View ArticlePubMedGoogle Scholar
- Binder M, Hibbett DS. Molecular systematics and biological diversification of Boletales. Mycologia. 2006;98(6):971–81.View ArticlePubMedGoogle Scholar
- Hibbett DS. A phylogenetic overview of the Agaricomycotina. Mycologia. 2006;98(6):917–25.View ArticlePubMedGoogle Scholar
- Zuccaro A, Lahrmann U, Guldener U, Langen G, Pfiffi S, Biedenkopf D, et al. Endophytic life strategies decoded by genome and transcriptome analyses of the mutualistic root symbiont Piriformospora indica. PLoS Pathog. 2011;7(10):e1002290.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin F, Aerts A, Ahren D, Brun A, Danchin EG, Duchaussoy F, et al. The genome of Laccaria bicolor provides insights into mycorrhizal symbiosis. Nature. 2008;452(7183):88–92.View ArticlePubMedGoogle Scholar
- Grigoriev IV, Nordberg H, Shabalov I, Aerts A, Cantor M, Goodstein D, et al. The genome portal of the Department of Energy Joint Genome Institute. Nucleic Acids Res. 2012;40(Database issue):D26–32.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Talavera G, Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007;56(4):564–77.View ArticlePubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992;8(3):275–82.PubMedGoogle Scholar
- Dean M, Annilo T. Evolution of the ATP-binding cassette (ABC) transporter superfamily in vertebrates. Annu Rev Genomics Hum Genet. 2005;6:123–42.View ArticlePubMedGoogle Scholar
- Wenzel JJ, Piehler A, Kaminski WE. ABC A-subclass proteins: gatekeepers of cellular phospho- and sphingolipid transport. Front Biosci. 2007;12:3177–93.View ArticlePubMedGoogle Scholar
- Gupta A, Chattoo BB. Functional analysis of a novel ABC transporter ABC4 from Magnaporthe grisea. FEMS Microbiol Lett. 2008;278(1):22–8.View ArticlePubMedGoogle Scholar
- Kuchler K, Sterne RE, Thorner J. Saccharomyces cerevisiae STE6 gene product: a novel pathway for protein export in eukaryotic cells. EMBO J. 1989;8(13):3973–84.PubMed CentralPubMedGoogle Scholar
- McGrath JP, Varshavsky A. The yeast STE6 gene encodes a homologue of the mammalian multidrug resistance P-glycoprotein. Nature. 1989;340(6232):400–4.View ArticlePubMedGoogle Scholar
- Christensen PU, Davey J, Nielsen O. The Schizosaccharomyces pombe mam1 gene encodes an ABC transporter mediating secretion of M-factor. Mol Gen Genet. 1997;255(2):226–36.View ArticlePubMedGoogle Scholar
- Kragl C, Schrettl M, Eisendle M, Haas H. The Aspergillus fumigatus ABC transporter AbcB is involved in excretion of siderophore breakdown product. Edinburgh, Scotland, UK: Proceedings of the 9th European Conference on Fungal Genetics; 2008. p. 130. 5-8 April 2008.Google Scholar
- von Dohren H. A survey of nonribosomal peptide synthetase (NRPS) genes in Aspergillus nidulans. Fungal Genet Biol. 2009;46 Suppl 1:S45–52.View ArticleGoogle Scholar
- Del Sorbo G, Schoonbeek H, de Waard MA. Fungal transporters involved in efflux of natural toxic compounds and fungicides. Fungal Genet Biol. 2000;30(1):1–15.View ArticlePubMedGoogle Scholar
- Kispal G, Csere P, Guiard B, Lill R. The ABC transporter Atm1p is required for mitochondrial iron homeostasis. FEBS Lett. 1997;418(3):346–50.View ArticlePubMedGoogle Scholar
- Kispal G, Csere P, Prohl C, Lill R. The mitochondrial proteins Atm1p and Nfs1p are essential for biogenesis of cytosolic Fe/S proteins. EMBO J. 1999;18(14):3981–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Ortiz DF, Kreppel L, Speiser DM, Scheel G, McDonald G, Ow DW. Heavy metal tolerance in the fission yeast requires an ATP-binding cassette-type vacuolar membrane transporter. EMBO J. 1992;11(10):3491–9.PubMed CentralPubMedGoogle Scholar
- Ortiz DF, Ruscitti T, McCue KF, Ow DW. Transport of metal-binding peptides by HMT1, a fission yeast ABC-type vacuolar membrane protein. J Biol Chem. 1995;270(9):4721–8.View ArticlePubMedGoogle Scholar
- Mendoza-Cozatl DG, Zhai Z, Jobe TO, Akmakjian GZ, Song WY, Limbo O, et al. Tonoplast-localized Abc2 transporter mediates phytochelatin accumulation in vacuoles and confers cadmium tolerance. J Biol Chem. 2010;285(52):40416–26.PubMed CentralView ArticlePubMedGoogle Scholar
- Young L, Leonhard K, Tatsuta T, Trowsdale J, Langer T. Role of the ABC transporter Mdl1 in peptide export from mitochondria. Science. 2001;291(5511):2135–8.View ArticlePubMedGoogle Scholar
- Chloupkova M, LeBard LS, Koeller DM. MDL1 is a high copy suppressor of ATM1: evidence for a role in resistance to oxidative stress. J Mol Biol. 2003;331(1):155–65.View ArticlePubMedGoogle Scholar
- Ortiz DF, St Pierre MV, Abdulmessih A, Arias IM. A yeast ATP-binding cassette-type protein mediating ATP-dependent bile acid transport. J Biol Chem. 1997;272(24):15358–65.View ArticlePubMedGoogle Scholar
- Gulshan K, Moye-Rowley WS. Vacuolar import of phosphatidylcholine requires the ATP-binding cassette transporter Ybt1. Traffic. 2011;12(9):1257–68.PubMed CentralView ArticlePubMedGoogle Scholar
- Li ZS, Szczypka M, Lu YP, Thiele DJ, Rea PA. The yeast cadmium factor protein (YCF1) is a vacuolar glutathione S-conjugate pump. J Biol Chem. 1996;271(11):6509–17.View ArticlePubMedGoogle Scholar
- Petrovic S, Pascolo L, Gallo R, Cupelli F, Ostrow JD, Goffeau A, et al. The products of YCF1 and YLL015w (BPT1) cooperate for the ATP-dependent vacuolar transport of unconjugated bilirubin in Saccharomyces cerevisiae. Yeast. 2000;16(6):561–71.View ArticlePubMedGoogle Scholar
- Klein M, Mamnun YM, Eggmann T, Schuller C, Wolfger H, Martinoia E, et al. The ATP-binding cassette (ABC) transporter Bpt1p mediates vacuolar sequestration of glutathione conjugates in yeast. FEBS Lett. 2002;520(1-3):63–7.View ArticlePubMedGoogle Scholar
- Katzmann DJ, Hallstrom TC, Voet M, Wysock W, Golin J, Volckaert G, et al. Expression of an ATP-binding cassette transporter-encoding gene (YOR1) is required for oligomycin resistance in Saccharomyces cerevisiae. Mol Cell Biol. 1995;15(12):6875–83.PubMed CentralView ArticlePubMedGoogle Scholar
- Cui Z, Hirata D, Tsuchiya E, Osada H, Miyakawa T. The multidrug resistance-associated protein (MRP) subfamily (Yrs1/Yor1) of Saccharomyces cerevisiae is important for the tolerance to a broad range of organic anions. J Biol Chem. 1996;271(25):14712–6.View ArticlePubMedGoogle Scholar
- Nagy Z, Montigny C, Leverrier P, Yeh S, Goffeau A, Garrigos M, et al. Role of the yeast ABC transporter Yor1p in cadmium detoxification. Biochimie. 2006;88(11):1665–71.View ArticlePubMedGoogle Scholar
- Hettema EH, van Roermund CW, Distel B, van den BM, Vilela C, Rodrigues-Pousada C, et al. The ABC transporter proteins Pat1 and Pat2 are required for import of long-chain fatty acids into peroxisomes of Saccharomyces cerevisiae. EMBO J. 1996;15(15):3813–22.PubMed CentralPubMedGoogle Scholar
- Dong J, Lai R, Nielsen K, Fekete CA, Qiu H, Hinnebusch AG. The essential ATP-binding cassette protein RLI1 functions in translation by promoting preinitiation complex assembly. J Biol Chem. 2004;279(40):42157–68.View ArticlePubMedGoogle Scholar
- Kispal G, Sipos K, Lange H, Fekete Z, Bedekovics T, Janaky T, et al. Biogenesis of cytosolic ribosomes requires the essential iron-sulphur protein Rli1p and mitochondria. EMBO J. 2005;24(3):589–98.PubMed CentralView ArticlePubMedGoogle Scholar
- Dong J, Lai R, Jennings JL, Link AJ, Hinnebusch AG. The novel ATP-binding cassette protein ARB1 is a shuttling factor that stimulates 40S and 60S ribosome biogenesis. Mol Cell Biol. 2005;25(22):9859–73.PubMed CentralView ArticlePubMedGoogle Scholar
- Garcia-Barrio M, Dong J, Ufano S, Hinnebusch AG. Association of GCN1-GCN20 regulatory complex with the N-terminus of eIF2alpha kinase GCN2 is required for GCN2 activation. EMBO J. 2000;19(8):1887–99.PubMed CentralView ArticlePubMedGoogle Scholar
- Kozak L, Gopal G, Yoon JH, Sauna ZE, Ambudkar SV, Thakurta AG, et al. Elf1p, a member of the ABC class of ATPases, functions as a mRNA export factor in Schizosacchromyces pombe. J Biol Chem. 2002;277(37):33580–9.View ArticlePubMedGoogle Scholar
- Kamath A, Chakraburtty K. Role of yeast elongation factor 3 in the elongation cycle. J Biol Chem. 1989;264(26):15423–8.PubMedGoogle Scholar
- Sarthy AV, McGonigal T, Capobianco JO, Schmidt M, Green SR, Moehle CM, et al. Identification and kinetic analysis of a functional homolog of elongation factor 3, YEF3 in Saccharomyces cerevisiae. Yeast. 1998;14(3):239–53.View ArticlePubMedGoogle Scholar
- Prasad R, Goffeau A. Yeast ATP-binding cassette transporters conferring multidrug resistance. Annu Rev Microbiol. 2012;66:39–63.View ArticlePubMedGoogle Scholar
- Liu HY, Chiang YC, Pan J, Chen J, Salvadore C, Audino DC, et al. Characterization of CAF4 and CAF16 reveals a functional connection between the CCR4-NOT complex and a subset of SRB proteins of the RNA polymerase II holoenzyme. J Biol Chem. 2001;276(10):7541–8.View ArticlePubMedGoogle Scholar