Rate and mode differences between nuclear and mitochondrial small-subunit rRNA genes in mushrooms

Species boundaries in the human pathogen Paracoccidioides

The use of molecular taxonomy for identifying recently diverged species has transformed the study of speciation in fungi. The pathogenic fungus Paracoccidioides spp has been hypothesized to be composed of five phylogenetic species, four of which compose the brasiliensis species complex. Nuclear gene genealogies support this divergence scenario, but mitochondrial loci do not; while all species from the brasiliensis complex are differentiated at nuclear coding loci, they are not at mitochondrial loci. We addressed the source of this incongruity using 11 previously published gene fragments, 10 newly-sequenced nuclear non-coding loci, and 10 microsatellites. We hypothesized and further demonstrated that the mito-nuclear incongruence in the brasiliensis species complex results from interspecific hybridization and mitochondrial introgression, a common phenomenon in eukaryotes. Additional population genetic analyses revealed possible nuclear introgression but much less than that seen in the mitochondrion. Our results are consistent with a divergence scenario of secondary contact and subsequent mitochondrial introgression despite the continued persistence of species boundaries. We also suggest that yeast morphology slightly-but significantly-differs across all five Paracoccidioides species and propose to elevate four of these phylogenetic species to formally described taxonomic species.

The evolution of sex: A new hypothesis based on mitochondrial mutational erosion: Mitochondrial mutational erosion in ancestral eukaryotes would favor the evolution of sex, harnessing nuclear recombination to optimize compensatory nuclear coadaptation

The evolution of sex in eukaryotes represents a paradox, given the "twofold" fitness cost it incurs. We hypothesize that the mutational dynamics of the mitochondrial genome would have favored the evolution of sexual reproduction. Mitochondrial DNA (mtDNA) exhibits a high-mutation rate across most eukaryote taxa, and several lines of evidence suggest that this high rate is an ancestral character. This seems inexplicable given that mtDNA-encoded genes underlie the expression of life's most salient functions, including energy conversion. We propose that negative metabolic effects linked to mitochondrial mutation accumulation would have invoked selection for sexual recombination between divergent host nuclear genomes in early eukaryote lineages. This would provide a mechanism by which recombinant host genotypes could be rapidly shuffled and screened for the presence of compensatory modifiers that offset mtDNA-induced harm. Under this hypothesis, recombination provides the genetic variation necessary for compensatory nuclear coadaptation to keep pace with mitochondrial mutation accumulation.

Whole genome and global gene expression analyses of the model mushroom Flammulina velutipes reveal a high capacity for lignocellulose degradation

Flammulina velutipes is a fungus with health and medicinal benefits that has been used for consumption and cultivation in East Asia. F. velutipes is also known to degrade lignocellulose and produce ethanol. The overlapping interests of mushroom production and wood bioconversion make F. velutipes an attractive new model for fungal wood related studies. Here, we present the complete sequence of the F. velutipes genome. This is the first sequenced genome for a commercially produced edible mushroom that also degrades wood. The 35.6-Mb genome contained 12,218 predicted protein-encoding genes and 287 tRNA genes assembled into 11 scaffolds corresponding with the 11 chromosomes of strain KACC42780. The 88.4-kb mitochondrial genome contained 35 genes. Well-developed wood degrading machinery with strong potential for lignin degradation (69 auxiliary activities, formerly FOLymes) and carbohydrate degradation (392 CAZymes), along with 58 alcohol dehydrogenase genes were highly expressed in the mycelium, demonstrating the potential application of this organism to bioethanol production. Thus, the newly uncovered wood degrading capacity and sequential nature of this process in F. velutipes, offer interesting possibilities for more detailed studies on either lignin or (hemi-) cellulose degradation in complex wood substrates. The mutual interest in wood degradation by the mushroom industry and (ligno-)cellulose biomass related industries further increase the significance of F. velutipes as a new model.

Finding single copy genes out of sequenced genomes for multilocus phylogenetics in non-model fungi

Historically, fungal multigene phylogenies have been reconstructed based on a small number of commonly used genes. The availability of complete fungal genomes has given rise to a new wave of model organisms that provide large number of genes potentially useful for building robust gene genealogies. Unfortunately, cross-utilization of these resources to study phylogenetic relationships in the vast majority of non-model fungi (i.e. “orphan” species) remains an unexamined question. To address this problem, we developed a method coupled with a program named “PHYLORPH” (PHYLogenetic markers for ORPHans). The method screens fungal genomic databases (107 fungal genomes fully sequenced) for single copy genes that might be easily transferable and well suited for studies at low taxonomic levels (for example, in species complexes) in non-model fungal species. To maximize the chance to target genes with informative regions, PHYLORPH displays a graphical evaluation system based on the estimation of nucleotide divergence relative to substitution type. The usefulness of this approach was tested by developing markers in four non-model groups of fungal pathogens. For each pathogen considered, 7 to 40% of the 10–15 best candidate genes proposed by PHYLORPH yielded sequencing success. Levels of polymorphism of these genes were compared with those obtained for some genes traditionally used to build fungal phylogenies (e.g. nuclear rDNA, β-tubulin, γ-actin, Elongation factor EF-1α). These genes were ranked among the best-performing ones and resolved accurately taxa relationships in each of the four non-model groups of fungi considered. We envision that PHYLORPH will constitute a useful tool for obtaining new and accurate phylogenetic markers to resolve relationships between closely related non-model fungal species.

Loss of the flagellum happened only once in the fungal lineage: phylogenetic structure of kingdom Fungi inferred from RNA polymerase II subunit genes

Background

At present, there is not a widely accepted consensus view regarding the phylogenetic structure of kingdom Fungi although two major phyla, Ascomycota and Basidiomycota, are clearly delineated. Regarding the lower fungi, Zygomycota and Chytridiomycota, a variety of proposals have been advanced. Microsporidia may or may not be fungi; the Glomales (vesicular-arbuscular mycorrhizal fungi) may or may not constitute a fifth fungal phylum, and the loss of the flagellum may have occurred either once or multiple times during fungal evolution. All of these issues are capable of being resolved by a molecular phylogenetic analysis which achieves strong statistical support for major branches. To date, no fungal phylogeny based upon molecular characters has satisfied this criterion.

Results

Using the translated amino acid sequences of the RPB1 and RPB2 genes, we have inferred a fungal phylogeny that consists largely of well-supported monophyletic phyla. Our major results, each with significant statistical support, are: (1) Microsporidia are sister to kingdom Fungi and are not members of Zygomycota; that is, Microsporidia and fungi originated from a common ancestor. (2) Chytridiomycota, the only fungal phylum having a developmental stage with a flagellum, is paraphyletic and is the basal lineage. (3) Zygomycota is monophyletic based upon sampling of Trichomycetes, Zygomycetes, and Glomales. (4) Zygomycota, Basidiomycota, and Ascomycota form a monophyletic group separate from Chytridiomycota. (5) Basidiomycota and Ascomycota are monophyletic sister groups.

Conclusion

In general, this paper highlights the evolutionary position and significance of the lower fungi (Zygomycota and Chytridiomycota). Our results suggest that loss of the flagellum happened only once during early stages of fungal evolution; consequently, the majority of fungi, unlike plants and animals, are nonflagellated. The phylogeny we infer from gene sequences is the first one that is congruent with the widely accepted morphology-based classification of Fungi. We find that, contrary to what has been published elsewhere, the four morphologically defined phyla (Ascomycota, Basidiomycota, Zygomycota and Chytridiomycota) do not overlap with one another. Microsporidia are not included within kingdom Fungi; rather they are a sister-group to the Fungi. Our study demonstrates the applicability of protein sequences from large, slowly-evolving genes to the derivation of well-resolved and highly supported phylogenies across long evolutionary distances.

Mitochondrial substitution rates are extraordinarily elevated and variable in a genus of flowering plants

Plant mitochondrial (mt) genomes have long been known to evolve slowly in sequence. Here we show remarkable departure from this pattern of conservative evolution in a genus of flowering plants. Substitution rates at synonymous sites vary substantially among lineages within Plantago. At the extreme, rates in Plantago exceed those in exceptionally slow plant lineages by approximately 4,000-fold. The fastest Plantago lineages set a new benchmark for rapid evolution in a DNA genome, exceeding even the fastest animal mt genome by an order of magnitude. All six mt genes examined show similarly elevated divergence in Plantago, implying that substitution rates are highly accelerated throughout the genome. In contrast, substitution rates show little or no elevation in Plantago for each of four chloroplast and three nuclear genes examined. These results, combined with relatively modest elevations in rates of nonsynonymous substitutions in Plantago mt genes, indicate that major, reversible changes in the mt mutation rate probably underlie the extensive variation in synonymous substitution rates. These rate changes could be caused by major changes in any number of factors that control the mt mutation rate, from the production and detoxification of oxygen free radicals in the mitochondrion to the efficacy of mt DNA replication and/or repair.

Secondary structure and molecular evolution of the mitochondrial small subunit ribosomal RNA in Agaricales (Euagarics clade, Homobasidiomycota)

The complete sequences and secondary structures of the mitochondrial small subunit (SSU) ribosomal RNAs of both mostly cultivated mushrooms Agaricus bisporus (1930 nt) and Lentinula edodes (2164 nt) were achieved. These secondary structures and that of Schizophyllum commune (1872 nt) were compared to that previously established for Agrocybe aegerita. The four structures are near the model established for Archae, Bacteria, plastids, and mitochondria; particularly the helices 23 and 37, described as specific to bacteria, are present. Within the four Agaricales (Homobasidiomycota), the SSU-rRNA "core" is conserved in size (966 to 1009 nt) with the exception of an unusual extension of 40 nt in the H17 helix of S. commune. The four core sequences possess 76% of conserved positions and a cluster of C in their 3' end, which could constitute a signal involved in the RNA maturation process. Among the nine putative variable domains, three (V3, V5, V7) do not show significant length variations and possess similar percentages of conserved positions (69%) than the core. The other six variable domains show important length variations, due to independent large size inserted/deleted sequences, and higher rates of nucleotide substitutions than the core (only 31% of conserved positions between the four species). Interestingly, the inserted/deleted sequences are located in few preferential sites (hot spots for insertion/deletion) where they seem to arise or disappear haphazardly during evolution. These sites are located on the surface of the tertiary structure of the 30S ribosomal subunit, at the beginning of hairpin loops; the insertions lead to a lengthening of existing hairpins or to branching loops bearing up to five additional helices.

Genetic relatedness among Filobasidiella species

The three accepted species of Filobasidiella, F. neoformans, F. depauperata, and F. lutea, are compared morphologically and by molecular analysis. Sequences of the internally transcribed spacer (ITS) and the small subunit (SSU) gene of the ribosomal RNA (rRNA) gene cluster were obtained, and analysed by Neighbor-joining and Maximum parsimony methods. The three species of Filobsidiella are shown to form a single monophyletic clade, rooted by Tremella mesenterica. F. lutea was recovered as a distinct, but closely related taxon with the Filobasidiella clade. This is the first report of DNA sequences from herbarium specimens of F. lutea.

Evolution of complex fruiting-body morphologies in homobasidiomycetes

The fruiting bodies of homobasidiomycetes include some of the most complex forms that have evolved in the fungi, such as gilled mushrooms, bracket fungi and puffballs ('pileate-erect') forms. Homobasidiomycetes also include relatively simple crust-like 'resupinate' forms, however, which account for ca. 13-15% of the described species in the group. Resupinate homobasidiomycetes have been interpreted either as a paraphyletic grade of plesiomorphic forms or a polyphyletic assemblage of reduced forms. The former view suggests that morphological evolution in homobasidiomycetes has been marked by independent elaboration in many clades, whereas the latter view suggests that parallel simplification has been a common mode of evolution. To infer patterns of morphological evolution in homobasidiomycetes, we constructed phylogenetic trees from a dataset of 481 species and performed ancestral state reconstruction (ASR) using parsimony and maximum likelihood (ML) methods. ASR with both parsimony and ML implies that the ancestor of the homobasidiomycetes was resupinate, and that there have been multiple gains and losses of complex forms in the homobasidiomycetes. We also used ML to address whether there is an asymmetry in the rate of transformations between simple and complex forms. Models of morphological evolution inferred with ML indicate that the rate of transformations from simple to complex forms is about three to six times greater than the rate of transformations in the reverse direction. A null model of morphological evolution, in which there is no asymmetry in transformation rates, was rejected. These results suggest that there is a 'driven' trend towards the evolution of complex forms in homobasidiomycetes.

Sequence and secondary structure of the mitochondrial small-subunit rRNA V4, V6, and V9 domains reveal highly species-specific variations within the genus Agrocybe

A comparative study of variable domains V4, V6, and V9 of the mitochondrial small-subunit (SSU) rRNA was carried out with the genus Agrocybe by PCR amplification of 42 wild isolates belonging to 10 species, Agrocybe aegerita, Agrocybe dura, Agrocybe chaxingu, Agrocybe erebia, Agrocybe firma, Agrocybe praecox, Agrocybe paludosa, Agrocybe pediades, Agrocybe alnetorum, and Agrocybe vervacti. Sequencing of the PCR products showed that the three domains in the isolates belonging to the same species were the same length and had the same sequence, while variations were found among the 10 species. Alignment of the sequences showed that nucleotide motifs encountered in the smallest sequence of each variable domain were also found in the largest sequence, indicating that the sequences evolved by insertion-deletion events. Determination of the secondary structure of each domain revealed that the insertion-deletion events commonly occurred in regions not directly involved in the secondary structure (i.e., the loops). Moreover, conserved sequences ranging from 4 to 25 nucleotides long were found at the beginning and end of each domain and could constitute genus-specific sequences. Comparisons of the V4, V6, and V9 secondary structures resulted in identification of the following four groups: (i) group I, which was characterized by the presence of additional P23-1 and P23-3 helices in the V4 domain and the lack of the P49-1 helix in V9 and included A. aegerita, A. chaxingu, and A. erebia; (ii) group II, which had the P23-3 helix in V4 and the P49-1 helix in V9 and included A. pediades; (iii) group III, which did not have additional helices in V4, had the P49-1 helix in V9 and included A. paludosa, A. firma, A. alnetorum, and A. praecox; and (iv) group IV, which lacked both the V4 additional helices and the P49-1 helix in V9 and included A. vervacti and A. dura. This grouping of species was supported by the structure of a consensus tree based on the variable domain sequences. The conservation of the sequences of the V4, V6, and V9 domains of the mitochondrial SSU rRNA within species and the high degree of interspecific variation found in the Agrocybe species studied open the way for these sequences to be used as specific molecular markers of the Basidiomycota.

Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the "euagarics." Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.

Mutational spectrometry without phenotypic selection: human mitochondrial DNA

By first separating mutant from nonmutant DNA sequences on the basis of their melting temperatures and then increasing the number of copies by high-fidelity DNA amplification, we have developed a method that allows observation of point mutations in biological samples at fractions at or above 10-6. Using this method, we have observed the hotspot point mutations that lie in 100 base pairs of the mitochondrial genome in samples of cultured cells and human tissues. To date, 19 mutants have been isolated, their fractions ranging from 4x10-4 down to the limit of detection. We performed specific tests to determine if the observed signals were artefacts arising from contamination, polymerase errors during PCR or DNA adducts created during the procedure. We also tested the possibilities that DNA replication mismatch intermediates, or endogenous DNA adducts that were originally present in the cells, were included with true mutants in our separation steps and converted to mutants during PCR. We show that while most of the mutants behave as double-stranded point mutants in the cells, some appear to arise at least in part from mismatch intermediates or cellular DNA adducts. This technology is therefore sufficient for the observation of the spectrum of point mutations in human mitochondrial DNA and is a tool for discovering the primary causes of these mutations.

The RecA protein as a model molecule for molecular systematic studies of bacteria: comparison of trees of RecAs and 16S rRNAs from the same species

The evolution of the RecA protein was analyzed using molecular phylogenetic techniques. Phylogenetic trees of all currently available complete RecA proteins were inferred using multiple maximum parsimony and distance matrix methods. Comparison and analysis of the trees reveal that the inferred relationships among these proteins are highly robust. The RecA trees show consistent subdivisions corresponding to many of the major bacterial groups found in trees of other molecules including the alpha, beta, gamma, delta, epsilon proteobacteria, cyanobacteria, high-GC gram-positives, and the Deinococcus-Thermus group. However, there are interesting differences between the RecA trees and these other trees. For example, in all the RecA trees the proteins from gram-positive species are not monophyletic. In addition, the RecAs of the cyanobacteria consistently group with those of the high-GC gram-positives. To evaluate possible causes and implications of these and other differences phylogenetic trees were generated for small-subunit rRNA sequences from the same (or closely related) species as represented in the RecA analysis. The trees of the two molecules using these equivalent species-sets are highly congruent and have similar resolving power for close, medium, and deep branches in the history of bacteria. The implications of the particular similarities and differences between the trees are discussed. Some of the features that make RecA useful for molecular systematics and for studies of protein evolution are also discussed.

Mitochondrial DNAs of the fungus Armillaria ostoyae: restriction map and length variation

A restriction-enzyme-site map is presented for the 147-kb mtDNA of North American Armillaria ostoyae. The locations of five structural genes, atp6, atp8, coxI, coxIII, and cob, along with the location and orientation of the large and small ribosomal RNA genes, were determined through Southern hybridizations with cloned genes from other fungal mtDNAs. Based on this map, the variation in mtDNA suggested geographic structure at two different levels. On a large geographic scale, 17 mtDNA types from North America were distinct, with respect to both size and restriction maps, from three mtDNA types from Europe. At the local scale, identical mtDNA types were evident among several different genetic individuals located no more than 1 km apart at a site in Michigan. No mtDNA type occurred more than once among genetic individuals from different regions of North America, although the occurrence of similar mtDNAs in isolates from distant regions suggested that this may occur at a low frequency with large sample sizes. Among the North American mtDNA types, analysis of discrete length variants was inconsistent with the hypothesis that the mtDNA of A. ostoyae evolves as a clonal lineage in which each length mutation represents a unique event. The two remaining hypotheses, that similar mutational events have occurred independently and that genetic exchange and recombination occurs among mtDNAs in natural populations of this species, remain to be tested.