Environmental drivers and cryptic biodiversity hotspots define endophytes in Earth's largest terrestrial biome.

Understanding how symbiotic associations differ across environmental gradients is key to predicting the fate of symbioses as environments change, and it is vital for detecting global reservoirs of symbiont biodiversity in a changing world.1,2,3 However, sampling of symbiotic partners at the full-biome scale is difficult and rare. As Earth's largest terrestrial biome, boreal forests influence carbon dynamics and climate regulation at a planetary scale. Plants and lichens in this biome host the highest known phylogenetic diversity of fungal endophytes, which occur within healthy photosynthetic tissues and can influence hosts' resilience to stress.4,5 We examined how communities of endophytes are structured across the climate gradient of the boreal biome, focusing on the dominant plant and lichen species occurring across the entire south-to-north span of the boreal zone in eastern North America. Although often invoked for understanding the distribution of biodiversity, neither a latitudinal gradient nor mid-domain effect5,6,7 can explain variation in endophyte diversity at this trans-biome scale. Instead, analyses considering shifts in forest characteristics, Picea biomass and age, and nutrients in host tissues from 46° to 58° N reveal strong and distinctive signatures of climate in defining endophyte assemblages in each host lineage. Host breadth of endophytes varies with climate factors, and biodiversity hotspots can be identified at plant-community transitions across the boreal zone at a global scale. Placed against a backdrop of global circumboreal sampling,4 our study reveals the sensitivity of endophytic fungi, their reservoirs of biodiversity, and their important symbiotic associations, to climate.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria.

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Phylogenetic structure of specialization: A new approach that integrates partner availability and phylogenetic diversity to quantify biotic specialization in ecological networks.

Biotic specialization holds information about the assembly, evolution, and stability of biological communities. Partner availabilities can play an important role in enabling species interactions, where uneven partner availabilities can bias estimates of biotic specialization when using phylogenetic diversity indices. It is therefore important to account for partner availability when characterizing biotic specialization using phylogenies. We developed an index, phylogenetic structure of specialization (PSS), that avoids bias from uneven partner availabilities by uncoupling the null models for interaction frequency and phylogenetic distance. We incorporate the deviation between observed and random interaction frequencies as weights into the calculation of partner phylogenetic α-diversity. To calculate the PSS index, we then compare observed partner phylogenetic α-diversity to a null distribution generated by randomizing phylogenetic distances among the same number of partners. PSS quantifies the phylogenetic structure (i.e., clustered, overdispersed, or random) of the partners of a focal species. We show with simulations that the PSS index is not correlated with network properties, which allows comparisons across multiple systems. We also implemented PSS on empirical networks of host-parasite, avian seed-dispersal, lichenized fungi-cyanobacteria, and hummingbird pollination interactions. Across these systems, a large proportion of taxa interact with phylogenetically random partners according to PSS, sometimes to a larger extent than detected with an existing method that does not account for partner availability. We also found that many taxa interact with phylogenetically clustered partners, while taxa with overdispersed partners were rare. We argue that species with phylogenetically overdispersed partners have often been misinterpreted as generalists when they should be considered specialists. Our results highlight the important role of randomness in shaping interaction networks, even in highly intimate symbioses, and provide a much-needed quantitative framework to assess the role that evolutionary history and symbiotic specialization play in shaping patterns of biodiversity. PSS is available as an R package at https://github.com/cjpardodelahoz/pss.

Ecological generalism drives hyperdiversity of secondary metabolite gene clusters in xylarialean endophytes.

Although secondary metabolites are typically associated with competitive or pathogenic interactions, the high bioactivity of endophytic fungi in the Xylariales, coupled with their abundance and broad host ranges spanning all lineages of land plants and lichens, suggests that enhanced secondary metabolism might facilitate symbioses with phylogenetically diverse hosts. Here, we examined secondary metabolite gene clusters (SMGCs) across 96 Xylariales genomes in two clades (Xylariaceae s.l. and Hypoxylaceae), including 88 newly sequenced genomes of endophytes and closely related saprotrophs and pathogens. We paired genomic data with extensive metadata on endophyte hosts and substrates, enabling us to examine genomic factors related to the breadth of symbiotic interactions and ecological roles. All genomes contain hyperabundant SMGCs; however, Xylariaceae have increased numbers of gene duplications, horizontal gene transfers (HGTs) and SMGCs. Enhanced metabolic diversity of endophytes is associated with a greater diversity of hosts and increased capacity for lignocellulose decomposition. Our results suggest that, as host and substrate generalists, Xylariaceae endophytes experience greater selection to diversify SMGCs compared with more ecologically specialised Hypoxylaceae species. Overall, our results provide new evidence that SMGCs may facilitate symbiosis with phylogenetically diverse hosts, highlighting the importance of microbial symbioses to drive fungal metabolic diversity.

Climate and seasonality drive the richness and composition of tropical fungal endophytes at a landscape scale.

Understanding how species-rich communities persist is a foundational question in ecology. In tropical forests, tree diversity is structured by edaphic factors, climate, and biotic interactions, with seasonality playing an essential role at landscape scales: wetter and less seasonal forests typically harbor higher tree diversity than more seasonal forests. We posited that the abiotic factors shaping tree diversity extend to hyperdiverse symbionts in leaves-fungal endophytes-that influence plant health, function, and resilience to stress. Through surveys in forests across Panama that considered climate, seasonality, and covarying biotic factors, we demonstrate that endophyte richness varies negatively with temperature seasonality. Endophyte community structure and taxonomic composition reflect both temperature seasonality and climate (mean annual temperature and precipitation). Overall our findings highlight the vital role of climate-related factors in shaping the hyperdiversity of these important and little-known symbionts of the trees that, in turn, form the foundations of tropical forest biodiversity.

Cyanolichen microbiome contains novel viruses that encode genes to promote microbial metabolism.

Lichen thalli are formed through the symbiotic association of a filamentous fungus and photosynthetic green alga and/or cyanobacterium. Recent studies have revealed lichens also host highly diverse communities of secondary fungal and bacterial symbionts, yet few studies have examined the viral component within these complex symbioses. Here, we describe viral biodiversity and functions in cyanolichens collected from across North America and Europe. As current machine-learning viral-detection tools are not trained on complex eukaryotic metagenomes, we first developed efficient methods to remove eukaryotic reads prior to viral detection and a custom pipeline to validate viral contigs predicted with three machine-learning methods. Our resulting high-quality viral data illustrate that every cyanolichen thallus contains diverse viruses that are distinct from viruses in other terrestrial ecosystems. In addition to cyanobacteria, predicted viral hosts include other lichen-associated bacterial lineages and algae, although a large fraction of viral contigs had no host prediction. Functional annotation of cyanolichen viral sequences predicts numerous viral-encoded auxiliary metabolic genes (AMGs) involved in amino acid, nucleotide, and carbohydrate metabolism, including AMGs for secondary metabolism (antibiotics and antimicrobials) and fatty acid biosynthesis. Overall, the diversity of cyanolichen AMGs suggests that viruses may alter microbial interactions within these complex symbiotic assemblages.

Turnover of Lecanoroid Mycobionts and Their Trebouxia Photobionts Along an Elevation Gradient in Bolivia Highlights the Role of Environment in Structuring the Lichen Symbiosis.

Shifts in climate along elevation gradients structure mycobiont-photobiont associations in lichens. We obtained mycobiont (lecanoroid Lecanoraceae) and photobiont (Trebouxia alga) DNA sequences from 89 lichen thalli collected in Bolivia from a ca. 4,700 m elevation gradient encompassing diverse natural communities and environmental conditions. The molecular dataset included six mycobiont loci (ITS, nrLSU, mtSSU, RPB1, RPB2, and MCM7) and two photobiont loci (ITS, rbcL); we designed new primers to amplify Lecanoraceae RPB1 and RPB2 with a nested PCR approach. Mycobionts belonged to Lecanora s.lat., Bryonora, Myriolecis, Protoparmeliopsis, the "Lecanora" polytropa group, and the "L." saligna group. All of these clades except for Lecanora s.lat. occurred only at high elevation. No single species of Lecanoraceae was present along the entire elevation gradient, and individual clades were restricted to a subset of the gradient. Most Lecanoraceae samples represent species which have not previously been sequenced. Trebouxia clade C, which has not previously been recorded in association with species of Lecanoraceae, predominates at low- to mid-elevation sites. Photobionts from Trebouxia clade I occur at the upper extent of mid-elevation forest and at some open, high-elevation sites, while Trebouxia clades A and S dominate open habitats at high elevation. We did not find Trebouxia clade D. Several putative new species were found in Trebouxia clades A, C, and I. These included one putative species in clade A associated with Myriolecis species growing on limestone at high elevation and a novel lineage sister to the rest of clade C associated with Lecanora on bark in low-elevation grassland. Three different kinds of photobiont switching were observed, with certain mycobiont species associating with Trebouxia from different major clades, species within a major clade, or haplotypes within a species. Lecanoraceae mycobionts and Trebouxia photobionts exhibit species turnover along the elevation gradient, but with each partner having a different elevation threshold at which the community shifts completely. A phylogenetically defined sampling of a single diverse family of lichen-forming fungi may be sufficient to document regional patterns of Trebouxia diversity and distribution.

Host availability drives distributions of fungal endophytes in the imperilled boreal realm.

Boreal forests represent the world's largest terrestrial biome and provide ecosystem services of global importance. Highly imperilled by climate change, these forests host Earth's greatest phylogenetic diversity of endophytes, a hyperdiverse group of symbionts that are defined by their occurrence within living, symptomless plant and lichen tissues. Endophytes shape the ecological and evolutionary trajectories of plants and are therefore key to the function and resilience of terrestrial ecosystems. A critical step in linking the ecological functions of endophytes with those of their hosts is to understand the distributions of these symbionts at the global scale; however, turnover in host taxa with geography and climate can confound insights into endophyte biogeography. As a result, global drivers of endophyte diversity and distributions are not known. Here, we leverage sampling from phylogenetically diverse boreal plants and lichens across North America and Eurasia to show that host filtering in distinctive environments, rather than turnover with geographical or environmental distance, is the main determinant of the community composition and diversity of endophytes. We reveal the distinctiveness of boreal endophytes relative to soil fungi worldwide and endophytes from diverse temperate biomes, highlighting a high degree of global endemism. Overall, the distributions of endophytes are directly linked to the availability of compatible hosts, highlighting the role of biotic interactions in shaping fungal communities across large spatial scales, and the threat that climate change poses to biological diversity and function in the imperilled boreal realm.

T-BAS Version 2.1: Tree-Based Alignment Selector Toolkit for Evolutionary Placement of DNA Sequences and Viewing Alignments and Specimen Metadata on Curated and Custom Trees.

The Tree-Based Alignment Selector (T-BAS) toolkit combines phylogenetic-based placement of DNA sequences with alignment and specimen metadata visualization tools in an integrative pipeline for analyzing microbial biodiversity. The release of T-BAS version 2.1 makes available reference phylogenies, supports multilocus sequence placements and permits uploading and downloading trees, alignments, and specimen metadata.

Contrasting Symbiotic Patterns in Two Closely Related Lineages of Trimembered Lichens of the Genus Peltigera.

Species circumscription is key to the characterization of patterns of specificity in symbiotic systems at a macroevolutionary scale. Here, a worldwide phylogenetic framework was used to assess the biodiversity and symbiotic patterns of association among partners in trimembered lichens from the genus Peltigera, section Chloropeltigera. We sequenced six loci of the main fungal partner and performed species discovery and validation analyses to establish putative species boundaries. Single locus phylogenies were used to establish the identity of both photobionts, Nostoc (cyanobacterium) and Coccomyxa (green alga). Distribution and specificity patterns were compared to the closely related clade, section Peltidea, which includes mainly Peltigera species with trimembered thalli. For section Chloropeltigera, eight fungal species (including five newly delimited putative species) were found in association with nine Nostoc phylogroups and two Coccomyxa species. In contrast, eight fungal species (including three newly delimited putative species) in section Peltidea were found in association with only four Nostoc phylogroups and the same two Coccomyxa species as for section Chloropeltigera. This difference in cyanobiont biodiversity between these two sections can potentially be explained by a significantly higher frequency of sexual reproductive structures in species from section Chloropeltigera compared to section Peltidea. Therefore, horizontal transmission of the cyanobiont might be more prevalent in Chloropeltigera species, while vertical transmission might be more common in Peltidea species. All Peltigera species in section Chloropeltigera are generalists in their association with Nostoc compared to more specialized Peltigera species in section Peltidea. Constrained distributions of Peltigera species that associate strictly with one species of green algae (Coccomyxa subellipsoidea) indicate that the availability of the green alga and the specificity of the interaction might be important factors limiting geographic ranges of trimembered Peltigera, in addition to constraints imposed by their interaction with Nostoc partners and by climatic factors.

Contemporaneous radiations of fungi and plants linked to symbiosis.

Interactions between fungi and plants, including parasitism, mutualism, and saprotrophy, have been invoked as key to their respective macroevolutionary success. Here we evaluate the origins of plant-fungal symbioses and saprotrophy using a time-calibrated phylogenetic framework that reveals linked and drastic shifts in diversification rates of each kingdom. Fungal colonization of land was associated with at least two origins of terrestrial green algae and preceded embryophytes (as evidenced by losses of fungal flagellum, ca. 720 Ma), likely facilitating terrestriality through endomycorrhizal and possibly endophytic symbioses. The largest radiation of fungi (Leotiomyceta), the origin of arbuscular mycorrhizae, and the diversification of extant embryophytes occurred ca. 480 Ma. This was followed by the origin of extant lichens. Saprotrophic mushrooms diversified in the Late Paleozoic as forests of seed plants started to dominate the landscape. The subsequent diversification and explosive radiation of Agaricomycetes, and eventually of ectomycorrhizal mushrooms, were associated with the evolution of Pinaceae in the Mesozoic, and establishment of angiosperm-dominated biomes in the Cretaceous.

Bioclimatic factors at an intrabiome scale are more limiting than cyanobiont availability for the lichen-forming genus Peltigera.

PREMISE OF THE STUDY: Factors shaping spatiotemporal patterns of associations in mutualistic systems are poorly understood. We used the lichen-forming fungi Peltigera and their cyanobacterial partners Nostoc to investigate the spatial structure of this symbiosis at an intrabiome scale and to identify potential factors shaping these associations. METHODS: Ninety-three thalli were sampled in Québec, Canada, along a south-north and an east-west transect of ~1300 km each. We identified the two main partners (Peltigera species and Nostoc phylogroups) using molecular markers and modeled the effects of environmental variables and partner occurrence on Peltigera-Nostoc distributions. KEY RESULTS: Peltigera species showed a high degree of specialization toward cyanobionts, whereas two Nostoc phylogroups dominated both transects by associating with several Peltigera species. Peltigera species had narrower ranges than these two main cyanobionts. Distributions of three Peltigera species were highly associated with precipitation and temperature variables, which was not detected for Nostoc phylogroups at this spatial scale. CONCLUSIONS: For these cyanolichens, factors driving patterns of symbiotic associations are scale dependent. Contrary to global-scale findings, generalist Peltigera species were not more widespread within the boreal biome than specialists. Nostoc availability was not the only driver of Peltigera species' geographic ranges; environmental factors also contributed to their intrabiome distributions. Climatic conditions (especially precipitation) limited the range of some Peltigera species more than the range of their cyanobacterial partners at an intrabiome (boreal) scale.

Conserved genomic collinearity as a source of broadly applicable, fast evolving, markers to resolve species complexes: A case study using the lichen-forming genus Peltigera section Polydactylon.

Synteny can be maintained for certain genomic regions across broad phylogenetic groups. In these homologous genomic regions, sites that are under relaxed purifying selection, such as intergenic regions, could be used broadly as markers for population genetic and phylogenetic studies on species complexes. To explore the potential of this approach, we found 125 Collinear Orthologous Regions (COR) ranging from 1 to >10kb across nine genomes representing the Lecanoromycetes and Eurotiomycetes (Pezizomycotina, Ascomycota). Twenty-six of these COR were found in all 24 eurotiomycete genomes surveyed for this study. Given the high abundance and availability of fungal genomes we believe this approach could be adopted for other large groups of fungi outside the Pezizomycotina. Asa proof of concept, we selected three Collinear Orthologous Regions (COR1b, COR3, and COR16), based on synteny analyses of several genomes representing three classes of Ascomycota: Eurotiomycetes, Lecanoromycetes, and Lichinomycetes. COR16, for example, was found across these three classes of fungi. Here we compare the resolving power of these three new markers with five loci commonly used in phylogenetic studies of fungi, using section Polydactylon of the cyanolichen-forming genus Peltigera (Lecanoromycetes) - a clade with several challenging species complexes. Sequence data were subjected to three species discovery and two validating methods. COR markers substantially increased phylogenetic resolution and confidence, and highly contributed to species delimitation. The level of phylogenetic signal provided by each of the COR markers was higher than the commonly used fungal barcode ITS. High cryptic diversity was revealed by all methods. As redefined here, most species represent lineages that have relatively narrower, and more homogeneous biogeographical ranges than previously understood. The scabrosoid clade consists of ten species, seven of which are new. For the dolichorhizoid clade, twenty-two new species were discovered for a total of twenty-nine species in this clade.

T-BAS: Tree-Based Alignment Selector toolkit for phylogenetic-based placement, alignment downloads and metadata visualization: an example with the Pezizomycotina tree of life.

Motivation: High-quality phylogenetic placement of sequence data has the potential to greatly accelerate studies of the diversity, systematics, ecology and functional biology of diverse groups. We developed the Tree-Based Alignment Selector (T-BAS) toolkit to allow evolutionary placement and visualization of diverse DNA sequences representing unknown taxa within a robust phylogenetic context, and to permit the downloading of highly curated, single- and multi-locus alignments for specific clades. Results: In its initial form, T-BAS v1.0 uses a core phylogeny of 979 taxa (including 23 outgroup taxa, as well as 61 orders, 175 families and 496 genera) representing all 13 classes of largest subphylum of Fungi-Pezizomycotina (Ascomycota)-based on sequence alignments for six loci (nr5.8S, nrLSU, nrSSU, mtSSU, RPB1, RPB2 ). T-BAS v1.0 has three main uses: (i) Users may download alignments and voucher tables for members of the Pezizomycotina directly from the reference tree, facilitating systematics studies of focal clades. (ii) Users may upload sequence files with reads representing unknown taxa and place these on the phylogeny using either BLAST or phylogeny-based approaches, and then use the displayed tree to select reference taxa to include when downloading alignments. The placement of unknowns can be performed for large numbers of Sanger sequences obtained from fungal cultures and for alignable, short reads of environmental amplicons. (iii) User-customizable metadata can be visualized on the tree. Availability and Implementation: T-BAS Version 1.0 is available online at http://tbas.hpc.ncsu.edu . Registration is required to access the CIPRES Science Gateway and NSF XSEDE's large computational resources. Contact: icarbon@ncsu.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Biological nitrogen fixation by alternative nitrogenases in boreal cyanolichens: importance of molybdenum availability and implications for current biological nitrogen fixation estimates.

Cryptogamic species and their associated cyanobacteria have attracted the attention of biogeochemists because of their critical roles in the nitrogen cycle through symbiotic and asymbiotic biological fixation of nitrogen (BNF). BNF is mediated by the nitrogenase enzyme, which, in its most common form, requires molybdenum at its active site. Molybdenum has been reported as a limiting nutrient for BNF in many ecosystems, including tropical and temperate forests. Recent studies have suggested that alternative nitrogenases, which use vanadium or iron in place of molybdenum at their active site, might play a more prominent role in natural ecosystems than previously recognized. Here, we studied the occurrence of vanadium, the role of molybdenum availability on vanadium acquisition and the contribution of alternative nitrogenases to BNF in the ubiquitous cyanolichen Peltigera aphthosa s.l. We confirmed the use of the alternative vanadium-based nitrogenase in the Nostoc cyanobiont of these lichens and its substantial contribution to BNF in this organism. We also showed that the acquisition of vanadium is strongly regulated by the abundance of molybdenum. These findings show that alternative nitrogenase can no longer be neglected in natural ecosystems, particularly in molybdenum-limited habitats.

Contributions of North American endophytes to the phylogeny, ecology, and taxonomy of Xylariaceae (Sordariomycetes, Ascomycota).

The Xylariaceae (Sordariomycetes) comprise one of the largest and most diverse families of Ascomycota, with at least 85 accepted genera and ca. 1343 accepted species. In addition to their frequent occurrence as saprotrophs, members of the family often are found as endophytes in living tissues of phylogenetically diverse plants and lichens. Many of these endophytes remain sterile in culture, precluding identification based on morphological characters. Previous studies indicate that endophytes are highly diverse and represent many xylariaceous genera; however, phylogenetic analyses at the family level generally have not included endophytes, such that their contributions to understanding phylogenetic relationships of Xylariaceae are not well known. Here we use a multi-locus, cumulative supermatrix approach to integrate 92 putative species of fungi isolated from plants and lichens into a phylogenetic framework for Xylariaceae. Our collection spans 1933 isolates from living and senescent tissues in five biomes across the continental United States, and here is analyzed in the context of previously published sequence data from described species and additional taxon sampling of type specimens from culture collections. We found that the majority of strains obtained in our surveys can be classified in the hypoxyloid and xylaroid subfamilies, although many also were found outside of these lineages (as currently circumscribed). Many endophytes were placed in lineages previously not known for endophytism. Most endophytes appear to represent novel species, but inferences are limited by potential gaps in public databases. By linking our data, publicly available sequence data, and records of ascomata, we identify many geographically widespread, host-generalist clades capable of symbiotic associations with diverse photosynthetic partners. Concomitant with such cosmopolitan host use and distributions, many xylariaceous endophytes appear to inhabit both living and non-living plant tissues, with potentially important roles as saprotrophs. Overall, our study reveals major gaps in the availability of multi-locus datasets and metadata for this iconic family, and provides new hypotheses regarding the ecology and evolution of endophytism and other trophic modes across the family Xylariaceae.

Interaction type influences ecological network structure more than local abiotic conditions: evidence from endophytic and endolichenic fungi at a continental scale.

Understanding the factors that shape community assembly remains one of the most enduring and important questions in modern ecology. Network theory can reveal rules of community assembly within and across study systems and suggest novel hypotheses regarding the formation and stability of communities. However, such studies generally face the challenge of disentangling the relative influence of factors such as interaction type and environmental conditions on shaping communities and associated networks. Endophytic and endolichenic symbioses, characterized by microbial species that occur within healthy plants and lichen thalli, represent some of the most ubiquitous interactions in nature. Fungi that engage in these symbioses are hyperdiverse, often horizontally transmitted, and functionally beneficial in many cases, and they represent the diversification of multiple phylogenetic groups. We evaluated six measures of ecological network structure for >4100 isolates of endophytic and endolichenic fungi collected systematically from five sites across North America. Our comparison of these co-occurring interactions in biomes ranging from tundra to subtropical forest showed that the type of interactions (i.e., endophytic vs. endolichenic) had a much more pronounced influence on network structure than did environmental conditions. In particular, endophytic networks were less nested, less connected, and more modular than endolichenic networks in all sites. The consistency of the network structure within each interaction type, independent of site, is encouraging for current efforts devoted to gathering metadata on ecological network structure at a global scale. We discuss several mechanisms potentially responsible for such patterns and draw attention to knowledge gaps in our understanding of networks for diverse interaction types.

Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi.

We studied the evolutionary history of the Parmeliaceae (Lecanoromycetes, Ascomycota), one of the largest families of lichen-forming fungi with complex and variable morphologies, also including several lichenicolous fungi. We assembled a six-locus data set including nuclear, mitochondrial and low-copy protein-coding genes from 293 operational taxonomic units (OTUs). The lichenicolous lifestyle originated independently three times in lichenized ancestors within Parmeliaceae, and a new generic name is introduced for one of these fungi. In all cases, the independent origins occurred c. 24 million yr ago. Further, we show that the Paleocene, Eocene and Oligocene were key periods when diversification of major lineages within Parmeliaceae occurred, with subsequent radiations occurring primarily during the Oligocene and Miocene. Our phylogenetic hypothesis supports the independent origin of lichenicolous fungi associated with climatic shifts at the Oligocene-Miocene boundary. Moreover, diversification bursts at different times may be crucial factors driving the diversification of Parmeliaceae. Additionally, our study provides novel insight into evolutionary relationships in this large and diverse family of lichen-forming ascomycetes.

Determination of elemental baseline using peltigeralean lichens from Northeastern Canada (Québec): Initial data collection for long term monitoring of the impact of global climate change on boreal and subarctic area in Canada.

Northeastern Canada is mostly free of anthropogenic activities. The extent to which this territory has been impacted by anthropogenic atmospheric depositions remains to be studied. The main goal of our study was to establish background levels for metals in boreal muscicolous/terricolous macrolichens over non-urbanized areas of northeastern Canada (Québec). Concentrations of 18 elements (Na, Mg, Al, P, K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Mo, Cd, and Pb) were determined for three species of the genus Peltigera (Peltigera aphthosa (L.) Willd. s.l., Peltigera neopolydactyla (Gyeln.) Gyeln. s.l., Peltigera scabrosa Th. Fr. s.l.), and Nephroma arcticum (L.) Torss., along a 1080 km south-north transect and along a of 730 km west-east transect. We report that elemental contents in the sampled lichen thalli are very low and similar to background levels found in other studies performed in pristine places (high elevation or remote ecosystems) throughout the world. Overall, our results demonstrate that most of the boreal and subarctic zone of Québec (northeastern Canada) is still pristine. The elemental baseline established in these lichen populations will contribute to monitor metal pollution in boreal and sub-polar ecosystems due to global climate change and future industrial expansion.