Development and characterization of nuclear microsatellite markers in the endophytic fungus Epichloë festucae (Clavicipitaceae)1

• Premise of the study: Microsatellite primers were developed for the endophytic fungus Epichloë festucae, which is symbiotic with Festuca rubra, to study the population genetics of the species and to compare population structures between E. festucae and its host F. rubra. • Methods and Results: We developed 14 polymorphic markers using the unplaced genomic scaffold sequences of E. festucae from GenBank. The number of alleles per locus (A) varied from four to 16, and unbiased haploid diversity (h) was 0.717 in eight populations located in the Faroe Islands, Finland, and Spain. The Spanish populations possessed a higher number of alleles and haploid diversity (on average A = 5.1 and h = 0.591, respectively) compared to northern populations (on average A = 1.5 and h = 0.199, respectively). • Conclusions: These polymorphic markers will be used by grass breeders for uses including the improvement of commercial turfgrass cultivars, and by population geneticists to study different species of the Epichloë genus.


Ap Applicati tions ons in in Pl Plant t Scien Sciences ces
Endophytic fungi are very common and important components in plant microbiomes. They are capable of infecting their host plant's tissues without causing obvious symptoms to the host ( Hyde and Soytong, 2008 ). The Epichloë endophytes (Ascomycota: Clavicipitaceae) are one of the most studied systems of plant-endophyte associations because they form symbiotic relationships with several economically important turf and forage cool season grasses ( Clay and Schardl, 2002 ). These endophytes are important agents infl uencing the growth and persistence of host grasses, and a genetically compatible endophyte infection has been demonstrated to provide a selective advantage to the host ( Ahlholm et al., 2002 ;Saikkonen et al., 2006Saikkonen et al., , 2010a. The Epichloë endophytes are a group of fi lamentous fungi that comprise sexual Epichloë species and their asexual derivatives, Neotyphodium species, which have been recently classifi ed as Epichloë ( Leuchtmann et al., 2014 ). Epichloë festucae Leuchtm., Schardl & Siegel is a fungal epichloid systemic endophyte , which systematically and intercellularly colonizes aboveground tissues and seeds of Festuca rubra L. by means of haploid hyphae. Festuca rubra is a perennial grass with rapid expansion worldwide in a wide range of ecosystems, and it is one of the most important turfgrasses in temperate regions ( Inda et al., 2008 ). Previously, nuclear microsatellite markers have been developed for Epichloë species by Moon et al. (1999) and expressed sequence tag (EST)-derived simple sequence repeats (SSRs) for pasture grass endophytes by van Zijll de Jong et al. (2003) , and among those markers, four polymorphic nuclear microsatellite markers have been used in a population genetic study on E. festucae ( Wäli et al., 2007 ). Because highly polymorphic genomic microsatellites are effective tools for studying population genetic characteristics, in this study, we aim to develop additional polymorphic microsatellite primers for E. festucae .

METHODS AND RESULTS
The unplaced genomic scaffold sequences of E. festucae were downloaded (GenBank accession no. JH158803-JH158837) and searched for ≥ 10 monoand dinucleotide repeats, and for ≥ 8 tri-, tetra-, penta-, and hexanucleotide repeats by using MSATCOMMANDER ( Faircloth, 2008 ). The selected marker regions possessed the maximum length of the repeat motif, a minimum distance of 100,000 bp between the repeat motifs within the same accession, and the presence of appropriate fl anking sequences for primer design with the following criteria: primer length of 18-27 bp, GC content 40-60%, annealing temperature 55-58 ° C, and expected amplicon size of 100-300 bp. Twenty-four primer pairs were designed with Primer3 software ( Rozen and Skaletsky, 2000 ). The forward primers were labeled with fl uorescent dyes for automated electrophoresis, and the primers were obtained from Oligomer Oy (Helsinki, Finland).
To isolate E. festucae from endophyte-infected (E+) plants of F. rubra , three leaves were collected from each selected tiller from pots containing replanted F. rubra in the greenhouse of Ruissalo Botanical Garden, Finland. The plant material was kept at 4 ° C for 24 h, followed by surface sterilization, including a treatment in 75% ethanol for 30 s, 4% sodium hypochlorite for 3 min, and 75% http://www.bioone.org/loi/apps ethanol for 15 s. Then, a leaf was cut into fi ve segments and planted on autoclaved Petri dishes containing 5% potato dextrose agar (PDA). Agar plates were stored at room temperature until mycelium emerged from the plated leaf fragments, after which a small sample of mycelium was transferred to a new PDA plate on a piece of sterile cellophane (9 cm in diameter). Epichloë festucae identity was determined based on the morphological characteristics observed in cultures, and based on the sequences of the ITS1, 5.8S rRNA, and ITS2 region of a set of isolates from different locations, compared with GenBank resources (http://www.ncbi.nlm.nih.gov) using BLAST searches. PCR amplifi cation of the ITS1, 5.8S rRNA, and ITS2 region was performed using primers ITS1 and ITS5 (http://www.fungalbarcoding.org/). Epichloë festucae is the only systemic fungus described for F. rubra , and the systemic endophyte has always been E. festucae when we have previously sequenced fungi from F. rubra (e.g., Wäli et al., 2007 ). The risk that the mycelium growing on PDA corresponds to another related fungal endophyte different from E. festucae is marginal, especially as grasses are commonly assumed to be infected by only one systemic fungus. Replanted agar plates were stored at room temperature until the growth of the mycelium was suffi cient for DNA extraction. Mycelium growth was    Kit (Omega Bio-Tek, Norcross, Georgia, USA), and a NanoDrop Lite spectrophotometer (Thermo Fisher Scientifi c, Wilmington, Delaware, USA) was used to reveal the yield and purity of DNA. The PCR reactions were performed with a single microsatellite primer pair in a 10-μ L reaction mixture containing 5-10 ng genomic DNA, 1 × GoTaq Flexi Buffer, 1.0 mM MgCl 2 solution, 0.2 mM of each dNTP, 0.2 μ M of each primer, and 1.25 units GoTaq G2 HotStart Polymerase (Promega Corporation, Madison, Wisconsin, USA). PCR reactions were performed in a C1000 Thermal Cycler (Bio-Rad, Applied Biosystems, Foster City, California, USA) as follows: an initial denaturation at 95 ° C for 2 min; followed by 30 cycles at 95 ° C for 30 s, 57 ° C for 30 s, and 73 ° C for 30 s; and a fi nal extension at 73 ° C for 5 min. Each PCR product was amplifi ed singly. The amplifi cation success was controlled with a set of PCR products using 2% agarose gels (SeaKem LE Agarose; Lonza, Rockland, Maine, USA). The products were run on an ABI 3130xl Genetic Analyzer using the GeneScan 500 ROX Size Standard (Applied Biosystems) at the Institute of Biotechnology, University of Helsinki, Finland, and assigned to allelic sizes with Peak Scanner version 1 software (Applied Biosystems). The unbiased haploid diversity ( h ) and the number of alleles ( A ) per locus and population were calculated using GenAlEx version 6.5 Smouse, 2006 , 2012 ).
Initially, 24 individuals originating from populations in the Faroe Islands ( N = 8), Finland ( N = 8), and Spain ( N = 8) were screened to reveal the competence of the 24 primer pairs. Nine primer pairs out of 24 were rejected from the further analysis because of unclear patterns with multiple bands and allelic dropouts, whereas 15 primer pairs amplifi ed reliably and produced clearly interpretable single bands, and these were used in the further analyses ( Table 1 ) . Fifteen loci were screened for polymorphism using 70 individuals originating from four different populations in the Faroe Islands, one population from Finland, and three populations from Spain ( Tables 1 and 2 ) . One locus was monomorphic while 14 loci revealed polymorphism with altogether 123 alleles ( Table 1 ). The number of alleles per locus varied from four to 16 at the species level and from one to seven at the population level ( Table 2 ). The unbiased haploid diversity per locus varied from 0.435 to 0.874 at the species level, and from 0.000 to 0.933 at the population level ( Table 2 ). The Spanish populations possessed a considerably higher number of alleles and haploid diversity (on average A = 5.1 and h = 0.591, respectively) compared to northern populations (on average A = 1.5 and h = 0.199, respectively).

CONCLUSIONS
Because endophytes have both scientifi c relevance and applied importance, these new polymorphic microsatellite markers will be useful for grass breeders, e.g., to improve commercial turfgrass cultivars of F. rubra , and for researchers to study different aspects of grass endophyte evolution. These markers presumably cross-amplify within the genus Epichloë , which includes host-specifi c endophytes of several important forage grasses ( Leuchtmann et al., 2014 ).