Bayes or bootstrap? A simulation study comparing the performance of Bayesian Markov chain Monte Carlo sampling and bootstrapping in assessing phylogenetic confidence

Bayesian Markov chain Monte Carlo sampling has become increasingly popular in phylogenetics as a method for both estimating the maximum likelihood topology and for assessing nodal confidence. Despite the growing use of posterior probabilities, the relationship between the Bayesian measure of confidence and the most commonly used confidence measure in phylogenetics, the nonparametric bootstrap proportion, is poorly understood. We used computer simulation to investigate the behavior of three phylogenetic confidence methods: Bayesian posterior probabilities calculated via Markov chain Monte Carlo sampling (BMCMC-PP), maximum likelihood bootstrap proportion (ML-BP), and maximum parsimony bootstrap proportion (MP-BP). We simulated the evolution of DNA sequence on 17-taxon topologies under 18 evolutionary scenarios and examined the performance of these methods in assigning confidence to correct monophyletic and incorrect monophyletic groups, and we examined the effects of increasing character number on support value. BMCMC-PP and ML-BP were often strongly correlated with one another but could provide substantially different estimates of support on short internodes. In contrast, BMCMC-PP correlated poorly with MP-BP across most of the simulation conditions that we examined. For a given threshold value, more correct monophyletic groups were supported by BMCMC-PP than by either ML-BP or MP-BP. When threshold values were chosen that fixed the rate of accepting incorrect monophyletic relationship as true at 5%, all three methods recovered most of the correct relationships on the simulated topologies, although BMCMC-PP and ML-BP performed better than MP-BP. BMCMC-PP was usually a less biased predictor of phylogenetic accuracy than either bootstrapping method. BMCMC-PP provided high support values for correct topological bipartitions with fewer characters than was needed for nonparametric bootstrap.     

Widespread occurrence of spliceosomal introns in the rDNA genes of ascomycetes

Spliceosomal (pre-mRNA) introns have previously been found in eukaryotic protein-coding genes, in the small nuclear RNAs of some fungi, and in the small- and large-subunit ribosomal DNA genes of a limited number of ascomycetes. How the majority of these introns originate remains an open question because few proven cases of recent and pervasive intron origin have been documented. We report here the widespread occurrence of spliceosomal introns (69 introns at 27 different sites) in the small- and large-subunit nuclear-encoded rDNA of lichen-forming and free-living members of the Ascomycota. Our analyses suggest that these spliceosomal introns are of relatively recent origin, i.e., within the Euascomycetes, and have arisen through aberrant reverse-splicing (in trans) of free pre-mRNA introns into rRNAs. The spliceosome itself, and not an external agent (e.g., transposable elements, group II introns), may have given rise to these introns. A nonrandom sequence pattern was found at sites flanking the rRNA spliceosomal introns. This pattern (AG-intron-G) closely resembles the proto-splice site (MAG-intron-R) postulated for intron insertions in pre-mRNA genes. The clustered positions of spliceosomal introns on secondary structures suggest that particular rRNA regions are preferred sites for insertion through reverse-splicing.     

Diversity and phylogenetic affinities of foliar fungal endophytes in loblolly pine inferred by culturing and environmental PCR

We examined endophytic fungi in asymptomatic foliage of loblolly pine (Pinus taeda) in North Carolina, U.S.A., with four goals: (i) to evaluate morphotaxa, BLAST matches and groups based on sequence similarity as functional taxonomic units; (ii) to explore methods to maximize phylogenetic signal for environmental datasets, which typically contain many taxa but few characters; (iii) to compare culturing vs. culture-free methods (environmental PCR of surface sterilized foliage) for estimating endophyte diversity and species composition; and (iv) to investigate the relationships between traditional ecological indices (e.g. Shannon index) and phylogenetic diversity (PD) in estimating endophyte diversity and spatial heterogeneity. Endophytes were recovered in culture from 87 of 90 P. taeda leaves sampled, yielding 439 isolates that represented 24 morphotaxa. Sequence data from the nuclear ribosomal internal transcribed spacer (ITS) for 150 isolates revealed 59 distinct ITS genotypes that represented 24 and 37 unique groups based on 90% and 95% sequence similarity, respectively. By recoding ambiguously aligned regions to extract phylogenetic signal and implementing a conservative phylogenetic backbone constraint, we recovered well supported phylogenies based on ca. 600 bp of the nuclear ribosomal large subunit (LSUrDNA) for 72 Ascomycota and Basidiomycota, 145 cultured endophytes and 33 environmental PCR samples. Comparisons with LSUrDNA-delimited species showed that morphotaxa adequately estimated total species richness but rarely corresponded to biologically meaningful groups. ITS BLAST results were variable in their utility, but ITS genotype groups based on 90% sequence similarity were concordant with LSUrDNA-delimited species. Environmental PCR yielded more genotypes per sampling effort and recovered several distinct clades relative to culturing, but some commonly cultured clades were never found (Sordariomycetes) or were rare relative to their high frequency among cultures (Leotiomycetes). In contrast to traditional indices, PD demonstrated spatial heterogeneity in endophyte assemblages among P. taeda trees and study plots. Our results highlight the need for caution in designating taxonomic units based on gross cultural morphology or ITS BLAST matches, the utility of phylogenetic tools for extracting robust phylogenies from environmental samples, the complementarity of culturing and environmental PCR, the utility of PD relative to traditional ecological indices, and the remarkably high diversity of foliar fungal endophytes in this simplified temperate ecosystem.     

Phylogenetic relationships, host affinity, and geographic structure of boreal and arctic endophytes from three major plant lineages

Although associated with all plants, fungal endophytes (microfungi that live within healthy plant tissues) represent an unknown proportion of fungal diversity. While there is a growing appreciation of their ecological importance and human uses, little is known about their host specificity, geographic structure, or phylogenetic relationships. We surveyed endophytic Ascomycota from healthy photosynthetic tissues of three plant species (Huperzia selago, Picea mariana, and Dryas integrifolia, representing lycophytes, conifers, and angiosperms, respectively) in northern and southern boreal forest (Québec, Canada) and arctic tundra (Nunavut, Canada). Endophytes were recovered from all plant species surveyed, and were present in <1-41% of 2 mm2 tissue segments examined per host species. Sequence data from the nuclear ribosomal internal transcribed spacer region (ITS) were obtained for 280 of 558 isolates. Species-accumulation curves based on ITS genotypes remained non-asymptotic, and bootstrap analyses indicated that a large number of genotypes remain to be found. The majority of genotypes were recovered from only a single host species, and only 6% of genotypes were shared between boreal and arctic communities. Two independent Bayesian analyses and a neighbor-joining bootstrapping analysis of combined data from the nuclear large and small ribosomal subunits (LSUrDNA, SSUrDNA; 2.4 kb) showed that boreal and arctic endophytes represent Dothideomycetes, Sordariomycetes, Chaetothyriomycetidae, Leotiomycetes, and Pezizomycetes. Many well-supported phylotypes contained only endophytes despite exhaustive sampling of available sequences of Ascomycota. Together, these data demonstrate greater than expected diversity of endophytes at high-latitude sites and provide a framework for assessing the evolution of these poorly known but ubiquitous symbionts of living plants.     

Low genetic variation and no detectable population structure in aspergillus fumigatus compared to closely related Neosartorya species

Aspergillus fumigatus is an anamorphic euascomycete mold with a ubiquitous presence worldwide. Despite intensive work to understand its success as a pathogen infecting immunosuppressed patients, the population dynamics and recent evolutionary history of A. fumigatus remain understudied. We examined patterns of genetic variation at three intergenic loci for 70 natural isolates from Europe, North America, South America, Asia, Africa, and Australia. The same loci were used to analyze within-population genetic variation for 33 isolates obtained from five geographic locations. Neither data set detected evidence of population differentiation or found any association between the genetic and geographic distances among these isolates. No evidence for genetic differentiation within the two A. fumigatus mating types was detected. The genetic diversity of A. fumigatus, contrasted with that of its close teleomorphic relatives, Neosartorya fischeri and Neosartorya spinosa, is remarkably low.     

Photoautotrophic symbiont and geography are major factors affecting highly structured and diverse bacterial communities in the lichen microbiome

Although common knowledge dictates that the lichen thallus is formed solely by a fungus (mycobiont) that develops a symbiotic relationship with an alga and/or cyanobacterium (photobiont), the non-photoautotrophic bacteria found in lichen microbiomes are increasingly regarded as integral components of lichen thalli. For this study, comparative analyses were conducted on lichen-associated bacterial communities to test for effects of photobiont-types (i.e. green algal vs. cyanobacterial), mycobiont-types and large-scale spatial distances (from tropical to arctic latitudes). Amplicons of the 16S (SSU) rRNA gene were examined using both Sanger sequencing of cloned fragments and barcoded pyrosequencing. Rhizobiales is typically the most abundant and taxonomically diverse order in lichen microbiomes; however, overall bacterial diversity in lichens is shown to be much higher than previously reported. Members of Acidobacteriaceae, Acetobacteraceae, Brucellaceae and sequence group LAR1 are the most commonly found groups across the phylogenetically and geographically broad array of lichens examined here. Major bacterial community trends are significantly correlated with differences in large-scale geography, photobiont-type and mycobiont-type. The lichen as a microcosm represents a structured, unique microbial habitat with greater ecological complexity and bacterial diversity than previously appreciated and can serve as a model system for studying larger ecological and evolutionary principles.     

RECOGNIZING DINOFLAGELLATE SPECIES USING ITS rDNA SEQUENCES1

Dinoflagellate taxonomy is based primarily on morphology and morphometric data that can be difficult to obtain. In contrast, molecular data can be rapidly and cost‐effectively acquired, which has led to a rapid accumulation of sequence data in GenBank. Currently there are no systematic criteria for utilizing taxonomically unassigned sequence data to identify putative species that could in turn serve as a basis for testable hypotheses concerning the taxonomy, diversity, distribution, and toxicity of these organisms. The goal of this research was to evaluate whether simple, uncorrected genetic distances (p) calculated using ITS1/5.8S/ITS2 (ITS region) rDNA sequences could be used to develop criteria for recognizing putative species before formal morphological evaluation and classification. The current analysis used sequences from 81 dinoflagellate species belonging to 14 genera. For this diverse assemblage of dinoflagellate species, the within‐species genetic distances between ITS region copies (p=0.000–0.021 substitutions per site) were consistently less than those observed between species (p=0.042–0.580). Our results indicate that a between‐species uncorrected genetic distance of p≥0.04 could be used to delineate most free‐living dinoflagellate species. Recently evolved species, however, may have ITS p values <0.04 and would require more extensive morphological and genetic analyses to resolve. For most species, the sequence of the dominant ITS region allele has the potential to serve as a unique species‐specific “DNA barcode” that could be used for the rapid identification of dinoflagellates in field and laboratory studies.