Jessica Lee Abbate

Pierre Gladieux

Michael E Hood

Damien M De Vienne

Janis Antonovics

Alodie Snirc

Tatiana Giraud
email: tatiana.giraud@u-psud.fr

Damien Antonovics

Co-occurrence among three divergent

Keywords: Adaptation, Fungi, Host Parasite Interactions, Hybridization, Phylogeography, Population Genetics -Empirical population structure, endemicity, speciation, hybrid zones, fungi, Microbotryum violaceum, Silene maritima, biogeography, altitude

come   L'archive ouverte pluridisciplinaire

Introduction

The competitive exclusion principle postulates that two species occupying the same niche in sympatry cannot coexist indefinitely [START_REF] Gause | The Struggle For Existence[END_REF][START_REF] Hardin | The competitive exclusion principle[END_REF]. In pathogens, this means that species coexisting in sympatry in association with a common host species are likely to differ in other components of their ecological niches. The coexistence of multiple pathogens on the same crop species is in fact often associated with separation by time, resource use, abiotic conditions or geographical area [START_REF] Fitt | Coexistence of related pathogen species on arable crops in space and time[END_REF][START_REF] Giraud | Introduction: microbial local adaptation: insights from natural populations, genomics and experimental evolution[END_REF].

Closely related pathogens can also specialize on either the early or the late season part of a given plant life cycle [START_REF] Hamelin | Temporal niche differentiation of parasites sharing the same plant host: oak powdery mildew as a case study[END_REF][START_REF] Mailleret | From elaborate to compact seasonal plant epidemic models and back: is competitive exclusion in the details?[END_REF]. For example, two eyespot wheat pathogens, O. yallundae and O. acuformis, differ in temperature optima for growth and fungicide resistance [START_REF] Fitt | Coexistence of related pathogen species on arable crops in space and time[END_REF]. In oilseed rape, the pathogens Leptosphaeria maculans and L. biglobosa differ in their location in stem and leaf tissues and in their maturation rates under low temperatures [START_REF] Fitt | Coexistence of related pathogen species on arable crops in space and time[END_REF]. In malaria, variation in host resistance to different Plasmodium species has been invoked to explain pathogen coexistence [START_REF] Snounou | hooThe co-existence of Plasmodium: sidelights from falciparum and vivax malaria in Thailand[END_REF]). Furthermore, some Plasmodium species are specialized to attack only the youngest or oldest red blood cells, leading to complex within-host dynamics during co-infection (McQueen & McKenzie 2004a, b). Pathogens with different transmission modes, e.g. sexual versus non-sexual transmission, can also stably coexist [START_REF] Thrall | Polymorphism in sexual vs. non-sexual transmission[END_REF].

The goal of this study was to delineate and understand the coexistence of different species of anther-smut fungi on two closely related host plants, Silene vulgaris and S. uniflora. Microbotryum is a species complex of basidiomycete fungi responsible for anthersmut disease in many plants in the Caryophyllaceae. Microbotryum anther-smut fungi are obligate pathogens that sterilize their hosts. Infected plants produce fungal teliospores in place of pollen, and female structures do not mature; in dioecious or gynodioecious species, infected female plants also develop spore-bearing anthers. Within Microbotryum, there are numerous phylogenetically divergent lineages many of which have been given species status [START_REF] Denchev | Three new species of anthericolous smut fungi on Caryophyllaceae[END_REF][START_REF] Kemler | Implications of molecular characters for the phylogeny of the Microbotryaceae (Basidiomycota : Urediniomycetes)[END_REF]Le Gac et al. 2007a). These typically show intersterility (de Vienne et al. 2009b;Le Gac et al. 2007b) and strong host specificity (de Vienne et al. 2009a). There is usually only one endemic Microbotryum species per host plant species (Le Gac et al. 2007a;[START_REF] Refrégier | Cophylogeny of the anther smut fungi and their Caryophyllaceous hosts: Prevalence of host shifts and Importance of delimiting parasite species[END_REF], although transient host shifts are not uncommon [START_REF] Antonovics | The ecology and genetics of a host shift: Microbotryum as a model system[END_REF][START_REF] Gladieux | Maintenance of fungal pathogen species that are specialized to different hosts: allopatric divergence and introgression through secondary contact[END_REF][START_REF] Hood | Karyotypic similarity identifies multiple hostshifts of a pathogenic fungus in natural populations[END_REF].

A conspicuous exception to this pattern is the presence of three divergent Microbotryum species on S. vulgaris and its closely related sister species S. uniflora (Figure 1). These are: 1) M. lagerheimii [START_REF] Denchev | Microbotryum lagerheimii sp. nov. (Microbotryaceae)[END_REF]) (MvSv1 in Le Gac et al. 2007a), 2) M. silenes-inflatae (MvSv2 in Le Gac et al. 2007a), and 3) M. violaceo-irregulare [START_REF] Kemler | Implications of molecular characters for the phylogeny of the Microbotryaceae (Basidiomycota : Urediniomycetes)[END_REF][START_REF] Lutz | Anther smuts of Caryophyllaceae: Molecular analyses reveal further new species[END_REF]. These three Microbotryum species are not particularly closely related to one-another (Figure 1), and show differences in spore morphology (color and ornamentation; Figure 1; [START_REF] Denchev | Microbotryum lagerheimii sp. nov. (Microbotryaceae)[END_REF][START_REF] Gold | Within-host competitive exclusion among species of the anther smut pathogen[END_REF][START_REF] Lutz | Anther smuts of Caryophyllaceae: molecular characters indicate host-dependent species delimitation[END_REF][START_REF] Lutz | Anther smuts of Caryophyllaceae: Molecular analyses reveal further new species[END_REF].

Microbotryum lagerheimii has lighter spore color, while verrucose spore ornamentation is unique to M. violaceo-irregulare (Figure 1; [START_REF] Lutz | Anther smuts of Caryophyllaceae: molecular characters indicate host-dependent species delimitation[END_REF][START_REF] Vánky | European smut fungi Gustav[END_REF]). All three species castrate host flowers in the same way, replacing the pollen by their spores and aborting ovaries, and are sexually transmitted (the spores are dispersed by pollinators), suggesting that it is unlikely they can have distinct ecological niches within host individuals. None of the three species has ever been found to persist in a host species other than S. vulgaris or S. uniflora, indicating that they are maintained endemically on these host plants.

One explanation for coexistence of these three fungal species could be that they are allopatric, although they have been found to co-occur in the same sites [START_REF] Abbate | Elevational disease distribution in a natural plant-pathogen system: insights from changes across host populations and climate[END_REF][START_REF] Chung | Variation in resistance to multiple pathogen species: anther smuts of Silene uniflora[END_REF]Le Gac et al. 2007a); however sampling has been limited thus far

F o r R e v i e w O n l y 6
and it remains possible that contact zones of allopatric distributions were sampled. An alternative hypothesis is that the species display differential specialization towards S. vulgaris versus S. uniflora. Previous studies further suggest that this hypothesis is unlikely, given that both host species have been found to be parasitized in nature by at least two of the Microbotryum species also parasitizing S. vulgaris: M. lagerheimii and M. silenes-inflatae [START_REF] Chung | Variation in resistance to multiple pathogen species: anther smuts of Silene uniflora[END_REF][START_REF] Smith | Anther smuts of Silene acaulis and S. uniflora in the Outer Hebrides, including an assessment of ITS genotypes of Microbotryum silenes-acaulis[END_REF].

A more plausible explanation for the species coexistence could be adaptation to different abiotic conditions. All three species affecting S. vulgaris were found almost exclusively in host populations at higher elevations in the French Alps, where overall temperatures were lower and precipitation more consistently higher than conditions experienced by lower-elevation host populations [START_REF] Abbate | Elevational disease distribution in a natural plant-pathogen system: insights from changes across host populations and climate[END_REF]). Here we investigated the generality of this pattern with a larger and broader sampling in order to further understand their coexistence. Differentiation in disjunct refugia during the last glaciation in Europe may also have shaped differential adaptation to abiotic conditions and/or generated allopatry. Strong population differentiation has been found in other Microbotryum species that reflects footprints of persistence in glaciation refugia, e.g. in M. lychnidis-dioicae parasitizing S. latifolia [START_REF] Badouin | Widespread selective sweeps throughout the genome of model plant pathogenic fungi and identification of effector candidates[END_REF][START_REF] Feurtey | Strong phylogeographic co-structure between the anther-smut fungus and its white campion host[END_REF][START_REF] Gladieux | Maintenance of fungal pathogen species that are specialized to different hosts: allopatric divergence and introgression through secondary contact[END_REF][START_REF] Vercken | Glacial refugia in pathogens: European genetic structure of anther smut pathogens on Silene latifolia and Silene dioica[END_REF] and M. silenes-acaulis parasitizing S. acaulis [START_REF] Bueker | Distribution and population structure of the anther smut Microbotryum silenes-acaulis parasitizing an arctic-alpine plant[END_REF]). High selfing rates were inferred for previously studied Microbotryum species [START_REF] Bueker | Distribution and population structure of the anther smut Microbotryum silenes-acaulis parasitizing an arctic-alpine plant[END_REF][START_REF] Giraud | Patterns of within population dispersion and mating of the fungus Microbotryum violaceum parasitising the plant Silene latifolia[END_REF][START_REF] Gladieux | Maintenance of fungal pathogen species that are specialized to different hosts: allopatric divergence and introgression through secondary contact[END_REF][START_REF] Hood | Intratetrad mating, heterozygosity, and the maintenance of deleterious alleles in Microbotryum violaceum (=Ustilago violacea)[END_REF][START_REF] Hood | Mating within the meiotic tetrad and the maintenance of genomic heterozygosity[END_REF]), and such a closed mating system may contribute to generate and maintain strong population structure and genetic isolation between species [START_REF] Gibson | Sibling competition arena: selfing and a competition arena can combine to constitute a barrier to gene flow in sympatry[END_REF][START_REF] Vercken | Glacial refugia in pathogens: European genetic structure of anther smut pathogens on Silene latifolia and Silene dioica[END_REF].

To assess the wider geographical and elevational distributions of the three different Microbotryum species occurring on S. vulgaris and S. uniflora, we collected anther-smut with microsatellite markers. We then examined whether the three Microbotryum species 1) had contrasting geographical or elevational distributions, 2) showed different frequencies on the two sibling host species S. vulgaris and S. uniflora, 3) co-occurred in the same sampling sites, and if so, whether hybrids could be detected, 4) exhibited within-species geographical population subdivision, and 5) displayed high selfing rates like other Microbotryum species.

Materials and Methods

Teliospore collection

The Microbotryum samples analyzed in this study were collected across western Europe (Supplementary Figure S1, Supplementary Table S1) from 51 sites with S. vulgaris (n = 407), and 12 sites with S. uniflora (n = 59). A sampling site was defined as a location in which samples were collected a few meters apart. Sampling was performed by numerous investigators throughout Europe and was carried out as part of general disease surveys, other field trips, or as part of more formal investigations [START_REF] Chung | Variation in resistance to multiple pathogen species: anther smuts of Silene uniflora[END_REF]Le Gac et al. 2007a, and unpublished studies). Within sampling sites, a single diseased flower was collected per individual plant, stored in silica gel, and DNA from the diploid teliospores was extracted for genetic analyses. Multiple infections by different genotypes are not uncommon in Silene-Microbotryum systems [START_REF] Bueker | Distribution and population structure of the anther smut Microbotryum silenes-acaulis parasitizing an arctic-alpine plant[END_REF][START_REF] Buono | Influence of multiple infection and relatedness on virulence: disease dynamics in an experimental plant population and its castrating parasite[END_REF][START_REF] Chung | Variation in resistance to multiple pathogen species: anther smuts of Silene uniflora[END_REF][START_REF] Hood | Dynamics of multiple infection and within-host competition by the anthersmut pathogen[END_REF]López-Villavicencio et al. 2011), but teliospores within a single flower originate from a single diploid genotype [START_REF] Gold | Within-host competitive exclusion among species of the anther smut pathogen[END_REF][START_REF] Lopez-Villavicencio | Multiple infections by the anther smut pathogen are frequent and involve related strains[END_REF]. Teliospores sampled from a given flower were therefore assumed to be from a particular pathogen individual. Because of the difficulty of identifying the different Microbotryum species in the

F o r R e v i e w O n l y
8 field, sampling was carried out without knowledge of fungal species identity, so our sample corresponds to an unbiased representation of the Microbotryum species frequencies in the sampled regions. We obtained permits for collections in national parks and all samples were collected before 2014, so not falling under the Nagoya protocol. Dried spores from each sample are freely available upon request.

Microsatellite, ITS genotyping and spore morphology

DNA was extracted and genotyped using 11 microsatellite loci as described previously (SVG1, SVG2, SVG5, SVG6, SVG8, SVG15, SN2, SN5, SN11, DC2, and DC5; Giraud 2004; [START_REF] Giraud | Isolation of 60 polymorphic microsatellite loci in EST libraries of four sibling species of the phytopathogenic fungal complex Microbotryum[END_REF]. Microsatellite scoring was performed as described previously [START_REF] Vercken | Glacial refugia in pathogens: European genetic structure of anther smut pathogens on Silene latifolia and Silene dioica[END_REF]. Analyses using GENCLONE [START_REF] Arnaud-Haond | GENCLONE: a computer program to analyse genotypic data, test for clonality and describe spatial clonal organization[END_REF] showed that the eleven markers were sufficient to discriminate multilocus genotypes in the dataset as a plateau was reached in terms of number of genotypes detected (Supplementary Figure S2).

In order to assign the genetic clusters from microsatellite analyses to the three Microbotryum species, we sequenced the ribosomal internal transcribed spacer (ITS), spanning ITS1, 5.8S and ITS2 fragments, in 27 individuals using the fungus-specific primers ITS1 and ITS4 [START_REF] White | Analysis of phylogenetic relationships by amplification and direct sequencing of ribosomal RNA genes[END_REF]) and matched these to published ITS sequences from M. silenes-inflatae, M. lagerheimii, M. violaceo-irregulare [START_REF] Kemler | Implications of molecular characters for the phylogeny of the Microbotryaceae (Basidiomycota : Urediniomycetes)[END_REF][START_REF] Lutz | Anther smuts of Caryophyllaceae: Molecular analyses reveal further new species[END_REF].

Additionally, we verified in 15 individuals spore color (light versus dark purple) under the binocular and spore ornamentation (i.e. verrucose versus reticulate) under light microscopy.

The samples chosen for ITS sequencing and for spore morphology assessment were picked at random within each cluster, and there was complete agreement in the species assignments between ITS, spore morphology and genetic clusters. A maximum likelihood phylogeny of ITS region sequenced was inferred using PhyML 3.0 [START_REF] Guindon | New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0[END_REF].

Map. The sample map (Supplementary Figure 2) was generated using the GGMAP package in R environment [START_REF] Kahle | ggmap: Spatial visualization with ggplot2[END_REF].

Descriptive statistics. Microsatellite variability was quantified by the unbiased gene diversity (H d ), allelic richness (A r ), and the fixation index (F IS ) using FSTAT 2.9.3.2 [START_REF] Goudet | FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3)[END_REF].

Observed and unbiased expected heterozygosities (H O and H E ) were computed using a custom python script. Analyses of variance (ANOVA) and Kruskall-Wallis tests were performed using JMP v7.0 (SAS Institute).

Population subdivision. We assessed population subdivision using Bayesian clustering implemented in STRUCTURE 2.3.4 [START_REF] Falush | Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies[END_REF][START_REF] Pritchard | Inference of population structure using multilocus genotype data[END_REF]. The method partitions multilocus genotypes into clusters while minimizing departure from expected frequencies and linkage equilibrium among loci. For this analysis, we used a haploid setting because Microbotryum individuals were highly homozygous (Table 1), which could bias structure inferences based on Hardy-Weinberg expectations in a diploid setting. We conducted independent runs with different number of clusters (K = 1 to 15). Each run of 400,000 iterations and 100,000 burn-in iterations was repeated ten times, using a model allowing for admixture and correlated allele frequencies. We used CLUMPP 1.1.1 (Jakobsson Because M. lagerheimii and M. violaceo-irregulare were minority components of the data set, they could not be differentiated in the overall analysis. We therefore also ran STRUCTURE on a sub-dataset encompassing only genotypes assigned to these species based on ITS sequencing. We retained for this sub-dataset the genotypes with total memberships q≥0.8 in the cluster corresponding to both M. lagerheimii and M. violaceo-irregulare in the analysis on the whole dataset (blue cluster at K=5, see results). The 0.8 threshold is arbitrary but is typically used for assigning individuals, and the genotypes with main assignments below 0.8 had assignment probabilities balanced among the different clusters, indicating they likely corresponded to individuals with low assignment power (see results).

Hybrid detection

To assess whether individuals with intermediate assignment probability to clusters were high-confidence inter-cluster hybrids or may be simply genotypes with low power of assignment, possibly carrying alleles frequent in all species, we used a two-step procedure. In the first step, we used the membership proportions of individual genomes in different clusters (q values) inferred using the admixture model of STRUCTURE to identify 'pure' representatives of the different species or of clusters within species. Based on the STRUCTURE outputs of the full dataset at K=5 (identified as the most relevant subdivision, see results), pure representatives of a given species were identified as genotypes for which the sum of q values

F o r R e v i e w O n l y
11 in the different clusters assigned to this species was above 0.8; genotypes not meeting this requirement were classified as putative hybrids. Pure representatives of the various clusters within a given species were defined as genotypes displaying q values above 0.8 within one of the clusters identified within this species; genotypes not meeting this requirement were considered as putative admixed genotypes among clusters.

In the second step, we used pure representatives of the species or clusters as 'learning samples', allowing estimation of the membership proportions for the putative hybrid samples in 'learning clusters', using STRUCTURE, options USEPOPINFO and UPDATEPFROMPOPFLAGONLY. In assignment tests of putative hybrids between M. silenes-inflatae and either M. lagerheimii or M. violaceo-irregulare, genotypes were considered as hybrids if they had q values above 0.2 and below 0.8 in learning clusters. For M. lagerheimii and M. violaceo-irregulare, genotypes were considered as hybrids if they had a membership proportion q above 0.2 and below 0.8 in the two learning clusters corresponding to M. lagerheimii and M. violaceo-irregulare for differentiating these two species. Because clustering with STRUCTURE can be strongly influenced by unbalanced sample sizes [START_REF] Kalinowski | The computer program STRUCTURE does not reliably identify the main genetic clusters within species: simulations and implications for human population structure[END_REF], and because sample size for M. silenes-inflatae was one order of magnitude higher than for the other species, datasets of pure representatives were randomly resampled ten times to the size of the smallest species sample (n=7, M. violaceoirregulare), and STRUCTURE analyses were run on each resampled dataset.

In order to check the validity of the genetic subdivision retrieved from the STRUCTURE analysis, using a method not assuming outcrossing or lack of genetic disequilibrium, we also carried out a principal component analysis (PCA) on the microsatellite data using the 'dudi.pca' function in the R-package ADE4 [START_REF] Dray | The ade4 package: implementing the duality diagram for ecologists[END_REF]. We then used a discriminant analysis of the principal components (DAPC), as implemented in the R-package [START_REF] Jombart | adegenet: a R package for the multivariate analysis of genetic markers[END_REF][START_REF] Jombart | Discriminant analysis of principal components: a new method for the analysis of genetically structured populations[END_REF]. The 30 first principal components were retained, representing 90% of total genetic variation, as well as the six first discriminant functions. The membership probabilities of genotypes were computed using DAPC for numbers of clusters K ranging from 1 to 25.

Field Code Changed Field Code Changed F o r R e v i e w O n l y 12 ADEGENET

Host specificity, geographical and climatic distributions

To assess host specificity, we examined the differential occurrence of the pathogen species and clusters on S. vulgaris and S. uniflora. We also plotted the geographical distribution of the Microbotryum species and their genetic clusters to assess their degree of allopatry. Because there was large variance in sample size at each site, with sites often having only one individual sample, we used analyses based on presence and absence of a species or cluster at a site. Analyses based on raw numbers or relative frequencies often violated model assumptions. The associations between host species and pathogen species or genetic clusters were tested using the 'CATMOD' procedure in the SAS Statistical Computing Environment (SAS Institute Inc. 2013). Associations between species or genetic clusters and geographic regions were tested using Fisher's exact tests with the SAS procedure 'FREQ' due to small sample size for some categories.

To assess environmental effects, we compared the elevational distributions of the Microbotryum species, and of the genetic clusters within M. silenes-inflatae and M. lagerheimii. In addition to elevation, we incorporated in our analysis a suite of bioclimatic variables for each of the sampling locations in order to assess whether components of climatic variation could help explain the observed distributions. Nineteen bioclimatic variables (Supplementary Table S2) were obtained from the WorldClim database (means across years ). We used a principal components analysis to analyze and reduce the variation to three major axes using the SAS procedure PRINCOMP (based on the correlation matrix). We then carried out a logistic regression to examine the relationship between species presence and the first three components generated by the PCA, using the 'glm' function in the R stats package with a quasi-binomial link.

Results

Species identification, species distributions and population subdivision

The STRUCTURE analysis on the whole dataset, modeling K=2 to 15 clusters, revealed several levels of strong genetic differentiation, with well-defined Microbotryum clusters appearing up to K=5 (i.e., clusters encompassing multiple genotypes with membership coefficients near the maximum; Figure 2 and Supplementary Figure S3). The genetic clusters were overall broadly distributed among sampled regions, although a few clusters appeared geographically more restricted (Figure 2 and Supplementary Figure S3). Some genetic differentiation was observed among samples collected on S. uniflora and S. vulgaris, with some genetic clusters being specific to either host species, but also with other genetic clusters occurring on both host species.

In order to assign the genetic clusters to the three previously described Microbotryum species, we sequenced ITS in 27 randomly chosen strains and built a maximum likelihood Supplementary Figure S4), indicating that this cluster corresponds to M. silenes-inflatae.

Analyses at K>2 only further subdivided M. silenes-inflatae into several clusters (Figure 2).

The ITS sequences falling in the fully supported clades containing the published M. lagerheimii (n=5) and M. violaceo-irregulare (n=2) sequences all belonged to genotypes in the second (blue) cluster of the K=2 STRUCTURE model (Supplementary Figure S4).

The log-likelihood plateaued at K>5 with only a single clustering solution (Supplementary Figure S5A); the additional clusters appearing at K values above 5 mainly introduced additional subdivision within M. silenes-inflatae, produced many genotypes with intermediate assignment probability and never clearly separated the genotypes assigned to M. lagerheimii and M. violaceo-irregulare using ITS into distinct clusters. K=5 was therefore considered as the most relevant K value, with the blue cluster corresponding to both M. lagerheimii and M. violaceo-irregulare, and the four other clusters (green, orange, pink and yellow) to different lineages within M. silenes-inflatae.

The species M. lagerheimii and M. violaceo-irregulare could not be separated into distinct clusters at any K values using the full dataset, despite their large phylogenetic distance (Figure 1). The lack of ability to distinguish these species is likely due to their small sample sizes compared to that collected for M. silenes-inflatae. Such difficulties in recovering species differentiation in cases of unbalanced sample sizes have been shown by simulations [START_REF] Neophytou | Bayesian clustering analyses for genetic assignment and study of hybridization in oaks: effects of asymmetric phylogenies and asymmetric sampling schemes[END_REF]. We therefore ran STRUCTURE on a sub-dataset with the (dark blue). Spore morphology (color and ornamentation) of 15 individuals, chosen at random within each species, provided further evidence for cluster assignment to species: all five strains assigned to M. silenes-inflatae had dark purple spores with reticulate ornamentation, all five strains assigned to M. lagerheimii had light purple spores with reticulate ornamentation, and all strains assigned to M. violaceo-irregulare had dark purple spores with verrucose ornamentation (Supplementary Table S3).

The principal component analysis (PCA) on the full dataset separated the three species, with M. lagerheimii and M. silenes-inflatae in separate corners of the parameter space and M. violaceo-irregulare appearing intermediate between them (Figure 4A). The ITS sequence variants particular to M. silenes-inflatae were not confined to any of its four genetic clusters (Supplementary Figure S3; Figures 4A and4B), providing evidence that they constitute within-species genetic subdivision. The M. lagerheimii and M. violaceo-irregulare species and two M. silenes-inflatae clusters were found on both S. vulgaris and S. uniflora, while the two other M. silenes-inflatae clusters occurred only on S. vulgaris (Figure 2). A PCA run on the sub-dataset of the two less frequent Microbotryum species also discriminated them well (Figures 4C and4D).

Clustering using DAPC supported the clustering patterns found in STRUCTURE analyses (Supplementary Figure S7 and Figure 2). The DAPC Bayesian information criterion reached a minimum value at K=28 (Supplementary Figure S8). However, visual inspection of

F o r R e v i e w O n l y
16 membership probability patterns across K values indicated that DAPC and STRUCTURE identified the same five main clusters (Supplementary Figure S7 and Figure 2). The DAPC on the full dataset could not distinguish M. lagerheimii and M. violaceo-irregulare into different clusters either, while the DAPC on the genotypes from these two species separately produced a pattern of membership proportions largely similar to the one obtained with STRUCTURE (Supplementary Figure S9).

Hybridization and admixture checking

We performed assignment tests using the option in STRUCTURE that uses pure reference genotypes as learning samples to determine whether individuals with intermediate assignment probabilities could simply be genotypes with low assignment power instead of hybrids or admixed genotypes. Within M. silenes-inflatae, we identified 64, 40, 75 and 79 genotypes representatives of the yellow, orange, pink and green clusters, respectively, that we used as learning samples. In assignment tests for within-species admixture, all individuals with intermediate assignment probabilities had comparable q values in multiple learning samples (q values ranging from 0.193 to 0.349, i.e. approximately one divided by the number of clusters), suggesting that all genotypes with intermediate assignment probabilities could be genotypes with low assignment power instead of genuine admixed individuals (Supplementary Table S4). Within M. lagerheimii, we identified 24 and 12 pure M. lagerheimii and M. silenes-inflatae, respectively, that we used as learning samples. In assignment tests for inter-specific hybrid identification, all genotypes had substantial membership in all learning samples considered. Therefore, no genotype could be identified as a high-confidence hybrid, neither between species nor between lineages within M. silenesinflatae (Supplementary Table S5), and they may instead represent genotypes with low assignment power. This inference was supported by the heterozygosity values of the genotypes with intermediate assignment probabilities, whose mean (H O =0.18) was not higher than that in other samples (Table 1), in contrast to what would be expected for hybrids.

Host specificity, geographical and climatic distributions

Pathogen species and several clusters within M. silenes-inflatae and M. lagerheimii were differentially distributed among host species. Because pathogens were sampled without prior knowledge of their genotype, their differential distribution was not due to sampling bias. The rarest species, M. violaceo-irregulare, was found only on S. vulgaris in the western and central Alps (Figure 5H and Supplementary Figure S10). The two other species, M. lagerheimii and M. silenes-inflatae, were both found in the Alps, the Pyrenees and the UK Figure 5H and Supplementary Figure S10). However, the two genetic clusters within M. lagerheimii were well-separated both geographically and with regard to host species, with the light blue cluster occurring only on S. uniflora in the UK and the dark blue cluster only on S. vulgaris in the Alps and the Pyrenees (Figure 5, F andG). While all Microbotryum species were found on S. vulgaris, M. silenes-inflatae was found in 90% of S. vulgaris populations but just 50% of S. uniflora populations, the difference in frequency being significant (Fisher's exact test, p = 0.004). Clusters within M. silenes-inflatae had largely overlapping

F o r R e v i e w O n l y
18 distributions on S. vulgaris, except the pink cluster which occurred only in the Alps and the Pyrenees (Figure 5 A-E). The green cluster of M. silenes-inflatae was found commonly on S. uniflora in the UK (Fig 5 B), while only single records of the yellow and orange clusters were identified there. This bias towards occurrence in S. vulgaris populations was statistically significant only for the yellow cluster (Fisher's exact test, p = 0.025). In the single site where host plants were identified as S. vulgaris var. prostrata (Site 22), the high elevation form of the host, all strains were assigned to M. violaceo-irregulare.

Different species and clusters within species often occurred in the same sampling site. Among the 45 sites with more than one assigned sample, nine displayed sympatry of Microbotryum species (M. silenes-inflatae and either M. lagerheimii or M.violaceoirregulare). Within M. silenes-inflatae, 11 sampling sites out of 32 displayed sympatry among at least two genetic clusters (Supplementary Figure S10). Genotypes with intermediate assignment probabilities in STRUCTURE were not more frequently found in sites with cooccurring clusters (χ 2 = 0.36, df = 1, p = 0.54), supporting the interpretation that they represent genotypes with low assignment power rather than hybrid individuals. Moreover, the genotypes with intermediate assignment probabilities in different species were present as frequently in sites where only a single species was found as in sites where two or more species were found (χ 2 = 1.17, df = 1, p = 0.26). Genotypes with intermediate assignment probabilities in different M. silenes-inflatae clusters were in 33 out of 40 (82%) sites with a single cluster, and the genotypes with intermediate assignment probabilities in different M. lagerheimii clusters were in 13 out of 17 (76%) sites with a single cluster (grey genotypes in Supplementary Figure S10).

Where our sampling was most extensive (Pyrenees, western Alps and central Alps), we investigated whether the distributions of the species (51 sites) and of clusters within M. 5 and Supplementary Figure 10). Within M. silenes-inflatae, the green cluster showed significant differences in occurrence among regions (Fisher's exact test, p = 0.025), being more common in the western Alps. The analysis approached significance for the orange cluster (Fisher's exact test, p = 0.052) which was more common in the Pyrenees (Figure 5 and Supplementary Figure 10). The other two clusters showed no significant differences in occurrence among regions (Fisher's exact tests, p = 0.634 for pink and p = 0.408 for yellow).

The three species occurred at significantly different elevations (Figure 6; ANOVA, df = 2, Sum of squares = 4486478, F-ratio = 8.7, p = 0.0002). The elevational distribution was not normal, but ANOVAs are robust to such assumption violations [START_REF] Lindman | Analysis of variance in complex experimental designs W. H. Freeman & Co[END_REF]) and a non-parametric test also indicated significant differences in elevation among species (Kruskall-Wallis non-parametric test, df = 2, χ 2 = 16.7, p = 0.002). The species M. violaceoirregulare was found only at high elevations, M. lagerheimii at the lowest elevations, while M. silenes-inflatae was found across a greater range of elevations (Figure 6). The two clusters within M. lagerheimii were also found at significantly different elevations (Figure 6; ANOVA, df = 1, Sum of squares = 21620792, F-ratio = 167.9, p < 0.0001; Kruskall-Wallis test, df = 1, χ 2 = 24.1, p < 0.0001), as well as the four clusters within M. silenes-inflatae (Figure 6; ANOVA, df = 3, Sum of squares = 10011015, F-ratio = 24.3, p < 0.0001; Kruskall-Wallis test, df = 2, χ 2 = 56.6, p < 0.0001). In M. lagerheimii, the contrast in elevation between S2) were reduced to three axes using PCA. All PC axis loadings were very similar for both species and cluster analyses (Supplementary Tables S2 andS6). Axis 1 (explaining 63 and 66% of the variation for species and cluster, respectively) represented decreasing temperature and increasing precipitation; axis 2 (19.3 and 19.4% of the variation) was weighted heavily on increased seasonality of precipitation and more moderate interactions between temperature and rainfall; axis 3 (8.3 and 7.6% of the variation) was weighted positively on more mild winter temperatures and lower temperature seasonality (Supplementary Tables S2 andS6).

Due to the high degree of collinearity among variables, effects of climate, elevation and geographic location were then analyzed in a second PCA that incorporated latitude, longitude, elevation, as well as scores for the first three principal components of the bioclimatic analysis. Three major axes summarized the variation, similarly for species and S7). The first axis (accounting for 66.3 and 48.1% of the variation in the species and cluster analyses, respectively) was weighted heavily on geographical position (latitude and longitude) and precipitation seasonality (PC axis 2 in the bioclimate analysis). The second axis (33.0 and 32.6% of the variation) was weighted heavily on elevation and the three climatic PC axes (especially axis 1 for temperature and rainfall) and was almost completely independent of location. Axis 3 (17.0 and 17.1%) was relatively independent of both location and elevation, but weighted strongly on climatic axes representing fewer extremes of cold and dry periods.

The logistic regression analysis on species (Supplementary Table S8) showed a significant negative association between the occurrence of M. lagerheimii and location and seasonality in precipitation (PC axis 1), suggesting that either climate, location, or both contributed to explain this species' distribution. For the logistic regression analysis on M. silenes-inflatae genetic clusters (Supplementary Table S8), occurrence of the orange cluster was significantly associated with PC axis 2 (reflecting increased elevation and its characteristic colder and wetter climate), and the yellow cluster approached significant negative association with PC3 (indicating greater chance of occurrence at sites with more extreme cold and dry periods, independent of geography or elevation) (Supplementary Table S8).

Levels of genetic variation and mating system

Summary statistics were computed for clusters identified at K=5 in analyses of the full dataset for M. silenes-inflatae, and at K=3 in analyses of the M. lagerheimii-M. violaceo-irregulare dataset (Table 1). The genotypes with intermediate assignment probabilities in STRUCTURE outputs were excluded. Low heterozygosity levels were found in all three species, leading to

F o r R e v i e w O n l y
22 high levels of inferred selfing rates (mean of 0.67, with range of 0.33-0.92, Table 1). Higher selfing rates were estimated in M. violaceo-irregulare (0.92) and M. lagerheimii (0.80-0.86) than in M. silenes-inflatae (0. 33-0.73).

Discussion

This study aimed to understand the coexistence of three distant species of anther-smut fungi, all castrating pathogens of the sister host species, Silene vulgaris and S. uniflora. The three recognized Microbotryum species that are specific to these hosts were found in our sampling, although with contrasting frequencies: M. silenes-inflatae was by far the most prevalent, followed by M. lagerheimii and then M. violaceo-irregulare. These three phylogenetically distant species, previously distinguished based on morphological characters, were unambiguously identified based on microsatellite markers complemented with ITS sequences and spore morphology. We had to analyze a restricted microsatellite dataset to discriminate between M. lagerheimii and M. violaceo-irregulare, likely because of their low abundances compared to M. silenes-inflatae [START_REF] Neophytou | Bayesian clustering analyses for genetic assignment and study of hybridization in oaks: effects of asymmetric phylogenies and asymmetric sampling schemes[END_REF]. Despite co-occurring in the same host species and sites, no high-confidence inter-specific hybrids were identified between the pathogen lineages, in agreement with the strong post-zygotic isolation generally found among distant Microbotryum species (de Vienne et al. 2009b;Le Gac et al. 2007b;[START_REF] Petit | Co-occurrence and hybridization of anther-smut pathogens specialized on Dianthus hosts[END_REF]).

However, no crossing studies have experimentally investigated post-zygotic isolation among these particular Microbotryum species, and we cannot exclude the possibility that our microsatellite markers lacked power to identify hybrids with sufficient confidence, in particular later-generation backcrossed genotypes. This suggests that the differentiation between the two M. lagerheimii clusters found on S. vulgaris and S. uniflora, respectively, may result from past separation in disjunct glacial refugia with contrasting climates rather than host adaptation. Subdivisions within M. silenes-inflatae had less obvious biological interpretation as the clusters were broadly overlapping in geography. Indeed, all four M. silenes-inflatae clusters were found in the Alps, the UK and the Pyrenees on S. vulgaris. However, only the green cluster was found regularly on S. uniflora in the UK. Additionally, this green cluster was more frequent in the western Alps than in the Pyrenees, while the opposite was found for the orange cluster, suggesting possible allopatric or host-related differentiation of the M. silenes-inflatae genetic clusters. The geographic co-occurrence of the four M. silenes-inflatae clusters in the Alps may result from recent movements and secondary contact. Footprints of a glacial refugia in Europe have been found in M. lychnidis-dioicae parasitizing S. latifolia [START_REF] Badouin | Widespread selective sweeps throughout the genome of model plant pathogenic fungi and identification of effector candidates[END_REF][START_REF] Feurtey | Strong phylogeographic co-structure between the anther-smut fungus and its white campion host[END_REF][START_REF] Gladieux | Maintenance of fungal pathogen species that are specialized to different hosts: allopatric divergence and introgression through secondary contact[END_REF][START_REF] Vercken | Glacial refugia in pathogens: European genetic structure of anther smut pathogens on Silene latifolia and Silene dioica[END_REF]) and M. silenes-acaulis parasitizing S. acaulis [START_REF] Bueker | Distribution and population structure of the anther smut Microbotryum silenes-acaulis parasitizing an arctic-alpine plant[END_REF]. The genetic clusters within M. silenes-inflatae have not homogenised despite sympatry, which is likely due to the selfing mating system of these fungi. All three Microbotryum species indeed had high estimated lagerheimii, as has been found for other Microbotryum species previously studied (e.g., those on S. latifolia, S. dioica, S. acaulis, Saponaria officinalis and Dianthus spp.; [START_REF] Bucheli | Differences in population structure of the anther smut fungus Microbotryum violaceum on two closely related host species, Silene latifolia and S. dioica[END_REF][START_REF] Bueker | Distribution and population structure of the anther smut Microbotryum silenes-acaulis parasitizing an arctic-alpine plant[END_REF][START_REF] Delmotte | Host and parasite population structure in a natural plant-pathogen system[END_REF][START_REF] Fortuna | Polymorphic microsatellite markers for the tetrapolar anther-smut fungus Microbotryum saponariae based on genome sequencing[END_REF][START_REF] Giraud | Patterns of within population dispersion and mating of the fungus Microbotryum violaceum parasitising the plant Silene latifolia[END_REF][START_REF] Gladieux | Maintenance of fungal pathogen species that are specialized to different hosts: allopatric divergence and introgression through secondary contact[END_REF][START_REF] Hood | Intratetrad mating, heterozygosity, and the maintenance of deleterious alleles in Microbotryum violaceum (=Ustilago violacea)[END_REF][START_REF] Hood | Mating within the meiotic tetrad and the maintenance of genomic heterozygosity[END_REF][START_REF] Petit | Co-occurrence and hybridization of anther-smut pathogens specialized on Dianthus hosts[END_REF]. Such high selfing rates may be the result of mating among sporidia from one diploid individual (a spore load from a single infected flower). Additionally, it may be the result of mating between sporidia from single teliospores, a process also promoted by intra-tetrad mating that is present in many Microbotryum species as a result of ordered segregation of the mating-type loci [START_REF] Hood | Intratetrad mating, heterozygosity, and the maintenance of deleterious alleles in Microbotryum violaceum (=Ustilago violacea)[END_REF][START_REF] Hood | Breaking linkage between mating compatibility factors: Tetrapolarity in Microbotryum[END_REF]. Interestingly, M. silenes-inflatae had lower inferred selfing rates than the other species, and it would be interestinginformative to investigate in future studies if this is due to ecological factors, such as more frequent codispersal of spores from different fungal individuals by pollinators, or to intrinsic factors, such as a lower propensity of intra-tetrad mating.

The three Microbotryum species had large range overlaps on S. vulgaris, with frequent co-occurrence within sampling sites, indicating that allopatry may not be a major factor currently allowing the persistence of these distinct pathogen species on the same host. The distribution of M. violaceo-irregulare and M. lagerheimii was nevertheless significantly different between the local regions in the Pyrenees and Alps, suggesting they may have a history of past allopatry. However, their current broad range overlap indicates that their coexistence must either be transient or promoted by niche specialization.

The co-existence of the three Microbotryum species is clearly not maintained by host specialization, as all three species occurred on S. vulgaris. Moreover, it has been shown that resistance of S. uniflora to M. lagerheimii and M. silenes-inflatae was significantly positively [START_REF] Chung | Variation in resistance to multiple pathogen species: anther smuts of Silene uniflora[END_REF], which indicates that the plant resistance is not specific to any one Microbotryum species. This further supports the inference that host specialization on S. uniflora versus S. vulgaris is not a major factor contributing to the cooccurrence of these Microbotryum species. However, because the host was not genotyped in our study and no infection trials were carried out, we cannot be ruled out thatthe influence of different host genotypes influenceon Microbotryum species distribution and coexistence.

The situation of the anther-smut species co-occurring on S. vulgaris is in some respects similar to that of anther-smut disease on the genus Dianthus where multiple Microbotryum species co-occur in the same sampling sites on multiple Dianthus species [START_REF] Petit | Co-occurrence and hybridization of anther-smut pathogens specialized on Dianthus hosts[END_REF]. In this case, however, the Microbotryum species on Dianthus hosts are closely related phylogenetically and hybrids, likely facilitated by the proximal relationships between the Microbotryum species, are frequently found [START_REF] Petit | Co-occurrence and hybridization of anther-smut pathogens specialized on Dianthus hosts[END_REF]. The situation on Dianthus spp. and the S. vulgaris/S. uniflora sister pair contrasts with the pattern found in all Caryophyllaceae studied to date, where each host species has been found to harbor only one endemic Microbotryum species in natural populations (Le Gac et al. 2007a;[START_REF] Refrégier | No evidence of reproductive character displacement between two sister fungal species causing anther smut disease in Silene[END_REF]. Silene vulgaris also occasionally harbors Microbotryum species as a result of transient cross-species disease transmission from other Silene hosts, such as S. latifolia and S. dioica [START_REF] Antonovics | The ecology and genetics of a host shift: Microbotryum as a model system[END_REF][START_REF] Hood | Karyotypic similarity identifies multiple hostshifts of a pathogenic fungus in natural populations[END_REF]. It is unclear if this is due to a particular susceptibility of S. vulgaris to anther smut or to its very wide geographic and environmental ranges that allow contact with the disease on many other Silene hosts.

The three anther-smut species on S. vulgaris occurred at significantly different elevations, suggesting that they may occupy different ecological niches in terms of abiotic conditions. This was supported by some association of M. lagerheimiii occurrence with climatic variables related to temperature and precipitation. Within M. silenes-inflatae, the where Site 22 is located showed that none of the three fungal species were confined to any specific host ecotype [START_REF] Abbate | Elevational disease distribution in a natural plant-pathogen system: insights from changes across host populations and climate[END_REF]. Furthermore, phenotypic traits including those that mark differences between the ecotypes change gradually over elevation, and the differentiation between populations is relative to their geographic distance [START_REF] Berardi | Elevational divergence and clinal variation in floral color and leaf chemistry in Silene vulgaris[END_REF]. The much wider sampling in the present study confirms that Microbotryum on S. vulgaris is largely confined to higher elevations (79% of the samples, 88% of the sites, were found above 1300m), even though populations of the host plant are abundant at lower elevations across its range. More generally, elevation and abiotic conditions have been shown to impact the presence of fungal pathogens in other systems [START_REF] Cordier | The composition of phyllosphere fungal assemblages of European beech (Fagus sylvatica) varies significantly along an elevation gradient[END_REF][START_REF] Desprez-Loustau | Relating germination dynamics of Melampsora pinitorqua teliospores to temperature and rainfall during overwintering[END_REF][START_REF] Enjalbert | Genetic evidence of local adaptation of wheat yellow rust (Puccinia striiformis f. sp tritici) within France[END_REF][START_REF] Gange | Rapid and recent changes in fungal fruiting patterns[END_REF][START_REF] Mboup | Genetic structure and local adaptation of European wheat yellow rust populations: the role of temperature-specific adaptation[END_REF][START_REF] Vacher | Distribution of parasitic fungal species richness: influence of climate versus host species diversity[END_REF]). This study thus begs important questions for future work with Microbotryum pathogens investigating how the genetics of adaptation to elevation in fungi interacts with the adaptation to, or of, host plants. For example, one question is whether Microbotryum is rare at low elevations because the fungus cannot adapt to low-elevation climatic conditions, or because the costs of resistance in the host to the pathogen decline under those conditions. In conclusion, using a molecular ecology approach we showed that we can assign samples to the three Microbotryum anther-smut pathogens parasitizing S. vulgaris and S. uniflora using microsatellite markers. Furthermore, we revealed the existence of four genetic clusters within M. silenes-inflatae. We also showed that the different species and genetic clusters had large range overlaps and were all found on S. vulgaris. These findings suggest that the coexistence of multiple species and genetic clusters of Microbotryum pathogens on S. vulgaris and S. uniflora cannot be explained by differences in host or geographic distributions alone, although recent secondary contacts cannot be entirely excluded. Instead, we found that the different species and genetic clusters appeared to occupy different ecological niches in terms of abiotic conditions, in particular elevation, temperature and precipitation. This study highlights the importance of including not only species differences and geography but also abiotic conditions when considering the distribution and co-existence of pathogenic fungi on similar hosts. violaceo-irregulare were determined using the K=3 output from the sub-dataset excluding M.

silenes-inflatae.

  host populations in different regions in western Europe, and genotyped them

  ) to identify potential distinct solutions among the results of independent replicate runs for each K ('Greedy' algorithm; 100 random input sequences; G' statistic), and to average individual assignment probabilities (q) over replicated runs with identical clustering solution.

  ) recorded from a network of weather stations and interpolated over a 30-arc second resolution grid (www.worldclim.org)[START_REF] Hijmans | Very high resolution interpolated climate surfaces for global land areas[END_REF]. For each sampling site, values of the 19 variables were extracted using nearest-neighbor interpolation with the Spatial Analyst Tool in ArcMap 9.2 software (ESRI (Environmental Systems Resource Institute) 2009

  including published ITS sequences from M. silenes-inflatae, M. lagerheimii, M. violaceo-irregulare and other Microbotryum species (Supplementary Figure S4; Genbank accession numbers MH491551-MH491577). Out of the 20 sequences in the fully supported clade containing the published M. silenes-inflatae ITS sequences, 19 were genotypes assigned to the largest STRUCTURE cluster in analyses assuming K=2 clusters (pink cluster;

  M. lagerheimii and M. violaceo-irregulare only. At K=3 in this sub-dataset (Figure 3; Supplementary Figures S5B and S6), one cluster (brown) corresponded to M. violaceo-irregulare genotypes while the other two clusters (dark and light blue) represented M. lagerheimii genotypes, mainly separating isolates collected on S. uniflora in the UK and Netherlands (light blue) versus genotypes collected on S. vulgaris in the Pyrenees and Alps

  representatives of the dark and light blue clusters, respectively; the only individual with intermediate assignment probabilities in the M. lagerheimii full dataset (q=0.465 in the blue cluster; q=0.535 in the light blue cluster) had high membership in the blue cluster (q=0.731) in the M. violaceo-irregulare/M. lagerheimii subset analysis.

  38 and 400 pure genotype representatives of M. violaceo-irregulare,

  (37 sites with unambiguous assignments) were associated with local geography, elevation and/or climate. The occurrence of M. silenes-inflatae did not differ between the Pyrenees, western Alps or central Alps (Fisher's exact test, p = 0.168). The occurrence of M. violaceo-irregulare and M. lagerheimii differed significantly among regions (Fisher's exact tests, p = 0.031 and p = 0.001 respectively), with the former being absent in the Pyrenees and the latter absent in the Western Alps (Figure

  was associated with the differences in host species, as described above, i.e. coastal S. uniflora in the UK versus broad elevational distribution of S. vulgaris elsewhere.The Figure2also illustrates the differences in elevation between species and clusters within each geographical region: clusters remain well separated when genotypes are sorted by elevation, in particular at K=5. Differences in Microbotryum distributions among both elevations and geographic regions could indicate impacts of climate. To test whether species or clusters were associated with different elevations and/or climates among the most heavily sampled regions (S. vulgaris populations in the Pyrenees, western Alps, and central Alps), we carried out a two-step analysis to first reduce climatic variables to uncorrelated principal components (PCs), and then test for the relative ability of climate, elevation and geographical location to explain species and cluster distributions. For both species and cluster datasets, the 19 bioclimatic variables (Supplementary Table

  subdivision within species, with two clusters in M. lagerheimii and four clusters in M. silenes-inflatae. Given their weak genetic differentiation (as shown by ITS sequences and distribution in the PCA), the genetic clusters within M. silenes-inflatae likely do not represent further cryptic species. Within M. lagerheimii, the genetic subdivision corresponded to separation in geography, host and environmental conditions and it remains to be determined if these might be cryptic species. Silene vulgaris and S. uniflora were allopatric during the last glaciation[START_REF] Prentice | Mosaic variation in allozyme and plastid DNA markers in the European ranges of Silene vulgaris and its partially sympatric relative S. uniflora (Caryophyllaceae)[END_REF], with S. uniflora likely persisting in northern refugia from the Baltic/Scandinavian region and S. vulgaris in southern European refugia.

  being highly homozygous, especially M. violaceo-irregulare and M.

  some of the genetic clusters also varied with elevation and climatic factors. These effects may reflect environmental tolerances of the fungal species and clusters, or ecotypic differences among host populations. It is well known that the host S. vulgaris has different ecotypes adapted to contrasting elevations (Marsden-Jones & Turrill 1957), and these may have different anther-smut specificities. In the present study, just one host population was identified as S. vulgaris var. prostrata, the high elevation form of the host (Site 22, in which all infections were caused by M. violaceo-irregulare). However, a previous study of Microbotryum on S. vulgaris in the small region of the south-eastern French Alps

  grant awarded to Benjamin Roche. JA acknowledges NSF grant DEB-1115765 and NIH grant R01GM122061. Sequencing was funded by the ''Consortium National de Recherche en Génomique'' and the ''service de systématique moléculaire'' of the Muséum National d'Histoire Naturelle (CNRS IFR 101). It is part of the agreement n_2005/ between the Genoscope, the Muséum National d'Histoire Naturelle the INRA and the CEA (Genoscope), on the projects ''Macrophylogeny of life'', "SPEED-ID" and "Bibliothèque du Vivant" directed by Guillaume Lecointre and Jean-Yves Rasplus.

F

Figure 1 .

 1 Figure LegendsFigure 1. Maximum-likelihood tree showing the phylogenetic placement of Microbotryum lagerheimii (HQ832090), M. violaceo-irregulare (AY588104), and M. silenes-inflatae (JN223404), all parasitizing Silene vulgaris, among other Microbotryum species, based on internal transcribed spacer (ITS) sequences. Microbotryum from Polygonum bistorta was used as an outgroup. Bootstraps are indicated. Pictures show differences in teliospore color between the dark purple of M. silenes-inflatae and M. violaceo-irregulare versus the lighter lavender of M. lagerheimii, as well as differences in teliospore ornamentation (verrucose in M. violaceo-irregulare, fully reticulated in M. silenes-inflatae and M. lagerheimii) under 1000x light-microscope magnification.

Figure 2 .

 2 Figure 2. Proportions of ancestry in K (from 2 to 10) clusters of Microbotryum spp. genotypes inferred with the STRUCTURE program. Each genotype is represented by a vertical bar, partitioned into K segments representing the amount of ancestry of its genome in K clusters. When several clustering solutions (''modes'') were found within replicate runs, only the major mode is shown with its corresponding proportion of runs. IS: Iceland, NL: Netherlands, UK: United Kingdom, CH: Switzerland, IT: Italy, F: France. For each region, genotypes are ordered by sampling elevation (represented at bottom). See Supplementary Figure S3 for a sorting by membership coefficient.

Figure 3 .

 3 Figure 3. Proportions of ancestry in K (from 2 to 5) clusters of the Microbotryum violaceoirregulare and M. lagerheimii genotypes inferred with the STRUCTURE program, using a sub-dataset excluding M. silenes-inflatae. In the sub-dataset we kept only the genotypes assigned to the blue cluster at K=5 in the output of the Structure analysis on the whole dataset (Figure 4). A single clustering solution (''mode'') was found among replicate runs. Each genotype is represented by a vertical bar, partitioned into K segments representing the amount

Figure 4 .

 4 Figure 4. Principal component analysis (PCA) on multilocus microsatellite genotypes for the dataset including the three Microbotryum species (A: principal component 1 vs principal component 2; B: principal component 1 vs principal component 3) for the sub-dataset including only M. lagerheimii and M. violaceo-irregulare (C: principal component 1 vs principal component 2; D: principal component 1 vs principal component 3). Scatterplots for the first three principal components are shown using a colour labelling of genotypes defined according to the assignment of multilocus genotypes to three species and six clusters within species using Bayesian clustering analyses. Colors indicate clusters as in Figures 2 and 3.Black outlines indicate the samples for which the ITS region has been sequenced.

Figure 5 .

 5 Figure 5. Distribution maps of the different genetic clusters and species of anther-smut Microbotryum fungi on Silene vulgaris and S. uniflora, and maps of their co-occurrence. Pie charts are proportional to the number of samples. Colors correspond to those used in other figures.

Figure 6 .

 6 Figure 6. Violin plots showing elevational distribution of species and clusters within species. All samples belonging to Microbotryum silenes-inflatae at K=5 (including all genotypes with low assignment power to the green, orange, pink and yellow M. silenes-inflatae clusters) are included in the leftmost M. silenes-inflatae column. Otherwise, only samples with highconfidence cluster assignment are shown. Samples belonging to M. lagerheimii and M.

Table 1 .

 1 Summary statistics of genetic variation and estimates of selfing rates in the identified Microbotryum species and clusters. , number of pure representatives of the species/clusters, as determined using Bayesian clustering analyses; H d , gene diversity; A r , allelic richness, corrected for sample size using rarefaction (standardized sample size of two diploid individuals); H O /H E , observed/expected heterozygosity; F IS , fixation index; s, selfing rate estimated as 2F IS /(1+F IS ).

	o
	r
	R
	e
	v i
	e
	w
	O n
	l
	y

N

Acknowledgements

We thank the following people for help with strain collection: Mickael Le Gac, Jacqui Shykoff, Anne Thomas, Michel Bartoli, Alex Widmer, Lucie Salvaudon, Erhart, Marcel Dorken, Marc-André Selosse, Rémy Relli, Jean-Pierre Thoirain, Fabienne van Rossum, Amanda Gibson, Alexis Matamoro, and the French 'Parc National des Pyrénées'. We thank Ricardo Rodriguez de la Vega for help with figures. T.G. acknowledges the grant ANR-12-ADAP-0009 (Gandalf project). M.E.H. acknowledges the grant NSF-DEB 0747222 and 1115765. J.L.A. acknowledges the National Science Foundation award DEB-0910212, ANR grant 10-PDOC-017-01 'SPATEVOLEPID' awarded to Sébastien Lion, and ANR JC grant F o r R e v i e w O n l y (A) Yellow cluster of M. silenes-in atae (B) Green cluster of M. silenes-in atae (C) Pink cluster of M. silenes-in atae (D) Orange cluster of M. silenes-in atae (E) Coexistence of M. silenes-in atae clusters M. silenes-inflatae M. lagerheimii M. violaceo-irregulare M. silenes-inflatae (E) Blue cluster of M. lagerheimii (F) Light blue cluster of M. lagerheimii (G) M. violaceo-irregulare (H) Coexistence of M. silenes-in atae and M. lagerheimii (I) Coexistence of M. silenes-in atae and M. violaceo-irregulare

Data accessibility

The genotypes and elevation data are available on Dryad : doi:10.5061/dryad.627fm05.

GenBank accession numbers for ITS sequences: MH491551-MH491577.

The genotypes and elevation data will be submitted to Dryad upon acceptance.


Martina Sandroni

Favereau Ludovic

Aurélien Planchat

Huriye Akdas-Kilig

Nadine Szuwarski

Yann Pellegrin

Errol Blart

Hubert Le Bozec

Mohammed Boujtita

Heteroleptic copper(I

Since 1991 1 and the discovery of DSC (Grätzel cells), many attempts to replace the costly and toxic (albeit remarkably efficient) rutheniumpolypyridine complexes have been reported. 2,3 Copper(I)-bis(diimine) 15 complexes have early shown promising results in this field. 4,5,6 Lately, the use of heteroleptic copper(I) complexes has afforded significant PCEs thanks to an improved extinction coefficient in the visible and electron transfer vectorialization. 6 The latter point is an essential criterion to fulfil in the design of efficient sensitizers for TiO2. Indeed, each ligand is set to 20 play one (or more) well-defined role such as anchoring, passivation of the surface and assisting charge injection. Accordingly, ligands differ by their molecular structures and therefore by their electronic natures. In the course of our program on heteroleptic bis-diimine copper(I) complexes, 7,8,9 prepared according to the HETPHEN concept developed by Schmittel 25 and colleagues, 10 we have prepared and studied four new stable heteroleptic copper(I) complexes [CuL 0 L n ] + hereafter named Cn (n = 1-4, Figure 1). The anchoring ligand L 0 (6,6'-dimesityl-2,2'-bipyridine-4,4'dicarboxylic acid) is based on the classical 4,4'-dicarboxylic acid bipyridine onto which were attached two mesityl groups in positions 6 and 6', providing the necessary steric bulk to avoid the formation of homoleptic complexes. The ligands completing the coordination sphere of 35 the copper(I) ion belong to the family of 4,4'-bis(styrylphenyl)-2,2'bipyridines, derivatized with electron releasing moieties of different strength. Methyl groups in α of the chelating nitrogen atoms confer rigidity to the final scaffold, preserving the excited state from exciplex quenching and excessive flattening upon excitation, to a certain extent. 40 Three complexes C2, C3 and C4, bearing respectively alkoxy, N,Ndiethylamine and N,N-diphenylamine moieties were thus isolated. For the sake of comparison, a fourth model complex [CuL 0 L 1 ] + (C1) was synthesized, with L 1 = 2,2',4,4'-tetramethylbipyridine.

The syntheses of all ligands are reported in ESI. The 45 HETPHEN modus operandi was used to isolate C1-4 and started with the synthesis of the Cu(L 0 )] + intermediate in DMF. An equivalent of L n was subsequently added dropwise, entailing an immediate colour change of the medium from yellow to deep red. Impurities were removed by size exclusion chromatography. A similar protocol was used to isolate the 50 dimethyl-ester forms of each complex (named hereafter Cnester, n=1-4, synthesis given in ESI).

The electronic absorption spectra of the complexes were recorded in solution and on nanocrytalline TiO2 films (Figures 2 andS3). All the complexes featured the classical MLCT absorption band at ca. 500 nm 55 (Table 1 and Figure 2). 11 The increased conjugation of the π system on both L 0 and L n (n = 2-4) induces a stabilization of the π* orbitals, explaining the red-shift of this transition compared to the benchmark bisneocuproine Cu(I) complex C5 (Figure S7). 9 One notices that the MLCT bands are more intense as well, because of the increased ground state 60 dipolar moment generated by the combination of electron poor L 0 and electron rich L n . The complexes C3 and C4 present higher light harvesting efficiency in the visible than C1 and C2 because of an intense additional intraligand charge transfer transition (ILCT), located at the edge of the visible around 420 nm. 9

65 λ (nm) [ε (M -1 .cm -1 ] E (V)* [ΔE (mV)] C1 477 [4.7•10 3 ] 0.94 [96] ─ C2 500 [9.8•10 3 ] 0.91 [96] ─ C3 504 [1.3•10 4 ] 1.08 a [-] 0.80 a [-] C4 502 [1.4•10 4 ] 1.03 b [-] 0.95 b [-]
Table 1. UV-Visible and electrochemical data for C1-4. *data collected with the methyl ester forms of C1-4.

This very intense ILCT transition corresponds to a shift of the electron 70 density from the electron rich amine moieties to the electron poor pyridine. Such band does exist for C2 too, but is significantly blue-shifted compared to C3 and C4 because of the poorer electron donating power of L 2 . Spectra recorded on TiO2 transparent electrodes (Figure S3) feature the same patterns than those recorded in solution phase (Figure 2). 75 Overall, the complexes displayed a rather broad and intense absorption over a large wavelength frame (λonset~620 nm), revealing their potentials No luminescence was detected upon excitation in the MLCT band, regardless the conditions. This could be due to cis-trans isomerization of the vinyl double bond 12 or the lesser rigidity of the Cu-bpy coordination 20 cage compared to Cu-Phen, facilitating the deleterious exciplex quenching.

To record better resolved cyclic and pulsed voltammograms (no adsorption of the dye on the electrode), these measurements were performed on the diester forms of complexes. The latter featured the 25 expected, reversible copper-centred oxidation around 1 V vs. SCE (Table 1). The voltammograms of C3ester and C4ester displayed an additional oxidation wave at 0.95 and 0.80 V vs. SCE respectively, corresponding to the removal of one electron from the NR2 (R = ethyl or phenyl) amine moieties. Only differential pulse voltammetry allowed discriminating the 30 two close oxidation steps for complex C4ester. The higher Cu II /Cu I potentials displayed by C3ester and C4ester likely originate in the coulombic repulsion between the copper cation and the electrogenerated hole on the amine fragment.

The combination of electrochemical and UV-Vis data allowed 35 evaluating the Gibbs energies associated to the various charge transfer processes. In all cases, both charge injection and dye regeneration are exergonic (ca. 300 meV, see Table S1 in SI). Energy-wise, C1-4 feature roughly the same behaviours. TiO2 electrodes were dipped while still hot for two days in ethanolic solutions of C1-4 and the photovoltaic devices 40 were then assembled with a platinum counter-electrode, sealed with a hotmelt polymer frame and their performances along those of the reference benchmark N719 were evaluated under AM 1.5 calibrated artificial sunlight ( Spin coating a 0.1M CDCA (chenodeoxycholic acid) ethanolic solution onto the photo-electrodes prior to the final sealing is anticipated 70 to eliminate the deleterious self-quenching process induced by aggregation. Rewardingly, unprecedented improvements in the power conversion efficiencies (PCEs) of all DSCs were observed upon such CDCA-surface treatment. First of all, an increase of the photopotential was observed for all DSCs. C1-4 based devices exhibited a 50-60 mV rise 75 of the Voc, together with a decrease of the dark current. This improvement was therefore assigned to the higher electron concentration in the CB and to a passivation of recombination sites by the co-adsorbent molecules. In the case of C2, the octyl chains may provide a built-in, efficient protection of titanium dioxide's surface, thus explaining the 80 lesser increase of the Voc (ca. 30 mV). The electron lifetime (τn) and mean transit time (τtr) of photoinjected charge carriers were then recorded by intensity-modulated photovoltage spectroscopy (IMVS) and intensitymodulated photocurrent spectroscopy (IMPS). However, these measurements reveal that no significant improvement of both τn and τd 85 was observed when CDCA was added in the preparation of the series of solar cells (see SI). As a result, the charge collection efficiency ηcoll measured as a function of the illumination intensity is quite similar for all the dyes (see SI). the Jsc for all devices. One calculates a 70% increase of the photocurrent for C1 and C4 and 45 and 50% for C2 and C3 based photovoltaic devices, respectively. For a better understanding, IPCEs of DSCs with and without CDCA treatment were compared (Figure S5). A significant increase of the IPCE is observed for CDCA-treated cells, regardless the 5 dye (Figure 3), together with a broadening of the signals. Several reasons can be invoked to rationalize this important result. First, CDCA molecules release protons and this bends the conduction band downwards, increasing thus the driving force of the electron injection into the semi-conductor, and consequently improving the electron injection 10 yield. 13 Second, a noticeable increase of the absorbance of C1-4 based photo-electrodes was monitored upon CDCA treatment, along with a slightly broadened MLCT transition (Figure S4). These subtle changes in the absorption spectra of the chemisorbed complexes are in line with the IPCE, and are probably grounded in a reorganization of the dye 15 monolayer upon CDCA adsorption. The role of CDCA is often associated with the disruption of dye aggregates and certainly comes into play here, especially due to the presence of organic styryl branches on the complexes C2-4. Based on the effect of CDCA on both IPCE measurements and on the Jsc enhancement, we conclude that the main 20 role of CDCA with these complexes is certainly to decrease the aggregation on TiO2 surface leading to higher LHE and injection quantum yield. In these conditions, DSCs provided a maximum PCE of 4.66% for C3-based device. This is to date the highest PCE ever reported for a DSC based on a copper(I) complex sensitizer, and hold great promises for the 25 future of these cheap solar cells. Most highly performing dyes, including ruthenium complexes, are neutral species, while these first series of copper(I) complexes are positively charged. This is certainly one weak point of these dyes, which can be overcome by using new ancillary ligands 30

Conclusions

We successfully isolated four stable heteroleptic copper(I)polypyridine complexes, using the HETPHEN concept. Through a careful choice of ligands, unprecedented PCE were measured, reaching 4.66%. The new anchoring ligand L 0 paves the route to prepare other sensitizers 35 as it certainly forms stable heteroleptic copper(I) with many unhindered diimine ligands. This contribution brings further credit to these molecular complexes as efficient sensitizers for DSCs, en route for a cheap and less toxic substitute to ruthenium dyes.

A c c e

  p t e d m a n u s c r i t p t Heteroleptic copper(I)-polypyridine complexes as efficient sensitizers for dye sensitized solar cells Martina Sandroni, ac Ludovic Favereau, a Aurelien Planchat, a Huriye Akdas-Kilig, b Nadine Szuwarski, a Yann Pellegrin, a Errol Blart, a Hubert Le Bozec, b Mohammed Boujtita* a and Fabrice Odobel* a The synthesis and the physico-chemical characterizations of HETPHEN based heteroleptic copper(I)-bis(diimine) complexes are reported. In TiO 2 based dye sensitized solar cells (DSCs), the latter display impressive photoconversion 10 efficiencies (PCEs), unprecedented for first row transition metal coordination complexes.

Figure 1 .

 1 Figure 1. Molecular structures of L n and Cn (n=1-4) 30

1 A 10 Figure 2 .

 1102 Figure 2. UV-Visible spectra of complexes C1 (dash-dot), C2 (plain), C3 15

  photocurrent and photovoltage. The latter is assigned to a lower light harvesting efficiency (LHE) and probably to an exacerbated charge recombination with the electrolyte. Indeed, the positive charge of C1-4 entails a coulombic repulsion between them on the surface of TiO2, increasing the number of unoccupied adsorption sites and thus 55

Figure 3 .

 3 Figure 3. IPCE for DSCs sensitized with C1 (dash-dot), C2 (dot), C3 (plain), and C4 (dash) recorded with CDCA.

A c c e p t e d m a n u s c r i t p t

  40

Table 2

 2

	45
		Voc (mV)	Jsc (mA.cm -2 )	FF (%)	PCE (%)
	C1 a	475	2.20	72.80	0.76
	C2 a	535	2.89	72.54	1.12
	C3 a	545	7.51	71.52	2.93
	C4 a	565	6.70	73.32	2.77
	N719	635	16.87	68.69	7.36
	C1 b	525	3.76	74.64	1.47
	C2 b	565	4.99	72.39	2.04
	C3 b	605	10.86	70.97	4.66
	C4 b	625	10.13	69.76	4.42

and SI for details).

Table 2 .

 2 Photovoltaic data for DSCs

Acknowledgements

ANR agency is gratefully acknowledged for the financial support of these researches through program HeteroCop (n° ANR-09-BLAN-0183-01).

Notes and references