Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2005, p. 6584-6589, Vol. 71, No. 11
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.11.6584-6589.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
National Research Council, Plant Genetics Institute, Perugia Division, Via Madonna Alta 130, I-06128 Perugia, Italy,1 National Research Council, Plant Genetics Institute, Florence Division, Via Madonna del Piano, I-50019 Sesto Fiorentino, Florence, Italy2
Received 9 December 2004/ Accepted 13 July 2005
|
|
|---|
|
|
|---|
The complexity of the truffle life cycle, the difficulties of growing these fungi under controlled conditions, and the lack of reliable phenotypic markers to differentiate morphologically similar species have been major obstacles to understanding the distribution, propagation, and fructification of these hypogeous fungi. Molecular markers have been developed to type most of the economically important truffle species (2, 10, 15, 19, 23, 24, 25, 33, 34), but studies of the environmental and molecular determinants for different portions of the life cycle generally are lacking. The spatial distribution and ecological requirements of these symbiotic fungi vary by species. Some species are widely distributed and have pronounced morphological and molecular variability (9, 11, 22, 25, 27, 32), while others have a more restricted distribution and little intraspecific polymorphism in either morphological or genetic traits (3, 4, 8, 9, 26). Tuber magnatum Pico, which produces white ascomata, is harvested only in Italy and some countries on the Balkan Peninsula.
In the related species T. melanosporum, there is little or no genetic differentiation between the sampled populations, and the species is regarded as one that self-fertilizes for sexual reproduction and that has an ascocarp consisting primarily of diploid (dikaryotic) hyphae (3, 4, 21). The few studies of genetic structure and mode of reproduction of T. magnatum found very limited intraspecific allozyme variation and a complete absence of heterozygous individuals (8). Whether T. magnatum is normally dikaryotic and reproduces homothallically or by self-fertilization or if it is an outbreeding ascomycete remains to be determined.
Our objectives in this study were to use microsatellite markers to determine (i) if T. magnatum populations are genetically homogeneous, (ii) if the fungus is dikaryotic, and (iii) the degree of outbreeding in field populations. Understanding the life cycle and genetic differences between truffle populations is important for commercial truffle growth management and marketing strategies.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Name and location of the studied populations and estimated diversity parameters
|
Data analysis.
The mean number of alleles and the allelic richness, over all loci, were calculated for each population. Allelic richness, computed for each population, was weighted to the smallest sample size (n = 5) using a rarefaction method (29).
Two-locus linkage disequilibrium analysis was performed with GENEPOP, version 3.2a (30), within each population and each geographic group and across all populations. The index of association (IA) used to test for multilocus linkage disequilibrium was calculated as described by Maynard Smith et al. (18) with the program MULTILOCUS, version 1.2 (1). The same software was used to determine if IA differed significantly from zero, by performing 1,000 randomizations of the data set. The randomization procedure was first performed by shuffling alleles among all individuals. To determine if linkage disequilibrium results from the reproductive isolation of geographically distant populations, randomizations by shuffling alleles among individuals within populations and within geographical groups also were performed (17).
To test the effects of sample size reduction on linkage disequilibrium (LD) values, two locus and multilocus analyses also were performed by considering a data set of 275 individuals pooled as a single population and obtained by excluding the southernmost and northwesternmost populations (3, 4, 25, and 26).
The extent of population subdivision was evaluated by calculating Wright's FST (39) and the related RST (37) indexes, with the latter taking explicit account of the mutation process that can occur at microsatellite loci. These two indexes were estimated across populations and between pairs of populations, according to Weir and Cockerham (38) and Michalakis and Excoffier (20), by using the program SPAGEDI, version 1.1 (14).
A genetic distance matrix of pairwise FST values also was used to perform a hierarchical analysis of molecular variance (AMOVA) (7) with the ARLEQUIN program (version 2.0; Genetics and Biometry Laboratory, University of Geneva, Switzerland [http://lgb.unige.ch/arlequin/]).
The possible presence of geographic structure of genetic variation in nuclear microsatellites in truffles was evaluated with four tests. First, we tested for the presence of phylogeographic structure by comparing RST estimates to RST values computed after 10,000 random permutations of allele types among alleles (13) were calculated using the SPAGEDI program. If RST is >RST (permuted), then there is phylogeographic structure, i.e., on average, phylogenetically similar alleles are found in the same population more often than are randomly chosen alleles.
Second, we tested for a pattern of isolation by distance (31). A Mantel test with 1,000 random permutations was performed with the matrix of pairwise genetic differentiation between populations, using RST/(1 RST), and a matrix of the ln(geographic distance) with the software GENALEX (version 5.1; Australian National University, Canberra, Australia [http://anu.edu.au/BoZo/GenAlEx/]).
Third, a spatial autocorrelation analysis was performed with SPAGEDI. This analysis used RST and evaluated 18 spatial distance classes with similar sample sizes. The 95% confidence intervals were estimated for each distance class by 10,000 random permutations.
Finally, a simulated annealing procedure implemented in the spatial analysis of molecular variance (SAMOVA) algorithm (6) was used to define groups of populations that are geographically homogeneous and maximally differentiated from each other. The program repeatedly seeks the composition of a user-defined number, K, of groups of geographically adjacent populations that maximizes FCT, the proportion of total genetic variance due to differences among groups of populations. The program was run for 10,000 iterations for K
{2, ..., 12} from each of 250 random initial conditions. The simulated annealing process was repeated 250 times to ensure that the final configuration of the K groups is not affected by a given initial configuration (6). For each K, the configuration with the largest FCT values after the 250 independent simulated annealing processes were completed was retained as the best grouping of populations.
|
|
|---|
There also was a significant departure from linkage equilibrium (IA = 0.109, P < 0.001) in the multilocus analysis. Linkage disequilibrium can result from reduced gene flow due to geographic separation. To assess the potential geographic bias in our combined analysis, these tests also were performed with samples in smaller geographic groups. When the two-locus LD was calculated for each of the 13 geographic groups, only one of the two-locus combinations had a significant LD, and no significant two-locus LD was detected when the 26 local populations were considered. Similarly, IA was not significant when the randomizations were performed by shuffling alleles among samples within the 13 geographical groups and within the 26 local populations (P values of 0.065 and 0.119, respectively).
To evaluate the effect of sample size on LD values, the two-locus LD and IA tests also were calculated for a data set of 275 individuals obtained by excluding from the whole data set the populations 3, 4, 25, and 26 that SAMOVA proved to be genetically differentiated (see below) and thus were expected to contribute largely to LD. In this situation, only 5 out of 21 two-locus combinations resulted in significant LD (P = 0.05), and the IA was not significant (IA = 0.036, P = 0.062), suggesting that the decrease in LD values when populations are analyzed separately is not due to the reduced sample size.
Patterns within and among population variability.
The seven SSR loci had 2 to 18 alleles per locus. Within each of the 351 truffles analyzed, only a single allele per locus was detected. Each multilocus genotype was then treated as a haplotype on the assumption that T. magnatum is a haploid fungus (see Discussion). The 26 populations had a mean number of alleles of 3.03 ± 0.12, ranging from 1.86 ± 0.34 to 4.43 ± 0.97, while the average allelic richness was 2.67 ± 0.06, ranging from 1.74 ± 0.26 to 3.30 ± 0.65 (Table 1).
Based on AMOVA, most of the genetic variability is located at the within-population level, which accounts for about 84% of the total molecular variance, while 6% and 10% of the genetic variation were attributable to differences among populations within geographical groups and among the 13 geographic groups, respectively (Table 2). Overall, 263 haplotypes were detected within the 316 samples. Some haplotypes were found at multiple sites in the same region, and some were found in multiple populations or groups (see supplemental Table SA1).
|
View this table: [in a new window] |
TABLE 2. AMOVA based on 316 truffles using seven SSR loci
|
Most of the alleles were not unique to a population, but some alleles could be used to characterize truffles according to their geographic origin because of their differential frequencies in some groups of populations. For example, the frequency of allele 17 at the MA7 locus was high in southern populations 23, 24, 25, and 26 (0.2, 0.4, 0.4, and 1.0, respectively), while alleles 14, 24, 34, and 35 of MA4 were found only in individuals from populations 3, 4, and 5 from Piemonte, the northwesternmost area of sampling.
A marked and highly significant difference between populations was suggested by the FST and RST values, which were 0.15 and 0.22, respectively. A phylogeographic signal also was detected with permutation procedures in the total sample [RST of 0.22 versus RST(permuted) of 0.12, P = 0.0005]. The test of isolation by distance showed that among-population differentiation increased significantly with ln(geographical distance) (Mantel test; P < 0.006), although the regression explains only 5.4% of the total variance (R2 = 0.054) (Fig. 1A).
![]() View larger version (19K): [in a new window] |
FIG. 1. (A) Correlation between genetic and geographical pairwise distances. (B) Average genetic distance (RST) for 18 spatial distance classes. Broken lines indicate the boundaries of the 95% confidence intervals estimated by permutation.
|
![]() View larger version (25K): [in a new window] |
FIG. 2. Sampling localities and geographic distribution of genetic variability. Numbers refer to populations listed in Table 1. Diameters of the circles are proportional to the sample sizes. Differently shaded circles show genetically different population groups according to SAMOVA based on RST index and a K of 4.
|
|
|
|---|
We used codominant markers (polymorphic microsatellite loci) and strains sampled across the entire known range of the species to infer the genetic structure, population dynamics, and reproductive mode of T. magnatum, the most prestigious of the edible species in the genus Tuber. We found only a single allele per locus in each of the 351 truffles analyzed. These results are consistent with those of Bertault et al. (3) for T. melanosporum. They previously interpreted the absence of heterozygotes to mean that these fungi are functionally diploid with a strictly self-fertilizing reproductive system. We interpret these data to mean that the ascocarps of T. magnatum are haploid maternal tissue, as found for many ascomycetes. Two-locus and multilocus linkage disequilibrium analyses indicate that within both the T. magnatum geographical groups and the populations, extensive genetic exchange occurs. The relatively high number of multilocus genotypes within populations suggests that these fungi are not highly inbred. Within the same truffle grounds, differences between haplotypes usually occurred at multiple loci rather than at a single locus (see supplemental Table SA1), a pattern consistent with recombination in addition to mutation as a source for the variation observed in the haplotypes.
The truffle ascocarps are presumed to be primarily maternal tissue and thus would be haploid. Paternal DNA might not be as easily recoverable since it would be present only in the ascospores, and these spores are not usually broken during the DNA extraction process that we use. Thus, haploid maternal tissue will be the dominant component of the truffle and presumably represents the genotypes we analyzed in this study.
Significant linkage disequilibrium was found only when the entire truffle data set was evaluated as a whole and could be explained if there was significant geographic population subdivision. The FST and RST statistics and the correlation between genetic and geographic distance are consistent with genetic differentiation of populations, while the results of the RST permutation analysis are consistent with phylogeographic structure. If individual populations are removed from the analysis, then removal of the southernmost and two northwesternmost populations allows the rest of the populations to be grouped together.
The geographic distribution of T. magnatum tracks the postglacial expansion of the species with which this fungus can establish a mutualistic symbiosis. Refugia for some T. magnatum host species (e.g., Quercus, Corylus, Tilia, and Carpinus) were located in central and southern Italy and spread northward as the ice receded (5, 28). A refugium for T. magnatum in central Italy, from which the northernmost and southernmost populations originated, could explain the results observed.
The genetic variation and geographic distribution of the SSR markers could enable the tracing and identification of T. magnatum populations. The geographic origin of truffles often is important to associations of truffle harvesters and local governments which are promoting economic and social development of rural and marginal areas. The genetic structure of populations of T. magnatum also has important ecological implications. Artificial fungal plantations often are established in naturally productive areas to counterbalance a sharp decline in wild truffle harvests. Potential problems related to microbial competition are largely unexplored with respect to truffle cultivation, and the consequences of the deliberate introduction of nonnative strains into an indigenous population have been serious for other edible mushrooms (16, 36). Hence, the profitable commercial practice of artificial fungal plantations may cause a loss of fungal biodiversity, and local strains with unique commercial characteristics might become extinct. Thus, when white truffle cultivation is planned in naturally productive areas, ecologically adapted truffle strains that are similar to those already present in the region should be used. Similar cautions also apply to T. melanosporum (21). The SSR markers that we used could be used to type strains of T. magnatum that are inoculated on commercially produced seedlings as well as for assessing the impact the introduced strains might have on native ectomycorrhizal populations.
The presence of genetic structure in populations of both of the most prestigious truffle species argues against the hypothesis that phenotypic and organoleptic differences among truffles of different origin are due solely to environmental factors (3). If the SSR loci examined are linked to the genes controlling these commercially important traits, then the SSR markers could be used to select strains with desirable characteristics for future commercial use. As our data suggest that truffles are not strictly self-fertilizing, the possibility of improving these mushrooms through conventional breeding protocols needs to be evaluated with more care and in greater detail. This study provides the first evidence suggesting that truffles are not strictly self-fertilizing organisms and should result in a significant reevaluation of the life cycle and biology of T. magnatum and other Tuber spp.
This research was supported in part by a grant from the Regione Umbria/Comunità Montane Umbre/CNR-IGV Sezione di Perugia, contribution no. 58 from the Institute of Plant Genetics, Division of Perugia.
Andrea Rubini and Francesco Paolocci contributed equally to this work.
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»