Comparative Analysis of the Complete Plastid Genome Sequence of the Red Alga Gracilaria tenuistipitata var. liui Provides Insights into the Evolution of Rhodoplasts and Their Relationship to Other Plastids

We sequenced to completion the circular plastid genome of the red alga Gracilaria tenuistipitata var. liui. This is the first plastid genome sequence from the subclass Florideophycidae (Rhodophyta). The genome is composed of 183,883 bp and contains 238 predicted genes, including a single copy of the ribosomal RNA operon. Comparisons with the plastid genome of Porphyra pupurea reveal strong conservation of gene content and order, but we found major genomic rearrangements and the presence of coding regions that are specific to Gracilaria. Phylogenetic analysis of a data set of 41 concatenated proteins from 23 plastid and two cyanobacterial genomes support red algal plastid monophyly and a specific evolutionary relationship between the Florideophycidae and the Bangiales. Gracilaria maintains a surprisingly ancient gene content in its plastid genome and, together with other Rhodophyta, contains the most complete repertoire of plastid genes known in photosynthetic eukaryotes.Supplementary material () is available for this article.[Reviewing Editor: Dr. W. Ford Doolittle]     

Evolution of Red Algal Plastid Genomes: Ancient Architectures,Introns, Horizontal Gene Transfer,and Taxonomic Utility of Plastid Markers

Red algae have the most gene-rich plastid genomes known, but despite their evolutionary importance these genomes remain poorly sampled. Here we characterize three complete and one partial plastid genome from a diverse range of florideophytes. By unifying annotations across all available red algal plastid genomes we show they all share a highly compact and slowly-evolving architecture and uniquely rich gene complements. Both chromosome structure and gene content have changed very little during red algal diversification, and suggest that plastid-to nucleus gene transfers have been rare. Despite their ancient character, however, the red algal plastids also contain several unprecedented features, including a group II intron in a tRNA-Met gene that encodes the first example of red algal plastid intron maturase – a feature uniquely shared among florideophytes. We also identify a rare case of a horizontally-acquired proteobacterial operon, and propose this operon may have been recruited for plastid function and potentially replaced a nucleus-encoded plastid-targeted paralogue. Plastid genome phylogenies yield a fully resolved tree and suggest that plastid DNA is a useful tool for resolving red algal relationships. Lastly, we estimate the evolutionary rates among more than 200 plastid genes, and assess their usefulness for species and subspecies taxonomy by comparison to well-established barcoding markers such as cox1 and rbcL. Overall, these data demonstrates that red algal plastid genomes are easily obtainable using high-throughput sequencing of total genomic DNA, interesting from evolutionary perspectives, and promising in resolving red algal relationships at evolutionarily-deep and species/subspecies levels.     

Hypothesis: Gene‐rich plastid genomes in red algae may be an outcome of nuclear genome reduction

Red algae (Rhodophyta) putatively diverged from the eukaryote tree of life >1.2 billion years ago and are the source of plastids in the ecologically important diatoms, haptophytes, and dinoflagellates. In general, red algae contain the largest plastid gene inventory among all such organelles derived from primary, secondary, or additional rounds of endosymbiosis. In contrast, their nuclear gene inventory is reduced when compared to their putative sister lineage, the Viridiplantae, and other photosynthetic lineages. The latter is thought to have resulted from a phase of genome reduction that occurred in the stem lineage of Rhodophyta. A recent comparative analysis of a taxonomically broad collection of red algal and Viridiplantae plastid genomes demonstrates that the red algal ancestor encoded ~1.5× more plastid genes than Viridiplantae. This difference is primarily explained by more extensive endosymbiotic gene transfer (EGT) in the stem lineage of Viridiplantae, when compared to red algae. We postulate that limited EGT in Rhodophytes resulted from the countervailing force of ancient, and likely recurrent, nuclear genome reduction. In other words, the propensity for nuclear gene loss led to the retention of red algal plastid genes that would otherwise have undergone intracellular gene transfer to the nucleus. This hypothesis recognizes the primacy of nuclear genome evolution over that of plastids, which have no inherent control of their gene inventory and can change dramatically (e.g., secondarily non‐photosynthetic eukaryotes, dinoflagellates) in response to selection acting on the host lineage.     

Reappraising plastid markers of the red algae for phylogenetic community ecology in the genomic era

Selection of appropriate genetic markers to quantify phylogenetic diversity is crucial for community ecology studies. Yet, systematic evaluation of marker genes for this purpose is scarcely done. Recently, the combined effort of phycologists has produced a rich plastid genome resource with taxonomic representation spanning all of the major lineages of the red algae (Rhodophyta). In this proof‐of‐concept study, we leveraged this resource by developing and applying a phylogenomic strategy to seek candidate plastid markers suitable for phylogenetic community analysis. We ranked the core genes of 107 published plastid genomes based on various sequence‐derived properties and their tree distance to plastid genome phylogenies. The resulting ranking revealed that the most widely used marker, rbcL, is not necessarily the optimal marker, while other promising markers might have been overlooked. We designed and tested PCR primers for several candidate marker genes, and successfully amplified one of them, rpoC1, in a taxonomically broad set of red algal specimens. We suggest that our general marker identification methodology and the rpoC1 primers will be useful to the phycological community for investigating the biodiversity and community ecology of the red algae.     

Nucleomorph and plastid genome sequences of the chlorarachniophyte Lotharella oceanica: convergent reductive evolution and frequent recombination in nucleomorph-bearing algae

Background

Nucleomorphs are residual nuclei derived from eukaryotic endosymbionts in chlorarachniophyte and cryptophyte algae. The endosymbionts that gave rise to nucleomorphs and plastids in these two algal groups were green and red algae, respectively. Despite their independent origin, the chlorarachniophyte and cryptophyte nucleomorph genomes share similar genomic features such as extreme size reduction and a three-chromosome architecture. This suggests that similar reductive evolutionary forces have acted to shape the nucleomorph genomes in the two groups. Thus far, however, only a single chlorarachniophyte nucleomorph and plastid genome has been sequenced, making broad evolutionary inferences within the chlorarachniophytes and between chlorarachniophytes and cryptophytes difficult. We have sequenced the nucleomorph and plastid genomes of the chlorarachniophyte Lotharella oceanica in order to gain insight into nucleomorph and plastid genome diversity and evolution.

Results

The L. oceanica nucleomorph genome was found to consist of three linear chromosomes totaling ~610 kilobase pairs (kbp), much larger than the 373 kbp nucleomorph genome of the model chlorarachniophyte Bigelowiella natans. The L. oceanica plastid genome is 71 kbp in size, similar to that of B. natans. Unexpectedly long (~35 kbp) sub-telomeric repeat regions were identified in the L. oceanica nucleomorph genome; internal multi-copy regions were also detected. Gene content analyses revealed that nucleomorph house-keeping genes and spliceosomal intron positions are well conserved between the L. oceanica and B. natans nucleomorph genomes. More broadly, gene retention patterns were found to be similar between nucleomorph genomes in chlorarachniophytes and cryptophytes. Chlorarachniophyte plastid genomes showed near identical protein coding gene complements as well as a high level of synteny.

Conclusions

We have provided insight into the process of nucleomorph genome evolution by elucidating the fine-scale dynamics of sub-telomeric repeat regions. Homologous recombination at the chromosome ends appears to be frequent, serving to expand and contract nucleomorph genome size. The main factor influencing nucleomorph genome size variation between different chlorarachniophyte species appears to be expansion-contraction of these telomere-associated repeats rather than changes in the number of unique protein coding genes. The dynamic nature of chlorarachniophyte nucleomorph genomes lies in stark contrast to their plastid genomes, which appear to be highly stable in terms of gene content and synteny.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-374) contains supplementary material, which is available to authorized users.     

Molecular phylogenetics supports a clade of red algal parasites retaining native plastids: taxonomy and terminology revised

Parasitism is a life strategy that has repeatedly evolved within the Florideophyceae. Historically, the terms adelphoparasite and alloparasite have been used to distinguish parasites based on the relative phylogenetic relationship of host and parasite. However, analyses using molecular phylogenetics indicate that nearly all red algal parasites infect within their taxonomic family, and a range of relationships exist between host and parasite. To date, all investigated adelphoparasites have lost their plastid, and instead, incorporate a host‐derived plastid when packaging spores. In contrast, a highly reduced plastid lacking photosynthesis genes was sequenced from the alloparasite Choreocolax polysiphoniae. Here we present the complete Harveyella mirabilis plastid genome, which has also lost genes involved in photosynthesis, and a partial plastid genome from Leachiella pacifica. The H. mirabilis plastid shares more synteny with free‐living red algal plastids than that of C. polysiphoniae. Phylogenetic analysis demonstrates that C. polysiphoniae, H. mirabilis, and L. pacifica form a robustly supported clade of parasites, which retain their own plastid genomes, within the Rhodomelaceae. We therefore transfer all three genera from the exclusively parasitic family, Choreocolacaceae, to the Rhodomelaceae. Additionally, we recommend applying the terms archaeplastic parasites (formerly alloparasites), and neoplastic parasites (formerly adelphoparasites) to distinguish red algal parasites using a biological framework rather than taxonomic affiliation with their hosts.     

Analysis of the complete plastid genome of the unicellular red alga Porphyridium purpureum

We determined the complete nucleotide sequence of the plastid genome of the unicellular marine red alga Porphyridium purpureum strain NIES 2140, belonging to the unsequenced class Porphyridiophyceae. The genome is a circular DNA composed of 217,694 bp with the GC content of 30.3 %. Twenty-nine of the 224 protein-coding genes contain one or multiple intron(s). A group I intron was found in the rpl28 gene, whereas the other introns were group II introns. The P. purpureum plastid genome has one non-coding RNA (ncRNA) gene, 29 tRNA genes and two nonidentical ribosomal RNA operons. One rRNA operon has a tRNA^Ala(UGC) gene between the rrs and the rrl genes, whereas another has a tRNA^Ile(GAU) gene. Phylogenetic analyses suggest that the plastids of Heterokontophyta, Cryptophyta and Haptophyta originated from the subphylum Rhodophytina. The order of the genes in the ribosomal protein cluster of the P. purpureum plastid genome differs from that of other Rhodophyta and Chromalveolata. These results suggest that a large-scale rearrangement occurred in the plastid genome of P. purpureum after its separation from other Rhodophyta.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

The Mitochondrial Genome of an Aquatic Plant,Spirodela polyrhiza

Background

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

Methods

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

Conclusions

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.     

The Complete Chloroplast and Mitochondrial Genomes of the Green Macroalga Ulva sp. UNA00071828 (Ulvophyceae,Chlorophyta)

Sequencing mitochondrial and chloroplast genomes has become an integral part in understanding the genomic machinery and the phylogenetic histories of green algae. Previously, only three chloroplast genomes (Oltmannsiellopsis viridis, Pseudendoclonium akinetum, and Bryopsis hypnoides) and two mitochondrial genomes (O. viridis and P. akinetum) from the class Ulvophyceae have been published. Here, we present the first chloroplast and mitochondrial genomes from the ecologically and economically important marine, green algal genus Ulva. The chloroplast genome of Ulva sp. was 99,983 bp in a circular-mapping molecule that lacked inverted repeats, and thus far, was the smallest ulvophycean plastid genome. This cpDNA was a highly compact, AT-rich genome that contained a total of 102 identified genes (71 protein-coding genes, 28 tRNA genes, and three ribosomal RNA genes). Additionally, five introns were annotated in four genes: atpA (1), petB (1), psbB (2), and rrl (1). The circular-mapping mitochondrial genome of Ulva sp. was 73,493 bp and follows the expanded pattern also seen in other ulvophyceans and trebouxiophyceans. The Ulva sp. mtDNA contained 29 protein-coding genes, 25 tRNA genes, and two rRNA genes for a total of 56 identifiable genes. Ten introns were annotated in this mtDNA: cox1 (4), atp1 (1), nad3 (1), nad5 (1), and rrs (3). Double-cut-and-join (DCJ) values showed that organellar genomes across Chlorophyta are highly rearranged, in contrast to the highly conserved organellar genomes of the red algae (Rhodophyta). A phylogenomic investigation of 51 plastid protein-coding genes showed that Ulvophyceae is not monophyletic, and also placed Oltmannsiellopsis (Oltmannsiellopsidales) and Tetraselmis (Chlorodendrophyceae) closely to Ulva (Ulvales) and Pseudendoclonium (Ulothrichales).     

The Complete Plastid Genome Sequence of the Haptophyte Emiliania huxleyi: a Comparison to Other Plastid Genomes

The complete nucleotide sequence of the plastid genome of thehaptophyte Emiliania huxleyi has been determined. E. huxleyiis the most abundant coccolithophorid and has a key role inthe carbon cycle. It is also implicated in the production ofdimethylsulphide (DMS), which is involved in cloud nucleationand may affect the global climate. Here, we report the plastidgenome sequence of this ecologically and economically importantspecies and compare its gene content and arrangement to otherknown plastid genomes. The genome is circular and consists of105,309 bp with an inverted repeat of 4,841 bp. In terms ofboth genome size and gene content E. huxleyi cpDNA is substantiallysmaller than any other from the red plastid lineage. The geneticinformation is densely packed, with 86.8% of the genome specifying110 identified protein-coding genes, 9 open reading frames,28 different tRNAs, and 3 rRNAs. A detailed comparison to otherplastid genomes, based on gene content, gene function, and genecluster analysis is discussed. These analyses suggest a closerelationship of the E. huxleyi cpDNA to the chlorophyll c-containingplastids from heterokonts and cryptophytes, and they supportthe origin of the chromophyte plastids from the red algal lineage.     

Complete Plastid Genome Sequence of the Brown Alga Undaria pinnatifida

In this study, we fully sequenced the circular plastid genome of a brown alga, Undaria pinnatifida. The genome is 130,383 base pairs (bp) in size; it contains a large single-copy (LSC, 76,598 bp) and a small single-copy region (SSC, 42,977 bp), separated by two inverted repeats (IRa and IRb: 5,404 bp). The genome contains 139 protein-coding, 28 tRNA, and 6 rRNA genes; none of these genes contains introns. Organization and gene contents of the U. pinnatifida plastid genome were similar to those of Saccharina japonica. There is a co-linear relationship between the plastid genome of U. pinnatifida and that of three previously sequenced large brown algal species. Phylogenetic analyses of 43 taxa based on 23 plastid protein-coding genes grouped all plastids into a red or green lineage. In the large brown algae branch, U. pinnatifida and S. japonica formed a sister clade with much closer relationship to Ectocarpus siliculosus than to Fucus vesiculosus. For the first time, the start codon ATT was identified in the plastid genome of large brown algae, in the atpA gene of U. pinnatifida. In addition, we found a gene-length change induced by a 3-bp repetitive DNA in ycf35 and ilvB genes of the U. pinnatifida plastid genome.     

Complete Plastid Genome Sequencing of Four Tilia Species (Malvaceae): A Comparative Analysis and Phylogenetic Implications

Tilia is an ecologically and economically important genus in the family Malvaceae. However, there is no complete plastid genome of Tilia sequenced to date, and the taxonomy of Tilia is difficult owing to frequent hybridization and polyploidization. A well-supported interspecific relationships of this genus is not available due to limited informative sites from the commonly used molecular markers. We report here the complete plastid genome sequences of four Tilia species determined by the Illumina technology. The Tilia plastid genome is 162,653 bp to 162,796 bp in length, encoding 113 unique genes and a total number of 130 genes. The gene order and organization of the Tilia plastid genome exhibits the general structure of angiosperms and is very similar to other published plastid genomes of Malvaceae. As other long-lived tree genera, the sequence divergence among the four Tilia plastid genomes is very low. And we analyzed the nucleotide substitution patterns and the evolution of insertions and deletions in the Tilia plastid genomes. Finally, we build a phylogeny of the four sampled Tilia species with high supports using plastid phylogenomics, suggesting that it is an efficient way to resolve the phylogenetic relationships of this genus.     

Gene-rich plastid genomes of two parasitic red algal species,Laurencia australis and L. verruciformis (Rhodomelaceae,Ceramiales), and a taxonomic revision of Janczewskia

Parasitic red algae are an interesting system for investigating the genetic changes that occur in parasites. These parasites have evolved independently multiple times within the red algae. The functional loss of plastid genomes can be investigated in these multiple independent examples, and fine-scale patterns may be discerned. The only plastid genomes from red algal parasites known so far are highly reduced and missing almost all photosynthetic genes. Our study assembled and annotated plastid genomes from the parasites Janczewskia tasmanica and its two Laurencia host species (Laurencia elata and one unidentified Laurencia sp. A25) from Australia and Janczewskia verruciformis, its host species (Laurencia catarinensis), and the closest known free-living relative (Laurencia obtusa) from the Canary Islands (Spain). For the first time we show parasitic red algal plastid genomes that are similar in size and gene content to free-living host species without any gene loss or genome reduction. The only exception was two pseudogenes (moeB and ycf46) found in the plastid genome of both isolates of J. tasmanica, indicating potential for future loss of these genes. Further comparative analyses with the three highly reduced plastid genomes showed possible gene loss patterns, in which photosynthetic gene categories were lost followed by other gene categories. Phylogenetic analyses did not confirm monophyly of Janczewskia, and the genus was subsumed into Laurencia. Further investigations will determine if any convergent small-scale patterns of gene loss exist in parasitic red algae and how these are applicable to other parasitic systems.     

Phylogenomics and multigene phylogenies decipher two new cryptic marine algae from California,Gelidium gabrielsonii and G. kathyanniae (Gelidiales,Rhodophyta)

Molecular surveys are leading to the discovery of many new cryptic species of marine algae. This is particularly true for red algal intertidal species, which exhibit a high degree of morphological convergence. DNA sequencing of recent collections of Gelidium along the coast of California, USA, identified two morphologically similar entities that differed in DNA sequence from existing species. To characterize the two new species of Gelidium and to determine their evolutionary relationships to other known taxa, phylogenomic, multigene analyses, and morphological observations were performed. Three complete mitogenomes and five plastid genomes were deciphered, including those from the new species candidates and the type materials of two closely related congeners. The mitogenomes contained 45 genes and had similar lengths (24,963–24,964 bp). The plastid genomes contained 232 genes and were roughly similar in size (175,499–177,099 bp). The organellar genomes showed a high level of gene synteny. The two Gelidium species are diminutive, turf‐forming, and superficially resemble several long established species from the Pacific Ocean. The phylogenomic analysis, multigene phylogeny, and morphological evidence confirms the recognition and naming of two new species, describe herein as G. gabrielsonii and G. kathyanniae. On the basis of the monophyly of G. coulteri, G. gabrielsonii, G. galapagense, and G. kathyanniae, we suggest that this lineage likely evolved in California. Organellar genomes provide a powerful tool for discovering cryptic intertidal species and they continue to improve our understanding of the evolutionary biology of red algae and the systematics of the Gelidiales.     

When Less is More: Red Algae as Models for Studying Gene Loss and Genome Evolution in Eukaryotes

Genome evolution is usually viewed through the lens of growth in size and complexity over time, exemplified by plants and animals. In contrast, genome reduction is associated with a narrowing of ecological potential, such as in parasites and endosymbionts. But, can nuclear genome reduction also occur in, and potentially underpin a major radiation of free-living eukaryotes? An intriguing example of this phenomenon is provided by the red algae (Rhodophyta) that have lost many conserved pathways such as for flagellar motility, macroautophagy regulation, and phytochrome based light sensing. This anciently diverged, species-rich, and ecologically important algal lineage has undergone at least two rounds of large-scale genome reduction during its >1 billion-year evolutionary history. Here, using recent analyses of genome data, we review knowledge about the evolutionary trajectory of red algal nuclear and organelle gene inventories and plastid encoded autocatalytic introns. We compare and contrast Rhodophyta genome evolution to Viridiplantae (green algae and plants), both of which are members of the Archaeplastida, and highlight their divergent paths. We also discuss evidence for the speculative hypothesis that reduction in red algal plastid genome size through endosymbiotic gene transfer is counteracted by ongoing selection for compact nuclear genomes in red algae. Finally, we describe how the spliceosomal intron splicing apparatus provides an example of “evolution in action” in Rhodophyta and how the overall constraints on genome size in this lineage has left significant imprints on this key step in RNA maturation. Our review reveals the red algae to be an exciting, yet under-studied model that offers numerous novel insights as well as many unanswered questions that remain to be explored using modern genomic, genetic, and biochemical methods. The fact that a speciose lineage of free-living eukaryotes has spread throughout many aquatic habitats after having lost about 25% of its primordial gene inventory challenges us to elucidate the mechanisms underlying this remarkable feat.     

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.