Fucose as a nutrient ligand for Dikarya and a building block of early diverging lineages
IMA Fungus volume 14, Article number: 17 (2023)
Fucose is a deoxyhexose sugar present and studied in mammals. The process of fucosylation has been the primary focus in studies relating to fucose in animals due to the presence of fucose in Lewis antigens. Very few studies have reported its presence in Fungi, mostly in Mucoromycotina. The constitution of 25% and 12% of this sugar in the carbohydrates of cell wall in the respective Umbelopsis and Mucorales strains boosts the need to bridge the gap of knowledge on fucose metabolism across the fungal tree of life. In the absence of a network map involving fucose proteins, we carried out an in-silico approach to construct the fucose metabolic map in Fungi. We analyzed the taxonomic distribution of 85 protein families in Fungi including diverse early diverging fungal lineages. The expression of fucose-related protein-coding genes proteins was validated with the help of transcriptomic data originating from representatives of early diverging fungi. We found proteins involved in several metabolic activities apart from fucosylation such as synthesis, transport and binding. Most of the identified protein families are shared with Metazoa suggesting an ancestral origin in Opisthokonta. However, the overall complexity of fucose metabolism is greater in Metazoa than in Fungi. Massive gene loss has shaped the evolutionary history of these metabolic pathways, leading to a repeated reduction of these pathways in most yeast-forming lineages. Our results point to a distinctive mode of utilization of fucose among fungi belonging to Dikarya and the early diverging lineages. We speculate that, while Dikarya used fucose as a source of nutrients for metabolism, the early diverging group of fungi depended on fucose as a building block and signaling compound.
Fucose (6-deoxy hexose in the L-configuration) is a relatively less abundant carbohydrate compared to other simple sugars. It can bind to glycans, glycoproteins and glycolipids (Ma et al. 2006). All domains of life contain fucosylated molecules either on the surface of the cell in the form of cell wall polysaccharides in fungi and plants, lipopolysaccharides in bacteria, ABO group antigens in mammals or within the cell. The fucosylation process and its role, as well as the general fucose metabolism, are best described in mammals (Kizuka et al. 2017). However, studies on fucose biology have also been carried out for Caenorhabditis (Yan et al. 2015, 2018) and Drosophila (Ishio et al. 2015). To date, three main fucosylation mechanisms can be distinguished. Terminal fucosylation and core fucosylation of N-linked glycans are common and well-described types of fucose modification, in contrast to O-fucosylation that has been well-characterized only in mammals and in Helicobacter pylori (Schneider et al. 2017). L-Fucose can be bound to serine and threonine residues via O-fucosylation and act as a signaling molecule. Fucose can also be linked to N-glycans via α-linkage to galactose or N-acetylglucosamine (GlcNAc) residues. Such structures are involved in various recognition processes, including fertilization and development, immunological response, cancer antigens hiding or during apoptosis (Staudacher et al. 1999).
The greatest emphasis has been put on studying the role of fucose in mammals. Some fucosyltransferases are known to be essential for the proper gastrulation of mice (Du et al. 2010; Schneider et al. 2017). Fucose is also an important compound in the human gastrointestinal tract as part of mucus glycans and blood group ABO and Lewis antigens (Garber et al. 2021). Fucosylation plays a role in the regulation of immunity and cancer development, synapse formation along with fertilization and host-microbe interactions (Schneider et al. 2017). The best-known interaction mediated by this kind of modification is the colonization of human gastric mucosa by the bacterium H. pylori. This microbe can mimic human Lewis antigen to evade the host immune response (Moran 2008).
So far a few reports on fungal proteins acting on fucose (and perhaps arabinose (Ogawa et al. 2001)) substrates describe their fucose-binding lectins as players in animal recognition and mycorrhiza formation (Lebreton et al. 2021). Lectins are widely distributed oligosaccharide-binding proteins with diverse specificities. They take part in cell–cell interactions and, especially in plants, are believed to be involved in defense against pathogens. In fungi, fucose-binding lectins play a role in Aspergillus fumigatus human infection (Sakai et al. 2019), nematode infection by Arthrobotrys oligospora (Tian et al. 2020) and antifungal activity of Aleuria aurantia L-fucose-specific lectin against Mucor racemosum (Amano et al. 2012).
Fucose was long known as a constituent of the cell wall of Mucoromycota (Bartnicki-Garcia 1968), but the structural underpinnings were described only recently in Mucorales (Lecointe et al. 2019; Mélida et al. 2015). Fucose constitutes up to 25% of all carbohydrates in the cell wall of two analyzed Umbelopsis strains (Muszewska et al. 2021). Remaining subphyla of Mucoromycota, namely Glomeromycotina and Mortierellomycotina, as well as non-Mucorales orders within Mucoromycotina: Umbelopsidales and Endogonales remain understudied. Studies have also shown fucose to comprise from 9 to 25% of the EPS layer in several species of Rhizopus and Mucor (Lecointe et al. 2019; De Ruiter et al. 1992). Furthermore, Mucorales secrete exopolysaccharides with significant amounts of L-fucose (Miyazaki and Irino 1972).
Fucosylated proteins were reported in fungi only once (Grass et al. 2011). Enzymes responsible for fucosylation in fungi remain elusive. Although α-1,6-fucose residues in fungal glycans were confirmed, until today there was no evidence of fungal enzymes responsible for this linkage (Martínez-Duncker et al. 2014).
In previous studies, we identified protein O-α-fucosyltransferase as one of the protein families forming fungal adaptasome, the dynamically evolving part of protein-coding set of genes(Orłowska and Muszewska 2022) and confirmed a high content of fucose in Umbelopsis cell wall (Muszewska et al. 2021). These surprising discoveries and the lack of knowledge on fucose metabolism in fungi led us to ask questions about the occurrence and distribution of proteins associated with fucose in the fungal tree of life (FToL). Hence, we provide a comprehensive screening of proteins involved in fucose metabolism across the fungal kingdom and propose a hypothetical metabolic map.
The starting point for this study was to create the initial list of all possible proteins associated with fucose metabolism across the whole tree of life. We searched the following databases using key words “fucose,” “fucosy,” “fucos”: BRENDA (Chang et al. 2021), UniProt (UniProt Consortium 2021), MetaCyc (Caspi et al. 2016), Pfam (Mistry et al. 2021), TCDB (Saier et al. 2021), UniLectins (Bonnardel et al. 2019). As a result, 75 protein families were identified (38 enzymes, 18 transporters and 19 lectins). For each record from the initial list, we picked reference sequences that were used as a query for searching against the NCBI NR protein sequence database (Sayers et al. 2022) (Additional file 3: Table S1). Reference sequences were chosen based on UniProt annotation score. Prepared queries were used to perform blastp v. 2.8.1 + (Altschul et al. 1990) against NCBI NR database last updated December, 2nd, 2021 limited to fungal taxa with e-value < = 1e−5. Fasta sequences for obtained accession numbers were downloaded using NCBI Entrez (Sayers et al. 2022). For each set of proteins, we performed clustering in CLANS with p-value thresholds of 1e−5, 1e−10, 1e−15 and 1e−20, expulsive force exponent set to 2 and remaining parameters set as default values (Frickey and Lupas 2004). To eliminate false positives, we analyzed resulting clusters in terms of putative protein function based on domain architecture predicted by PfamScan with default parameters against Pfam database v. 33 and v. 35 (Mistry et al. 2007) and InterProScan 5.60–92.0 (Jones et al. 2014). Fifteen enzymes, 2 groups of transporters and 9 lectins were confirmed to be present in fungal protein sets derived from genomic sequencing (proteomes). Logos representing sequence alignment were prepared in WebLogo 3.7.12 (Crooks et al. 2004). Metabolic maps were built with Inkscape (Inkscape Project 2020) based on pathways stored in MetaCyc database (Caspi et al. 2016). Cellular localization and the presence of signal peptides were predicted by using WoLF PSORT v. 0.2 (organism type—fungi, other parameters set with defaults) (Horton et al. 2007) and SignalP v 5.0 (organism type—eukaryotes, other parameters set with defaults) (Almagro Armenteros et al. 2019). Alignments were calculated using MAFFT v7.490 (–localpair) (Katoh et al. 2002), and trimmed with TrimAl 1.2rev59 (-automated1) (Capella-Gutiérrez, et al. 2009). Phylogenetic trees were constructed with IQTree 2.0.7 (-m MFP -alrt 0 -t RANDOM -nt AUTO) (Minh et al. 2020). Phylogenetic trees included as an outgroup fucosyltransferases identified with the same query list Additional file 3: Table S1) in Capsaspora owczarzaki (GCF_000151315.2) (Suga et al. 2013), Salpingoeca arctica (GCF_001186125.1) (Ruiz-Trillo et al. 2007), Fonticula alba (GCF_000388065.1) (Ruiz-Trillo et al. 2007), Salpingoeca rosetta (GCF_000188695.1) (Fairclough et al. 2013), Monosiga brevicollis (GCF_000002865.3) (King et al. 2008), Mus musculus (GCF_000001635.27) (Church et al. 2011), Drosophila melanogaster (GCF_000001215.4) (Matthews et al. 2015) and Homo sapiens (GCF_000001405.40_GRCh38.p14) (Lander et al. 2001).
To verify the binding affinity of identified fungal Fringe homologs (O-fucosylpeptide 3-beta-N-acetylglucosaminyltransferase), we clustered them with reference sequences (Uniprot: O09008 and Q96EU7) using CLANS (same parameters as above) (Frickey and Lupas 2004). We performed structure prediction with ColabFold v 1.4.0 (default parameters) (Mirdita et al. 2022; Jumper et al. 2021). Structures of proteins and ligands were prepared for molecular docking with Open Babel 3.0.1 (O’Boyle et al. 2011). Docking was performed with AutoDock Vina 1.2.0 (Eberhardt et al. 2021) via PyRx 0.9.2 (Dallakyan and Olson 2015) with default parameters. As ligands, we used α-L-fucosyl-L-threonyl residue and O-(3-O-D-galactosyl-N-acetyl-β-D-galactosaminyl)-L-serine.
In order to verify the expression of our predicted proteins, we made use of publicly available transcriptomic datasets of different species of early-diverging fungi (EDFs). The data found were whole transcriptome datasets obtained from pure culture as well from stress conditions. They were collected from the ENA server as fastq files. They were quality checked using FastQC (v0.11.8) (Andrews and Others 2010), followed by trimming of adapters using fastp (v0.19.6) (S. Chen et al. 2018) with default parameters, respectively. The reference genomes of the EDF species were retrieved from Ensembl Fungi (release 53) and aligned with the trimmed reads using Hisat2 (Sirén, Välimäki and Mäkinen 2014). The SAM files obtained from alignment were compressed into binary file format using samtools (v1.10). Finally, the aligned reads were counted for their expression using featureCounts (v1.6.3) (Liao et al. 2014) taking their respective annotations in the form of GTF files from Ensembl (Yates et al. 2022).
The TPM (transcripts per million reads) values and coefficient of variation of TPM were calculated to determine protein expression from whole transcriptome data (Stanton et al. 2017), while differential gene expression was carried out using DESeq2 R package (Love et al. 2014) for protein-coding gene expression from stress transcriptomic datasets. The reads were then mapped with fucose-related protein-coding genes. The log2 fold change criteria of [downregulation < 0 > upregulation] were used to determine their expression profiles during stress and coefficient of variation was used as a determinant of expression from whole transcriptome analysis.
To construct the fungal fucose metabolic map and to uncover possible fucose-related functions in fungal cells, we analyzed known eukaryotic and prokaryotic pathways. We used 75 protein families—38 enzymes, 18 transporters and 19 lectins from the whole tree of life as references (Additional file 3: Table S1, Table 1). Out of the aforementioned 15 enzymes, 2 groups of transporters and 9 lectins were confirmed to be present in fungal protein sequence sets (Fig. 1). None of the bacterial fucose-related proteins had fungal homologs. Most of the identified fungal homologs belonged to protein families previously described in metazoa. Only a few were reported from fungi, for instance, lectins, fucose synthase and GT10 fucosyltransferases. Seven main functions in which the studied proteins are involved can be distinguished: 1. Nutrition; 2. Fucose biosynthesis; 3. Fucose recycling; 4. Fucosylation; 5. Fucose degradation; 6. Fucose binding (lectins); 7. Transport of fucose. Here, we navigate the fucose-related metabolic map based on the protein functions.
Fucose as a simple hexose can be an easy source of energy. The use of fucose as a nutrient is possible due to degradation of complex compounds by ɑ-L-fucosidase (EC: 184.108.40.206 and 220.127.116.11) (Paper et al. 2013). α-L-fucosidases belong to either CAZy glycoside hydrolases family GH29 and are characterized by domains Alpha_L_fucos (PF01120) and Fucosidase_C (PF16757) or they belong to family GH95 and have a Glyco_hyd_65N_2 (PF14498) domain. Fucosidases from family GH29 act on fucosyl linkages α-1,2/3/4/6, while GH95 homologs act only on α-1,2 linkage (Paper et al. 2013). Identified proteins from both families are predicted to have extracellular localization.
The ability to treat fucose as a nutrient inferred from the presence of ɑ-L-fucosidase coding genes seems to be limited to Dikarya. We found homologs of this enzyme only in the representatives of Agaricomycotina, Ustilagomycotina (Basidiomycota) and Pezizomycotina (Ascomycota). Notably many of the taxa encoding ɑ-L-fucosidase genes rely on plant organic matter either as pathogens of saprotrophs or belong to lineages evolutionarily linked to plant hosts. We found only one homolog of this family in Mucoromycotina which suggests the possibility of contamination.
Fucose biosynthesis and recycling
Fucose can be a building block of complex carbohydrates and be added to diverse molecules. As such, it needs to be acquired or synthesized in the fungal cell. GDP-L-fucose is an activated nucleotide form of L-fucose and is required for fucosylation carried out by fucosyltransferases. GDP-L-fucose may be synthesized in two ways: via de-novo and salvage pathways (Fig. 2). Components of the salvage pathway (also of the fucose recycling process)—fucokinase FUK (EC: 18.104.22.168) and fucose-1-phosphate guanylyltransferase GFPP (EC: 22.214.171.124) are present only in Mucoromycota—we found only single homologs of FUK and GFPP outside this taxon in some chytrids. The enzymes involved in the de novo synthesis pathway (GDP-mannose 4,6-dehydratase (EC: 126.96.36.199) and GDP-L-fucose synthase (EC: 188.8.131.521)) are present in all sequenced Mortierellomycotina but not limited to this taxon. Representatives of Dikarya, Zoopagomycota, Chytridiomycota and Blastocladiomycota can probably synthesize GDP-L-fucose de novo. However, the conservation of this pathway is incomplete. Several Ascomycota (Pezizomycotina, n = 11) only have the second enzyme in the pathway the GDP-L-fucose synthase, whereas Ustilagomycotina have GDP-mannose 4,6-dehydratase alone.
Fungi possess five families of enzymes involved in fucose-related modification of proteins, lipids and polysaccharides (Fig. 3). These fucosyltransferases are homologous to animal transferases (Fig. 4). These enzymes belong to two Pfam clans GT-A and GT-B, both possessing a Rossman fold and grouping diverse glycosyltransferases. O-fucosylpeptide 3-beta-N-acetylglucosaminyltransferase (LFNG) (EC: 184.108.40.206) defined by Fringe domain (PF02434) is a transferase belonging to Pfam clan GT-A (CL0110). This protein is not able to fucosylate, but is interesting for this work as it operates on fucosyl peptides.
Fringe proteins act on fucosylated and galactosylated peptides (Ju and Cummings 2002). These two groups are highly similar in terms of sequence and structure. We were unable to recognize homologs of Fringe-domain containing proteins that act exclusively on fucosyl peptides (neither based on sequence and structure nor by molecular docking of simplified ligands, see Additional file 1: for clustering results).
The remaining four enzymes are representatives of the GT-B clan (CL0113) (Fig. 5). Homologs of CAZY glycosyl transferase 10 family (GT10, EC: 220.127.116.11 / 18.104.22.168 / 22.214.171.124) have a very similar distribution to proteins with a XG_FTase (PF03254) domain. Both families are present in Pucciniomycotina, Agaricomycotina, Mucoromycota and Chytridiomycota. Blastocladiomycota fungi seem to have only GT10 homologs in terms of fucosyltransferases in general. Interestingly, Basidiobolus is the only Zoopagomycota lineage encoding a rich repertoire of fucosyltransferase genes (GT10, XG_FTase and O-FucT). Proteins with O-FucT domain (PF10250) (peptide-O-fucosyltransferase, EC: 126.96.36.199) make up the most diverse group of fucosyltransferases. We recognize eight subfamilies of these transferases mapping on the Pfam O-FucT domain but with distinguishable sequence motifs and taxonomic distribution (Figs. 5 and 6). Two subfamilies are plant-specific and the remaining ones are fungal subfamilies with the highest occurrence in Agaricomycotina, Pucciniomycotina and Mortierellomycotina. Representatives of O-Fuct are present in the earliest fungal lineages including Rozella and Olpidium. Still, they show a repeated loss of this trait in yeast-forming fungi (Ustilagomycotina, Saccharomycotina, and just a single copy in Taphrinomycotina). The characteristic sequence motifs of each new O-FucT subfamily corresponding to conserved sequences of POFUTs fucosyltransferases (Chen et al. 2012) are summarized with alignment logos (Fig. 6). Almost all novel subfamilies have well-preserved key POFUTs residues; we analyzed the conservation of catalytic amino acids residues—E54 and R294, as well as conservation of amino acids critical for GDP-fucose binding—P53, D371, S387, T388 and F389 (Chen et al. 2012). R294 residue is well preserved among all identified POFUTs. In the case of E54, this residue is conserved only in C and G, but the surroundings keep a similar, aromatic character in most POFUTs—in A, B, F. P53 have a very similar distribution. In the case of D and E, we can observe well-conserved G55 residue, but without following aromatic residues. D371 and S387 are strongly conserved among all groups. T388 and F389 are not very well conserved but are substituted by amino acids with similar properties.
Subfamily G encompasses characterized fucosyltransferases involved in strobilurin polyketide modification (Nofiani et al. 2018), cell-cycle regulation coupled with cell signaling and differentiation. O-FucT subfamilies present in Dikarya A and B are surrounded by secondary and cellular metabolism genes acting as a regulator of autophagy with its role in the formation of autosomes and in vacuolar protein sorting. The subfamilies E and F are also involved in secondary metabolism since the subfamily E genes co-occur with polyketide synthetase, ketopantoate reductase as well as MFS1 class transporters found ubiquitously in SMC. Genes encoding fucosyltransferases from subfamilies A and B occur mostly in a housekeeping neighborhood related to transcription, DNA repair and lipid metabolism. Subfamily C present mostly in EDFs may have an intracellular role linked to dolichol metabolism due to the occurrence of decaprenyl diphosphate synthase-like and undecaprenyl diphosphate synthase domains.
We did not find fungal homologs of animal FUT1 and FUT2 transferases, which belong to the CAZY GT11 family. FUT1,2 catalyzes the transfer of fucose to the terminal galactose through ɑ-1,2 linkage (Kudo and Narimatsu 2008).
Enzymes with FUT8_N_cat domain (PF19745) (glycoprotein 6-alpha-L-fucosyltransferase, EC: 188.8.131.52) are mostly present in Mortierellomycotina and Chytridiomycotina. In Dikarya, they are present in Pezizomycotina (only a few homologs), Pucciniomycotina, Taphrinomycotina and Saccharomycotina. The Dikarya homologs are mostly present in entomopathogenic fungi and form a separate clade in sequence clustering (Fig. 5). Mortierellomycotina sequences group into a clade distinct from the Dikarya.
RNA-Seq data analysis identified the expression of four O-FucT subfamilies. While Mucoromycotina and Glomeromycotina expressed protein-coding genes from subfamilies C and E, Mortierellomycotina expressed subfamilies C, F and G from the whole transcriptome data. The analysis of stress response data displayed up and downregulation of fucose protein-coding genes in Glomeromycotina. All of the genes expressed in Mucoromycotina were upregulated. We did not find any publicly available transcriptomic stress experiments in organisms belonging to Mortierellomycotina.
In fungi—in addition to alpha-L-fucosidase GH29 involved in nutrition—three enzymes are responsible for the degradation of fucose and fucose-containing polysaccharides (Fig. 7). The most numerous is the family of L-fuconate dehydratase (EC: 184.108.40.206). These proteins are built from the domains MR_MLE_C (PF13378) and MR_MLE_N (PF02746). They are present in Agaricomycotina, Pezizomycotina, Mortierellomycotina and Mucoromycotina. The same taxonomic distribution was observed for L-fucose mutarotase (EC: 220.127.116.11), defined by the RbsD_FucU (PF05025) domain. 1,2-alpha-L-fucosidase (EC: 18.104.22.168) homologs play a role in the degradation of polysaccharides by cleaving fucose residue. They have Glyco_hyd_65N_2 (PF14498) domain and occur in Ustilagomycotina, Agaricomycotina, Pezizomycotina. A single sequence was also found in Umbelopsis vinacea. The relatively narrow distribution of both enzymes opens the question of how other fucose-containing lineages degrade this substrate.
The MetaCyc database describes two L-fucose degradation pathways (named I and II). Fungi probably do not have L-fucose degradation pathway I as we did not detect L-fucose isomerase. We found only two homologs of this enzyme in Hortaea werneckii and Mortierella sp. GBA39. These sequences show extreme similarity to bacterial sequences from the Paenibacillus genus that point to contamination of fungal genomic sequences. The function of other enzymes of pathway—L-fuculokinase and L-fuculose-phosphate aldolase—may be handled by fungal homologs with the same domain architectures described as xylulose kinase and unannotated aldolases, but the question remains if these enzymes have appropriate substrate specificity. In the case of degradation pathway II, we do not find homologs of specific enzymes described as pathways components in the MetaCyc database, but we speculate that the function of these enzymes may be handled by other aldolases, hydrolases and dehydrogenases. For example, TrAA12 dehydrogenase from Trichoderma reesei shows weak activity against uncommon carbohydrates including L-fucose (Turbe-Doan et al. 2019).
Fucose binding (lectins)
We identified nine families of proteins that can bind compounds with fucose or bind fucose directly. Fucose-binding lectins are most abundant in Dikarya, particularly in Pezizomycotina and Agaricomycotina. We found homologs of previously described fungal lectins like Aleuria aurantia lectin (AAL) group (Fungal_lectin (PF07938)) (Wimmerova et al. 2003), galectins (Gal-bind_lectin (PF00337)), e.g., Coprinus cinereus protein CGL2 (Cooper et al. 1997), Pholiota squarrosa lectin (PhoSL) (also called RSL when from Rhizopus stolonifer (Oda et al. 2003)), tectonins (PF19193), e.g., LbTec2 from Laccaria bicolor (Sommer et al. 2018), Agglutinin_C (PF14200) and RicinB_lectin_2 (PF14200) domains containing proteins like MOA lectins (M. oreades agglutinin) (Grahn et al. 2009) and RicinB_like family (CCL2 protein from C. cinerea. CCL2 (Schubert et al. 2012).
The other lectins found are known to date from other organisms but were not described in fungi—DUF346 (PF03984) lectins, PA-IIL (PF07472) and Lectin_legB (PF00139). Outside Dikarya, we found lectins in Glomeromycotina and Blastocladiomycota (Tectonin), Mucoromycotina (PhoSL and Lectin_legB) and Chytridiomycota (Lectin_legB).
Transport of fucose
TCDB classification was used to determine candidates for fucose transporters. All identified fungal proteins belong to 2.A.7 family, subfamilies 2.A.7.10.4 and 2.A.7.16. Fungal 2.A.7.10.4 family transporters form three groups, each with a different taxonomic distribution. The first group consists of Basidiomycota and EDF sequences, the second group includes only Ascomycota sequences and the third group is Basidiomycota specific. In the case of 2.A.7.16 fungal transporters, we also observed three groups. Two of them are specific for EDF (one consists only of Mortierellomycotina and Glomeromycotina homologs) and one group for both EDF and Dikarya sequences.
Fucose-related traits not confirmed in fungi
L-fucose is the main component of fucoidan—a complex polymer with sulfation and acetylation modifications produced mostly by brown algae (Li et al. 2008). We did not find homologs of enzymes capable of synthesizing fucosyl polymers and in consequence have no intermediate evidence for their existence in fungi. For instance, we have not found homologs of enzymes involved in fucose polysaccharide synthesis like TDP-fucosamine acetyltransferase, fucoidan-2-sulfatase and UDP-N-acetyl-L-fucosamine synthase. Fucose polymer degradation, if existent, has to be performed without homologs of bacterial fucoidanases. We did not identify any of the prokaryotic components of the fucose metabolism. Suppose fungi are able to synthesize secondary metabolites with fucosyl residues. In that case, the synthesis has to be independent of dTDP-fucopyranose mutase, aclacinomycin-T 2-deoxy-L-fucose transferase and dTDP-alpha-D-fucopyranose:NAD(P) + oxidoreductase. Bacterial antigen biosynthesis was also missing based on the absence of TDP-N-acetylfucosamine: lipid II N-acetylfucosaminyltransferase and GDP-beta-L-fucose: beta-D-Gal-(1- > 3)-alpha-D-GalNAc-(1- > 3)-alpha-D-GalNAc-diphospho-ditrans, octa-cis-undecaprenol alpha-1,2-fucosyltransferase.
Our aim in this work was to provide a generalized picture of fucose usage and metabolism on the fungal tree of life. To design the metabolic map, we have compared existing prokaryotic and eukaryotic proteins depending on or using fucose residues in any form. We have then detected their presence in diverse fungal lineages and reconstructed metabolic maps for the main fungal lineages taking into account an ecological and biochemical perspective. Most of the identified protein families are homologs of animal fucose-related proteins and none of the bacterial protein families was found in fungi. This led us to predict different patterns of fucose usage in EDF compared to Dikarya. The former seems to predominantly use fucose either as a building block for polymers or as a residue for macromolecule modification. The latter detects fucose in their hosts with lectins and possesses a set of putative enzymates that allows to use it as an energy source. Noteworthy, yeast-forming lineages of both Ascomycota (Saccharomycotina) and Basidiomycota (Ustilagomycotina) have lost most of the fucose-related protein-coding genes. Overall the evolutionary history of fucose metabolism in Fungi seems to be shaped by consecutive and parallel losses of its components, leading to its simplification or even complete loss. The overall fucose metabolic capabilities of fungi and animals remain similar. One of the fucosyltransferase families GT11 is apparently lost in all Fungi.
Most of the fucose-related proteins form a coherent group in terms of sequence diversity, confirming a common origin predating the radiation of animals and fungi. Despite the ubiquity of endohyphal bacteria among fungi, and documented co-evolution with some Mucoromycota taxa, we observed no transfer of fucose metabolic genes to fungal genomes. We have, however, detected a burst of diversity in the O-FucT fucosylase family with changes both in the core sequence motifs and occurrence in fungal lineages.
The complexity of the fucose map genes varies greatly between fungal lineages, with Mortierellomycotina possessing the most complete map among studied taxa. They are the only lineage to encode both fucose synthesis and recycling pathways as well as a wide range of fucosyltransferases including the newly described O-FucT subfamilies. In the case of Mucoromycotina, some enzymes (fucokinase, fucose-1-phosphate guanylyltransferase) are limited only to Endogonales. On the other hand, all yeast-forming lineages do not use fucose as a building block and hardly ever have fucose-detecting lectins. Most Dikarya lineages with a complex fucose map are related to a plant host either as saprotrophs or pathogens.
In our previous work, we showed that two Umbelopsis species have up to 25% of fucose in their cell wall (Muszewska et al. 2021). Mucoran is a polysaccharide described in the cell wall of Mucoromycota. It is composed of D-glucuronic acid, L-fucose and D-mannose in a molar ratio of 5:2:3 (Bartnicki-Garcia and Lindberg 1972). The composition of the Umbelopsis spp. cell wall differs from this ratio. However, we did not find in their genomes homologs of enzymes involved in the biosynthesis of fucose polymers. A question arises about the form of fucose present in Umbelopsis cell wall. Some of the remaining fucose may exist as a part of glycuronan—polysaccharide present in Mucorales and it is composed of fucose, mannose, galactose and glucuronic acid with a molar ratio of 5:1:1:6 (Lecointe et al. 2019). Mucoran can also be an acceptor of fucose residues (Camacho-Aguero, Balcazar-Orozco, and Flores-Carreo´n 1990). Furthermore, Mucorales have also exopolysaccharides (EPS) that contain fucose (De Ruiter et al. 1992) so maybe this is also the case for Umbelopsis spp. One might speculate that the detection of high amounts of fucose in the cell wall might be also attributed to the presence of EPS. It is also likely that fucose is localized in both EPS and the cell wall.
The presence of both GDP-L-fucose synthesis pathways—the de novo and the salvage pathway in M. alpina—was reported already in 2010 (Ren et al. 2010; Wang et al. 2016). According to our results, both pathways co-occur only in Mortierellomycotina. Our results suggest that the salvage pathway is limited to Mucoromycota. As the functional salvage pathway is present in animals (including mammals) (Schneider et al. 2017), we assume that genes encoding components of this pathway were lost in the evolution of fungi outside Mucoromycota. An additional argument supporting this thesis is the presence of the salvage pathway in the tree of life before the formation of Opisthokonta—its occurrence was described in Dictyostelium discoideum (Gonzalez-Yanes et al. 1989).
Enzymes involved in the de novo synthesis pathway are present in all fungal divisions. However, several lineages encode only one of the two enzymes needed which opens the question on the possibility to produce GDP-L-fucose without GDP-mannose 4,6-dehydratase and function replacement by another dehydratase. The presence of both pathways in Mortierellomycotina points to a peculiar complexity of fucose metabolism in this fungal subphylum.
More than 13 different fucosyltransferases were described in humans (FUTs and POFUTs). They play a role in posttranslational modification of proteins, glycans and lipids (Schneider et al. 2017; Shimizu et al. 2001). We identify five families of fungal fucosyltransferases, all present both in Dikarya and EDF. It suggests that fucosylation is an ancient trait common in Opisthokonta. The ancestral set of fucosyltransferases comprises.
CAZY GT10 family (with domains Glyco_transf_10 (PF00852) and Glyco_tran_10_N (PF17039)) homologs of human FUT3-7 and 9–11 involved in the synthesis of Lewis antigens, homologs of plant cell wall-related FUT1 protein (Perrin et al. 1999) characterized by XG_FTase domain (PF03254), homologs of O-FucT called POFUT1 and POFUT2—involved in Notch signaling (Stahl et al. 2008; Sparrow et al. 2006)). The Fungal O-FucT family harbors a few novel subfamilies with derived conservation of the active site motif suggesting a possible change of the protein substrate. This is particularly likely in the absence of EGF motives in fungal hyphae. One might speculate that Dikarya O-FucT subfamilies can be a part of candidate secondary metabolite clusters similar to the strobilurin polyketide. The subfamilies occurring in EDF might be involved in some housekeeping function or signaling based on the genomic neighborhood observations. The revealed complexity of fucosyltransferases is smaller than in Metazoa but greater than expected based on published data. The fucosyltransferase repertoire may work on a wide range of substrates and contribute to novel natural compounds, molecular biology tools or diagnostic markers.
The fungal transferases catalyze the transfer of L-fucose from GDP-L-fucose to the N-acetyl glucosamine via linkage ɑ-1,3 (FUT4,6,7,9) or ɑ-1,4 (FUT3,5) (Breton et al. 1998). POFUTs catalyze the reaction by attaching fucose to serine/threonine residues via O-glycosidic linkage (C.-I. Chen et al. 2012). We did not find fungal homologs of animal FUT1 and FUT2 transferases, which belong to CAZY GT11 family. FUT1,2 catalyzes the transfer of fucose to the terminal galactose through ɑ-1,2 linkage (Kudo and Narimatsu 2008). But we found homologs of plant FUT1 that operate on the same linkage (Sarria et al. 2001). We also found homologs of FUT8 that add fucose to N-glycans via ɑ-1,6 linkage (García-García et al. 2021). The possibility of transfer of N-acetyl-α-D-glucosaminyl residue to O-fucose-modified peptides/proteins by Fringe homologs opens up several questions regarding what could be the fucosylated peptide/protein and in which biological processes this mechanism would be involved. Taken together, fungi are able to process all kinds of fucosylation known in animals—it is especially true for Mortierellomycotina, as fungi from this taxa have homologs of all described transferases from GT-B clan.
In addition to glycoside hydrolase family 29 present in Dikarya and involved in nutrition, we found other hydrolases among Dikarya and EDF that are involved in the degradation and recycling of fucose within a cell. 1,2-alpha-L-fucosidase (EC: 22.214.171.124) is responsible for the degradation of polysaccharides to L-fucopyranose. This enzyme belongs to the GH95 family and was described before in Aspergillus niger (Bahl 1970).
MetaCyc database describes two L-fucose degradation pathways (named I and II). Fungi probably do not have L-fucose degradation pathway I as they lack L-fucose isomerase. We found only two homologs of this enzyme in Hortaea werneckii and Mortierella sp. GBA39. This sequence shows extreme similarity to bacterial sequences from the Paenibacillus genus, which point to contamination of fungal genomic sequences. The function of other enzymes of the pathway—L-fuculokinase and L-fuculose-phosphate aldolase—may be handled by fungal homologs with the same domain architectures described as xylulose kinase (Y. M. Chen et al. 1987) and unannotated aldolases, but the question remains if these enzymes have appropriate substrate specificity.
We hypothesize that fungi may have L-fucose degradation pathway II. However, we did not find homologs of L-fucose dehydrogenase, L-fuconolactonase, 2-dehydro-3-deoxy-L-fuconate-4-dehydrogenase and 2,4-didehydro-3-deoxy-L-fuconate hydrolase, but we speculate that function of these enzymes may be handled by different enzymes that possess the same protein domains. L-fucose dehydrogenase activity may be performed by broad specificity enzymes such as AA12 from T. reesei (Turbe-Doan et al. 2019).
Lectins are a diverse and widespread group of proteins, capable of binding oligosaccharides in a specific, non-catalytic way (Rini 1995), and are produced by all organisms (Kobayashi and Kawagishi 2014). In Fungi, these proteins are known to play a role in storage, growth, morphogenesis, parasitism, molecular recognition and mating (Kobayashi and Kawagishi 2014). Most of the described fungal lectins were from Dikarya (especially fruiting body forming). Fungal lectins are of particular biotechnological interest and numerous studies deciphered their structural properties allowing them to be exploited as a diagnostic tool for cancer or in targeted drug delivery (Varrot et al. 2013).
Lectin LecB (called also PA-IIL, same as domain PF07472) is one of two lectins synthesized by the opportunistic pathogen Pseudomonas aeruginosa and are important virulence factors. LecB plays a role, especially in Pseudomonas respiratory chronic infection (Tielker et al. 2005). This kind of lectins was never described in fungi before. We observed the occurrence of LecB homologs in a very narrow taxonomic context—they are present in two broad host range Clavicipitaceae enthomoparasitic genera, Metarhizium and Hirsutella. Obtained sequences show similarity to bacterial homologs, e.g., from Nostoc, Kitasatospora and Sorangium, which may suggest horizontal gene transfer from bacteria. However, as similarity ranks around 40%, and there is no ecological overlap between these bacteria and fungi, the evolutionary origin of these lectins is not clear.
We found homologs of fucose-binding lectin from Photorhabdus luminescens (PLL). These proteins are characterized by DUF346 domain (PF03984). P. luminescens is a bacteria known for its symbiotic relationship with nematodes from the family Heterorhabditidae and its pathogenicity toward insects (Forst et al. 1997). The biological role of PLL remains unclear. Some sources suggest that PLL has larvicidal activity (Eisemann et al. 1994), while others suggest that the role of PLL is limited to cell recognition and adhesion (Kumar et al. 2016).
Lectin_legB (PF00139) lectins are originally known from plants as legume lectins. Lotus tetragonolobus lectin (LTA) is a fucose-specific lectin with described diverse activities (Moreno et al. 2008), e.g., reactivity against blood group O (Konami et al. 1990). Legume lectins play a role in nodulation by Rhizobium (Hirsch 1999).
Almost twenty TCDB families of transporters are related to fucose transport. From them, only three subfamilies have homologs in fungi. Proteins from the drug/metabolite transporter (DMT) superfamily 2.A.7 were described among others as Efr and Gfr (ER (endoplasmic reticulum)/Golgi GDP-fucose transporter) in Drosophila (Ishikawa et al. 2010) and SLC35C1 in Homo sapiens (Ishida and Kawakita 2004). Transport of GDP-fucose into the ER and Golgi is necessary for fucosylation because GDP-fucose is synthesized in the cytosol but fucosylation takes place in the lumen of ER and Golgi (Ishikawa et al. 2010).
Also, ATP-binding cassette (ABC) superfamily transporters (3.A.1) are probably involved in fucose transport in fungi. Proteins from family 3.A.1.2.5 act like multiple sugar transporters with the ability to transport fucose among other sugars. The occurrence of such transporters was described in bacteria (Conners et al. 2005; Kemner et al. 1997). So far fucose transporters were not verified experimentally in fungi.
Here, we present a novel putative metabolic map in fungi with a distinction between non-Dikarya and Dikarya subtypes of fucose utilization. The former likely depends on fucose as a building and signaling compound while the latter mainly detects it and uses it as a carbon source. We have observed some clear patterns with most of the fucose-related proteins co-occurring within Opisthokonta and no HGT bacterial protein families present in the analyzed fungi. Fungi possess a limited subset of the metazoa fucose metabolic map, however, sufficient to be able to synthesize, degrade, bind, transfer and transport this compound. This ensemble provides a list of hypotheses to be tested experimentally.
Availability of data and materials
All accessions and assemblies are listed in Additional file 3: Table S1.
Fungal tree of life
Aleuria aurantia lectin
Transcripts per million reads
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215(3):403–410
Amano K, Katayama H, Saito A, Ando A, Nagata Y (2012) Aleuria aurantia lectin exhibits antifungal activity against Mucor racemosus. Biosci Biotechnol Biochem 76(5):967–970
Andrews S (2010) FastQC: a quality control tool for high throughput sequence data. Babraham Bioinformatics, Babraham Institute, Cambridge, United Kingdom
Armenteros A, Juan J, Tsirigos KD, Sønderby CK, Petersen TN, Winther O, Brunak S, von Heijne G, Nielsen H (2019) SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat Biotechnol 37(4):420–423
Bahl OP (1970) Glycosidases of Aspergillus niger. J Biol Chem 245(2):299–304
Bartnicki-Garcia S (1968) Cell wall chemistry, morphogenesis, and taxonomy of fungi. Annu Rev Microbiol 22:87–108
Bartnicki-Garcia S, Lindberg B (1972) Partial characterization of mucoran: the glucuronomannan component. Carbohyd Res 23(1):75–85
Bonnardel F, Mariethoz J, Salentin S, Robin X, Schroeder M, Perez S, Lisacek F, Imberty A (2019) UniLectin3D, a database of carbohydrate binding proteins with curated information on 3D structures and interacting ligands. Nucleic Acids Res 47(D1):D1236–D1244
Breton C, Oriol R, Imberty A (1998) Conserved structural features in eukaryotic and prokaryotic fucosyltransferases. Glycobiology 8(1):87–94
Camacho-Aguero S, Balcazar-Orozco R, Flores-Carreo A (1990) Biosynthesis of polyuronides in Mucor rouxii: partial characterization of fucosyl transferase. Exper Mycol 14(3):227–233
Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T (2009) trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25(15):1972-1973
Caspi R, Billington R, Ferrer L, Foerster H, Fulcher CA, Keseler IM, Kothari A et al (2016) The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res 44(D1):D471–D480
Chang A, Jeske L, Ulbrich S, Hofmann J, Koblitz J, Schomburg I, Neumann-Schaal M, Jahn D, Schomburg D (2021) BRENDA, the ELIXIR core data resource in 2021: new developments and updates. Nucleic Acids Res 49(D1):D498-508
Chen YM, Zhu Y, Lin EC (1987) The organization of the fuc regulon specifying L-fucose dissimilation in Escherichia coli K12 as determined by gene cloning. Mol Gen Genet MGG 210(2):331–337
Chen C-I, Keusch JJ, Klein D, Hess D, Hofsteenge J, Gut H (2012) Structure of human POFUT2: insights into thrombospondin type 1 repeat fold and O-fucosylation. EMBO J 31(14):3183–3197
Chen S, Zhou Y, Chen Y, Jia Gu (2018) Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34(17):i884–i890
Church DM, Schneider VA, Graves T, Auger K, Cunningham F, Bouk N, Chen H-C et al (2011) Modernizing reference genome assemblies. PLoS Biol 9(7):e1001091
Conners SB, Montero CI, Comfort DA, Shockley KR, Johnson MR, Chhabra SR, Kelly RM (2005) An expression-driven approach to the prediction of carbohydrate transport and utilization regulons in the hyperthermophilic bacterium Thermotoga maritima. J Bacteriol 187(21):7267–7282
Cooper DN, Boulianne RP, Charlton S, Farrell EM, Sucher A, Lu BC (1997) Fungal galectins, sequence and specificity of two isolectins from Coprinus cinereus. J Biol Chem 272(3):1514–1521
Crooks GE, Hon G, Chandonia J-M, Brenner SE (2004) WebLogo: a sequence logo generator. Genome Res. https://doi.org/10.1101/gr.849004
Dallakyan S, Olson AJ (2015) Small-molecule library screening by docking with PyRx. Methods Mol Biol 1263:243–250
Davis WJ, Amses KR, Benny GL, Carter-House D, Chang Y, Grigoriev I, Smith ME, Spatafora JW, Stajich JE, James TY (2019) Genome-scale phylogenetics reveals a monophyletic Zoopagales (Zoopagomycota, Fungi). Molec Phylogenet Evol 133:152-163
Du J, Takeuchi H, Leonhard-Melief C, Shroyer KR, Dlugosz M, Haltiwanger RS, Holdener BC (2010) O-fucosylation of thrombospondin type 1 repeats restricts epithelial to mesenchymal transition (EMT) and maintains epiblast pluripotency during mouse gastrulation. Dev Biol 346(1):25–38
Eberhardt J, Santos-Martins D, Tillack A, Forli S (2021) AutoDock Vina 1.2.0: new docking methods, expanded force field, and python bindings. https://doi.org/10.26434/chemrxiv.14774223.
Eisemann CH, Donaldson RA, Pearson RD, Cadogan LC, Vuocolo T, Tellam RL (1994) Larvicidal activity of lectins on Lucilia cuprina: mechanism of action. Entomol Exp Appl 72(1):1–10
Fairclough SR, Chen Z, Kramer E, Zeng Q, Young S, Robertson HM, Begovic E et al (2013) Premetazoan genome evolution and the regulation of cell differentiation in the choanoflagellate Salpingoeca rosetta. Genome Biol 14(2):R15
Forst S, Dowds B, Boemare N, Stackebrandt E (1997) Xenorhabdus and Photorhabdus spp.: bugs that kill bugs. Annu Rev Microbiol 51:47–72
Frickey T, Lupas A (2004) CLANS: a java application for visualizing protein families based on pairwise similarity. Bioinformatics 20(18):3702–3704
Garber JM, Hennet T, Szymanski CM (2021) Significance of fucose in intestinal health and disease. Mol Microbiol 115(6):1086–1093
García-García A, Serna S, Yang Z, Delso I, Taleb V, Hicks T, Artschwager R et al (2021) FUT8-directed core fucosylation of N-glycans is regulated by the glycan structure and protein environment. ACS Catal 11(15):9052–9065
Gonzalez-Yanes B, Mandell RB, Girard M, Henry S, Aparicio O, Gritzali M, Brown RD Jr, Erdos GW, West CM (1989) The spore coat of a fucosylation mutant in Dictyostelium discoideum. Dev Biol 133(2):576–587
Grahn EM, Winter HC, Tateno H, Goldstein IJ, Krengel U (2009) Structural characterization of a lectin from the mushroom Marasmius oreades in complex with the blood group B trisaccharide and calcium. J Mol Biol 390(3):457–466
Grass J, Pabst M, Kolarich D, Pöltl G, Léonard R, Brecker L, Altmann F (2011) Discovery and structural characterization of Fucosylated oligomannosidic n-glycans in mushrooms. J Biol Chem 286(8):5977–5984
Hirsch AM (1999) Role of lectins (and rhizobial exopolysaccharides) in legume nodulation. Curr Opin Plant Biol 2(4):320–326
Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K (2007) WoLF PSORT: protein localization predictor. Nucleic Acids Res 35:W585–W587
Inkscape Project. 2020. Inkscape (version 1.2). https://inkscape.org.
Ishida N, Kawakita M (2004) Molecular physiology and pathology of the nucleotide sugar transporter family (SLC35). Pflugers Arch 447(5):768–775
Ishikawa HO, Ayukawa T, Nakayama M, Higashi S, Kamiyama S, Nishihara S, Aoki K, Ishida N, Sanai Y, Matsuno K (2010) Two pathways for importing GDP-fucose into the endoplasmic reticulum lumen function redundantly in the O-fucosylation of notch in Drosophila. J Biol Chem 285(6):4122–4129
Ishio A, Sasamura T, Ayukawa T, Kuroda J, Ishikawa HO, Aoyama N, Matsumoto K et al (2015) O-fucose monosaccharide of drosophila notch has a temperature-sensitive function and cooperates with O-glucose glycan in notch transport and notch signaling activation. J Biol Chem 290(1):505–519. https://doi.org/10.1074/jbc.M114.616847
Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, McWilliam H et al (2014) InterProScan 5: genome-scale protein function classification. Bioinformatics 30(9):1236–1240
Ju T, Cummings RD (2002) A unique molecular chaperone cosmc required for activity of the mammalian core 1 beta 3-galactosyltransferase. Proc Natl Acad Sci USA 99(26):16613–16618
Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, Tunyasuvunakool K et al (2021) Highly accurate protein structure prediction with alphafold. Nature 596(7873):583–589
Katoh K, Misawa K, Kuma KI, Miyata T (2002) MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucl Acids Res 30(14):3059–3066
Kemner JM, Liang X, Nester EW (1997) The Agrobacterium tumefaciens virulence gene chvE is part of a putative ABC-type sugar transport operon. J Bacteriol 179(7):2452–2458
King N, Jody Westbrook M, Young SL, Kuo A, Abedin M, Chapman J, Fairclough S et al (2008) The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans. Nature 451(7180):783–788
Kizuka Y, Nakano M, Yamaguchi Y, Nakajima K, Oka R, Sato K, Ren C-T, Hsu T-L, Wong C-H, Taniguchi N (2017) An alkynyl-fucose halts hepatoma cell migration and invasion by inhibiting GDP-fucose-synthesizing enzyme FX, TSTA3. Cell Chem Biol 24(12):1467–78.e5
Kobayashi Y, Kawagishi H (2014) Fungal lectins: a growing family. Methods Mol Biol 1200:15–38
Konami Y, Yamamoto K, Osawa T (1990) The primary structure of the lotus tetragonolobus seed lectin. FEBS Lett 268(1):281–286
Kudo T, Narimatsu H (2008) Fucosyltransferase (α1,2/ α1,3/ α1,4-fucosyltransferases). In: Taniguchi N, Suzuki A, Ito Y, Narimatsu H, Kawasaki T, Hase S (eds) Experimental glycoscience: glycobiology. Springer Japan, Tokyo, pp 36–38
Kumar A, Sýkorová P, Demo G, Dobeš P, Hyršl P, Wimmerová M (2016) A novel fucose-binding lectin from Photorhabdus luminescens (PLL) with an unusual heptabladed β-propeller tetrameric structure. J Biol Chem 291(48):25032–25049
Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K et al (2001) Initial sequencing and analysis of the human genome. Nature 409(6822):860–921
Lebreton A, Bonnardel F, Dai Y-C, Imberty A, Martin FM, Lisacek F (2021) A comprehensive phylogenetic and bioinformatics survey of lectins in the fungal kingdom. J Fungi 7(6):453. https://doi.org/10.3390/jof7060453
Lecointe K, Cornu M, Leroy J, Coulon P, Sendid B (2019) Polysaccharides cell wall architecture of mucorales. Front Microbiol 10(March):469
Li Bo, Fei Lu, Wei X, Zhao R (2008) Fucoidan: structure and bioactivity. Molecules 13(8):1671–1695
Liao Y, Smyth GK, Shi W (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30(7):923–930
Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome Biol 15(12):550
Ma B, Simala-Grant JL, Taylor DE (2006) Fucosylation in prokaryotes and eukaryotes. Glycobiology 16(12):158R-184R
Martínez-Duncker I, Díaz-Jímenez DF, Mora-Montes HM (2014) Comparative analysis of protein glycosylation pathways in humans and the fungal pathogen candida albicans. Int J Microbiol 2014:267497
Matthews BB, Dos Santos G, Crosby MA, Emmert DB, St. Pierre SE, Gramates LS, Zhou P, Schroeder AJ, Falls K, Strelets V, Russo SM (2015) Gene model annotations for Drosophila melanogaster: impact of high-throughput data. G3: Genes Genomes Genet 5(8):1721–1736
Mélida H, Sain D, Stajich JE, Bulone V (2015) Deciphering the uniqueness of mucoromycotina cell walls by combining biochemical and phylogenomic approaches. Environ Microbiol 17(5):1649–1662
Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, Von Haeseler A, Lanfear R (2020) IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Molec Biol Evol 37(5):1530–1534
Mirdita M, Schütze K, Moriwaki Y, Heo L, Ovchinnikov S, Steinegger M (2022) ColabFold: making protein folding accessible to all. Nat Methods 19(6):679–682
Mistry J, Bateman A, Finn RD (2007) Predicting active site residue annotations in the pfam database. BMC Bioinform 8(August):298
Mistry J, Chuguransky S, Williams L, Qureshi M, Salazar GA, Sonnhammer ELL, Tosatto SCE et al (2021) Pfam: the protein families database in 2021. Nucleic Acids Res 49(D1):D412–D419
Miyazaki T, Irino T (1972) Studies on fungal polysaccharides. X. extracellular heteroglycans of Absidia cylindrospora and Mucor mucedo. Chem Pharm Bull 20(2):330–335
Moran AP (2008) Relevance of fucosylation and lewis antigen expression in the bacterial gastroduodenal pathogen Helicobacter pylori. Carbohyd Res 343(12):1952–1965
Moreno FB, Batista M, Maia T, de Oliveira D, Martil E, Viçoti MM, Bezerra GA, Abrego JRB, Cavada BS, Filgueira W, de Azevedo F Jr (2008) Identification of a new quaternary association for legume lectins. J Struct Biol 161(2):133–143
Muszewska A, Okrasińska A, Steczkiewicz K, Drgas O, Orłowska M, Perlińska-Lenart U, Aleksandrzak-Piekarczyk T et al (2021) Metabolic potential, ecology and presence of associated bacteria is reflected in genomic diversity of mucoromycotina. Front Microbiol 12:239
Nofiani R, de Mattos-Shipley K, Lebe KE, Han L-C, Iqbal Z, Bailey AM, Willis CL, Simpson TJ, Cox RJ (2018) Strobilurin biosynthesis in basidiomycete fungi. Nat Commun 9(1):3940
O’Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR (2011) Open babel: an open chemical toolbox. J Cheminformtics 3:33
Oda Y, Senaha T, Matsuno Y, Nakajima K, Naka R, Kinoshita M, Honda E, Furuta I, Kakehi K (2003) A new fungal lectin recognizing alpha(1–6)-linked fucose in the N-glycan. J Biol Chem 278(34):32439–32447
Ogawa S, Otta Y, Ando A, Nagata Y (2001) A lectin from an ascomycete mushroom, Melastiza chateri: no synthesis of the lectin in mycelial isolate. Biosci Biotechnol Biochem 65(3):686–689
Orłowska M, Muszewska A (2022) In silico predictions of ecological plasticity mediated by protein family expansions in early-diverging fungi. J Fungi 8(1):67. https://doi.org/10.3390/jof8010067
Paper JM, Scott-Craig JS, Cavalier D, Faik A, Wiemels RE, Borrusch MS, Bongers M, Walton JD (2013) α-Fucosidases with different substrate specificities from two species of Fusarium. Appl Microbiol Biotechnol 97(12):5371–5380
Perrin RM, DeRocher AE, Bar-Peled M, Zeng W, Norambuena L, Orellana A, Raikhel NV, Keegstra K (1999) Xyloglucan fucosyltransferase, an enzyme involved in plant cell wall biosynthesis. Science 284(5422):1976–1979
Ren Y, Perepelov AV, Wang H, Zhang H, Knirel YA, Wang L, Chen W (2010) Biochemical characterization of Gdp-L-Fucose de novo synthesis pathway in fungus Mortierella alpina. Biochem Biophys Res Commun 391(4):1663–1669
Rini JM (1995) Lectin structure. Annu Rev Biophys Biomol Struct 24:551–577
Ruiter De, Gerhard A, Josso SL, Colquhoun IJ, Voragen AGJ, Rombouts FM (1992) Isolation and characterization of β(1–4)-D-glucuronans from extracellular polysaccharides of moulds belonging to mucorales. Carbohyd Polym 18(1):1–7
Ruiz-Trillo I, Burger G, Holland PWH, Nicole King B, Lang F, Roger AJ, Gray MW (2007) The origins of multicellularity: a multi-taxon genome initiative. Trends Genet TIG 23(3):113–118
Saier MH, Reddy VS, Moreno-Hagelsieb G, Hendargo KJ, Zhang Y, Iddamsetty V, Lam KJK et al (2021) The transporter classification database (TCDB): 2021 update. Nucleic Acids Res 49(D1):D461–D467
Sakai K, Hiemori K, Tateno H, Hirabayashi J, Gonoi T (2019) Fucose-specific lectin of Aspergillus fumigatus: binding properties and effects on immune response stimulation. Med Mycol 57(1):71–83
Sarria R, Wagner TA, O’Neill MA, Faik A, Wilkerson CG, Keegstra K, Raikhel NV (2001) Characterization of a Family of arabidopsis genes related to xyloglucan fucosyltransferase1. Plant Physiol 127(4):1595
Sayers EW, Bolton EE, Rodney Brister J, Canese K, Chan J, Comeau DC, Connor R et al (2022) Database resources of the national center for biotechnology information. Nucleic Acids Res 50(D1):D20-26
Schneider M, Al-Shareffi E, Haltiwanger RS (2017) Biological functions of fucose in mammals. Glycobiology 27(7):601–618
Schubert M, Bleuler-Martinez S, Butschi A, Wälti MA, Egloff P, Stutz K, Yan S et al (2012) Plasticity of the β-trefoil protein fold in the recognition and control of invertebrate predators and parasites by a fungal defence system. PLoS Pathog 8(5):e1002706
Shimizu K, Chiba S, Saito T, Kumano K, Takahashi T, Hirai H (2001) Manic fringe and lunatic fringe modify different sites of the notch2 extracellular region, resulting in different signaling modulation. J Biol Chem 276(28):25753–25758
Sirén J, Välimäki N, Mäkinen V (2014) HISAT2-fast and sensitive alignment against general human population. IEEE/ACM Trans Comput Biol Bioinforma 11(2):375–388
Sommer R, Makshakova ON, Wohlschlager T, Hutin S, Marsh M, Titz A, Künzler M, Varrot A (2018) Crystal structures of fungal tectonin in complex with o-methylated glycans suggest key role in innate immune defense. Structure 26(3):391-402.e4
Sparrow DB, Chapman G, Wouters MA, Whittock NV, Ellard S, Fatkin D, Turnpenny PD, Kusumi K, Sillence D, Dunwoodie SL (2006) Mutation of the LUNATIC FRINGE gene in humans causes Spondylocostal dysostosis with a severe vertebral phenotype. Am J Hum Genet 78(1):28–37
Spatafora JW, Chang Y, Benny GL, Lazarus K, Smith ME, Berbee ML, Bonito G, Corradi N, Grigoriev I, Gryganskyi A, James TY (2016) A phylum-level phylogenetic classification of zygomycete fungi based ongenome-scale data. Mycologia 108(5):1028-1046
Stahl M, Uemura K, Ge C, Shi S, Tashima Y, Stanley P (2008) Roles of Pofut1 and O-fucose in mammalian notch signaling. J Biol Chem 283(20):13638–13651
Stanton KA, Edger PP, Puzey JR, Kinser T, Cheng P, Vernon DM, Forsthoefel NR, Cooley AM (2017) A whole-transcriptome approach to evaluating reference genes for quantitative gene expression studies: a case study in Mimulus. G3 Genes Genomes Genet 7(4):1085–1095
Staudacher E, Altmann F, Wilson IB, März L (1999) Fucose in N-glycans: from plant to man. Biochem Biophys Acta 1473(1):216–236
Suga H, Chen Z, de Mendoza A, Sebé-Pedrós A, Brown MW, Kramer E, Carr M et al (2013) The Capsaspora genome reveals a complex unicellular prehistory of animals. Nat Commun 4:2325
Tian D, Zhang Li, Zhang S, Kong X, Sheng K, Wang J, Zhang M, Wang Y (2020) Heterologous expression and molecular binding properties of aoflea, a fucose-specific lectin from nematophagous fungus Arthrobotrys oligospora. Int J Biol Macromol 149(April):1213–1221
Tielker D, Hacker S, Loris R, Strathmann M, Wingender J, Wilhelm S, Rosenau F, Jaeger K-E (2005) Pseudomonas aeruginosa lectin LecB is located in the outer membrane and is involved in biofilm formation. Microbiology 151(Pt 5):1313–1323
Turbe-Doan A, Record E, Lombard V, Kumar R, Levasseur A, Henrissat B, Garron ML (2019) Trichoderma reesei dehydrogenase, a pyrroloquinoline quinone-dependent member of auxiliary activity family 12 of the carbohydrate-active Enzymes database: Functional and structural characterization. Appl Environ Microbiol 85(24):e00964-19
UniProt Consortium (2021) UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res 49(D1):D480–D489
Varrot A, Basheer SM, Imberty A (2013) Fungal lectins: structure, function and potential applications. Curr Opin Struct Biol 23(5):678–685
Wang H, Zhang C, Chen H, Yang Q, Zhou X, Zhennan Gu, Zhang H, Chen W, Chen YQ (2016) Biochemical characterization of an isoform of GDP-D-mannose-4,6-dehydratase from Mortierella alpina. Biotech Lett 38(10):1761–1768
Wimmerova M, Mitchell E, Sanchez J-F, Gautier C, Imberty A (2003) Crystal structure of fungal lectin: six-bladed beta-propeller fold and novel fucose recognition mode for Aleuria aurantia lectin. J Biol Chem 278(29):27059–27067
Yan S, Jin C, Wilson IB, Paschinger K (2015) Comparisons of Caenorhabditis fucosyltransferase mutants reveal a multiplicity of isomeric N-glycan structures. J Proteome Res 14(12):5291–5305. https://doi.org/10.1021/acs.jproteome.5b00746
Yan S, Vanbeselaere J, Jin C, Blaukopf M, Wöls F, Wilson IBH, Paschinger K (2018) Core richness of N-glycans of Caenorhabditis elegans: a case study on chemical and enzymatic release. Anal Chem 90(1):935. https://doi.org/10.1021/acs.analchem.7b03898
Yates AD, Allen J, Amode RM, Azov AG, Barba M, Becerra A, Bhai J et al (2022) Ensembl genomes 2022: an expanding genome resource for non-vertebrates. Nucleic Acids Res 50(D1):D996-1003
This work was supported by National Science Centre grant #2021/41/B/NZ2/02426 to AM.
Ethics approval and consent to participate
Adherence to national and international regulations
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. CLANS clustering of GT-A clan peptidyltransferases present in Fungi with animal representatives of transferases acting on glucosylated, fucosylated and galactosylated peptides. After clustering of fungal Fringe homologs with 42 animal representatives of both LFNG and C1GLT transferases, we obtained four clusters with 4946 fungal sequences. However, only four fungal sequences (Rozella allomycis RKP19558.1,RKP21212.1, EPZ36320.1 and Fusarium oxysporum EXK84998.1) group together with the human LFNG (Q8NES3) protein sequence. Two major fungal clusters are equally distant from LFNG and C1GLT clans. Basidiobolus meristosporus ORX91553.1 and Batrachochytrium salamandrivorans KAH6567989.1, KAH6579126.1 and B. dendrobatidis EGF83417.1 sequences group together with the animal C1GLT sequences. We provide all Fringe-like accessions in Supplementary Table S1 despite unresolved specificity questions.
. Phylogenetic trees for fungal fucosyltrasferases.
. Supplementary Table S1.
About this article
Cite this article
Orłowska, M., Barua, D., Piłsyk, S. et al. Fucose as a nutrient ligand for Dikarya and a building block of early diverging lineages. IMA Fungus 14, 17 (2023). https://doi.org/10.1186/s43008-023-00123-8
- Fucose degradation
- Fungal proteins
- Fungal enzymes
- Fucose biosynthesis
- Protein family expansion
- Fucose metabolism