Landscape genomics: natural selection drives the evolution of mitogenome in penguins

Complete Mitochondrial Genome and Phylogenetic Analysis of Tarsiger indicus (Aves: Passeriformes: Muscicapidae)

Tarsiger indicus (Vieillot, 1817), the White-browed Bush Robin, is a small passerine bird widely distributed in Asian countries. Here, we successfully sequenced its mitogenome using the Illumina Novaseq 6000 platform (Illumina, San Diego, CA, USA) for PE 2 × 150 bp sequencing. Combined with other published mitogenomes, we conducted the first comprehensive comparative mitogenome analysis of Muscicapidae birds and reconstructed the phylogenetic relationships between Muscicapidae and related groups. The T. indicus mitogenome was 16,723 bp in size, and it possessed the typical avian mitogenome structure and organization. Most PCGs of T. indicus were initiated strictly with the typical start codon ATG, while COX1 and ND2 were started with GTG. RSCU statistics showed that CUA, CGA, and GCC were relatively high frequency in the T. indicus mitogenome. T. cyanurus and T. indicus shared very similar mitogenomic features. All 13 PCGs of Muscicapidae mitogenomes had experienced purifying selection. Specifically, ATP8 had the highest rate of evolution (0.13296), whereas COX1 had the lowest (0.01373). The monophylies of Muscicapidae, Turdidae, and Paradoxornithidae were strongly supported. The clade of ((Muscicapidae + Turdidae) + Sturnidae) in Passeriformes was supported by both Bayesian Inference and Maximum likelihood analyses. The latest taxonomic status of many passerine birds with complex taxonomic histories were also supported. For example, Monticola gularis, T. indicus, and T. cyanurus were allocated to Turdidae in other literature; our phylogenetic topologies clearly supported their membership in Muscicapidae; Paradoxornis heudei, Suthora webbiana, S. nipalensis, and S. fulvifrons were formerly classified into Muscicapidae; we supported their membership in Paradoxornithidae; Culicicapa ceylonensis was originally classified as a member of Muscicapidae; our results are consistent with a position in Stenostiridae. Our study enriches the genetic data of T. indicus and provides new insights into the molecular phylogeny and evolution of passerine birds.

The mitochondrial genome of the mountain wooly tapir, Tapirus pinchaque and a formal test of the effect of altitude on the adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates

BACKGROUND

The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments.

RESULTS

The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4.

CONCLUSION

The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

Natural selection drives the evolution of mitogenomes in Acrossocheilus

The mitochondrial genome plays a crucial role in the balance of energy and heat production in organisms and, thus, may be under natural selection due to its potential role in adaptive divergence and speciation. Here, we studied natural selection on the mitogenome of Acrossocheilus (Cypriniformes: Cyprinidae), a genus of fish that inhabits a broad latitudinal distribution ranging from the tropics and subtropics through temperate regions. Specifically, we used 25 published mitogenome sequences of Acrossocheilus species to investigate phylogenetic relationships in this genus and detected signals of positive selection on 13 protein-coding, mitochondrial genes. We found that relaxed purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, we also found evidence of diversifying selection on some codons, indicating episodes of positive selection. Additionally, we analyzed the mitogenomic data within an environmental modeling framework and found that the Ka/Ks ratio of ATP6 may correlated with a mean diurnal temperature range (p = 0.0449), while the Ka/Ks ratio of COX2 may correlated with precipitation during the driest month (p = 0.00761). These results suggest that the mitogenomes of Acrossocheilus species may be involved in evolutionary adaptations to different habitats. Based on this, we believe that our study provides a new insight into the role of the mitochondrial genome of Acrossocheilus species in adaptation to different environments. During our study, we also discovered several cases of paraphyly and polyphyly among accessions of species and their putative synonyms. Thus, our study suggests that a careful reassessment of the taxonomy of Acrossocheilus is using high-quality molecular data merited.

Complete Mitogenomes of Polypedates Tree Frogs Unveil Gene Rearrangement and Concerted Evolution within Rhacophoridae

Simple Summary

Duplicated control regions have been reported several times in the tree frog family Rhacophoridae, and previous studies have mostly relied on sequence analysis to reconstruct their evolution. This is the first study to employ a phylogenetic method to demonstrate the existence of concerted and parallel evolution succinctly and intuitively in the duplicated control regions of the family Rhacophoridae. Phylogenetic relationships were also used to illustrate the parallel evolution of ATP8 loss of function in the genus Polypedates. In general, this study elucidated the evolutionary patterns and pathways of mitochondrial gene rearrangement of the family Rhacophoridae from a phylogenetic perspective, which aids in understanding the evolutionary history of this fascinating tree frog taxon from a molecular evolution standpoint.

Abstract

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy “LTPF” tRNA gene cluster compared to the novel “TLPF” order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

The complete mitochondrial genome of the edible mushroom Pleurotus giganteus (Agaricales, Pleurotus) and insights into its phylogeny

Pleurotus giganteus (Berk.) Karunarathna & K.D. Hyde 2011 is one of the largest edible mushrooms integrating medicinal value and edible value. The complete mitochondrial genome of the edible fungus P. giganteus was published in this paper. It was determined using Pacbio and Illumina sequencing. The circular molecule is 102,950 bp in length, consisting of 30 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 24 transfer RNA (tRNA) genes. The base composition of the whole mitogenome is A (37.3%), T (37.7%), G (12.2%), and C (12.8%). The phylogenetic tree shows P. giganteus was the basal taxon in Pleurotus and closely related to Pleurotus citrinopileatus Singer 1990.

Positive selection over the mitochondrial genome and its role in the diversification of gentoo penguins in response to adaptation in isolation

Although mitochondrial DNA has been widely used in phylogeography, evidence has emerged that factors such as climate, food availability, and environmental pressures that produce high levels of stress can exert a strong influence on mitochondrial genomes, to the point of promoting the persistence of certain genotypes in order to compensate for the metabolic requirements of the local environment. As recently discovered, the gentoo penguins (Pygoscelis papua) comprise four highly divergent lineages across their distribution spanning the Antarctic and sub-Antarctic regions. Gentoo penguins therefore represent a suitable animal model to study adaptive processes across divergent environments. Based on 62 mitogenomes that we obtained from nine locations spanning all four gentoo penguin lineages, we demonstrated lineage-specific nucleotide substitutions for various genes, but only lineage-specific amino acid replacements for the ND1 and ND5 protein-coding genes. Purifying selection (dN/dS < 1) is the main driving force in the protein-coding genes that shape the diversity of mitogenomes in gentoo penguins. Positive selection (dN/dS > 1) was mostly present in codons of the Complex I (NADH genes), supported by two different codon-based methods at the ND1 and ND4 in the most divergent lineages, the eastern gentoo penguin from Crozet and Marion Islands and the southern gentoo penguin from Antarctica respectively. Additionally, ND5 and ATP6 were under selection in the branches of the phylogeny involving all gentoo penguins except the eastern lineage. Our study suggests that local adaptation of gentoo penguins has emerged as a response to environmental variability promoting the fixation of mitochondrial haplotypes in a non-random manner. Mitogenome adaptation is thus likely to have been associated with gentoo penguin diversification across the Southern Ocean and to have promoted their survival in extreme environments such as Antarctica. Such selective processes on the mitochondrial genome may also be responsible for the discordance detected between nuclear- and mitochondrial-based phylogenies of gentoo penguin lineages.

Phylogeny and evolution of Lasiopodomys in subfamily Arvivolinae based on mitochondrial genomics

The species of Lasiopodomys Lataste 1887 with their related genera remains undetermined owing to inconsistent morphological characteristics and molecular phylogeny. To investigate the phylogenetic relationship and speciation among species of the genus Lasiopodomys, we sequenced and annotated the whole mitochondrial genomes of three individual species, namely Lasiopodomys brandtii Radde 1861, L. mandarinus Milne-Edwards 1871, and Neodon (Lasiopodomys) fuscus Büchner 1889. The nucleotide sequences of the circular mitogenomes were identical for each individual species of L. brandtii, L. mandarinus, and N. fuscus. Each species contained 13 protein-coding genes (PCGs), 22 transfer RNAs, and 2 ribosomal RNAs, with mitochondrial genome lengths of 16,557 bp, 16,562 bp, and 16,324 bp, respectively. The mitogenomes and PCGs showed positive AT skew and negative GC skew. Mitogenomic phylogenetic analyses suggested that L. brandtii, L. mandarinus, and L. gregalis Pallas 1779 belong to the genus Lasiopodomys, whereas N. fuscus belongs to the genus Neodon grouped with N. irene. Lasiopodomys showed the closest relationship with Microtus fortis Büchner 1889 and M. kikuchii Kuroda 1920, which are considered as the paraphyletic species of genera Microtus. T_MRCA and niche model analysis revealed that Lasiopodomys may have first appeared during the early Pleistocene epoch. Further, L. gregalis separated from others over 1.53 million years ago (Ma) and then diverged into L. brandtii and L. mandarinus 0.76 Ma. The relative contribution of climatic fluctuations to speciation and selection in this group requires further research.

Revealing the biodiversity of Chilean birds through the COI barcode approach

The mitochondrial cytochrome c oxidase subunit I (COI) gene is an effective molecular tool for the estimation of genetic variation and the identification of bird species. This molecular marker is used to differentiate among Chilean bird species by analyzing barcodes for 76 species (197 individuals), comprising 28 species with no previous barcode data and 48 species with sequences retrieved from the BOLD and GenBank databases. The DNA barcodes correctly identified 94.7% of the species analyzed (72 of 76 species). Mean intraspecific K2P distance was 0.3% (range 0-8.7%). Within the intraspecific divergence range, three species, Phrygilus gayi, Sephanoides sephanoides and Curaeus curaeus, showed relatively high intraspecific divergence (1.5-8.7%), possibly due to the presence of a species complex or geographic isolation of sub-populations. Mean interspecific K2P distance was 24.7% (range 1.3-43.5%). Consequently, the intraspecific K2P distance showed limited overlap with interspecific K2P distance. The mean intraspecific divergence in our study was similar to that found in temperate regions of South America (0.24%). However, it was approximately one order of magnitude lower than values reported for bird species in tropical regions of northern South America (1.8-2.13%). This result suggests that bird species from Chile show low levels of genetic structure and divergence. The small overlap between intra- and inter-specific distances implies that COI barcodes could be used as an effective tool to identify nearly all the Chilean bird species analyzed.

Mitogenomes of Two Phallus Mushroom Species Reveal Gene Rearrangement, Intron Dynamics, and Basidiomycete Phylogeny

Phallus indusiatus and Phallus echinovolvatus are edible bamboo mushrooms with pharmacological properties. We sequenced, assembled, annotated, and compared the mitogenomes of these species. Both mitogenomes were composed of circular DNA molecules, with sizes of 89,139 and 50,098 bp, respectively. Introns were the most important factor in mitogenome size variation within the genus Phallus. Phallus indusiatus, P. echinovolvatus, and Turbinellus floccosus in the subclass Phallomycetidae have conservative gene arrangements. Large-scale gene rearrangements were observed in species representing 42 different genera of Basidiomycetes. A variety of intron position classes were found in the 44 Basidiomycete species analyzed. A novel group II intron from the P. indusiatus mitogenome was compared with other fungus species containing the same intron, and we demonstrated that the insertion sites of the intron had a base preference. Phylogenetic analyses based on combined gene datasets yielded well-supported Bayesian posterior probability (BPP = 1) topologies. This indicated that mitochondrial genes are reliable molecular markers for analyzing the phylogenetic relationships of the Basidiomycetes. This is the first study of the mitogenome of the genus Phallus, and it increases our understanding of the population genetics and evolution of bamboo mushrooms and related species.

The role of selection in the evolution of marine turtles mitogenomes

Sea turtles are the only extant chelonian representatives that inhabit the marine environment. One key to successful colonization of this habitat is the adaptation to different energetic demands. Such energetic requirement is intrinsically related to the mitochondrial ability to generate energy through oxidative phosphorylation (OXPHOS) process. Here, we estimated Testudines phylogenetic relationships from 90 complete chelonian mitochondrial genomes and tested the adaptive evolution of 13 mitochondrial protein-coding genes of sea turtles to determine how natural selection shaped mitochondrial genes of the Chelonioidea clade. Complete mitogenomes showed strong support and resolution, differing at the position of the Chelonioidea clade in comparison to the turtle phylogeny based on nuclear genomic data. Codon models retrieved a relatively increased dN/dS (ω) on three OXPHOS genes for sea turtle lineages. Also, we found evidence of positive selection on at least three codon positions, encoded by NADH dehydrogenase genes (ND4 and ND5). The accelerated evolutionary rates found for sea turtles on COX2, ND1 and CYTB and the molecular footprints of positive selection found on ND4 and ND5 genes may be related to mitochondrial molecular adaptation to stress likely resulted from a more active lifestyle in sea turtles. Our study provides insight into the adaptive evolution of the mtDNA genome in sea turtles and its implications for the molecular mechanism of oxidative phosphorylation.

Complete mitochondrial genome of Iniistius trivittatus and unique variation in two observed inserts between rRNA and tRNA genes in wrasses

Background

The family Labridae made up of 519 species in the world. The functional evolution of the feeding-related jaws leaded to differentiation of species, and the pharyngeal jaw apparatus evolved independently, but evolutionary mechanism still remain unaddressed in wrasses. Mitogenomes data can be used to infer genetic diversification and investigate evolutionary history of wrasses, whereas only eight complete mitogenomes in this family have been sequenced to date. Here, we sequenced the complete mitogenomes of Iniistius trivittatus to investigate genetic differentiation among wrasse species.

Results

We sequenced the complete mitogenomes of I. trivittatus using a novel PCR strategy. The I. trivittatus mitogenomes is 16,820 bp in length and includes 13 protein -coding genes, 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region. Compared to eight known mitochondrial genome, 2 additional noncoding regions (lengths of 121 and 107 bp), or so-called inserts, are found in the intergenic regions 12S rRNA - tRNA^Val - 16S rRNA. The presumed origin of the two rare inserts is from tRNA- related retrotransposons. Compared with cytochrome b gene, the two insert sequences are highly conserved at the intraspecies level, but they showed significant variation and low similarity (< 70%) at the interspecies level. The insert events were only observed in I. trivittatus by checking the phylogenetic trees based on the complete mitogenomes of Labrida species. This finding provides evidence that in the mitogenomes, retrotransposon inserts result in intraspecific homoplasmy and interspecific heteroplasmy by natural selection and adaptation to various environments.

Conclusions

This study found additional mitogenome inserts limited in wrasse species. The rRNA genes with inserts might have experienced a selective pressure for adaptation to feeding modes. Such knowledge can enable a better understanding of molecular mechanism underlying morphological evolution in wrasses.

High-coverage genomes to elucidate the evolution of penguins

BACKGROUND

Penguins (Sphenisciformes) are a remarkable order of flightless wing-propelled diving seabirds distributed widely across the southern hemisphere. They share a volant common ancestor with Procellariiformes close to the Cretaceous-Paleogene boundary (66 million years ago) and subsequently lost the ability to fly but enhanced their diving capabilities. With ∼20 species among 6 genera, penguins range from the tropical Galápagos Islands to the oceanic temperate forests of New Zealand, the rocky coastlines of the sub-Antarctic islands, and the sea ice around Antarctica. To inhabit such diverse and extreme environments, penguins evolved many physiological and morphological adaptations. However, they are also highly sensitive to climate change. Therefore, penguins provide an exciting target system for understanding the evolutionary processes of speciation, adaptation, and demography. Genomic data are an emerging resource for addressing questions about such processes.

RESULTS

Here we present a novel dataset of 19 high-coverage genomes that, together with 2 previously published genomes, encompass all extant penguin species. We also present a well-supported phylogeny to clarify the relationships among penguins. In contrast to recent studies, our results demonstrate that the genus Aptenodytes is basal and sister to all other extant penguin genera, providing intriguing new insights into the adaptation of penguins to Antarctica. As such, our dataset provides a novel resource for understanding the evolutionary history of penguins as a clade, as well as the fine-scale relationships of individual penguin lineages. Against this background, we introduce a major consortium of international scientists dedicated to studying these genomes. Moreover, we highlight emerging issues regarding ensuring legal and respectful indigenous consultation, particularly for genomic data originating from New Zealand Taonga species.

CONCLUSIONS

We believe that our dataset and project will be important for understanding evolution, increasing cultural heritage and guiding the conservation of this iconic southern hemisphere species assemblage.

Molecular genetic studies of Chilean avifauna: an overview about current progress

The biodiversity of Chilean avifauna has been shaped by geological events, geographical factors and extreme climatic regimens. Molecular genetic studies developed to date have begun to provide valuable information to elucidate the effect of these factors on the evolutionary history of Chilean birds. We conducted a search in two bibliographic databases (PubMed and Scopus) over the period between 1998 and 2018 to identify the number and class of articles published in this research field. The quantification of specific research categories was also performed to assess progress in different molecular genetic issues of Chilean birds. Representative studies of the more active research areas: population structure, phylogeography, systematics and taxonomy and genetic diversity, were briefly commented on. Studies on genetic divergence analysis of cryptic species and the underlying genetics at molecular level of adaptive phenotypes in populations, were very scarce. Several types of mitochondrial and nuclear markers have been used, with COI plus D-loop and microsatellite loci, respectively, being the most widely-used markers. However, high-throughput analyses based on complete mitochondrial genomes have also been published. More than one hundred species have been analyzed; the most studied groups being penguins, followed by New World swallows. Although, to date, molecular genetic studies of Chilean birds still require further analysis, the studies available on some species or groups have provided important data to improve our understanding of the origin and evolution of this avifauna.

The complete mitochondrial genome of a wild edible mushroom, Russula griseocarnosa

Russula griseocarnosa is a wild edible ectomycorrhizal mushroom in southern China. In this study, we assembled the complete mitochondrial genome of R. griseocarnosa. Its total length was 60995 bp with a GC content of 21% and contained a total of 52 genes, including 14 standard protein-coding genes, two rRNA genes, 21 tRNA genes and 15 free-standing open reading frames (ORFs). Phylogenetic analysis reflected that the evolutionary processes between R. griseocarnosa and some agaricomycetes.

Mitogenomes Uncover Extinct Penguin Taxa and Reveal Island Formation as a Key Driver of Speciation

The emergence of islands has been linked to spectacular radiations of diverse organisms. Although penguins spend much of their lives at sea, they rely on land for nesting, and a high proportion of extant species are endemic to geologically young islands. Islands may thus have been crucial to the evolutionary diversification of penguins. We test this hypothesis using a fossil-calibrated phylogeny of mitochondrial genomes (mitogenomes) from all extant and recently extinct penguin taxa. Our temporal analysis demonstrates that numerous recent island-endemic penguin taxa diverged following the formation of their islands during the Plio-Pleistocene, including the Galápagos (Galápagos Islands), northern rockhopper (Gough Island), erect-crested (Antipodes Islands), Snares crested (Snares) and royal (Macquarie Island) penguins. Our analysis also reveals two new recently extinct island-endemic penguin taxa from New Zealand's Chatham Islands: Eudyptes warhami sp. nov. and a dwarf subspecies of the yellow-eyed penguin, Megadyptes antipodes richdalei ssp. nov. Eudyptes warhami diverged from the Antipodes Islands erect-crested penguin between 1.1 and 2.5 Ma, shortly after the emergence of the Chatham Islands (∼3 Ma). This new finding of recently evolved taxa on this young archipelago provides further evidence that the radiation of penguins over the last 5 Ma has been linked to island emergence. Mitogenomic analyses of all penguin species, and the discovery of two new extinct penguin taxa, highlight the importance of island formation in the diversification of penguins, as well as the extent to which anthropogenic extinctions have affected island-endemic taxa across the Southern Hemisphere's isolated archipelagos.

Contrasting phylogeographic pattern among Eudyptes penguins around the Southern Ocean

Since at least the middle-Miocene, the Antarctic Polar Front (APF) and the Subtropical Front (STF) appear to have been the main drivers of diversification of marine biota in the Southern Ocean. However, highly migratory marine birds and mammals challenge this paradigm and the importance of oceanographic barriers. Eudyptes penguins range from the Antarctic Peninsula to subantarctic islands and some of the southernmost subtropical islands. Because of recent diversification, the number of species remains uncertain. Here we analyze two mtDNA (HVRI, COI) and two nuclear (ODC, AK1) markers from 13 locations of five putative Eudyptes species: rockhopper (E. filholi, E. chrysocome, and E. moseleyi), macaroni (E. chrysolophus) and royal penguins (E. schlegeli). Our results show a strong phylogeographic structure among rockhopper penguins from South America, subantarctic and subtropical islands supporting the recognition of three separated species of rockhopper penguins. Although genetic divergence was neither observed among macaroni penguins from the Antarctic Peninsula and sub-Antarctic islands nor between macaroni and royal penguins, population genetic analyses revealed population genetic structure in both cases. We suggest that the APF and STF can act as barriers for these species. While the geographic distance between colonies might play a role, their impact/incidence on gene flow may vary between species and colonies.

Characterization and comparative mitogenomic analysis of six newly sequenced mitochondrial genomes from ectomycorrhizal fungi (Russula) and phylogenetic analysis of the Agaricomycetes

In this study, the mitochondrial genomes of six Russula species were sequenced using next generation sequencing. The six mitogenomes were all composed of circular DNA molecules, with lengths ranging from 40,961 bp to 69,423 bp. The length and number of protein coding genes (PCGs), GC content, AT skew, and GC skew varied among the six mitogenomes. The increased number and total size of introns likely contributed to the size expansion of mitogenomes in some Russula species. Gene synteny analysis revealed some gene rearrangements among the six mitochondrial genomes. The nad4L gene had the lowest K2P genetic distance of the 15 core PCGs among the six Russula species, indicating that this gene was highly conserved. The Ka/Ks values for all 15 core PCGs were <1, suggesting that they were all subject to purifying selection. Phylogenetic analyses based on two gene datasets (15 core PCGs, and 15 core PCGs + rnl + rns) recovered identical and well-supported trees. In addition, cox1 was identified as a potential single-gene molecular marker for the phylogenetic analysis of relationships among Agaricomycetes species. This study provides the first report of mitogenomes from the Russulaceae family and facilitates the investigation of population genetics and evolution of other ectomycorrhizal fungi.