Marie Vasse

Hermann Voglmayr

Veronika Mayer

Cecile Gueidan

Maximilian Nepel

Leandro Moreno

Sybren G De Hoog

Marc-Andre Selosse

Doyle Mckey

Rumsais Blatrix
email: rumsais.blatrix@cefe.cnrs.fr

Cécile Gueidan

Marc-André Selosse

Rumsaïs Blatrix

A phylogenetic perspective on the association between ants (Hymenoptera: Formicidae) and black yeasts (Ascomycota: Chaetothyriales)

Keywords: Symbiosis, mutualism, ant-plant, Chaetothyriales, phylogeny. RIM3substitutions3B3site

published or not. The documents may come

Introduction

Fungi of the order Chaetothyriales (Ascomycetes), also referred to as 'black yeasts and relatives', are mostly known from a range of oligotrophic or extreme environments, such as bare rocks, plant surfaces, indoor surfaces of buildings and substrates contaminated with aromatic hydrocarbons [START_REF] Gueidan | A reappraisal of orders and families within the subclass Chaetothyriomycetidae (Eurotiomycetes, Ascomycota)[END_REF][START_REF] Prenafeta-Boldú | Fungi growing on aromatic hydrocarbons: biotechnology's unexpected encounter with biohazard?[END_REF][START_REF] Badali | Cladophialophora psammophila, a novel species of Chaetothyriales with a potential use in the bioremediation of volatile aromatic hydrocarbons[END_REF][START_REF] Chomnunti | Phylogeny of Chaetothyriaceae in northern Thailand including three new species[END_REF]. Some Chaetothyriales are also opportunistic human pathogens causing various diseases ranging from skin to neurotropic infections [START_REF] De Hoog | Black fungi: clinical and pathogenic approaches[END_REF]. Their extremophilic adaptation-high resistance to a physical or chemical stress-is accompanied by very low competitive ability in mild, buffered environments. Several strains of Chaetothyriales have recently been found in symbiotic association with ants, but these remain poorly investigated.

Most of them appear to be species new to science. As ants produce a large diversity of secondary metabolites, it has been suggested that Chaetothyriales might be predisposed to occupy ant nests because of their particular capacity to metabolise various chemical substances, including aromatic hydrocarbons and other secondary metabolites [START_REF] Nascimento | Arthrocladium, an unexpected human opportunist in Trichomeriaceae (Chaetothyriales)[END_REF][START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF]. As ants are expected to shape their microbial environment by the production of multiple antimicrobial substances, investigating these associations will contribute to elucidating the processes underlying the evolution of insect-microbe symbioses, for which convergence is a pervasive theme in insect diversification [START_REF] Mueller | The evolution of agriculture in insects[END_REF][START_REF] Gibson | Extraordinarily widespread and fantastically complex: comparative biology of endosymbiotic bacterial and fungal mutualists of insects[END_REF]. Ant-associated Chaetothyriales can be classified into three main functional groups: [START_REF] Gueidan | A reappraisal of orders and families within the subclass Chaetothyriomycetidae (Eurotiomycetes, Ascomycota)[END_REF] strains associated with fungus-growing attine ants [START_REF] Little | Symbiotic complexity: discovery of a fifth symbiont in the attine ant-microbe symbiosis[END_REF][START_REF] Attili-Angelis | Novel Phialophora species from leaf-cutting ants (tribe Attini)[END_REF][START_REF] Duarte | Leaf-cutting ants: an unexpected microenvironment holding human opportunistic black fungi[END_REF]; [START_REF] Prenafeta-Boldú | Fungi growing on aromatic hydrocarbons: biotechnology's unexpected encounter with biohazard?[END_REF] strains involved in 'carton' structures built by ants [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF][START_REF] Ruiz-González | Specific, non-nutritional association between an ascomycete fungus and Allomerus plant-ants[END_REF][START_REF] Schlick-Steiner | Specificity and transmission mosaic of ant nest-wall fungi[END_REF][START_REF] Nepel | High diversity and low specificity of Chaetothyrialean fungi in carton galleries in a neotropical ant-plant association[END_REF]; and (3) strains associated with ant-plant symbioses [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF][START_REF] Mayer | Current issues in the evolutionary ecology of ant-plant symbioses[END_REF].

The first group, Chaetothyriales associated with attine ants, are not the focus of our study. However, because our analysis allows a comparison of their phylogenetic placement, they are briefly presented here. The nature of the association between Chaetothyriales and attine ants is mostly unknown. It was first proposed that Chaetothyriales were detrimental to the symbiosis between attine ants and their cultivated fungi because they were shown to lower the growth of a specific bacterium (Pseudonocardia) known to protect the cultivated fungus through the production of antibiotics [START_REF] Little | Black yeast symbionts compromise the efficiency of antibiotic defenses in fungus-growing ants[END_REF]. However, the proposition that black yeasts inhibit Pseudonocardia was based on in vitro experiments confronting these two organisms on artificial medium. Since then, the hypothesis of coevolution between attine ants and Pseudonocardia has been challenged and the system is now viewed as a complex community of bacteria recruited by the ants from the environment [START_REF] Mueller | Symbiont recruitment versus ant-symbiont co-evolution in the attine ant-microbe symbiosis[END_REF]. The role of Chaetothyriales in this context has not yet been re-investigated. Nevertheless, evidence of the presence of Chaetothyriales on the cuticle of attine ants is accumulating [START_REF] Attili-Angelis | Novel Phialophora species from leaf-cutting ants (tribe Attini)[END_REF][START_REF] Duarte | Leaf-cutting ants: an unexpected microenvironment holding human opportunistic black fungi[END_REF].

Several ant species make a cardboard-like material ('carton') to build their entire nest, partition the space available in natural cavities, cover and protect their hemipteran trophobionts, or construct runways and galleries for protection or for capturing large prey.

Fungal mycelium growing within ant carton-and cultured by the ants themselves-has long been known to contribute to structural support, stabilizing the construction material used as carton components (i.e. plant debris or fibres) [START_REF] Lagerheim | Über Lasius fuliginosus (Latr.) und seine Pilzzucht[END_REF][START_REF] Maschwitz | Der Kartonnestbau bei Lasius fuliginosus Latr. (Hym. Formicidae)[END_REF][START_REF] Weissflog | Freinestbau von Ameisen (Hymenoptera, Formicidae) in der Kronenregion feuchttropischer Wälder Südostasiens[END_REF][START_REF] Dejean | Arboreal ants build traps to capture prey[END_REF][START_REF] Mayer | Mycelial carton galleries of Azteca brevis (Formicidae) as a multi-species network[END_REF]. However, the taxonomic affiliation of these fungi was unknown until recently and the frequency and extent of associations with ants were underestimated. Fungi of various orders (e.g. Botryosphaeriales, Eurotiales, Capnodiales, Chaetothyriales, Hypocreales, Pleosporales) can be isolated from carton. Only strains of the orders Chaetothyriales and Capnodiales are isolated from pieces of hyphae and are thus considered to have a structural function in carton construction [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF][START_REF] Ruiz-González | Specific, non-nutritional association between an ascomycete fungus and Allomerus plant-ants[END_REF][START_REF] Schlick-Steiner | Specificity and transmission mosaic of ant nest-wall fungi[END_REF][START_REF] Nepel | High diversity and low specificity of Chaetothyrialean fungi in carton galleries in a neotropical ant-plant association[END_REF]. The other strains have been isolated only from spores and are most likely contaminants. Numerous strains of Chaetothyriales are usually mixed within a single piece of carton, and show low specificity in regard to the identity of the carton-making ant [START_REF] Nepel | High diversity and low specificity of Chaetothyrialean fungi in carton galleries in a neotropical ant-plant association[END_REF][START_REF] Mayer | Mycelial carton galleries of Azteca brevis (Formicidae) as a multi-species network[END_REF]. Carton-making is widespread among tropical ants, but also occurs in temperate species [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF][START_REF] Schlick-Steiner | Specificity and transmission mosaic of ant nest-wall fungi[END_REF][START_REF] Lagerheim | Über Lasius fuliginosus (Latr.) und seine Pilzzucht[END_REF]. Some grounddwelling ants construct carton, but carton-making is most widespread in arboreal ants, with some ant genera being particularly renowned for building carton structures (e.g.

Crematogaster, Azteca, Camponotus, Polyrhachis). The capacity to build carton is scattered throughout the ant phylogeny. It has either evolved and/or been lost several times independently, or is facultative or sporadic for all ants.

Chaetothyriales are astonishingly abundant in tropical ant-plant symbioses. In these symbioses, the plants host ants in specialized hollow structures called domatia. Such symbioses are restricted to the tropics, but widespread throughout tropical regions. The ant colony usually lives within the host plant throughout its whole lifespan. Domatia have various anatomical origins that vary across taxa; these include twigs, leaves and stipules [START_REF] Davidson | The evolutionary ecology of symbiotic ant-plant relationships[END_REF]. In addition to providing ants with nesting sites, ant-plants also often contribute to ant nutrition, either directly through the production of extrafloral nectar and food bodies, or indirectly through sap-sucking Hemiptera that ants rear in domatia for honeydew and 'meat' [START_REF] Heil | Protective ant-plant interactions as model systems in ecological and evolutionary research[END_REF]. Antplants benefit from the presence of ants in three ways: (1) protection against herbivores (vertebrates and/or invertebrates) and pathogens, (2) reduced competition through pruning of the surrounding vegetation, and (3) uptake of nitrogen and possibly other nutrients from ant activity [START_REF] Mayer | Current issues in the evolutionary ecology of ant-plant symbioses[END_REF][START_REF] Davidson | The evolutionary ecology of symbiotic ant-plant relationships[END_REF][START_REF] Rosumek | Ants on plants: a meta-analysis of the role of ants as plant biotic defenses[END_REF]. While the occurrence of fungi within ant-occupied domatia has long been recorded, the nature of these fungi and of their relationship with ant-plant symbioses has only recently been investigated [START_REF] Mayer | Current issues in the evolutionary ecology of ant-plant symbioses[END_REF][START_REF] Nepel | Ant-cultivated Chaetothyriales in hollow stems of myrmecophytic Cecropia sp. trees -diversity and patterns[END_REF]. The inner wall of each domatium has a dark, thin patch where the fungal symbionts are located. Fungal patches within domatia are mainly constituted of a few strains of Chaetothyriales that initially appeared to be more or less specific to the ant species [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF][START_REF] Nepel | Ant-cultivated Chaetothyriales in hollow stems of myrmecophytic Cecropia sp. trees -diversity and patterns[END_REF][START_REF] Blatrix | Repeated evolution of fungal cultivar specificity in independently evolved ant-plant-fungus symbioses[END_REF][START_REF] Kokolo | Pattern of specificity in the tripartite symbiosis between Barteria plants, ants and Chaetothyriales fungi[END_REF]. Fungal patches also contain nematodes, bacteria and spores of opportunistic fungi [START_REF] Blatrix | Ant-plant mutualisms should be viewed as symbiotic communities[END_REF][START_REF] Maschwitz | Nematode associates and bacteria in ant-tree symbioses[END_REF][START_REF] Esquivel | A new species of nematode, Sclerorhabditis neotropicalis sp. n. (Rhabditida), associated with Azteca ants in Cecropia obtusifolia[END_REF]. Detailed functional investigation of a few systems revealed that domatia patches are manured and used as a food source by the ants [START_REF] Defossez | Plant-ants feed their host plant, but above all a fungal symbiont to recycle nitrogen[END_REF][START_REF] Blatrix | Plant-ants use symbiotic fungi as a food source: new insight into the nutritional ecology of ant-plant interactions[END_REF] in a way reminiscent of fungal culture by attine ants. Ants and plants living in symbiosis are scattered throughout their respective phylogenies and have evolved many times independently, at least 40 and 158 times in their respective lineages [START_REF] Davidson | The evolutionary ecology of symbiotic ant-plant relationships[END_REF][START_REF] Chomicki | Phylogenetics and molecular clocks reveal the repeated evolution of ant-plants after the late Miocene in Africa and the early Miocene in Australasia and the Neotropics[END_REF]. Ant-plant symbioses evolved recently (probably not before the Miocene) and are rarely associated with substantial speciation of plant or ant lineages [START_REF] Chomicki | Phylogenetics and molecular clocks reveal the repeated evolution of ant-plants after the late Miocene in Africa and the early Miocene in Australasia and the Neotropics[END_REF][START_REF] Chomicki | Macroevolutionary assembly of ant/plant symbioses: Pseudomyrmex ants and their ant-housing plants in the Neotropics[END_REF].

A previous phylogenetic study indicated that ant-associated Chaetothyriales form a few monophyletic clades, and that each domatia strain was associated with only one ant species, leading the authors to postulate that ant-associated Chaetothyriales from domatia were ant-species specific and specialized to this lifestyle [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF]. The rapid rate at which antassociated species of Chaetothyriales are being discovered suggests that ant-Chaetothyriales associations are extremely frequent. A reappraisal of the phylogenetic distribution of antassociated strains is needed because apparent ant-specificity could result from the undersampling of this poorly known order of fungi. Here, we reconstruct a new multi-gene phylogeny of the Chaetothyriales based on one mitochondrial and three nuclear genes. In addition, we reconstruct phylogenies of subsets of strains using three additional, more variable nuclear markers. Our phylogenetic approach includes ant-associated strains from antplant associations that had not been analysed in previous studies (about one-third of the strains are new). The present study also includes a larger number of Chaetothyriales strains not associated with ants, improving the ability of our phylogeny to clarify the phylogenetic position of ant-associated strains.

Methods

Sampling

Our data set included a total of 242 strains of Chaetothyriales (Table S1): 46 strains from ant-made carton (among which seven are new isolates), 55 strains from domatia patches (among which 25 are new isolates), and 141 strains of various origins intended to cover the broadest possible range of functional and taxonomic diversity in the order Chaetothyriales.

Our choice of strains not associated with ants was constrained by the availability of already published sequences. We cannot exclude the possibility that additional strains would change the topology of the phylogenetic trees of this study. As outgroups, we used five taxa in the order Verrucariales according to the phylogeny of Gueidan et al. [START_REF] Gueidan | A rock-inhabiting ancestor for mutualistic and pathogen-rich fungal lineages[END_REF].

Seven gene regions were used as phylogenetic markers: the small and the large subunits of the nuclear ribosomal RNA gene (nSSU and nLSU, respectively), the small subunit of the mitochondrial ribosomal RNA gene (mtSSU), the fragments AC and DF of the largest subunit of the RNA polymerase II (RPB1), the internal transcribed spacer region (ITS) composed of ITS1, 5.8S and ITS2, the partial beta tubulin gene (Bt2) and the translation elongation factor 1 alpha (EF1). Sequences were produced de novo using the protocols described below or retrieved from GenBank. Many strains with missing information were included in the phylogenetic analyses, in particular strains not associated with ants. We decided to include these strains because of their potential importance for our study according to their taxonomic position or ecology. For more details on missing data, see Table S1 and alignment matrices at TreeBase. Collection of samples of ant-made carton and domatia fungal patches, isolation of fungal strains and DNA extraction from pure cultures followed the protocols described in Voglmayr et al. [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF]. Strains and GenBank accession numbers used in this study are listed in Table S1.

DNA amplification and sequencing

Amplification by Polymerase Chain Reaction (PCR) was performed using either the Qiagen multiplex kit (Qiagen, Venlo, Netherlands) or the Sigma REDExtract-N-Amp PCR Ready Mix (Sigma-Aldrich). Conditions for PCR first followed manufacturer's instructions, and were then adjusted if amplification failed. Amplification used 25 µl of solution containing 12.5 µl of PCR mix (from either the Qiagen or Sigma kits), 5 µl of Q solution or water (Qiagen and Sigma respectively), 2.5 µl of each primer at 5 µM, and 2.5 µl of DNA extract.

Primers used for amplification and sequencing are in Table S2. The PCR program was: 15 min at 95 °C, 35 cycles of 1 min at 94 °C, 1 min at 45-55 °C, 1.5-3 min at 72 °C (depending on expected sequence length), and a final elongation of 20 min at 60 °C. Sanger dideoxy sequencing of the PCR products was performed at the Genoscope (Evry, France).

Sequences were edited using CodonCode Aligner (CodonCode Corporation, Dedham, MA, USA), and contigs were built from forward and reverse sequences generated for each gene.

Conflicting base calls were coded using the ambiguity code.

Phylogenetic analyses

The sequences were aligned manually using Mesquite v.3.04 [START_REF] Maddison | Mesquite: a modular system for evolutionary analysis[END_REF]. We delimited and excluded from the alignments introns and all ambiguously aligned regions following the protocol of Lutzoni et al. [START_REF] Lutzoni | Integrating ambiguously aligned regions of DNA sequences in phylogenetic analyses without violating positional homology[END_REF].

We first used the four least variable markers (nSSU, nLSU, mtSSU, RPB1 region A-C and region D-F) to build the order-level phylogeny with all strains. To test for congruence between these markers, we ran maximum likelihood analyses on each marker matrix with RAxML HPC2 v.8.2.4 [START_REF] Stamatakis | RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models[END_REF] on the Cipres Web Portal [START_REF] Miller | Creating the CIPRES Science Gateway for inference of large phylogenetic trees[END_REF] using a GTR + I + G model and 1000 bootstrap pseudo-replicates. RPB1 regions were partitioned according to the three positions of the codons. Bootstrap consensus trees of each marker were compared by eye for sets of taxa supported by bootstrap values equal to or above 70 %. We excluded the mtSSU sequences of strains CTeY1 and CTeY9 due to incongruent placement. We then combined the four markers into a concatenated matrix and ran a maximum likelihood (ML) analysis to investigate the phylogenetic relationships among the taxa. This ML analysis was conducted on the Cipres Web Portal with RAxML HPC2 v.8.2.4 using a GTR + I + G model on 6 partitions (nSSU, nLSU, mtSSU, first, second and third codon positions of RPB1) as estimated by the Akaike Information Criterion in jModelTest2 v.2.1.9 [START_REF] Darriba | jModelTest 2: more models, new heuristics and parallel computing[END_REF][START_REF] Guindon | A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood[END_REF], and for 1000 bootstrap pseudo-replicates. In addition, we ran a Bayesian analysis on the concatenated matrix with MrBayes v.3.2.6 [START_REF] Ronquist | MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space[END_REF] using the same model and partitions as in the ML analysis.

Two analyses of four chains were run for 10 000 000 generations, sampling trees every 500 generations with 25 % burn-in samples discarded for each run.

Then, to better resolve the phylogenetic relationships among closely related strains, we built phylogenies of subsets of taxa based on three markers, ITS, Bt2 and EF1, containing parts that are more variable than the markers used for the order-level phylogeny. In each subset we included additional strains that were not in the order-level phylogeny, and for which we retrieved sequences from GenBank. We did not remove introns and indels in the alignments of these markers. For each subset, we ran ML analysis on each matrix as described above. After checking the congruence between the markers, we concatenated the matrices and ran maximum likelihood analysis. For consistency we used the same parameters and models as for the main tree. The ITS region is considered to be the most relevant DNA barcode for fungi [START_REF] Schoch | Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi[END_REF]. In order to detect Chaetothyriales species that could potentially occur in association with several ant species, and/or on several continents, we delineated Molecular Operational Taxonomic Units (MOTUs) in ant-associated strains using a cut-off similarity of 99 % between ITS sequences. As our concern was to group strains that were highly likely to belong to a single species, we chose a cut-off value in the upper range observed for Eurotiomycetes [START_REF] Schoch | Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi[END_REF]. A lower value increases the risk of false conclusion regarding the range of potential hosts and the extent of geographic distribution of a particular species. The distribution curve of pairwise ITS similarity for the 38 ant-associated strains in the family Cyphellophoraceae (which proved to form a monophyletic clade, see below) shows a drop around the 99 % value, but not around the 97 % value, indicating that the traditional cut-off value of 97 % is less consistent with discontinuity in sequence similarity than the 99 % value used in our study (Fig. S1). The drop of the distribution curve is interpreted as the transition between the average intraspecific and the average interspecific sequence similarity.

Results

For this study, we generated 93 mtSSU, 27 RPB1, 88 nSSU, 81 nLSU, 79 ITS, 61 EF1 and 58

Bt2 new sequences. GenBank accessions of new sequences are given in the section on data accessibility and in Table S1. The dataset of the order-level phylogeny based on four markers included 6292 characters (1290 for nLSU, 1580 for nSSU, 596 for mtSSU and 2826 for RBP1). Bayesian and maximum likelihood analyses yielded similar trees: we present only the tree from the Bayesian analysis (Fig. 1). The accession URL of the TreeBASE project with alignments and phylogenetic trees is indicated in the section on data accessibility. The families Trichomeriaceae, Herpotrichiellaceae and Cyphellophoraceae, in which most antassociated strains are located, formed well-supported clades.

Ant-associated strains were scattered throughout the phylogeny and branched at various depths in the phylogeny (Fig. 1). Most strains clustered according to whether they originated from carton or from domatia, although a few clades in Trichomeriaceae included both functional types. Ant-associated strains occurred in four out of the five formally described families of Chaetothyriales, although most of them belonged to the Trichomeriaceae, followed by the Cyphellophoraceae, and to a lesser extent by the Herpotrichiellaceae. Various types of association with ants were found in each of these three families. Strains from both domatia patches and ant-made carton occurred in Trichomeriaceae and Herpotrichiellaceae. The Cyphellophoraceae contained most of the strains from domatia, but also strains isolated from the cuticle of attine ant workers and gynes in previous studies [START_REF] Little | Symbiotic complexity: discovery of a fifth symbiont in the attine ant-microbe symbiosis[END_REF][START_REF] Attili-Angelis | Novel Phialophora species from leaf-cutting ants (tribe Attini)[END_REF]. Although there were several independent origins of each functional type, most of the ant-associated strains appeared to cluster in a few clades. For instance, the Trichomeriaceae contained a monophyletic clade of carton strains (Fig. 1) and a large clade with strong support in the Cyphellophoraceae was composed only of domatia strains (Fig. 2a).

Our order-level phylogeny clearly showed the absence of phylogenetic clustering of the ant-associated strains according to their geographic origin (Fig. 1). As illustrated in several clades, strains from two or three continents were sometimes closely related. Similarly, there was no clustering of strains according to the taxonomic relationships between host ants or between host plants. A highly illustrative example is a monophyletic clade in the Trichomeriaceae composed of closely related strains from Central American Azteca (Dolichoderinae) plant-ants on Cecropia (Urticaceae) and Tetrathylacium (Salicaceae) (domatia and carton fungi), south-east Asian Cladomyrma (Formicinae) plant-ants on Saraca (Fabaceae) (domatia fungi) and African Crematogaster (Myrmicinae) (carton fungi) (Fig. 1).

Phylogenies based on sub-sets of strains (Fig. 2) and computation of ITS sequence similarity revealed seven cases -two of which had already been noticed by Nepel et al. [START_REF] Nepel | High diversity and low specificity of Chaetothyrialean fungi in carton galleries in a neotropical ant-plant association[END_REF] -of a single MOTU (> 99 % ITS sequence similarity) comprised of strains isolated from several ant species and/or on several continents. One MOTU (strain T394 from Nepel et al.

[15] and MCRE12) was found in carton nests of Crematogaster sp. in Thailand and in carton galleries of Azteca brevis in Costa Rica (for illustration see Nepel et al. [START_REF] Nepel | High diversity and low specificity of Chaetothyrialean fungi in carton galleries in a neotropical ant-plant association[END_REF]). Two MOTUs from Azteca nests in Costa Rica, Cyph1 (Fig. 2a) and Tric1 (Fig. 2b

Discussion

Our multilocus phylogeny of the Chaetothyriales based on a sample of taxa representative of the entire diversity of the order shows a pattern of phylogenetic scatter of ant-associated strains. These strains do not cluster according to their geographic origin nor to the identity of their species of host ant or host plant. Ant-associated fungal taxa are not necessarily hostspecific, but can occur in niches other than ant nests, and can be distributed worldwide. Many aspects of the evolutionary history of ant-associated Chaetothyriales are still unclear. More data on the nature of the interaction between ants and particular Chaetothyriales taxa (for example, distinguishing mutualistic interactions from others), and more comprehensive sampling (for example, to infer ancestral habitats of ant species-specific taxa) will be required to resolve open questions.

Studies on the pattern of association between ants and Chaetothyriales at both local and larger spatial scales can contribute greatly to understanding the ecology of ant-associated strains. In Cameroon, we investigated Chaetothyriales associated with three ant-plant symbioses co-occurring at the same site: Petalomyrmex phylax on Leonardoxa africana subsp. africana, Tetraponera aethiops on Barteria fistulosa and Crematogaster margaritae on Keetia hispida. We found that each ant-plant symbiosis had its own set of domatia fungal strains [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF][START_REF] Blatrix | Repeated evolution of fungal cultivar specificity in independently evolved ant-plant-fungus symbioses[END_REF]. In Central America, four Azteca species living with three different Cecropia species were investigated along a transect of a few kilometers length. Although some strains were shared across Azteca species, some others were ant-species specific [START_REF] Nepel | Ant-cultivated Chaetothyriales in hollow stems of myrmecophytic Cecropia sp. trees -diversity and patterns[END_REF]. In contrast, the present study shows that several MOTUs are composed of strains isolated from various ant species in various continents, suggesting both low ant-specificity and long-distance dispersal of ant-associated Chaetothyriales. Moreover, strains found in domatia of Tachigali (in French Guiana) and Triplaris (in Costa Rica) occupied by different Pseudomyrmex ant species both correspond to Exophiala oligosperma (99.4 % ITS sequence similarity), a fungus isolated from various low-nutrient environmental substrates, reported to be an opportunistic human pathogen and known to metabolize toluene [START_REF] De Hoog | Species diversity and polymorphism in the Exophiala spinifera clade containing opportunistic black yeast-like fungi[END_REF]47]. In addition, a recently described species, Arthrocladium tardum (Trichomeriaceae), has been isolated both from a domatium of Leonardoxa africana occupied by the ant Petalomyrmex phylax in Cameroon and from decaying coconut shells in Brazil [START_REF] Nascimento | Arthrocladium, an unexpected human opportunist in Trichomeriaceae (Chaetothyriales)[END_REF]. Although some Chaetothyriales strains are clearly mutualistic symbionts of ants [START_REF] Blatrix | Repeated evolution of fungal cultivar specificity in independently evolved ant-plant-fungus symbioses[END_REF][START_REF] Defossez | Plant-ants feed their host plant, but above all a fungal symbiont to recycle nitrogen[END_REF][START_REF] Blatrix | Plant-ants use symbiotic fungi as a food source: new insight into the nutritional ecology of ant-plant interactions[END_REF], our study indicates that ant nests can be colonized by opportunistic species and that Chaetothyriales taxa might display exceptional dispersal capacities and surprising ubiquity. In the future, it will be most useful to investigate the degree of specialization of the interactions between particular ant and fungal partners in order to interpret the phylogenetic pattern in evolutionary terms.

The order-level phylogeny shows that the strains we isolated from ant-made carton or from domatia occupied by plant-associated ants are scattered throughout the phylogeny of the order Chaetothyriales. Among the ant-associated strains added by this study to the phylogeny published by Voglmayr et al. [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF], two (MACrh12a and Pet5a) form long and isolated branches occurring basally in their respective families (Trichomeriaceae and Epibryaceae).

Epiphyllous strains recently described in the genus Trichomerium (Trichomeriaceae) [START_REF] Chomnunti | Phylogeny of Chaetothyriaceae in northern Thailand including three new species[END_REF] are placed within a clade that previously was known to contain only strains from ant-made carton [START_REF] Voglmayr | The diversity of ant-associated black yeasts: Insights into a newly discovered world of symbiotic interactions[END_REF]. In addition to a better knowledge of the functional ecology of ant-Chaetothyriales interactions, we need a more comprehensive sampling of strains in the phylogeny to properly infer the evolutionary history of cases of mutualistic symbioses between ants and Chaetothyriales. Our results suggest that focussing on communities of epiphylls should be the next step to elucidate the functional ecology of black yeasts and their roles in the biology of ants.

Conclusion

The molecular characterisation of fungi associated with tropical arboreal ants opens a window on an unsuspected aspect of the evolution and ecology of the Chaetothyriales. In addition, it reveals complex patterns of association with ants. The functional and evolutionary ecology underlying these patterns remains to be understood. Accomplishing this will require screening the environments surrounding nests of tropical arboreal ants-especially the plant hosts of these ants-for Chaetothyriales, and investigating the mode of transmission of the fungal symbiont. Preliminary results suggest that founding queens may carry fungi from their mother colony in their infrabuccal pockets [48]. Tropical arboreal ants aside, Chaetothyriales have also been found on the cuticles of attine ants and in carton structures of ground-dwelling ants of temperate latitudes. Are these fungi universal symbionts of ants? We may have just barely seen the tip of the iceberg. Each label of an ant-associated strain corresponding to an undescribed species from carton and domatia is composed of the name of the strain, the name of the ant and the name of the host plant (if any). Names of strains from Africa, America and South-East Asia are represented in green, orange and magenta, respectively. Strains from carton, domatia and attine ants have blue, red and yellow branches, respectively. an ant-associated strain corresponding to an undescribed species from carton or domatia is composed of the name of the strain, the name of the ant and the name of the host plant (if any). Names of strains from Africa, America and South-East Asia are represented in green, orange and magenta, respectively. Strains from carton, domatia and attine ants have blue, red and yellow branches, respectively.

Supplementary material

  ), were found respectively in domatia of Keetia hispida and a closely related Rubiaceae occupied by Crematogaster margaritae in Cameroon, and in carton nests of Monomorium sp. in Malaysia. Three MOTUs (Tric2, Herp1 and Herp2, Fig. 2b, c) were found in carton nests of Crematogaster sp. and/or Pheidole sp. in Cameroon and in carton galleries of Azteca brevis in Costa Rica. Finally, one MOTU (Herp3, Fig. 2d), which corresponds to the species Exophiala oligosperma, was found in domatia of Tachigali sp. (Fabaceae) occupied by Pseudomyrmex malignus or P. penetrator in French Guiana, and in domatia of Triplaris sp. (Polygonaceae) occupied by Pseudomyrmex sp. in Costa Rica.

Figure 1 :

 1 Figure 1: Phylogenetic reconstruction of the Chaetothyriales based on a Bayesian approach and four markers (nLSU, nSSU, mtSSU, RPB1), showing the position of strains isolated from ant nests (ant-made carton, ant-occupied domatia of plants symbiotically associated with ants, or attine ants). Strong support values (100 % bootstrap [BS] and posterior probabilities [PP]) are indicated by black dots. Other values are indicated as BS/PP. Dashes replace non-significant values (< 70 % for BS and < 95 % for PP).

Figure 2 :

 2 Figure 2: Phylogenetic trees of subsets of Chaetothyriales strains generated by maximum likelihood analysis based on three markers (ITS, Bt2, EF1). (a) Subset for the family Cyphellophoraceae. (b) Subset for the family Trichomeriaceae. (c) and (d) Subsets for the family Herpotrichiellaceae. Clades corresponding to MOTUs (> 99 % ITS similarity) are highlighted in grey. Those MOTUs that contain strains from various ant genera or continents are named. Bootstrap values above 70% are added to supported branches. Strong support values (100 % BS) are indicated by black dots. Each label of

Figure S1 :

 S1 Figure S1: Distribution of pairwise similarity at ITS for 38 strains of ant-associated Cyphellophoraceae, showing a drop of values at 0.99 but not at 0.97.

  RB, VM, HV and MN carried out field collection and isolation of fungal strains; MV, RB, HV, MN, VM and LM carried out the molecular lab work; MV and CG carried out sequence alignments; MV, RB and CG participated in data analysis; RB, VM, MAS, DM, CG and SDH conceived and designed the study; MV and RB drafted the manuscript. All authors improved the manuscript and gave final approval for publication. Veiga MC, Kennes C. 2005 Biodegradation of toluene by the new fungal isolates Paecilomyces variotii and Exophiala oligosperma. J. Ind. Microbiol. Biotechnol.

	47. Estévez E, 32, 33-37. (doi:10.1007/s10295-004-0203-0)
	48.	Baker CC. 2015 Complexity in mutualisms: indirect interactions with multiple parties.
	PhD	thesis.	Harvard	University,	Cambridge,	MA,	USA.
	https://dash.harvard.edu/handle/1/23845506
	Figure captions

Table S1 :

 S1 List of the fungal strains used in this study, with Genbank accession numbers for DNA sequences.

Table S2 :

 S2 List of the primer pairs used for DNA amplification and sequencing.Figure S1. Distribution of pairwise similarity at ITS for 38 strains of ant-associated Cyphellophraceae, showing a drop of values at 0.99 but not at 0.97.

		:aS3Cyphellophoraceae					:bS3Trichomeriaceae
		OM
										U+
		gP	g+								gN
										gO
		25									gU gU		TricZ
		OM	OM							gg
			gN
			gP
			gg
		20
		Og
									RIM3substitutions3B3site	gN		TricM
				UU	Og
	Frequency	15		OU				CyphM	:cS3Herpotrichiellaceae gg gZ
				gg
									OZ
		10		gL gO					Og
									gg
					FG
										OU			HerpZ
										gg
													HerpM
		5
									:dS3Herpotrichiellaceae
		Colour3code3for3text3and3branchesk					gg
		Strains3from3Africa
		Strains3from3America						gg
		Strains3from3South,East3Asia						UU
		33333333Strains3from3attine3ants 33333333Strains3from3ant,made3carton 33333333Strains3from3domatia 0					gU	gR gR			HerpN
			0.85	0.86	0.87	0.88	0.89	0.90	0.91	0.92	0.93	0.94	0.95	0.96	0.97	0.98	0.99	1.00
											ITS similarity

Funding

This project was supported by the network "Bibliothèque du Vivant" funded by the CNRS, the Muséum National d'Histoire Naturelle, the INRA and the CEA (Centre National de Séquençage).

Data accessibility

DNA sequences: Genbank accessions KX822178-KX822578. Phylogenetic data, including alignments: TreeBASE accession URL http://purl.org/phylo/treebase/phylows/study/TB2:S19475.

Competing interests

We have no competing interests.

Authors' contributions