Comparative Analysis of Boletales Mitogenomes

Base compositions of the two Coniophora mitogenomes and other Boletales mitogenomes were calculated using the DNASTAR Lasergene v7.1 (http://www.dnastar.com/). Strand asymmetries of the 6 Boletales mitogenomes were assessed based on the following formulas: AT skew = [A - T] / [A + T], and GC skew = [G - C] / [G + C] [37] (link). The 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 6 Boletales mitogenomes were calculated using MEGA v6.06 [38] based on the Kimura-2-parameter (K2P) substitution model. The DnaSP v6.10.01 software [39] (link) was used to calculate synonymous (Ks) and nonsynonymous substitution rates (Ka) for core PCGs in the 6 Boletales mitogenomes. Gene collinearity analysis of 6 Boletales species was conducted by using Mauve v2.4.0 [40] (link).