Genomic analysis of ant domatia-associated melanized fungi (Chaetothyriales, Ascomycota)

Several species of melanized (“black yeast-like”) fungi in the order Chaetothyriales live in symbiotic association with ants inhabiting plant cavities (domatia) or with ants that use carton-like material for the construction of nests and tunnels. To investigate the genetic basis and evolution of this lifestyle, genomes of a lineage of four chaetothyrialean strains were sequenced and the genes annotated. While many members of Chaetothyriales have a large ensemble of enzymes enhancing tolerance of extreme or toxic environmental conditions, such as soil polluted with benzene, toluene, ethylbenzene, and xylene (BTEX), members of the domatia symbiont clade experienced significant gene family contractions. This includes enzymes involved in detoxification processes such as cytochrome P450s, transporters, and alcohol dehydrogenases. Overall, the genomes of domatia-associated species are relatively small compared to other Chaetothyriales, with low numbers of protein-coding genes and with a high content of repetitive elements. Biosynthetic clusters involved in the production of secondary metabolites and with potential antimicrobial activities are overrepresented in the genomes of these fungi. We speculate that, despite the reduction of several protein families, members of the domatia-associated clade might tolerate, and perhaps even metabolize, toxic compounds produced from exocrine glands of the ants as defense against microbes. In addition, in this symbiotic association, the plant and the ants could benefit from the production of secondary metabolites by the Chaetothyriales that participate in this tripartite association. We consider a new ecological classification for Chaetothyriales based on genomic features: (i) derived species with high abundance of paralogs colonizing habitats rich in polyaromatic and (ii) potential producers of secondary metabolites with antimicrobial activities, beneficial for symbiotic interactions, occupying specific micro-habitats such as ant domatia.


Introduction
Black yeast-like fungi in the order Chaetothyriales are known for their preference for toxic (Seyedmousavi et al. 2014) or hostile climatic (Selbmann et al. 2015) conditions-lifestyles that may have enhanced their ability to cause infection in human hosts (Moreno et al. 2018d). In addition to less stressful habitats, such as soil and decaying plat material, several members of Chaetothyriales are polyextremophiles, colonizing a wide range of hostile environments. Some of the black yeast-like fungi living in symbiotic association with ants, however, seem to occupy restricted ecological niches and might be specialists ( Fig. 1) (Voglmayr et al. 2011;Vasse et al. 2017).
Ant-associated Chaetothyriales can be classified in three ecological categories based on the nature of the fungal-ant interaction: (a) living in ant-occupied domatia in plants, (b) colonizing ant-made carton structures, and (c) living as parts of fungus-gardens of leaf-cutter ants (Vasse et al. 2017). Each type of ant-plant symbiosis seems to be associated with a particular (series of) chaetothyrialean species. Domatia are mostly inhabited by only a single, or rarely up to three species (Voglmayr et al. 2011;Nepel et al. 2016). As hundreds of ant species are known to live in symbiosis with plants (Chomicki and Renner 2015), and as the ant species investigated thus far often live with black yeast-like fungi (Mayer et al. 2014), many new ant-plant associated Chaetothyriales remain to be d i s co v e r e d . I n ge n e r a l , t he do m at i a -a s s o c i a t e d Chaetothyriales are slow-growing and mesophilic. A recent phylogenetic classification revealed that many of the melanized ant-fungi cluster in a single clade, but some species are scattered across the phylogeny of the Chaetothyriales (Voglmayr et al. 2011;Vasse et al. 2017). Some of the generalist extremotolerant species, such as a fungus close to (99.4% ITS similarity) the hydrocarbon degrader Exophiala oligosperma, have also been found colonizing ant-plant domatia (Vasse et al. 2017), suggesting that hydrocarbon degradation may be an essential factor in symbiosis with ants.
Fungi of all three categories of ant-association described above experience pronounced environmental challenges. One Fig. 1 Chaetothyrialean domatia symbionts of the Barteria fistulosa/ Tetraponera aethiops (a-e) and Cecropia obtusifolia/Azteca constructor (f, g) ant/plant symbioses. a Opened stem of Barteria fistulosa with Tetraponera aethiops; black patch containing chaetothyrialean symbiont fungi. b, c Strain CBS 134920 on 2% MEA showing conidiophores with unicellular pyriform conidia (b) and surface mycelium from the colony center with yellow granular contents (c). d, e Strain CBS 135597 on 2% MEA showing conidiophores with bicellular cyclindrical conidia. f Opened stem of Cecropia obtusifolia with small Azteca constructor workers and ring-shaped fungal patches. g Strain CBS 132003, showing conidiophores and unicellular clavate conidia. Bar: b-e, g 10 μm; f 2 cm of these is mechanical stress due to chewing by the ants, and another is the presence of toxic exhaust to keep the nests free of pathogens (Defossez et al. 2009). On the other hand, plant domatia may confer a certain degree of protection against abiotic stress, preventing fungi from dehydration, UV radiation exposure, and temperature variation. The presence of aromatic acids, mostly having antimicrobial activity, could be a main intrinsic characteristic of this habitat. Ants possess numerous exocrine glands producing a large diversity of compounds (Hölldobler and Wilson 1990), and this characteristic has been proposed as a key selective agent responsible for dominance of Chaetothyriales in domatia (Voglmayr et al. 2011). Interestingly, the major class of chemicals in the metapleural gland secretion of some attine ants are phenylacetic acids (PAA). Some Chaetothyriales are able to metabolize PAA during the breakdown of xenobiotics, via homogentisic acid, possibly as an alternative nutrient source (Cox et al. 1993;Cox et al. 1996;Prenafeta-Boldu et al. 2006;Rustler et al. 2008;Moreno et al. 2018a).
In this context, determination of genome and gene content of Chaetothyriales living in association with ants may help to unravel the functional capacities of these fungi and how they emerged during the evolution of fungal symbiosis with ants. To address this question, we sequenced the whole genome of a selected number of domatia-inhabiting members of the order Chaetothyriales. Their genomic landscape was compared with 24 previously sequenced black yeasts of the chaetothyrialean families Herpotrichiellaceae, Cyphellophoraceae, and Trichomeriaceae. We analyzed key gene families (cytochrome P450, transporters, alcohol dehydrogenase, CAZymes) thought to be essential for the success of these organisms in colonizing extreme and toxic habitats, and which might also play a role in recurrent opportunism in the most derived family, Herpotrichiellaceae.

Strains, growth conditions, and DNA extraction
For extracting the genomic DNA, four strains (CBS 134916, CBS 135597, CBS 132003, CBS 134920; Table 1) were obtained from the reference collection of Centraalbureau voor Schimmelcultures housed at Westerdijk Fungal Biodiversity Institute and cultivated in malt extract broth (MEB) for 7 days, with shaking at 150 rpm at 25°C. The mycelia were separated using 0.2-μm filters and frozen by liquid nitrogen followed by grinding with a mortar and pestle. DNA extraction was p e r f o r m e d f r o m t h e f i n e f r o z e n p o w d e r v i a a cetyltrimethylammonium bromide (CTAB)-based method followed by phenol-chloroform/isoamyl alcohol purification (Moller et al. 1992).

Gene prediction
For gene prediction, the software Augustus (Stanke and Waack 2003) was initially used for generating training gene sets and therefore used as input the models of each organism produced for Genemark-ES v4.30 (Lomsadze et al. 2005). The trained models were used again by Augustus (Stanke and Waack 2003) for the gene prediction in the masked genomic sequences. In addition to protein-coding genes, tRNAs were predicted by tRNAscan-SE v1.23 (Lowe and Eddy 1997). The protein sequences were functionally annotated by InterProScan (database version 5.27-66.0) and BLAST against Swiss-Prot database.

Functional annotation
To predict CYP450s, we identified proteins carrying the PFAM domain PF00067 and compared, by using BLAST, their sequences against the CYP previously described by Teixeira et al. (2017). Putative CYP450 genes were classified into families and subfamilies as recommended by the International P450 Nomenclature Committee (Nelson 2006) using the BLAST identity cutoff of 35% for families and 55% for subfamilies. CYP proteins that could not be assigned to families or subfamilies based on the International P450 Nomenclature Committee were aligned and subjected to phylogenetic analyses, as described elsewhere (Chen et al. 2014a).
Carbohydrate-Active Enzymes (CAZymes) were identified using the dbCAN2 web server (Zhang et al. 2018). For comparative studies, we repeated the same analyses using the 24 previously sequenced black yeasts (Teixeira et al. 2017).
Protein sequences corresponding to the MAT1-1 or MAT1-2 genes, as well as the flanking genes APN2, SLA2, APC5, and COX13, commonly found conserved among Eurotiomycetes, were obtained from Teixeira et al. (2017) and compared by BLAST searches against the newly sequenced genomes.
In order to compare the protein content among the black yeasts, we built clusters of orthologs using the OrthoMCL pipeline v1.4 with a Markov index of 1.5 and maximum evalue of 1 × 10 -5 . The protein sequences of 20 black yeasts included in this analysis were downloaded NCBI: To assess phylogenetic relationships of the newly sequenced species, the amino acid sequences of each single-copy orthologs (SCO) were aligned with MUSCLE (Edgar 2004) and poorly aligned regions were automatically removed using TRIMAL v1. To identify and classify clusters of genes involved in the production of secondary metabolites, the genomes of the domatia-associated species were mined by the antiSMASH web server (Weber et al. 2015), fungal version 4.1.0 (https:// fungismash.secondarymetabolites.org/).
Membrane transport proteins were identified and classified as described elsewhere (Teixeira et al. 2017). Briefly, protein sequences were compared using BLASTP (E-value threshold 1e−05) against the transporter sequences available at Transporter Classification Database (TCDB) (Saier et al. 2014).
To predict changes in size of protein families, functional d o m a i n s w e r e c o u n t e d a n d c o m p i l e d f o r 2 1 Herpotrichiellaceae and in the domatia-associated species. Significant contractions of specific InterPro domain were determined by Fisher's exact test using the Python library SciPy.

Code availability
Custom codes used in this paper as well as raw data produced by the prediction of CAZy and InterProScan are available in http://www.github.com/leandrofmoreno/ant-associated_ black_yeasts.

Sequencing, assembly, and annotation
Genome sequences of four Chaetothyriales colonizing ant domatia were determined by Illumina high-throughput sequencing and de novo assembled. Assembly size varied from 20.6 Mbp, in the fungal species isolated from plant-ant domatium of Barteria fistulosa-Tetraponera aethiops, to 26.3 Mbp in the species associated with Saraca thaipingensis-Cladomyrma petalae (Table 1 and Fig. 2). The GC content ranged from 50.55 to 54%, being consistent with other black yeast species (Teixeira et al. 2017). The content of repetitive elements was considered high and unique among the order Chaetothyriales, varying from 4.09 to 16.32%. The majority of the transposons found in these species belonged to the subclasses of retrotransposons LINE and LTR. We assessed the completeness and evaluated our assemblies quantifying the presence of the Core Eukaryotic Genes (CEGs). This analysis revealed that the genomes sequenced in this study are highly complete including at least 93% of 248 CEGs (Table 1). The protein-coding gene composition was determined by ab initio gene prediction method (without experimental evidence). We observed a phylogenetic correlation between genome size and total gene number across the order Chaetothyriales (Fig. 2). Highest counts of genes were found in the isolate CBS 134916 (6719 gene modes), although this number is lower than any other Chaetothyriales sequenced so far (Teixeira et al. 2017). Even more extreme values were found in the isolates CBS 134920, CBS 132003, and CBS 135597 (Table 1), in which compact proteomes containing less than 6000 proteins were estimated. It is worth noting that the genome assembly of CBS 134916 was more fragmented, which could lead to an increase in the number of gene counts caused by coding regions being split.
We compared 323,280 protein sequences from 21 Herpotrichiellaceae, 2 Cyphellophoraceae and 1 Trichomeriaceae in order to group protein homologs into clusters. We obtained 29,551 orthologous clusters, being 679 clusters composed by single-copy orthologues (one copy of gene from each species). Maximum likelihood tree was constructed using concatenated alignments of all the single-copy orthologous identified (Fig. 2), revealing that the ant domatiaassociated black fungi are placed at an intermediated position to Cyphellophoraceae and Trichomeriaceae.
Mitochondrial genomes of the four sequenced species were assembled into single contigs and are AT-rich with G + C content ranging between 23.8% in CBS 134916 to 26.4% in CBS 134920. The sizes of the mitochondrial genomes vary from 40,056 to 66,938 Mbp in CBS 132033 and CBS 134920, respectively. Genome size variation between the different species is correlated to the variable number of introns. All species analyzed in this study possess introns in the following proteincoding genes: cob, nad5, and cox1. Additionally, we observed that the strains CBS 134920 and CBS 134916 hold introns located in the gene coding for a cytochrome c oxidase (cox13).

Gene family expansions and contractions
By comparing the content of protein functional domains identified in the Chaetothyriales colonizing ant domatia against the family Herpotrichiellaceae, we verified 23 significant (Fisher's exact test, q value < 0.05) gene family contractions (Table 2). Such contractions are unique among the Chaetothyriales and provide specific genomic signatures of black domatia-associated fungi. Several contracted families belong to classes of enzymes thought to confer the metabolic plasticity previously described in Herpotrichiellaceae, including cytochrome P450s, transcription factors, alcohol dehydrogenases, and membrane transporters of the superfamily major facilitator, such as sugar transporter-like (Table 2). Protein domains associated with the metabolism of volatile organic compounds (VOC), such as the S-hydroxymethylglutathione (IPR006913), as well as flavin monooxygenase-like (IPR020946), which constitute a family of xenobioticmetabolizing enzymes, were also found depleted, but not absent, in this group of fungi.

Cytochrome P450s (CYPs)
The plethora of CYPs was assessed and revealed that the domatia-associated species possess the smallest CYP repertory among the Chaetothyriales analyzed so far, ranging from 36 to 46 genes in CBS 134926 and CBS 135597, respectively. The entire repertoire of CYPs of the four sequenced species contains 115 CYPs classified into 24 families and 36 subfamilies (Table S1). In addition, several CYPs are found in the domatia-colonizing species, with homologs in non-related fungi, but not other Chaetothyriales (Table S1). Four CYPs could not be assigned to any already defined CYP family or subfamily following the International P450 Nomenclature Committee. Non-classified CYPs were considered as speciesspecific CYPs (Table S1). Seven CYP families are shared by all four analyzed fungi: CYP617, CYP532, CYP505, CYP545, CYP504 (phenylacetate 2-hydroxylase), CYP628, and CYP53 (benzoate 4-hydroxylase).

Carbohydrate-active enzymes (CAZymes)
The number of CAZyme families is low compared to Chaetothyriales in other families, with the exception of the dermatitidis clade (Fig. 3). Species associated with the ant host Tetraponera aethiops (CBS 134920 and CBS 135597) possess the pectate lyase subfamilies 4 (PL1_4), g2833.t1 and g3705.t1, respectively, responsible for degradation and modification of plant-derived polysaccharides. The presence of pectate lyase in the order Chaetothyriales is exceptional, thus far having been reported in members of the family Cyphellophoraceae and in Capronia coronata only (Teixeira et al. 2017). Comparative analysis with distantly related ascomycetes showed that highly conserved homologs of pectate lyase are shared with some important fungal plant pathogens, such as Diaporthe helianthi and Grosmannia clavigera (Fig. 3). Our analyses confirmed that other polysaccharide lyase families, such as PL3, PL4, PL7, PL9, and PL10, have been lost completely in black yeasts associated to ant domatia. Enzymes linked to degradation of cellulose, such as the glycoside hydrolase subfamilies GH5_5 and GH5_22, were completely absent in the domatia-associated species. Such enzymes were previously found among other black yeast-like fungi, but as yet their activity has not been tested (Table S2) (Teixeira et al. 2017). Similarly, enzymes involved in the breakdown of many polysaccharides are absent, such as the polymerics xylan (GH10, GH11, GH62, GH67, GH115, CE1, and CE15), galactomannan (GH5_7, GH26, GH27, GH36, and GH134), and xyloglucan (GH12, GH29, GH74, and GH95). While an apparent reduction was noted of enzymes related to the degradation of cellulose, the abundance of enzymes in chitin catabolism (chitinases and N- Fig. 2 Genomic landscape of the order Chaetothyriales. Phylogenomic tree was built based on the amino acid sequences of 679 single-copy orthologs shared by all the studied species. Arrows in the node of the phylogenetic tree represent gene family expansion (green) and contraction (red) events that occurred during evolution of black yeasts. The abundance of protein domains corresponding to the major facilitator superfamily (MFS), alcohol dehydrogenase (ADH), cytochrome P450 (P450), and transcription factor (TRP) are represented in each species by bars colored in light purple, dark purple, light orange, and dark orange, respectively. Segments that are not clustered in the MAT locus are indicated by double slash acetylglucosaminidases) was consistent across black yeasts, with CAZy families GH18 (chitinases), GH20 (Nacetylglucosaminidases), GH75 (chitosanases), GH2 (glucosaminidase), and CE4 (chitin deacetylases) (Table S2).

Secondary metabolism
The genomes were screened for the prediction of clusters of genes involved in the production of secondary metabolites. On average, the species possess 19 biosynthetic clusters, with type I PKS (t1pks) being the most prevalent cluster type ranging from 10 copies in CBS 134916 to 14 copies in CBS 132003 (Table S3). Several of these genes were homologous with enzymes for the production of antimicrobial compounds (Table 3). Type III PKS cluster (t3pks), previously reported in Herpotrichiellaceae, Cyphellophoraceae, and Trichomeriaceae, was not observed in the domatia-associated species. The hybrid cluster III PKS/I PKS was found only in CBS 135597 and CBS 134920 that both lived in symbiotic relationship with the ant host Tetraponera aethiops and the plant host Barteria fistulosa (Table S3).

Asexual and sexual reproduction
The sexual identity of the sequenced strains was investigated by homology with the mating-type locus (MAT) structure previously described in Chaetothyriales. Strains proved to be heterothallic, since either MAT1-1 or MAT1-2 genes were found in a single scaffold, in addition to the flanking genes SLA2-APN2-COX13. CBS 134916 and CBS 134920 had MAT1-1 genes, while CBS 135597 and CBS 132003 harbored MAT1-2 (Fig. 2). In CBS 134920 and CBS 132003 the SLA2 gene was located within the MAT locus, while in CBS 134916 and CBS 135597 this gene was not clustered but was found in a separate scaffold. Regarding protein family domains associated with the parasexual cycle, the heterokaryon incompatibility protein (HET-PF06985) ranged from 10 domains in the fungus CBS 134920 to 31 domains in CBS 134916 (Table S4).

Discussion
We fully sequenced the genome of four ant domatiaassociated strains that represent as yet undescribed species, using high-throughput genome sequencing. We observed that the content of repetitive elements is inversely proportional to the genome size of these species, which have the smallest  (Teixeira et al. 2017). Despite their reduced genome size, these organisms possess the highest repetitive element content in the order Chaetothyriales (Table 1). Similarly to what has been published in Exophiala dermatitidis and Rhinocladiella mackenziei, the mitochondrial genomes of black yeasts colonizing ant domatia are AT-rich and their sizes correlate with the variable number of introns found on mitochondrial genes (Moreno et al. 2015;Teixeira et al. 2017;Vicente et al. 2017;Moreno et al. 2018c). We compared the protein domains of the four black yeasts c o l o n i z i n g a n t d o m a t i a a g a i n s t m e m b e r s o f Herpotrichiellaceae, a derived family in Chaetothyriales, to examine differences in functional content, using hypergeometric test. We observed that the domatiaassociated species underwent contraction in gene families previously described to be expanded in other Chaetothyriales (Chen et al. 2014b;Teixeira et al. 2017;Vicente et al. 2017;Moreno et al. 2018a). At a gross level, the reduced genome size of these species might be explained by a large number of gene family contractions (Fig. 2). Depletion of gene families in this group of fungi seems to be an extremely rare evolutionary event, having been described in the fungus Rhinocladiella mackenziei only (Chen et al. 2014b;Teixeira et al. 2017;Vicente et al. 2017;Moreno et al. 2018a). Among the significantly contracted functional domains, CYPs, transporters, and alcohol dehydrogenases are thought to play an important role for the survival of black yeasts under hostile conditions, conferring, e.g., tolerance to toxic xenobiotics. Expansion/ contraction of transcription factors in Chaetothyriales has not yet been discussed in detail, although these are essential regulators of gene expression. Their expansion correlates with the number of protein-coding genes and seems to be linked to gene duplication (Shelest 2017). This suggests that the increased num ber of duplicated genes in derived Chaetothyriales are still functional and under fine-tuned regulation. The expansion of gene families has been associated with the success of free-living black yeasts in colonizing extreme habitats. Although the environment of the domatium occupied by black yeasts can be harsh and toxic because of the antimicrobial chemicals produced by ants (Fernández-Marín et al. 2006Penick et al. 2018), the chemical stress induced is predictable and the number of compounds relatively low. This could explain the lower number of isoforms for specific gene families in Chaetothyriales specialized in colonizing ant-occupied domatia. As the nature of chemical compounds produced by ants is species-specific (Attygalle and Morgan 1984;Hölldobler and Wilson 1990), we may expect that gene family contraction is correlated with ant-host specialization of Chaetothyriales. Further studies determining the exact chemical nature and content of antimicrobial chemicals within domatia will shed light upon the origin of resistance to toxic compounds in black yeasts.
Due to the abundance and high diversity of CYPs in Chaetothyriales (Teixeira et al. 2017;Vicente et al. 2017), this family of enzymes has been suggested to play an essential role in the adaptation to extreme environments. For example, in the hydrocarbon degrading fungus Cladophialora immunda, CYPs are among the most overrepresented protein domains and are upregulated when the fungus was grown in the presence of toluene (Blasi et al. 2017). Indeed, recent studies have concluded that some families of CYPs are associated with the initial oxidation of the methyl group of toluene being one of the key reactions for the degradation of pollutants (Luykx et al. 2003;Blasi et al. 2017). CYPs are essential for the metabolism of polycyclic aromatic hydrocarbons from multiple anthropogenic sources, such as oil spills and incomplete combustion of organic materials (e.g., coal, oil, petrol, and wood) (Abdel-Shafy and Mansour 2016). It is noteworthy that, despite the observed reductions, frequent CYP families, involved in the degradation of xenobiotics, were found in all analyzed genomes, from ant-associated fungi to human opportunists in Herpotrichiellaceae. A remarkable example is CYP504 (coding for phenylacetate 2-hydroxylase), which is highly duplicated across the black yeasts (Teixeira et al. 2017) but was found as a single copy in the domatia-associated species. This enzyme is overexpressed in fungi when they are grown in the presence of phenylacetic acid (PAA), and it is responsible for the conversion of this chemical into 2hydroxyphenylacetate (Harris et al. 2009). Moreover, PAA is the most abundant active compound secreted by the metapleural gland of attine ants to avoid pathogens entering the fungus gardens (Do Nascimento et al. 1996;Fernandez-Marin et al. 2015). The presence of CYP504 in ant-associated species and its conservation during evolution of the black yeasts suggests that the breakdown of chemicals via PAA by phenylacetate 2-hydroxylase may represent one of the ancient mechanisms for degrading xenobiotics in this group of fungi.
We have mined the genomes of the black yeasts living in symbiotic association with ants in order to identify clusters of genes coding for enzymes related to the production of secondary metabolites. Remarkably, these species harbor the highest number of biosynthetic clusters predicted in Chaetothyriales so far (Table S3). While the average of clusters in the families Herpotrichiellaceae and Cyphellophoraceae was predicted to be 11 and 9, respectively (Teixeira et al. 2017), domatiaassociated species possess on average 19 clusters. In particular, the I PKS cluster is well represented in all strains isolated from domatia. The gene composition of the I PKS clusters identified in these species have homology with other biosynthetic pathways already described in several microorganisms, for example, the antibiotic/antiviral viridicatumtoxin/brefeldin biosynthetic gene cluster, identified in CBS 134920 and CBS 135597, and the pestheic acid biosynthetic gene cluster found in CBS 134916 and CBS 132003. We speculate that, in addition to the production of compounds with antimicrobial activity by the ants (Do Nascimento et al. 1996;Voglmayr et al. 2011), the antibiotic-producing black yeasts have an important role in keeping the domatia free of pathogens. Therefore, black yeasts colonizing ant domatia might benefit from the abundance of aromatic hydrocarbons produced by the ants, which could be used as carbon source via PAA degradation, while helping in the protection of the ant nests by producing antibiotics.
The composition of CAZyme families in fungi seems to be responsible, to some extent, for increased metabolic plasticity with respect to carbon utilization (Benoit et al. 2015). For example, the generalist lifestyle of some fungi is linked with the ability of degrading a diversity of polysaccharides, particularly those present in plant material. Therefore, species possessing a large arsenal of these enzymes seems to be more capable of thriving in a wide range of substrates, being not restricted to a specific ecological niche (de Vries et al. 2017). Comparative genomic analyses revealed that the black yeasts and relatives own a complex set of CAZymes, but experienced important gene losses in the families associated with the breakdown of pectin. The absence of such enzymes in the genome of Chaetothyriales suggests that these organisms, including the domatia-colonizing species, might have limited ability to obtain nutrients from plant material. The composition of CAZymes found in the black yeasts associated to ant domatia suggests that these organisms are unlikely to be efficient polysaccharide decomposers; this contrasts with the classic view of mutualistic fungus-ant relationship, where the fungal partner degrades plant material incorporated by the ants (De Fine Licht et al. 2010).
Upregulation of membrane transporters and increased gene dosage of this family of proteins are well-known mechanisms to acquire resistance against a variety of stressors (Sanglard et al. 2009;Morschhauser 2010;Yew et al. 2016). In the genomes of fungi colonizing ant domatia, the number of membrane proteins is reduced compared to other Chaetothyriales (Blasi et al. 2017;Teixeira et al. 2017). This finding is consistent with the contraction of the major facilitator, sugar transporter-like domain IPR005828 represented in Fig. 2. Despite this reduction, members of the investigated lineage have most of the conserved transporters reported previously in the order Chaetothyriales; these are likely involved in niche adaptation to hostile environments, although with limited gene dosage. For example, the toxin efflux pumps of the class trichothecene (TC 2.A.1.3.47) are highly abundant in most of the Chaetothyriales (i.e., 65 gene copies in Cladophialophora immunda) (Blasi et al. 2017;Teixeira et al. 2017), while the domatia-associated lineages harbor 9 to 10 gene copies coding for this transporter. These data suggest that, similar to CYPs and alcohol dehydrogenases, transporters are vastly duplicated in the derived Chaetothyriales which are not associated with ant domatia; this may be related to their lifestyle in environments where nutrient access requires complex metabolic plasticity, such as soil polluted with benzene, toluene, ethylbenzene, and xylene (BTEX) and shells of the babassu coconut palm tree containing a rich diversity of terpenoids and fatty acids (Badali et al. 2011;Nascimento et al. 2017).
By assessing the MAT locus of the sequenced strains, we propose heterothallism as prevalent reproductive mode, because a single mating type idiomorph (either MAT1-1 or MAT1-2) was found per haploid genome. In contrast to the MAT locus in the family Herpotrichiellaceae, in the ant-associated Chaetothyriales, the syntenic flanking genes are organized in way that resembles that of black fungi in the families Cyphellophoraceae and Trichomeriaceae, and of other Eurotiomycetes: APC5-COX13-APN2 (Fig. 2) (Coppin et al. 1997;Fraser et al. 2007;Teixeira et al. 2017;Moreno et al. 2018b). An interesting pattern within the genomes of several previously sequenced black yeasts was the high abundance of hypothetical protein coding genes and HET protein domains (Teixeira et al. 2017), which was not observed in CBS 134916, CBS 134920, CBS 135597, and CBS 132003. The parasexual cycle is an important source of genetic diversity in fungi, in which the HET locus confers non-self-discrimination during vegetative growth, and has been proposed to function as defense systems against competitors and pathogens (Glass and Dementhon 2006). Another incompatibility system characterized in fungi includes the participation of the WD-repeat domain (Espagne et al. 2002;Paoletti and Saupe 2009). Contrary to the abundance of HET domains that vary considerably across the black yeasts, the number of WD-repeat domains seems to be constant in Chaetothyriales.
Overall, our findings suggest that Chaetothyriales living in domatia together with ants possess more compact genomes, and despite their lower diversity of enzymes involved in the metabolism of toxic compounds, this lineage is expected to tolerate and eventually utilize xenobiotics as carbon source.
Funding Acquisition of fungal genomic data was partially supported through funds from the Faculty of Life Sciences, University of Vienna, to VM.
Open Access This article is distributed under the terms of the Creative Comm ons 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.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.