Trichoderma koningiopsis (Hypocreaceae) has the smallest mitogenome of the genus Trichoderma

Introduction

Fungal mitogenomes exhibit remarkable variation in conformation, size, gene content, arrangement and expression, including their intergenic spacers and introns.

Methods

The complete mitochondrial genome sequence of the mycoparasitic fungus Trichoderma koningiopsis was determined using the Illumina next-generation sequencing technology. We used data from our recent Illumina NGS-based project of T. koningiopsis genome sequencing to study its mitochondrial genome. The mitogenome was assembled, annotated, and compared with other fungal mitogenomes.

Results

T. koningiopsis strain POS7 mitogenome is a circular molecule of 27,560 bp long with a GC content of 27.80%. It harbors the whole complement of the 14 conserved mitochondrial protein-coding genes (PCG) such as atp6, atp8, atp9, cox1, cox2, cox3, cob, nad1, nad2, nad3, nad4, nad4L, nad5, and nad6, also found in the same gene order to other Hypocreales. The mitogenome also contains 26 transfer RNA genes (tRNAs), 5 of them with more than one copy. Other genes also present in the assembled mitochondrial genome are a small rRNA subunit and a large rRNA subunit containing ribosomal protein S3 gene. Despite the small genome size, two introns were detected in the T. koningiopsis POS7 mitogenome, one of them in cox3 gene and the other in rnl gene, accounting 7.34% of this mitogenome with a total size of 2,024 bp. A phylogenetic analysis was done using the 14 PCGs genes of T. koningiopsis strain POS7 mitogenome to compare them with those from other fungi of the Subphyla Pezizomycotina and Saccharomycotina. T. koningiopsis strain POS7 was clustered together with other representatives of Trichoderma lineage, within the Hypocreales group, which is also supported by previous phylogenetic studies based on nuclear markers.

Discussion

The mitochondrial genome of T. koningiopsis POS7 will allow further investigations into the taxonomy, phylogenetics, conservation genetics, and evolutionary biology of this important genus as well as other closely related species.

Genomic and Transcriptomic Insight of Giant Sclerotium Formation of Wood-Decay Fungi

Many fungi form persistent and dormant sclerotia with compact hardened mycelia during unfavorable circumstances. While most of these sclerotia are small in size, Wolfiporia cocos, a wood-decay fungus, grows into giant sclerotia, which are mainly composed of polysaccharides of linear (1→3)-β-D-glucans. To explore the underlying mechanism of converting sophisticated wood polysaccharides for biosynthesis of highly homogenized glucans in W. cocos, we sequenced and assembled the genome of a cultivated W. cocos strain (WCLT) in China. The 62-Mb haploid genome contains 44.2% repeat sequences, of which, 48.0% are transposable elements (TEs). Contrary to the genome of W. cocos from North America, WCLT has independently undergone a partial genome duplication (PGD) event. The large-scale TE insertion and PGD occurrence overlapped with an archeological Pleistocene stage of low oxygen and high temperature, and these stresses might have induced the differences in sclerotium due to geographical distribution. The wood decomposition enzymes, as well as sclerotium-regulator kinases, aquaporins, and highly expanded gene families such as NAD-related families, together with actively expressed 1,3-β-glucan synthase for sclerotium polysaccharides, all have contributed to the sclerotium formation and expansion. This study shall inspire further exploration on how fungi convert wood into simple glucans in the sclerotium of W. cocos.

Comparative mitochondrial genome analysis reveals intron dynamics and gene rearrangements in two Trametes species

Trametes species are efficient wood decomposers that are widespread throughout the world. Mitogenomes have been widely used to understand the phylogeny and evolution of fungi. Up to now, two mitogenomes from the Trametes genus have been revealed. In the present study, the complete mitogenomes of two novel Trametes species, Trametes versicolor and T. coccinea, were assembled and compared with other Polyporales mitogenomes. Both species contained circular DNA molecules, with sizes of 67,318 bp and 99,976 bp, respectively. Comparative mitogenomic analysis indicated that the gene number, length and base composition varied between the four Trametes mitogenomes we tested. In addition, all of the core protein coding genes in Trametes species were identified and subjected to purifying selection. The mitogenome of T. coccinea contained the largest number of introns among the four Trametes species tested, and introns were considered the main factors contributing to size variations of Polyporales. Several novel introns were detected in the Trametes species we assembled, and introns identified in Polyporales were found to undergo frequent loss/gain events. Large-scale gene rearrangements were detected between closely related Trametes species, including gene inversions, insertions, and migrations. A well-supported phylogenetic tree for 77 Basidiomycetes was obtained based on the combined mitochondrial gene set using 2 phylogenetic inference methods. The results showed that mitochondrial genes are effective molecular markers for understanding the phylogeny of Basidiomycetes. This study is the first to report the mitogenome rearrangement and intron dynamics of Trametes species, which shed light on the evolution of Trametes and other related species.

The complete mitochondrial genome of medicinal fungus Taiwanofungus camphoratus reveals gene rearrangements and intron dynamics of Polyporales

Taiwanofungus camphoratus is a highly valued medicinal mushroom that is endemic to Taiwan, China. In the present study, the mitogenome of T. camphoratus was assembled and compared with other published Polyporales mitogenomes. The T. camphoratus mitogenome was composed of circular DNA molecules, with a total size of 114,922 bp. Genome collinearity analysis revealed large-scale gene rearrangements between the mitogenomes of Polyporales, and T. camphoratus contained a unique gene order. The number and classes of introns were highly variable in 12 Polyporales species we examined, which proved that numerous intron loss or gain events occurred in the evolution of Polyporales. The Ka/Ks values for most core protein coding genes in Polyporales species were less than 1, indicating that these genes were subject to purifying selection. However, the rps3 gene was found under positive or relaxed selection between some Polyporales species. Phylogenetic analysis based on the combined mitochondrial gene set obtained a well-supported topology, and T. camphoratus was identified as a sister species to Laetiporus sulphureus. This study served as the first report on the mitogenome in the Taiwanofungus genus, which will provide a basis for understanding the phylogeny and evolution of this important fungus.

Characterization of Two Mitochondrial Genomes and Gene Expression Analysis Reveal Clues for Variations, Evolution, and Large-Sclerotium Formation in Medical Fungus Wolfiporia cocos

Wolfiporia cocos, a precious mushroom with a long history as an edible food and Asian traditional medicine, remains unclear in the genetic mechanism underlying the formation of large sclerotia. Here, two complete circular mitogenomes (BL16, 135,686 bp and MD-104 SS10, 124,842 bp, respectively) were presented in detail first. The salient features in the mitogenomes of W. cocos include an intron in the tRNA (trnQ-UUG²), and an obvious gene rearrangement identified between the two mitogenomes from the widely geographically separated W. cocos strains. Genome comparison and phylogenetic analyses reveal some variations and evolutional characteristics in W. cocos. Whether the mitochondrion is functional in W. cocos sclerotium development was investigated by analyzing the mitogenome synteny of 10 sclerotium-forming fungi and mitochondrial gene expression patterns in different W. cocos sclerotium-developmental stages. Three common homologous genes identified across ten sclerotium-forming fungi were also found to exhibit significant differential expression levels during W. cocos sclerotium development. Most of the mitogenomic genes are not expressed in the mycelial stage but highly expressed in the sclerotium initial or developmental stage. These results indicate that some of mitochondrial genes may play a role in the development of sclerotium in W. cocos, which needs to be further elucidated in future studies. This study will stimulate new ideas on cytoplasmic inheritance of W. cocos and facilitate the research on the role of mitochondria in large sclerotium formation.

Comparison of the Mitochondrial Genome Sequences of Six Annulohypoxylon stygium Isolates Suggests Short Fragment Insertions as a Potential Factor Leading to Larger Genomic Size

Mitochondrial DNA (mtDNA) is a core non-nuclear genetic material found in all eukaryotic organisms, the size of which varies extensively in the eumycota, even within species. In this study, mitochondrial genomes of six isolates of Annulohypoxylon stygium (Lév.) were assembled from raw reads from PacBio and Illumina sequencing. The diversity of genomic structures, conserved genes, intergenic regions and introns were analyzed and compared. Genome sizes ranged from 132 to 147 kb and contained the same sets of conserved protein-coding, tRNA and rRNA genes and shared the same gene arrangements and orientation. In addition, most intergenic regions were homogeneous and had similar sizes except for the region between cytochrome b (cob) and cytochrome c oxidase I (cox1) genes which ranged from 2,998 to 8,039 bp among the six isolates. Sixty-five intron insertion sites and 99 different introns were detected in these genomes. Each genome contained 45 or more introns, which varied in distribution and content. Introns from homologous insertion sites also showed high diversity in size, type and content. Comparison of introns at the same loci showed some complex introns, such as twintrons and ORF-less introns. There were 44 short fragment insertions detected within introns, intergenic regions, or as introns, some of them located at conserved domain regions of homing endonuclease genes. Insertions of short fragments such as small inverted repeats might affect or hinder the movement of introns, and these allowed for intron accumulation in the mitochondrial genomes analyzed, and enlarged their size. This study showed that the evolution of fungal mitochondrial introns is complex, and the results suggest short fragment insertions as a potential factor leading to larger mitochondrial genomes in A. stygium.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

Whole genome sequencing of Rhodotorula mucilaginosa isolated from the chewing stick (Distemonanthus benthamianus): insights into Rhodotorula phylogeny, mitogenome dynamics and carotenoid biosynthesis

In industry, the yeast Rhodotorula mucilaginosa is commonly used for the production of carotenoids. The production of carotenoids is important because they are used as natural colorants in food and some carotenoids are precursors of retinol (vitamin A). However, the identification and molecular characterization of the carotenoid pathway/s in species belonging to the genus Rhodotorula is scarce due to the lack of genomic information thus potentially impeding effective metabolic engineering of these yeast strains for improved carotenoid production. In this study, we report the isolation, identification, characterization and the whole nuclear genome and mitogenome sequence of the endophyte R. mucilaginosa RIT389 isolated from Distemonanthus benthamianus, a plant known for its anti-fungal and antibacterial properties and commonly used as chewing sticks. The assembled genome of R. mucilaginosa RIT389 is 19 Mbp in length with an estimated genomic heterozygosity of 9.29%. Whole genome phylogeny supports the species designation of strain RIT389 within the genus in addition to supporting the monophyly of the currently sequenced Rhodotorula species. Further, we report for the first time, the recovery of the complete mitochondrial genome of R. mucilaginosa using the genome skimming approach. The assembled mitogenome is at least 7,000 bases larger than that of Rhodotorula taiwanensis which is largely attributed to the presence of large intronic regions containing open reading frames coding for homing endonuclease from the LAGLIDADG and GIY-YIG families. Furthermore, genomic regions containing the key genes for carotenoid production were identified in R. mucilaginosa RIT389, revealing differences in gene synteny that may play a role in the regulation of the biotechnologically important carotenoid synthesis pathways in yeasts.

Escaping introns in COI through cDNA barcoding of mushrooms: Pleurotus as a test case

DNA barcoding involves the use of one or more short, standardized DNA fragments for the rapid identification of species. A 648-bp segment near the 5' terminus of the mitochondrial cytochrome c oxidase subunit I (COI) gene has been adopted as the universal DNA barcode for members of the animal kingdom, but its utility in mushrooms is complicated by the frequent occurrence of large introns. As a consequence, ITS has been adopted as the standard DNA barcode marker for mushrooms despite several shortcomings. This study employed newly designed primers coupled with cDNA analysis to examine COI sequence diversity in six species of Pleurotus and compared these results with those for ITS. The ability of the COI gene to discriminate six species of Pleurotus, the commonly cultivated oyster mushroom, was examined by analysis of cDNA. The amplification success, sequence variation within and among species, and the ability to design effective primers was tested. We compared ITS sequences to their COI cDNA counterparts for all isolates. ITS discriminated between all six species, but some sequence results were uninterpretable, because of length variation among ITS copies. By comparison, a complete COI sequences were recovered from all but three individuals of Pleurotus giganteus where only the 5' region was obtained. The COI sequences permitted the resolution of all species when partial data was excluded for P. giganteus. Our results suggest that COI can be a useful barcode marker for mushrooms when cDNA analysis is adopted, permitting identifications in cases where ITS cannot be recovered or where it offers higher resolution when fresh tissue is. The suitability of this approach remains to be confirmed for other mushrooms.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Complete mitochondrial genome of the ectomycorrhizal fungus Tricholoma matsutake

Tricholoma matsutake is an ectomycorrhizal fungus belonging to the class Agaricomycetes and the phylum Basidiomycota. Here, we report the complete mitochondrial genome sequence of T. matsutake. The mitochondrial DNA (mtDNA) of T. matsutake seems to be concatenated, and harbors 76,037 nucleotides, with an overall GC content of 20.69%. The mtDNA comprises 15 protein-coding genes, 28 tRNA genes, and 2 rRNA genes. A comparative mitogenomic analysis between Agaricomycetes and other classes (Exobasidiomycetes, Pucciniomycetes, Tremellomycetes, and Ustilaginomycetes) revealed that species belonging to the class Agaricomycetes preferably harbor the codon UUA (coding for leucine) in their mtDNAs, whereas other classes prefer the codon CUA. This bias in the usage of the leucine codon could be the result of evolutionary divergence between the Agaricomycetes and other classes. The T. matsutake mtDNA sequence thus provides insight into the evolutionary characteristics of the Agaricomycetes mtDNAs.

Intronic and plasmid-derived regions contribute to the large mitochondrial genome sizes of Agaricomycetes

Sizes of mitochondrial genomes vary extensively between fungal species although they typically contain a conserved set of core genes. We have characterised the mitochondrial genome of the conifer root rot pathogen Heterobasidion irregulare and compared the size, gene content and structure of 20 Basidiomycete mitochondrial genomes. The mitochondrial genome of H. irregulare was 114, 193 bp and contained a core set of 15 protein coding genes, two rRNA genes and 26 tRNA genes. In addition, we found six non-conserved open reading frames (ORFs) and four putative plasmid genes clustered in three separate regions together with 24 introns and 14 intronic homing endonuclease genes, unequally spread across seven of the core genes. The size differences among the 20 Basidiomycetes can largely be explained by length variation of intergenic regions and introns. The Agaricomycetes contained the nine largest mitochondrial genomes in the Basidiomycete group and Agaricomycete genomes are significantly (p < 0.001) larger than the other Basidiomycetes. A feature of the Agaricomycete mitochondrial genomes in this study was the simultaneous occurrence of putative plasmid genes and non-conserved ORFs, with Cantharellus cibarius as only exception, where no non-conserved ORF was identified. This indicates a mitochondrial plasmid origin of the non-conserved ORFs or increased mitochondrial genome dynamics of species harbouring mitochondrial plasmids. We hypothesise that two independent factors are the driving forces for large mitochondrial genomes: the homing endonuclease genes in introns and integration of plasmid DNA.

Mitochondrial genome of Phlebia radiata is the second largest (156 kbp) among fungi and features signs of genome flexibility and recent recombination events

Mitochondria are eukaryotic organelles supporting individual life-style via generation of proton motive force and cellular energy, and indispensable metabolic pathways. As part of genome sequencing of the white rot Basidiomycota species Phlebia radiata, we first assembled its mitochondrial genome (mtDNA). So far, the 156 348 bp mtDNA is the second largest described for fungi, and of considerable size among eukaryotes. The P. radiata mtDNA assembled as single circular dsDNA molecule containing genes for the large and small ribosomal RNAs, 28 transfer RNAs, and over 100 open reading frames encoding the 14 fungal conserved protein subunits of the mitochondrial complexes I, III, IV, and V. Two genes (atp6 and tRNA-IleGAU) were duplicated within 6.1 kbp inverted region, which is a unique feature of the genome. The large mtDNA size, however, is explained by the dominance of intronic and intergenic regions (sum 80% of mtDNA sequence). The intergenic DNA stretches harness short (≤ 200 nt) repetitive, dispersed and overlapping sequence elements in abundance. Long self-splicing introns of types I and II interrupt eleven of the conserved genes (cox1,2,3; cob; nad1,2,4,4L,5; rnl; rns). The introns embrace a total of 57 homing endonucleases with LAGLIDADGD and GYI-YIG core motifs, which makes P. radiata mtDNA to one of the largest known reservoirs of intron-homing endonucleases. The inverted duplication, intergenic stretches, and intronic features are indications of dynamics and genetic flexibility of the mtDNA, not fully recognized to this extent in fungal mitochondrial genomes previously, thus giving new insights for the evolution of organelle genomes in eukaryotes.

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.

Mitochondrial genome of the basidiomycetous yeast Jaminaea angkorensis

Jaminaea angkorensis is an anamorphic basidiomycetous yeast species originally isolated from decaying leaves in Cambodia. Taxonomically, J. angkorensis is affiliated with Microstromatales (Exobasidiomycetes, Ustilaginomycotina, Basidiomycota) and represents a basal phylogenetic lineage of this fungal order. To perform a comparative analysis of J. angkorensis with other basidiomycetes, we determined and analyzed its complete mitochondrial DNA sequence. The mitochondrial genome is represented by 29,999 base pairs long, circular DNA containing 32 % guanine and cytosine residues. Its genetic organization is relatively compact and comprises typical genes for 15 conserved proteins involved in oxidative phosphorylation (atp6, 8, and 9; cob; cox1, 2, and 3; and nad1, 2, 3, 4, 4L, 5, and 6) and translation (rps3), two ribosomal RNAs (rnl and rns) and twenty-two transfer RNAs (trnA-Y). Although the gene content is similar to other basidiomycetes, the gene orders in the examined species exhibit only a limited synteny, reflecting their phylogenetic distances and extensive genome rearrangements. In addition, a comparative analysis of basidiomycete mitochondrial genomes indicates that stop-to-tryptophan reassignment of the UGA codon was accompanied by structural alterations of tRNA-Trp(CCA). These results provide an insight into the evolution of the genetic code in fungal mitochondria.

Evidence for inter-specific recombination among the mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex

Background

The availability of mitochondrial genomes has allowed for the resolution of numerous questions regarding the evolutionary history of fungi and other eukaryotes. In the Gibberella fujikuroi species complex, the exact relationships among the so-called “African”, “Asian” and “American” Clades remain largely unresolved, irrespective of the markers employed. In this study, we considered the feasibility of using mitochondrial genes to infer the phylogenetic relationships among Fusarium species in this complex. The mitochondrial genomes of representatives of the three Clades (Fusarium circinatum, F. verticillioides and F. fujikuroi) were characterized and we determined whether or not the mitochondrial genomes of these fungi have value in resolving the higher level evolutionary relationships in the complex.

Results

Overall, the mitochondrial genomes of the three species displayed a high degree of synteny, with all the genes (protein coding genes, unique ORFs, ribosomal RNA and tRNA genes) in identical order and orientation, as well as introns that share similar positions within genes. The intergenic regions and introns generally contributed significantly to the size differences and diversity observed among these genomes. Phylogenetic analysis of the concatenated protein-coding dataset separated members of the Gibberella fujikuroi complex from other Fusarium species and suggested that F. fujikuroi (“Asian” Clade) is basal in the complex. However, individual mitochondrial gene trees were largely incongruent with one another and with the concatenated gene tree, because six distinct phylogenetic trees were recovered from the various single gene datasets.

Conclusion

The mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex are remarkably similar to those of the previously characterized Fusarium species and Sordariomycetes. Despite apparently representing a single replicative unit, all of the genes encoded on the mitochondrial genomes of these fungi do not share the same evolutionary history. This incongruence could be due to biased selection on some genes or recombination among mitochondrial genomes. The results thus suggest that the use of individual mitochondrial genes for phylogenetic inference could mask the true relationships between species in this complex.

Complete mitochondrial genome of the medicinal mushroom Ganoderma lucidum

Ganoderma lucidum is one of the well-known medicinal basidiomycetes worldwide. The mitochondrion, referred to as the second genome, is an organelle found in most eukaryotic cells and participates in critical cellular functions. Elucidating the structure and function of this genome is important to understand completely the genetic contents of G. lucidum. In this study, we assembled the mitochondrial genome of G. lucidum and analyzed the differential expressions of its encoded genes across three developmental stages. The mitochondrial genome is a typical circular DNA molecule of 60,630 bp with a GC content of 26.67%. Genome annotation identified genes that encode 15 conserved proteins, 27 tRNAs, small and large rRNAs, four homing endonucleases, and two hypothetical proteins. Except for genes encoding trnW and two hypothetical proteins, all genes were located on the positive strand. For the repeat structure analysis, eight forward, two inverted, and three tandem repeats were detected. A pair of fragments with a total length around 5.5 kb was found in both the nuclear and mitochondrial genomes, which suggests the possible transfer of DNA sequences between two genomes. RNA-Seq data for samples derived from three stages, namely, mycelia, primordia, and fruiting bodies, were mapped to the mitochondrial genome and qualified. The protein-coding genes were expressed higher in mycelia or primordial stages compared with those in the fruiting bodies. The rRNA abundances were significantly higher in all three stages. Two regions were transcribed but did not contain any identified protein or tRNA genes. Furthermore, three RNA-editing sites were detected. Genome synteny analysis showed that significant genome rearrangements occurred in the mitochondrial genomes. This study provides valuable information on the gene contents of the mitochondrial genome and their differential expressions at various developmental stages of G. lucidum. The results contribute to the understanding of the functions and evolution of fungal mitochondrial DNA.

Establishment and interpretation of the genome sequence of the phytopathogenic fungus Rhizoctonia solani AG1-IB isolate 7/3/14

Anastomosis group AG1-IB isolates of the anamorphic basidiomycetous fungus Rhizoctonia solani Kühn affect various agricultural and horticultural important crops including bean, rice, soybean, figs, hortensia, cabbage and lettuce. To gain insights into the genome structure and content, the first draft genome sequence of R. solani AG1-IB isolate 7/3/14 was established. Four complete runs on the Genome Sequencer (GS) FLX platform (Roche Applied Science) yielding approx. a 25-fold coverage of the R. solani genome were accomplished. Assembly of the sequence reads by means of the gsAssembler software version 2.6 applying the heterozygotic mode resulted in numerous contigs and scaffolds and a predicted size of 87.1 Mb for the diploid status of the genome. 'Contig-length vs. read-count' analysis revealed that the assembled contigs can be classified into five different groups. Detailed BLAST-analysis revealed that most contigs of group II feature high-scoring matches to other contigs of the same group suggesting that distinguishable allelic variants exist for many genes. Due to the supposed diploid and heterokaryotic nature of R. solani AG1-IB 7/3/14, this result has been anticipated. However, the heterokaryotic character of the isolate is not really supported by sequencing data obtained for the isolate R. solani AG1-IB 7/3/14. Coverage of group III contigs is twice as high as for group II contigs which can also be explained by the diploid status of the genome and indistinguishable alleles on homologous chromosomes. Assembly of sequence data led to the identification of the rRNA unit (group V contigs) and the mitochondrial (mt) genome (group IV contigs) which is a circular molecule of 162,751 bp in size featuring a GC-content of 36.4%. The R. solani 7/3/14 mt-genome is one of the largest fungal mitochondrial genomes known to date. Its large size essentially is due to the presence of numerous non-conserved hypothetical ORFs and introns. Gene prediction for the R. solani AG1-IB 7/3/14 genome was conducted by the Augustus Gene Prediction Software for Eukaryotes (version 2.6.) applying the parameter set for the fungus Coprinopsis cinerea okayama 7#130. Gene prediction and annotation provided first insights into the R. solani AG1-IB 7/3/14 gene structure and content. In total, 12,422 genes were predicted. The average number of exons per gene is five. Exons have a mean length of 214 bp, whereas introns on average are 66 bp in length. Annotation of the genome revealed that 4169 of 12,422 genes could be assigned to KOG functional categories.

The genome and transcriptome of the pine saprophyte Ophiostoma piceae, and a comparison with the bark beetle-associated pine pathogen Grosmannia clavigera

BACKGROUND

Ophiostoma piceae is a wood-staining fungus that grows in the sapwood of conifer logs and lumber. We sequenced its genome and analyzed its transcriptomes under a range of growth conditions. A comparison with the genome and transcriptomes of the mountain pine beetle-associated pathogen Grosmannia clavigera highlights differences between a pathogen that colonizes and kills living pine trees and a saprophyte that colonizes wood and the inner bark of dead trees.

RESULTS

We assembled a 33 Mbp genome in 45 scaffolds, and predicted approximately 8,884 genes. The genome size and gene content were similar to those of other ascomycetes. Despite having similar ecological niches, O. piceae and G. clavigera showed no large-scale synteny. We identified O. piceae genes involved in the biosynthesis of melanin, which causes wood discoloration and reduces the commercial value of wood products. We also identified genes and pathways involved in growth on simple carbon sources and in sapwood, O. piceae's natural substrate. Like the pathogen, the saprophyte is able to tolerate terpenes, which are a major class of pine tree defense compounds; unlike the pathogen, it cannot utilize monoterpenes as a carbon source.

CONCLUSIONS

This work makes available the second annotated genome of a softwood ophiostomatoid fungus, and suggests that O. piceae's tolerance to terpenes may be due in part to these chemicals being removed from the cells by an ABC transporter that is highly induced by terpenes. The data generated will provide the research community with resources for work on host-vector-fungus interactions for wood-inhabiting, beetle-associated saprophytes and pathogens.

Complete mitochondrial genome of the aluminum-tolerant fungus Rhodotorula taiwanensis RS1 and comparative analysis of Basidiomycota mitochondrial genomes

The complete mitochondrial genome of Rhodotorula taiwanensis RS1, an aluminum-tolerant Basidiomycota fungus, was determined and compared with the known mitochondrial genomes of 12 Basidiomycota species. The mitochondrial genome of R. taiwanensis RS1 is a circular DNA molecule of 40,392 bp and encodes the typical 15 mitochondrial proteins, 23 tRNAs, and small and large rRNAs as well as 10 intronic open reading frames. These genes are apparently transcribed in two directions and do not show syntenies in gene order with other investigated Basidiomycota species. The average G+C content (41%) of the mitochondrial genome of R. taiwanensis RS1 is the highest among the Basidiomycota species. Two introns were detected in the sequence of the atp9 gene of R. taiwanensis RS1, but not in that of other Basidiomycota species. Rhodotorula taiwanensis is the first species of the genus Rhodotorula whose full mitochondrial genome has been sequenced; and the data presented here supply valuable information for understanding the evolution of fungal mitochondrial genomes and researching the mechanism of aluminum tolerance in microorganisms.

Comparing COI and ITS as DNA barcode markers for mushrooms and allies (Agaricomycotina)

DNA barcoding is an approach to rapidly identify species using short, standard genetic markers. The mitochondrial cytochrome oxidase I gene (COI) has been proposed as the universal barcode locus, but its utility for barcoding in mushrooms (ca. 20,000 species) has not been established. We succeeded in generating 167 partial COI sequences (~450 bp) representing ~100 morphospecies from ~650 collections of Agaricomycotina using several sets of new primers. Large introns (~1500 bp) at variable locations were detected in ~5% of the sequences we obtained. We suspect that widespread presence of large introns is responsible for our low PCR success (~30%) with this locus. We also sequenced the nuclear internal transcribed spacer rDNA regions (ITS) to compare with COI. Among the small proportion of taxa for which COI could be sequenced, COI and ITS perform similarly as a barcode. However, in a densely sampled set of closely related taxa, COI was less divergent than ITS and failed to distinguish all terminal clades. Given our results and the wealth of ITS data already available in public databases, we recommend that COI be abandoned in favor of ITS as the primary DNA barcode locus in mushrooms.

A biologist's guide to de novo genome assembly using next-generation sequence data: A test with fungal genomes

We offer a guide to de novo genome assembly using sequence data generated by the Illumina platform for biologists working with fungi or other organisms whose genomes are less than 100Mb in size. The guide requires no familiarity with sequencing assembly technology or associated computer programs. It defines commonly used terms in genome sequencing and assembly; provides examples of assembling short-read genome sequence data for four strains of the fungus Grosmannia clavigera using four assembly programs; gives examples of protocols and software; and presents a commented flowchart that extends from DNA preparation for submission to a sequencing center, through to processing and assembly of the raw sequence reads using freely available operating systems and software.

The Agaricus bisporus cox1 gene: the longest mitochondrial gene and the largest reservoir of mitochondrial group i introns

In eukaryotes, introns are located in nuclear and organelle genes from several kingdoms. Large introns (up to 5 kbp) are frequent in mitochondrial genomes of plant and fungi but scarce in Metazoa, even if these organisms are grouped with fungi among the Opisthokonts. Mitochondrial introns are classified in two groups (I and II) according to their RNA secondary structure involved in the intron self-splicing mechanism. Most of these mitochondrial group I introns carry a “Homing Endonuclease Gene” (heg) encoding a DNA endonuclease acting in transfer and site-specific integration (“homing”) and allowing intron spreading and gain after lateral transfer even between species from different kingdoms. Opposed to this gain mechanism, is another which implies that introns, which would have been abundant in the ancestral genes, would mainly evolve by loss. The importance of both mechanisms (loss and gain) is matter of debate. Here we report the sequence of the cox1 gene of the button mushroom Agaricus bisporus, the most widely cultivated mushroom in the world. This gene is both the longest mitochondrial gene (29,902 nt) and the largest group I intron reservoir reported to date with 18 group I and 1 group II. An exhaustive analysis of the group I introns available in cox1 genes shows that they are mobile genetic elements whose numerous events of loss and gain by lateral transfer combine to explain their wide and patchy distribution extending over several kingdoms. An overview of intron distribution, together with the high frequency of eroded heg, suggests that they are evolving towards loss. In this landscape of eroded and lost intron sequences, the A. bisporus cox1 gene exhibits a peculiar dynamics of intron keeping and catching, leading to the largest collection of mitochondrial group I introns reported to date in a Eukaryote.