Next Article in Journal
The HMGA1 Pseudogene 7 Induces miR-483 and miR-675 Upregulation by Activating Egr1 through a ceRNA Mechanism
Previous Article in Journal
Identification of the Caprine Keratin-Associated Protein 20-2 (KAP20-2) Gene and Its Effect on Cashmere Traits
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 8
Issue 11
10.3390/genes8110325
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of Novel Polymorphic EST-SSR Markers in Bailinggu (Pleurotus tuoliensis) for Crossbreeding

by
Yueting Dai
1,
Wenying Su
1,
Chentao Yang
2,
Bing Song
1,
Yu Li
1,* and
Yongping Fu
1,*
1
Engineering Research Center of Chinese Ministry of Education for Edible and Medicinal Fungi, Jilin Agricultural University, Changchun 130118, China
2
China National GeneBank, Environmental Genomics, Beijing Genomics Institute, Shenzhen 518083, China
*
Authors to whom correspondence should be addressed.
Genes 2017, 8(11), 325; https://doi.org/10.3390/genes8110325
Submission received: 16 October 2017 / Revised: 16 October 2017 / Accepted: 8 November 2017 / Published: 17 November 2017
(This article belongs to the Section Microbial Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Identification of monokaryons and their mating types and discrimination of hybrid offspring are key steps for the crossbreeding of Pleurotus tuoliensis (Bailinggu). However, conventional crossbreeding methods are troublesome and time consuming. Using RNA-seq technology, we developed new expressed sequence tag-simple sequence repeat (EST-SSR) markers for Bailinggu to easily and rapidly identify monokaryons and their mating types, genetic diversity and hybrid offspring. We identified 1110 potential EST-based SSR loci from a newly-sequenced Bailinggu transcriptome and then randomly selected 100 EST-SSRs for further validation. Results showed that 39, 43 and 34 novel EST-SSR markers successfully identified monokaryons from their parent dikaryons, differentiated two different mating types and discriminated F1 and F2 hybrid offspring, respectively. Furthermore, a total of 86 alleles were detected in 37 monokaryons using 18 highly informative EST-SSRs. The observed number of alleles per locus ranged from three to seven. Cluster analysis revealed that these monokaryons have a relatively high level of genetic diversity. Transfer rates of the EST-SSRs in the monokaryons of closely-related species Pleurotus eryngii var. ferulae and Pleurotus ostreatus were 72% and 64%, respectively. Therefore, our study provides new SSR markers and an efficient method to enhance the crossbreeding of Bailinggu and closely-related species.

1. Introduction

Pleurotus tuoliensis [1], also known as Bailinggu, is one of the many significant, choice edible mushrooms [2]. Bailinggu was first cultivated successfully in China in 1986 [3]. In recent years, more research has been devoted to its genetic diversity, breeding and cultivation [1,2,3,4,5]. Presently, however, only six varieties are used for commercial cultivation. Therefore, strengthening the genetics and breeding of Bailinggu is urgently needed.
Crossbreeding plays an important role in mushroom breeding [6,7], but it is seldom done with Bailinggu. Conventional procedures of mushroom crossbreeding include monokaryon isolation and identification from dikaryotic strains, determination of their mating types, cross-mating two monokaryotic mycelia, microscopic observation of clamp connections at the contact line and hybrid offspring detection [2,6,7,8]. However, these conventional methods are troublesome and time consuming. Simpler, more efficient and precise methods are urgently needed. One such method is molecular-assisted selection, which would benefit the application of mushroom genetics and breeding strategies.
Nowadays, molecular markers are key tools that support and expedite genetics and breeding programs [9,10,11]. Simple sequence repeat (SSR) markers are becoming more effective markers, largely due to their codominance, easy amplification, high reproducibility and multi-allelic variation [12]. Transcriptome sequencing is an efficient, cost-effective way to develop expressed sequence tag-simple sequence repeat (EST-SSR) markers. EST-SSR markers have been useful for the assessment of genetic diversity surveys, construction of genetic maps, cultivar identification and other related studies [5,13,14,15,16]. However, research related to the use of molecular markers for identifying monokaryons and their corresponding genetic diversity, mating types and hybrid offspring in Bailinggu breeding is limited [5,16].
In this study, we developed EST-SSR markers from a newly-sequenced transcriptome of Bailinggu using RNA-seq technology. The specific objectives of this study were to: (1) identify SSR loci in the transcriptome of Bailinggu and analyze their distribution, frequency and size; (2) develop novel EST-SSR markers for Bailinggu that can be used to differentiate monokaryons and their parent dikaryons, identify monokaryon mating types and genetic diversity and detect crossings of F1 and F2 hybrids with monokaryons; and (3) test the transferability of the novel EST-SSR markers in closely-related species Pleurotus eryngii var. ferulae and Pleurotus ostreatus. The newly-identified EST-SSR markers could become a valuable resource for further crossbreeding of Bailinggu and closely-related species.

2. Materials and Methods

2.1. Isolation of Protoplast-Derived Monokaryons and Determination of Mating Types

Five dikaryotic strains of Bailinggu (CCMJ814, 973, 980, 1077 and 1123) were used to isolate protoplast-derived monokaryons (Table 1). The growth medium for the dikaryotic mycelia was a liquid MYG medium (10 g maltose, 5 g yeast extract powder, 5 g glucose and 1000 mL water). Incubation using static cultures was performed at 24 °C for seven to 11 days. Mycelia were then harvested by filtration, washed with 0.6 M D-mannitol (Dingguo, Beijing, China) as an osmotic stabilizer and incubated for 150 min at 27 °C in 2.0% lywallzyme (Guangdong Institute of Microbiology, Guangzhou, China) [9,17]. The obtained protoplasts were filtrated using cotton with 0.6 M D-mannitol and then resuspended in 6 M D-mannitol. The final concentration was approximately 1 × 106 protoplasts/mL. Next, we cultured the diluted protoplasts on regeneration media (liquid maltose yeast glucose (MYG) medium added to 15 g agar and 0.6 M D-mannitol per liter) at 24 °C. The regenerated colonies from the protoplasts were examined under a microscope with fluorescent staining using 4’,6-diamidino-2-phenylindole (DAPI, Sigma-Aldrich, St. Louis, MO, USA) to observe clamp connections and the number of nuclei [18]. The monokaryon should only have one nucleus without clamp connections.
To determine the mating types of monokaryons from the same parent strain, the identified monokaryons were subjected to a monokaryon-monokaryon (mono-mono) mating test at 24 °C on a solid maltose yeast glucose medium for two weeks using somatic incompatibility tests [3,9]. Mating mycelia at the contact zone were selected to observe the presence of clamp connections using a microscope. The presence of clamp connections indicates that the mating types of two monokaryons are different, while the absence of clamp connections indicates that the mating types of two monokaryons are the same.

2.2. Isolation of Basidiospore-Derived Monokaryons

Seven dikaryotic strains of Bailinggu, including CCMJ967, 968, 973 (CCMJ973 was used to obtain protoplast- and basidiospore-derived monokaryons at the same time), 974, 1001, 1002 and 1044, were used to isolate monokaryons from basidiospore suspensions, performed as previously described by Larraya et al. [19]. The obtained monokaryons were also examined under a microscope with fluorescent staining with DAPI to observe clamp connections and the number of nuclei. Then, the identified monokaryons were subjected to a mono-mono mating test [2]. For the across-taxa transfer of the novel EST-SSR markers in the related species, P. eryngii var. ferulae and P. ostreatus, the basidiospore-derived monokaryons were derived from dikaryotic strains of CCMJ970 (P. eryngii var. ferulae) and CCMJ1080 (P. ostreatus) (Table 2).

2.3. Crossbreeding of Monokaryons

In order to obtain Bailinggu hybrid dikaryotic strains, we used the monokaryons obtained from 11 dikaryotic strains for intra-strain crossing. We then examined the presence of clamp connections in the obtained hybrid dikaryotic mycelia under a microscope. Among them, 72 hybridized combinations were successfully cross-mated and the hybrid F1 mycelia obtained. We then selected a hybridized combination of M28 × M33 to obtain F1 fruiting bodies and F2 population. Monokaryon M28 was obtained from the dikaryotic strain CCMJ1123 (D11), and monokaryon M33 was obtained from the dikaryotic strain CCMJ1077 (D2). To test fruiting of the F1 hybrid dikaryotic mycelia from M28 × M33, the mycelia were transferred to a new MYG plate and incubated at 24 °C for two weeks [2]. We next used plastic bottles as containers for the cultivation substrate, which contained 780 g of medium, including 35% sawdust, 5% corncob, 24% wheat bran, 10% maize powder and 4.5% soybean meal. These bottles were incubated at 24 °C for 60 days. Bottles were then processed by scraping off a part of the mycelia surface and moving them to a cropping room at 13–15 °C, with fluorescent lighting of 600–800 lx for 10 h and humidity of 93–95%. We harvested 13 F1 fruiting bodies within 30–50 days in the cropping room and gathered their spores. Single-spore isolation from F1 fruiting bodies and hybridization were performed to acquire the F2 population. All strains were maintained in the Engineering Research Center of the Chinese Ministry of Education for Edible and Medicinal Fungi of the Jilin Agriculture University, China.

2.4. Development of Novel Expressed Sequence Tag-Simple Sequence Repeat Markers Using Transcriptome Sequencing

The dikaryotic mycelia samples of Bailinggu (CCMJ1044) used for transcriptome analysis were kindly provided by the Hengdaxing mushroom factory (Beijing, China). The materials were handled with the same method described by Fu et al. [16]. To perform transcriptome analysis, the mycelia from four treatment stages were used including the control sample that was grown at 25 °C and three cold stage samples that were grown at −3 °C for 2, 6 and 10 d. All samples were immediately frozen in liquid nitrogen and stored at −80 °C. We used a TRIzol reagent (Life technologies, Grand Island, NY, USA) to extract the total RNA. We further evaluated the samples to ensure their eligibility by examining their purity, concentration and integrity with a Nanodrop 2000 UV-Vis spectrophotometer (Thermo Scientific, Waltham, MA, USA), a Qubit 3.0 Fluorometer (Life Technologies, Grand Island, NY, USA) and an Agilent Technologies 2100 Bioanalyzer (Santa Clara, CA, USA), respectively. We constructed 12 cDNA libraries (each stage had three replicates) and performed paired-end sequencing reactions on an Illumina HiSeq 2500 platform at Biomarker Technologies Co., Ltd. (Beijing, China). Sequencing data are available in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA), under Accession Number SRP090421. Raw reads were trimmed with adapters, and low-quality reads were filtered. We then used Trinity (http://trinityrnaseq.sf.net) [16] software to assemble clean, short reads from the 12 libraries, yield transcripts and further cluster into unigenes. Filtration was done after assembly, and the parameter set was relaxed slightly compared with our previous analysis [16] in order to obtain more useful data.
MIcroSAtellite (MISA) (http://pgrc.ipk-gatersleben.de/misa/) [16] was applied to search for SSRs in 31,139 assembled unigenes. The criteria for the default parameters of the SSR loci were as follows: mono-, di-, tri-, tetra-, penta- and hexa-nucleotide motifs with a minimum of ten, six, five, five, five and five repeats, respectively. The maximum number of bases that interrupted two SSRs in a compound locus was 100. Because of unigene redundancy, we filtered SSRs by integrating those that were located in the same gene, but corresponded to different isoforms. Based on the above procedures, we investigated the characteristics of SSRs in the transcriptome of Bailinggu, including density (number of loci per Kbp), motif (type, length and repeat times) and location (coding sequence, non-coding sequence or untranslated regions).

2.5. Validation of Novel Expressed Sequence Tag-Simple Sequence Repeat Markers

Primer 3 software (http://www.premierbiosoft.com) [16] was used to design primer pairs for EST-SSRs, with the following criteria: polymerase chain reaction (PCR) product size was 100–400 bp, primer length was 18–22 nucleotides, GC content was 40–70% and annealing temperature was 60 °C. These markers were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China). We then randomly chose 100 primer pairs to test whether they could successfully amplify with eight Bailinggu monokaryons (four from columnar-shaped strains and four from palmate-shaped strains).
The EST-SSR primer pairs that successfully obtained stabilized bands were used to differentiate monokaryons from their dikaryotic strains, identify monokaryon mating type, and distinguish F1 and F2 hybrids. Next, we investigated the genetic diversity of 37 Bailinggu monokaryons. First, the genomic DNA from fresh mycelia was extracted with the Plant DNA Mini Kit (Kangwei, Beijing, China), following the manufacturer’s instructions. Second, the PCR were carried out in 25-μL final reaction volumes, containing 20 ng template DNA, 1× buffer and 0.2 mM dNTPs, 0.2 μM primer pairs and 0.5 U Taq DNA polymerase (Dingguo, Beijing, China). PCR was performed using a Thermal Cycler (Bio-Rad, Brockport, NY, USA), with the following thermal cycling protocol: 5 min at 94 °C, then the samples were processed through 30 cycles for 30 s at 94 °C, 30 s at 60 °C and 30 s at 72 °C. All amplifications were completed with a 10-min extension step at 72 °C. The PCR products were separated by 8% (w/v) non-denaturing polyacrylamide gel electrophoresis in 1× tris-borate-EDTA (TBE) buffer. The polymorphic bands were visualized by silver staining [20]. Amplification products were scored as present (1) or absent (0). We used GenAIEx6.502 software [21] and NTSYSpc2.1 software (East Setauket, NY, USA) to calculate the genetic diversity parameters and construct an unweighted pair group method with arithmetic mean (UPGMA) dendrogram, respectively [22].

2.6. Across-Taxa Transferability of Novel Expressed Sequence Tag-Simple Sequence Repeat Markers in Monokaryons of Related Species of the Genus Pleurotus

Two closely-related species of the genus Pleurotus (P. eryngii var. ferulae and P. ostreatus) were selected to measure the across-taxa transfer rate of the new EST-SSR markers in their six single-spore-derived monokaryons. All 100 EST-SSR markers mentioned above were used to measure the transfer rate, which was calculated as the proportion of SSR primers that showed successful amplification in each species.

3. Results and Discussion

3.1. Isolation of Protoplast-Derived Monokaryons

We optimized the protocol for obtaining the monokaryons from Bailinggu by using the protoplast-derived method. We found that culturing liquid mycelia for 11 days was more effective than the traditional seven days. The rate of the protoplast-derived monokaryon development from regenerated colonies was more than 93% using 11-day-old liquid mycelia, compared to about 60% with a seven-day-old liquid. Microscopic analysis showed that a total of 162 protoplast-derived monokaryons were obtained from five dikaryotic strains of Bailinggu. Among them, 36, 4, 50, 50 and 22 protoplast-derived monokaryons were obtained from dikaryotic strains CCMJ973, 980, 1077, 1123 and 814, respectively. We used these monokaryons for further EST–SSR marker validation, genetic diversity estimation and Bailinggu hybridization.
Monokaryons obtained from protoplasts were vegetative progeny taken directly from the nutrient hypha without undergoing meiosis. Therefore, we usually could obtain two different mating types (AxBx and AyBy) of protoplast-derived monokaryons. However, we found that there is an imbalance in the recovery of the two mating types from a dikaryon after protoplasting in this study. The mono-mono mating-type tests revealed that the four dikaryotic strains CCMJ973, 980, 1077 and 1123 produced monokaryotic strains that had the two mating types, while the dikaryotic strain CCMJ814 produced monokaryotic strains with only one of the mating types. Previous fungi studies also demonstrated that one or two mating type monokaryotic strains could be obtained from their parent dikaryotic strains through protoplast regeneration technology [9,23]. In addition, the ratio of two different mating types of monokaryons from the four dikaryotic strains CCMJ973, 980, 1077 and 1123 were 3:1, 3:1, 1.63:1 and 2.57:1, respectively. Theoretically, if protoplast-derived monokaryons are regenerated from the same dikaryon, the two mating types will be equally distributed [23]. However, in our study, the ratio of two different mating types ranged from 1.63:1–3:1. Lin and Zhang also demonstrated that the proportion of the two mating types is not always 1:1 in many types of fungi [24]. We infer that our imbalanced ratios might be caused by different protoplast regeneration rates and different mycelia growth rates of the diverse parent dikaryons. Therefore, we urgently need useful molecular markers to improve the identification efficiency of monokaryons and their mating types due to the imbalance that occurred in protoplast-derived monokaryons.

3.2. Novel EST-SSR Markers Developed Using Transcriptome Sequencing

Most research on the genetics and molecular biology of Bailinggu has concentrated on genetic diversity analyses of dikaryotic germplasm resources [1,5,16]. Consequently, effective molecular markers for use in monokaryon identification are limited, which in turn hinders crossbreeding of Bailinggu. Therefore, in order to broaden the range of molecular tools for Bailinggu, we developed a large number of EST-SSR loci using Bailinggu transcriptome data. With the Illumina HiSeq 2500, we employed a paired-end sequencing strategy on samples of Bailinggu mycelia. This technology has a high throughput and is an efficient, rapid and cost-effective tool for SSR mining, particularly for less well-studied species [25,26,27]. In this study, we obtained substantially new sequence data and SSR loci information for Bailinggu. We found that 2211 (5.66%) out of 31,139 unigenes harbored SSR loci and that 165 unigenes have more than one SSR locus. In total, 1110 EST-SSR loci were identified in these unigenes (Table 3). The density of EST-SSR loci in the transcriptome was 0.018 loci per Kbp, which was lower than Grifola frondosa (0.034 per Kbp) [28], Agaricus subrufescens (0.024 per Kbp) [27] and Auricularia polytricha (0.089 per Kbp) [29]. The discrepancy between SSR frequencies may be partially related to differences in species, SSR screening criteria or the amount of detected sequences.
The number of SSR loci identified in this Bailinggu study was approximately twice that of the previous report [16] because we relaxed the filtration parameters slightly after assembly with Trinity. In order not to affect our downstream analysis, the identified SSR loci were then filtered—those located in the same gene, but corresponding to different isoforms were integrated—because some unigenes were redundant. As a result, we increased the possibility for developing effective markers that can be used for further studies.
The trinucleotide repeats (TNR) were the most common type of SSR motif (565, 50.90%), followed by dinucleotide repeats (DNR) (232, 20.90%), hexanucleotide repeats (HNR) (107, 9.64%), mononucleotide repeats (MNR) (94, 8.49%), pentanucleotide repeats (PNR) (44, 3.96%) and tetranucleotide repeats (TTNR) (30, 2.70%). In addition, 38 EST-SSRs were present in compound formations (Figure 1). The finding of TNR as the most abundant repeat motif in the Bailinggu transcriptome is similar to studies of many other fungi, such as A. polytricha [29], Lentinula edodes [30], Boletus edulis [31], Volvariella volvacea [32] and Phellinus linteus [14]. Therefore, a high frequency of tri-nucleotide repeats might be a common feature in EST-derived SSRs. Previous studies suggested that the TNR type of microsatellites should be generated from DNA replication slippage and would be less likely to cause frame-shift mutations in coding regions, resulting in a relatively lower impact on gene expressions [33]. However, an aberrance in other motifs types (except HNR) can lead to a variation of open reading frames, which can cause a frameshift mutation and lead to the formation of shorter or completely different proteins [34]. Meanwhile, the high frequency of DNR may be explained by their abundance in several codons with different nucleotide arrangements. Some DNR, such as (AG)n/(CT)n, are not selectively neutral and may have functional roles [35].
We identified more than 234 motif sequence types (including complementary sequences). The G/C motif appeared 58 times in MNR and was more prevalent than the A/T motif. Five tandem repeats were the most common repeat number (260, 46.02%) in TNR, followed by six tandem repeats in DNR (126, 54.31%) and TNR (107, 18.94%). In addition, the mean length of all 1110 SSRs was 21.76 bp, and the mean lengths of the different repeat types, in increasing order, were DNR (15.34 bp), MNR (16.48 bp), TNR (21.14 bp), TTNR (26.40 bp), PNR (23.98 bp) and HNR (27.70 bp). The largest number of repetitions was 56 for DNR and 39, 26, 21, 16 and 9 for MNR, TNR, TTNR, PNR and HNR, respectively.
We next investigated the distributional characteristics of SSR loci in the unigenes from the Bailinggu transcriptome. Of these 1110 EST-SSRs, 131 and 490 were located in the coding sequence regions (CDS) and untranslated regions (UTR) respectively, and 192 were located over CDS such that each SSR locus was located in both the CDS and UTR regions. Due to the lack of enough information to delimit the CDS region in some unigenes harboring SSR loci, 297 EST-SSRs were not ascertained for their distributional characteristics. The density of SSRs in the UTR was three-times more than the density in the CDS. A similar pattern was also observed in previous transcriptome studies [35,36]. For SSR present in the CDS, TNR (63, 48.10%) was the most abundant followed by DNR (33, 25.19%), MDR (10, 7.63%), HDR (8, 6.11%), PDR (7, 5.343%) and TTDR (6, 4.58%). A high frequency of these SSRs located in transcribed regions, especially in the untranslated portions, could be attributed to mutation and selection pressure of amino acids [33,37]. For breeding, because of the association with important economic and agronomic traits, SSR loci located in the transcribed regions may be of more interest.

3.3. Validation of Novel Expressed Sequence Tag-Simple Sequence Repeats

We successfully designed 317 primer pairs using Primer 3. To identify the efficiency of these primer pairs and to screen for high-quality EST-SSR markers useful for crossbreeding of Bailinggu, we selected a set of 100 primer pairs (Supplementary Materials, Table S1). These primer pairs included 34 DNR, 54 TNR, 1 TTNR, 1 PNR, 3 HNR and 7 multiple repeats. Among these primers, 86 pairs (86.00%) showed a clear amplification in Bailinggu monokaryons, whereas the remaining 14 failed to amplify. These 86 EST-SSR markers will be useful in constructing a genetic linkage map for Bailinggu. We also developed 39 EST-SSRs (blgSSR1–6, 10–42) and 43 EST-SSRs (blgSSR1–38, 43–47) and used them to quickly and efficiently differentiate monokaryons from dikaryons and identify the mating types of monokaryons, respectively (Figure 2A).
Next, due to the codominance of SSR markers, we found 34 EST-SSRs (blgSSR9–19, 24–27, 36–41, 45–57) useful for identifying F1 and F2 hybrid offspring (Figure 2B,C). We also found that all individuals of the F1 generation had the same genotype and similar phenotype, which corresponds with Mendel’s law of uniformity. Moreover, the genetic separation ratio of the F2 generation using EST-SSR markers equaled 1:2:1 (homozygous like parent A: heterozygous: homozygous like parent B), which matches the F2 generation of plants [38]. Therefore, these newly-developed EST-SSR markers are effective and can be used for further genetic and breeding studies of Bailinggu.
Conventional methods for crossbreeding of Bailinggu are troublesome and time-consuming, involving monokaryon isolation, antagonistic activity tests, mono-mono mating tests, microscopic observation and isozyme electrophoresis [39,40,41]. Additionally, to identify hybrids, the traditional detection method of esterase isozyme electrophoresis is quite complicated and is influenced by external conditions [42]. Therefore, we need to explore a new appraisal method that is simpler, more efficient and precise. The development of molecular markers is an efficient method to enhance crossbreeding of Bailinggu. Previous reports have shown that SSR markers are simple, highly stable, more efficient and more precise in the identification and discrimination of breeding materials [10,43]. For our study, using traditional methods to validate the 162 monokaryons and hybrids would have required more than 20–30 days. In contrast, using our newly-developed SSR markers for DNA preparation, PCR amplification and polyacrylamide gel electrophoresis detection, this same validation can be completed within 1–2 days. Compared with the traditional method of microscopic examination and isozyme electrophoresis for identifying Bailinggu’s hybrid offspring, SSR markers are faster and more accurate. SSR markers may also be applied to track related genotypes within the offspring. Above all, our newly-developed SSR markers will benefit and facilitate breeding research of Bailinggu.

3.4. Genetic Diversity Analysis of Monokaryons and Dikaryons Using 18 EST-SSRs

To analyze the genetic diversity of Bailinggu monokaryons, we selected 28 single spores derived monokaryons from seven dikaryotic strains and nine protoplast-derived monokaryons from five dikaryotic strains (Supplementary Materials Table S2) using 18 polymorphic, clear and stable EST-SSRs. The genetic parameters of these SSR markers are shown in Table 4. In total, 86 alleles were detected across all 37 monokaryons. The percentage of polymorphic loci was 100%. The observed number of alleles per locus (Na) ranged from three (blgSSR30) to seven (blgSSR24), with a mean value of 4.778 per locus. Shannon’s information index (I) value ranged from 0.616–1.587, with a mean of 1.138. The polymorphism information content (PIC) of the 18 EST-SSRs ranged from 0.321–0.766, with a mean of 0.584, which was higher than previously reported means for Cordyceps militaris (0.191) [44], Flammulina velutipes (0.420) [45] and Auricularia auricula-judae (0.470) [12], but lower than Agaricus bisporus (0.618) [46]. The PIC value was greater than 0.500, depicting a high discrimination power. Among these 18 EST-SSRs, 14 (blgSSR11, 13, 14, 16, 24, 26, 27, 28, 29, 32, 37, 38, 43 and 56) were high polymorphic loci with PIC values of more than 0.500; none were low polymorphic loci with PIC values less than 0.250. The higher PIC values for our EST-SSRs indicate that these polymorphic EST-SSR markers, derived from the mycelia transcriptome, have suitable power for discriminating among the monokaryons of Bailinggu. Thus, our EST-SSR markers will help broaden the gene pool and improve Bailinggu crossbreeding programs aimed at reducing costs and conserving energy. The maximum and the minimum genetic distance between the monokaryons was 1.022 and 0.000, respectively, with a mean of 0.641. These results indicated that the 37 monokaryons of Bailinggu had a relatively high genetic diversity. In further breeding programs, these polymorphic and stable EST-SSR markers can be used for construction of a genetic linkage map and quantitative trait locus (QTL) mapping for Bailinggu.
According to the cluster analysis, based on the Dice similarity coefficient and UPGMA algorithm, the genetic similarity coefficient among all 37 monokaryons ranged from 0.19–1.00 (Figure 3A). Monokaryons derived from one dikaryotic strain were clustered into the same group, except monokaryons (M1–M9) derived from the dikaryotic strains CCMJ974 (D4) and CCMJ968 (D5). We discovered that these two wild dikaryotic strains were collected from Yumin, China, in 2012 and inferred that they may have the same genetic background. Combined analysis of monokaryons and dikaryons allows the possibility of more information about the markers, so we also analyzed the genetic relationship of the dikaryotic strains using the same 18 EST-SSR markers. We found that the genetic similarity coefficient between the 11 dikaryotic strains ranged from 0.29–0.95 (Figure 3B). This result was similar to the range of the clustered monokaryons. However, these dikaryotic strains have more heterozygosity associated with each locus.
The above results showed that these Bailinggu strains from natural populations have a relatively high genetic diversity. We hypothesized that the natural genetic recombination may also contribute to the relatively high genetic diversity of Bailinggu populations. Therefore, we used these 18 EST-SSR markers to assess the possibility of natural genetic recombination occurring in these 11 dikaryotic strains of Bailinggu. Based on the PCR amplification products and non-denaturing PAGE results (Figure S1), four out of 18 EST-SSR markers (blgSSR 14, 30, 31, 56) could identified as having the natural genetic recombination in some of the 11 dikaryotic strains of Bailinggu, although they were from natural populations. These results indicated that not only does natural gene mutation contribute to the relatively high genetic diversity of Bailinggu populations, but also natural genetic recombination.

3.5. Across-Taxa Transferability of Novel Expressed Sequence Tag-Simple Sequence Repeat Markers in the Monokaryons of Two Species of the Genus Pleurotus

To evaluate the possibility of the EST-SSR markers from the Bailinggu transcriptome used in monokaryons from closely-related species of the genus Pleurotus, we carried out across-taxa transferability analysis for these 100 novel EST-SSR markers in the monokaryons from P. eryngii var. ferulae and P. ostreatus. In the monokaryons of P. eryngii var. ferulae, 72 primers (blgSSR 2–14, 18, 20–23, 26–28, 31, 33, 34, 36, 38, 40–44, 46, 48–49, 52–54, 56–58, 60–82, 88–90, 93–98) were amplified products. In monokaryons of P. ostreatus, 64 primers (blgSSR 4–13, 15, 20–22, 25–27, 34, 36–37, 40–46, 48–54, 57, 59–82, 88, 95–98) were amplified products. Therefore, the overall transferability rates of these 100 novel EST-SSRs were 72% in P. eryngii var. ferulae and 64% in P. ostreatus.
Previous studies showed the significant or moderate conservation of EST-SSR markers in the related species within the same genus in fungi and plants [47]. We also found relatively high transferability rates (82%) of EST-SSR markers developed from Bailinggu in the dikaryons of P. eryngii var. ferulae [5]. In addition, the across-taxa transferability rate in monokaryons of P. eryngii var. ferulae (72%) was higher than in P. ostreatus (64%). A likely reason for this difference is that P. eryngii var. ferulae has a closer genetic relationship with Bailinggu than P. ostreatus, resulting in moderate cross-species conservation of SSR-containing genes [1,5]. In fact, we found some unigenes known to encode conserved proteins in fungi contained in our SSR loci, such as unigenes that code hydrophobin (blgSSR31), ATP-dependent helicase (blgSSR10), PH_Scd1 and vWFA domain-containing protein (blgSSR15), the bZIP domain of ATF-2 and similar proteins (blgSSR52), the catalytic domain of cell division control protein 7-like serine/threonine kinases (blgSSR43) and some hypothetical proteins (blgSSR53, 63, 84). Therefore, the transferability analysis indicates that our EST-SSR markers are also potentially useful for correlational research in closely-related species of the genus Pleurotus, especially species without sequenced genomes or molecular markers such as P. eryngii var. ferulae. Furthermore, compared to the genomic SSRs, EST-SSRs have a higher transfer rate in related species [47,48]. The likely reason for this difference is that most of the EST-derived SSRs are originated from transcribed regions, and therefore, they are more conserved across a number of related species within the same genus. Therefore, EST-SSR markers are not only a valuable tool for mushroom breeding, but also a better choice for cross-species phylogenetic studies.

4. Conclusions

We developed a new set of EST-SSR markers, using Bailinggu transcriptome data, which can be used for Bailinggu crossbreeding. Our research contributes the following resources and information: (1) 39, 43 and 34 novel EST-SSR markers that successfully identified monokaryons, differentiated two different monokaryon mating types and discriminated F1 and F2 hybrid offspring, respectively; (2) monokaryons of Bailinggu have a relatively high genetic diversity; (3) novel EST-SSR markers with relatively high transferability rates in monokaryons of P. eryngii var. ferulae and P. ostreatus; and (4) newly-developed EST-SSR markers that can be used to quickly, easily and precisely differentiate monokaryons from parent dikaryons and identify monokaryon mating types and hybrid dikaryotic strains. Consequently, the newly-developed EST-SSR markers could replace conventional methods for crossbreeding. Ultimately, our study has the potential to significantly enhance breeding research in Bailinggu and other closely-related species.

Supplementary Materials

The following are available online at www.mdpi.com/2073-4425/8/11/325/s1. Table S1: Characteristics of 100 simple sequence repeat (SSR) markers of Bailinggu used in this study. Table S2: Information for the 37 monokaryons used in the genetic diversity analysis. Figure S1: The natural genetic recombination found in the dikaryotic strains of Bailinggu, based on the PCR amplification products using four EST-SSR primers and non-denaturing PAGE.

Acknowledgments

This study was financially supported by the National Basic Research Program of China (No. 2014CB138305), the National Natural Science Foundation of China (31471926) and the University S & T Innovation Platform of Jilin Province for Economic Fungi (#2014B-1)). We are grateful to Yuning He and Dingxuan Huang for kindly providing the sequenced samples of Bailinggu and to Peng Wei for kindly providing the wild strains of Bailinggu. We also thank Linda R. Klein for valuable writing advice and editing the manuscript.

Author Contributions

Yu Li and Yongping Fu conceptualized the study and designed the research framework. Yueting Dai and Yongping Fu wrote the manuscript and performed the data analyses. Yueting Dai worked on the SSR identification experiments. Wenying Su screened the SSR markers. Chentao Yang performed assembly and annotation of the transcriptome. Bing Song designed the SSR primers. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare that they have no competing financial interests.

References

  1. Zhao, M.R.; Zhang, J.X.; Chen, Q.; Wu, X.L.; Gao, W.; Deng, W.Q.; Huang, C.Y. The famous cultivated mushroom Bailinggu is a separate species of the Pleurotus eryngii species complex. Sci. Rep. 2016, 6, 33066. [Google Scholar] [CrossRef] [PubMed]
  2. Kawai, G.; Babasaki, K.; Neda, H. Taxonomic position of a Chinese Pleurotus “Bai-Ling-Gu”, it belongs to Pleurotus eryngii (DC., Fr.) Quél. and evolved independently in China. Mycoscience 2008, 49, 75–87. [Google Scholar] [CrossRef]
  3. Mou, C.J.; Cao, Y.Q.; Ma, J.L. A new variety of Pleurotus eryngii and its cultural characters. Acta Mycol. Sin. 1987, 6, 153–156. [Google Scholar]
  4. Zervakis, G.I.; Ntougias, S.; Gargano, M.L.; Besi, M.I.; Polemis, E.; Typas, M.A.; Venturella, G. A reappraisal of the Pleurotus eryngii complex—New species and taxonomic combinations based on the application of a polyphasic approach, and an identification key to Pleurotus taxa associated with Apiaceae plants. Fungal Biol. 2014, 118, 814–834. [Google Scholar] [CrossRef] [PubMed]
  5. Fu, Y.P.; Yueting Dai, Y.T.; Yang, C.T.; Wei, P.; Song, B.; Yang, Y.; Sun, L.; Zhang, Z.W.; Li, Y. Comparative transcriptome analysis identified candidate genes related to Bailinggu mushroom formation and genetic markers for genetic analyses and breeding. Sci. Rep. 2017, 7, 9266. [Google Scholar] [CrossRef] [PubMed]
  6. Kim, M.K.; Ryu, J.-S.; Lee, Y.-H.; Kim, H.-R. Breeding of a long shelf-life strain for commercial cultivation by mono–mono crossing in Pleurotus eryngii. Sci. Hortic. 2013, 162, 265–270. [Google Scholar] [CrossRef]
  7. Liu, S.R.; Ke, B.B.; Zhang, W.R.; Liu, X.H.; Wu, X.P. Breeding of new Ganoderma lucidum strains simultaneously rich in polysaccharides and triterpenes by mating basidiospore-derived monokaryons of two commercial cultivars. Sci. Hortic. 2017, 216, 58–65. [Google Scholar] [CrossRef]
  8. Kim, K.-H.; Kang, Y.M.; Im, C.H.; Ali, A.; Kim, S.Y.; Je, H.-J.; Kim, M.-K.; Rho, H.S.; Lee, H.S.; Kong, W.-S.; et al. Identification and functional analysis of pheromone and receptor genes in the B3 mating locus of Pleurotus eryngii. PLoS ONE 2014, 9, e104693. [Google Scholar]
  9. Li, D.; Liu, Y.; Wang, P.; Ma, Y.; Wang, S.; Zhao, S.; Xu, F. Development of SCAR markers to determine the mating types of Lepista nuda protoplast monokaryons. Curr. Microbiol. 2014, 68, 536–542. [Google Scholar] [CrossRef] [PubMed]
  10. Liu, X.B.; Feng, B.; Li, J.; Yan, C.; Yang, Z.L. Genetic diversity and breeding history of winter mushroom (Flammulina velutipes) in China uncovered by genomic SSR markers. Gene 2016, 591, 227–235. [Google Scholar] [CrossRef] [PubMed]
  11. Kim, H.J.; Jung, J.; Kim, M.-S.; Lee, J.M.; Choi, D.; Yeam, I. Molecular marker development and genetic diversity exploration by RNA-seq in Platycodon grandiflorum. Genome 2015, 58, 1–11. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, R.Y.; Hu, D.D.; Gu, J.G.; Hu, Q.X.; Zuo, X.M.; Wang, H.X. Development of SSR markers for typing cultivars in the mushroom Auricularia auricula-judae. Mycol. Prog. 2012, 11, 587–592. [Google Scholar] [CrossRef]
  13. Yan, H.D.; Zhang, Y.; Zeng, B.; Yin, G.H.; Zhang, X.Q.; Ji, Y.; Huang, L.K.; Jiang, X.M.; Liu, X.C.; Peng, Y.; et al. Genetic diversity and association of EST-SSR and SCoT markers with rust traits in orchardgrass (Dactylis glomerata L.). Molecules 2016, 11, 66. [Google Scholar] [CrossRef] [PubMed]
  14. Huang, Y.; Wu, X.; Jian, D.; Zhan, Y.; Fan, G. De novo transcriptome analysis of a medicinal fungi Phellinus linteus and identification of SSR markers. Biotechnol. Biotechnol. Equip. 2015, 29, 395–403. [Google Scholar] [CrossRef] [PubMed]
  15. Kim, T.S.; Raveendar, S.; Suresh, S.; Lee, G.A.; Lee, J.R.; Cho, J.H.; Lee, S.Y.; Ma, K.H.; Cho, G.T.; Chung, J.W. Transcriptome analysis of two Vicia sativa subspecies: Mining molecular markers to enhance genomic resources for vetch improvement. Genes 2015, 6, 1164–1182. [Google Scholar] [CrossRef] [PubMed]
  16. Fu, Y.P.; Liang, Y.; Dai, Y.T.; Yang, C.T.; Duan, M.Z.; Zhang, Z.; Hu, S.N.; Zhang, Z.W.; Li, Y. De novo sequencing and transcriptome analysis of Pleurotus eryngii subsp. tuoliensis (Bailinggu) mycelia in response to cold stimulation. Molecules 2016, 21, 560. [Google Scholar]
  17. Wang, M.Z.; Liu, S.S.; Li, Y.Y.; Xu, R.; Lu, C.H.; Shen, Y.M. Protoplast mutation and genome shuffling induce the endophytic fungus Tubercularia sp. TF5 to produce new compounds. Curr. Microbiol. 2010, 61, 254–260. [Google Scholar] [CrossRef] [PubMed]
  18. Oliveira, L.N.; Casaletti, L.; Bao, S.N.; Borges, C.L.; Lima, P.D.S.; Soares, C.M.D.A. Characterizing the nuclear proteome of Paracoccidioides spp. Fungal Biol. 2016, 120, 1209e1224. [Google Scholar] [CrossRef] [PubMed]
  19. Larraya, L.; Peñas, M.M.; Pérez, G.; Santos, C.; Ritter, E.; Pisabarro, A.G.; Ramírez, L. Identification of incompatibility alleles and characterisation of molecular markers genetically linked to the A incompatibility locus in the white rot fungus Pleurotus ostreatus. Curr. Genet. 1999, 34, 486–493. [Google Scholar] [CrossRef] [PubMed]
  20. Chen, X.M.; Li, J.K.; Xiao, S.J.; Liu, X.D. De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile. Gene 2016, 576, 537–543. [Google Scholar] [CrossRef] [PubMed]
  21. Peakall, R.; Smouse, P.E. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [PubMed]
  22. Parekha, M.J.; Kumar, S.; Zalaa, H.N.; Fougat, R.S.; Patel, C.B.; Bosamia, T.C.; Kulkarni, K.S.; Parihar, A. Development and validation of novel fiber relevant dbEST-SSR markers and their utility in revealing genetic diversity in diploid cotton (Gossypium herbaceum and G. arboreum). Ind. Crop. Prod. 2016, 83, 620–629. [Google Scholar] [CrossRef]
  23. Zhao, J.; Zhang, S.T. Monokaryotization by protoplasting heterothallic species of edible mushrooms. World J. Microbiol. Biotechnol. 1993, 9, 538–543. [Google Scholar] [CrossRef] [PubMed]
  24. Lin, F.C.; Zhang, S.T. Analysis of incompatibility factors in cultivated strains of Lentinus edodes in China. J. Huazhong Agric. Univ. 1995, 14, 459–466. [Google Scholar]
  25. Seong, J.; Kang, S.W.; Patnaik, B.B.; Park, S.Y.; Hwang, H.J.; Chung, J.M.; Song, D.K.; Noh, M.Y.; Park, S.-H.; Jeon, G.J.; et al. Transcriptome analysis of the tadpole shrimp (Triops longicaudatus) by Illumina paired-end sequencing: Assembly, annotation, and marker discovery. Genes 2016, 7, 114. [Google Scholar] [CrossRef] [PubMed]
  26. Vatanparast, M.; Shetty, P.; Chopra, R.; Doyle, J.J.; Sathyanarayana, N.; Egan, A.N. Transcriptome sequencing and marker development in winged bean (Psophocarpus tetragonolobus; Leguminosae). Sci. Rep. 2016, 6, 29070. [Google Scholar] [CrossRef] [PubMed]
  27. Foulongne-Oriol, M.; Lapalu, N.; Ferandon, C.; Spataro, C.; Ferrer, N.; Amselem, J.; Savoie, J.M. The first set of expressed sequence tags (EST) from the medicinal mushroom Agaricus subrufescens delivers resource for gene discovery and marker development. Appl. Microbiol. Biol. 2014, 98, 7879–7892. [Google Scholar] [CrossRef] [PubMed]
  28. Yuan, W.D.; Lu, N.; Chen, Q.; Song, J.L.; Wang, W.K. Analysis on transcriptome sequenced of maitake. J. Fudan Univ. 2015, 54, 673–678. [Google Scholar]
  29. Zhou, Y.; Chen, L.F.; Fan, X.Z.; Bian, Y.B. De novo assembly of Auricularia polytricha transcriptome using Illumina sequencing for gene discovery and SSR marker identification. PLoS ONE 2014, 9, e91740. [Google Scholar] [CrossRef] [PubMed]
  30. Lin, F.X.; Cheng, S.M.; Li, A.Z.; Lin, F.C. Compose and distribution of EST-SSR in genomes of Lentinula edodes. J. Microbiol. Chin. 2007, 34, 438–442. [Google Scholar]
  31. Wang, Y.; Chen, M.J.; Wang, H.; Bao, D.P. Distribution pattern analysis of SSR loci in the genome of Boletus edulis. Mycosystema 2015, 34, 204–214. [Google Scholar]
  32. Wang, Y.; Chen, M.; Wang, H.; Wang, J.F.; Bao, D. Microsatellites in the genome of the edible mushroom, Volvariella volvacea. Biol. Med. Res. Int. 2014, 2014, 281912. [Google Scholar] [CrossRef] [PubMed]
  33. Huang, X.; Yan, H.D.; Zhang, X.Q.; Zhang, J.; Frazier, T.P.; Huang, D.J.; Lu, L.; Huang, L.K.; Liu, W.; Peng, Y.; et al. De novo transcriptome analysis and molecular marker development of two Hemarthria species. Front. Plant Sci. 2016, 7, 496. [Google Scholar] [CrossRef] [PubMed]
  34. Metzgar, D.; Bytof, J.; Wills, C. Selection against frameshift mutations limits microsatellite expansion in coding DNA. Genome Res. 2000, 10, 72. [Google Scholar] [PubMed]
  35. Hamarsheh, O.; Amro, A. Characterization of simple sequence repeats (SSRs) from Phlebotomus papatasi (Diptera: Psychodidae) expressed sequence tags (ESTs). Parasite Vector 2011, 4, 189. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, F.; Hu, Z.; Liu, W.; Li, J.; Wang, W.; Liang, Z.; Wang, F.; Sun, X. Distribution, function and evolution characterization of microsatellite in Sargassum thunbergii (Fucales, Phaeophyta) transcriptome and their application in marker development. Sci. Rep. 2016, 6, 18947. [Google Scholar] [CrossRef] [PubMed]
  37. Singh, R.; Pandey, B.; Danishuddin, M.; Sheoran, S.; Sharma, P.; Chatrath, R. Mining and survey of simple sequence repeats in wheat rust Puccinia sp. Bioinformation 2011, 7, 291–295. [Google Scholar] [CrossRef] [PubMed]
  38. You, C.Y.; Chen, Q.J.; Qu, Y.U.; Ma, L. Construction and evaluation on crossing population from Gossypium hirsutum and G. barbadense. J. Xinjiang Agric. Univ. 2007, 30, 68–71. [Google Scholar]
  39. Zhao, Y.; Lin, F.; Yan, S.Y.; Wang, H.; Chen, M.J. Establishment and application of an efficient cross breeding method assisted by molecular markers of Volvariella volvacea. Microbiol. Chin. 2015, 42, 1165–1174. [Google Scholar]
  40. Xie, B.G.; Yi, H.; Huang, Z.L.; Lin, Y.Z.; Xie, F.Q.; Iiang, Y.J. Cross-breeding and isoenzyme analysis of Volvariella volvacea. J. Fujian Agric. Univ. 2001, 30, 372–376. [Google Scholar]
  41. Bian, Y.B.; Wu, K.Y.; Wu, C.S. Genetic analysis of isozyme loci of intraspecific hybrid in Auricularia auricula. Acta Genet. Sin. 2003, 30, 76–80. [Google Scholar] [PubMed]
  42. Bian, Y.B.; Luo, X.C.; Zhou, Q. Esterase isozyme zymogram polymorphisms of cultivated strains in Auricularia auricula. Mycosystema 2000, 19, 87–90. [Google Scholar]
  43. Wu, P.; Zhang, L.J.; Zhang, D.; Shang, X.D.; Tan, Q.; Song, C.Y. Identification of Xianggu (Lentinula edodes) monokaryons and hybrid progenies using SSR markers. Microbiol. Chin. 2016, 43, 444–455. [Google Scholar]
  44. Yuan, F.; Yu, H.; Zuo, S.M.; Adams, A. Development of EST-SSR for preliminary analysis of genetic diversity of Cordyceps militaris. Biochem. Syst. Ecol. 2015, 58, 126–131. [Google Scholar] [CrossRef]
  45. Zhang, R.; Hu, D.; Zhang, J.; Zuo, X.; Jiang, R.; Wang, H.; Ng, T.B. Development and characterization of simple sequence repeat (SSR) markers for the mushroom Flammulina velutipes. J. Biosci. Bioeng. 2010, 110, 273–275. [Google Scholar] [CrossRef] [PubMed]
  46. Fu, Y.P.; Wang, X.X.; Li, D.; Liu, Y.; Song, B.; Zhang, C.L.; Wang, Q.; Chen, M.Y.; Zhang, Z.W.; Li, Y. Identification of resistance to wet bubble disease and genetic diversity in wild and cultivated strains of Agaricus bisporus. Int. J. Mol. Sci. 2016, 17, 1568. [Google Scholar] [CrossRef] [PubMed]
  47. Li, X.Y.; Shangguan, L.F.; Song, C.N.; Wang, C.; Gao, Z.H.; Yu, H.P.; Fang, J.G. Analysis of expressed sequence tags from Prunus mume flower and fruit and development of simple sequence repeat markers. BMC Genet. 2010, 11, 66. [Google Scholar] [CrossRef] [PubMed]
  48. Castillo, A.; Budak, H.; Varshney, R.V.; Dorado, G.; Graner, A.; Hernandez, P. Transferability and polymorphism of barley EST-SSR markers used for phylogenetic analysis in Hordeum chilense. BMC Plant Biol. 2008, 8, 97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 08 00325 g001 550
Figure 1. The distribution of SSR motif types in Bailinggu transcriptome.
Figure 1. The distribution of SSR motif types in Bailinggu transcriptome.
Genes 08 00325 g001
Genes 08 00325 g002 550
Figure 2. Examples of expressed sequence tag-simple sequence repeat markers applied in Bailinggu. (A) Polymorphic bands between monokaryons and their parent dikaryons using blgSSR24, blgSSR13, blgSSR25 and blgSSR15 (1, 6, 11, 16: parent dikaryons CCMJ973 (D10); 2–3, 7–8, 12–13, 17–18: monokaryon M29; 4–5, 9–10, 14–15, 19–20: monokaryon M30; the dikaryons CCMJ973 (D10) had two bands, while monokaryons M29 and M30 only had one band at the corresponding position); (B) polymorphic bands between the monokaryotic strains and their F1 hybrid using blgSSR12, blgSSR22 and blgSSR21 (1–2, 7–8 and 13–14: monokaryon M28; 3–4, 9–10 and 15–16: monokaryon M33; 5–6, 11–12 and 17–18: F1 hybrid; monokaryon M28 and monokaryon M33 only had one band, while their F1 hybrid was the combination of M28 and M33); (C) the monokaryotic strains and their F2 hybrid using blgSSR13 (1: monokaryon M28; 2: monokaryon M33; 3–5: F2 homozygous-like monokaryon M28; 6–8: F2 homozygous-like monokaryon M33; 9–14: F2 heterozygous).
Figure 2. Examples of expressed sequence tag-simple sequence repeat markers applied in Bailinggu. (A) Polymorphic bands between monokaryons and their parent dikaryons using blgSSR24, blgSSR13, blgSSR25 and blgSSR15 (1, 6, 11, 16: parent dikaryons CCMJ973 (D10); 2–3, 7–8, 12–13, 17–18: monokaryon M29; 4–5, 9–10, 14–15, 19–20: monokaryon M30; the dikaryons CCMJ973 (D10) had two bands, while monokaryons M29 and M30 only had one band at the corresponding position); (B) polymorphic bands between the monokaryotic strains and their F1 hybrid using blgSSR12, blgSSR22 and blgSSR21 (1–2, 7–8 and 13–14: monokaryon M28; 3–4, 9–10 and 15–16: monokaryon M33; 5–6, 11–12 and 17–18: F1 hybrid; monokaryon M28 and monokaryon M33 only had one band, while their F1 hybrid was the combination of M28 and M33); (C) the monokaryotic strains and their F2 hybrid using blgSSR13 (1: monokaryon M28; 2: monokaryon M33; 3–5: F2 homozygous-like monokaryon M28; 6–8: F2 homozygous-like monokaryon M33; 9–14: F2 heterozygous).
Genes 08 00325 g002
Genes 08 00325 g003 550
Figure 3. Unweighted pair group method with arithmetic mean (UPGMA) dendrograms of Bailinggu based on 18 novel EST-SSR markers. (A) UPGMA dendrogram of 37 monokaryons of Bailinggu based on 18 novel EST-SSR markers. M1–M37 represented monokaryons obtained from parent dikaryons. D1–11 represented the parent dikaryotic strains: D1-CCMJ814, D2-CCMJ1077, D3-CCMJ967, D4-CCMJ974, D5-CCMJ968, D6-CCMJ1044, D7-CCMJ980, D8-CCMJ1001, D9-CCMJ1002, D10-CCMJ973, D11-CCMJ1123. (B) UPGMA dendrogram of 11 parent dikaryotic strains of Bailinggu based on the same 18 novel EST-SSR markers.
Figure 3. Unweighted pair group method with arithmetic mean (UPGMA) dendrograms of Bailinggu based on 18 novel EST-SSR markers. (A) UPGMA dendrogram of 37 monokaryons of Bailinggu based on 18 novel EST-SSR markers. M1–M37 represented monokaryons obtained from parent dikaryons. D1–11 represented the parent dikaryotic strains: D1-CCMJ814, D2-CCMJ1077, D3-CCMJ967, D4-CCMJ974, D5-CCMJ968, D6-CCMJ1044, D7-CCMJ980, D8-CCMJ1001, D9-CCMJ1002, D10-CCMJ973, D11-CCMJ1123. (B) UPGMA dendrogram of 11 parent dikaryotic strains of Bailinggu based on the same 18 novel EST-SSR markers.
Genes 08 00325 g003
Table
Table 1. Information about the five strains of Pleurotus tuoliensis used for isolating protoplast-derived monokaryons in this study.
Table 1. Information about the five strains of Pleurotus tuoliensis used for isolating protoplast-derived monokaryons in this study.
No.Strain NumberStrain NameStrain TypeOrigin
1CCMJ814P. tuoliensisWild strainXinjiang, China
2CCMJ973P. tuoliensisWild strainXinjiang, China
3CCMJ980P. tuoliensisWild strainXinjiang, China
4CCMJ1077P. tuoliensisCultivated strainXinjiang, China
5CCMJ1123P. tuoliensisCultivated strainShandong, China
Table
Table 2. Information about the nine Pleurotus strains used for isolating basidiospore-derived monokaryons in this study.
Table 2. Information about the nine Pleurotus strains used for isolating basidiospore-derived monokaryons in this study.
No.Strain NumberStrain NameStrain TypeOrigin
1CCMJ967P. tuoliensisWild strainXinjiang, China
2CCMJ968P. tuoliensisWild strainXinjiang, China
3CCMJ973P. tuoliensisWild strainXinjiang, China
4CCMJ974P. tuoliensisWild strainXinjiang, China
5CCMJ1001P. tuoliensisWild strainXinjiang, China
6CCMJ1002P. tuoliensisWild strainXinjiang, China
7CCMJ1044P. tuoliensisCultivated strainBeijing, China
8CCMJ970P. eryngii var. ferulaeWild strainXinjiang, China
9CCMJ1080P. ostreatusCultivated strainChangchun, China
Table
Table 3. General summary of Bailinggu transcriptome sequencing.
Table 3. General summary of Bailinggu transcriptome sequencing.
ItemsNumber
Total raw data199.81 million
Total clean data177.98 million
Total number of assembled unigenes31,139
Total size of examined sequences (bp)60,265,141
Total number of identified SSRs1110
Number of sequences containing SSRs2211
Number of sequences containing more than one SSR165
SSR: simple sequence repeat.
Table
Table 4. Characteristics of 18 polymorphic SSR markers used for the phylogenetic analysis of the 37 monokaryons of Bailinggu.
Table 4. Characteristics of 18 polymorphic SSR markers used for the phylogenetic analysis of the 37 monokaryons of Bailinggu.
LocusMotifPrimer Sequence (5′-3′)NaNePICI
blgSSR11(CCAGGA)11F-CTTGAATATTGACGGGAGCC
R-AGGCTGTGGATAAGAACCCC		42.6030.5551.011
blgSSR13(CCA)5(CAC)6(CCT)6F-TCCATAATTGTCATCTCCCCA
R-AAATAATTACGACGGTGGGC		53.8020.7191.457
blgSSR14(AT)7F-GTCAAACTGCCCAAATTCGT
R-CACGGTTGCTCGATTATTCTG		62.9010.6121.034
blgSSR16(AC)7F-ACTGCGTCGTGGTGTACAAG
R-GTCGCTAAGTATCGGTGGGA		42.0420.6070.616
blgSSR20(AG)7F-CGGTGCCAGTGTTCCTTATT
R-TGATTCCTGCCCTTTCATTC		52.3440.4130.708
blgSSR24(GA)7F-ATCGGGAATGGCAATCAATA
R-CTCGAGTCCCGAAACCAATA		74.8720.7661.586
blgSSR26(TCC)6F-CCCCGTCCTCTAGTTCATCA
R-AATACGGGTCGTCAGATTCG		53.0440.6821.166
blgSSR27(GGA)8F-ATCGAAAAATGATTACGCCG
R-GGCGAATTTCCTCTTTAGGG		42.2010.6020.765
blgSSR28(GCG)6F-ACTTTCCAACACCAACTGCC
R-GACCACGAGAGTATCGGAGC		43.2760.6501.208
blgSSR29(GAG)7F-TGTTGCAAGCAATGGAATGT
R-GTGAATTCCAGCGGTTGTTT		53.2930.6461.247
blgSSR30(ACC)6F-GAAGAGTCGCAGACCTCCAC
R-GATGCTTTGGTGGACTTGGT		31.9950.4310.83
blgSSR31(CGC)5F-AACATGTTCTCCAAGGCCAC
R-ATGCGTGCTAACTGAACGTG		42.2220.4940.844
blgSSR32(ATC)6F-GCGGATAACCTACTCGTGGA
R-ACGCCGAAAATTTTGATGTC		53.0060.6121.184
blgSSR37(CT)6F-GAGCCCAAGTGACGTTCCTA
R-CTAGGTGGGACTCCGAAACC		52.3770.6030.993
blgSSR38(AGC)6F-ATAACAAGGCATGTTTCCGC
R-ACAGTCAGGCTCTGGGAGAA		53.5510.6901.309
blgSSR43(TCC)5F-CTCAGCCCCAAATTGAACAT
R-CTAGTGGCGGGAAGACTGAG		42.3330.5841.035
blgSSR56(CCA)6F-GTGTGGGCAGTCGTAGCATA
R-CCCTTCCGTTGCTTTCATTA		6209870.5251.024
blgSSR57(TCC)6F-CTCGACTCCACGAAAGAAGG
R-AACACCCGAAATACGAATGC		41.7340.3210.635
Mean4.7782.8050.5841.036
Na: observed number of alleles per locus; Ne: effective number of alleles per locus; PIC: polymorphic information content; I: Shannon’s information index.

Share and Cite

MDPI and ACS Style

Dai, Y.; Su, W.; Yang, C.; Song, B.; Li, Y.; Fu, Y. Development of Novel Polymorphic EST-SSR Markers in Bailinggu (Pleurotus tuoliensis) for Crossbreeding. Genes 2017, 8, 325. https://doi.org/10.3390/genes8110325

AMA Style

Dai Y, Su W, Yang C, Song B, Li Y, Fu Y. Development of Novel Polymorphic EST-SSR Markers in Bailinggu (Pleurotus tuoliensis) for Crossbreeding. Genes. 2017; 8(11):325. https://doi.org/10.3390/genes8110325

Chicago/Turabian Style

Dai, Yueting, Wenying Su, Chentao Yang, Bing Song, Yu Li, and Yongping Fu. 2017. "Development of Novel Polymorphic EST-SSR Markers in Bailinggu (Pleurotus tuoliensis) for Crossbreeding" Genes 8, no. 11: 325. https://doi.org/10.3390/genes8110325

APA Style

Dai, Y., Su, W., Yang, C., Song, B., Li, Y., & Fu, Y. (2017). Development of Novel Polymorphic EST-SSR Markers in Bailinggu (Pleurotus tuoliensis) for Crossbreeding. Genes, 8(11), 325. https://doi.org/10.3390/genes8110325

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop