Climate acts as an environmental filter to plant pathogens.

Climate shapes the distribution of plant-associated microbes such as mycorrhizal and endophytic fungi. However, the role of climate in plant pathogen community assembly is less understood. Here, we explored the role of climate in the assembly of Phytophthora communities at >250 sites along a latitudinal gradient from Spain to northern Sweden and an altitudinal gradient from the Spanish Pyrenees to lowland areas. Communities were detected by ITS sequencing of river filtrates. Mediation analysis supported the role of climate in the biogeography of Phytophthora and ruled out other environmental factors such as geography or tree diversity. Comparisons of functional and species diversity showed that environmental filtering dominated over competitive exclusion in Europe. Temperature and precipitation acted as environmental filters at different extremes of the gradients. In northern regions, winter temperatures acted as an environmental filter on Phytophthora community assembly, selecting species adapted to survive low minimum temperatures. In southern latitudes, a hot dry climate was the main environmental filter, resulting in communities dominated by drought-tolerant Phytophthora species with thick oospore walls, a high optimum temperature for growth, and a high maximum temperature limit for growth. By taking a community ecology approach, we show that the establishment of Phytophthora plant pathogens in Europe is mainly restricted by cold temperatures.

Multiple rearrangements and low inter- and intra-species mitogenome sequence variation in the Heterobasidion annosum s.l. species complex.

Introduction: Mitochondria are essential organelles in the eukaryotic cells and responsible for the energy production but are also involved in many other functions including virulence of some fungal species. Although the evolution of fungal mitogenomes have been studied at some taxonomic levels there are still many things to be learned from studies of closely related species. Methods: In this study, we have analyzed 60 mitogenomes in the five species of the Heterobasidion annosum sensu lato complex that all are necrotrophic pathogens on conifers. Results and Discussion: Compared to other fungal genera the genomic and genetic variation between and within species in the complex was low except for multiple rearrangements. Several translocations of large blocks with core genes have occurred between the five species and rearrangements were frequent in intergenic areas. Mitogenome lengths ranged between 108 878 to 116 176 bp, mostly as a result of intron variation. There was a high degree of homology of introns, homing endonuclease genes, and intergenic ORFs among the five Heterobasidion species. Three intergenic ORFs with unknown function (uORF6, uORF8 and uORF9) were found in all five species and was located in conserved synteny blocks. A 13 bp long GC-containing self-complementary palindrome was discovered in many places in the five species that were optional in presence/absence. The within species variation is very low, among 48 H. parviporum mitogenomes, there was only one single intron exchange, and SNP frequency was 0.28% and indel frequency 0.043%. The overall low variation in the Heterobasidion annosum sensu lato complex suggests a slow evolution of the mitogenome.

A Reliable and Simple Method for the Production of Viable Pycnidiospores of the Pine Pathogen Diplodia sapinea and a Spore-Based Infection Assay on Scots Pine.

Diplodia sapinea is a globally distributed opportunistic fungal pathogen of conifers that causes severe production losses in forestry. The fungus frequently colonizes pine trees as an endophyte without causing visible symptoms but can become pathogenic when the host plant is weakened by stress, such as drought or heat. Forest damage might therefore further increase due to the effects of climate change. The future development of control strategies depends on a better understanding of the fungus' biology, which requires experimental methods for its investigation in the laboratory. An efficient, standardized protocol for the production and storage of highly viable pycnidiospores was developed, and a spore-based infection method was devised. We compared infection rates of dormant and actively growing, wounded, or nonwounded Scots pine seedlings inoculated with in vitro-produced spores and mycelium from agar-plugs. Spores were a much more efficient inoculum for causing disease symptoms on wounded plants than the conventional agar plug. The application of spores on nonwounded plants lead to high rates of asymptomatic infection, suggesting endophytic fungal development. These methods enable standardized spore infection and virulence assays and promote D. sapinea as a model organism for studying the switch from endophytic to pathogenic life styles of forest pathogens.

Interaction of drought- and pathogen-induced mortality in Norway spruce and Scots pine.

Pathogenic diseases frequently occur in drought-stressed trees. However, their contribution to the process of drought-induced mortality is poorly understood. We combined drought and stem inoculation treatments to study the physiological processes leading to drought-induced mortality in Norway spruce (Picea abies) and Scots pine (Pinus sylvestris) saplings infected with Heterobasidion annosum s.s. We analysed the saplings' water status, gas exchange, nonstructural carbohydrates (NSCs) and defence responses, and how they related to mortality. Saplings were followed for two growing seasons, including an artificially induced 3-month dormancy period. The combined drought and pathogen treatment significantly increased spruce mortality; however, no interaction between these stressors was observed in pine, although individually each stressor caused mortality. Our results suggest that pathogen infection decreased carbon reserves in spruce, reducing the capacity of saplings to cope with drought, resulting in increased mortality rates. Defoliation, relative water content and the starch concentration of needles were predictors of mortality in both species under drought and pathogen infection. Infection and drought stress create conflicting needs for carbon to compartmentalize the pathogen and to avoid turgor loss, respectively. Heterobasidion annosum reduces the functional sapwood area and shifts NSC allocation patterns, reducing the capacity of trees to cope with drought.

Community composition of aquatic fungi across the thawing Arctic.

Thermokarst activity at permafrost sites releases considerable amounts of ancient carbon to the atmosphere. A large part of this carbon is released via thermokarst ponds, and fungi could be an important organismal group enabling its recycling. However, our knowledge about aquatic fungi in thermokarstic systems is extremely limited. In this study, we collected samples from five permafrost sites distributed across circumpolar Arctic and representing different stages of permafrost integrity. Surface water samples were taken from the ponds and, additionally, for most of the ponds also the detritus and sediment samples were taken. All the samples were extracted for total DNA, which was then amplified for the fungal ITS2 region of the ribosomal genes. These amplicons were sequenced using PacBio technology. Water samples were also collected to analyze the chemical conditions in the ponds, including nutrient status and the quality and quantity of dissolved organic carbon. This dataset gives a unique overview of the impact of the thawing permafrost on fungal communities and their potential role on carbon recycling.

Declining fungal diversity in Arctic freshwaters along a permafrost thaw gradient.

Climate change-driven permafrost thaw has a strong influence on pan-Arctic regions, via, for example, the formation of thermokarst ponds. These ponds are hotspots of microbial carbon cycling and greenhouse gas production, and efforts have been put on disentangling the role of bacteria and archaea in recycling the increasing amounts of carbon arriving to the ponds from degrading watersheds. However, despite the well-established role of fungi in carbon cycling in the terrestrial environments, the interactions between permafrost thaw and fungal communities in Arctic freshwaters have remained unknown. We integrated data from 60 ponds in Arctic hydro-ecosystems, representing a gradient of permafrost integrity and spanning over five regions, namely Alaska, Greenland, Canada, Sweden, and Western Siberia. The results revealed that differences in pH and organic matter quality and availability were linked to distinct fungal community compositions and that a large fraction of the community represented unknown fungal phyla. Results display a 16%-19% decrease in fungal diversity, assessed by beta diversity, across ponds in landscapes with more degraded permafrost. At the same time, sites with similar carbon quality shared more species, aligning a shift in species composition with the quality and availability of terrestrial dissolved organic matter. We demonstrate that the degradation of permafrost has a strong negative impact on aquatic fungal diversity, likely via interactions with the carbon pool released from ancient deposits. This is expected to have implications for carbon cycling and climate feedback loops in the rapidly warming Arctic.

Killing two enemies with one stone? Genomics of resistance to two sympatric pathogens in Norway spruce.

Trees must cope with the attack of multiple pathogens, often simultaneously during their long lifespan. Ironically, the genetic and molecular mechanisms controlling this process are poorly understood. The objective of this study was to compare the genetic component of resistance in Norway spruce to Heterobasidion annosum s.s. and its sympatric congener Heterobasidion parviporum. Heterobasidion root- and stem-rot is a major disease of Norway spruce caused by members of the Heterobasidion annosum species complex. Resistance to both pathogens was measured using artificial inoculations in half-sib families of Norway spruce trees originating from central to northern Europe. The genetic component of resistance was analysed using 63,760 genome-wide exome-capture sequenced SNPs and multitrait genome-wide associations. No correlation was found for resistance to the two pathogens; however, associations were found between genomic variants and resistance traits with synergic or antagonist pleiotropic effects to both pathogens. Additionally, a latitudinal cline in resistance in the bark to H. annosum s.s. was found; trees from southern latitudes, with a later bud-set and thicker stem diameter, allowed longer lesions, but this was not the case for H. parviporum. In summary, this study detects genomic variants with pleiotropic effects which explain multiple disease resistance from a genic level and could be useful for selection of resistant trees to both pathogens. Furthermore, it highlights the need for additional research to understand the evolution of resistance traits to multiple pathogens in trees.

Comparative analyses of the Hymenoscyphus fraxineus and Hymenoscyphus albidus genomes reveals potentially adaptive differences in secondary metabolite and transposable element repertoires.

BACKGROUND: The dieback epidemic decimating common ash (Fraxinus excelsior) in Europe is caused by the invasive fungus Hymenoscyphus fraxineus. In this study we analyzed the genomes of H. fraxineus and H. albidus, its native but, now essentially displaced, non-pathogenic sister species, and compared them with several other members of Helotiales. The focus of the analyses was to identify signals in the genome that may explain the rapid establishment of H. fraxineus and displacement of H. albidus. RESULTS: The genomes of H. fraxineus and H. albidus showed a high level of synteny and identity. The assembly of H. fraxineus is 13 Mb longer than that of H. albidus', most of this difference can be attributed to higher dispersed repeat content (i.e. transposable elements [TEs]) in H. fraxineus. In general, TE families in H. fraxineus showed more signals of repeat-induced point mutations (RIP) than in H. albidus, especially in Long-terminal repeat (LTR)/Copia and LTR/Gypsy elements. Comparing gene family expansions and 1:1 orthologs, relatively few genes show signs of positive selection between species. However, several of those did appeared to be associated with secondary metabolite genes families, including gene families containing two of the genes in the H. fraxineus-specific, hymenosetin biosynthetic gene cluster (BGC). CONCLUSION: The genomes of H. fraxineus and H. albidus show a high degree of synteny, and are rich in both TEs and BGCs, but the genomic signatures also indicated that H. albidus may be less well equipped to adapt and maintain its ecological niche in a rapidly changing environment.

Unravelling hybridization in Phytophthora using phylogenomics and genome size estimation.

The genus Phytophthora comprises many economically and ecologically important plant pathogens. Hybrid species have previously been identified in at least six of the 12 phylogenetic clades. These hybrids can potentially infect a wider host range and display enhanced vigour compared to their progenitors. Phytophthora hybrids therefore pose a serious threat to agriculture as well as to natural ecosystems. Early and correct identification of hybrids is therefore essential for adequate plant protection but this is hampered by the limitations of morphological and traditional molecular methods. Identification of hybrids is also important in evolutionary studies as the positioning of hybrids in a phylogenetic tree can lead to suboptimal topologies. To improve the identification of hybrids we have combined genotyping-by-sequencing (GBS) and genome size estimation on a genus-wide collection of 614 Phytophthora isolates. Analyses based on locus- and allele counts and especially on the combination of species-specific loci and genome size estimations allowed us to confirm and characterize 27 previously described hybrid species and discover 16 new hybrid species. Our method was also valuable for species identification at an unprecedented resolution and further allowed correct naming of misidentified isolates. We used both a concatenation- and a coalescent-based phylogenomic method to construct a reliable phylogeny using the GBS data of 140 non-hybrid Phytophthora isolates. Hybrid species were subsequently connected to their progenitors in this phylogenetic tree. In this study we demonstrate the application of two validated techniques (GBS and flow cytometry) for relatively low cost but high resolution identification of hybrids and their phylogenetic relations.

Transcriptional responses in developing lesions of European common ash (Fraxinus excelsior) reveal genes responding to infection by Hymenoscyphus fraxineus.

BACKGROUND: With the expanding ash dieback epidemic that has spread across the European continent, an improved functional understanding of the disease development in afflicted hosts is needed. The study investigated whether differences in necrosis extension between common ash (Fraxinus excelsior) trees with different levels of susceptibility to the fungus Hymenoscyphus fraxineus are associated with, and can be explained by, the differences in gene expression patterns. We inoculated seemingly healthy branches of each of two resistant and susceptible ash genotypes with H. fraxineus grown in a common garden. RESULTS: Ten months after the inoculation, the length of necrosis on the resistant genotypes were shorter than on the susceptible genotypes. RNA sequencing of bark samples collected at the border of necrotic lesions and from healthy tissues distal to the lesion revealed relatively limited differences in gene expression patterns between susceptible and resistant genotypes. At the necrosis front, only 138 transcripts were differentially expressed between the genotype categories while 1082 were differentially expressed in distal, non-symptomatic tissues. Among these differentially expressed genes, several genes in the mevalonate (MVA) and iridoid pathways were found to be co-regulated, possibly indicating increased fluxes through these pathways in response to H. fraxineus. Comparison of transcriptional responses of symptomatic and non-symptomatic ash in a controlled greenhouse experiment revealed a relatively small set of genes that were differentially and concordantly expressed in both studies. This gene-set included the rate-limiting enzyme in the MVA pathway and a number of transcription factors. Furthermore, several of the concordantly expressed candidate genes show significant similarity to genes encoding players in the abscisic acid- or Jasmonate-signalling pathways. CONCLUSIONS: A set of candidate genes, concordantly expressed between field and greenhouse experiments, was identified. The candidates are associated with hormone signalling and specialized metabolite biosynthesis pathways indicating the involvement of these pathways in the response of the host to infection by H. fraxineus.

Combining transcriptomics and genetic linkage based information to identify candidate genes associated with Heterobasidion-resistance in Norway spruce.

The Heterobasidion annosum s.l species complex comprises the most damaging forest pathogens to Norway spruce. We revisited previously identified Quantitative Trait Loci (QTLs) related to Heterobasidion-resistance in Norway spruce to identify candidate genes associated with these QTLs. We identified 329 candidate genes associated with the resistance QTLs using a gene-based composite map for Pinaceae. To evaluate the transcriptional responses of these candidate genes to H. parviporum, we inoculated Norway spruce plants and sequenced the transcriptome of the interaction at 3 and 7 days post inoculation. Out of 298 expressed candidate genes 124 were differentially expressed between inoculation and wounding control treatment. Interestingly, PaNAC04 and two of its paralogs in the subgroup III-3 of the NAC family transcription factors were found to be associated with one of the QTLs and was also highly induced in response to H. parviporum. These genes are possibly involved in the regulation of biosynthesis of flavonoid compounds. Furthermore, several of the differentially expressed candidate genes were associated with the phenylpropanoid pathway including a phenylalanine ammonia-lyase, a cinnamoyl-CoA reductase, a caffeoyl-CoA O-methyltransferase and a PgMYB11-like transcription factor gene. Combining transcriptome and genetic linkage analyses can help identifying candidate genes for functional studies and molecular breeding in non-model species.

Optimized metabarcoding with Pacific biosciences enables semi-quantitative analysis of fungal communities.

Recent studies have questioned the use of high-throughput sequencing of the nuclear ribosomal internal transcribed spacer (ITS) region to derive a semi-quantitative representation of fungal community composition. However, comprehensive studies that quantify biases occurring during PCR and sequencing of ITS amplicons are still lacking. We used artificially assembled communities consisting of 10 ITS-like fragments of varying lengths and guanine-cytosine (GC) contents to evaluate and quantify biases during PCR and sequencing with Illumina MiSeq, PacBio RS II and PacBio Sequel I technologies. Fragment length variation was the main source of bias in observed community composition relative to the template, with longer fragments generally being under-represented for all sequencing platforms. This bias was three times higher for Illumina MiSeq than for PacBio RS II and Sequel I. All 10 fragments in the artificial community were recovered when sequenced with PacBio technologies, whereas the three longest fragments (> 447 bases) were lost when sequenced with Illumina MiSeq. Fragment length bias also increased linearly with increasing number of PCR cycles but could be mitigated by optimization of the PCR setup. No significant biases related to GC content were observed. Despite lower sequencing output, PacBio sequencing was better able to reflect the community composition of the template than Illumina MiSeq sequencing.

Functional Ecology of Forest Disease.

Global change is pressing forest pathologists to solve increasingly complex problems. We argue that understanding interactive effects between forest pathogens and global warming, globalization, and land-use changes may benefit from a functional ecology mindset. Traits can be more informative about ecological functions than species inventories and may deliver a more mechanistic description of forest disease. Myriad microbes with pathogenic potential interact with forest ecosystems at different organizational levels. Elucidation of functional traits may enable the microbial complexity to be reduced into manageable categories with predictive power. In this review, we propose guidelines that allow the research community to develop a functional forest pathology approach. We suggest new angles by which functional questions can be used to resolve burning issues on tree disease. Building up functional databases for pathogenicity is key to implementing these approaches.

Association genetics identifies a specifically regulated Norway spruce laccase gene, PaLAC5, linked to Heterobasidion parviporum resistance.

It is important to improve the understanding of the interactions between the trees and pathogens and integrate this knowledge about disease resistance into tree breeding programs. The conifer Norway spruce (Picea abies) is an important species for the forest industry in Europe. Its major pathogen is Heterobasidion parviporum, causing stem and root rot. In this study, we identified 11 Norway spruce QTLs (Quantitative trait loci) that correlate with variation in resistance to H. parviporum in a population of 466 trees by association genetics. Individual QTLs explained between 2.1 and 5.2% of the phenotypic variance. The expression of candidate genes associated with the QTLs was analysed in silico and in response to H. parviporum hypothesizing that (a) candidate genes linked to control of fungal sapwood growth are more commonly expressed in sapwood, and; (b) candidate genes associated with induced defences are respond to H. parviporum inoculation. The Norway spruce laccase PaLAC5 associated with control of lesion length development is likely to be involved in the induced defences. Expression analyses showed that PaLAC5 responds specifically and strongly in close proximity to the H. parviporum inoculation. Thus, PaLAC5 may be associated with the lignosuberized boundary zone formation in bark adjacent to the inoculation site.

The conifer root rot pathogens Heterobasidion irregulare and Heterobasidion occidentale employ different strategies to infect Norway spruce.

Heterobasidion irregulare and H. occidentale are two closely related conifer root rot pathogens in the H. annosum sensu lato (s.l.) species complex. The two species H. irregulare and H. occidentale have different host preference with pine and non-pine tree species favored, respectively. The comparison of transcriptomes of H. irregulare and H. occidentale growing in Norway spruce bark, a susceptible host non-native to North America, showed large differences in gene expression. Heterobasidion irregulare induced more genes involved in detoxification of host compounds and in production of secondary metabolites, while the transcriptome induced in H. occidentale was more oriented towards carbohydrate degradation. Along with their separated evolutionary history, the difference might be driven by their host preferences as indicated by the differentially expressed genes enriched in particular Gene Ontology terms.

Occurrence and pathogenicity of Corinectria spp. - an emerging canker disease of Abies sibirica in Central Siberia.

During recent years, a new disease of Siberian fir (A. sibirica) emerged in Central Siberia, exhibiting symptoms of stem/branch deformation, cambium necrosis, and dieback of branches and twigs, the causal agent remaining unknown. The aim was to identify agent of the disease and to investigate its pathogenicity to A. sibirica and Norway spruce (Picea abies). Symptomatic tissues of fir were subjected to pure culture isolation of anticipated pathogen(s). Obtained isolates were subjected to molecular identification, phylogenetic analyses, and pathogenicity tests with A. sibirica saplings, and seeds and seedlings of A. sibirica and P. abies. The study demonstrated that, (i) most commonly isolated fungus from canker wounds of A. sibirica exhibited Acremonium-like anamorphs; (ii) phylogeny demonstrated that investigated fungi belong to genus Corinectria, but are genetically well separated from other worldwide known Corinectria spp.; (iii) one species of isolated fungi has the capacity to cause the disease and kill A. sibirica saplings and seedlings, but also seedlings of P. abies. Guidelines for future research were defined in order to generate needed information on species description, its origin and ecology, and estimation of potential risks upon the eventual invasion of the pathogen to new geographic areas, in particular of Europe.

Genetic Variation Explains Changes in Susceptibility in a Naïve Host Against an Invasive Forest Pathogen: The Case of Alder and the Phytophthora alni Complex.

Predicting whether naïve tree populations have the potential to adapt to exotic pathogens is necessary owing to the increasing rate of invasions. Adaptation may occur as a result of natural selection when heritable variation in terms of susceptibility exists in the naïve population. We searched for signs of selection on black alder (Alnus glutinosa) stands growing on riverbanks invaded by two pathogens differing in aggressiveness, namely, Phytophthora uniformis (PU) and Phytophthora × alni (PA). We compared the survival and heritability measures from 72 families originating from six invaded and uninvaded (naïve) sites by performing in vitro inoculations. The results from the inoculations were used to assess the relative contribution of host genetic variation on natural selection. We found putative signs of natural selection on alder exerted by PU but not by PA. For PU, we found a higher survival in families originating from invaded sites compared with uninvaded sites. The narrow sense heritability of susceptibility to PU of uninvaded populations was significantly higher than to PA. Simulated data supported the role of heritable genetic variation on natural selection and discarded a high aggressiveness of PA decreasing the transmission rate as an alternative hypothesis for a slow natural selection. Our findings expand on previous attempts of using heritability as a predictor for the likelihood of natural adaptation of naïve tree populations to invasive pathogens. Measures of genetic variation can be useful for risk assessment purposes or when managing Phytophthora invasions.

Genotypic variation in Norway spruce correlates to fungal communities in vegetative buds.

The taxonomically diverse phyllosphere fungi inhabit leaves of plants. Thus, apart from the fungi's dispersal capacities and environmental factors, the assembly of the phyllosphere community associated with a given host plant depends on factors encoded by the host's genome. The host genetic factors and their influence on the assembly of phyllosphere communities under natural conditions are poorly understood, especially in trees. Recent work indicates that Norway spruce (Picea abies) vegetative buds harbour active fungal communities, but these are hitherto largely uncharacterized. This study combines internal transcribed spacer sequencing of the fungal communities associated with dormant vegetative buds with a genome-wide association study (GWAS) in 478 unrelated Norway spruce trees. The aim was to detect host loci associated with variation in the fungal communities across the population, and to identify loci correlating with the presence of specific, latent, pathogens. The fungal communities were dominated by known Norway spruce phyllosphere endophytes and pathogens. We identified six quantitative trait loci (QTLs) associated with the relative abundance of the dominating taxa (i.e., top 1% most abundant taxa). Three additional QTLs associated with colonization by the spruce needle cast pathogen Lirula macrospora or the cherry spruce rust (Thekopsora areolata) in asymptomatic tissues were detected. The identification of the nine QTLs shows that the genetic variation in Norway spruce influences the fungal community in dormant buds and that mechanisms underlying the assembly of the communities and the colonization of latent pathogens in trees may be uncovered by combining molecular identification of fungi with GWAS.

Invasive forest pathogens in Europe: Cross-country variation in public awareness but consistency in policy acceptability.

Political action can reduce introductions of diseases caused by invasive forest pathogens (IPs) and public support is important for effective prevention. The public's awareness of IP problems and the acceptability of policies aiming to combat these pathogens were surveyed in nine European countries (N = 3469). Although awareness of specific diseases (e.g., ash dieback) varied, problem awareness and policy acceptability were similar across countries. The public was positive towards policies for informational measures and stricter standards for plant production, but less positive towards restricting public access to protected areas. Multilevel models, including individual and country level variables, revealed that media exposure was positively associated with awareness of IP problems, and strengthened the link between problem awareness and policy acceptability. Results suggest that learning about IPs through the media and recognizing the associated problems increase policy acceptability. Overall, the study elaborates on the anthropogenic dimension of diseases caused by IPs.

Contrasting distribution patterns between aquatic and terrestrial Phytophthora species along a climatic gradient are linked to functional traits.

Diversity of microbial organisms is linked to global climatic gradients. The genus Phytophthora includes both aquatic and terrestrial plant pathogenic species that display a large variation of functional traits. The extent to which the physical environment (water or soil) modulates the interaction of microorganisms with climate is unknown. Here, we explored the main environmental drivers of diversity and functional trait composition of Phytophthora communities. Communities were obtained by a novel metabarcoding setup based on PacBio sequencing of river filtrates in 96 river sites along a geographical gradient. Species were classified as terrestrial or aquatic based on their phylogenetic clade. Overall, terrestrial and aquatic species showed contrasting patterns of diversity. For terrestrial species, precipitation was a stronger driver than temperature, and diversity and functional diversity decreased with decreasing temperature and precipitation. In cold and dry areas, the dominant species formed resistant structures and had a low optimum temperature. By contrast, for aquatic species, temperature and water chemistry were the strongest drivers, and diversity increased with decreasing temperature and precipitation. Within the same area, environmental filtering affected terrestrial species more strongly than aquatic species (20% versus 3% of the studied communities, respectively). Our results highlight the importance of functional traits and the physical environment in which microorganisms develop their life cycle when predicting their distribution under changing climatic conditions. Temperature and rainfall may be buffered differently by water and soil, and thus pose contrasting constrains to microbial assemblies.