Peng Wu
Tian Yao2
Yuanhang Ren
Jinghua Ye
Qiang Li- 1Yunnan Plateau Characteristic Agricultural Industry Research Institute, Yunnan Agricultural University, Kunming, China
- 2School of Food and Biological Engineering, Chengdu University, Chengdu, China
- 3College of Information Science and Engineering, Chengdu University, Chengdu, China
- 4Panxi Featured Crops Research and Utilization Key Laboratory of Sichuan Province, Xichang University, Xichang, China
The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.
Introduction
Pisolithus species are ectomycorrhizal fungi with wide geographic distribution and host ranges, such as pine, oak and so on (Marx, 1977; Hitchcock et al., 2011). The symbiosis of Pisolithus species with trees promotes the carbon and nutrient cycling of the forest ecosystem (Cairney et al., 1989). In addition, the Pisolithus symbionts can promote the growth of plants, the absorption of mineral elements, and enhance the ability of plants to cope with abiotic stress, including drought, heavy metals and nutritional deficiencies (Jourand et al., 2014; Shi et al., 2017; Sebastiana et al., 2018). Up to now, more than ten species have been described in the genus Pisolithus. The limited morphological characteristics caused confusion in the classification of subspecies, varieties and species complex in the genus Pisolithus (Junghans et al., 1998; Hitchcock et al., 2003). The mitochondrial genome is a powerful tool to analyze the phylogenetic relationship of basidiomycetes (Basse, 2010). Up to now, six mitogenomes from the order Boletales have been reported, including four ectomycorrhizal Boletales species, Rhizopogon salebrosus, R. vinicolor (Li et al., 2019a), Paxillus involutus, P. rubicundulus (Li et al., 2020b), and two saprophytic Boletales species, Coniophora olivacea (MT375015), and C. puteana (MT375016). However, no complete mitogenome of Pisolithus species has been published.
Mitochondria play an important role in environmental adaptation, stress regulation and growth of eukaryotes (Ernster and Schatz, 1981). As the energy factory of cells, the genome of mitochondria mainly encodes respiratory chain related proteins, components of the ATPase and genes required for translation, such as rRNA and tRNA genes (McBride et al., 2006). Mutations in the mitogenome could lead to aging, disease, and even death of eukaryotic species (Ernster and Schatz, 1981; McBride et al., 2006; Gu et al., 2019). The mitogenome of basidiomycetes varies greatly in gene arrangement, genome size, intron type and number, repeat sequence content, and gene composition (Li et al., 2018d, 2020a; Sandor et al., 2018). Some of these characteristics can be used to infer the phylogenetic relationship of basidiomycetes and understand the origin or evolution of basidiomycetes (Richards et al., 1998; Li et al., 2015). In addition, the mitogenome of basidiomycetes contains several available molecular markers, which are widely used in phylogeny (Richards et al., 1998; Li et al., 2019c). So far, several nuclear genes, such as the internal transcribed spacer (ITS), the second largest subunit of RNA polymerase II (RPB2), beta-tubulin sequences, and translation elongation factor 1-α (tef1-α), have been widely used in the taxonomic classification and phylogeny of basidiomycetes (James et al., 2006). With more and more fungal nuclear genomes sequenced, the construction of fungal phylogeny based on whole genome sequence has become a trend (Li et al., 2018a). Mitogenome provides more genetic information than single nuclear gene molecular markers, and is easier to be obtained than the whole nuclear genome. So the mitogenome has become an effective supplementary tool to analyze fungal phylogeny. Compared with the mitogenome of animals, the mitogenome of fungi contains introns, and their types and numbers vary greatly among different fungal species (Kolesnikova et al., 2019; Li et al., 2020c). The number of available fungal mitogenomes, particularly basidiomycete mitogenomes, is far less than that of animals, or even less than the available nuclear genomes of fungi1. So far, two Pisolithus nuclear genomes have been reported, including Pisolithus microcarpus, and Pisolithus tinctorius (Kohler et al., 2015), but their mitogenomes have not been resolved.
In the present study, mitogenomes of two Pisolithus species, including P. microcarpus and P. tinctorius, were assembled and compared with other Boletales mitogenomes to assess the variations and conservation based on their gene arrangement, genome size, gene content and gene structure. Phylogenetic relationships of 79 basidiomycetes were revealed based on the combined mitochondrial datasets. Intron dynamics in cox1 genes of 79 basidiomycetes were also detected to understand the origin and evolution of introns in the phylum Basidiomycota. This study served as the first report on mitogenomes of the genus Pisolithus and the family Pisolithaceae, which will promote the understanding of the phylogeny, evolution and population genetics of Pisolithus species.
Materials and Methods
Pisolithus Mitogenomes Assembly and Annotation
The raw sequencing data of P. microcarpus and P. tinctorius used for mitogenomes assembly were downloaded from the Sequence Read Archive (SRA) database under the accession numbers of SRR3927183 and SRR3927210, respectively, which were submitted by the DOE Joint Genome Institute (Kohler et al., 2015). We conducted several quality control steps to generate clean data from the raw sequencing data, including removing adapter reads by ngsShoRT (Chen et al., 2014) and filtering reads with low-quality values by AdapterRemoval v 2 (Schubert et al., 2016). The two Pisolithus mitogenomes were assembled using the SPAdes 3.9.0 software based on the obtained clean reads (Bankevich et al., 2012). Gaps between contigs were filled using the MITObim V1.9 software (Hahn et al., 2013). The complete mitogenomes of P. microcarpus and P. tinctorius were obtained and further annotated according to the previously described methods (Li et al., 2018b). The protein-coding genes (PCGs), open reading frames (ORFs), introns, rRNA genes, and tRNA genes were first annotated by MFannot (Valach et al., 2014) and MITOS (Bernt et al., 2013). Then the annotated ORFs were verified or modified using the NCBI Open Reading Frame Finder (Coordinators, 2017), and further annotated by BLASTP searches against the NCBI non-redundant protein sequence database (Bleasby and Wootton, 1990). Intron-exon borders of PCGs were determined by the exonerate v2.2 software (Slater and Birney, 2005). tRNA genes in the two Pisolithus mitogenomes were also predicted using tRNAscan-SE v1.3.1 (Lowe and Chan, 2016). The graphical maps of the two Pisolithus mitogenomes were drawn using OGDraw v1.2 (Lohse et al., 2013).
Sequence Analysis of Pisolithus Mitogenomes
We calculated base compositions of the two Pisolithus mitogenomes and six Boletales mitogenomes using the DNASTAR Lasergene v7.12. Strand asymmetries of all the eight Boletales mitogenomes were assessed using the following formulas: AT skew = [A − T]/[A + T], and GC skew = [G − C]/[G + C]. Pairwise genetic distances between each pair of the 15 core PCGs (atp6, atp8, atp9, cob, cox1, cox2, cox3, nad1, nad2, nad3, nad4, nad4L, nad5, nad6, and rps3) in the eight Boletales mitogenomes tested were calculated using MEGA v6.06 (Caspermeyer, 2016) based on the Kimura-2-parameter (K2P) substitution model. We also detected the synonymous substitution rates (Ks) and non-synonymous substitution rates (Ka) for core PCGs in the eight Boletales mitogenomes using the DnaSP v6.10.01 software (Rozas et al., 2017). Gene collinearity analysis of the eight Boletales species was conducted using Mauve v2.4.0 software (Darling et al., 2004).
Repetitive Element Analysis
We conducted BLASTN searches (Chen et al., 2015) of the two Pisolithus mitogenomes against themselves to detect any interspersed repeats or intra-genomic duplications in the two Pisolithus mitogenomes, based on an E-value of <10–10. Tandem repeats (>10 bp in length) in the two Pisolithus mitogenomes were detected using the Tandem Repeats Finder (Benson, 1999). We also performed BlastN searches of the two Pisolithus mitogenomes against their published nuclear genomes (GCA_000827275.1; GCA_000827335.1) to detect if there were any gene fragments that natural transferred between the nuclear and mitochondrial genomes of the two Pisolithus species.
Comparative Mitogenomic Analysis and Intron Analysis
In the present study, we compared the genome sizes, gene numbers, intron numbers, base compositions, and gene arrangements of the eight Boletales mitogenomes reported to assess variations and conservativeness between different Boletales species. We assigned introns of cox1 genes in 79 published Basidiomycota mitogenomes into different position classes (Pcls) according to the method described by Ferandon et al. (2010). We first aligned the cox1 genes of 79 Basidiomycota mitogenomes with the reference cox1 gene of the medical fungus G. calidophilum (Li et al., 2019c) using Clustal W (Thompson et al., 1994). Introns were named according to the insert sites (nt) in corresponding reference gene. Each Pcl was constituted by introns inserted at the same position of corresponding reference cox1 gene (Cheng et al., 2021). Pcls with the same insert sites were considered as orthologous introns and had high sequence similarities (Ferandon et al., 2010).
Phylogenetic Analysis
In order to investigate the phylogenetic relationships of Boletales species and other Basidiomycota species, we constructed a phylogenetic tree of 79 species based on the combined mitochondrial gene datasets (15 core PCGs) (Li et al., 2021b). The mitogenome of Annulohypoxylon stygium from the phylum Ascomycota was used as an outgroup (Deng et al., 2018). We first aligned individual mitochondrial genes using MAFFT v7.037 software (Katoh et al., 2019), and then concatenated them into a combined mitochondrial gene dataset using SequenceMatrix v1.7.8 (Vaidya et al., 2011). We conducted partition homogeneity test to detect if there was potential phylogenetic conflict between different mitochondrial genes. PartitionFinder v2.1.1 software (Lanfear et al., 2017) was used to determine partitioning schemes for the mitochondrial gene set and best-fit models for phylogeny. The phylogenetic analysis was conducted using both maximum likelihood (ML) and bayesian inference (BI) methods (Li et al., 2018c). We used RAxML v 8.0.0 (Stamatakis, 2014) to conduct the ML analysis, and used MrBayes v3.2.6 (Ronquist et al., 2012) to perform the BI analysis. When conducted the ML analysis, bootstrap values (BS) were assessed through an ultrafast bootstrap approach with 10,000 replicates. BI analysis was conducted according to the previously described methods (Wang et al., 2020; Ye et al., 2020).
Data Availability
The complete mitogenomes of P. microcarpus and P. tinctorius were deposited in the GenBank database under the accession numbers of MT577034 and MT577035, respectively.
Results
Characterization and PCGs of Pisolithus Mitogenomes
The mitogenomes of P. microcarpus and P. tinctorius were all circular, with sizes of 43,990 and 44,054 bp, respectively (Figure 1). The GC contents of the P. microcarpus and P. tinctorius mitogenomes were 23.34 and 22.18%, respectively (Supplementary Table 1). Both the two Pisolithus species had negative AT skews and positive GC skews. There were 25 and 22 PCGs present in the P. microcarpus and P. tinctorius mitogenomes, respectively. Both the two mitogenomes contained a set of core PCGs, including atp6, atp8, atp9, cob, cox1, cox2, cox3, nad1, nad2, nad3, nad4, nad4L, nad5, nad6, and rps3 (Supplementary Table 2). In addition, we detected two plasmid-derived genes in the mitogenome of P. tinctorius, which encoded DNA polymerase and RNA polymerase. We also detected five and two PCGs with unknown functions in the P. microcarpus and P. tinctorius mitogenomes, respectively. Four and three introns were detected in the P. microcarpus and P. tinctorius mitogenomes, respectively, which were distributed in cox1 and cob genes of Pisolithus mitogenomes. Five and three intronic ORFs were detected in the P. microcarpus and P. tinctorius mitogenomes, respectively, which encoded two families of homing endonucleases, namely GIY YIG and LAGLIDADG endonucleases. The first intron of cox1 gene in P. microcarpus contained two intronic ORFs.
Figure 1. Circular maps of the two Pisolithus mitogenomes. Genes are represented by different colored blocks. Colored blocks outside each ring indicate that the genes are on the direct strand, while colored blocks within the ring indicates that the genes are located on the reverse strand. The inner graph is the GC content of mitochondrial sequences, and the circle inside the GC content graph marks the 50% threshold.
RNA Genes in Pisolithus Mitogenomes
Both the P. microcarpus and P. tinctorius mitogenomes contained two rRNA genes, namely the small subunit ribosomal RNA gene (rns) and the large subunit ribosomal RNA gene (rnl) (Supplementary Table 2). The average lengths of rns and rnl genes in the two Pisolithus mitogenomes were 1,973 and 3,768 bp, respectively. The mitogenome of P. microcarpus contained shorter rns and longer rnl genes than that of P. tinctorius.
We detected 26 and 25 tRNA genes in the mitogenomes of P. microcarpus and P. tinctorius, respectively (Supplementary Table 2). The mitogenome of P. microcarpus contained an additional trnI compared with the P. tinctorius mitogenome. All tRNAs detected in the two Pisolithus mitogenomes were folded into classical cloverleaf structures (Figure 2). The length of these tRNAs ranged from 71 bp to 86 bp, and the trnS gene was the largest. The extra stem/loop regions were found to be the main factors contributing to size variations of tRNA genes in Pisolithus. Ten of the 26 tRNAs shared in the two Pisolithus mitogenomes contained sites that varied between the two Pisolithus mitogenomes. We detected a total of 15 variable sites between the two Pisolithus mitogenomes, of which seven occurred on the D arm.
Figure 2. Putative secondary structures of tRNA genes in the two Pisolithus mitogenomes. tRNA genes with green color represent tRNAs shared by the two Pisolithus species, while tRNA genes with blue color represent the unique tRNAs in Pisolithus microcarpus. Residues conserved across the two mitogenomes are shown in green, while variable sites are shown in red. All genes are shown in order of occurrence in the mitogenome of Pisolithus microcarpus, starting from trnM.
Repeats of Pisolithus itogenomes
A total of two and four intra-genomic duplications were detected in the mitogenomes of P. microcarpus and P. tinctorius, respectively (Supplementary Table 3). The size of the four intra-genomic duplications ranged from 39 to 52 bp, and the pair-wise nucleotide similarities were between 79.49 and 100%. The largest intra-genomic duplications were observed in the coding regions of intronic ORFs, orf349 and orf325, in the P. tinctorius mitogenome. Intra-genomic duplications in the P. microcarpus mitogenomes were found located in the protein coding regions of atp6, as well as in the intergenic region between neighboring genes atp6 and trnM. Intra-genomic duplications accounted for 0.23 and 0.41% of the P. microcarpus and P. tinctorius mitogenomes, repectively. A total of 14 and 36 tandem repeats were detected in the mitogenomes of P. microcarpus and P. tinctorius, repectively, (Supplementary Table 4). The largest tandem repeat was found located in the intergenic region between nad6 and nad2 in the P. tinctorius mitogenome, with a size of 144 bp. Individual tandem repeat in the two Pisolithus mitogenomes contained 2 to 33 copies, with the most copy number observed in the P. tinctorius mitogenome. Tandem repeats accounted for 1.58 and 4.96% of the P. microcarpus and P. tinctorius mitogenomes, repectively.
We conducted BlastN searches of the two Pisolithus mitogenomes against their published nuclear genomes to detect any gene fragments that transferred between the mitochondrial and nuclear genomes. There were 4 and 22 aligned fragments found between mitochondrial and nuclear genomes of P. microcarpus and P. tinctorius, respectively (Supplementary Table 5). The length of these repetitive fragments ranged from 31 to 7,709 bp, and the pair-wise nucleotide similarities were between 93.10 and 100%. The largest aligned fragment was found located between rnl and trnM genes in the mitogenome of P. tinctorius. A total of 129 and 36,944 bp aligned fragments were detected in the mitogenomes of P. microcarpus and P. tinctorius, respectively. Large aligned fragments found between mitochondrial and nuclear genomes of the two Pisolithus species indicated that gene fragment transfer events may have occurred in the evolution of Pisolithus species.
Mitochondrial Gene Arrangement in Boletales Species
To date, six complete mitogenomes from the order Boletales have been reported (Li et al., 2019a). In the present study, we compared gene arrangements of 15 core PCGs and 2 rRNA genes in the eight Boletales species. Five out of the eight Boletales species had identical gene arrangements, including Rhizopogon vinicolor, Paxillus involutus, Paxillus rubicundulus, P. microcarpus, and P. tinctorius, which may represent the gene arrangement of the common ancestor of ectomycorrhizal Boletales species (Figure 3). However, we detected large-scale gene rearrangements in the other three Boletales species, including C. olivacea, Coniophora puteana, and Rhizopogon salebrosus. Compared with the gene arrangement of the putative ancestor of ectomycorrhizal Boletales, the three Boletales species were found underwent large-scale gene rearrangements, involving gene migrations, insertions, and inversions.
Figure 3. Mitochondrial gene arrangement analyses of eight Boletales species. All genes are shown in order of occurrence in the mitochondrial genome, starting from cox1. Fifteen core protein coding genes and two rRNA genes were included in the gene arrangement analysis. Genes with a green background indicate that they are conservative in Boletales species. The genes with red background indicate that they have gene rearrangements. The phylogenetic positions of the eight Boletales species were established using the Bayesian inference (BI) method and Maximum Likelihood (ML) method based on concatenated mitochondrial genes. Coli, C. olivacea; Cput, C. puteana; Pinv, P. involutus; Prub, P. rubicundulus; Pmic, P. microcarpus; Ptin, P. tinctorius; Rsal, R. salebrosus; Rvin, R. vinicolor.
A total of 28 homologous regions were detected between the 8 Boletales species using Mauve (Figure 4). The mitogenomes of P. microcarpus, P. tinctorius, and R. vinicolor were found containing the homologous region Y, which have never been detected in other Boletales species. The Y region contained some intergenetic sequences, which may be obtained from the same ancestor. Homologous region X was unique homologous region in C. puteana and P. tinctorius mitogenomes. The two Coniophora contained unique homologous regions G, H, I, J, K, M, and N, which were involved in orf276, orf113, and orf516 encoding DNA polymerases, and orf416, orf400, and orf302 encoding proteins with unknown functions. The results showed that the saprophytic Boletales species (Coniophora) differentiated from ectomycorrhizal Boletales species in gene content. Collinearity analysis indicated that Pisolithus, Paxillus, and Rhizopogon species showed a high degree of collinearity, while large-scale gene rearrangements were found between Coniophora species, indicating gene arrangement differentiated in Coniophora species.
Figure 4. Gene collinearity analysis of eight Boletales species using Mauve v2.4.0. Color blocks of the same color represent homologous regions between different mitogenomes, which are represented by different capital letters. The schematic diagram of Coniophora olivacea’s mitogenome is shown at the top of the picture.
Evolution of Core PCGs in Boletales
The largest K2P genetic distance was found in rps3 gene among the 15 core PCGs detected, which indicted the rps3 gene varied greatly between the 8 Boletales species (Figure 5). Large gene differentiations were also observed in the cox3 and nad2 genes between the 8 Boletales species tested. The lowest mean K2P genetic distances were observed in the atp8 and atp9 genes among the 15 core PCGs, indicating that the two genes were highly conserved between the 8 Boletales species.
Figure 5. Genetic analysis of 15 core protein coding genes conserved in the 8 Boletales mitogenomes. K2P, pairwise genetic distances between each pair of the 15 core PCGs in the 8 Boletales mitogenomes based on the Kimura-2-parameter model. Ka, non-synonymous substitution site; Ks, synonymous substitution site. Species and NCBI accession number used for genetic analysis in the present study are listed in Supplementary Table 6.
The highest mean non-synonymous substitution rate (Ka) was observed in the rps3 gene among the 15 core PCGs detected, followed by the cox3 gene (Figure 5), while the lowest mean Ka value was observed in the atp9 gene among the 15 core PCGs from the eight Boletales species. The cox1 genes in the 8 Boletales species were found exhibited the highest mean synonymous substitution rate (Ks), while the rps3 gene exhibited the lowest mean Ks value among the 15 PCGs detected. We found the Ka/Ks values for 14 of the 15 core PCGs were less than 1, indicating that the 14 core genes were subjected to purifying selection in evolution. However, we found the Ka/Ks value of the rps3 was >1 between some Boletales species, including between P. microcarpus and Paxillus involutus, between P. microcarpus and Paxillus rubicundulus, between P. tinctorius and C. olivacea, between P. tinctorius and C. puteana, as well as between Rhizopogon vinicolor and R. salebrosus, indicating that the rps3 gene may be under strong positive selection pressure in some Boletales species.
Intron Dynamics in cox1 Genes of Basidiomycota Species
Intron dynamics in cox1 genes of 79 Basidiomycota species were analyzed. Species included in the intron analysis accounted for over 2/3 of Basidiomycota mitogenomes available in the NCBI database. We detected 1,071 introns in the 79 Basidiomycota species, and most of them belonged to the group I. The number of intron in individual Basidiomycota species ranged from 0 to 46, and the cox1 gene was found containing the largest number of introns among all the host genes detected in Basidiomycota. A total of 488 introns were found located in cox1 genes of the 79 Basidiomycota species, which occupied 45.56% of the introns found in the mitogenomes of the 79 examined basidiomycetes.
We further classified introns of cox1 genes into different position classes (Pcls). In the present study, we detected 47 Pcls in the cox1 genes of 79 Basidiomycota species (Figure 6). Pcls present in more than 1/5 of Basidiomycota species were considered as common Pcls in Basidiomycota, while introns detected in less than 1/5 of Basidiomycota species were considered as rare introns. In the present study, 12 common Pcls and 35 rare Pcls were detected in the 79 Basidiomycota species. The P383 was the most common introns in Basidiomycota, which distributed in 40 of the 79 basidiomycete species. P706 and P1107 were also widely distributed in Basidiomycota, which could be detected in 35 and 34 basidiomycete species, respectively. However, 14 Pcls could only be detected in only one of the 79 basidiomycete species, including P166, P193, and P218, P309, P318, etc. A total of 20 Pcls were detected in the 8 Boletales species, of which 10 were common introns in Basidiomycota. The quantity and class of introns in Boletales species varied greatly, indicating intron loss/gain occurred in the evolution of Boletales. The Pcls of the two novel mitogenomes we assembled in this study varied, and P717 present in P. tinctorius has not been detected in other Boletales species. However, the novel Pcl (P717) in Boletales was also detected in Lyophyllum shimeji (Li et al., 2019b), Ganoderma calidophilum (Li et al., 2019c), and Cryptococcus gattii (Gomes et al., 2018) from other orders, indicating potential horizontal gene transfer events.
Figure 6. Position class (Pcl) information of cox1 genes in the 79 Basidiomycota species. Introns in cox1 genes of 79 published Basidiomycota mitogenomes were classified into different position classes (Pcls) using the cox1 gene of G. calidophilum as the reference. Each Pcl was constituted by introns inserted at the same position of corresponding cox1 gene and named according to its insertion site in the aligned corresponding reference sequence (nt). The Pcls present in more than 1/5 of Basidiomycota species were considered as common Pcls in Basidiomycota, while introns detected in less than 1/5 of Basidiomycota species were considered to be rare introns. Phylogenetic positions of the 79 Basidiomycota species were established using the Bayesian inference (BI) method and Maximum Likelihood (ML) method based on concatenated mitochondrial genes. Species information is shown in Supplementary Table 6.
Comparative Mitogenome Analysis and Phylogenetic Analysis
Comparative mitogenomic analysis showed that the two Pisolithus mitogenomes were smaller than the four mitogenomes from the genera Coniophora and Rhizopogon, and larger than the two mitogenomes from the genus Paxillus (Supplementary Table 1). The average GC content of the two Pisolithus mitogenomes (22.76%) was lower that of the two species from the genus Coniophora, but larger than that of the other four Boletales species. All the 8 Boletales species contained negative AT skews and positive GC skews. Comparative mitogenome analysis indicated that the number of non-intronic ORFs, plasmid-derived genes, PCGs with unknown functions, introns, intronic ORFs, and tRNA genes varied between Boletales species. All the eight Boletales species contained two rRNA genes. The two Pisolithus mitogenomes were found containing the lowest intra-genomic repeats in their mitogenomes than in the other six Boletales species. The proportions of introns, protein coding region, intergenic region, and RNA region were also variable in different Boletales species. The correlation analysis of mitogenome size and intron number of 79 basidiomycetes showed that the mitogenome size had a significant correlation with the number of introns (p < 0.01), with the pearson correlation coefficient 0.808 and the spearman correlation coefficient 0.840 (Figure 7).
Figure 7. Pearson and spearman correlation analyses of mitogenome size and intron number of 79 basidiomycetes.
An identical and well-supported phylogenetic tree for 79 basidiomycete species was obtained using both Maximum likelihood (ML) and Bayesian inference (BI) methods (Figure 8). The combined mitochondrial gene dataset (15 core PCGs) was used to conduct the phylogenetic analysis. We found all major clades within the phylogenetic tree had good supported values (BPP ≥ 0.96; BS ≥ 98). The 79 basidiomycete species could be divided into 15 major clades in the phylogenetic tree, corresponding to the orders Agaricales, Boletales, Cantharellales, Gomphales, Hymenochaetales, Microbotryales, Microstromatales, Polyporales, Pucciniales, Russulales, Sporidiobolales, Tilletiales, Tremellales, Trichosporonales, and Ustilaginales (Supplementary Table 6). The eight Boletales species could be assigned into two groups, wherein the first group comprised two Coniophora species, and the second group comprised six species from several genera, including Rhizopogon, Paxillus, and Pisolithus. The results indicated that the ectomycorrhizal Boletales species may have common origin. The results also showed that the two Pisolithus species exhibited a close relationship with Paxillus species.
Figure 8. Molecular phylogeny of 79 Basidiomycota species based on Bayesian inference (BI) and Maximum likelihood (ML) analysis of 15 protein coding genes and two rRNA genes. Support values are Bayesian posterior probabilities (BPP, before slash) and bootstrap values (BS, after slash). The asterisk indicates that the BPP value is one and the BS value is 100 of the branch. Species and NCBI accession numbers for mitogenomes used in the phylogenetic analysis are provided in Supplementary Table 6.
Discussion
Variations of Gene Content in Boletales Mitogenomes
It has been reported that the mitogenome of eukaryotes is obtained from bacteria through endosymbiosis (Lang et al., 1999; Munoz-Gomez et al., 2017). In the process of evolution, most of the mitochondrial genes have been transferred to the nuclear genomes (Adams et al., 2002; Adams and Palmer, 2003). However, most animal and fungal mitogenomes retain a whole set of core PCGs in the evolutionary process, which is believed to facilitate the local control of oxidative phosphorylation (Allen, 2015; Bjorkholm et al., 2015). In the present study, we found that these core PCGs varied in length and base composition between different Boletales species. The rps3 gene differentiated greatly between different Boletales species, while atp8 and atp9 genes were quite conservative between Boletales species, indicating that different mitochondrial genes had different evolution rates in Boletales. The rps3 gene of basidiomycetes is an ancient gene which come from alphaproteobacteria and was relocated to another region of the mitogenome (Korovesi et al., 2018). Previous studies found that the rps3 gene was probably implicated in the assembly of the mitochondrial small (37S) ribosomal subunit (Seif et al., 2005). In the present study, we found that the rps3 gene may have been subjected to positive selection between some Boletales species. In addition to core PCGs, we also detected a series of non-conserved PCGs in Boletales species, including plasmid-derived genes, intronic ORFs and PCGs with unknown functions. The number and species of these non-conserved PCGs varied greatly in different Boletales species. Plasmid-derived genes were thought to be obtained from mitochondrial plasmids (Himmelstrand et al., 2014), which were considered as self-replicating genetic elements. Mitochondrial plasmids were free-standing in some Basidiomycota species (Yui et al., 1988; Nakai et al., 2000), and integrated into mitogenomes in some other Basidiomycota species (Formighieri et al., 2008; Ferandon et al., 2013). The dynamic changes of mitochondrial plasmids also play an important role in the evolution of fungal mitogenome (Himmelstrand et al., 2014; Li et al., 2020a; Wang et al., 2020). In addition, different numbers of PCGs with unknown functions have been detected in the eight Boletales mitogenomes. Further studies are needed to fully understand the role of these unknown functional PCGs in species differentiation and ecological adaptation of Boletales. We also detected 7 out of 15 variable sites in the D arm of tRNA genes, which showed that the D arm was highly variable in the two Pisolithus mitogenomes.
Gene Rearrangements in Boletales Mitogenomes
Mitochondrial gene arrangements provide abundant information for understanding the evolutionary status of species (Beaudet et al., 2013; Zhang et al., 2017; Li et al., 2018b). In this study, we found that five out of the eight Boletales species had identical mitochondrial gene arrangement, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. The identical mitochondrial gene arrangement was found distributed in species from the genera, Rhizopogon (Li et al., 2019a), Paxillus (Li et al., 2020b), and Pisoliths, indicating close phylogenetic relationships between these species. Large-scale gene rearrangements were detected in the two saprophytic Boletales species (Coniophora) and the ectomycorrhizal fungus R. salebrosus. A total of 4, 14 and 5 out of the 17 genes rearranged in mitogenomes of C. olivacea, C. puteana, and R. salebrosus, respectively, compared with the mitochondrial gene arrangement of putative ectomycorrhizal Boletales ancestors. Up to now, the mechanism of mitochondrial gene rearrangement in fungi has not been fully revealed. Previous studies have shown that the accumulation of repetitive sequences in fungal mitogenomes may contribute to the mitochondrial gene rearrangements in fungi (Aguileta et al., 2014). However, in this study, we found weak correlations between repeat content and gene rearrangement in Boletales. Few repeats were detected in the mitogenome of R. salebrosus (Li et al., 2019a), but a large-scale gene rearrangements were found in it. Therefore, there may be other models or mechanisms for Boletales mitochondrial genome rearrangement.
Introns Dynamics in cox1 Genes of Basidiomycota Mitogenomes
As mobile genetic elements in fungal mitogenome, intron dynamics play an important role in altering organization or size of fungal mitogenomes (Himmelstrand et al., 2014; Gomes et al., 2018; Fan et al., 2019; Li et al., 2021a). In the present study, we found the number and class of basidiomycete introns varied greatly, even between closely related basidiomycete species, and a significant correlation between the mitogenome size and intron number of 79 basidiomycetes were detected. A total of 1,071 introns were detected in the 79 basidiomycete species, with each species containing 0–46 introns. The introns of basidiomycetes were unevenly distributed in the core PCGs or rRNA genes of basidiomycetes, of which cox1 gene was the largest host gene of basidiomycete introns. A total of 47 Pcls were detected in the 79 basidiomycete species, and the P383 was the most widely distributed intron class. We found the classes and numbers of introns in the eight Boletales species tested varied greatly, which indicated intron loss/gain events occurred in the evolution of Boletales. Some introns of the Boletales were lost in the evolutionary process, resulting in the relatively small size of Boletales in Basidiomycota. The intron P717 intron was only found in Pisoliths tinctorius of the eight Boletales species, while homologous introns were detected in distant species (Li et al., 2019c), indicating potential intron transfer events.
Reconstruction of Basidiomycete Phylogeny Based on Mitochondrial Gene
Basidiomycetes, the largest mushroom-forming fungal group on the earth, play an important role in ecosystem, industry and medicine (Cahan and Kennell, 2005; Maciel et al., 2010; Martinkova et al., 2016). Understanding the origin, evolutionary and phylogenetic relationships of basidiomycetes will lay a foundation for the effective utilization of basidiomycetes. However, it is difficult to classify or identify basidiomycetes only by morphology, which is mainly due to the changeable and overlapping morphological characteristics or convergent evolution of basidiomycetes (Li et al., 2019b). The introduction of molecular markers promotes the understanding of basidiomycete origin and evolution. Up to now, the rRNA ITS sequence is the most commonly used molecular marker in basidiomycete phylogeny (Schoch et al., 2012). Nuclear genomes of basidiomycetes were also used for phylogeny of basidiomycetes (Cornell et al., 2007; Liu et al., 2017). The mitogenome of basidiomycetes can provide more genetic information than ITS sequence, and it is easier to be obtained than nuclear genome. Therefore, mitogenome of basidiomycetes may become a universal molecular marker to analyze the origin and phylogeny of basidiomycetes (Sandor et al., 2018; Wang et al., 2018; Li et al., 2020c; Wang and Xu, 2020). In the present study, more than 2/3 of basidiomycete mitogenomes available in the public database (see text footnote 1) were included in the phylogenetic study. A well-supported phylogenetic tree for the 79 basidiomycetes was obtained based on mitochondrial gene dataset, which showed that mitochondrial genome is an effective tool to analyze phylogenetic relationship of basidiomycetes. More mitogenomes of basidiomycetes need to be resolved to reconstruct the origin and evolution of basidiomycetes in the future.
Data Availability Statement
The complete mitogenomes of P. microcarpus and P. tinctorius were deposited in the GenBank database under the accession numbers of MT577034 and MT577035, respectively.
Author Contributions
MG and QL: conceived and designed experiments. PW, TY, JY, YQ, and YR: analyzed the data. PW and QL: wrote and review the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the project of protection, breeding and high value utilization of Plateau Characteristic edible fungi resources (2019ZG00906 and 2019ZG00906-04).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.583129/full#supplementary-material
Footnotes
References
Adams, K. L., and Palmer, J. D. (2003). Evolution of mitochondrial gene content: gene loss and transfer to the nucleus. Mol. Phylogenet. Evol. 29, 380–395. doi: 10.1016/s1055-7903(03)00194-5
Adams, K. L., Qiu, Y. L., Stoutemyer, M., and Palmer, J. D. (2002). Punctuated evolution of mitochondrial gene content: high and variable rates of mitochondrial gene loss and transfer to the nucleus during angiosperm evolution. Proc. Natl. Acad. Sci. U.S.A. 99, 9905–9912. doi: 10.1073/pnas.042694899
Aguileta, G., de Vienne, D. M., Ross, O. N., Hood, M. E., Giraud, T., Petit, E., et al. (2014). High variability of mitochondrial gene order among fungi. Genome Biol. Evol. 6, 451–465. doi: 10.1093/gbe/evu028
Allen, J. F. (2015). Why chloroplasts and mitochondria retain their own genomes and genetic systems: colocation for redox regulation of gene expression. Proc. Natl. Acad. Sci. U.S.A. 112, 10231–10238. doi: 10.1073/pnas.1500012112
Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021
Basse, C. W. (2010). Mitochondrial inheritance in fungi. Curr. Opin. Microbiol. 13, 712–719. doi: 10.1016/j.mib.2010.09.003
Beaudet, D., Terrat, Y., Halary, S., de la Providencia, I. E., and Hijri, M. (2013). Mitochondrial genome rearrangements in glomus species triggered by homologous recombination between distinct mtDNA haplotypes. Genome Biol. Evol. 5, 1628–1643. doi: 10.1093/gbe/evt120
Benson, G. (1999). Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573–580. doi: 10.1093/nar/27.2.573
Bernt, M., Donath, A., Juhling, F., Externbrink, F., Florentz, C., Fritzsch, G., et al. (2013). MITOS: improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol. 69, 313–319. doi: 10.1016/j.ympev.2012.08.023
Bjorkholm, P., Harish, A., Hagstrom, E., Ernst, A. M., and Andersson, S. G. (2015). Mitochondrial genomes are retained by selective constraints on protein targeting. Proc. Natl. Acad. Sci. U.S.A. 112, 10154–10161. doi: 10.1073/pnas.1421372112
Bleasby, A. J., and Wootton, J. C. (1990). Construction of validated, non-redundant composite protein sequence databases. Protein Eng. 3, 153–159. doi: 10.1093/protein/3.3.153
Cahan, P., and Kennell, J. C. (2005). Identification and distribution of sequences having similarity to mitochondrial plasmids in mitochondrial genomes of filamentous fungi. Mol. Genet. Genomics 273, 462–473. doi: 10.1007/s00438-005-1133-x
Cairney, J. W. G., Ashford, A. E., and Allaway, W. G. (1989). Distribution of photosynthetically fixed carbon within root systems of Eucalyptus pilularis plants ectomycorrhizal with Pisolithus tinctorius. New Phytol. 112, 495–500. doi: 10.1111/j.1469-8137.1989.tb00343.x
Caspermeyer, J. (2016). MEGA evolutionary software re-engineered to handle Today’s big data demands. Mol. Biol. Evol. 33:1887.
Chen, C., Khaleel, S. S., Huang, H., and Wu, C. H. (2014). Software for pre-processing Illumina next-generation sequencing short read sequences. Source Code Biol. Med. 9:8.
Chen, Y., Ye, W., Zhang, Y., and Xu, Y. (2015). High speed BLASTN: an accelerated MegaBLAST search tool. Nucleic Acids Res. 43, 7762–7768. doi: 10.1093/nar/gkv784
Cheng, J., Luo, Q., Ren, Y. H., Luo, Z., Liao, W. L., Wang, X., et al. (2021). Panorama of intron dynamics and gene rearrangements in the phylum Basidiomycota as revealed by the complete mitochondrial genome of Turbinellus floccosus. Appl. Microbiol. Biotechnol. 105, 2017–2032. doi: 10.1007/s00253-021-11153-w
Coordinators, N. R. (2017). Database resources of the National Center for Biotechnology Informati. Nucleic Acids Res. 47, D23–D28.
Cornell, M. J., Alam, I., Soanes, D. M., Wong, H. M., Hedeler, C., Paton, N. W., et al. (2007). Comparative genome analysis across a kingdom of eukaryotic organisms: specialization and diversification in the fungi. Genome Res. 17, 1809–1822. doi: 10.1101/gr.6531807
Darling, A. C., Mau, B., Blattner, F. R., and Perna, N. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403. doi: 10.1101/gr.2289704
Deng, Y., Hsiang, T., Li, S., Lin, L., Wang, Q., Chen, Q., et al. (2018). Comparison of the mitochondrial genome sequences of six Annulohypoxylon stygium isolates suggests short fragment insertions as a potential factor leading to larger genomic size. Front. Microbiol. 9:2079. doi: 10.3389/fmicb.2018.02079
Fan, W. W., Zhang, S., and Zhang, Y. J. (2019). The complete mitochondrial genome of the Chan-hua fungus Isaria cicadae: a tale of intron evolution in Cordycipitaceae. Environ. Microbiol. 21, 864–879. doi: 10.1111/1462-2920.14522
Ferandon, C., Moukha, S., Callac, P., Benedetto, J. P., Castroviejo, M., and Barroso, G. (2010). The Agaricus bisporus cox1 Gene: the longest mitochondrial gene and the largest reservoir of mitochondrial group I Introns. PLoS One 5:e14048. doi: 10.1371/journal.pone.0014048
Ferandon, C., Xu, J., and Barroso, G. (2013). The 135 kbp mitochondrial genome of Agaricus bisporus is the largest known eukaryotic reservoir of group I introns and plasmid-related sequences. Fungal Genet. Biol. 55, 85–91. doi: 10.1016/j.fgb.2013.01.009
Formighieri, E. F., Tiburcio, R. A., Armas, E. D., Medrano, F. J., Shimo, H., Carels, N., et al. (2008). The mitochondrial genome of the phytopathogenic basidiomycete Moniliophthora perniciosa is 109 kb in size and contains a stable integrated plasmid. Mycol. Res. 112, 1136–1152. doi: 10.1016/j.mycres.2008.04.014
Gomes, F., Arantes, T. D., Fernandes, J. A. L., Ferreira, L. C., Romero, H., Bosco, S. M. G., et al. (2018). Polymorphism in Mitochondrial Group I Introns among Cryptococcus neoformans and Cryptococcus gattii genotypes and its association with drug susceptibility. Front. Microbiol. 9:86.
Gu, X., Kang, X., and Liu, J. (2019). Mutation signatures in germline mitochondrial genome provide insights into human mitochondrial evolution and disease. Hum. Genet. 138, 613–624. doi: 10.1007/s00439-019-02009-5
Hahn, C., Bachmann, L., and Chevreux, B. (2013). Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads–a baiting and iterative mapping approach. Nucleic Acids Res. 41:e129. doi: 10.1093/nar/gkt371
Himmelstrand, K., Olson, A., Brandstrom Durling, M., Karlsson, M., and Stenlid, J. (2014). Intronic and plasmid-derived regions contribute to the large mitochondrial genome sizes of Agaricomycetes. Curr. Genet. 60, 303–313. doi: 10.1007/s00294-014-0436-z
Hitchcock, C. J., Chambers, S. M., Anderson, I. C., and Cairney, J. W. (2003). Development of markers for simple sequence repeat-rich regions that discriminate between Pisolithus albus and P. microcarpus. Mycol. Res. 107, 699–706. doi: 10.1017/s0953756203007871
Hitchcock, C. J., Chambers, S. M., and Cairney, J. W. (2011). Genetic population structure of the ectomycorrhizal fungus Pisolithus microcarpus suggests high gene flow in south-eastern Australia. Mycorrhiza 21, 131–137. doi: 10.1007/s00572-010-0317-3
James, T. Y., Kauff, F., Schoch, C. L., Matheny, P. B., Hofstetter, V., Cox, C. J., et al. (2006). Reconstructing the early evolution of Fungi using a six-gene phylogeny. Nature 443, 818–822.
Jourand, P., Hannibal, L., Majorel, C., Mengant, S., Ducousso, M., and Lebrun, M. (2014). Ectomycorrhizal Pisolithus albus inoculation of Acacia spirorbis and Eucalyptus globulus grown in ultramafic topsoil enhances plant growth and mineral nutrition while limits metal uptake. J. Plant Physiol. 171, 164–172. doi: 10.1016/j.jplph.2013.10.011
Junghans, D. T., Gomes, E. A., Guimaraes, W. V., Barros, E. G., and Araujo, E. F. (1998). Genetic diversity of the ectomycorrhizal fungus Pisolithus tinctorius based on RAPD-PCR analysis. Mycorrhiza 7, 243–248. doi: 10.1007/s005720050187
Katoh, K., Rozewicki, J., and Yamada, K. D. (2019). MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. 20, 1160–1166. doi: 10.1093/bib/bbx108
Kohler, A., Kuo, A., Nagy, L. G., Morin, E., Barry, K. W., Buscot, F., et al. (2015). Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. Nat. Genet. 47, 410–415. doi: 10.1038/ng.3223
Kolesnikova, A. I., Putintseva, Y. A., Simonov, E. P., Biriukov, V. V., Oreshkova, N. V., Pavlov, I. N., et al. (2019). Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species. BMC Genomics 20:351. doi: 10.1186/s12864-019-5732-z
Korovesi, A. G., Ntertilis, M., and Kouvelis, V. N. (2018). Mt-rps3 is an ancient gene which provides insight into the evolution of fungal mitochondrial genomes. Mol. Phylogenet. Evol. 127, 74–86. doi: 10.1016/j.ympev.2018.04.037
Lanfear, R., Frandsen, P. B., Wright, A. M., Senfeld, T., and Calcott, B. (2017). PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol. 34, 772–773.
Lang, B. F., Gray, M. W., and Burger, G. (1999). Mitochondrial genome evolution and the origin of eukaryotes. Annu. Rev. Genet. 33, 351–397. doi: 10.1146/annurev.genet.33.1.351
Li, H., Shao, R., Song, N., Song, F., Jiang, P., Li, Z., et al. (2015). Higher-level phylogeny of paraneopteran insects inferred from mitochondrial genome sequences. Sci. Rep. 5:8527.
Li, H., Wu, S., Ma, X., Chen, W., Zhang, J., Duan, S., et al. (2018a). The Genome Sequences of 90 Mushrooms. Sci. Rep. 8:9982.
Li, Q., Chen, C., Xiong, C., Jin, X., Chen, Z., and Huang, W. (2018b). Comparative mitogenomics reveals large-scale gene rearrangements in the mitochondrial genome of two Pleurotus species. Appl. Microbiol. Biotechnol. 102, 6143–6153. doi: 10.1007/s00253-018-9082-6
Li, Q., He, X., Ren, Y., Xiong, C., Jin, X., Peng, L., et al. (2020a). Comparative Mitogenome analysis reveals mitochondrial genome differentiation in Ectomycorrhizal and Asymbiotic Amanita Species. Front. Microbiol. 11:1382. doi: 10.3389/fmicb.2020.01382
Li, Q., Li, L. J., Feng, H. Y., Tu, W. Y., Bao, Z. J., Xiong, C., et al. (2021a). Characterization of the complete mitochondrial genome of Basidiomycete yeast Hannaella oryzae: intron evolution, gene rearrangement and its phylogeny. Front. Microbiol. 12:646567. doi: 10.3389/fmicb.2021.646567
Li, Q., Liao, M., Yang, M., Xiong, C., Jin, X., Chen, Z., et al. (2018c). Characterization of the mitochondrial genomes of three species in the ectomycorrhizal genus Cantharellus and phylogeny of Agaricomycetes. Int. J. Biol. Macromol. 118, 756–769. doi: 10.1016/j.ijbiomac.2018.06.129
Li, Q., Ren, Y., Shi, X., Peng, L., Zhao, J., Song, Y., et al. (2019a). Comparative mitochondrial genome analysis of two Ectomycorrhizal fungi (Rhizopogon) reveals dynamic changes of intron and phylogenetic relationships of the Subphylum Agaricomycotina. Int. J. Mol. Sci. 20:5167. doi: 10.3390/ijms20205167
Li, Q., Ren, Y., Xiang, D., Shi, X., Zhao, J., Peng, L., et al. (2020b). Comparative mitogenome analysis of two ectomycorrhizal fungi (Paxillus) reveals gene rearrangement, intron dynamics, and phylogeny of basidiomycetes. IMA Fungus 11:12.
Li, Q., Wang, Q., Chen, C., Jin, X., Chen, Z., Xiong, C., et al. (2018d). Characterization and comparative mitogenomic analysis of six newly sequenced mitochondrial genomes from ectomycorrhizal fungi (Russula) and phylogenetic analysis of the Agaricomycetes. Int. J. Biol. Macromol. 119, 792–802. doi: 10.1016/j.ijbiomac.2018.07.197
Li, Q., Wang, Q., Jin, X., Chen, Z., Xiong, C., Li, P., et al. (2019b). Characterization and comparison of the mitochondrial genomes from two Lyophyllum fungal species and insights into phylogeny of Agaricomycetes. Int. J. Biol. Macromol. 121, 364–372. doi: 10.1016/j.ijbiomac.2018.10.037
Li, Q., Wu, P., Li, L., Feng, H., Tu, W., Bao, Z., et al. (2021b). The first eleven mitochondrial genomes from the ectomycorrhizal fungal genus (Boletus) reveal intron loss and gene rearrangement. Int. J. Biol. Macromol. 172, 560–572. doi: 10.1016/j.ijbiomac.2021.01.087
Li, Q., Xiang, D., Wan, Y., Wu, Q., Wu, X., Ma, C., et al. (2019c). The complete mitochondrial genomes of five important medicinal Ganoderma species: features, evolution, and phylogeny. Int. J. Biol. Macromol. 139, 397–408. doi: 10.1016/j.ijbiomac.2019.08.003
Li, Q., Yang, L., Xiang, D., Wan, Y., Wu, Q., Huang, W., et al. (2020c). The complete mitochondrial genomes of two model ectomycorrhizal fungi (Laccaria): features, intron dynamics and phylogenetic implications. Int. J. Biol. Macromol. 145, 974–984. doi: 10.1016/j.ijbiomac.2019.09.188
Liu, X. B., Zhang, Q. Q., Xia, X. Y., Liu, X. Y., Ge, L., and Yang, L. (2017). Phylogenetic relationships of ascomycetes and basidiomycetes based on comparative genomics analysis. Genes Genomics 39, 1307–1316. doi: 10.1007/s13258-017-0595-7
Lohse, M., Drechsel, O., Kahlau, S., and Bock, R. (2013). OrganellarGenomeDRAW–a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 41, W575–W581.
Lowe, T. M., and Chan, P. P. (2016). tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 44, W54–W57.
Maciel, M. J. M., Silva, A. C. E., and Ribeiro, H. C. T. (2010). Industrial and biotechnological applications of ligninolytic enzymes of the basidiomycota: a review. Electron. J. Biotechnol. 13, 14–15.
Martinkova, L., Kotik, M., Markova, E., and Homolka, L. (2016). Biodegradation of phenolic compounds by Basidiomycota and its phenol oxidases: a review. Chemosphere 149, 373–382. doi: 10.1016/j.chemosphere.2016.01.022
Marx, D. H. (1977). Tree host range and world distribution of the extomycorrhizal fungus Pisolithus tinctorius. Can. J. Microbiol. 23, 217–223. doi: 10.1139/m77-033
McBride, H. M., Neuspiel, M., and Wasiak, S. (2006). Mitochondria: more than just a powerhouse. Curr. Biol. 16, R551–R560.
Munoz-Gomez, S. A., Wideman, J. G., Roger, A. J., and Slamovits, C. H. (2017). The origin of mitochondrial cristae from Alphaproteobacteria. Mol. Biol. Evol. 34, 943–956.
Nakai, R., Sen, K., Kurosawa, S., and Shibai, H. (2000). Basidiomycetous fungus Flammulina velutipes harbors two linear mitochondrial plasmids encoding DNA and RNA polymerases. FEMS Microbiol. Lett. 190, 99–102. doi: 10.1111/j.1574-6968.2000.tb09269.x
Richards, M. B., Macaulay, V. A., Bandelt, H. J., and Sykes, B. C. (1998). Phylogeography of mitochondrial DNA in western Europe. Ann. Hum. Genet. 62, 241–260. doi: 10.1046/j.1469-1809.1998.6230241.x
Ronquist, F., Teslenko, M., van der Mark, P., Ayres, D. L., Darling, A., Hohna, S., et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542. doi: 10.1093/sysbio/sys029
Rozas, J., Ferrer-Mata, A., Sanchez-DelBarrio, J. C., Guirao-Rico, S., Librado, P., Ramos-Onsins, S. E., et al. (2017). DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299–3302. doi: 10.1093/molbev/msx248
Sandor, S., Zhang, Y., and Xu, J. (2018). Fungal mitochondrial genomes and genetic polymorphisms. Appl. Microbiol. Biotechnol. 102, 9433–9448. doi: 10.1007/s00253-018-9350-5
Schoch, C. L., Seifert, K. A., Huhndorf, S., Robert, V., Spouge, J. L., Levesque, C. A., et al. (2012). Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc. Natl. Acad. Sci. U.S.A. 109, 6241–6246.
Schubert, M., Lindgreen, S., and Orlando, L. (2016). AdapterRemoval v2: rapid adapter trimming, identification, and read merging. BMC Res. Notes 9:88. doi: 10.1186/s13104-016-1900-2
Sebastiana, M., da Silva, A. B., Matos, A. R., Alcantara, A., Silvestre, S., and Malho, R. (2018). Ectomycorrhizal inoculation with Pisolithus tinctorius reduces stress induced by drought in cork oak. Mycorrhiza 28, 247–258. doi: 10.1007/s00572-018-0823-2
Seif, E., Leigh, J., Liu, Y., Roewer, I., Forget, L., and Lang, B. F. (2005). Comparative mitochondrial genomics in zygomycetes: bacteria-like RNase P RNAs, mobile elements and a close source of the group I intron invasion in angiosperms. Nucleic Acids Res. 33, 734–744. doi: 10.1093/nar/gki199
Shi, L., Wang, J., Liu, B., Nara, K., Lian, C., Shen, Z., et al. (2017). Ectomycorrhizal fungi reduce the light compensation point and promote carbon fixation of Pinus thunbergii seedlings to adapt to shade environments. Mycorrhiza 27, 823–830. doi: 10.1007/s00572-017-0795-7
Slater, G. S., and Birney, E. (2005). Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics 6:31. doi: 10.1186/1471-2105-6-31
Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033
Thompson, J. D., Higgins, D. G., and Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673–4680. doi: 10.1093/nar/22.22.4673
Vaidya, G., Lohman, D. L., and Meier, R. (2011). SequenceMatrix: concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics 27, 171–180. doi: 10.1111/j.1096-0031.2010.00329.x
Valach, M., Burger, G., Gray, M. W., and Lang, B. F. (2014). Widespread occurrence of organelle genome-encoded 5S rRNAs including permuted molecules. Nucleic Acids Res. 42, 13764–13777. doi: 10.1093/nar/gku1266
Wang, L., Zhang, S., Li, J. H., and Zhang, Y. J. (2018). Mitochondrial genome, comparative analysis and evolutionary insights into the entomopathogenic fungus Hirsutella thompsonii. Environ. Microbiol. 20, 3393–3405. doi: 10.1111/1462-2920.14379
Wang, X., Song, A., Wang, F., Chen, M., Li, X., Li, Q., et al. (2020). The 206 kbp mitochondrial genome of Phanerochaete carnosa reveals dynamics of introns, accumulation of repeat sequences and plasmid-derived genes. Int. J. Biol. Macromol. 162, 209–219. doi: 10.1016/j.ijbiomac.2020.06.142
Wang, Y., and Xu, J. (2020). Mitochondrial genome polymorphisms in the human Pathogenic Fungus Cryptococcus neoformans. Front. Microbiol. 11:706. doi: 10.3389/fmicb.2020.00706
Ye, J., Cheng, J., Ren, Y., Liao, W., and Li, Q. (2020). The First Mitochondrial Genome for Geastrales (Sphaerobolus stellatus) reveals intron dynamics and large-scale gene rearrangements of Basidiomycota. Front. Microbiol. 11:1970. doi: 10.3389/fmicb.2020.01970
Yui, Y., Katayose, Y., and Shishido, K. (1988). Two linear plasmid-like DNA elements simultaneously maintained in Pleurotus ostreatus. Biochim. Biophys. Acta 951, 53–60. doi: 10.1016/0167-4781(88)90024-3
Keywords: mitochondrial genome, plasmid, intron, gene rearrangement, evolution, phylogenetic analysis
Citation: Wu P, Yao T, Ren Y, Ye J, Qing Y, Li Q and Gui M (2021) Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes. Front. Microbiol. 12:583129. doi: 10.3389/fmicb.2021.583129
Received: 14 July 2020; Accepted: 16 June 2021;
Published: 05 July 2021.
Edited by:
Daniela Minerdi, University of Turin, ItalyReviewed by:
Georg Hausner, University of Manitoba, CanadaZhengfei Wang, Yancheng Teachers University, China
Huei-Mien Ke, Biodiversity Research Center, Academia Sinica, Taiwan
Copyright © 2021 Wu, Yao, Ren, Ye, Qing, Li and Gui. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yuan Qing, 64899867@qq.com; Qiang Li, leeq110@126.com; Mingying Gui, Guimingying@126.com