Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Polymorphic Microsatellite Markers for the Tetrapolar Anther-Smut Fungus Microbotryum saponariae Based on Genome Sequencing

  • Taiadjana M. Fortuna ,

    taiadjanafortuna@gmail.com

    Affiliation Laboratoire d’Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, 91400, Orsay, France

  • Alodie Snirc,

    Affiliation Laboratoire d’Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, 91400, Orsay, France

  • Hélène Badouin,

    Affiliation Laboratoire d’Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, 91400, Orsay, France

  • Jérome Gouzy,

    Affiliations INRA, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR441, Castanet-Tolosan, F-31326, France, CNRS, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR2594, Castanet-Tolosan, F-31326, France

  • Sophie Siguenza,

    Affiliations INRA, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR441, Castanet-Tolosan, F-31326, France, CNRS, Laboratoire des Interactions Plantes-Microorganismes (LIPM), UMR2594, Castanet-Tolosan, F-31326, France

  • Diane Esquerre,

    Affiliation GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet-Tolosan, F-31326, France

  • Stéphanie Le Prieur,

    Affiliation Laboratoire d’Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, 91400, Orsay, France

  • Jacqui A. Shykoff,

    Affiliation Laboratoire d’Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, 91400, Orsay, France

  • Tatiana Giraud

    Affiliation Laboratoire d’Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, 91400, Orsay, France

Abstract

Background

Anther-smut fungi belonging to the genus Microbotryum sterilize their host plants by aborting ovaries and replacing pollen by fungal spores. Sibling Microbotryum species are highly specialized on their host plants and they have been widely used as models for studies of ecology and evolution of plant pathogenic fungi. However, most studies have focused, so far, on M. lychnidis-dioicae that parasitizes the white campion Silene latifolia. Microbotryum saponariae, parasitizing mainly Saponaria officinalis, is an interesting anther-smut fungus, since it belongs to a tetrapolar lineage (i.e., with two independently segregating mating-type loci), while most of the anther-smut Microbotryum fungi are bipolar (i.e., with a single mating-type locus). Saponaria officinalis is a widespread long-lived perennial plant species with multiple flowering stems, which makes its anther-smut pathogen a good model for studying phylogeography and within-host multiple infections.

Principal Findings

Here, based on a generated genome sequence of M. saponariae we developed 6 multiplexes with a total of 22 polymorphic microsatellite markers using an inexpensive and efficient method. We scored these markers in fungal individuals collected from 97 populations across Europe, and found that the number of their alleles ranged from 2 to 11, and their expected heterozygosity from 0.01 to 0.58. Cross-species amplification was examined using nine other Microbotryum species parasitizing hosts belonging to Silene, Dianthus and Knautia genera. All loci were successfully amplified in at least two other Microbotryum species.

Significance

These newly developed markers will provide insights into the population genetic structure and the occurrence of within-host multiple infections of M. saponariae. In addition, the draft genome of M. saponariae, as well as one of the described markers will be useful resources for studying the evolution of the breeding systems in the genus Microbotryum and the evolution of specialization onto different plant species.

Introduction

Anther-smut fungi belonging to the genus Microbotryum sterilize their host plant by aborting ovaries and replacing pollen by fungal spores. Sibling Microbotryum species are highly specialized on their respective host plants and show strong post-mating reproductive isolation [13]. Microbotryum fungi have become model species in ecology and evolution of plant pathogenic fungi [4,5]. Although anther-smut disease occurs in hundreds of plant species in the Caryophyllaceae [6] and at least nine other plant families [7,8], most studies to date have focused on M. lychnidis-dioicae, which parasitizes the white campion Silene latifolia. This is the first fungal species in which distinct mating types were described [9], it has been used in classical genetic analyses for decades (e.g., [1012]), and more recently as a model in disease transmission and metapopulation dynamics (e.g., [13,14]).

More recently, population genetics studies using microsatellite markers allowed elucidation of the genetic population structure of M. lychnidis-dioicae in Europe [1517], patterns of post-glaciation expansion [18], and the invasion history of M. lychnidis-dioicae in North America [5,19]. The use of microsatellite markers also revealed a high prevalence of multiple infections by different fungal genotypes, which segregated in different stems of the same Si. latifolia host [20]. The co-infecting genotypes also revealed a high level of relatedness, suggesting competitive exclusion between unrelated fungal genotypes [2123]. In addition, a reference genome of M. lychnidis-dioicae is now available [24,25], which provides opportunities for examining genomic changes within the genus. In particular, in this species, the two independent mating-type loci typical of basidiomycetes are trapped within a very large region of suppressed recombination on what have thus become the sex chromosomes, spanning ca. 90% of the chromosome length. These genomic data are consistent with previous findings based on marker segregation analyses and optical maps [2628]. This large non-recombining region explains why M. lychnidis-dioicae is bipolar, with a single mating-type locus. Bipolarity is the rule for most Microbotryum species studied so far, suggesting similar linkage between mating-type loci [29].

Many resources are therefore available for the species M. lychnidis-dioicae, while the other anther-smut fungi remain poorly studied. Other Microbotryum species, however, deserve further investigations [7,30], and in particular M. saponariae, parasitizing Saponaria ocymoides and Sa. officinalis [2,8]. The common soapwort Sa. officinalis, with a broad distribution in Europe, is a large long-lived perennial plant with multiple stems, which makes it a particularly interesting subject for the study of multiple infections (i.e., multiple pathogen genotypes coexisting within a host) and phylogeography. In addition, M. saponariae has been recently identified as one of the few tetrapolar species (i.e., maintaining the independence of the two mating-type loci), recently identified, within the anther-smut Microbotryum fungi [29]. Of the available microsatellite markers for various Microbotryum fungi, very few cross-amplify on M. saponariae [15,17,31], and only four were polymorphic. Therefore, our aim here was to develop polymorphic microsatellite markers by generating a draft genome of M. saponariae using mate-pair Illumina sequencing. A similar approach to develop SNP markers has been used in other plant pathogenic fungi for studying multiple infections [32].

Microsatellite loci have been widely used as genetic markers due to their ubiquity, reproducibility, neutrality and high level of polymorphism [33]. These features make them valuable molecular markers for inference of evolutionary and demographic parameters [34]. However, for each new species, microsatellites often need to be isolated de novo and the characterization process of each locus can be time consuming and expensive given the need for numerous, locus-specific fluorescent primers. Recent methods using more simple and inexpensive approaches allow the isolation of polymorphic markers for large-scale multiplex assays. Here we isolated microsatellite markers using a high resolution and rapid approach with generic vector-specific fluorescent primers and unlabelled locus-specific primers [35], which makes the costs considerably lower.

The isolated microsatellite markers will allow studies for unravelling patterns of infection (i.e. single or multiple genotypes) within host, population structure and phylogeography of M. saponariae. In addition, the molecular and genomic resources obtained here will contribute to increasing our knowledge on the evolution of reproductive systems and host specialization of anther-smut fungi.

Material and Methods

Genome Sequencing

The Microbotryum saponariae genome sequenced in this study was originated from a fungal individual (strain 1053) collected on a host plant Sa. officinalis growing on the University Paris-Sud campus (Orsay, France) in 2012 and is available upon request. Diploid teliospores produced in the anther of a diseased plant were collected and cultivated on a nutrient-rich medium, as Potato Dextrose Agar (PDA) at 23°C for a few days. On nutritive media, fungal teliospores germinate and undergoes meiosis, each producing four haploid products, called sporidia, that replicate clonally. The sporidia harvested from PDA plates represent therefore the haploid products of many independent meiosis events of a single diploid fungal genotype. After adequate growth on PDA, sporidia were collected and stored in silica gel at -20°C until use. Fungal DNA was extracted using a hydraulic press (Carver, Inc., Wabash, USA) to lyse the cell walls and the Genomic-Tip 100/G kit (Qiagen®, CA, USA) following the manufacturer’s protocol. DNA purity was assessed by measuring 260/280 and 260/230 absorbance ratios with a NanoDrop® spectrophotometer (2000; Thermo Scientific, Wilmington, USA), and DNA concentration was measured using a Qubit® fluorometer (2.0; Thermo Scientific). Preparation of DNA libraries was performed by the GeT-Place INRA platform, Toulouse, France. Mate-pair libraries with insert sizes of 3kb were prepared with Illumina NexteraTM Mate Pair Sample Prep kits, and sequencing was performed on a HiSeq 2000 IlluminaTM sequencer at a depth of coverage of 100X on average. The genome sequence of M. saponariae is available in GenBank under the accession number PRJEB11435.

Detection of Microsatellite Loci in the Genome and Primer Design

Simple repeats were detected with RepeatMasker (v4.0.5; [36]) on the assembled genome. Only perfect di- and tri-nucleotide repeats of 8 to 12 repetitions were considered. Flanking regions of 200 bp were selected on each side of the nucleotide repeat. Thus, repeats located less than 200 bp from the contig extremity were discarded. Primers were designed using Primer3 [37] on 100 to 210 bp long fragments containing the target microsatellite sequences (S1 Table). For each microsatellite locus, up to five primer pairs were mapped to the target genome with GMAP (version 2014-10-16; [38]) using default parameters, and only pairs mapping to a single locus were further considered.

Polymorphism Screening of Microsatellite Markers

For a first screening of polymorphism, a sample subset of eleven M. saponariae individuals were selected from different geographical locations across Europe (Fig 1, S2 Table), and DNA was extracted with the NucleoSpin® Soil kit (Macherey-Nagel, Germany). To detect polymorphic microsatellite markers, we used a simple and inexpensive approach that avoids generating fluorescent labels for all tested markers [35]. The eleven DNA samples, each at 5 ng.μL-1, were pooled as a template DNA for downstream applications. PCR amplifications were performed for each microsatellite marker with unlabelled specific primers using a high fidelity Taq DNA Polymerase Taq Pfu® (Promega, USA). We used a PCR touchdown program with 94°C for 5 min, 35 cycles of 94°C for 30 s, 62°C for 30°C with a decrease of 1°C at each cycle during 12 cycles and 50°C during 23 cycles, 72°C for 45 s, followed by 72°C for 5 min. Successful PCR amplifications were confirmed by 2% (w/v) agarose gel electrophoresis. PCR products for multiple individual microsatellite loci, amplified from the same DNA pool, were combined and purified using NucleoSpin® 96 PCR Clean-up (Macherey-Nagel, Germany), before ligation into pDrive vector (Qiagen®, USA). Ligation products were diluted 1/10 with deionised water and amplified again for each marker using the specific ‘reverse’ microsatellite primer and a universal fluorescently labelled primer, M13 forward (-40; 5’ GTTTTCCCAGTCACGAC 3’), targeting the plasmid insert flanking region. We used the same PCR touchdown program as described above. Labelled amplicons were diluted 1/225 and fractionated by capillary electrophoresis on an ABI PRISM X3730XL at the GENTYANE INRA platform, Clermont, France. Alleles of each microsatellite locus were identified and their size scored using GENEMAPPER 5.0 (Applied BiosystemsTM, Foster city, USA).

thumbnail
Fig 1. Map of the geographical location of Microbotryum saponariae individuals sampled in 97 populations of the host plants: 1) Saponaria officinalis (red circles, n = 92), including those used for the first polymorphism screening (green circles, n = 11), and 2) Saponaria ocymoides (blue circles, n = 5) in Europe.

Fungal DNA was extracted from a sporulated anther of a diseased plant and was used to test the polymorphism of the 22 microsatellite markers.

https://doi.org/10.1371/journal.pone.0165656.g001

PCR Multiplex Assembling

For the assignment of PCR multiplexes we defined four amplification size ranges: 1) 100 to 115 bp; 2) 125 to 150 bp; 3) 155 to 170 bp; and 4) 180 to 210 bp. These groups were defined to separate the microsatellite markers by size and thus avoid interference between markers within each multiplex. For each multiplex, compatibility of the primer pairs for each locus was tested with the software Multiplex Manager (1.2; [39]). As only four different colours of dye are available we set the maximum number of loci per reaction to four and a complementary threshold equal to seven (default parameter). Primer pairs were then successfully combined into six multiplexes, ranging from three to four markers per multiplex (Table 1).

thumbnail
Table 1. Summary statistics for the 22 microsatellite markers in Microbotryum saponariae (N = 97).

Repeat: Motifs and repeat numbers of the microsatellite markers in the individual genome from which the markers were developed. Range: allelic size range in number of base pairs (bp). Multiplex: number indicating to which multiplex group the marker belongs. Multiplex 3 a and b were kept separate during PCR amplification to avoid the amplification interference that occurred when all four markers were amplified together, but they were pooled for genotyping in the end. Dye: forward primers of the selected markers were labelled with the following fluorescent dyes: Yakima Yellow equivalent to VIC, ATTO 550 equivalent to NED, ATTO 565 equivalent to PET. Alleles: number of alleles observed in the sample. HO: observed heterozygosity. HE: expected heterozygosity; Symbols (*) and (+) show significant deficit or excess in heterozygosity compared to Hardy-Weinberg expectations (* P < 0.05; ** P < 0.01; *** P < 0.001). FIS: Fixation index. SE: standard error.

https://doi.org/10.1371/journal.pone.0165656.t001

Polymorphism of the Markers in Natural Populations

We further tested the polymorphism of the 22 microsatellite markers retained after the first screening described above, using 97 fungal individuals collected from 92 populations of Sa. officinalis and 5 populations of Sa. ocymoides across Europe (Fig 1, S2 Table). In addition, 24 individuals of a population of M. saponariae (population 1188, in Bures-sur-Yvette near the University Paris-Sud Campus, S2 Table) were also genotyped in order to assess within-population deviation from Hardy-Weinberg equilibrium. Diploid teliospores were harvested from one anther of a diseased sporulating flower. Fungal DNA was extracted in 100 μL using Chelex (Bio-Rad®, USA) protocol previously described [40], and ten-fold diluted for PCR amplification. Fungal samples were genotyped using the 6 multiplexes of microsatellite markers developed above. PCR reactions were performed separately for each multiplex in 12,5 μL volume containing 3 μL of DNA, ¼ of the recommended volume of multiplex PCR Kit (Qiagen®, USA), and 1.25 μL of the primer mix. The primer mix included 2 μM of unlabelled forward and reverse primer and 0.5 μM or 0.75 μM of labelled forward primer depending on the dye label. FAM and Yakima Yellow dyes were at 0.75 μM, while ATTO 550 and ATTO 565 dyes were at 0.5 uM. The PCR cycling program used was the same as described above. After checking for successful PCR amplifications on 2% (w/v) agarose gel electrophoresis, PCR amplicons were diluted and fractionated by capillary electrophoresis at the GENTYANE INRA platform as previously described. Alleles of each locus were scored using GENEMAPPER 5.0 (Applied BiosystemsTM, USA).

Cross-Species Amplification

To test the microsatellite markers developed for M. saponariae on other Microbotryum species, we sampled fungal individuals (S2 Table) and extracted DNA using the methods described for testing marker polymorphism in natural populations. The following species were tested: M. lychnidis-dioicae parasitizing Si. latifolia, M. silenes-dioicae parasitizing Si. dioicae, M. silene-acaulis parasitizing Si. acaulis, M. lagerheimii parasitizing Si. vulgaris, M. violaceum s.s. parasitizing Si. nutans, M. violaceum s.l. parasitizing Si. caroliniana, M. violaceum s.l. parasitizing Si. flos-cuculi, M. shykoffianum parasitizing Dianthus pavonius, M. carthusianorum parasitizing D. superbus, M. dianthorum parasitizing D. seguieri and M. scabiosae parasitizing Knautia arvensis (S2 Table). The success of amplification was checked on 2% (w/v) agarose gel electrophoresis and PCR amplicons were submitted to capillary electrophoresis at GENTYANE platform. For each locus, alleles were scored as previously described.

Descriptive Statistics

The observed and expected heterozygosity (HO and HE), departure from Hardy-Weinberg expectations (FIS), linkage disequilibrium (LD) among locus pairs and genotypic disequilibrium were calculated using GENEPOP [41] on the web at http://genepop.curtin.edu.au using the sample of one individual per population. Furthermore, FIS was calculated within the population of M. saponariae for which all individuals were genotyped. The nominal P-value of 0.05 was adjusted for multiple comparisons using a Bonferroni correction (i.e., 2.3x10-4 based on 219 tests).

Results and Discussion

The search for microsatellite loci in the M. saponariae draft genome yielded 188 loci meeting our criteria. We chose 96 loci located on different scaffolds and with different amplification sizes for designing multiplexes, with half of di- and half of tri-nucleotide repeats. Out of the 96 tested markers, four yielded no amplification. We thus screened six multiplexes of three to four microsatellite markers (Table 1, S1 Table) for polymorphism based on pooled DNA of a sample subset. We selected among those showing the highest level of polymorphism and the most easily readable peaks on the capillary electrophoresis chromatograms (S1 Fig).

We further tested the polymorphism of the markers retained after the first screening described above. A total of 22 microsatellite markers were found polymorphic and easy to score (Table 1). The number of alleles per locus varied from 2 to 11, for a total of 51 alleles observed. Descriptive statistics on the polymorphism of M. saponariae and on deviations from Hardy-Weinberg equilibrium are shown in Table 1. The mean observed (HO) and mean expected heterozygosity (HE) values ranged from 0 to 0.97 and from 0.01 to 0.58, respectively. Most markers (20 out of 22, i.e., 91%) displayed significantly lower levels of heterozygosity than expected under Hardy-Weinberg equilibrium. The fixation index (FIS) ranged from -0.95 to 1 for a mean multilocus value of 0.61 (±0.10). Microbotryum saponariae exhibited only 11% of heterozygous genotypes while 21% were expected under Hardy-Weinberg equilibrium, which suggests high selfing rates as previously reported in other Microbotryum species [16,40,42]. Only one marker, out of the 22, Msap_146, displayed higher levels of heterozygosity than expected (extreme FIS value of -0.95) with most individuals being heterozygous. Further analysis on the marker position in the M. saponariae genome revealed that it is located on the contig carrying the pheromone receptor gene. Most likely this marker is situated in the non-recombining region of the mating-type chromosome, which leads to permanent heterozygosity and high differentiation between the two mating types [25,29]. This marker will be useful for studies on mating types and mating system, but should not be used in analyses assuming Hardy-Weinberg equilibrium. Deficits in heterozygotes may also results from a Wahlund effect, which would be expected because we genotyped here only one individual per population and populations were sampled at a large geographical scale. If allele frequencies vary among populations, our pooled sample would necessarily have a deficit of heterozygotes [43]. In fact, the genotyping of a whole population showed little deviation from Hardy-Weinberg expectations (mean FIS across loci, excluding SAP_146, of -0.03). On the other hand, despite a likely Wahlund effect, the heterozygote deficit observed here was lower than those previously reported within populations of M. lychnidis-dioicae (e.g., mean 0.96±0.06 in M. lychnidis-dioicae; [19]), suggesting lower selfing rates in M. saponariae than in M. lychnidis-dioicae. This is consistent with M. lychnidis-dioicae being bipolar and M. saponariae tetrapolar [44]. Bipolarity is indeed thought to evolve under selfing mating systems while tetrapolarity is beneficial under outcrossing, increasing discrimination against self [44].

After Bonferroni correction, significant linkage disequilibrium was observed only between six pairs of markers out of the 219 comparisons: Msap_11 and Msap_94, Msap_11 and Msap_31, Msap_4 and Msap_96, Msap_4 and Msap_27, Msap_96 and Msap_27, Msap_90 and Msap_77. This linkage disequilibrium may be due to genetic linkage between these markers or to population structure. Indeed, we pooled individuals from multiple populations and differences in alleles frequencies between populations will lead to apparent linkage disequilibrium when analysed as a single population.

Cross-species amplification was further examined among nine other Microbotryum species and was successful at 6 to 19 loci (Table 2). Some markers were amplified in all species, as for instance Msap_146 and Msap_73, while others, e.g. Msap_11 and Msap_4, were only amplified in the species most closely related phylogenetically to M. saponariae, such as the Microbotryum species on Dianthus hosts [3]. In general, the size of the alleles of each marker was different in the other species, except for nine markers where the alleles had identical size compared to those identified in M. saponariae. In some species, such as M. scabiosae, few markers could be amplified, unsurprisingly, given the high genetic distance from M. saponariae species [8]. Microsatellite markers in fungi are known to be difficult to amplify in other species, even closely related ones [45]. Among Microbotryum species, the divergence between the sister species M. lychnidis-dioicae and M. silenes-dioicae has been estimated to be 420,000 years old [46]. All the other species pairs studied here had older divergence times, as shown by phylogenies [3,7,30,46,47].

thumbnail
Table 2. Results of cross-species amplification of the 22 microsatellite loci in nine additional Microbotryum species.

No amplification (-); successful amplification is given by the allelic profile of homozygous (colour scale) or heterozygous individuals (white). The name M. violaceum s.l. refers to the species complex name and it is indicated when the fungal species has not yet been described. The colour scale illustrates the allele size: darker colour corresponds to larger amplification products.

https://doi.org/10.1371/journal.pone.0165656.t002

Conclusion

The method used for screening polymorphism at microsatellite loci [35] was efficient and powerful. Using it enabled us to isolate 22 polymorphic markers grouped in 6 multiplexes without needing to develop the 96 specific fluorescent primers. The developed microsatellite markers will be very useful for future studies on the population structure and phylogeography of M. saponariae and for comparison with patterns in M. lychnidis-dioicae [18,40,48]. Furthermore, knowing the genome sequence and population structure of M. saponariae will allow the comparison of polymorphism at genes evolving under particular selection regime, such as the mating-type genes [49,50]. Microsatellites are also rapidly evolving markers that may reveal cryptic genetic subdivisions among populations parasitizing the two closely related host species Sa. officinalis and Sa. ocymoides, whose relationship remains unclear. It will also be interesting to assess whether M. saponariae is more outbreeding than the other anther-smut fungi studied so far [51]. Tetrapolar species are indeed usually outcrossing [44], and M. saponariae may reveal the first known tetrapolar selfing species [29].

The microsatellite markers will also be useful for studying the occurrence of multiple infections in Sa. officinalis, for assessing whether the results found in Si. latifolia [21,22] and in Si. acaulis [34] are general. In particular, we are interested in knowing whether the large long-lived Sa. officinalis plants carry multiple genotypes of M. saponariae, if they segregate spatially in different stems, and if genotypes co-occurring within the same host plant are significantly related.

In addition, the genomic resources released here will be useful for studies of fungal comparative genomics of anther-smut fungi and even more broadly. Fungi are indeed very useful models for tackling questions of genomic evolution and adaptation [52]. The published genome will help in particular studying mating-type chromosome evolution in the Microbotryum genus [25,26,29] and host specialization [1,30] in plant pathogenic fungi.

Supporting Information

S1 Fig. Capillary electrophoresis chromatograms with peaks of heterozygous individuals of Microbotryum saponariae at the polymorphic microsatellite loci 102 (Msap_102).

Heterozygous individuals have the following allelic profile: a) 173 and 187; b) 173 and 196; c) 173 and 184.

https://doi.org/10.1371/journal.pone.0165656.s001

(PDF)

S1 Table. Template sequence used to design the primers for the 22 microsatellite markers in Microbotryum saponariae.

Microsatellite repeats are highlighted in bold blue.

https://doi.org/10.1371/journal.pone.0165656.s002

(PDF)

S2 Table. Information on the Microbotryum individuals sampled from different host plants across Europe and analysed in this study.

Individual fungal DNA was extracted using Chelex protocol (Biorad, USA) and was used to test the polymorphism of the 22 microsatellite markers in M. saponariae and in the other Microbotryum species used for cross-species amplification. Symbols (*) refers to DNA samples that were also extracted with the Nucleospin Soil kit (Macherey-Nagel, Germany) and were pooled as a template DNA for downstream applications in the first screening of markers polymorphism.

https://doi.org/10.1371/journal.pone.0165656.s003

(PDF)

Acknowledgments

Concerning the M. saponariae samples, we would like to thank all collectors that helped us building such a large collection. We thank the GeT-Place INRA platform for sequencing and the GENTYANE INRA platform for assistance with fungal genotyping. We thank Antoine Branca for help in statistical analyses and map elaboration. The manuscript benefited from helpful comments from Michael Hood.

Author Contributions

  1. Conceptualization: TG JAS.
  2. Formal analysis: HB JG SS TMF AS.
  3. Funding acquisition: TG JAS.
  4. Investigation: HB DE SLP TMF AS.
  5. Methodology: TG JAS.
  6. Project administration: TG JAS.
  7. Resources: TG JAS.
  8. Supervision: TG JAS.
  9. Validation: TMF AS.
  10. Visualization: TMF.
  11. Writing – original draft: TMF AS TG.
  12. Writing – review & editing: TMF HB TG JAS.

References

  1. 1. Le Gac M, Hood ME, Giraud T (2007) Evolution of reproductive isolation within a parasitic fungal species complex. Evolution 61: 1781–1787. pmid:17598756
  2. 2. Lutz M, Göker M, Piatek M, Kemler M, Begerow D, Oberwinkler F (2005) Anther smuts of Caryophyllaceae: molecular characters indicate host-dependent species delimitation. Mycological progress 4: 225–238.
  3. 3. Refrégier G, Le Gac M, Jabbour F, Widmer A, Shykoff JA, Yockteng R, et al. (2008) Cophylogeny of the anther smut fungi and their caryophyllaceous hosts: Prevalence of host shifts and importance of delimiting parasite species for inferring cospeciation. BMC Evol Biol 8: 1–17.
  4. 4. Antonovics J, Hood M, Partain J (2002) The ecology and genetics of a host shift: Microbotryum as a model system. Am Nat 160: S40–S53. pmid:18707452
  5. 5. Gladieux P, Feurtey A, Hood ME, Snirc A, Clavel J, Dutech C, et al. (2015) The population biology of fungal invasions. Mol Ecol 24: 1969–1986. pmid:25469955
  6. 6. Hood ME, Mena-Alí JI, Gibson AK, Oxelman B, Giraud T, Yockteng R, et al. (2010) Distribution of the anther-smut pathogen Microbotryum on species of the Caryophyllaceae. New Phytol 187: 217–229. pmid:20406409
  7. 7. Kemler M, Lutz M, Göker M, Oberwinkler F, Begerow D (2009) Hidden diversity in the non-caryophyllaceous plant-parasitic member of Microbotryum (Pucciniomycotina: Microbotryales). Systematics and Biodiversity 7: 297–306.
  8. 8. Kemler M, Göker M, Oberwinkler F, Begerow D (2006) Implications of molecular characters for the phylogeny of the Microbotryaceae (Basidiomycota: Urediniomycetes). BMC Evol Biol 6: 1–15.
  9. 9. Kniep H (1919) Untersuchungen über den Antherenbrand (Ustilago violacea Pers.). Ein Beitrag zum Sexualitätsproblem. Z Bot 11: 257–284.
  10. 10. Garber E, Ruddat M (2002) Transmission genetics of Microbotryum violaceum (Ustilago violacea): A case history. Adv Appl Microbiol 51: 107–127. pmid:12236055
  11. 11. Day AW, Jones JK (1968) The production and characteristics of diploids in Ustilago violacea. Genet Res 1968 11: 63.
  12. 12. Day AW (1976) Somatic nuclear division in the sporidia of Ustilago violacea. III. Ultrastructural observations. Can J Microbiol 22: 495. pmid:1260542
  13. 13. Thrall PH, Biere A, Uyenoyama MK (1995) Frequency-dependent disease transmission and the dynamics of the Silene-Ustilago host-pathogen system. Am Nat 73: S1249–S1258.
  14. 14. Thrall PH, Burdon JJ (1997) Host-pathogen dynamics in a metapopulation context: The ecological and evolutionary consequences. J Ecol 85: 743–753.
  15. 15. Giraud T, Fournier E, Vautrin D, Solignac M, Shykoff JA (2002) Isolation of 44 polymorphic microsatellite loci in three host races of the phytopathogenic fungus Microbotryum violaceum. Mol Ecol Notes 2: 142–146.
  16. 16. Bucheli E, Gautschi B, Shykoff JA (2001) Differences in population structure of the anther smut fungus Microbotryum violaceum on two closely related host species, Silene latifolia and S. dioica. Mol Ecol 10: 285–294. pmid:11298945
  17. 17. Bucheli E, Gautschi B, Shykoff JA (2000) Host-specific differentiation in the anther smut fungus Microbotryum violaceum as revealed by microsatellites. J Evol Biol 13: 188–198.
  18. 18. Vercken E, Fontaine MC, Gladieux P, Hood ME, Jonot O, Giraud T (2010) Glacial refugia in pathogens: European genetic structure of anther smut pathogens on Silene latifolia and Silene dioica. PLoS Pathog 6: e1001229. pmid:21187901
  19. 19. Fontaine MC, Gladieux P, Hood ME, Giraud T (2013) History of the invasion of the anther smut pathogen on Silene latifolia in North America. New Phytol 198: 946–956. pmid:23406496
  20. 20. Hood ME (2003) Dynamics of multiple infection and within‐host competition by the anther‐amut pathogen. Am Nat 162: 122–133. pmid:12856241
  21. 21. López-Villavicencio M, Jonot O, Coantic A, Hood ME, Enjalbert J, Giraud T (2007) Multiple infections by the anther smut pathogen are frequent and involve related strains. PLoS Pathog 3: e176. pmid:18020704
  22. 22. López-Villavicencio M, Courjol F, Gibson AK, Hood ME, Jonot O, Shykoff JA, et al. (2011) Competition, cooperation among kin, and virulence in multiple infections. Evolution 65: 1357–1366. pmid:21121914
  23. 23. Koskella B, Giraud T, Hood ME, Bronstein AEJL, Whitlock EMC (2006) Pathogen relatedness affects the prevalence of within‐host competition. Am Nat 168: 121–126. pmid:16874619
  24. 24. Perlin MH, Amselem J, Fontanillas E, Toh SS, Chen Z, Goldberg J, et al. (2015) Sex and parasites: genomic and transcriptomic analysis of Microbotryum lychnidis-dioicae, the biotrophic and plant-castrating anther smut fungus. BMC Genomics 16: 461. pmid:26076695
  25. 25. Badouin H, Hood ME, Gouzy J, Aguileta G, Siguenza S, Perlin MH, et al. (2015) Chaos of rearrangements in the mating-type chromosomes of the anther-smut fungus Microbotryum lychnidis-dioicae. Genetics 200: 1275–1284. pmid:26044594
  26. 26. Fontanillas E, Hood ME, Badouin H, Petit E, Barbe V, Gouzy J, et al. (2015) Degeneration of the nonrecombining regions in the mating-type chromosomes of the anther-smut fungi. Mol Biol Evol 32: 928–943. pmid:25534033
  27. 27. Hood ME (2002) Dimorphic mating-type chromosomes in the fungus Microbotryum violaceum. Gentics 160: 457–461.
  28. 28. Hood ME, Petit E, Giraud T (2013) Extensive divergence between mating-type chromosomes of the anther-smut fungus. Genetics 193: 309–315. pmid:23150606
  29. 29. Hood ME, Scott M, Hwang M (2015) Breaking linkage between mating compatibility factors: Tetrapolarity in Microbotryum. Evolution 69: 2561–2572. pmid:26339889
  30. 30. Lutz M, Piątek M, Kemler M, Chlebicki A, Oberwinkler F (2008) Anther smuts of Caryophyllaceae: Molecular analyses reveal further new species. Mycol Res 112: 1280–1296. pmid:18951773
  31. 31. Giraud T, Yockteng R, Marthey S, Chiapello H, Jonot O, Lopez-Villavicencio M, et al. (2008) Isolation of 60 polymorphic microsatellite loci in EST libraries of four sibling species of the phytopathogenic fungal complex Microbotryum. Mol Ecol Resour 8: 387–392. pmid:21585800
  32. 32. Tollenaere C, Susi H, Nokso-Koivisto J, Koskinen P, Tack A, Auvinen P, et al. (2012) SNP design from 454 sequencing of Podosphaera plantaginis transcriptome reveals a genetically diverse pathogen metapopulation with high levels of mixed-genotype infection. PLoS ONE 7: e52492. pmid:23300684
  33. 33. Jarne P, Lagoda PJL (1996) Microsatellites, from molecules to populations and back. Trends Ecol Evol 11: 424–429. pmid:21237902
  34. 34. Büker B, Eberlein C, Gladieux P, Schaefer A, Snirc A, Bennett DJ, et al. (2016) Distribution and population structure of the anther smut Microbotryum silenes-acaulis parasitizing an arctic–alpine plant. Mol Ecol 25: 811–824. pmid:26671732
  35. 35. Cryer NC, Butler DR, Wilkinson MJ. (2005) High throughput, high resolution selection of polymorphic microsatellite loci for multiplex analysis. Plant Methods 1: 3–5. pmid:16270924
  36. 36. Smit A, Hubley R, Green P. RepeatMasker Open-4.0. 2013–2015. Available: http://www.repeatmasker.org. 2015.
  37. 37. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. (2012) Primer3—new capabilities and interfaces. Nucleic Acids Res 40: e115. pmid:22730293
  38. 38. Wu TD, Watanabe CK (2005) GMAP: A Genomic Mapping and Alignment Program for mRNA and EST Sequences. Bioinformatics 21: 1859–75. pmid:15728110
  39. 39. Holleley C, Geerts P (2009) Multiplex Manager 1.0: a cross-platform computer program that plans and optimizes multiplex PCR. BioTechniques 46: 511–517. pmid:19594450
  40. 40. Giraud T (2004) Patterns of within population dispersal and mating of the fungus Microbotryum violaceum parasitising the plant Silene latifolia. Heredity 93: 559–565. pmid:15292913
  41. 41. Raymond M, Rousset F (1995) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered 86: 248–249.
  42. 42. Giraud T, Jonot O, Shykoff JA (2005) Selfing propensity under choice conditions in a parasitic fungus, Microbotryum violaceum, and parameters influencing infection success in artificial inoculations. Int J Plant Sci 166: 649–657.
  43. 43. Hartl DL, Clark AG. Principles of population genetics, 4th ed. Sunderland, Massachusetts: Sinauer Associates, Inc. Publishers; 2007.
  44. 44. Nieuwenhuis BPS, Billiard S, Vuilleumier S, Petit E, Hood ME, Giraud T (2013) Evolution of uni- and bifactorial sexual compatibility systems in fungi. Heredity 111: 445–455. pmid:23838688
  45. 45. Dutech C, Enjalbert J, Fournier E, Delmotte F, Barrès B, Carlier J, et al. (2007) Challenges of microsatellite isolation in fungi. Fungal Genet Biol 44: 933–949. pmid:17659989
  46. 46. Gladieux P, Vercken E, Fontaine MC, Hood ME, Jonot O, Couloux A, Giraud T (2011) Maintenance of fungal pathogen species that are specialized to different hosts: allopatric divergence and introgression through secondary contact. Mol Biol Evol 28: 459–71. pmid:20837605
  47. 47. Le Gac M, Hood ME, Fournier E, Giraud T (2007) Phylogenetic evidence of host-specific cryptic species in the anther smut fungus. Evolution 61: 15–26. pmid:17300424
  48. 48. Gladieux P, Heitman J, Giraud T (2010) Epidemiology and evolution of fungal pathogens, in plants and animals. Genet Evol Infect Dis: 59–106.
  49. 49. Xu L, Petit E, Hood ME (2016) Variation in mate-recognition pheromones of the fungal genus Microbotryum. Heredity 116: 44–51. pmid:26306729
  50. 50. Devier B, Aguileta G, Hood ME, Giraud T (2009) Ancient trans-specific polymorphism at pheromone receptor genes in basidiomycetes. Genetics 181: 209–223. pmid:19001292
  51. 51. Hood ME, Antonovics J (2000) Intratetrad mating, heterozygosity, and the maintenance of deleterious alleles in Microbotryum violaceum (= Ustilago violacea). Heredity 85: 231–241. pmid:11012726
  52. 52. Gladieux P, Ropars J, Badouin H, Branca A, Aguileta G, Vienne DM, et al. (2014) Fungal evolutionary genomics provides insight into the mechanisms of adaptive divergence in eukaryotes. Mol Ecol 23: 753–773. pmid:24341913