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.].

Comparative transcriptomics of fungal endophytes in co-culture with their moss host Dicranum scoparium reveals fungal trophic lability and moss unchanged to slightly increased growth rates.

Mosses harbor fungi whose interactions within their hosts remain largely unexplored. Trophic ranges of fungal endophytes from the moss Dicranum scoparium were hypothesized to encompass saprotrophism. This moss is an ideal host to study fungal trophic lability because of its natural senescence gradient, and because it can be grown axenically. Dicranum scoparium was co-cultured with each of eight endophytic fungi isolated from naturally occurring D. scoparium. Moss growth rates, and gene expression levels (RNA sequencing) of fungi and D. scoparium, were compared between axenic and co-culture treatments. Functional lability of two fungal endophytes was tested by comparing their RNA expression levels when colonizing living vs dead gametophytes. Growth rates of D. scoparium were unchanged, or increased, when in co-culture. One fungal isolate (Hyaloscyphaceae sp.) that promoted moss growth was associated with differential expression of auxin-related genes. When grown with living vs dead gametophytes, Coniochaeta sp. switched from having upregulated carbohydrate transporter activity to upregulated oxidation-based degradation, suggesting an endophytism to saprotrophism transition. However, no such transition was detected for Hyaloscyphaceae sp. Individually, fungal endophytes did not negatively impact growth rates of D. scoparium. Our results support the long-standing hypothesis that some fungal endophytes can switch to saprotrophism.

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.

Community dynamics of soil-borne fungal communities along elevation gradients in neotropical and palaeotropical forests.

Because of their steep gradients in abiotic and biotic factors, mountains offer an ideal setting to illuminate the mechanisms that underlie patterns of species distributions and community assembly. We compared the composition of taxonomically and functionally diverse fungal communities in soils along five elevational gradients in mountains of the Neo- and Palaeotropics (northern Argentina, southern Brazil, Panama, Malaysian Borneo and Papua New Guinea). Both the richness and composition of soil fungal communities reflect environmental factors, particularly temperature and soil pH, with some shared patterns among neotropical and palaeotropical regions. Community dynamics are characterized by replacement of species along elevation gradients, implying a relatively narrow elevation range for most fungi, which appears to be driven by contrasting environmental preferences among both functional and taxonomic groups. For functional groups dependent on symbioses with plants (especially ectomycorrhizal fungi), the distribution of host plants drives richness and community composition, resulting in important differences in elevational patterns between neotropical and palaeotropical montane communities. The pronounced compositional and functional turnover along elevation gradients implies that tropical montane forest fungi will be sensitive to climate change, resulting in shifts in composition and functionality over time.

Molybdenum threshold for ecosystem scale alternative vanadium nitrogenase activity in boreal forests.

Biological nitrogen fixation (BNF) by microorganisms associated with cryptogamic covers, such as cyanolichens and bryophytes, is a primary source of fixed nitrogen in pristine, high-latitude ecosystems. On land, low molybdenum (Mo) availability has been shown to limit BNF by the most common form of nitrogenase (Nase), which requires Mo in its active site. Vanadium (V) and iron-only Nases have been suggested as viable alternatives to countering Mo limitation of BNF; however, field data supporting this long-standing hypothesis have been lacking. Here, we elucidate the contribution of vanadium nitrogenase (V-Nase) to BNF by cyanolichens across a 600-km latitudinal transect in eastern boreal forests of North America. Widespread V-Nase activity was detected (∼15-50% of total BNF rates), with most of the activity found in the northern part of the transect. We observed a 3-fold increase of V-Nase contribution during the 20-wk growing season. By including the contribution of V-Nase to BNF, estimates of new N input by cyanolichens increase by up to 30%. We find that variability in V-based BNF is strongly related to Mo availability, and we identify a Mo threshold of ∼250 ng·glichen-1 for the onset of V-based BNF. Our results provide compelling ecosystem-scale evidence for the use of the V-Nase as a surrogate enzyme that contributes to BNF when Mo is limiting. Given widespread findings of terrestrial Mo limitation, including the carbon-rich circumboreal belt where global change is most rapid, additional consideration of V-based BNF is required in experimental and modeling studies of terrestrial biogeochemistry.

Exploring syntenic conservation across genomes for phylogenetic studies of organisms subjected to horizontal gene transfers: A case study with Cyanobacteria and cyanolichens.

Understanding the evolutionary history of symbiotic Cyanobacteria at a fine scale is essential to unveil patterns of associations with their hosts and factors driving their spatiotemporal interactions. As for bacteria in general, Horizontal Gene Transfers (HGT) are expected to be rampant throughout their evolution, which justified the use of single-locus phylogenies in macroevolutionary studies of these photoautotrophic bacteria. Genomic approaches have greatly increased the amount of molecular data available, but the selection of orthologous, congruent genes that are more likely to reflect bacterial macroevolutionary histories remains problematic. In this study, we developed a synteny-based approach and searched for Collinear Orthologous Regions (COR), under the assumption that genes that are present in the same order and orientation across a wide monophyletic clade are less likely to have undergone HGT. We searched sixteen reference Nostocales genomes and identified 99 genes, part of 28 COR comprising three to eight genes each. We then developed a bioinformatic pipeline, designed to minimize inter-genome contamination and processed twelve Nostoc-associated lichen metagenomes. This reduced our original dataset to 90 genes representing 25 COR, which were used to infer phylogenetic relationships within Nostocales and among lichenized Cyanobacteria. This dataset was narrowed down further to 71 genes representing 22 COR by selecting only genes part of one (largest) operon per COR. We found a relatively high level of congruence among trees derived from the 90-gene dataset, but congruence was only slightly higher among genes within a COR compared to genes across COR. However, topological congruence was significantly higher among the 71 genes part of one operon per COR. Nostocales phylogenies resulting from concatenation and species tree approaches based on the 90- and 71-gene datasets were highly congruent, but the most highly supported result was obtained when using synteny, collinearity, and operon information (i.e., 71-gene dataset) as gene selection criteria, which outperformed larger datasets with more genes.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata.

BACKGROUND: Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution. RESULTS: A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting. CONCLUSIONS: The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

Differential gene expression associated with fungal trophic shifts along the senescence gradient of the moss Dicranum scoparium.

Bryophytes harbour microbiomes, including diverse communities of fungi. The molecular mechanisms by which perennial mosses interact with these fungal partners along their senescence gradients are unknown, yet this is an ideal system to study variation in gene expression associated with trophic state transitions. We investigated differentially expressed genes of fungal communities and their host Dicranum scoparium across its naturally occurring senescence gradient using a metatranscriptomic approach. Higher activity of fungal nutrient-related (carbon, nitrogen, phosphorus and sulfur) transporters and Carbohydrate-Active enZyme (CAZy) genes was detected toward the bottom, partially decomposed, layer of the moss. The most prominent variation in the expression levels of fungal nutrient transporters was from inorganic nitrogen-related transporters, whereas the breakdown of organonitrogens was detected as the most enriched gene ontology term for the host D. scoparium, for those transcripts having higher expression in the partially decomposed layer. The abundance of bacterial rRNA transcripts suggested that more living members of Cyanobacteria are associated with the photosynthetic layer of D. scoparium, while members of Rhizobiales are detected throughout the gametophytes. Plant genes for specific fungal-plant communication, including defense responses, were differentially expressed, suggesting that different genetic pathways are involved in plant-microbe crosstalk in photosynthetic tissues compared to partially decomposed tissues.

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.

RNA-based analyses reveal fungal communities structured by a senescence gradient in the moss Dicranum scoparium and the presence of putative multi-trophic fungi.

Diverse plant-associated fungi are thought to have symbiotrophic and saprotrophic states because they can be isolated from both dead and living plant tissues. However, such tissues often are separated in time and space, and fungal activity at various stages of plant senescence is rarely assessed directly in fungal community studies. We used fungal ribosomal RNA metatranscriptomics to detect active fungal communities across a natural senescence gradient within wild-collected gametophytes of Dicranum scoparium (Bryophyta) to understand the distribution of active fungal communities in adjacent living, senescing and dead tissues. Ascomycota were active in all tissues across the senescence gradient. By contrast, Basidiomycota were prevalent and active in senescing and dead tissues. Several fungi were detected as active in living and dead tissues, suggesting their capacity for multi-trophy. Differences in community assembly detected by metatranscriptomics were echoed by amplicon sequencing of cDNA and compared to culture-based inferences and observation of fungal fruit bodies in the field. The combination of amplicon sequencing of cDNA and metatranscriptomics is promising for studying symbiotic systems with complex microbial diversity, allowing for the simultaneous detection of their presence and activity.

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.

Stable isotope analyses reveal previously unknown trophic mode diversity in the Hymenochaetales.

PREMISE OF THE STUDY: The Hymenochaetales are dominated by lignicolous saprotrophic fungi involved in wood decay. However, the group also includes bryophilous and terricolous taxa, but their modes of nutrition are not clear. Here, we investigate patterns of carbon and nitrogen utilization in numerous non-lignicolous Hymenochaetales and provide a phylogenetic context in which these non-canonical ecological guilds arose. METHODS: We combined stable isotope analyses of δ13 C and δ15 N and phylogenetic analyses to explore assignment and evolution of nutritional modes. Clustering procedures and statistical tests were performed to assign trophic modes to Hymenochaetales and test for differences between varying ecologies. Genomes of Hymenochaetales were mined for presence of enzymes involved in plant cell wall and lignin degradation and sucrolytic activity. KEY RESULTS: Three different trophic clusters were detected - biotrophic, saprotrophic, and a second biotrophic cluster including many bryophilous Hymenochaetales and mosses. Non-lignicolous Hymenochaetales are generally biotrophic. All lignicolous Hymenochaetales clustered as saprotrophic and most terricolous Hymenochaetales clustered as ectomycorrhizal. Overall, at least 15 species of Hymenochaetales are inferred as biotrophic. Bryophilous species of Rickenella can degrade plant cell walls and lignin, and cleave sucrose to glucose consistent with a parasitic or endophytic life style. CONCLUSIONS: Most non-lignicolous Hymenochaetales are biotrophic. Stable isotope values of many bryophilous Hymenochaetales cluster as ectomycorrhizal or in a biotrophic cluster indicative of parasitism or an endophytic life style. Overall, trophic mode diversity in the Hymenochaetales is greater than anticipated, and non-lignicolous ecological traits and biotrophic modes of nutrition are evolutionarily derived features.

The adaptive radiation of lichen-forming Teloschistaceae is associated with sunscreening pigments and a bark-to-rock substrate shift.

Adaptive radiations play key roles in the generation of biodiversity and biological novelty, and therefore understanding the factors that drive them remains one of the most important challenges of evolutionary biology. Although both intrinsic innovations and extrinsic ecological opportunities contribute to diversification bursts, few studies have looked at the synergistic effect of such factors. Here we investigate the Teloschistales (Ascomycota), a group of >1,000 lichenized species with variation in species richness and phenotypic traits that hinted at a potential adaptive radiation. We found evidence for a dramatic increase in diversification rate for one of four families within this order--Teloschistaceae--which occurred ∼ 100 Mya (Late Cretaceous) and was associated with a switch from bark to rock and from shady to sun-exposed habitats. This adaptation to sunny habitats is likely to have been enabled by a contemporaneous key novel phenotypic innovation: the production in both vegetative structure (thallus) and fruiting body (apothecia) of anthraquinones, secondary metabolites known to protect against UV light. We found that the two ecological factors (sun exposure and rock substrate) and the phenotypic innovation (anthraquinones in the thallus) were all significant when testing for state-dependent shifts in diversification rates, and together they seem likely to be responsible for the success of the Teloschistaceae, one of the largest lichen-forming fungal lineages. Our results support the idea that adaptive radiations are driven not by a single factor or key innovation, but require a serendipitous combination of both intrinsic biotic and extrinsic abiotic and ecological factors.

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.

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.

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.

Phylogenetic analyses of eurotiomycetous endophytes reveal their close affinities to Chaetothyriales, Eurotiales, and a new order - Phaeomoniellales.

Symbiotic fungi living in plants as endophytes, and in lichens as endolichenic fungi, cause no apparent symptoms to their hosts. They are ubiquitous, ecologically important, hyperdiverse, and represent a rich source of secondary compounds for new pharmaceutical and biocontrol products. Due in part to the lack of visible reproductive structures and other distinctive phenotypic traits for many species, the diversity and phylogenetic affiliations of these cryptic fungi are often poorly known. The goal of this study was to determine the phylogenetic placement of representative endophytes within the Eurotiomycetes (Pezizomycotina, Ascomycota), one of the most diverse and evolutionarily dynamic fungal classes, and to use that information to infer processes of macroevolution in trophic modes. Sequences of a single locus marker spanning the nuclear ribosomal internal transcribed spacer region (nrITS) and 600 base pairs at the 5' end of the nuclear ribosomal large subunit (nrLSU) were obtained from previous studies of >6000 endophytic and endolichenic fungi from diverse biogeographic locations and hosts. We conducted phylum-wide phylogenetic searches using this marker to determine which fungal strains belonged to Eurotiomycetes and the results were used as the basis for a class-wide, seven-locus phylogenetic study focusing on endophytic and endolichenic Eurotiomycetes. Our cumulative supermatrix-based analyses revealed that representative endophytes within Eurotiomycetes are distributed in three main clades: Eurotiales, Chaetothyriales and Phaeomoniellales ord. nov., a clade that had not yet been described formally. This new order, described herein, is sister to the clade including Verrucariales and Chaetothyriales. It appears to consist mainly of endophytes and plant pathogens. Morphological characters of endophytic Phaeomoniellales resemble those of the pathogenic genus Phaeomoniella. This study highlights the capacity of endophytic and endolichenic fungi to expand our understanding of the ecological modes associated with particular clades, and provides a first estimation of their phylogenetic relationships in the Eurotiomycetes.

Is vanadium a biometal for boreal cyanolichens?

Molybdenum (Mo) nitrogenase has long been considered the predominant isoenzyme responsible for dinitrogen fixation worldwide. Recent findings have challenged the paradigm of Mo hegemony, and highlighted the role of alternative nitrogenases, such as the vanadium-nitrogenase. Here, we first characterized homeostasis of vanadium (V) along with other metals in situ in the dinitrogen fixing cyanolichen Peltigera aphthosa. These lichens were sampled in natural sites exposed to various levels of atmospheric metal deposition. These results were compared with laboratory experiments where Anabaena variabilis, which is also hosting the V-nitrogenase, and a relatively close relative of the lichen cyanobiont Nostoc, was subjected to various levels of V. We report here that V is preferentially allocated to cephalodia, specialized structures where dinitrogen fixation occurs in tri-membered lichens. This specific allocation is biologically controlled and tightly regulated. Vanadium homeostasis in lichen cephalodia exposed to various V concentrations is comparable to the one observed in Anabaena variabilis and other dinitrogen fixing organisms using V-nitrogenase. Overall, our findings support current hypotheses that V could be a more important factor in mediating nitrogen input in high latitude ecosystems than previously recognized. They invite the reassessment of current theoretical models linking metal dynamics and dinitrogen fixation in boreal and subarctic ecosystems.

A multigene phylogenetic synthesis for the class Lecanoromycetes (Ascomycota): 1307 fungi representing 1139 infrageneric taxa, 317 genera and 66 families.

The Lecanoromycetes is the largest class of lichenized Fungi, and one of the most species-rich classes in the kingdom. Here we provide a multigene phylogenetic synthesis (using three ribosomal RNA-coding and two protein-coding genes) of the Lecanoromycetes based on 642 newly generated and 3329 publicly available sequences representing 1139 taxa, 317 genera, 66 families, 17 orders and five subclasses (four currently recognized: Acarosporomycetidae, Lecanoromycetidae, Ostropomycetidae, Umbilicariomycetidae; and one provisionarily recognized, 'Candelariomycetidae'). Maximum likelihood phylogenetic analyses on four multigene datasets assembled using a cumulative supermatrix approach with a progressively higher number of species and missing data (5-gene, 5+4-gene, 5+4+3-gene and 5+4+3+2-gene datasets) show that the current classification includes non-monophyletic taxa at various ranks, which need to be recircumscribed and require revisionary treatments based on denser taxon sampling and more loci. Two newly circumscribed orders (Arctomiales and Hymeneliales in the Ostropomycetidae) and three families (Ramboldiaceae and Psilolechiaceae in the Lecanorales, and Strangosporaceae in the Lecanoromycetes inc. sed.) are introduced. The potential resurrection of the families Eigleraceae and Lopadiaceae is considered here to alleviate phylogenetic and classification disparities. An overview of the photobionts associated with the main fungal lineages in the Lecanoromycetes based on available published records is provided. A revised schematic classification at the family level in the phylogenetic context of widely accepted and newly revealed relationships across Lecanoromycetes is included. The cumulative addition of taxa with an increasing amount of missing data (i.e., a cumulative supermatrix approach, starting with taxa for which sequences were available for all five targeted genes and ending with the addition of taxa for which only two genes have been sequenced) revealed relatively stable relationships for many families and orders. However, the increasing number of taxa without the addition of more loci also resulted in an expected substantial loss of phylogenetic resolving power and support (especially for deep phylogenetic relationships), potentially including the misplacements of several taxa. Future phylogenetic analyses should include additional single copy protein-coding markers in order to improve the tree of the Lecanoromycetes. As part of this study, a new module ("Hypha") of the freely available Mesquite software was developed to compare and display the internodal support values derived from this cumulative supermatrix approach.