Fine-scale geographic difference of the endangered Big-headed Turtle (Platysternon megacephalum) fecal microbiota, and comparison with the syntopic Beale's Eyed Turtle (Sacalia bealei).

BACKGROUND: Studies have elucidated the importance of gut microbiota for an organism, but we are still learning about the important influencing factors. Several factors have been identified in helping shape the microbiome of a host, and in this study we focus on two factors-geography and host. We characterize the fecal microbiota of the Big-headed Turtle (Platysternon megacephalum) and compare across a relatively fine geographic scale (three populations within an 8-km radius) and between two syntopic hosts (P. megacephalum and Sacalia bealei). Both species are endangered, which limits the number of samples we include in the study. Despite this limitation, these data serve as baseline data for healthy, wild fecal microbiotas of two endangered turtle species to aid in conservation management. RESULTS: For geography, the beta diversity of fecal microbiota differed between the most distant sites. The genus Citrobacter significantly differs between sites, which may indicate a difference in food availability, environmental microbiota, or both. Also, we identify the common core microbiome for Platysternon across Hong Kong as the shared taxa across the three sites. Additionally, beta diversity differs between host species. Since the two species are from the same site and encounter the same environmental microbiota, we infer that there is a host effect on the fecal microbiota, such as diet or the recruitment of host-adapted bacteria. Lastly, functional analyses found metabolism pathways (KEGG level 1) to be the most common, and pathways (KEGG level 3) to be statistically significant between sites, but statistically indistinguishable between species at the same site. CONCLUSIONS: We find that fecal microbiota can significantly differ at a fine geographic scale and between syntopic hosts. Also, the function of fecal microbiota seems to be strongly affected by geographic site, rather than species. This study characterizes the identity and function of the fecal microbiota of two endangered turtle species, from what is likely their last remaining wild populations. These data of healthy, wild fecal microbiota will serve as a baseline for comparison and contribute to the conservation of these two endangered species.

Diversity and abundance of human-pathogenic fungi associated with pigeon faeces in urban environments.

Pathogenic fungi are a growing health concern worldwide, particularly in large, densely populated cities. The dramatic upsurge of pigeon populations in cities has been implicated in the increased incidence of invasive fungal infections. In this study, we used a culture-independent, high-throughput sequencing approach to describe the diversity of clinically relevant fungi (CRF) associated with pigeon faeces and map the relative abundance of CRF across Seoul, Korea. In addition, we tested whether certain geographical, sociological and meteorological factors were significantly associated with the diversity and relative abundance of CRF. Finally, we compared the CRF diversity of fresh and old pigeon faeces to identify the source of the fungi and the role of pigeons in dispersal. Our results demonstrated that both the composition and relative abundance of CRF are unevenly distributed across Seoul. The green area ratio and the number of multiplex houses were positively correlated with species diversity, whereas wind speed and number of households were negatively correlated. The number of workers and green area ratio were positively correlated with the relative abundance of CRF, whereas wind speed was negatively correlated. Because many CRF were absent in fresh faeces, we inferred that most species cannot survive the gastrointestinal tract of pigeons and instead are likely transmitted through soil or air and use pigeon faeces as a substrate for proliferation.

Diversity and phylogeography of Northeast Asian brown frogs allied to Rana dybowskii (Anura, Ranidae).

We investigated the species diversity and phylogeography of the Northeast Asian brown frogs allied to Rana dybowskii (the R. dybowskii species complex: R. dybowskii, R. pirica, and R. uenoi) using four mitochondrial and three nuclear loci. Phylogenetic analyses confirmed the existence of three distinct species in this complex; using extensive molecular data, we confirm the validity of Rana uenoi recognized as a distinct species, and infer R. dybowskii and R. pirica to be sister species. Also, we included populations from previously unsampled regions in Northeast China, and identified them to be R. dybowskii. While many species in Northeast Asia diverged due to Pleistocene glaciation, divergence-dating analyses inferred older, Miocene speciation in the R. dybowskii species complex. Ancestral area reconstruction identified the orogenic movement of the Changbai Mountain Range and the opening of the Sea of Japan/East Sea being major events influencing allopatric speciation.

Influence of geology and human activity on the genetic structure and demography of the Oriental fire-bellied toad (Bombina orientalis).

The Oriental fire-bellied toad (Bombina orientalis) is a commonly used study organism, but knowledge of its evolutionary history is incomplete. We analyze sequence data from four genetic markers (mtDNA genes encoding cytochrome c oxidase subunit I, cytochrome b, and 12S-16S rRNA; nuDNA gene encoding recombination activating gene 2) from 188 individuals across its range in Northeast Asia to elucidate phylogeographic patterns and to identify the historic events that shaped its evolutionary history. Although morphologically similar across its range, B. orientalis exhibits phylogeographic structure, which we infer was shaped by geologic, climatic, and anthropogenic events. Phylogenetic and divergence-dating analyses recover four genetically distinct groups of B. orientalis: Lineage 1-Shandong Province and Beijing (China); Lineage 2-Bukhan Mountain (Korea); Lineage 3-Russia, Northeast China, and northern South Korea; and Lineage 4-South Korea. Lineage 2 was previously unknown. Additionally, we discover an area of secondary contact on the Korean Peninsula, and infer a single dispersal event as the origin of the insular Jeju population. Skyline plots estimate different population histories for the four lineages: Lineages 1 and 2 experienced population decreases, Lineage 3 remained stable, while Lineage 4 experienced a sharp increase during the Holocene. The timing of the population expansion of Lineage 4 coincides with the advent of rice cultivation, which may have facilitated the increase in population size by providing additional breeding habitat.

Penicillium jejuense sp. nov., isolated from the marine environments of Jeju Island, Korea.

Three strains of an unidentified Penicillium species were isolated during a fungal diversity survey of marine environments in Korea. These strains are described here as a new species following a multigene phylogenetic analyses of nuc rDNA internal transcribed spacer barcodes (ITS1-5.8S-ITS2), genes for ß-tubulin, calmodulin and RNA polymerase II second largest subunit, and observation of macro-and micromorphological characters. Phylogenetic analyses revealed that the three strains formed a strongly supported monophyletic group distinct from previously reported species of section Aspergilloides. Morphologically this species can be distinguished from its sister species, P. crocicola, by the reverse color on Czapek yeast autolysate agar, abundant production of sclerotia on malt extract agar and colony characters on yeast extract sucrose agar. We name this new species P. jejuense, after the locality where it was discovered. At 25 C for 7 d, P. jejuense colonies grew to 55-60 mm on CYA, 45-48 mm on MEA, 48-52 mm on YES and 23-26 mm on CREA. Conidia (2.2-3.4 × 2.0-2.6 µm) and sclerotia (160-340 × 125-210 µm) were globose to ellipsoidal.

Trichoderma songyi sp. nov., a new species associated with the pine mushroom (Tricholoma matsutake).

A new species, Trichoderma songyi, was found to be associated with the pine mushroom (Tricholoma matsutake) in Korea. This species was isolated from three different substrates: Tricholoma matsutake basidiomata, as well as roots of Pinus densiflora and soil in the fairy ring. Based on its molecular and phenotypic characteristics, we demonstrate that Trichoderma songyi is unique and distinguishable from closely related species. We performed phylogenetic analyses based on two molecular markers, the genes for both translation elongation factor 1-alpha and the second largest subunit of RNA polymerase II. Phylogenetic analyses showed that Trichoderma songyi is closely related to Trichoderma koningii aggregate and Trichoderma caerulescens. Morphologically, Trichoderma songyi can be distinguished from these closely related taxa by its growth rates, colony morphology on PDA in darkness, and coconut-like odour. Due to the economic importance of the pine mushroom, the relationship between Trichoderma songyi and Tricholoma matsutake should be studied further.

Marine-derived Penicillium in Korea: diversity, enzyme activity, and antifungal properties.

The diversity of marine-derived Penicillium from Korea was investigated using morphological and multigene phylogenetic approaches, analyzing sequences of the internal transcribed spacer region, ß-tubulin gene, and RNA polymerase subunit II gene. In addition, the biological activity of all isolated strains was evaluated. We tested for the extracellular enzyme activity of alginase, endoglucanase, and ß-glucosidase, and antifungal activity against two plant pathogens (Colletotrichum acutatum and Fusarium oxysporum). A total of 184 strains of 36 Penicillium species were isolated, with 27 species being identified. The most common species were Penicillium polonicum (19.6 %), P. rubens (11.4 %), P. chrysogenum (11.4 %), and P. crustosum (10.9 %). The diversity of Penicillium strains isolated from soil (foreshore soil and sand) and marine macroorganisms was higher than the diversity of strains isolated from seawater. While many of the isolated strains showed alginase and ß-glucosidase activity, no endoglucanase activity was found. More than half the strains (50.5 %) showed antifungal activity against at least one of the plant pathogens tested. Compared with other strains in this study, P. citrinum (strain SFC20140101-M662) showed high antifungal activity against both plant pathogens. The results reported here expand our knowledge of marine-derived Penicillium diversity. The relatively high proportion of strains that showed antifungal and enzyme activity demonstrates that marine-derived Penicillium have great potential to be used in the production of natural bioactive products for pharmaceutical and/or industrial use.

Are fecal samples an appropriate proxy for amphibian intestinal microbiota?

The intestinal microbiota, an invisible organ supporting a host's survival, has essential roles in metabolism, immunity, growth, and development. Since intestinal microbiota influences a host's biology, application of such data to wildlife conservation has gained interest. There are standard protocols for studying the human intestinal microbiota, but no equivalent for wildlife. A major challenge is sampling the intestinal microbiota in an effective, unbiased way. Fecal samples are a popular proxy for intestinal microbiota because collection is non-invasive and allows for longitudinal sampling. Yet it is unclear whether the fecal microbiota is representative of the intestinal microbiota. In wildlife studies, research on the sampling methodology is limited. In this study focusing on amphibians, we characterize and compare the microbiota (small intestine, large intestine, and feces) of two Hong Kong stream-dwelling frog species: Lesser Spiny Frog (Quasipaa exilispinosa) and Hong Kong Cascade Frog (Amolops hongkongensis). We found that the microbiota of both species are similar at the level of phylum and family, but diverge at the level of genus. When we assessed the performance of fecal microbiota in representing the intestinal microbiota in these two species, we found that (1) the microbiota of the small and large intestine differs significantly, (2) feces are not an appropriate proxy of either intestinal sections, and (3) a set of microbial taxa significantly differs between sample types. Our findings raise caution equating fecal and intestinal microbiota in stream-dwelling frogs. Sampling feces can avoid sacrifice of an animal, but researchers should avoid over-extrapolation and interpret results carefully.

Genetic evidence for a high diversity and wide distribution of endemic strains of the pathogenic chytrid fungus Batrachochytrium dendrobatidis in wild Asian amphibians.

Population declines and extinctions of amphibians have been attributed to the chytrid fungus Batrachochytrium dendrobatidis (Bd), especially one globally emerging recombinant lineage ('Bd-GPL'). We used PCR assays that target the ribosomal internal transcribed spacer region (ITS) of Bd to determine the prevalence and genetic diversity of Bd in South Korea, where Bd is widely distributed but is not known to cause morbidity or mortality in wild populations. We isolated Korean Bd strains from native amphibians with low infection loads and compared them to known worldwide Bd strains using 19 polymorphic SNP and microsatellite loci. Bd prevalence ranged between 12.5 and 48.0%, in 11 of 17 native Korean species, and 24.7% in the introduced bullfrog Lithobates catesbeianus. Based on ITS sequence variation, 47 of the 50 identified Korean haplotypes formed a group closely associated with a native Brazilian Bd lineage, separated from the Bd-GPL lineage. However, multilocus genotyping of three Korean Bd isolates revealed strong divergence from both Bd-GPL and the native Brazilian Bd lineages. Thus, the ITS region resolves genotypes that diverge from Bd-GPL but otherwise generates ambiguous phylogenies. Our results point to the presence of highly diversified endemic strains of Bd across Asian amphibian species. The rarity of Bd-GPL-associated haplotypes suggests that either this lineage was introduced into Korea only recently or Bd-GPL has been outcompeted by native Bd strains. Our results highlight the need to consider possible complex interactions among native Bd lineages, Bd-GPL and their associated amphibian hosts when assessing the spread and impact of Bd-GPL on worldwide amphibian populations.

Re-description of the genus Heterobiantes and its type species H. geniculatus (Opiliones: Laniatores: Epedanidae) with additional genital and external morphology data.

The genus Heterobiantes Roewer, 1912 and its type species Epedanus geniculatus Pocock, 1903 are redescribed based on the type specimen currently housed within the British Museum of Natural History (BMNH), London, along with new specimens captured in Hong Kong, providing updated morphological data and male genital description. Closer inspections reveal how the species is different from other similar genera within Epedaninae and rationalizes its placement as a monotypic genus within the sub-family. Genetic data are also included.

Revisiting the genetic diversity and population structure of the endangered Green Sea Turtle (Chelonia mydas) breeding populations in the Xisha (Paracel) Islands, South China Sea.

The Green Sea Turtle (Chelonia mydas) is an umbrella species in the South China Sea, a Chinese national first-level protected wild animal, and the only sea turtle that nests in waters around China. The largest C. mydas nesting ground is distributed in the Xisha (Paracel) Islands, which plays a vital role in the survival of sea turtle populations in the region. This study reveals the genetic diversity and population structure of the breeding population of C. mydas in the Xisha (Paracel) Islands using three mitochondrial markers. A total of 15 D-loop, five Cytochrome b (Cyt b), and seven Cytochrome C Oxidase subunit I (COI) haplotypes were identified in the breeding population of C. mydas in the Xisha (Paracel) Islands. D-loop haplotypes are distributed in clades III, IV, and VIII of the C. mydas mitochondrial control region. It is the first time that one haplotype from Clade IV was found in this C. mydas population, and five new D-loop haplotypes were also identified. The haplotype and nucleotide diversity were calculated for each marker: D-loop (0.415 haplotype diversity, 0.00204 nucleotide diversity), Cyt b (0.140, 0.00038) and COI (0.308, 0.00083). The average genetic distance (p) of each molecular marker was less than 0.01. Neutral detection and nucleotide mismatch analysis suggested that the breeding population of C. mydas in the Xisha (Paracel) Islands did not experience a population expansion event in recent history. It is recommended that a sea turtle protection area be established in the Xisha (Paracel) Islands area to strengthen protection and effectively protect the uniqueness and sustainability of the breeding population of C. mydas in the South China Sea.

Using eDNA techniques to find the endangered big-headed turtle (Platysternon megacephalum).

We evaluate the efficacy of environmental DNA (eDNA) techniques to locate wild populations and estimate the population size of the endangered big-headed turtle (Platysternon megacephalum) in Hong Kong. The results from this study are important for identifying priority sites for protection and further research. Additionally, we assess the impact of two environmental variables (temperature and pH) on eDNA quantity. We surveyed 34 streams for three years, sampling four times each year. Four new populations were first identified with eDNA analysis, and then verified by field surveys. Our multi-year survey highlights that eDNA detection can be inconsistent over time, even in streams with known populations. There was no significant relationship between eDNA quantity and the environmental variables tested. Lastly, our results suggest that eDNA methods remain promising to estimate population size, since number of positive detections were positively correlated with population size in streams with known populations. We conclude that eDNA methods are powerful, but care must be taken when interpreting field results as they are affected by species ecology and environmental conditions.

Evaluating phylogenetic informativeness and data-type usage for new protein-coding genes across Vertebrata.

As a resource for vertebrate phylogenetics, we developed 75 new protein-coding genes using a combination of expressed sequence tags (ESTs) available in Genbank, and targeted amplification of complementary DNA (cDNA). In addition, we performed three additional analyses in order to assess the utility of our approach. First, we profiled the phylogenetic informativeness of these new markers using the online program PhyDesign. Next, we compared the utility of four different data-types used in phylogenetics: nucleotides (NUCL), amino acids (AA), 1st and 2nd codon positions only (N12), and modified sequences to account for codon degeneracy (DEGEN1; Regier et al., 2010). Lastly, we use these new markers to construct a vertebrate phylogeny and address the uncertain relationship between higher-level mammal groups: monotremes, marsupials, and placentals. Our results show that phylogenetic informativeness of the 75 new markers varies, both in the amount of phylogenetic signal and optimal timescale. When comparing the four data-types, we find that the NUCL data-type, due to the high level of phylogenetic signal, performs the best across all divergence times. The remaining three data-types (AA, N12, DEGEN1) are less subject to homoplasy, but have greatly reduced levels of phylogenetic signal relative to NUCL. Our phylogenetic inference supports the Theria hypothesis of mammalian relationships, with marsupials and placentals being sister groups.

Taxonomic Revision of the Genus Lactifluus (Russulales, Basidiomycota) of South Korea.

Lactifluus (Pers.) Roussel is an ectomycorrhizal genus that was recently recognized to be distinct from the genus Lactarius. To date, 226 Lactifluus species have been reported worldwide. Misidentification of Lactifluus species is common because of intraspecific morphological variation, cryptic diversity, and the limited number of taxonomic keys available. Molecular data are indispensable for species delimitation; a multilocus phylogenetic analysis showed that most Asian Lactifluus species are not conspecific with morphologically similar species present on other continents. In particular, Korea has misused European and North American Lactifluus names. In this study, we evaluated the taxonomy of Lactifluus in Korea using both morphological and multilocus molecular (ITS, nrLSU, rpb1, and rpb2) data. We examined 199 Lactifluus specimens collected between 1980 and 2016, and a total of 24 species across the four Lactifluus subgenera were identified. All Korean species are distinct and clearly separated from European and North American species. Five taxa corresponded to previously described species from Asia and the remaining 19 taxa are confirmed as new species. Herein, we provide keys to the Korean Lactifluus species within their subgenera, molecular phylogenies, a summary of diversity, and detailed description of the new species.

Phylogeographic study of the Bufo gargarizans species complex, with emphasis on Northeast Asia.

We conduct a phylogeographic and population genetic study of the Asiatic toad (Bufo gargarizans) to understand its evolutionary history, and the influence of geology and climate. A total of 292 individuals from 94 locations were genotyped for two mitochondrial loci (cytb, ND2) and five nuclear introns (Sox9-2, Rho-3, CCNB2-3, UCH-2, and DBI-2), and we performed a suite of phylogenetic, population genetic, and divergence dating analyses. The phylogenetic trees constructed using mitochondrial loci inferred B. gargarizans being divided into two major groups: China mainland and Northeast Asia (Northeast China, Russia, and Korean Peninsula). As with previous studies of this species, we recover population genetic structure not tied to geographic region. Additionally, we discover a new genetic clade restricted to Northeast Asia that points towards the Korean Peninsula being a glacial refugium during the Pleistocene. The weak phylogeographic pattern of B. gargarizans is likely the result of multiple biological, anthropogenic, and historical factors - robust dispersal abilities as a consequence of physiological adaptations, human translocation, geologic activity, and glacial cycles of the Pleistocene. We highlight the complex geologic and climatic history of Northeast Asia and encourage further research to understand its impact on the biodiversity in the region.

Comparative Analysis of the Fecal Microbiota of Wild and Captive Beal's Eyed Turtle (Sacalia bealei) by 16S rRNA Gene Sequencing.

The Beal's eyed turtle (Sacalia bealei) is threatened with extinction due to hunting for large-scale trade. In Hong Kong, there are some of the world's remaining wild populations of S. bealei, as well as a breeding colony. This breeding colony is at the core of conservation efforts (captive breeding, reintroduction programs). Therefore, we would like to know how captivity, in particular diet, affects the gut microbiota. Using high-throughput 16S rRNA gene sequencing, we comparatively analyzed the fecal microbiota of wild and captive S. bealei. We found that wild S. bealei have higher alpha diversity than captive S. bealei, but the difference was not significant. Significant differences were found in ß-diversity; at the phylum level, wild S. bealei have higher relative abundances of Proteobacteria and captive S. bealei have higher relative abundances of Firmicutes. At the genus level, Cetobacterium and Citrobacter are more abundant in wild S. bealei, while Clostridium spp. are significantly more abundant in captive S. bealei. These results suggest conditions in captivity, with diet being a major factor, influence the gut microbiota of S. bealei. The connection between diet and health has always been considered for captive animals, and in this study we use the gut microbiota as an another tool to assess health.

Large-Scale Hybridisation as an Extinction Threat to the Suweon Treefrog (Hylidae: Dryophytes suweonensis).

Amphibians are in the midst of a sixth mass extinction, and human activities play a major role in pushing species towards extinction. Landscape anthropisation has impacts that indirectly threaten species, in addition to the obvious destruction of natural habitats. For instance, land modification may bring human-commensal species in contact with sister-clades from which they were previously isolated. The species in these new contact zones are then able to hybridise to the point of reaching lineage fusion, through which the gene pool of the two species merges and one of the parental lineages becomes extirpated. Here, we documented the patterns of hybridisation between the spatially restricted D. suweonensis and the widespread D. japonicus. On the basis of the analysis of Cytochrome c oxidase subunit I mitochondrial DNA sequences (404 individuals from 35 sites) and six polymorphic microsatellites (381 individuals from 34 sites), we revealed a generalised, bi-directional, and geographically widespread hybridisation between the two species. Evidence of fertile back-crosses is provided by relatively high numbers of individuals in cyto-nuclear disequilibrium, as well as the presence of hybrid individuals further south than the species distribution limit, determined on the basis of call properties. Hybridisation is an additional threat to the endangered D. suweonensis.

Integrating and updating wildlife conservation in China.

China has about 11% of the world's total wildlife species, so strengthening China's wildlife conservation is of great significance to global biodiversity. Despite some successful cases and conservation efforts, 21.4% of China's vertebrate species are threatened by human activities. The booming wildlife trade in China has posed serious threat to wildlife in China and throughout the world, while leading to a high risk of transmission of infectious zoonotic diseases. China's wildlife conservation has faced a series of challenges, two of which are an impractical, separated management of wildlife and outdated protected species lists. Although the Wildlife Protection Law of China was revised in 2016, the issues of separated management remain, and the protected species lists are still not adequately revised. These issues have led to inefficient and overlapping management, waste of administrative resources, and serious obstacles to wildlife protection. In this article, we analyze the negative effects of current separated management of wildlife species and outdated protected species lists, and provide some suggestions for amendment of the laws and reform of wildlife management system.

Influence of Season and Soil Properties on Fungal Communities of Neighboring Climax Forests (Carpinus cordata and Fraxinus rhynchophylla).

Trees in forest ecosystems constantly interact with the soil fungal community, and this interaction plays a key role in nutrient cycling. The diversity of soil fungal communities is affected by both environmental factors and host tree species. We investigated the influence of both of these factors by examining the total fungal communities in the rhizospheric soil of climax tree species that have similar ecological roles (Carpinus cordata, an ectomycorrhizal [ECM] tree, and Fraxinus rhynchophylla, an arbuscular mycorrhizal [AM] tree) in temperate forests with continental climates of Mt. Jeombong, South Korea. Fungal communities were assessed by Illumina-MiSeq sequencing the internal transcribed spacer (ITS) region of environmental DNA, and comparing their environmental factors (season and soil properties). We found that soil fungi of the two forest types differed in terms of community structure and ecological guild composition. The total fungal community composition changed significantly with seasons and soil properties in the F. rhynchophylla forest, but not in the C. cordata forest. However, potassium and carbon were significantly correlated with fungal diversity in both forests, and a positive correlation was found only between symbiotrophs of C. cordata and the carbon to nitrogen (C/N) ratio. Thus, the effects of environmental factors on soil fungal communities depended on the host trees, but some factors were common in both forests. Our results indicate that individual tree species should be considered when anticipating how the fungal communities will respond to environmental change.

Successional Change of the Fungal Microbiome Pine Seedling Roots Inoculated With Tricholoma matsutake.

The pine mushroom (Tricholoma matsutake; Agaricales, Tricholomataceae) is an ectomycorrhizal fungus that produces a commercially valuable, edible mushrooms. Attempts to artificially cultivate T. matsutake has so far been unsuccessful. One method used to induce T. matsutake to produce fruiting bodies of in the wild is shiro (mycelial aggregations of T. matsutake) transplantation. In vitro ectomycorrhization of T. matsutake with seedlings of Pinus densiflora has been successful, but field trials showed limited production of fruiting bodies. Few studies have been done to test what happens after transplantation in the wild, whether T. matsutake persists on the pine seedling roots or gets replaced by other fungi. Here, we investigated the composition and the interaction of the root fungal microbiome of P. densiflora seedlings inoculated with T. matsutake over a 3 year period after field transplantation, using high-throughput sequencing. We found a decline of T. matsutake colonization on pine roots and succession of mycorrhizal fungi as P. densiflora seedlings grew. Early on, roots were colonized by fast-growing, saprotrophic Ascomycota, then later replaced by early stage ectomycorrhiza such as Wilcoxina. At the end, more competitive Suillus species dominated the host roots. Most of the major OTUs had negative or neutral correlation with T. matsutake, but several saprotrophic/plant pathogenic/mycoparasitic species in genera Fusarium, Oidiodendron, and Trichoderma had positive correlation with T. matsutake. Four keystone species were identified during succession; two species (Fusarium oxysporum, and F. trincintum) had a positive correlation with T. matsutake, while the other two had a negative correlation (Suillus granulatus, Cylindrocarpon pauciseptatum). These findings have important implications for further studies on the artificial cultivation of T. matsutake.