Arctic root‐associated fungal community composition reflects environmental filtering

There is growing evidence that root‐associated fungi have important roles in Arctic ecosystems. Here, we assess the diversity of fungal communities associated with roots of the ectomycorrhizal perennial herb Bistorta vivipara on the Arctic archipelago of Svalbard and investigate whether spatial separation and bioclimatic variation are important structuring factors of fungal community composition. We sampled 160 plants of B. vivipara from 32 localities across Svalbard. DNA was extracted from entire root systems, and 454 pyrosequencing of ITS1 amplicons was used to profile the fungal communities. The fungal communities were predominantly composed of Basidiomycota (55% of reads) and Ascomycota (35%), with the orders Thelephorales (24%), Agaricales (13.8%), Pezizales (12.6%) and Sebacinales (11.3%) accounting for most of the reads. Plants from the same site or region had more similar fungal communities to one another than plants from other sites or regions, and sites clustered together along a weak latitudinal gradient. Furthermore, a decrease in per‐plant OTU richness with increasing latitude was observed. However, no statistically significant spatial autocorrelation between sites was detected, suggesting that environmental filtering, not dispersal limitation, causes the observed patterns. Our analyses suggest that while latitudinal patterns in community composition and richness might reflect bioclimatic influences at global spatial scales, at the smaller spatial scale of the Svalbard archipelago, these changes more likely reflect varied bedrock composition and associated edaphic factors. The need for further studies focusing on identifying those specific bioclimatic and edaphic factors structuring root‐associated fungal community composition at both global and local scales is emphasized.     

Experimental warming increases herbivory by leaf‐chewing insects in an alpine plant community

Climate warming is predicted to affect species and trophic interactions worldwide, and alpine ecosystems are expected to be especially sensitive to changes. In this study, we used two ongoing climate warming (open‐top chambers) experiments at Finse, southern Norway, to examine whether warming had an effect on herbivory by leaf‐chewing insects in an alpine Dryas heath community. We recorded feeding marks on the most common vascular plant species in warmed and control plots at two experimental sites at different elevations and carried out a brief inventory of insect herbivores. Experimental warming increased herbivory on Dryas octopetala and Bistorta vivipara. Dryas octopetala also experienced increased herbivory at the lower and warmer site, indicating an overall positive effect of warming, whereas B. vivipara experienced an increased herbivory at the colder and higher site indicating a mixed effect of warming. The Lepidoptera Zygaena exulans and Sympistis nigrita were the two most common leaf‐chewing insects in the Dryas heath. Based on the observed patterns of herbivory, the insects life cycles and feeding preferences, we argue that Z. exulans is the most important herbivore on B. vivipara, and S. nigrita the most important herbivore on D. octopetala. We conclude that if the degree of insect herbivory increases in a warmer world, as suggested by this study and others, complex interactions between plants, insects, and site‐specific conditions make it hard to predict overall effects on plant communities.     

Temporal variation of Bistorta vivipara‐associated ectomycorrhizal fungal communities in the High Arctic

Ectomycorrhizal (ECM) fungi are important for efficient nutrient uptake of several widespread arctic plant species. Knowledge of temporal variation of ECM fungi, and the relationship of these patterns to environmental variables, is essential to understand energy and nutrient cycling in Arctic ecosystems. We sampled roots of Bistorta vivipara ten times over two years; three times during the growing‐season (June, July and September) and twice during winter (November and April) of both years. We found 668 ECM OTUs belonging to 25 different ECM lineages, whereof 157 OTUs persisted throughout all sampling time‐points. Overall, ECM fungal richness peaked in winter and species belonging to Cortinarius, Serendipita and Sebacina were more frequent in winter than during summer. Structure of ECM fungal communities was primarily affected by spatial factors. However, after accounting for spatial effects, significant seasonal variation was evident revealing correspondence with seasonal changes in environmental conditions. We demonstrate that arctic ECM richness and community structure differ between summer (growing‐season) and winter, possibly due to reduced activity of the core community, and addition of fungi adapted for winter conditions forming a winter‐active fungal community. Significant month × year interactions were observed both for fungal richness and community composition, indicating unpredictable between‐year variation. Our study indicates that addressing seasonal changes requires replication over several years.     

A genome scan of diversifying selection in Ophiocordyceps zombie‐ant fungi suggests a role for enterotoxins in co‐evolution and host specificity

Identification of the genes underlying adaptation sheds light on the biological functions targeted by natural selection. Searches for footprints of positive selection, in the form of rapid amino acid substitutions, and the identification of species‐specific genes have proved to be powerful approaches to identifying the genes involved in host specialization in plant‐pathogenic fungi. We used an evolutionary comparative genomic approach to identify genes underlying host adaptation in the ant‐infecting genus Ophiocordyceps, which manipulates ant behaviour. A comparison of the predicted genes in the genomes of species from three species complexes—O. unilateralis, O. australis and O. subramanianii—revealed an enrichment in pathogenesis‐associated functions, including heat‐labile enterotoxins, among species‐specific genes. Furthermore, these genes were overrepresented among those displaying significant footprints of positive selection. Other categories of genes suspected to be important for virulence and pathogenicity in entomopathogenic fungi (e.g., chitinases, lipases, proteases, core secondary metabolism genes) were much less represented, although a few candidate genes were found to evolve under positive selection. An analysis including orthologs from other entomopathogenic fungi in a broader context showed that positive selection on enterotoxins was specific to the ant‐infecting genus Ophiocordyceps. Together with previous studies reporting the overexpression of an enterotoxin during behavioural manipulation in diseased ants, our findings suggest that heat‐labile enterotoxins are important effectors in host adaptation and co‐evolution in the Ophiocordyceps entomopathogenic fungi.     

Aphid–willow interactions in a high Arctic ecosystem: responses to raised temperature and goose disturbance

Recently, there have been several studies using open top chambers (OTCs) or cloches to examine the response of Arctic plant communities to artificially elevated temperatures. Few, however, have investigated multitrophic systems, or the effects of both temperature and vertebrate grazing treatments on invertebrates. This study investigated trophic interactions between an herbivorous insect (Sitobion calvulum, Aphididae), a woody perennial host plant (Salix polaris) and a selective vertebrate grazer (barnacle geese, Branta leucopsis). In a factorial experiment, the responses of the insect and its host to elevated temperatures using open top chambers (OTCs) and to three levels of goose grazing pressure were assessed over two summer growing seasons (2004 and 2005). OTCs significantly enhanced the leaf phenology of Salix in both years and there was a significant OTC by goose presence interaction in 2004. Salix leaf number was unaffected by treatments in both years, but OTCs increased leaf size and mass in 2005. Salix reproduction and the phenology of flowers were unaffected by both treatments. Aphid densities were increased by OTCs but unaffected by goose presence in both years. While goose presence had little effect on aphid density or host plant phenology in this system, the OTC effects provide interesting insights into the possibility of phenological synchrony disruption. The advanced phenology of Salix effectively lengthens the growing season for the plant, but despite a close association with leaf maturity, the population dynamics of the aphid appeared to lack a similar phenological response, except for the increased population observed.     

Host specificity of two pollinating seed‐consuming fly species is not related to soil moisture of host plant in the high Himalayas

Studying the drivers of host specificity can contribute to our understanding of the origin and evolution of obligate pollination mutualisms. The preference–performance hypothesis predicts that host plant choice of female insects is related mainly to the performance of their offspring. Soil moisture is thought to be particularly important for the survival of larvae and pupae that inhabit soil. In the high Himalayas, Rheum nobile and R. alexandrae differ in their distribution in terms of soil moisture; that is, R. nobile typically occurs in scree with well‐drained soils, R. alexandrae in wetlands. The two plant species are pollinated by their respective mutualistic seed‐consuming flies, Bradysia sp1. and Bradysia sp2. We investigated whether soil moisture is important for regulating host specificity by comparing pupation and adult emergence of the two fly species using field and laboratory experiments. Laboratory experiments revealed soil moisture did have significant effects on larval and pupal performances in both fly species, but the two fly species had similar optimal soil moisture requirements for pupation and adult emergence. Moreover, a field reciprocal transfer experiment showed that there was no significant difference in adult emergence for both fly species between their native and non‐native habitats. Nevertheless, Bradysia sp1., associated with R. nobile, was more tolerant to drought stress, while Bradysia sp2., associated with R. alexandrae, was more tolerant to flooding stress. These results indicate that soil moisture is unlikely to play a determining role in regulating host specificity of the two fly species. However, their pupation and adult emergence in response to extremely wet or dry soils are habitat‐specific.     

Neighboring trees regulate the root‐associated pathogenic fungi on the host plant in a subtropical forest

Root‐associated fungi and host‐specific pathogens are major determinants of species coexistence in forests. Phylogenetically related neighboring trees can strongly affect the fungal community structure of the host plant, which, in turn, will affect the ecological processes. Unfortunately, our understanding of the factors influencing fungal community composition in forests is still limited. In particular, investigation of the relationship between the phytopathogenic fungal community and neighboring trees is incomplete. In the current study, we tested the host specificity of members of the root‐associated fungal community collected from seven tree species and determined the influence of neighboring trees and habitat variation on the composition of the phytopathogenic fungal community of the focal plant in a subtropical evergreen forest. Using high‐throughput sequencing data with respect to the internal transcribed spacer (ITS) region, we characterized the community composition of the root‐associated fungi and found significant differences with respect to fungal groups among the seven tree species. The density of conspecific neighboring trees had a significantly positive influence on the relative abundance of phytopathogens, especially host‐specific pathogens, while the heterospecific neighbor density had a significant negative impact on the species richness of host‐specific pathogens, as well as phytopathogens. Our work provides evidence that the root‐associated phytopathogenic fungi of a host plant depend greatly on the tree neighbors of the host plant.     

Divergence in gene expression related to variation in host specificity of an ectomycorrhizal fungus

Ectomycorrhizae are formed by mutualistic interactions between fungi and the roots of woody plants. During symbiosis the two organisms exchange carbon and nutrients in a specific tissue that is formed at the contact between a compatible fungus and plant. There is considerable variation in the degree of host specificity among species and strains of ectomycorrhizal fungi. In this study, we have for the first time shown that this variation is associated with quantitative differences in gene expression, and with divergence in nucleotide sequences of symbiosis-regulated genes. Gene expression and sequence evolution were compared in different strains of the ectomycorrhizal fungus Paxillus involutus; the strains included Nau, which is not compatible with birch and poplar, and the two compatible strains Maj and ATCC200175. On a genomic level, Nau and Maj were very similar. The sequence identity was 98.9% in the 16 loci analysed, and only three out of 1075 genes analysed by microarray-based hybridizations had signals indicating differences in gene copy numbers. In contrast, 66 out of the 1075 genes were differentially expressed in Maj compared to Nau after contact with birch roots. Thirty-seven of these symbiosis-regulated genes were also differentially expressed in the ATCC strain. Comparative analysis of DNA sequences of the symbiosis-regulated genes in different strains showed that two of them have evolved at an enhanced rate in Nau. The sequence divergence can be explained by a decreased selection pressure, which in turn is determined by lower functional constraints on these proteins in Nau as compared to the compatible strains.     

Substantial compositional turnover of fungal communities in an alpine ridge‐to‐snowbed gradient

The main gradient in vascular plant, bryophyte and lichen species composition in alpine areas, structured by the topographic gradient from wind‐exposed ridges to snowbeds, has been extensively studied. Tolerance to environmental stress, resulting from wind abrasion and desiccation towards windswept ridges or reduced growing season due to prolonged snow cover towards snowbeds, is an important ecological mechanism in this gradient. The extent to which belowground fungal communities are structured by the same topographic gradient and the eventual mechanisms involved are less well known. In this study, we analysed variation in fungal diversity and community composition associated with roots of the ectomycorrhizal plant Bistorta vivipara along the ridge‐to‐snowbed gradient. We collected root samples from fifty B. vivipara plants in ten plots in an alpine area in central Norway. The fungal communities were analysed using 454 pyrosequencing analyses of tag‐encoded ITS1 amplicons. A distinct gradient in the fungal community composition was found that coincided with variation from ridge to snowbeds. This gradient was paralleled by change in soil content of carbon, nitrogen and phosphorus. A large proportion (66%) of the detected 801 nonsingleton operational taxonomic units (OTUs) were ascomycetes, while basidiomycetes dominated quantitatively (i.e. with respect to number of reads). Numerous fungal OTUs, many with taxonomic affinity to Sebacinales, Cortinarius and Meliniomyces, showed distinct affinities either to ridge or to snowbed plots, indicating habitat specialization. The compositional turnover of fungal communities along the gradient was not paralleled by a gradient in species richness.     

Community composition and diversity of Neotropical root‐associated fungi in common and rare trees

Interactions between plants and root‐associated fungi can affect the assembly, diversity, and relative abundances of tropical plant species. Host–symbiont compatibility and some degree of host specificity are prerequisites for these processes to occur, and these prerequisites may vary with host abundance. However, direct assessments of whether specificity of root‐associated fungi varies with host abundance are lacking. Here, in a diverse tropical forest in Los Tuxtlas, Mexico, we couple DNA metabarcoding with a sampling design that controls for host phylogeny, host age, and habitat variation, to characterize fungal communities associated with the roots of three confamilial pairs of host species that exhibit contrasting (high and low) relative abundances. We uncovered a functionally and phylogenetically diverse fungal community composed of 1,038 OTUs (operational taxonomic units with 97% genetic similarity), only 14 of which exhibited host specificity. Host species was a significant predictor of fungal community composition only for the subset of OTUs composed of putatively pathogenic fungi. We found no significant difference in the number of specialists associating with common versus rare trees, but we found that host abundance was negatively correlated with the diversity of root fungal communities. This latter result was significant for symbiotrophs (mostly arbuscular mycorrhizal fungi) and, to a lesser extent, for pathotrophs (mostly plant pathogens). Thus, root fungal communities differ between common and rare trees, which may impact the strength of conspecific negative density dependence. Further studies from other tropical sites and host lineages are warranted, given the role of root‐associated fungi in biodiversity maintenance.     

Contrasting the seasonal and elevational prevalence of generalist avian haemosporidia in co‐occurring host species

Understanding the ecology and evolution of parasites is contingent on identifying the selection pressures they face across their infection landscape. Such a task is made challenging by the fact that these pressures will likely vary across time and space, as a result of seasonal and geographical differences in host susceptibility or transmission opportunities. Avian haemosporidian blood parasites are capable of infecting multiple co‐occurring hosts within their ranges, yet whether their distribution across time and space varies similarly in their different host species remains unclear. Here, we applied a new PCR method to detect avian haemosporidia (genera Haemoproteus, Leucocytozoon, and Plasmodium) and to determine parasite prevalence in two closely related and co‐occurring host species, blue tits (Cyanistes caeruleus, N = 529) and great tits (Parus major, N = 443). Our samples were collected between autumn and spring, along an elevational gradient in the French Pyrenees and over a three‐year period. Most parasites were found to infect both host species, and while these generalist parasites displayed similar elevational patterns of prevalence in the two host species, this was not always the case for seasonal prevalence patterns. For example, Leucocytozoon group A parasites showed inverse seasonal prevalence when comparing between the two host species, being highest in winter and spring in blue tits but higher in autumn in great tits. While Plasmodium relictum prevalence was overall lower in spring relative to winter or autumn in both species, spring prevalence was also lower in blue tits than in great tits. Together, these results reveal how generalist parasites can exhibit host‐specific epidemiology, which is likely to complicate predictions of host–parasite co‐evolution.     

Effects of host genetics and environment on egg‐associated microbiotas in brown trout (Salmo trutta)

Recent studies found fish egg‐specific bacterial communities that changed over the course of embryogenesis, suggesting an interaction between the developing host and its microbiota. Indeed, single‐strain infections demonstrated that the virulence of opportunistic bacteria is influenced by environmental factors and host immune genes. However, the interplay between a fish embryo host and its microbiota has not been studied yet at the community level. To test whether host genetics affects the assemblage of egg‐associated bacteria, adult brown trout (Salmo trutta) were sampled from a natural population. Their gametes were used for full‐factorial in vitro fertilizations to separate sire from dam effects. In total, 2520 embryos were singly raised under experimental conditions that differently support microbial growth. High‐throughput 16S rRNA amplicon sequencing was applied to characterize bacterial communities on milt and fertilized eggs across treatments. Dam and sire identity influenced embryo mortality, time until hatching and composition of egg‐associated microbiotas, but no link between bacterial communities on milt and on fertilized eggs could be found. Elevated resources increased embryo mortality and modified bacterial communities with a shift in their putative functional potential. Resource availability did not significantly affect any parental effects on embryo performance. Sire identity affected bacterial diversity that turned out to be a significant predictor of hatching time: embryos associated with high bacterial diversity hatched later. We conclude that both host genetics and the availability of resources define diversity and composition of egg‐associated bacterial communities that then affect the life history of their hosts.     

Amplicon‐pyrosequencing‐based detection of compositional shifts in bryophyte‐associated fungal communities along an elevation gradient

Although bryophytes are a dominant vegetation component of boreal and alpine ecosystems, little is known about their associated fungal communities. HPLC assays of ergosterol (fungal biomass) and amplicon pyrosequencing of the ITS2 region of rDNA were used to investigate how the fungal communities associated with four bryophyte species changed across an elevational gradient transitioning from conifer forest to the low‐alpine. Fungal biomass and OTU richness associated with the four moss hosts did not vary significantly across the gradient (P > 0.05), and both were more strongly affected by host and tissue type. Despite largely constant levels of fungal biomass, distinct shifts in community composition of fungi associated with Hylocomium, Pleurozium and Polytrichum occurred between the elevation zones of the gradient. This likely is a result of influence on fungal communities by major environmental factors such as temperature, directly or indirectly mediated by, or interacting with, the response of other components of the vegetation (i.e. the dominant trees). Fungal communities associated with Dicranum were an exception, exhibiting spatial autocorrelation between plots, and no significant structuring by elevation. Nevertheless, the detection of distinct fungal assemblages associated with a single host growing in different elevation zones along an elevational gradient is of particular relevance in the light of the ongoing changes in vegetation patterns in boreal and alpine systems due to global climate warming.     

Cryptic diversity,high host specificity and reproductive synchronization in army ant‐associated Vatesus beetles

Army ants and their arthropod symbionts represent one of the most species‐rich animal associations on Earth, and constitute a fascinating example of diverse host–symbiont interaction networks. However, despite decades of research, our knowledge of army ant symbionts remains fragmentary due to taxonomic ambiguity and the inability to study army ants in the laboratory. Here, we present an integrative approach that allows us to reliably determine species boundaries, assess biodiversity, match different developmental stages and sexes, and to study the life cycles of army ant symbionts. This approach is based on a combination of community sampling, DNA barcoding, morphology and physiology. As a test case, we applied this approach to the staphylinid beetle genus Vatesus and its different Eciton army ant host species at La Selva Biological Station, Costa Rica. DNA barcoding led to the discovery of cryptic biodiversity and, in combination with extensive community sampling, revealed strict host partitioning with no overlap in host range. Using DNA barcoding, we were also able to match the larval stages of all focal Vatesus species. In combination with studies of female reproductive physiology, this allowed us to reconstruct almost the complete life cycles of the different beetle species. We show that Vatesus beetles are highly adapted to the symbiosis with army ants, in that their reproduction and larval development are synchronized with the stereotypical reproductive and behavioural cycles of their host colonies. Our approach can now be used to study army ant‐symbiont communities more broadly, and to obtain novel insights into co‐evolutionary and ecological dynamics in species‐rich host–symbiont systems.     

Seasonal variation in epifaunal communities associated with giant kelp (Macrocystis pyrifera) at an upwelling‐dominated site

Kelp forests are highly productive and species‐rich benthic ecosystems in temperate regions that provide biogenic habitat for numerous associated species. Diverse epifaunal communities inhabit kelp sporophytes and are subject to variations in the physical environment and to changes experienced by the kelp habitat itself. We assessed seasonal variations in epifaunal invertebrate communities inhabiting giant kelps, Macrocystis pyrifera, and their effects on this seaweed. Six seasonal samplings were conducted over a year at an upwelling‐dominated site in northern‐central Chile where physical conditions are known to fluctuate temporally. More than 30 taxa were identified, among which peracarid crustaceans stood out in both diversity and abundance. Species richness and abundance differed among sporophyte sections (holdfast and fronds) and throughout the year. The frond community was dominated by two grazers (the amphipod Peramphithoe femorata and the isopod Amphoroidea typa), while suspension feeders, grazers, and omnivores (the amphipod Aora typica, the isopod Limnoria quadripunctata, and polychaetes) dominated the holdfasts. Abundances of the dominant species fluctuated throughout the year but patterns of variation differed among species. The most abundant grazer (P. femorata) had highest densities in summer, while the less abundant grazer (A. typa) reached its peak densities in winter. Interestingly, the area of kelp damaged by grazers was highest in autumn and early winter, suggesting that grazing impacts accumulate during periods of low kelp growth, which can thus be considered as ‘vestiges of herbivory past.’ Among the factors determining the observed seasonal patterns, strong variability of environmental conditions, reproductive cycles of associated fauna, and predation by fishes vary in importance. Our results suggest that during spring and early summer, bottom‐up processes shape the community structure of organisms inhabiting large perennial seaweeds, whereas during late summer and autumn, top‐down processes are more important.