Bacterial communities associated with an island radiation of lichen-forming fungi.

Evolutionary radiations are one of the most striking processes biologists have studied in islands. A radiation is often sparked by the appearance of ecological opportunity, which can originate in processes like trophic niche segregation or the evolution of key innovations. Another recently proposed mechanism is facilitation mediated by the bacterial communities associated with the radiating species. Here we explore the role of the bacterial communities in a radiation of lichen-forming fungi endemic to Macaronesia. Bacterial diversity was quantified by high throughput sequencing of the V1-V2 hyper-variable region of 172 specimens. We characterized the taxonomic and phylogenetic diversity of the bacterial communities associated with the different species, tested for compositional differences between these communities, carried out a functional prediction, explored the relative importance of different factors in bacterial community structure, searched for phylosymbiosis and tried to identify the origin of this pattern. The species of the radiation differed in the composition of their bacterial communities, which were mostly comprised of Alphaproteobacteria and Acidobacteriia, but not in the functionality of those communities. A phylosimbiotic pattern was detected, but it was probably caused by environmental filtering. These findings are congruent with the combined effect of secondary chemistry and mycobiont identity being the main driver of bacterial community structure. Altogether, our results suggest that the associated bacterial communities are not the radiation's main driver. There is one possible exception, however, a species that has an abnormally diverse core microbiome and whose bacterial communities could be subject to a specific environmental filter at the functional level.

Ocean crossers: A tale of disjunctions and speciation in the dwarf-fruticose Lichina (lichenized Ascomycota).

Lichens thrive in rocky coastal areas in temperate and cold regions of both hemispheres. Species of the genus Lichina, which form characteristic black fruiting thalli associated with cyanobacteria, often create distinguishable bands in the intertidal and supralittoral zones. The present study uses a comprehensive specimen dataset and four gene loci to (1) delineate and discuss species boundaries in this genus, (2) assess evolutionary relationships among species, and (3) infer the most likely causes of their current geographic distribution in the Northern and Southern hemispheres. A dated phylogeny describes the time frame in which extant disjunctions of species and populations were established. The results showed that the genus is integrated by four species, with Lichina pygmaea, L. confinis and the newly described L. canariensis from rocky seashores in the Canary Islands, occurring in the Northern Hemisphere, whereas L. intermedia is restricted to the Southern Hemisphere. Lichina intermedia hosted a much higher intraspecific genetic diversity than the other species, with subclades interpreted as species-level lineages by the different species delimitation approaches. However, a conservative taxonomic approach was adopted. This species showed a striking disjunct distribution between Australasia and southern South America. The timing for the observed interspecific and intraspecific divergences and population disjunctions postdated continental plate movements, suggesting that long-distance dispersal across body waters in the two hemispheres played a major role in shaping the current species distributions. Such ocean crossings were, as in L. canariensis, followed by speciation. New substitution rates for the nrITS of the genus Lichina were inferred using a tree spanning the major Ascomycota lineages calibrated using fossils. In conclusion, this work lays the foundation for a better understanding of the evolution through time and space of maritime lichens.

Biological Microbial Interactions from Cooccurrence Networks in a High Mountain Lacustrine District.

A fundamental question in biology is why some species tend to occur together in the same locations, while others are never observed coexisting. This question becomes particularly relevant for microorganisms thriving in the highly diluted waters of high mountain lakes, where biotic interactions might be required to make the most of an extreme environment. We studied a high-throughput gene data set of alpine lakes (>220 Pyrenean lakes) with cooccurrence network analysis to infer potential biotic interactions, using the combination of a probabilistic method for determining significant cooccurrences and coexclusions between pairs of species and a conceptual framework for classifying the nature of the observed cooccurrences and coexclusions. This computational approach (i) determined and quantified the importance of environmental variables and spatial distribution and (ii) defined potential interacting microbial assemblages. We determined the properties and relationships between these assemblages by examining node properties at the taxonomic level, indicating associations with their potential habitat sources (i.e., aquatic versus terrestrial) and their functional strategies (i.e., parasitic versus mixotrophic). Environmental variables explained fewer pairs in bacteria than in microbial eukaryotes for the alpine data set, with pH alone explaining the highest proportion of bacterial pairs. Nutrient composition was also relevant for explaining association pairs, particularly in microeukaryotes. We identified a reduced subset of pairs with the highest probability of species interactions ("interacting guilds") that significantly reached higher occupancies and lower mean relative abundances in agreement with the carrying capacity hypothesis. The interacting bacterial guilds could be more related to habitat and microdispersal processes (i.e., aquatic versus soil microbes), whereas for microeukaryotes trophic roles (osmotrophs, mixotrophs, and parasitics) could potentially play a major role. Overall, our approach may add helpful information to guide further efforts for a mechanistic understanding of microbial interactions in situ. IMPORTANCE A fundamental question in biology is why some species tend to occur together in the same locations, while others are never observed to coexist. This question becomes particularly relevant for microorganisms thriving in the highly diluted waters of high mountain lakes, in which biotic interactions might be required to make the most of an extreme environment. Microbial metacommunities are too often only studied in terms of their environmental niches and geographic barriers since they show inherent difficulties to quantify biological interactions and their role as drivers of ecosystem functioning. Our study highlights that telling apart potential interactions from both environmental and geographic niches may help for the initial characterization of organisms with similar ecologies in a large scope of ecosystems, even when information about actual interactions is partial and limited. The multilayered statistical approach carried out here offers the possibility of going beyond taxonomy to understand microbiological behavior in situ.

A global microbiome survey of vineyard soils highlights the microbial dimension of viticultural terroirs.

The microbial biodiversity found in different vitivinicultural regions is an important determinant of wine terroir. It should be studied and preserved, although it may, in the future, be subjected to manipulation by precision agriculture and oenology. Here, we conducted a global survey of vineyards' soil microbial communities. We analysed soil samples from 200 vineyards on four continents to establish the basis for the development of a vineyard soil microbiome's map, representing microbial biogeographical patterns on a global scale. This study describes vineyard microbial communities worldwide and establishes links between vineyard locations and microbial biodiversity on different scales: between continents, countries, and between different regions within the same country. Climate data correlates with fungal alpha diversity but not with prokaryotes alpha diversity, while spatial distance, on a global and national scale, is the main variable explaining beta-diversity in fungal and prokaryotes communities. Proteobacteria, Actinobacteria and Acidobacteria phyla, and Archaea genus Nitrososphaera dominate prokaryotic communities in soil samples while the overall fungal community is dominated by the genera Solicoccozyma, Mortierella and Alternaria. Finally, we used microbiome data to develop a predictive model based on random forest analyses to discriminate between microbial patterns and to predict the geographical source of the samples with reasonable precision.

Network Properties of Local Fungal Communities Reveal the Anthropogenic Disturbance Consequences of Farming Practices in Vineyard Soils.

Agroecosystems are human-managed ecosystems subject to generalized ecological rules. Understanding the ecology behind the assembly and dynamics of soil fungal communities is a fruitful way to improve management practices and plant productivity. Thus, monitoring soil health would benefit from the use of metrics that arise from ecological explanations that can also be informative for agricultural management. Beyond traditional biodiversity descriptors, community-level properties have the potential of informing about particular ecological situations. Here we assess the impact of different farming practices in a survey of 350 vineyard soils from the United States and Spain by estimating network properties based on spatial associations. Our observations using traditional approaches show results concurring with previous literature: the influence of geographic and climatic factors on sample distributions, or different operational taxonomic unit (OTU) compositions depending on agricultural managements. Furthermore, using network properties, we observe that fungal communities ranged from dense arrangements of associations to a sparser structure of associations, indicating differential levels of niche specialization. We detect fungal arrangements capable of thriving in wider or smaller ranges of temperature, revealing that niche specialization may be a critical soil process impacting soil health. Low-intervention practices (organic and biodynamic managements) promoted densely clustered networks, describing an equilibrium state based on mixed collaborative communities. In contrast, conventionally managed vineyards had highly modular sparser communities, supported by a higher coexclusion proportion. Thus, we hypothesize that network properties at the community level may help to understand how the environment and land use can affect community structure and ecological processes in agroecosystems.IMPORTANCE Soil fungal communities play a key role in agroecosystem sustainability. The complexity of fungal communities, at both taxonomic and functional levels, makes it difficult to find clear patterns connecting community composition with ecosystem function and to understand the impact of biotic (interspecies interactions) and abiotic (e.g., climate or anthropogenic disturbances) factors on it. Here we combine network analysis methods and properties, proposing a novel analytical approach: to infer ecological properties from local networks, which we apply to the study of fungal communities in vineyard soils. We conclude that different levels of farming intensification may lead to different ecological strategies in soil fungal communities settled by particular association arrangements.

Tuning up microbiome analysis to monitor WWTPs' biological reactors functioning.

Wastewater treatment plants (WWTPs) are necessary to protect ecosystems quality and human health. Their function relies on the degradation of organic matter and nutrients from a water influent, prior to the effluent release into the environment. In this work we studied the bacterial community dynamics of a municipal WWTP with a membrane bioreactor through 16S rRNA gene sequencing. The main phyla identified in the wastewater were Proteobacteria, Bacteroidetes, Chloroflexi, Planctomycetes and Actinobacteria. The WWTP is located in Spain and, like other studied WWTP in temperate climate zones, the temperature played a major role in community assembly. Seasonal community succession is observed along the two years sampling period, in addition to a continual annual drift in the microbial populations. The core community of the WWTP bioreactor was also studied, where a small fraction of sequence variants constituted a large fraction of the total abundance. This core microbiome stability along the sampling period and the likewise dissimilarity patterns along the temperature gradient makes this feature a good candidate for a new process control in WWTPs.

Regional community assembly drivers and microbial environmental sources shaping bacterioplankton in an alpine lacustrine district (Pyrenees, Spain).

Microbial communities in natural ecosystems are subject to strong ecological rules. The study of local communities along a regional metacommunity can reveal patterns of community assembly, and disentangle the underlying ecological processes. In particular, we seek drivers of community assembly at the regional scale using a large lacustrine dataset (>300 lakes) along the geographical, limnological and physico-chemical gradients in the Pyrenees. By using high throughput amplicon sequencing of the 16S rRNA gene, and inferring environmental sources of bacterial immigrants, we showed that surface aquatic bacterial assemblages were strongly influenced by terrestrial populations from soil, biofilms or sediments, and primarily selected by a pH-alkalinity gradient. Indeed, source proportions explained 27% of the community variation, and chemistry 15% of the total variation, half of it shared with the sources. Major taxonomic groups such as Verrucomicrobia, Actinobacteria and Bacteroidetes showed higher aquatic affinities than Parcubacteria, Gammaproteobacteria, Alphaproteobacteria or Betaproteobacteria, which may be recruited and selected through different hydrographic habitats. A regional fingerprint was observed with lower alpha diversity and higher beta diversity in the central Pyrenees than in both ends. We suggest an ecological succession process, likely influenced by complex interactions of environmental source dispersal and environmental filtering along the mountain range geography.

A Randomized Trait Community Clustering approach to unveil consistent environmental thresholds in community assembly.

Similarities and differences of phenotypes within local co-occurring species hold the key to inferring the contribution of stochastic or deterministic processes in community assembly. Developing both phylogenetic-based and trait-based quantitative methods to unravel these processes is a major aim in community ecology. We developed a trait-based approach that: (i) assesses if a community trait clustering pattern is related to increasing environmental constraints along a gradient; and (ii) determines quantitative thresholds for an environmental variable along a gradient to interpret changes in prevailing community assembly drivers. We used a regional set of natural shallow saline ponds covering a wide salinity gradient (0.1-40% w/v). We identify a consistent discrete salinity threshold (ca. 5%) for microbial community assembly drivers. Above 5% salinity a strong environmental filtering prevailed as an assembly force, whereas a combination of biotic and abiotic factors dominated at lower salinities. This method provides a conceptual approach to identify consistent environmental thresholds in community assembly and enables quantitative predictions for the ecological impact of environmental changes.

High planktonic diversity in mountain lakes contains similar contributions of autotrophic, heterotrophic and parasitic eukaryotic life forms.

A rich eukaryotic planktonic community exists in high-mountain lakes despite the diluted, oligotrophic and cold, harsh prevailing conditions. Attempts of an overarching appraisal have been traditionally hampered by observational limitations of small, colorless, and soft eukaryotes. We aimed to uncover the regional eukaryotic biodiversity of a mountain lakes district to obtain general conclusions on diversity patterns, dominance, geographic diversification, and food-web players common to oligotrophic worldwide distributed freshwater systems. An unprecedented survey of 227 high-altitude lakes comprising large environmental gradients was carried out using Illumina massive tag sequencing of the 18S rRNA gene. We observed a large Chrysophyceae dominance in richness, abundance and novelty, and unveiled an unexpected richness in heterotrophic phagotrophs and parasites. In particular, Cercozoa and Chytridiomycota showed diversity features similar to the dominant autotrophic groups. The prominent beta-dispersion shown by parasites suggests highly specific interactions and a relevant role in food webs. Interestingly, the freshwater Pyrenean metacommunity contained more diverse specific populations than its closest marine oligotrophic equivalent, with consistently higher beta-diversity. The relevance of unseen groups opens new perspectives for the better understanding of planktonic food webs. Mountain lakes, with remarkable environmental idiosyncrasies, may be suitable environments for the genetic diversification of microscopic eukaryotic life forms.

Consistent changes in the taxonomic structure and functional attributes of bacterial communities during primary succession.

Ecologists have long studied primary succession, the changes that occur in biological communities after initial colonization of an environment. Most of this work has focused on succession in plant communities, laying the conceptual foundation for much of what we currently know about community assembly patterns over time. Because of their prevalence and importance in ecosystems, an increasing number of studies have focused on microbial community dynamics during succession. Here, we conducted a meta-analysis of bacterial primary succession patterns across a range of distinct habitats, including the infant gut, plant surfaces, soil chronosequences, and aquatic environments, to determine whether consistent changes in bacterial diversity, community composition, and functional traits are evident over the course of succession. Although these distinct habitats harbor unique bacterial communities, we were able to identify patterns in community assembly that were shared across habitat types. We found an increase in taxonomic and functional diversity with time while the taxonomic composition and functional profiles of communities became less variable (lower beta diversity) in late successional stages. In addition, we found consistent decreases in the rRNA operon copy number and in the high-efficient phosphate assimilation process (Pst system) suggesting that reductions in resource availability during succession select for taxa adapted to low-resource conditions. Together, these results highlight that, like many plant communities, microbial communities also exhibit predictable patterns during primary succession.

High occurrence of Pacearchaeota and Woesearchaeota (Archaea superphylum DPANN) in the surface waters of oligotrophic high-altitude lakes.

We carried out a regional survey on the archaea composition from surface waters of > 300 high-altitude Pyrenean lakes (average altitude 2300 m, pH range 4.4-10.1) by 16S rRNA gene tag sequencing. Relative Archaea abundances ranged between 0% and 6.3% of total prokaryotes amplicons in the polymerase chain reaction (PCR) mixture, and we detected 769 operational taxonomic units (OTUs; grouped at 97% identity) that split into 13 different lineages, with altitude and pH having a significant effect on the community composition. Woesearchaeota and Pacearchaeota (formerly Euryarchaeota DHVEG-6 cluster) dominated the data set (83% of total OTUS), showed a high occurrence (presence in c. 75% of the lakes) and had relative abundances significantly and positively correlated with the phylogenetic diversity of bacterial communities. Micrarchaeota-Diapherotrites (formerly Euryarchaeota MEG cluster), Methanomicrobia, Thermoplasmata and ammonia-oxidizing thaumarchaeota (AOA) showed relative abundances between 1% and 3% and occurrences between 14% and 26%. Minor lineages were SM1K20, Aenigmarchaeota (formerly Euryarchaeota DSEG cluster), Methanobacteria, Bathyarchaeota and SCG. Environmental preferences substantially differed among lineages, with Aenigmarchaeota and Methanomicrobia having the largest habitat breadth, and Thermoplasmata, AOA and Micrarchaeota having the smallest. Pacearchaeota and Woesearchaeota had been mostly reported from saline habitats and sediments, but surface waters of oligotrophic alpine lakes are suitable environments for such ecologically spread and genetically diverse archaeal lineages.

Ecological Specialization of Two Photobiont-Specific Maritime Cyanolichen Species of the Genus Lichina.

All fungi in the class Lichinomycetes are lichen-forming and exclusively associate with cyanobacteria. Two closely related maritime species of the genus Lichina (L. confinis and L. pygmaea) show similar distribution ranges in the Northeast Atlantic, commonly co-occurring at the same rocky shores but occupying different littoral zones. By means of 16S rRNA and phycocyanin operon markers we studied a) the phylogenetic relationships of cyanobionts associated with these species, b) the match of divergence times between both symbionts, and c) whether Lichina species differ in photobiont association and in how geography and ecology affect selectivity. The cyanobionts studied are closely related to both marine and freshwater strains of the genus Rivularia. We found evidence of a high specificity to particular cyanobiont lineages in both species: Lichina pygmaea and L. confinis incorporate specific lineages of Rivularia that do not overlap at the haplotype nor the OTU levels. Dating divergences of the fungal and cyanobacterial partners revealed an asynchronous origin of both lineages. Within each fungal species, selectivity varied across the studied area, influenced by environmental conditions (both atmospheric and marine), although patterns were highly correlated between both lichen taxa. Ecological speciation due to the differential association of photobionts to each littoral zone is suspected to have occurred in marine Lichina.

Lichen myco- and photobiont diversity and their relationships at the edge of life (McMurdo Dry Valleys, Antarctica).

Lichen-forming fungi are among the most diverse group of organisms in Antarctica. Being poikilohydric, lichens are able to cope with harsh environmental conditions that exclude other organisms like vascular plants. The McMurdo Dry Valleys (Victoria Land, Continental Antarctica) are a hyperarid cold desert where macroscopic life is reduced to a few lichens and bryophyte species. We investigated the diversity of lichen-forming fungi and their associated photobionts in three valleys (Garwood, Marshall, and Miers). Correct identification of lichen-forming fungi from extreme ecosystems is complicated by the presence of numerous sterile and extremely modified thalli. To overcome this problem, we used a combined approach for the identification of the species present in the area, the first involving identification by means of standard characters and the second, two DNA-based (ITS region) species delimitation methods (General Mixed Yule-Coalescent model and genetic distances). In addition, we also used ITS sequences for the identification of the photobionts associated with the mycobionts. We studied the relationships between both bionts and assessed the degree of selectivity and specificity found in those associations. We also looked for landscape level spatial patterns in these associations. The two DNA-based methods performed quite differently, but 27 species of lichen-forming fungi and five putative species of photobionts were found in the studied area. Although there was a general trend for low selectivity in the relationships, high specificity was found in some associations and differential selectivity was observed in some lichen-forming fungi. No spatial structure was detected in the distribution of photobionts in the studied area.