Deforestation impacts soil biodiversity and ecosystem services worldwide.

Deforestation poses a global threat to biodiversity and its capacity to deliver ecosystem services. Yet, the impacts of deforestation on soil biodiversity and its associated ecosystem services remain virtually unknown. We generated a global dataset including 696 paired-site observations to investigate how native forest conversion to other land uses affects soil properties, biodiversity, and functions associated with the delivery of multiple ecosystem services. The conversion of native forests to plantations, grasslands, and croplands resulted in higher bacterial diversity and more homogeneous fungal communities dominated by pathogens and with a lower abundance of symbionts. Such conversions also resulted in significant reductions in carbon storage, nutrient cycling, and soil functional rates related to organic matter decomposition. Responses of the microbial community to deforestation, including bacterial and fungal diversity and fungal guilds, were predominantly regulated by changes in soil pH and total phosphorus. Moreover, we found that soil fungal diversity and functioning in warmer and wetter native forests is especially vulnerable to deforestation. Our work highlights that the loss of native forests to managed ecosystems poses a major global threat to the biodiversity and functioning of soils and their capacity to deliver ecosystem services.

Complete Mitochondrial Genomes of Nedyopus patrioticus: New Insights into the Color Polymorphism of Millipedes.

There has been debate about whether individuals with different color phenotypes should have different taxonomic status. In order to determine whether the different color phenotypes of Nedyopus patrioticus require separate taxonomic status or are simply synonyms, here, the complete mitochondrial genomes (mitogenomes) of two different colored N. patrioticus, i.e., red N. patrioticus and white N. patrioticus, are presented. The two mitogenomes were 15,781 bp and 15,798 bp in length, respectively. Each mitogenome contained 13 PCGs, 19 tRNAs, 2 rRNAs, and 1 CR, with a lack of trnI, trnL2, and trnV compared to other Polydesmida species. All genes were located on a single strand in two mitogenomes. Mitochondrial DNA analyses revealed that red N. patrioticus and white N. patrioticus did not show clear evolutionary differences. Furthermore, no significant divergence was discovered by means of base composition analysis. As a result, we suggest that white N. patrioticus might be regarded as a synonym for red N. patrioticus. The current findings confirmed the existence of color polymorphism in N. patrioticus, which provides exciting possibilities for future research. It is necessary to apply a combination of molecular and morphological methods in the taxonomy of millipedes.

Transcriptome-wide analysis of a superior xylan degrading isolate Penicillium oxalicum 5-18 revealed active lignocellulosic degrading genes.

Lignocellulose biomass raw materials have a high value in energy conversion. Recently, there has been growing interest in using microorganisms to secret a series of enzymes for converting low-cost biomass into high-value products such as biofuels. We previously isolated a strain of Penicillium oxalicun 5-18 with promising lignocellulose-degrading capability. However, the mechanisms of lignocellulosic degradation of this fungus on various substrates are still unclear. In this study, we performed transcriptome-wide profiling and comparative analysis of strain 5-18 cultivated in liquid media with glucose (Glu), xylan (Xyl) or wheat bran (WB) as sole carbon source. In comparison to Glu culture, the number of differentially expressed genes (DEGs) induced by WB and Xyl was 4134 and 1484, respectively, with 1176 and 868 genes upregulated. Identified DEGs were enriched in many of the same pathways in both comparison groups (WB vs. Glu and Xly vs. Glu). Specially, 118 and 82 CAZyme coding genes were highly upregulated in WB and Xyl cultures, respectively. Some specific pathways including (Hemi)cellulose metabolic processes were enriched in both comparison groups. The high upregulation of these genes also confirmed the ability of strain 5-18 to degrade lignocellulose. Co-expression and co-upregulated of genes encoding CE and AA CAZy families, as well as other (hemi)cellulase revealed a complex degradation strategy in this strain. Our findings provide new insights into critical genes, key pathways and enzyme arsenal involved in the biomass degradation of P. oxalicum 5-18.

Fungal necromass is reduced by intensive drought in subsoil but not in topsoil.

The frequency and intensity of droughts worldwide are challenging the conservation of soil organic carbon (SOC) pool. Microbial necromass is a key component of SOC, but how it responds to drought at specific soil depths remains largely unknown. Here, we conducted a 3-year field experiment in a forest plantation to investigate the impacts of drought intensities under three treatments (ambient control [CK], moderate drought [30% throughfall removal], and intensive drought [50% throughfall removal]) on soil microbial necromass pools (i.e., bacterial necromass carbon, fungal necromass carbon, and total microbial necromass carbon). We showed that the effects of drought on microbial necromass depended on microbial groups, soil depth, and drought intensity. While moderate drought increased total (+9.1% ± 3.3%) and fungal (+13.5% ± 4.9%) necromass carbon in the topsoil layer (0-15 cm), intensive drought reduced total (-31.6% ± 3.7%) and fungal (-43.6% ± 4.0%) necromass in the subsoil layer (15-30 cm). In contrast, both drought treatments significantly increased the BNC in the topsoil and subsoil. Our results suggested that the effects of drought on the microbial necromass of the subsoil were more pronounced than those of the topsoil. This study highlights the complex responses of microbial necromass to drought events depending on microbial community structure, drought intensity and soil depth with global implications when forecasting carbon cycling under climate change.

A High-Quality Genome Sequence of the Penicillium oxalicum 5-18 Strain Isolated from a Poplar Plantation Provides Insights into Its Lignocellulose Degradation.

Studies on the degradation of plant cell wall polysaccharides by fungal extracellular enzymes have attracted recent attention from researchers. Xylan, abundant in hemicellulose, that play great role in connection between cellulose and lignin, has seen interest in its hydrolytic enzymatic complex. In this study, dozens of fungus species spanning genera were isolated from rotting leaves based on their ability to decompose xylan. Among these isolates, a strain with strong xylanase-producing ability was selected for further investigation by genome sequencing. Based on phylogenetic analysis of ITS (rDNA internal transcribed spacer) and LSU (Large subunit 28S rDNA) regions, the isolate was identified as Penicillium oxalicum. Morphological analysis also supported this finding. Xylanase activity of this isolated P. oxalicum 5-18 strain was recorded to be 30.83 U/mL using the 3,5-dinitro-salicylic acid (DNS) method. Further genome sequencing reveals that sequenced reads were assembled into a 30.78 Mb genome containing 10,074 predicted protein-encoding genes. In total, 439 carbohydrate-active enzymes (CAZymes) encoding genes were predicted, many of which were associated with cellulose, hemicellulose, pectin, chitin and starch degradation. Further analysis and comparison showed that the isolate P. oxalicum 5-18 contains a diverse set of CAZyme genes involved in degradation of plant cell wall components, particularly cellulose and hemicellulose. These findings provide us with valuable genetic information about the plant biomass-degrading enzyme system of P. oxalicum, facilitating a further exploration of the repertoire of industrially relevant lignocellulolytic enzymes of P. oxalicum 5-18.

Phosphorus additions imbalance terrestrial ecosystem C:N:P stoichiometry.

Carbon (C):nitrogen (N):phosphorus (P) stoichiometry in plants, soils, and microbial biomass influences productivity and nutrient cycling in terrestrial ecosystems. Anthropogenic inputs of P to ecosystems are increasing; however, our understanding of the impacts of P addition on terrestrial ecosystem C:N:P ratios remains elusive. By conducting a meta-analysis with 1413 paired observations from 121 publications, we showed that P addition significantly decreased plant, soil, and microbial biomass N:P and C:P ratios, but had negligible effects on C:N ratios. The reductions in N:P and C:P ratios became more evident as the P application rates and experimental duration increased. The P addition effects on terrestrial ecosystem C:N:P stoichiometry did not vary with ecosystem types or climates. Moreover, the responses of N:P and C:P ratios in soil and microbial biomass were associated with the responses of soil pH and fungi:bacteria ratios. Additionally, P additions increased net primary productivity, microbial biomass, soil respiration, N mineralization, and N nitrification, but decreased ammonium and nitrate contents. Decreases in plant N:P and C:P ratios were both negatively correlated to net primary productivity and soil respiration, but positively correlated to ammonium and nitrate contents; microbial biomass, soil respiration, ammonium contents, and nitrate contents all increased with declining soil N:P and C:P ratios. Our findings highlight that P additions could imbalance C:N:P stoichiometry and potentially impact the terrestrial ecosystem functions.

Characterization of the complete mitochondrial genome of Drawida gisti (Metagynophora, Moniligastridae) and comparison with other Metagynophora species.

Here, the complete mitochondrial genome (mitogenome) of Drawida gisti was sequenced and compared with the mitogenomes of other Metagynophora species. The circular mitogenome was 14,648 bp in length and contained two ribosomal RNA genes (rRNAs), 13 protein-coding genes (PCGs), and 22 transfer RNA genes (tRNAs). The types of constitutive genes and the direction of the coding strand that appeared in Drawida mitogenome were identical to those observed in other Metagynophora species, except for a missing lengthy non-coding region. The conservative relationships between Drawida species were supported by the overall analyses of 13 PCGs, two rRNAs, and 22 tRNAs. A comparison of the Metagynophora mitogenomes revealed that the ATP8 gene possessed the highest polymorphism among the 13 PCGs and two rRNAs. Phylogenetic analysis suggested that the Moniligastridae contained Drawida, which is a primitive Metagynophora group. Our study provides a step forward toward elucidating the evolutionary linkages within Drawida and even Metagynophora.

Identification, evolution, expression, and docking studies of fatty acid desaturase genes in wheat (Triticum aestivum L.).

BACKGROUNDS: Fatty acid desaturases (FADs) introduce a double bond into the fatty acids acyl chain resulting in unsaturated fatty acids that have essential roles in plant development and response to biotic and abiotic stresses. Wheat germ oil, one of the important by-products of wheat, can be a good alternative for edible oils with clinical advantages due to the high amount of unsaturated fatty acids. Therefore, we performed a genome-wide analysis of the wheat FAD gene family (TaFADs). RESULTS: 68 FAD genes were identified from the wheat genome. Based on the phylogenetic analysis, wheat FADs clustered into five subfamilies, including FAB2, FAD2/FAD6, FAD4, DES/SLD, and FAD3/FAD7/FAD8. The TaFADs were distributed on chromosomes 2A-7B with 0 to 10 introns. The Ka/Ks ratio was less than one for most of the duplicated pair genes revealed that the function of the genes had been maintained during the evolution. Several cis-acting elements related to hormones and stresses in the TaFADs promoters indicated the role of these genes in plant development and responses to environmental stresses. Likewise, 72 SSRs and 91 miRNAs in 36 and 47 TaFADs have been identified. According to RNA-seq data analysis, the highest expression in all developmental stages and tissues was related to TaFAB2.5, TaFAB2.12, TaFAB2.15, TaFAB2.17, TaFAB2.20, TaFAD2.1, TaFAD2.6, and TaFAD2.8 genes while the highest expression in response to temperature stress was related to TaFAD2.6, TaFAD2.8, TaFAB2.15, TaFAB2.17, and TaFAB2.20. Furthermore, docking simulations revealed several residues in the active site of TaFAD2.6 and TaFAD2.8 in close contact with the docked oleic acid that could be useful in future site-directed mutagenesis studies to increase the catalytic efficiency of them and subsequently improve agronomic quality and tolerance of wheat against environmental stresses. CONCLUSIONS: This study provides comprehensive information that can lead to the detection of candidate genes for wheat genetic modification.

Development and characterization of microsatellite markers in the earthworm Drawida gisti Michaelsen, 1931 and cross-amplification in two other congeners.

The earthworm (Drawida gisti) is an ecologically important sentinel species for soils that is widely distributed throughout Eastern Asia; however, the molecular tools required for genetic diversity studies of this earthworm are still rare. The aim of the study was to develop and characterize microsatellite markers in D. gisti and to evaluate their transferability to other Drawida species. We employed a RAD-seq approach to develop 12 microsatellite markers for D. gisti. The characterization and analysis of loci was achieved using 24 individuals of D. gisti from a natural population. The number of alleles per locus ranged from four to eleven, with an average of 6.5. Observed and expected heterozygosities varied from 0.708 to 0.958, and from 0.568 to 0.883, respectively. No loci presented significant deviations from the Hardy-Weinberg equilibrium, while linkage disequilibrium was detected between three loci. Cross-species amplification tests suggested that the transferability of ten loci was positive for the two congeners D. japonica and D. ghilarovi. This set of microsatellite markers may be used to evaluate the genetic diversity and population structures of D. gisti and related species in the future.

Feasible Green Strategy for the Quantitative Bioaccumulation of Heavy Metals by Lemna minor: Application of the Self-Thinning Law.

This study involved the development of mathematical linear regression models to describe the relationships between mean plant biomass (M) and population density (D), M and frond diameter (L), frond numbers (N) and L of Lemna minor under different initial population densities (3200, 4450, and 6400 plants/m2), respectively, from the perspective of the self-thinning law. Our results revealed that the value of the allometric exponents for M and D were - 3/2. Further, the concentrations of Zn, Pb, Cu, Fe, and Ni accumulated in L. minor plants were 0.86, 0.32, 0.36, 0.62, and 0.39 mg/kg, respectively. Based on these developed equations and the heavy metal accumulations by L. minor, the phytoremediation capacity of L. minor was quantified via its frond diameters. Overall, the present study provides a cost-effective green method for managing the phytoremediation of heavy metal-contaminated aquatic environments.

Comparative study on the biodegradation of chlorpyrifos-methyl by Bacillus megaterium CM-Z19 and Pseudomonas syringae CM-Z6.

The strains CM-Z19 and CM-Z6, which are capable of highly degrading chlorpyrifos-methyl, were isolated from soil. They were identified as Bacillus megaterium CM-Z19 and Pseudomonas syringae CM-Z6, respectively, based on the 16S rRNA and an analysis of their morphological, physiological and biochemical characteristics. The strain CM-Z19 showed 92.6% degradation of chlorpyrifos-methyl (100 mg/L) within 5 days of incubation, and the strain CM-Z6 was 99.1% under the same conditions. In addition, the degradation characteristics of the two strains were compared and studied, and the results showed that the strain CM-Z19 had higher phosphoesterase activity and ability to degrade the organophosphorus pesticide than did the strain CM-Z6. However, the strain CM-Z19 could not degrade its first hydrolysis metabolite 3,5,6-trichloro-2-pyridinol (TCP) and could not completely degrade chlorpyrifos-methyl. The strain CM-Z6 could effectively degrade TCP and could degrade chlorpyrifos-methyl more quickly than strain CM-Z19.

Strategies to Increase On-Target and Reduce Off-Target Effects of the CRISPR/Cas9 System in Plants.

The CRISPR/Cas9 system (clustered regularly interspaced short palindromic repeat-associated protein 9) is a powerful genome-editing tool in animals, plants, and humans. This system has some advantages, such as a high on-target mutation rate (targeting efficiency), less cost, simplicity, and high-efficiency multiplex loci editing, over conventional genome editing tools, including meganucleases, transcription activator-like effector nucleases (TALENs), and zinc finger nucleases (ZFNs). One of the crucial shortcomings of this system is unwanted mutations at off-target sites. We summarize and discuss different approaches, such as dCas9 and Cas9 paired nickase, to decrease the off-target effects in plants. According to studies, the most effective method to reduce unintended mutations is the use of ligand-dependent ribozymes called aptazymes. The single guide RNA (sgRNA)/ligand-dependent aptazyme strategy has helped researchers avoid unwanted mutations in human cells and can be used in plants as an alternative method to dramatically decrease the frequency of off-target mutations. We hope our concept provides a new, simple, and fast gene transformation and genome-editing approach, with advantages including reduced time and energy consumption, the avoidance of unwanted mutations, increased frequency of on-target changes, and no need for external forces or expensive equipment.

Analyzing the impact of climate and management factors on the productivity and soil carbon sequestration of poplar plantations.

It is crucial to investigate how climate and management factors impact poplar plantation production and soil carbon sequestration interactively. We extracted above-ground net primary production (ANPP), climate and management factors from peer-reviewed journal articles and analyzed impact of management factor and climate on the mean annual increment (MAI) of poplar ANPP statistically. Previously validated mechanistic model (ED) is used to perform case simulations for managed poplar plantations under different harvesting rotations. The meta-analysis indicate that the dry matter MAI was 6.3 Mg ha(-1) yr(-1) (n=641, sd=4.9) globally, and 5.1 (n=292, sd=4.0), 8.1 (n=224, sd=4.7) and 4.4 Mg ha(-1) yr(-1) (n=125, sd=3.2) in Europe, the US and China, respectively. Poplar MAI showed a significant response to GDD, precipitation and planting density and formed a quadratic relationship with stand age. The low annual production for poplar globally was probably caused by suboptimal water availability, rotation length and planting density. SEM attributes the variance of poplar growth rate more to climate than to management effects. Case simulations indicated that longer rotation cycle significantly increased soil carbon storage. Findings of this work suggests that management factor of rotation cycle alone could have dramatic impact on the above ground growth, as well as on the soil carbon sequestration of poplar plantations and will be helpful to quantify the long-term carbon sequestration through short rotation plantation. The findings of this study are useful in guiding further research, policy and management decisions towards sustainable poplar plantations.

Decolorization characteristics of a newly isolated salt-tolerant Bacillus sp. strain and its application for azo dye-containing wastewater in immobilized form.

Strain CICC 23870 capable of decolorization of various azo dyes under high saline conditions was isolated from saline-alkali soil. The oxygen-insensitive azoreductase in crude extracts exhibited a wide substrate adaptively in the presence of NADH as a cofactor. The decolorization process by free cells followed first-order kinetics, with a high Methyl Orange (MO) tolerance concentration up to 100 mg l(-1) estimated by Haldane model. The average decolorization rate of free cell system was 26.30 mg g(-1) h(-1) at initial MO concentration of 32.7 mg l(-1). However, the values for the systems of immobilized cells (4 mm) in alginate, alginate and nano-TiO2, and alginate and powered activated carbon (PAC) were 6.83, 4.64, and 11.34 mg g(-1) h(-1), respectively. The effective diffusion factors in the tree different matrices were calculated by diffusion-based mathematic model. The diffusion step controls the overall decolorization rate, and the effective diffusion coefficients varied with internal structure of the bead matrices. The diffusion coefficients were increased from 4.98 × 10(-9) to 2.25 × 10(-8) cm(2) s(-1) when PAC was added, but decreased to 6.62 × 10(-10) cm(2) s(-1) when nano-TiO2 was added. The immobilized matrices could be reused for at least three cycles but with a decreased decolorization rate, possibly due to the breakage of beads at the end of each cycle, which led to the loss of immobilized bacteria.

Implications of carbon catabolite repression for Aspergillus-based cell factories: A review.

Carbon catabolite repression (CCR) is a global regulatory mechanism that allows organisms to preferentially utilize a preferred carbon source (usually glucose) by suppressing the expression of genes associated with the utilization of nonpreferred carbon sources. Aspergillus is a large genus of filamentous fungi, some species of which have been used as microbial cell factories for the production of organic acids, industrial enzymes, pharmaceuticals, and other fermented products due to their safety, substrate convenience, and well-established post-translational modifications. Many recent studies have verified that CCR-related genetic alterations can boost the yield of various carbohydrate-active enzymes (CAZymes), even under CCR conditions. Based on these findings, we emphasize that appropriate regulation of the CCR pathway, especially the expression of the key transcription factor CreA gene, has great potential for further expanding the application of Aspergillus cell factories to develop strains for industrial CAZymes production. Further, the genetically modified CCR strains (chassis hosts) can also be used for the production of other useful natural products and recombinant proteins, among others. We here review the regulatory mechanisms of CCR in Aspergillus and its direct application in enzyme production, as well as its potential application in organic acid and pharmaceutical production to illustrate the effects of CCR on Aspergillus cell factories.

Biodegradation of Low Density Polyethylene by the Fungus Cladosporium sp. Recovered from a Landfill Site.

Low density polyethylene (LDPE) has been widely used commercially for decades; however, as a non-degradable material, its continuous accumulation has contributed to serious environmental issues. A fungal strain, Cladosporium sp. CPEF-6 exhibiting a significant growth advantage on MSM-LDPE (minimal salt medium), was isolated and selected for biodegradation analysis. LDPE biodegradation was analyzed by weight loss percent, change in pH during fungal growth, environmental scanning electron microscopy (ESEM), and Fourier transformed infrared spectroscopy (FTIR). Inoculation with the strain Cladosporium sp. CPEF-6 resulted in a 0.30 ± 0.06% decrease in the weight of untreated LDPE (U-LDPE). After heat treatment (T-LDPE), the weight loss of LDPE increased significantly and reached 0.43 ± 0.01% after 30 days of culture. The pH of the medium was measured during LDPE degradation to assess the environmental changes caused by enzymes and organic acids secreted by the fungus. The fungal degradation of LDPE sheets was characterized by ESEM analysis of topographical alterations, such as cracks, pits, voids, and roughness. FTIR analysis of U-LDPE and T-LDPE revealed the appearance of novel functional groups associated with hydrocarbon biodegradation as well as changes in the polymer carbon chain, confirming the depolymerization of LDPE. This is the first report demonstrating the capacity of Cladosporium sp. to degrade LDPE, with the expectation that this finding can be used to ameliorate the negative impact of plastics on the environment.

A Method to Reduce off-Targets in CRISPR/Cas9 System in Plants.

One of the strategies to reduce the off-target mutations in CRISPR/Cas9 system is to use the temperature-independent gene transformation method. Mesoporous silica nanoparticles (MSNs)-gene delivery system is temperature-independent; thus, it can transfer the interesting plasmid (pDNA) to the target plant at different temperatures, including 37 °C. Due to the high activity of SpCas9 at 37 °C compared to lower temperatures, on-target mutagenesis increases at 37 °C. Therefore, we describe the synthesis of the functionalized MSNs with the particle size of less than 40 nm, binding pDNA to the MSNs, and transferring of the pDNA-MSNs into the target plants.

The Complete Mitochondrial Genome of Spirobolus bungii (Diplopoda, Spirobolidae): The First Sequence for the Genus Spirobolus.

Millipedes (Diplopoda) comprise one of the most important groups of large soil arthropods in terrestrial ecosystems; however, their phylogenetic relationships are poorly understood. Herein, the mitochondrial genome (mitogenome) of Spirobolus bungii was sequenced and annotated, which was 14,879 bp in size and included 37 typical mitochondrial genes (13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs)). Most of the 13 PCGs had ATN (AT/A/T/G) as the start codon except for COX1, which used CGA, and most PCGs ended with the T end codon. By comparing the gene arrangements of the mitogenomes among Diplopoda species, rearrangement occurred between and within orders. In contrast to Narceus annularus, the mitogenome genes of S. bungii had consistent orders but were transcribed in completely opposite directions, which was a novel finding in Spirobolidae. Moreover, the phylogenetic relationships within Diplopoda, which were based on the sequences of 13 PCGs, showed that S. bungii was clustered with N. annularus, followed by Abacion magmun. This indicated that there might be a close relationship between Callipodida and Spirobolida. These results could contribute to further studies on the genetics and evolutionary processes of S. bungii and other Diplopoda species.

'biogeom': An R package for simulating and fitting natural shapes.

Many natural objects exhibit radial or axial symmetry in a single plane. However, a universal tool for simulating and fitting the shapes of such objects is lacking. Herein, we present an R package called 'biogeom' that simulates and fits many shapes found in nature. The package incorporates novel universal parametric equations that generate the profiles of bird eggs, flowers, linear and lanceolate leaves, seeds, starfish, and tree-rings, and three growth-rate equations that generate the profiles of ovate leaves and the ontogenetic growth curves of animals and plants. 'biogeom' includes several empirical datasets comprising the boundary coordinates of bird eggs, fruits, lanceolate and ovate leaves, tree rings, seeds, and sea stars. The package can also be applied to other kinds of natural shapes similar to those in the datasets. In addition, the package includes sigmoid curves derived from the three growth-rate equations, which can be used to model animal and plant growth trajectories and predict the times associated with maximum growth rate. 'biogeom' can quantify the intra- or interspecific similarity of natural outlines, and it provides quantitative information of shape and ontogenetic modification of shape with important ecological and evolutionary implications for the growth and form of the living world.

Soil nutrients and plant diversity affect ectomycorrhizal fungal community structure and functional traits across three subalpine coniferous forests.

The symbiotic relationship between ectomycorrhizal fungi (EMF) and the roots of host plants is significantly important in regulating the health and stability of ecosystems, especially of those such as the climate warming affected subalpine forest ecosystems. Therefore, from the coniferous forest systems located in the Southern Qinghai-Tibetan Plateau, root tips from three forest tree species: Pinus wallichiana, Abies spectabilis and Picea spinulosa, were collected to look for the local causes of EMF community composition and diversity patterns. The EMF colonization rate, diversity and taxonomic community structure were determined by morphotyping and sanger sequencing of the fungal ITS gene from the root tip samples. Soil exploration types were identified based on the morphologies of the ectomycorrhizas, coupled with soil properties analysis and plant diversity survey. Contrasting patterns of EMF community and functional diversity were found across the studied three forests types dominated by different coniferous tree species. In terms of associations between soil and EMF properties, the total phosphorus (TP) and nitrate (NO3 -) contents in soil negatively correlated with the colonization rate and the Shannon diversity index of EMF in contrast to the positive relationship between TP and EMF richness. The soil total nitrogen (TN), ammonium (NH4 +) and plant diversity together caused 57.6% of the total variations in the EMF taxonomic community structure at the three investigated forest systems. Whereas based on the soil exploration types alone, NH4 + and TN explained 74.2% of variance in the EMF community structures. Overall, the findings of this study leverage our understanding of EMF dynamics and local influencing factors in coniferous forests dominated by different tree species within the subalpine climatic zone.