Comparative genomics of Mortierella elongata and its bacterial endosymbiont Mycoavidus cysteinexigens.

Endosymbiosis of bacteria by eukaryotes is a defining feature of cellular evolution. In addition to well-known bacterial origins for mitochondria and chloroplasts, multiple origins of bacterial endosymbiosis are known within the cells of diverse animals, plants and fungi. Early-diverging lineages of terrestrial fungi harbor endosymbiotic bacteria belonging to the Burkholderiaceae. We sequenced the metagenome of the soil-inhabiting fungus Mortierella elongata and assembled the complete circular chromosome of its endosymbiont, Mycoavidus cysteinexigens, which we place within a lineage of endofungal symbionts that are sister clade to Burkholderia. The genome of M. elongata strain AG77 features a core set of primary metabolic pathways for degradation of simple carbohydrates and lipid biosynthesis, while the M. cysteinexigens (AG77) genome is reduced in size and function. Experiments using antibiotics to cure the endobacterium from the host demonstrate that the fungal host metabolism is highly modulated by presence/absence of M. cysteinexigens. Independent comparative phylogenomic analyses of fungal and bacterial genomes are consistent with an ancient origin for M. elongata - M. cysteinexigens symbiosis, most likely over 350 million years ago and concomitant with the terrestrialization of Earth and diversification of land fungi and plants.

Metatranscriptomic analysis of ectomycorrhizal roots reveals genes associated with Piloderma-Pinus symbiosis: improved methodologies for assessing gene expression in situ.

Ectomycorrhizal (EM) fungi form symbiotic associations with plant roots that regulate nutrient exchange between forest plants and soil. Environmental metagenomics approaches that employ next-generation sequencing show great promise for studying EM symbioses; however, metatranscriptomic studies have been constrained by the inherent difficulties associated with isolation and sequencing of RNA from mycorrhizae. Here we apply an optimized method for combined DNA/RNA extraction using field-collected EM fungal-pine root clusters, together with protocols for taxonomic identification of expressed ribosomal RNA, and inference of EM function based on plant and fungal metatranscriptomics. We used transcribed portions of ribosomal RNA genes to identify several transcriptionally dominant fungal taxa associated with loblolly pine including Amphinema, Russula and Piloderma spp. One taxon, Piloderma croceum, has a publically available genome that allowed us to identify patterns of gene content and transcript abundance. Over 1500 abundantly expressed Piloderma genes were detected from mycorrhizal roots, including genes for protein metabolism, cell signalling, electron transport, terpene synthesis and other extracellular activities. In contrast, Piloderma gene encoding an ammonia transporter showed highest transcript abundance in soil samples. Our methodology highlights the potential of metatranscriptomics to identify genes associated with symbiosis and ecosystem function using field-collected samples.

High diversity and widespread occurrence of mitotic spore mats in ectomycorrhizal Pezizales.

Fungal mitospores may function as dispersal units and/ or spermatia and thus play a role in distribution and/or mating of species that produce them. Mitospore production in ectomycorrhizal (EcM) Pezizales is rarely reported, but here we document mitospore production by a high diversity of EcM Pezizales on three continents, in both hemispheres. We sequenced the internal transcribed spacer (ITS) and partial large subunit (LSU) nuclear rDNA from 292 spore mats (visible mitospore clumps) collected in Argentina, Chile, China, Mexico and the USA between 2009 and 2012. We collated spore mat ITS sequences with 105 fruit body and 47 EcM root sequences to generate operational taxonomic units (OTUs). Phylogenetic inferences were made through analyses of both molecular data sets. A total of 48 OTUs from spore mats represented six independent EcM Pezizales lineages and included truffles and cup fungi. Three clades of seven OTUs have no known meiospore stage. Mitospores failed to germinate on sterile media, or form ectomycorrhizas on Quercus, Pinus and Populus seedlings, consistent with a hypothesized role of spermatia. The broad geographic range, high frequency and phylogenetic diversity of spore mats produced by EcM Pezizales suggests that a mitospore stage is important for many species in this group in terms of mating, reproduction and/or dispersal.

Phylogenetic lineages in Entomophthoromycota.

Entomophthoromycota is one of six major phylogenetic lineages among the former phylum Zygomycota. These early terrestrial fungi share evolutionarily ancestral characters such as coenocytic mycelium and gametangiogamy as a sexual process resulting in zygospore formation. Previous molecular studies have shown the monophyly of Entomophthoromycota, thus justifying raising the taxonomic status of these fungi to a phylum. Multi-gene phylogenies have identified five major lineages of Entomophthoromycota. In this review we provide a detailed discussion about the biology and taxonomy of these lineages: I) Basidiobolus (Basidiobolomycetes: Basidiobolaceae; primarily saprobic); II) Conidiobolus (Entomophthoromycetes, Ancylistaceae; several clades of saprobes and invertebrate pathogens), as well as three rapidly evolving entomopathogenic lineages in the family Entomophthoraceae centering around; III) Batkoa; IV) Entomophthora and allied genera; and V) the subfamily Erynioideae which includes Zoophthora and allied genera. Molecular phylogenic analysis has recently determined the relationships of several taxa that were previously unresolved based on morphology alone: Eryniopsis, Macrobiotophthora, Massospora, Strongwellsea and two as yet undescribed genera of Basidiobolaceae.

The Populus holobiont: dissecting the effects of plant niches and genotype on the microbiome.

BACKGROUND: Microorganisms serve important functions within numerous eukaryotic host organisms. An understanding of the variation in the plant niche-level microbiome, from rhizosphere soils to plant canopies, is imperative to gain a better understanding of how both the structural and functional processes of microbiomes impact the health of the overall plant holobiome. Using Populus trees as a model ecosystem, we characterized the archaeal/bacterial and fungal microbiome across 30 different tissue-level niches within replicated Populus deltoides and hybrid Populus trichocarpa × deltoides individuals using 16S and ITS2 rRNA gene analyses. RESULTS: Our analyses indicate that archaeal/bacterial and fungal microbiomes varied primarily across broader plant habitat classes (leaves, stems, roots, soils) regardless of plant genotype, except for fungal communities within leaf niches, which were greatly impacted by the host genotype. Differences between tree genotypes are evident in the elevated presence of two potential fungal pathogens, Marssonina brunnea and Septoria sp., on hybrid P. trichocarpa × deltoides trees which may in turn be contributing to divergence in overall microbiome composition. Archaeal/bacterial diversity increased from leaves, to stem, to root, and to soil habitats, whereas fungal diversity was the greatest in stems and soils. CONCLUSIONS: This study provides a holistic understanding of microbiome structure within a bioenergy relevant plant host, one of the most complete niche-level analyses of any plant. As such, it constitutes a detailed atlas or map for further hypothesis testing on the significance of individual microbial taxa within specific niches and habitats of Populus and a baseline for comparisons to other plant species.

A New Perspective on Sustainable Soil Remediation-Case Study Suggests Novel Fungal Genera Could Facilitate in situ Biodegradation of Hazardous Contaminants.

Deciding upon a cost effective and sustainable method to address soil pollution is a challenge for many remedial project managers. High pressure to quickly achieve cleanup goals pushes for energy-intensive remedies that rapidly address the contaminants of concern with established technologies, often leaving little room for research and development especially for slower treatment technologies, such as bioremediation, for the more heavily polluted sites. In the present case study, new genomic approaches have been leveraged to assess fungal biostimulation potential in soils polluted with particularly persistent hydrophobic contaminants. This new approach provides insights into the genetic functions available at a given site in a way never before possible. In particular, this article presents a case study where next generation sequencing (NGS) has been used to categorize fungi in soils from the Atlantic Wood Industries Superfund site in Portsmouth, Virginia. Data suggest that original attempts to harness fungi for bioremediation may have focused on fungal genera poorly suited to survive under heavily polluted site conditions, and that more targeted approaches relying on native indigenous fungi which are better equipped to survive under site specific conditions may be more appropriate.

Polymorphism at the ribosomal DNA spacers and its relation to breeding structure of the widespread mushroom Schizophyllum commune.

The common split-gilled mushroom Schizophyllum commune is found throughout the world on woody substrates. This study addresses the dispersal and population structure of this fungal species by studying the phylogeny and evolutionary dynamics of ribosomal DNA (rDNA) spacer regions. Extensive sampling (n = 195) of sequences of the intergenic spacer region (IGS1) revealed a large number of unique haplotypes (n = 143). The phylogeny of these IGS1 sequences revealed strong geographic patterns and supported three evolutionarily distinct lineages within the global population. The same three geographic lineages were found in phylogenetic analysis of both other rDNA spacer regions (IGS2 and ITS). However, nested clade analysis of the IGS1 phylogeny suggested the population structure of S. commune has undergone recent changes, such as a long distance colonization of western North America from Europe as well as a recent range expansion in the Caribbean. Among all spacer regions, variation in length and nucleotide sequence was observed between but not within the tandem rDNA repeats (arrays). This pattern is consistent with strong within-array and weak among-array homogenizing forces. We present evidence for the suppression of recombination between rDNA arrays on homologous chromosomes that may account for this pattern of concerted evolution.

Widespread occurrence and phylogenetic placement of a soil clone group adds a prominent new branch to the fungal tree of life.

Fungi are one of the most diverse groups of Eukarya and play essential roles in terrestrial ecosystems as decomposers, pathogens and mutualists. This study unifies disparate reports of unclassified fungal sequences from soils of diverse origins and anchors many of them in a well-supported clade of the Ascomycota equivalent to a subphylum. We refer to this clade as Soil Clone Group I (SCGI). We expand the breadth of environments surveyed and develop a taxon-specific primer to amplify 2.4kbp rDNA fragments directly from soil. Our results also expand the known range of this group from North America to Europe and Australia. The ancient origin of SCGI implies that it may represent an important transitional form among the basal Ascomycota groups. SCGI is unusual because it currently represents the only major fungal lineage known only from sequence data. This is an important contribution towards building a more complete fungal phylogeny and highlights the need for further work to determine the function and biology of SCGI taxa.

Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species.

Detailed restriction analyses of many samples often require substantial amounts of time and effort for DNA extraction, restriction digests, Southern blotting, and hybridization. We describe a novel approach that uses the polymerase chain reaction (PCR) for rapid simplified restriction typing and mapping of DNA from many different isolates. DNA fragments up to 2 kilobase pairs in length were efficiently amplified from crude DNA samples of several pathogenic Cryptococcus species, including C. neoformans, C. albidus, C. laurentii, and C. uniguttulatus. Digestion and electrophoresis of the PCR products by using frequent-cutting restriction enzymes produced complex restriction phenotypes (fingerprints) that were often unique for each strain or species. We used the PCR to amplify and analyze restriction pattern variation within three major portions of the ribosomal DNA (rDNA) repeats from these fungi. Detailed mapping of many restriction sites within the rDNA locus was determined by fingerprint analysis of progressively larger PCR fragments sharing a common primer site at one end. As judged by PCR fingerprints, the rDNA of 19 C. neoformans isolates showed no variation for four restriction enzymes that we surveyed. Other Cryptococcus spp. showed varying levels of restriction pattern variation within their rDNAs and were shown to be genetically distinct from C. neoformans. The PCR primers used in this study have also been successfully applied for amplification of rDNAs from other pathogenic and nonpathogenic fungi, including Candida spp., and ought to have wide applicability for clinical detection and other studies.

Ancient and recent patterns of geographic speciation in the oyster mushroom Pleurotus revealed by phylogenetic analysis of ribosomal DNA sequences.

Evidence from molecular systematic studies suggests that many mushroom species may be quite ancient. Gene phylogenies were developed to examine the relationship between reproductive isolation, genetic divergence, and biogeography in oyster mushrooms (Pleurotus). Sequence data were obtained for two regions of DNA from populations belonging to eight intersterility groups (biological species). Phylogenetic analysis of sequences from the 5' portion of the nuclear encoded large subunit rDNA demonstrates an ancient origin for four intersterility groups of broad geographic distribution (world-wide), with a more recent radiation of several intersterility groups that are restricted to the Northern Hemisphere. An expanded analysis using sequence data from the more variable rDNA internal transcribed spacer region also reveals a phylogenetically based pattern of genetic divergence associated with allopatric speciation among populations from different continents in the Northern Hemisphere. The ability of rDNA sequences to resolve phylogenetic relationships among geographically isolated populations within intersterility groups illustrates the importance of biogeography for understanding speciation in Pleurotus. Patterns of geographic distribution among intersterility groups suggest that several species lineages evolved quite early, with recently evolved groups restricted to the Northern Hemisphere and older lineages occurring throughout the world. Based on phylogenetic evidence, analysis of historical biogeography using area cladograms shows that multiple dispersal and vicariance events are responsible for patterns of speciation observed.

Phylogenetic relationships in the mushroom genus Coprinus and dark-spored allies based on sequence data from the nuclear gene coding for the large ribosomal subunit RNA: divergent domains, outgroups, and monophyly.

Phylogenetic relationships were investigated in the mushroom genus Coprinus based on sequence data from the nuclear encoded large-subunit rDNA gene. Forty-seven species of Coprinus and 19 additional species from the families Coprinaceae, Strophariaceae, Bolbitiaceae, Agaricaceae, Podaxaceae, and Montagneaceae were studied. A total of 1360 sites was sequenced across seven divergent domains and intervening sequences. A total of 302 phylogenetically informative characters was found. Ninety-eight percent of the average divergence between taxa was located within the divergent domains, with domains D2 and D8 being most divergent and domains D7 and D10 the least divergent. An empirical test of phylogenetic signal among divergent domains also showed that domains D2 and D3 had the lowest levels of homoplasy. Two equally most parsimonious trees were resolved using Wagner parsimony. A character-state weighted analysis produced 12 equally most parsimonious trees similar to those generated by Wagner parsimony. Phylogenetic analyses employing topological constraints suggest that none of the major taxonomic systems proposed for subgeneric classification is able to completely reflect phylogenetic relationships in Coprinus. A strict consensus integration of the two Wagner trees demonstrates the problematic nature of choosing outgroups within dark-spored mushrooms. The genus Coprinus is found to be polyphyletic and is separated into three distinct clades. Most Coprinus taxa belong to the first two clades, which together form a larger monophyletic group with Lacrymaria and Psathyrella in basal positions. A third clade contains members of Coprinus section Comati as well as the genus Leucocoprinus, Podaxis pistillaris, Montagnea arenaria, and Agaricus pocillator. This third clade is separated from the other species of Coprinus by members of the families Strophariaceae and Bolbitiaceae and the genus Panaeolus.

Abundance and diversity of Schizophyllum commune spore clouds in the Caribbean detected by selective sampling.

Selective spore trapping and molecular genotyping methods were employed to examine potential long-distance gene flow among Caribbean populations of the common mushroom Schizophyllum commune. Spore-trap samples from five locations were analysed using restriction fragment polymorphisms of five enzymatically amplified gene regions. Successful trappings suggested S. commune spores to be abundant in the air, with an estimated sedimentation rate of approximately 18 spores/m2/h. High levels of genetic diversity characterized the spore-trap samples, with as many as 12 alleles observed at a single locus (chitin synthase) over all samples. In addition, spore-trap samples showed significant among sample heterogeneity including geographical population substructure. The ribosomal DNA (rDNA) intergenic spacer displayed the greatest allele frequency differences among samples, clearly separating the samples into those possessing only a South American-type allele and those segregating for both North and South American-type alleles. The molecular variation provided no clear evidence for dispersal over large, aquatic barriers within the Caribbean region, and instead suggested that spore-trapping experiments are primarily reflective of the local, established population.

Extensive genetic divergence associated with speciation in filamentous fungi.

Little is known about genetic differentiation during speciation in fungi. The Collybia dryophila complex (Basidiomycetes: Tricholomataceae) contains several closely related groups of fungi at various levels of evolutionary divergence. Mating compatibility studies show there are several intersterile groups within the complex, three of which are distributed over two or more continents. Genetic relationships among five intersterility groups were compared by the method of DNA.DNA hybridization. Clustering techniques were used to reconstruct phylogenetic relationships of these fungi. Genic identity based on DNA.DNA hybridization corresponds well with mating compatibility. Reduced genic identity was observed between geographically isolated populations of a single mating group. This evidence suggests an allopatric mode of evolution for the C. dryophila group. These results indicate extensive genic divergence is associated with the appearance of intersterility barriers in basidiomycetous fungi. The levels of divergence observed in these and in other fungi are significantly higher than that observed in many other eukaryotes, and this suggests that the rate of molecular or morphological evolution in fungi may differ from that found in other organisms.

Phylogenetic relationships of the liverworts (Hepaticae), a basal embryophyte lineage, inferred from nucleotide sequence data of the chloroplast gene rbcL.

Sequence data from the chloroplast-encoded gene rbcL were obtained for 24 liverworts, a basal group of embryophytes. Maximum likelihood and parsimony analyses of these data, along with data from other major green plant lineages, confirm hypotheses based on morphological data, such as the paraphyly of bryophytes, and the basal position of liverworts. Molecular data corroborate the deep separation between the complex thalloid and leafy/simple thalloid liverworts implied by morphological data, but the monophyly of liverworts could not be rejected. The effects of accounting for site-to-site rate heterogeneity in these data were examined using maximum likelihood methods. Comparison of trees obtained with and without rate heterogeneity showed that simply allowing for heterogeneity had a greater improvement on likelihood score than optimization of transition/transversion bias. Incorporation of site-to-site rate heterogeneity in the larger analysis, however, did not necessarily change which topology was favored. Properties of rbcL sequences from the two liverwort groups were compared. Significantly different substitution rates were found between leafy/simple thalloid and complex thalloid liverwort taxa, with rates of rbcL sequence evolution in leafy/simple thalloid taxa being higher and more indicative of those of vascular plants, and with those of complex thalloid taxa (such as Marchantia) being slower. Codon usage in rbcL in complex thalloid liverworts was biased toward NNU and NNA, compared to the leafy/simple thalloid liverworts. Although base composition and relative substitution rates differed between the two groups, no significant differences were detected within each of the two groups of liverworts. The signal present in first and second codon sites versus third codon sites was compared. While the third codon positions in rbcL across this taxon sampling are highly variable (with only 15 constant sites of 439), the trees obtained were in general agreement with trees from the entire data set and with trees obtained from independent sources of data. The presence of signal in third codon positions across greater than 400 MY of plant evolution means that definitions of saturation based on pair-wise comparisons of sequences inadequately assess phylogenetic signal.

PCR-restriction fragment length polymorphism (RFLP) analyses reveal both extensive clonality and local genetic differences in Candida albicans.

Using a polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method to obtain genotypes for the diploid pathogenic yeast, Candida albicans, we analysed 204 C. albicans isolates from three populations of the Duke University community: two from clinical sources [one from patients infected with human immunodeficiency virus (HIV) and the other from patients without HIV infection], and the third from healthy student volunteers. The results indicated: (i) extensive evidence for clonality within and between populations of C. albicans; and (ii) greater genotypic and gene diversities in the nonclinical population than those derived from clinical specimens, regardless of HIV status. The two clinical populations were genetically more similar to each other than either was to the population consisting of isolates from healthy people. Within each population sample there was a general lack of heterozygotes, and random associations of alleles within and between loci were found in less than 50% of the loci or pairs of loci. These findings were consistent between the two sets of samples analysed: those including all isolates and those including only clone-corrected isolates. Possible mechanisms are presented to explain the observed patterns of genetic variation within and between C. albicans populations.

Lack of genetic differentiation between two geographically diverse samples of Candida albicans isolated from patients infected with human immunodeficiency virus.

The patterns of genetic variation of samples of Candida albicans isolated from patients infected with human immunodeficiency virus in Durham, N.C., and Vitória, Brazil, were compared. Genotypes for 126 strains were obtained at 16 polymorphic restriction sites distributed on nine PCR fragments. The results indicated evidence of clonality both within and between these two geographically diverse samples. The samples are genetically very similar, with little evidence of genetic differentiation.

Analysis of genes coding for small-subunit rRNA sequences in studying phylogenetics of dematiaceous fungal pathogens.

Because of their ability to display yeast-like growth forms in various environmental conditions, dematiaceous (melanized) hyphomycetes of the form-genera Exophiala, Rhinocladiella, and Wangiella have been informally termed "black yeasts." Cladistic analysis of 1,050 bp of the genes coding for small-subunit rRNA (SSU rDNA) supported a close relationship among species of these black yeasts with other dematiaceous hyphomycetes in the form-genera Fonsecaea, Phialophora, and Ramichloridium. The conventional categories of these fungi based on asexual states are not supported by phylogenetic analysis of SSU rDNA sequences. Isolates exhibiting annellidic modes of blastic conidiogenesis (e.g., Exophiala spp.) were not monophyletic and were placed as sister taxa to isolates that produce phialides or sympodulae. The results indicated very close relationships between isolates of Wangiella dermatitidis and Exophiala mansonii and between Rhinocladiella aquaspersa and Exophiala jeanselmei. This clade of dematiaceous hyphomycetes was a sister group to a clade comprising members of two orders of cleistothecial ascomycetes, Eurotiales and Onygenales. The etiological agents of chromoblastomycosis were found to be a closely related group (clade), while the agents of phaeohyphomycosis displayed a broader distribution on the SSU rDNA tree.

Development and characterization of a genetic linkage map of Cryptococcus neoformans var. neoformans using amplified fragment length polymorphisms and other markers.

A segregating population of single basidiospore isolates from a sexual cross was used to generate the first moderately dense genetic linkage map of Cryptococcus neoformans var. neoformans (Serotype D). Polymorphic DNA markers were developed using amplified fragment length polymorphisms, random amplified polymorphic DNA, and gene-encoding sequences. These markers were used to analyze 100 meiotic progeny. All markers were tested for distorted segregation with a goodness of fit test. Of the total of 181 markers, 148 showed balanced (1:1) segregation ratios. Segregation distortion was observed for 33 markers. Based on all the markers, a linkage map was generated that consists of 14 major linkage groups with 127 markers, several small linkage groups, and 2 linkage groups that consist only of highly skewed markers. The genetic distance of the linkage map is 1356.3 cM. The estimated total haploid genome size for C. neoformans var. neoformans was calculated using Hulberts method and yielded a map size of 1917 cM. The number of major linkage groups correlates well with the proposed number of 13 chromosomes for C. neoformans var. neoformans. Several genes, including CAP64, CnLAC, and the mating-type locus, were mapped, and their associations were consistent with published data. To date, 6 linkage groups have been assigned to their corresponding chromosomes. This linkage map should provide a framework for the ongoing genome sequencing project and will be a useful tool for studying the genetics and pathogenicity of this important medical yeast.

Colony size can be used to determine the MIC of fluconazole for pathogenic yeasts.

This report describes a new statistical method for estimating the MIC of fluconazole for yeasts pathogenic to humans. This method is based on comparison of the colony sizes on solid media containing different concentrations of fluconazole. By this method, the MICs of fluconazole for 10 yeast strains were comparable to results obtained by the standard method recommended by the National Committee for Clinical Laboratory Standards. This method is simple to perform and easy to interpret. The turnaround time is faster than other methods. The method should be applicable to the determination MICs of other antifungal drugs for yeasts.