Recombination and genetic differentiation among natural populations of the ectomycorrhizal mushroom Tricholoma matsutake from southwestern China

Exploring Mitochondrial Heterogeneity and Evolutionary Dynamics in Thelephora ganbajun through Population Genomics

Limited exploration in fungal mitochondrial genetics has uncovered diverse inheritance modes. The mitochondrial genomes are inherited uniparentally in the majority of sexual eukaryotes, our discovery of persistent mitochondrial heterogeneity within the natural population of the basidiomycete fungus Thelephora ganbajun represents a significant advance in understanding mitochondrial inheritance and evolution in eukaryotes. Here, we present a comprehensive analysis by sequencing and assembling the complete mitogenomes of 40 samples exhibiting diverse cox1 heterogeneity patterns from various geographical origins. Additionally, we identified heterogeneous variants in the nad5 gene, which, similar to cox1, displayed variability across multiple copies. Notably, our study reveals a distinct prevalence of introns and homing endonucleases in these heterogeneous genes. Furthermore, we detected potential instances of horizontal gene transfer involving homing endonucleases. Population genomic analyses underscore regional variations in mitochondrial genome composition among natural samples exhibiting heterogeneity. Thus, polymorphisms in heterogeneous genes, introns, and homing endonucleases significantly influence mitochondrial structure, structural variation, and evolutionary dynamics in this species. This study contributes valuable insights into mitochondrial genome architecture, population dynamics, and the evolutionary implications of mitochondrial heterogeneity in sexual eukaryotes.

Genetic diversity and population differentiation in Earliella scabrosa, a pantropical species of Polyporales

Earliella scabrosa is a pantropical species of Polyporales (Basidiomycota) and well-studied concerning its morphology and taxonomy. However, its pantropical intraspecific genetic diversity and population differentiation is unknown. We initiated this study to better understand the genetic variation within E. scabrosa and to test if cryptic species are present. Sequences of three DNA regions, the nuclear ribosomal internal transcribed spacer (ITS), the large subunit ribosomal DNA (LSU), and the translation elongation factor (EF1α) were analysed for 66 samples from 15 geographical locations. We found a high level of genetic diversity (haplotype diversity, Hd = 0.88) and low nucleotide diversity (π = 0.006) across the known geographical range of E. scabrosa based on ITS sequences. The analysis of molecular variance (AMOVA) indicates that the genetic variability is mainly found among geographical populations. The results of Mantel tests confirmed that the genetic distance among populations of E. scabrosa is positively correlated with the geographical distance, which indicates that geographical isolation is an important factor for the observed genetic differentiation. Based on phylogenetic analyses of combined dataset ITS-LSU-EF1α, the low intraspecific divergences (0-0.3%), and the Automated Barcode Gap Discovery (ABGD) analysis, E. scabrosa can be considered as a single species with five different geographical populations. Each population might be in the process of allopatric divergence and in the long-term they may evolve and become distinct species.

Genetic Differentiation and Widespread Mitochondrial Heteroplasmy among Geographic Populations of the Gourmet Mushroom Thelephora ganbajun from Yunnan, China

The mitochondrial genomes are generally considered non-recombining and homoplasmic in nature. However, our previous study provided the first evidence of extensive and stable mitochondrial heteroplasmy in natural populations of the basidiomycete fungus Thelephora ganbajun from Yunnan province, China. The heteroplasmy was characterized by the presence of two types of introns residing at adjacent but different sites in the cytochrome oxidase subunits I (cox1) gene within an individual strain. However, the frequencies of these two introns among isolates from different geographical populations and the implications for the genetic structure in natural populations have not been investigated. In this study, we analyzed DNA sequence variation at the internal transcribed spacer (ITS) regions of the nuclear ribosomal RNA gene cluster among 489 specimens from 30 geographic locations from Yunnan and compared that variation with distribution patterns of the two signature introns in the cox1 gene that are indicative of heteroplasmy in this species. In our samples, evidence for gene flow, abundant genetic diversity, and genotypic uniqueness among geographic samples in Yunnan were revealed by ITS sequence variation. While there was insignificant positive correlation between geographic distance and genetic differentiation among the geographic samples based on ITS sequences, a moderate significant correlation was found between ITS sequence variation, geographical distance of sampling sites, and distribution patterns of the two heteroplasmic introns in the cox1 gene. Interestingly, there was a significantly negative correlation between the copy numbers of the two co-existing introns. We discussed the implications of our results for a better understanding of the spread of stable mitochondrial heteroplasmy, mito-nuclear interactions, and conservation of this important gourmet mushroom.

Dynamic evolution of eukaryotic mitochondrial and nuclear genomes: a case study in the gourmet pine mushroom Tricholoma matsutake

Fungi, as eukaryotic organisms, contain two genomes, the mitochondrial genome and the nuclear genome, in their cells. How the two genomes evolve and correlate to each other is debated. Herein, taking the gourmet pine mushroom Tricholoma matsutake as an example, we performed comparative mitogenomic analysis using samples collected from diverse locations and compared the evolution of the two genomes. The T. matsutake mitogenome encodes 49 genes and is rich of repetitive and non-coding DNAs. Six genes were invaded by up to 11 group I introns, with one cox1 intron cox1P372 showing presence/absence dynamics among different samples. Bioinformatic analyses suggested limited or no evidence of mitochondrial heteroplasmy. Interestingly, hundreds of mitochondrial DNA fragments were found in the nuclear genome, with several larger than 500 nt confirmed by PCR assays and read count comparisons, indicating clear evidence of transfer of mitochondrial DNA into the nuclear genome. Nuclear DNA of T. matsutake showed a higher mutation rate than mitochondrial DNA. Furthermore, we found evidence of incongruence between phylogenetic trees derived from mitogenome and nuclear DNA sequences. Together, our results reveal the dynamic genome evolution of the gourmet pine mushroom.

Exploring the Species Diversity of Edible Mushrooms in Yunnan, Southwestern China, by DNA Barcoding

Yunnan Province, China, is famous for its abundant wild edible mushroom diversity and a rich source of the world's wild mushroom trade markets. However, much remains unknown about the diversity of edible mushrooms, including the number of wild edible mushroom species and their distributions. In this study, we collected and analyzed 3585 mushroom samples from wild mushroom markets in 35 counties across Yunnan Province from 2010 to 2019. Among these samples, we successfully obtained the DNA barcode sequences from 2198 samples. Sequence comparisons revealed that these 2198 samples likely belonged to 159 known species in 56 different genera, 31 families, 11 orders, 2 classes, and 2 phyla. Significantly, 51.13% of these samples had sequence similarities to known species at lower than 97%, likely representing new taxa. Further phylogenetic analyses on several common mushroom groups including 1536 internal transcribed spacer (ITS) sequences suggested the existence of 20 new (cryptic) species in these groups. The extensive new and cryptic species diversity in wild mushroom markets in Yunnan calls for greater attention for the conservation and utilization of these resources. Our results on both the distinct barcode sequences and the distributions of these sequences should facilitate new mushroom species discovery and forensic authentication of high-valued mushrooms and contribute to the scientific inventory for the management of wild mushroom markets.

Genetic diversity and population structure of Meretrix petechialis in China revealed by sequence-related amplified polymorphism markers

Genetic variation in nine stocks of Meretrix petechialis collected from China was analyzed using sequence-related amplified polymorphism (SRAP) markers. Eight primer pairs produced 132 polymorphic loci with an average of 16.5 loci per primer pair. A population from Jiangsu had the highest percentage of polymorphic loci at 27.27%, suggesting that these resources had a rich genetic diversity. The Nei’s gene diversity of the nine populations ranged from 0.0647 to 0.0793; a population from Shandong was the lowest and a population from North Korea the highest. The Shannon’s information index was between 0.1023 and 0.1202, with the lowest in the Shandong population and the highest in the Jiangsu population. The Nei’s unbiased genetic distance between the nine populations was 0.0243–0.0570 and the genetic similarity was 0.9446–0.9760; the genetic distance between Guangxi and Shandong populations was the furthest (0.0570) and the genetic distance between Shandong and Jiangsu populations was the closest (0.0243). Nei’s gene diversity analysis indicated that the genetic variance was mainly found within individual geographical populations, and the analysis of molecular variance revealed low but significant genetic differentiation among local and regional populations. The limited gene flow (Nm = 0.555) was inferred as a major reason for the extent of genetic differentiation in M. petechialis. The results obtained here indicated that M. petechialis have high degree of genetic diversity and the potential of further breeding with excellent germplasm resources. This study provides a scientific basis for the protection of germplasm resources and the breeding of M. petechialis.

Identification of the Pol Gene as a Species-Specific Diagnostic Marker for Qualitative and Quantitative PCR Detection of Tricholoma matsutake

Tricholoma matsutake is a rare, precious, and wild edible fungus that could not be cultivated artificially until now. This situation has given way to the introduction of fake T. matsutake commodities to the mushroom market. Among the methods used to detect food adulteration, amplification of species-specific diagnostic marker is particularly important and accurate. In this study, the Pol gene is reported as a species-specific diagnostic marker to identify three T. matsutake varieties and 10 other types of edible mushrooms through qualitative and quantitative PCR. The PCR results did not reveal variations in the amplified region, and the detection limits of qualitative and quantitative PCR were found to be 8 ng and 32 pg, respectively. Southern blot showed that the Pol gene exists as a single copy in the T. matsutake genome. The method that produced the purest DNA of T. matsutake in this study was also determined, and the high-concentration salt precipitation method was confirmed to be the most suitable among the methods tested. The assay proposed in this work is applicable not only to the detection of raw materials but also to the examination of processed products containing T. matsutake.

SNP and SCAR Markers for Specific Discrimination of Antler-Shaped Ganoderma lucidum

In this study we identified single nucleotide polymorphism (SNP) and sequence characteristic amplification region (SCAR) markers for specific identification of antler-shaped Ganoderma lucidum strains. When the partial mitochondrial SSU rDNA gene sequence of various antler- and kidney-shaped G. lucidum strains were analyzed and aligned, an SNP was found only in the antler-shaped G. lucidum strain at position 456 bp. In addition, this SNP of antler-shaped strains was digested by HinfI restriction enzyme. We further analyzed the polymorphism of various G. lucidum strains by random amplified polymorphic DNA (RAPD) analysis. In RAPD analysis, we isolated and sequenced a fragment, specific for antler-shaped G. lucidum strains. Based on this specific fragment sequence, two sets of specific primer pairs for antler-shaped G. lucidum strains were designed. PCR analysis revealed that two specific bands were observed only from antler-shaped strains. These two molecular markers will be helpful for identification of morphological characteristics of G. lucidum.

Genetic diversity and distribution of Senegalia senegal (L.) Britton under climate change scenarios in West Africa

Climate change is predicted to impact species’ genetic diversity and distribution. We used Senegalia senegal (L.) Britton, an economically important species distributed in the Sudano-Sahelian savannah belt of West Africa, to investigate the impact of climate change on intraspecific genetic diversity and distribution. We used ten nuclear and two plastid microsatellite markers to assess genetic variation, population structure and differentiation across thirteen sites in West Africa. We projected suitable range, and potential impact of climate change on genetic diversity using a maximum entropy approach, under four different climate change scenarios. We found higher genetic and haplotype diversity at both nuclear and plastid markers than previously reported. Genetic differentiation was strong for chloroplast and moderate for the nuclear genome. Both genomes indicated three spatially structured genetic groups. The distribution of Senegalia senegal is strongly correlated with extractable nitrogen, coarse fragments, soil organic carbon stock, precipitation of warmest and coldest quarter and mean temperature of driest quarter. We predicted 40.96 to 6.34 per cent of the current distribution to favourably support the species’ ecological requirements under future climate scenarios. Our results suggest that climate change is going to affect the population genetic structure of Senegalia senegal, and that patterns of genetic diversity are going to influence the species’ adaptive response to climate change. Our study contributes to the growing evidence predicting the loss of economically relevant plants in West Africa in the next decades due to climate change.

Local-Level Genetic Diversity and Structure of Matsutake Mushroom (Tricholoma matsutake) Populations in Nagano Prefecture, Japan, Revealed by 15 Microsatellite Markers

The annual yield of matsutake mushrooms (Tricholoma matsutake) has consistently decreased in Japan over the past few decades. We used 15 polymorphic and codominant simple sequence repeat (SSR) markers, developed using next-generation sequencing, to carry out genetic analyses of 10 populations in Nagano, Japan. Using the SSRs, we identified 223 genotypes, none of which was observed in more than one population. The mean expected heterozygosity and standardized allelic richness values were 0.67 and 4.05, respectively. Many alleles appeared in only one of the 10 populations; 34 of these private alleles were detected with a mean number per population of 3.4. The fixation index (F_ST) and standardized genetic differentiation (G′_ST) values were 0.019 and 0.028, respectively. Analysis of molecular variance (AMOVA) showed that the contribution of among population, among genets within a population, and within genets variation to the total variation was 2.91%, 11.62%, and 85.47%, respectively, with genetic differentiation being detected for all sources. Twenty-eight of 45 pairwise F_ST values were significantly larger than zero, and no pattern of isolation by distance was detected among the 10 populations. Bayesian-based clustering did not show clear differences among populations. These results suggest that reestablishment of a colony would be best accomplished by transplantation within a field; if this is not possible, then transplantation from within several dozen kilometers will cause little damage to the original population genetic structure.

Prediction of the potential geographic distribution of the ectomycorrhizal mushroom Tricholoma matsutake under multiple climate change scenarios

Effective conservation and utilization strategies for natural biological resources require a clear understanding of the geographic distribution of the target species. Tricholoma matsutake is an ectomycorrhizal (ECM) mushroom with high ecological and economic value. In this study, the potential geographic distribution of T. matsutake under current conditions in China was simulated using MaxEnt software based on species presence data and 24 environmental variables. The future distributions of T. matsutake in the 2050s and 2070s were also projected under the RCP 8.5, RCP 6, RCP 4.5 and RCP 2.6 climate change emission scenarios described in the Special Report on Emissions Scenarios (SRES) by the Intergovernmental Panel on Climate Change (IPCC). The areas of marginally suitable, suitable and highly suitable habitats for T. matsutake in China were approximately 0.22 × 10⁶ km², 0.14 × 10⁶ km², and 0.11 × 10⁶ km², respectively. The model simulations indicated that the area of marginally suitable habitats would undergo a relatively small change under all four climate change scenarios; however, suitable habitats would significantly decrease, and highly suitable habitat would nearly disappear. Our results will be influential in the future ecological conservation and management of T. matsutake and can be used as a reference for studies on other ectomycorrhizal mushroom species.

A qPCR assay that specifically quantifies Tricholoma matsutake biomass in natural soil

Tricholoma matsutake is an ectomycorrhizal basidiomycete that produces prized, yet uncultivable, "matsutake" mushrooms along densely developed mycelia, called "shiro," in the rhizosphere of coniferous forests. Pinus densiflora is a major host of this fungus in Japan. Measuring T. matsutake biomass in soil allows us to determine the kinetics of fungal growth before and after fruiting, which is useful for analyzing the conditions of the shiro and its surrounding mycorrhizosphere, predicting fruiting timing, and managing forests to obtain better crop yields. Here, we document a novel method to quantify T. matsutake mycelia in soil by quantifying a single-copy DNA element that is uniquely conserved within T. matsutake but is absent from other fungal species, including close relatives and a wide range of ectomycorrhizal associates of P. densiflora. The targeted DNA region was amplified quantitatively in cultured mycelia that were mixed with other fungal species and soil, as well as in an in vitro co-culture system with P. densiflora seedlings. Using this method, we quantified T. matsutake mycelia not only from shiro in natural environments but also from the surrounding soil in which T. matsutake mycelia could not be observed by visual examination or distinguished by other means. It was demonstrated that the core of the shiro and its underlying area in the B horizon are predominantly composed of fungal mycelia. The fungal mass in the A or A₀ horizon was much lower, although many white mycelia were observed at the A horizon. Additionally, the rhizospheric fungal biomass peaked during the fruiting season.

Evidence for Inbreeding and Genetic Differentiation among Geographic Populations of the Saprophytic Mushroom Trogia venenata from Southwestern China

During the past 40 years, more than 400 Sudden Unexplained Deaths (SUDs) have occurred in Yunnan, southwestern China. Epidemiological and toxicological analyses suggested that a newly discovered mushroom called Trogia venenata was the leading culprit for SUDs. At present, relatively little is known about the genetics and natural history of this mushroom. In this study, we analyzed the sequence variation at four DNA fragments among 232 fruiting bodies of T. venenata collected from seven locations. Our ITS sequence analyses confirmed that all the isolates belonged to the same species. The widespread presence of sequence heterozygosity within many strains at each of three protein-coding genes suggested that the fruiting bodies were diploid, dikaryotic or heterokaryotic. Within individual geographic populations, we found significant deviations of genotype frequencies from Hardy-Weinberg expectations, with the overall observed heterozygosity lower than that expected under random mating, consistent with prevalent inbreeding within local populations. The geographic populations were overall genetically differentiated. Interestingly, while a positive correlation was found between population genetic distance and geographic distance, there was little correlation between genetic distance and barium concentration difference for the geographic populations. Our results suggest frequent inbreeding, geographic structuring, and limited gene flow among geographic populations of T. venenata from southwestern China.

Genetic variability and bottleneck detection of four Tricholoma matsutake populations from northeastern and southwestern China

The excessive commercial collection of matsutake mushrooms can lead to extreme reduction of population size, which may cause genetic bottleneck and decrease genetic diversity of Tricholoma matsutake. Here, six polymorphic microsatellite loci markers were used to examine the genetic diversity of four natural T. matsutake populations from two main producing regions of China. The minimum combinations of four loci were able to discriminate total 86 sampled individuals with distinctive multilocus genotypes. Our analysis of molecular variance (AMOVA) revealed that about 80% and 20% of the overall genetic variation were respectively partitioned within and among populations. The principal-coordinate analyses (PCA) distinguished the four tested populations into three genetic clusters, each of which was correlated with respective endemic host plants on a geographical basis. The AMOVA, PCA and pairwise population FST estimates consistently displayed the same genetic divergence patterns and spatial structure of T. matsutake mediated by host plants in China. The significant heterozygosity excesses demonstrated that a recent genetic bottleneck occurred in each population tested. The complementary M-ratio test indicated past genetic bottleneck events over longer periods. Only four individuals were identified as putative first generation migrants within northeastern China, which implies restricted interpopulation gene flow in T. matsutake. We discuss that the significant genetic differentiation among populations of T. matsutake is most likely a function of host adaptation, host specificity, genetic bottleneck, limited dispersal and habitat fragmentation.

The host ranges of conifer-associated Tricholoma matsutake, Fagaceae-associated T. bakamatsutake and T. fulvocastaneum are wider in vitro than in nature

Tricholoma matsutake is the most commercially important edible mushroom in pine forests in Japan. Tricholoma bakamatsutake and T. fulvocastaneum, species closely related to T. matsutake, occur in Fagaceae forests. We examined ectomycorrhizal (EM) formation by these Tricholoma species by in vitro synthesis among seven strains (two of T. matsutake, four of T. bakamatsutake, one of T. fulvocastaneum) and axenic plants of pine (Pinus densiflora) and oak (Quercus serrata, Q. phillyraeoides). All strains, except for one of T. matsutake, formed EM associations with both pine and oak. Plant growth and mycelial development were differently affected by EM formation depending on the plant-fungus combination.

Little to no genetic structure in the ectomycorrhizal basidiomycete Suillus spraguei (Syn. S. pictus) across parts of the Northeastern USA

Population genetic studies of ectomycorrhizal fungi at scales larger than 100 km are still relatively rare with highly variable results. In this study, we determined the population genetic structure of the ectomycorrhizal basidiomycete Suillus spraguei over distances up to 600 km in northeastern USA forests. S. spraguei associates exclusively with five-needled pines and only with white pine (Pinus strobus) in the eastern USA. We used six microsatellite loci to assess the genetic structure between eight sites sampled in the Adirondack Park of New York and seven sites sampled in other forests of New York, Pennsylvania, and Massachusetts. Except for one site, little to no genetic differentiation was detected in pairwise comparisons of the sites (F ST = 0 to 0.05). Only one site was moderately differentiated from most other sites (F ST = 0.02 to 0.15). The Mantel test showed no significant correlation between genetic and geographic distances (isolation by distance; R (2) = 0.003, P = 0.3). The STRUCTURE analysis also supported the presence of a single cluster (K = 1).

Strategies of functional food for cancer prevention in human beings

Functional food for prevention of chronic diseases is one of this century's key global challenges. Cancer is not only the first or second leading cause of death in China and other countries across the world, but also has diet as one of the most important modifiable risk factors. Major dietary factors now known to promote cancer development are polished grain foods and low intake of fresh vegetables, with general importance for an unhealthy lifestyle and obesity. The strategies of cancer prevention in human being are increased consumption of functional foods like whole grains (brown rice, barley, and buckwheat) and by-products, as well some vegetables (bitter melon, garlic, onions, broccoli, and cabbage) and mushrooms (boletes and Tricholoma matsutake). In addition some beverages (green tea and coffee) may be protective. Southwest China (especially Yunnan Province) is a geographical area where functional crop production is closely related to the origins of human evolution with implications for anticancer influence.

Genetic variability and population structure of the mushroom Pleurotus eryngii var. tuoliensis

The genetic diversity of 123 wild strains of Pleurotus eryngii var. tuoliensis, which were collected from nine geographical locations in Yumin, Tuoli, and Qinghe counties in the Xinjiang Autonomous Region of China, was analysed using two molecular marker systems (inter-simple sequence repeat and start codon targeted). At the variety level, the percentage of polymorphic loci and Nei’s gene diversity index for P. eryngii var. tuoliensis was 96.32% and 0.238, respectively. At the population level, Nei’s gene diversity index ranged from 0.149 to 0.218 with an average of 0.186, and Shannon's information index ranged from 0.213 to 0.339 with an average of 0.284. These results revealed the abundant genetic variability in the wild resources of P. eryngii var. tuoliensis. Nei’s gene diversity analysis indicated that the genetic variance was mainly found within individual geographical populations, and the analysis of molecular variance revealed low but significant genetic differentiation among local and regional populations. The limited gene flow (Nm = 1.794) was inferred as a major reason for the extent of genetic differentiation of P. eryngii var. tuoliensis. The results of Mantel tests showed that the genetic distance among geographical populations of P. eryngii var. tuoliensis was positively correlated with the geographical distance and the longitudinal distances (r_Go = 0.789 and r_Ln = 0.873, respectively), which indicates that geographical isolation is an important factor for the observed genetic differentiation. Nine geographical populations of P. eryngii var. tuoliensis were divided into three groups according to their geographical origins, which revealed that the genetic diversity was closely related to the geographical distribution of this wild fungus.

Genetic diversity and recombination in natural populations of the nematode-trapping fungus Arthrobotrys oligospora from China

Nematophagous fungi can trap and capture nematodes and other small invertebrates. This unique ability has made them ideal organisms from which to develop biological control agents against plant- and animal-parasitic nematodes. However, effective application of biocontrol agents in the field requires a comprehensive understanding about the ecology and population genetics of the nematophagous fungi in natural environments. Here, we genotyped 228 strains of the nematode-trapping fungus Arthrobotrys oligospora using 12 single nucleotide polymorphic markers located on eight random DNA fragments. The strains were from different ecological niches and geographical regions from China. Our analyses identified that ecological niche separations contributed significantly, whereas geographic separation contributed relatively little to the overall genetic variation in our samples of A. oligospora. Interestingly, populations from stressful environments seemed to be more variable and showed more evidence for recombination than those from benign environments at the same geographic areas. We discussed the implications of our results to the conservation and biocontrol application of A. oligospora in agriculture and forestry.

Extensive gene flow over Europe and possible speciation over Eurasia in the ectomycorrhizal basidiomycete Laccaria amethystina complex

Biogeographical patterns and large-scale genetic structure have been little studied in ectomycorrhizal (EM) fungi, despite the ecological and economic importance of EM symbioses. We coupled population genetics and phylogenetic approaches to understand spatial structure in fungal populations on a continental scale. Using nine microsatellite markers, we characterized gene flow among 16 populations of the widespread EM basidiomycete Laccaria amethystina over Europe (i.e. over 2900 km). We also widened our scope to two additional populations from Japan (10(4) km away) and compared them with European populations through microsatellite markers and multilocus phylogenies, using three nuclear genes (NAR, G6PD and ribosomal DNA) and two mitochondrial ribosomal genes. European L. amethystina populations displayed limited differentiation (average F(ST) = 0.041) and very weak isolation by distance (IBD). This panmictic European pattern may result from effective aerial dispersal of spores, high genetic diversity in populations and mutualistic interactions with multiple hosts that all facilitate migration. The multilocus phylogeny based on nuclear genes confirmed that Japanese and European specimens were closely related but clustered on a geographical basis. By using microsatellite markers, we found that Japanese populations were strongly differentiated from the European populations (F(ST) = 0.416), more than expected by extrapolating the European pattern of IBD. Population structure analyses clearly separated the populations into two clusters, i.e. European and Japanese clusters. We discuss the possibility of IBD in a continuous population (considering some evidence for a ring species over the Northern Hemisphere) vs. an allopatric speciation over Eurasia, making L. amethystina a promising model of intercontinental species for future studies.

Analysis of genetic diversity of Lactarius hatsudake in South China

Lactarius hatsudake is a type of ectomycorrhizal fungus that significantly influences the growth of pine trees. It is widely prevalent in Asian countries and has a high economic value. Artificial cultivation of this fungus has not been achieved as yet; therefore, excessive manual harvesting may cause serious damages to the site of its production. In this study, we analyzed 41 samples of L. hatsudake from south China using internal transcribed spacer (ITS) sequences. By comparing the differences among ITS sequences to identify the haplotype diversity within each population, the relationships among local populations, the relationship between the level of genetic differentiation and geographical separation, and the contributions of local and regional geographical separations to the overall ITS haplotype variation were analyzed. Genetic analysis indicates that ITS sequences obtained from these 41 L. hatsudake samples could be identified as 18 haplotypes, of which 13 haplotypes were contained in only a single sample, whereas the remaining sequence types all were contained in two or more samples. The most common sequence type, haplotype 6, was found in 16 samples and was distributed across nearly every region. The Mantel test demonstrated that there is no significant linear relationship between geographical distance and the FST value of genetic difference. Results of this research illustrates that there exists a certain degree of genetic intermixing among natural populations of L. hatsudake. From the group genetic analysis, it appears that there exists genetic differentiation of lower frequencies in natural populations of L. hatsudake; however, the linear relationship between the degree of genetic differentiation and geographical distance is not distinctly apparent.

Population genetics of ectomycorrhizal fungi: from current knowledge to emerging directions

Ectomycorrhizal (EM) fungi are major microbial components of boreal, temperate and Mediterranean forests, as well as some tropical forest ecosystems. Nearly two decades of studies have clarified many aspects of their population biology, based on several model species from diverse lineages of fungi where the EM symbiosis evolved, i.e. among Hymenomycetes and, to a lesser extent, among Ascomycetes. In this review, we show how tools for individual recognition have changed, shifting from the use of somatic incompatibility reactions to dominant and non-specific markers (such as random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP)) and, more recently, to co-dominant and specific markers (such as microsatellites and single nucleotide polymorphisms (SNPs)). At the same time, the theoretical focus has also changed. In earlier studies, a major aim was the description of genet size and popul/ation strategy. For example, we show how some studies supported or challenged the simple, classical model of colonization of new forest stands by ruderal (R) species, propagating by spores and forming small genets, progressively replaced in older forests by more competitive (C) species, propagating by mycelial growth and forming larger genets. By contrast, more recent studies give insights into some genetic traits, such as partners' assortment (allo- versus autogamy), genetic structure of populations and gene flow that turn out to depend both on distance and on whether spores are animal- or wind-dispersed. We discuss the rising awareness that (i) many morphospecies contain cryptic biological species (often sympatric) and (ii) trans- and inter-continental species may often contain several biological species isolated by distance. Finally, we show the emergence of biogeographic approaches and call for some aspects to be developed, such as fine-scale and long-term population monitoring, analyses of subterranean populations of extra-radical mycelia, or more model species from the tropics, as well as from the Ascomycetes (whose genetic idiosyncrasies are discussed). With the rise of the '-omics' sciences, analysis of population structure for non-neutral genes is expected to develop, and forest management and conservation biology will probably profit from published and expected work.

Genetic population structure of the ectomycorrhizal fungus Pisolithus microcarpus suggests high gene flow in south-eastern Australia

Pisolithus are ectomycorrhizal fungi that associate with roots of numerous plant species in natural and plantation forests worldwide. Despite the fact that Pisolithus spp. are present in plantation forests in many countries, knowledge of the genetic population structure of Pisolithus spp. remains limited. In this study, we have tested the hypothesis that a propensity for long-distance spore dispersal in Pisolithus microcarpus, along with the widespread distribution of potential eucalypt and acacia plant hosts in south-eastern Australia facilitates gene flow that limits population differentiation. Five polymorphic simple sequence repeat markers were used to investigate the population structure of P. microcarpus. Isolates were grouped according to geographical origin and isolate genotypes were analysed among the geographical populations. Pairwise F (ST) estimates indicated limited genetic differentiation among the geographical populations. Analysis of molecular variance revealed that most of the genetic variation present was within geographical populations, with only 1.3% of the genetic variation among P. microcarpus geographical populations. This was particularly pronounced for four geographical populations within a ca 7,000 km(2) area New South Wales, which were each separated by < 100 km and appeared to be genetically homogeneous. The lack of population structure is suggested to be due to a high degree of gene flow, via basidiospores, between the New South Wales geographical populations.

Structure, function, and phylogeny of the mating locus in the Rhizopus oryzae complex

The Rhizopus oryzae species complex is a group of zygomycete fungi that are common, cosmopolitan saprotrophs. Some strains are used beneficially for production of Asian fermented foods but they can also act as opportunistic human pathogens. Although R. oryzae reportedly has a heterothallic (+/-) mating system, most strains have not been observed to undergo sexual reproduction and the genetic structure of its mating locus has not been characterized. Here we report on the mating behavior and genetic structure of the mating locus for 54 isolates of the R. oryzae complex. All 54 strains have a mating locus similar in overall organization to Phycomyces blakesleeanus and Mucor circinelloides (Mucoromycotina, Zygomycota). In all of these fungi, the minus (-) allele features the SexM high mobility group (HMG) gene flanked by an RNA helicase gene and a TP transporter gene (TPT). Within the R. oryzae complex, the plus (+) mating allele includes an inserted region that codes for a BTB/POZ domain gene and the SexP HMG gene. Phylogenetic analyses of multiple genes, including the mating loci (HMG, TPT, RNA helicase), ITS1-5.8S-ITS2 rDNA, RPB2, and LDH genes, identified two distinct groups of strains. These correspond to previously described sibling species R. oryzae sensu stricto and R. delemar. Within each species, discordant gene phylogenies among multiple loci suggest an outcrossing population structure. The hypothesis of random-mating is also supported by a 50:50 ratio of plus and minus mating types in both cryptic species. When crossed with tester strains of the opposite mating type, most isolates of R. delemar failed to produce zygospores, while isolates of R. oryzae produced sterile zygospores. In spite of the reluctance of most strains to mate in vitro, the conserved sex locus structure and evidence for outcrossing suggest that a normal sexual cycle occurs in both species.

In vitro ectomycorrhizal specificity between the Asian red pine Pinus densiflora and Tricholoma matsutake and allied species from worldwide Pinaceae and Fagaceae forests

Tricholoma matsutake produces commercially valuable, yet uncultivable, mushrooms (matsutake) in association with pines in the Far East and Scandinavia and with both pines and oaks in the foothills of Tibet. Other matsutake mushrooms, such as Tricholoma anatolicum from the Mediterranean regions and Tricholoma magnivelare and Tricholoma sp. from the North Pacific Coast area of Canada and North America as well as Mexico, respectively, are associated with pines or oaks in their natural habitats. Tricholoma bakamatsutake and Tricholoma fulvocastaneum from Asia produce moderately valuable matsutake mushrooms and are solely associated with Fagaceae in nature. In this study, we demonstrate for the first time that matsutake mushrooms from Scandinavia, Mediterranean regions, North America, and Tibet form ectomycorrhizae with Pinus densiflora similar to the Far East T. matsutake. In general, worldwide T. matsutake and the symbionts of Pinaceae colonize the rhizospheres of P. densiflora as well as T. matsutake isolated from the host plant. However, T. fulvocastaneum and T. bakamatsutake formed a discontinuous Hartig net and no Hartig net, respectively, and colonized to a lesser extent as compared to T. matsutake. The data suggest that conifer-associated matsutake mushrooms in their native habitat will associate symbiotically with the Asian red pine.

Genetic diversity of Dahongjun, the commercially important "Big Red Mushroom" from southern China

BACKGROUND

In southern China, a wild ectomycorrhizal mushroom commonly called "Dahongjun" or "Big Red Mushroom" by the local residents, has been harvested, consumed, and/or exported as an exotic food for many years. Although ecologically and economically important, very little is known about this mushroom, including its diversity and population structure.

METHODOLOGY AND PRINCIPAL FINDINGS

In this study, we analyzed 122 samples from five local populations representing the known distribution ranges of this mushroom in southern China. We investigated the genetic diversity and geographic structure of this mushroom using sequences from four DNA fragments. Our analyses identified that this mushroom contained at least three divergent lineages: one corresponds to a recently described species Russula griseocarnosa from southern China and the remaining two likely represent two novel species. While these lineages were prominently structured geographically based on ITS sequences, evidence for ancient and/or recent gene flow was also identified within individual lineages. In addition, a local population from Ailaoshan in central Yunnan Province where 85 of our 122 specimens came from showed clear evidence of recombination.

CONCLUSION AND SIGNIFICANCE

The ectomycorrhizal mushroom "Dahongjun" from southern China is a species complex with at least three divergent lineages. These lineages are largely geographically structured and there is evidence for recombination in nature. Our results indicate mature Dahongjun mushrooms with abundant basidiospores are important for the reproduction of this mushroom in nature and that individual populations of this species should be managed separately.

Genetic diversity of the endemic gourmet mushroom Thelephora ganbajun from south-western China

The ectomycorrhizal fungus Thelephora ganbajun is an endemic gourmet mushroom in Yunnan province, south-western China. However, despite its widespread consumer appeal, nutritional value and potential ecological role in natural forests, very little is known about its genetics, diversity and ecology. In this study, we investigated DNA sequence variation at the internal transcribed spacer (ITS) regions among 156 specimens collected from 23 sites of nine regions in Yunnan Province. Our analysis identified a total of 34 ITS haplotypes and these haplotypes were clustered into five distinct phylogenetic groups. The evolutionary divergences among these clades are similar to or greater than many known sister species pairs within the genus Thelephora and the closely related genus Tomentella. Among the 34 ITS haplotypes, 22 were represented by one specimen each and the remaining 12 were each shared by two or more specimens. The most common haplotype contained 68 specimens distributed in 21 of the 23 sites, a result consistent with gene flow among geographical populations. However, analysis of molecular variance (AMOVA) revealed low but significant genetic differentiation among local and regional populations. Interestingly, the Mantel test identified that the extent of genetic differentiation was not significantly correlated with geographical distance. Our study revealed significant genetic divergence within Th. ganbajun and limited but detectable gene flow among geographical populations of this endemic ectomycorrhizal gourmet mushroom.