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

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.

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.

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.

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.

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.

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.

Phylogenetic placement, species delimitation, and cyanobiont identity of endangered aquatic Peltigera species (lichen-forming Ascomycota, Lecanoromycetes).

• Premise of this study: Aquatic cyanolichens from the genus Peltigera section Hydrothyriae are subject to anthropogenic threats and, therefore, are considered endangered. In this study we addressed the phylogenetic placement of section Hydrothyriae within Peltigera. We delimited species within the section and identified their symbiotic cyanobacteria.• Methods: Species delimitation and population structure were explored using monophyly as a grouping criterion (RAxML) and Structurama based on three protein-coding genes in combination with two nuclear ribosomal loci. The 16S and rbcLX sequences for the cyanobionts were analyzed in the broad phylogenetic context of free-living and symbiotic cyanobacteria.• Key results: We confirm with high confidence the placement of section Hydrothyriae within the monophyletic genus Peltigera; however, its phylogenetic position within the genus remains unsettled. We recovered three distinct monophyletic groups corresponding to three species: P. hydrothyria, P. gowardii s.s., and P. aquatica Miadl. & Lendemer, the latter being formally introduced here. Each species was associated with an exclusive set of Nostoc haplotypes.• Conclusions: The ITS region alone provides sufficient genetic information to distinguish the three morphologically cryptic species within section Hydrothyriae. Section Hydrothyriae seems to be associated with a monophyletic lineage of Nostoc, that has not been found in symbiotic association with other members of Peltigera. Capsosira lowei should be transferred to the genus Nostoc. Potential threats to P. aquatica should be re-examined based on the recognition of two aquatic species in western North America.

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.

Assessing population structure and host specialization in lichenized cyanobacteria.

Coevolutionary theory predicts that the distribution of obligately symbiotic organisms will be determined by the dispersal ability and ecological range of both partners. We examined this prediction for lichen-forming fungi that form obligate symbioses with cyanobacteria. We compared genotypes of both partners of 250 lichens collected at multiple spatial scales in British Columbia, Canada. Multilocus sequence data collected from a subset of 128 of the specimens were used to determine the degree of recombination within the cyanobacterial populations. We found that six distinct clusters of cyanobacterial genotypes are distributed throughout the known global phylogeny of the genus Nostoc, and that each appears to be evolving clonally. Fungal specialization is high, with each species associating with either one or two of the cyanobacterial clusters, while cyanobacterial specialization varies, with clusters associating with between one and 12 different fungal species. Specialization also varies geographically, with some combinations restricted to a single site despite the availability of both partners elsewhere. Photobiont association patterns are determined by a combination of genetically based specificity, spatial population structure, and ecological factors and cannot be easily predicted by photobiont dispersal syndromes.

Molecular phylogeny of the Entomophthoromycota.

The Entomophthoromycota is a ubiquitous group of fungi best known as pathogens of a wide variety of economically important insect pests, and other soil invertebrates. This group of fungi also includes a small number of parasites of reptiles, vertebrates (including humans), macromycetes, fern gametophytes, and desmid algae, as well as some saprobic species. Here we report on recent studies to resolve the phylogenetic relationships within the Entomophthoromycota and to reliably place this group among other basal fungal lineages. Bayesian Interference (BI) and Maximum Likelihood (ML) analyses of three genes (nuclear 18S and 28S rDNA, mitochondrial 16S, and the protein-coding RPB2) as well as non-molecular data consistently and unambiguously identify 31 taxa of Entomophthoromycota as a monophyletic group distinct from other Zygomycota and flagellated fungi. Using the constraints of our multi-gene dataset we constructed the most comprehensive rDNA phylogeny yet available for Entomophthoromycota. The taxa studied here belong to five distinct, well-supported lineages. The Basidiobolus clade is the earliest diverging lineage, comprised of saprobe species of Basidiobolus and the undescribed snake parasite Schizangiella serpentis nom. prov. The Conidiobolus lineage is represented by a paraphyletic grade of trophically diverse species that include saprobes, insect pathogens, and facultative human pathogens. Three well supported and exclusively entomopathogenic lineages in the Entomophthoraceae center around the genera Batkoa, Entomophthora and Zoophthora, although several genera within this crown clade are resolved as non-monophyletic. Ancestral state reconstruction suggests that the ancestor of all Entomophthoromycota was morphologically similar to species of Conidiobolus. Analyses using strict, relaxed, and local molecular clock models documented highly variable DNA substitution rates among lineages of Entomophthoromycota. Despite the complications caused by different rates of molecular evolution among lineages, our dating analysis indicates that the Entomophthoromycota originated 405±90 million years ago. We suggest that entomopathogenic lineages in Entomophthoraceae probably evolved from saprobic or facultatively pathogenic ancestors during or shortly after the evolutionary radiation of the arthropods.

Host and geographic structure of endophytic and endolichenic fungi at a continental scale.

PREMISE OF THE STUDY: Endophytic and endolichenic fungi occur in healthy tissues of plants and lichens, respectively, playing potentially important roles in the ecology and evolution of their hosts. However, previous sampling has not comprehensively evaluated the biotic, biogeographic, and abiotic factors that structure their communities. METHODS: Using molecular data we examined the diversity, composition, and distributions of 4154 endophytic and endolichenic Ascomycota cultured from replicate surveys of ca. 20 plant and lichen species in each of five North American sites (Madrean coniferous forest, Arizona; montane semideciduous forest, North Carolina; scrub forest, Florida; Beringian tundra and forest, western Alaska; subalpine tundra, eastern central Alaska). KEY RESULTS: Endolichenic fungi were more abundant and diverse per host species than endophytes, but communities of endophytes were more diverse overall, reflecting high diversity in mosses and lycophytes. Endophytes of vascular plants were largely distinct from fungal communities that inhabit mosses and lichens. Fungi from closely related hosts from different regions were similar in higher taxonomy, but differed at shallow taxonomic levels. These differences reflected climate factors more strongly than geographic distance alone. CONCLUSIONS: Our study provides a first evaluation of endophytic and endolichenic fungal associations with their hosts at a continental scale. Both plants and lichens harbor abundant and diverse fungal communities whose incidence, diversity, and composition reflect the interplay of climatic patterns, geographic separation, host type, and host lineage. Although culture-free methods will inform future work, our study sets the stage for empirical assessments of ecological specificity, metabolic capability, and comparative genomics.

Reconstructing the early evolution of Fungi using a six-gene phylogeny.

The ancestors of fungi are believed to be simple aquatic forms with flagellated spores, similar to members of the extant phylum Chytridiomycota (chytrids). Current classifications assume that chytrids form an early-diverging clade within the kingdom Fungi and imply a single loss of the spore flagellum, leading to the diversification of terrestrial fungi. Here we develop phylogenetic hypotheses for Fungi using data from six gene regions and nearly 200 species. Our results indicate that there may have been at least four independent losses of the flagellum in the kingdom Fungi. These losses of swimming spores coincided with the evolution of new mechanisms of spore dispersal, such as aerial dispersal in mycelial groups and polar tube eversion in the microsporidia (unicellular forms that lack mitochondria). The enigmatic microsporidia seem to be derived from an endoparasitic chytrid ancestor similar to Rozella allomycis, on the earliest diverging branch of the fungal phylogenetic tree.

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.