Molecular phylogeny of Zygomycota based on EF-1alpha and RPB1 sequences: limitations and utility of alternative markers to rDNA

Earlier molecular phylogenetic analyses based on nuclear small subunit ribosomal DNA (nSSU rDNA) suggest that the Zygomycota are polyphyletic within the Chytridiomycota. However, these analyses failed to resolve almost all interordinal relationships among basal fungi (Chytridiomycota and Zygomycota), due to lack of sufficient characters within the nSSU rDNA. To further elucidate the higher-level phylogeny of Zygomycota, we have sequenced partial RPB1 (DNA dependent RNA polymerase II largest subunit) and EF-1alpha (translation elongation factor 1 alpha) genes from 10 and 3 zygomycete fungi, respectively. Independent molecular phylogenetic analyses were performed based on each sequence by distance and maximum likelihood methods. Although deep phylogenetic relationships among basal fungi still remain poorly resolved using either gene, the RPB1-based phylogeny identified a novel monophyletic clade consisting of the Dimargaritales, Harpellales, and Kickxellales. This result suggests that regularly formed septa (cross walls that divide hyphae into segments) with a lenticular cavity are plesiomorphic for this clade, and indicates the importance of septal pore ultrastructure in zygomycete phylogeny. In addition, a peculiar mucoralean genus Mortierella, which was considered to be distantly related to the other Mucorales based on previous nSSU rDNA analyses, was resolved as the basal most divergence within the Mucorales, consistent with traditional phenotypic-based taxonomy. Although the taxa included in our analysis are restricted, the monophyly of each order suggested by nSSU rDNA phylogeny is supported by the present RPB1-based analysis. These results support the potential use of RPB1 as an alternative marker for fungal phylogenetic studies. Conversely, the overall fungal phylogeny based on EF-1alpha sequence is poorly resolved. A comparison of numbers of observed substitutions versus inferred substitutions within EF-1alpha indicates that this gene is much more saturated than RPB1. This result suggests that the EF-1alpha gene is unsuitable for resolving higher-level phylogenetic relationships within the Fungi.     

DNA分子标记技术在真菌系统学研究中的应用及影响

DNA分子标记技术为真菌系统进化研究提供了许多新的方法,真菌分子系统学已成为一门成熟的学科。简述了真菌分子系统学的发展简史和代表性的研究方法以及对真菌系统学的主要贡献,包括将广义的真菌划分为3个类群,粘菌和卵菌不再属于真菌界成员。真菌生命之树项目的研究结果对真菌界高阶分类系统作出重大调整,将先前的4个门(壶菌门、接合菌门、子囊菌门和担子菌门)变为7个门(微孢子虫门、壶菌门、新丽鞭毛菌门、芽枝霉门、球囊菌门、子囊菌门和担子菌门)和4个亚门,并对真菌各类群概念作出修订。此外,DNA分子标记技术对真菌种概念的认识、有性型-无性型关联及分子生态学等研究领域产生了重要影响。     

Calculating the evolutionary rates of different genes: a fast, accurate estimator with applications to maximum likelihood phylogenetic analysis

In phylogenetic analyses with combined multigene or multiprotein data sets, accounting for differing evolutionary dynamics at different loci is essential for accurate tree prediction. Existing maximum likelihood (ML) and Bayesian approaches are computationally intensive. We present an alternative approach that is orders of magnitude faster. The method, Distance Rates (DistR), estimates rates based upon distances derived from gene/protein sequence data. Simulation studies indicate that this technique is accurate compared with other methods and robust to missing sequence data. The DistR method was applied to a fungal mitochondrial data set, and the rate estimates compared well to those obtained using existing ML and Bayesian approaches. Inclusion of the protein rates estimated from the DistR method into the ML calculation of trees as a branch length multiplier resulted in a significantly improved fit as measured by the Akaike Information Criterion (AIC). Furthermore, bootstrap support for the ML topology was significantly greater when protein rates were used, and some evident errors in the concatenated ML tree topology (i.e., without protein rates) were corrected. [Bayesian credible intervals; DistR method; multigene phylogeny; PHYML; rate heterogeneity.].     

ImOSM: intermittent evolution and robustness of phylogenetic methods

Among the criteria to evaluate the performance of a phylogenetic method, robustness to model violation is of particular practical importance as complete a priori knowledge of evolutionary processes is typically unavailable. For studies of robustness in phylogenetic inference, a utility to add well-defined model violations to the simulated data would be helpful. We therefore introduce ImOSM, a tool to imbed intermittent evolution as model violation into an alignment. Intermittent evolution refers to extra substitutions occurring randomly on branches of a tree, thus changing alignment site patterns. This means that the extra substitutions are placed on the tree after the typical process of sequence evolution is completed. We then study the robustness of widely used phylogenetic methods: maximum likelihood (ML), maximum parsimony (MP), and a distance-based method (BIONJ) to various scenarios of model violation. Violation of rates across sites (RaS) heterogeneity and simultaneous violation of RaS and the transition/transversion ratio on two nonadjacent external branches hinder all the methods recovery of the true topology for a four-taxon tree. For an eight-taxon balanced tree, the violations cause each of the three methods to infer a different topology. Both ML and MP fail, whereas BIONJ, which calculates the distances based on the ML estimated parameters, reconstructs the true tree. Finally, we report that a test of model homogeneity and goodness of fit tests have enough power to detect such model violations. The outcome of the tests can help to actually gain confidence in the inferred trees. Therefore, we recommend using these tests in practical phylogenetic analyses.     

Taxon-specific fungal primers reveal unexpectedly high diversity during leaf decomposition in a stream

Traditional techniques for studying the fungal community composition in streams favour the detection and identification of aquatic hyphomycetes. Our objective was to use molecular techniques to determine the presence and contributions of other fungal groups. We designed primers specific for the ITS regions in Ascomycota, Basidiomycota, Chytridiomycota, Zygomycota and Oomycota. The primers were used to amplify DNA from linden, maple, and beech leaves, and birch wood submerged in a stream for 4 weeks in summer, autumn, winter and spring. The amplification products were separated by denaturing gradient gel electrophoresis. Ascomycota were present in large phylotype numbers (up to 21) on all substrates and all dates and represented ≥ 75 % of the fungal biomass. Basidiomycota were the second most abundant group in summer and autumn (up to 13 % on wood) and were absent only on linden and maple in spring. There were consistently large numbers of phylotypes from Chytridiomycota and their relative contribution to the microbial community peaked in winter on all substrates. Oomycota were present in summer and abundant only on wood. Zygomycota were present in low numbers and their estimated contribution to fungal biomass was ≤ 1%. Using primers to target individual groups facilitates a more balanced approach to studying fungal diversity in freshwater ecosystems.     

Loss of the flagellum happened only once in the fungal lineage: phylogenetic structure of Kingdom Fungi inferred from RNA polymerase II subunit genes

Background  

At present, there is not a widely accepted consensus view regarding the phylogenetic structure of kingdom Fungi although two major phyla, Ascomycota and Basidiomycota, are clearly delineated. Regarding the lower fungi, Zygomycota and Chytridiomycota, a variety of proposals have been advanced. Microsporidia may or may not be fungi; the Glomales (vesicular-arbuscular mycorrhizal fungi) may or may not constitute a fifth fungal phylum, and the loss of the flagellum may have occurred either once or multiple times during fungal evolution. All of these issues are capable of being resolved by a molecular phylogenetic analysis which achieves strong statistical support for major branches. To date, no fungal phylogeny based upon molecular characters has satisfied this criterion.     

Exploring among-site rate variation models in a maximum likelihood framework using empirical data: effects of model assumptions on estimates of topology,branch lengths,and bootstrap support

We have investigated the effects of different among-site rate variation models on the estimation of substitution model parameters, branch lengths, topology, and bootstrap proportions under minimum evolution (ME) and maximum likelihood (ML). Specifically, we examined equal rates, invariable sites, gamma-distributed rates, and site-specific rates (SSR) models, using mitochondrial DNA sequence data from three protein-coding genes and one tRNA gene from species of the New Zealand cicada genus Maoricicada. Estimates of topology were relatively insensitive to the substitution model used; however, estimates of bootstrap support, branch lengths, and R-matrices (underlying relative substitution rate matrix) were strongly influenced by the assumptions of the substitution model. We identified one situation where ME and ML tree building became inaccurate when implemented with an inappropriate among-site rate variation model. Despite the fact the SSR models often have a better fit to the data than do invariable sites and gamma rates models, SSR models have some serious weaknesses. First, SSR rate parameters are not comparable across data sets, unlike the proportion of invariable sites or the alpha shape parameter of the gamma distribution. Second, the extreme among-site rate variation within codon positions is problematic for SSR models, which explicitly assume rate homogeneity within each rate class. Third, the SSR models appear to give severe underestimates of R-matrices and branch lengths relative to invariable sites and gamma rates models in this example. We recommend performing phylogenetic analyses under a range of substitution models to test the effects of model assumptions not only on estimates of topology but also on estimates of branch length and nodal support.     

A structural EM algorithm for phylogenetic inference.

A central task in the study of molecular evolution is the reconstruction of a phylogenetic tree from sequences of current-day taxa. The most established approach to tree reconstruction is maximum likelihood (ML) analysis. Unfortunately, searching for the maximum likelihood phylogenetic tree is computationally prohibitive for large data sets. In this paper, we describe a new algorithm that uses Structural Expectation Maximization (EM) for learning maximum likelihood phylogenetic trees. This algorithm is similar to the standard EM method for edge-length estimation, except that during iterations of the Structural EM algorithm the topology is improved as well as the edge length. Our algorithm performs iterations of two steps. In the E-step, we use the current tree topology and edge lengths to compute expected sufficient statistics, which summarize the data. In the M-Step, we search for a topology that maximizes the likelihood with respect to these expected sufficient statistics. We show that searching for better topologies inside the M-step can be done efficiently, as opposed to standard methods for topology search. We prove that each iteration of this procedure increases the likelihood of the topology, and thus the procedure must converge. This convergence point, however, can be a suboptimal one. To escape from such "local optima," we further enhance our basic EM procedure by incorporating moves in the flavor of simulated annealing. We evaluate these new algorithms on both synthetic and real sequence data and show that for protein sequences even our basic algorithm finds more plausible trees than existing methods for searching maximum likelihood phylogenies. Furthermore, our algorithms are dramatically faster than such methods, enabling, for the first time, phylogenetic analysis of large protein data sets in the maximum likelihood framework.     

贝加尔针茅草原土壤真菌群落结构对氮素和水分添加的响应

研究氮沉降和降雨变化对土壤真菌群落结构的互作效应,对未来预测多个气候变化因子对草地生态系统的交互作用具有重要意义。以施氮和灌溉模拟氮沉降和降雨增加,采用裂区设计,应用高通量测序技术,研究8个氮添加水平(0、15、30、50、100、150、200、300 kg N hm~(-2)a~(-1))和2个水分添加水平(不灌溉、模拟夏季增雨100 mm灌溉)对土壤真菌群落结构的影响。结果表明,氮素和水分添加后,土壤真菌群落中占优势的门类分别为接合菌门Zygomycota(22.0%—48.9%)、担子菌门Basidiomycota(7.8%—18.5%)、子囊菌门Ascomycota(9.4%—20.1%)、球囊菌门Glomeromycota(0.7%—3.1%)、壶菌门Chytridiomycota(0.1%—1.3%)。常规降雨条件下,随着氮添加水平升高,接合菌门相对丰度呈现出先升高后降低的趋势,N50处理最高;子囊菌门相对丰度在高氮添加时(N100—N300)呈升高趋势。而在氮素和水分同时添加条件下,随着氮添加水平升高,接合菌门相对丰度呈降低趋势,子囊菌门相对丰度变化则不明显。在相同的氮添加水平下,水分添加使接合菌门相对丰度增加,而担子菌门、子囊菌门、球囊菌门和壶菌门的相对丰度降低。在不同氮素和水分添加条件下,有5个土壤真菌门类11个真菌纲相对丰度变化显著。接合菌门的Mortierella属,担子菌门的Entolomataceae科和Geastrum属相对丰度变化极显著,可作为土壤真菌群落结构变化的指示种。PCo A分析结果也表明氮素和水分添加改变了土壤真菌群落结构。植物-土壤-微生物系统的结构方程模型结果表明,植物群落组成及植物物种丰富度的变化是土壤真菌群落结构发生变化的主要影响因素,土壤无机氮及p H的变化主要通过影响植物群落间接影响真菌群落,其对真菌群落的直接影响则较小。综上,氮素和水分添加改变了土壤真菌群落结构,且两者存在明显的互作效应,水分添加可改变氮添加对土壤真菌群落的影响。     

Combined nuclear and mitochondrial DNA sequences resolve generic relationships within the Cracidae (Galliformes,Aves)

The Cracidae is one of the most endangered and distinctive bird families in the Neotropics, yet the higher relationships among taxa remain uncertain. The molecular phylogeny of its 11 genera was inferred using 10,678 analyzable sites (5,412 from seven different mitochondrial segments and 5,266 sites from four nuclear genes). We performed combinability tests to check conflicts in phylogenetic signals of separate genes and genomes. Phylogenetic analysis showed that the unrooted tree of ((curassows, horned guan) (guans, chachalacas)) was favored by most data partitions and that different data partitions provided support for different parts of the tree. In particular, the concatenated mitochondrial DNA (mtDNA) genes resolved shallower nodes, whereas the combined nuclear sequences resolved the basal connections among the major clades of curassows, horned guan, chachalacas, and guans. Therefore, we decided that for the Cracidae all data should be combined for phylogenetic analysis. Maximum parsimony (MP), maximum likelihood (ML), and Bayesian analyses of this large data set produced similar trees. The MP tree indicated that guans are the sister group to (horned guan, (curassows, chachalacas)), whereas the ML and Bayesian analysis recovered a tree where the horned guan is a sister clade to curassows, and these two clades had the chachalacas as a sister group. Parametric bootstrapping showed that alternative trees previously proposed for the cracid genera are significantly less likely than our estimate of their relationships. A likelihood ratio test of the hypothesis of a molecular clock for cracid mtDNA sequences using the optimal ML topology did not reject rate constancy of substitutions through time. We estimated cracids to have originated between 64 and 90 million years ago (MYA), with a mean estimate of 76 MYA. Diversification of the genera occurred approximately 41-3 MYA, corresponding with periods of global climate change and other Earth history events that likely promoted divergences of higher level taxa.     

Does choice in model selection affect maximum likelihood analysis?

In order to have confidence in model-based phylogenetic analysis, the model of nucleotide substitution adopted must be selected in a statistically rigorous manner. Several model-selection methods are applicable to maximum likelihood (ML) analysis, including the hierarchical likelihood-ratio test (hLRT), Akaike information criterion (AIC), Bayesian information criterion (BIC), and decision theory (DT), but their performance relative to empirical data has not been investigated thoroughly. In this study, we use 250 phylogenetic data sets obtained from TreeBASE to examine the effects that choice in model selection has on ML estimation of phylogeny, with an emphasis on optimal topology, bootstrap support, and hypothesis testing. We show that the use of different methods leads to the selection of two or more models for approximately 80% of the data sets and that the AIC typically selects more complex models than alternative approaches. Although ML estimation with different best-fit models results in incongruent tree topologies approximately 50% of the time, these differences are primarily attributable to alternative resolutions of poorly supported nodes. Furthermore, topologies and bootstrap values estimated with ML using alternative statistically supported models are more similar to each other than to topologies and bootstrap values estimated with ML under the Kimura two-parameter (K2P) model or maximum parsimony (MP). In addition, Swofford-Olsen-Waddell-Hillis (SOWH) tests indicate that ML trees estimated with alternative best-fit models are usually not significantly different from each other when evaluated with the same model. However, ML trees estimated with statistically supported models are often significantly suboptimal to ML trees made with the K2P model when both are evaluated with K2P, indicating that not all models perform in an equivalent manner. Nevertheless, the use of alternative statistically supported models generally does not affect tests of monophyletic relationships under either the Shimodaira-Hasegawa (S-H) or SOWH methods. Our results suggest that although choice in model selection has a strong impact on optimal tree topology, it rarely affects evolutionary inferences drawn from the data because differences are mainly confined to poorly supported nodes. Moreover, since ML with alternative best-fit models tends to produce more similar estimates of phylogeny than ML under the K2P model or MP, the use of any statistically based model-selection method is vastly preferable to forgoing the model-selection process altogether.     

Property and efficiency of the maximum likelihood method for molecular phylogeny

Summary The maximum likelihood (ML) method for constructing phylogenetic trees (both rooted and unrooted trees) from DNA sequence data was studied. Although there is some theoretical problem in the comparison of ML values conditional for each topology, it is possible to make a heuristic argument to justify the method. Based on this argument, a new algorithm for estimating the ML tree is presented. It is shown that under the assumption of a constant rate of evolution, the ML method and UPGMA always give the same rooted tree for the case of three operational taxonomic units (OTUs). This also seems to hold approximately for the case with four OTUs. When we consider unrooted trees with the assumption of a varying rate of nucleotide substitution, the efficiency of the ML method in obtaining the correct tree is similar to those of the maximum parsimony method and distance methods. The ML method was applied to Brown et al.'s data, and the tree topology obtained was the same as that found by the maximum parsimony method, but it was different from those obtained by distance methods.     

Evolutionary relationships among basal fungi (Chytridiomycota and Zygomycota): Insights from molecular phylogenetics

Evolutionary relationships of the two basal fungal phyla Chytridiomycota and Zygomycota are reviewed in light of recent molecular phylogenetic investigation based on rDNA (nSSU, nLSU rDNA), entire mitochondrial genomes, and nuclear protein coding gene sequences (e.g., EF-1alpha, RPB1). Accumulated molecular evidence strongly suggests that the two basal fungal phyla are not monophyletic. For example, the chytridiomycete order Blastocladiales appears to be closely related to the zygomycete order Entomophthorales. Within the Zygomycota, a monophyletic clade, consisting of the Dimargaritales, Harpellales, and Kickxellales, which is characterized by a shared unique septal ultrastructure, was identified. Moreover, evidence for the exclusion of zygomycete orders Amoebidiales and Eccrinales from the Fungi, and their placement at the Animal-Fungi boundary has been clearly documented. Microsporidia, a group of amitochondriate organisms currently under intensive study, is not supported as derived within the Fungi, but a fungal affinity cannot be ruled out. Taking these molecular phylogenetic studies into account, we proposed a hypothetical evolutionary framework of basal fungi.     

Fungal community structure at pelagic and littoral sites in Lake Biwa determined with high-throughput sequencing

Fungi may play an important role in material cycling in lakes and oceans; however, only limited information is available on fungal community structure, especially in large lakes such as Lake Biwa. In this study, whole fungal communities were determined seasonally and spatially using a high-throughput sequencing technique. Water samples were collected from the epilimnion, 0–20 m depth, with a Van Dorn sampler at a pelagic site and from the surface at a littoral site in the north basin of Lake Biwa. All pelagic depth samples were combined into one sample. Sampling occurred on 24 April, 22 May, 10 July, and 16 September 2015. DNA was extracted from filtered samples. Metabarcoding analysis targeting fungi-specific internal transcribed spacer 2 regions was performed using an Illumina MiSeq platform. Epilimnetic fungal communities showed high diversity, with 479 operational taxonomic units (OTUs). The OTUs included 122 belonging to the phylum Ascomycota, 127 to Basidiomycota, 38 to Zygomycota, 45 to Chytridiomycota, 2 to Glomeromycota, and 145 were unclassified fungi. Fungal community structures varied seasonally and spatially. Few of the fungal OTUs overlapped between seasons and sites, and specific communities of fungi were detected on each sampling occasion. Results indicated that spatio-temporal variations in fungal communities were high and may be influenced by both internal factors and external factors, such as terrestrial inputs.     

Group I Introns from Zygomycota: Evolutionary Implications for the Fungal IC1 Intron Subgroup

The origins of fungal group I introns within nuclear small-subunit (nSSU) rDNA are enigmatic. This is partly because they have never been reported in basal fungal phyla (Zygomycota and Chytridiomycota), which are hypothesized to be ancestral to derived phyla (Ascomycota and Basidiomycota). Here we report group I introns from the nSSU rDNA of two zygomycete fungi, Zoophagus insidians (Zoopagales) and Coemansia mojavensis (Kickxellales). Secondary structure analyses predicted that both introns belong to the IC1 subgroup and that they are distantly related to each other, which is also suggested by different insertion sites. Molecular phylogenetic analyses indicated that the IC1 intron of Z. insidians is closely related to the IC1 intron inserted in the LSU rDNA of the basidiomycete fungus Clavicorona taxophila, which strongly suggests interphylum horizontal transfer. The IC1 intron of C. mojavensis has a low phylogenetic affinity to other fungal IC1 introns inserted into site 943 of nSSU rDNA (relative to E. coli 16S rDNA). It is noteworthy that this intron contains a putative ORF containing a His–Cys box motif in the antisense strand, a hallmark for nuclear-encoded homing endonucleases. Overall, molecular phylogenetic analyses do not support the placement of these two introns in basal fungal IC1 intron lineages. This result leads to the suggestion that fungal IC1 introns might have invaded or been transferred laterally after the divergence of the four major fungal phyla. Received: 8 February 2001 / Accepted: 1 November 2001     

Evolutionary relationships within the fungi: analyses of nuclear small subunit rRNA sequences.

Nucleotide sequences of the small subunit ribosomal RNA (18S) gene were used to investigate evolutionary relationships within the Fungi. The inferred tree topologies are in general agreement with traditional classifications in the following ways: (1) the Chytridiomycota and Zygomycota appear to be basal groups within the Fungi. (2) The Ascomycota and Basidiomycota are a derived monophyletic group. (3) Relationships within the Ascomycota are concordant with traditional orders and divide the hemi- and euascomycetes into distinct lineages. (4) The Basidiomycota is divided between the holobasidiomycetes and phragmobasidiomycetes. Conflicts with traditional classification were limited to weakly supported branches of the tree. Strongly supported relationships were robust to minor changes in alignment, method of analysis, and various weighting schemes. Weighting, either of transversions or by site, did not convincingly improve the status of poorly supported portions of the tree. The rate of variation at particular sites does not appear to be independent of lineage, suggesting that covariation of sites may be an important phenomenon in these genes.     

Microbial Pathogens in the Fungal Kingdom

The fungal kingdom is vast, spanning ~1.5 to as many as 5 million species diverse as unicellular yeasts, filamentous fungi, mushrooms, lichens, and both plant and animal pathogens. The fungi are closely aligned with animals in one of the six to eight supergroups of eukaryotes, the opisthokonts. The animal and fungal kingdoms last shared a common ancestor ~1 billion years ago, more recently than other groups of eukaryotes. As a consequence of their close evolutionary history and shared cellular machinery with metazoans, fungi are exceptional models for mammalian biology, but prove more difficult to treat in infected animals. The last common ancestor to the fungal/metazoan lineages is thought to have been unicellular, aquatic, and motile with a posterior flagellum, and certain extant species closely resemble this hypothesized ancestor. Species within the fungal kingdom were traditionally assigned to four phyla, including the basal fungi (Chytridiomycota, Zygomycota) and the more recently derived monophyletic lineage, the dikarya (Ascomycota, Basidiomycota). The fungal tree of life project has revealed that the basal lineages are polyphyletic, and thus there are as many as eight to ten fungal phyla. Fungi that infect vertebrates are found in all of the major lineages, and virulence arose multiple times independently. A sobering recent development involves the species Batrachochytrium dendrobatidis from the basal fungal phylum, the Chytridiomycota, which has emerged to cause global amphibian declines and extinctions. Genomics is revolutionizing our view of the fungal kingdom, and genome sequences for zygomycete pathogens (Rhizopus, Mucor), skin-associated fungi (dermatophytes, Malassezia), and the Candida pathogenic species clade promise to provide insights into the origins of virulence. Here we survey the diversity of fungal pathogens and illustrate key principles revealed by genomics involving sexual reproduction and sex determination, loss of conserved pathways in derived fungal lineages that are retained in basal fungi, and shared and divergent virulence strategies of successful human pathogens, including dimorphic and trimorphic transitions in form. The overarching conclusion is that fungal pathogens of animals have arisen repeatedly and independently throughout the fungal tree of life, and while they share general properties, there are also unique features to the virulence strategies of each successful microbial pathogen.     

退化高寒草地土壤真菌群落与土壤环境因子间相互关系

【目的】为探究祁连山高寒草地退化过程中土壤真菌群落分布特征与土壤环境因子间的相互关系。【方法】利用Illumina Miseq PE250高通量测序技术对轻度、中度和重度退化草地土壤真菌群落结构变化及其多样性进行分析,并对土壤真菌群落与土壤环境因子的相互关系进行冗余分析(RDA)。【结果】随着退化程度加剧,土壤pH呈现出升高趋势,电导率呈现出先升高后降低趋势,土壤含水量、有机碳、全氮、全磷和全钾含量均逐渐降低。高通量测序共得到750575条有效序列和5788个OTUs;各试验点样地中真菌群落Chao1指数和Shannon-Wiener指数变化各异。在门分类水平上,子囊菌门(Ascomycota)、担子菌门(Basidiomycota)、接合菌门(Zygomycota)、球囊菌门(Glomeromycota)和壶菌门(Chytridiomycota)是各草地土壤的优势类群。RDA分析表明,土壤速效钾、全氮、速效氮和有机碳是祁连山不同退化高寒草地土壤真菌群落分布的主要驱动因子。【结论】祁连山不同退化高寒草地土壤真菌群落间差异明显,土壤环境因子是影响土壤真菌群落分布的重要因素。     

Potential bias of fungal 18S rDNA and internal transcribed spacer polymerase chain reaction primers for estimating fungal biodiversity in soil

Four fungal 18S rDNA and internal transcribed spacer (ITS) polymerase chain reaction (PCR) primer pairs were tested for their specificity towards target fungal DNA in soil DNA extracts, and their ability to assess the diversity of fungal communities in a natural grassland soil was compared. Amplified PCR products were cloned, and approximately 50 clones from each library were sequenced. Phylogenetic analysis and database searches indicated that each of the sequenced cloned DNA fragments was of fungal origin for each primer pair, with the exception of the sequences generated using the 18S rDNA primers nu-SSU-0817 and nu-SSU-1196, where 35 of the 50 sequenced clones represented soil invertebrates. Although some of the primers have previously been suggested to be biased towards certain fungal taxonomic groups, the ratio of sequences representing each of the four main fungal phyla, Ascomycota, Basidiomycota, Chytridiomycota and Zygomycota, was similar for each of the primer pairs, suggesting that primer bias may be less significant than previously thought. Collector's curves were plotted to estimate the coverage obtained for each of the clone libraries after clustering the sequences into operational taxonomic units at a level of 99% sequence similarity. The curves indicated that good coverage of diversity was achieved, with the exception of the clone library constructed using primers nu-SSU-0817 and nu-SSU-1196, on account of the high number of non-fungal sequences obtained. The work demonstrates the usefulness of 18S rDNA and ITS PCR primers for assessing fungal diversity in environmental samples, and it also highlights some potential limitations of the approach with respect to PCR primer specificity and bias.