Drivers of interlineage variability in mitogenomic evolutionary rates in Platyhelminthes

Studies of forces driving interlineage variability in the evolutionary rates (both sequence and architecture) of mitochondrial genomes often produce contradictory results. Flatworms (Platyhelminthes) exhibit the fastest-evolving mitogenomic sequences among all bilaterian phyla. To test the effects of multiple factors previously associated with different aspects of mitogenomic evolution, we used mitogenomes of 223 flatworm species, phylogenetic multilevel regression models, and causal inference. Thermic host environment (endothermic vs. ectothermic) had nonsignificant impacts on both sequence evolution and mitogenomic size. Mitogenomic gene order rearrangements (GORR) were mostly positively correlated with mitogenomic size (R2 ≈ 20–30%). Longevity was not (negatively) correlated with sequence evolution in flatworms. The predominantly free-living “turbellaria” exhibited much shorter branches and faster-evolving mitogenomic architecture than parasitic Neodermata. As a result, “parasitism” had a strong explanatory power on the branch length variability (>90%), and there was a negative correlation between GORR and branch length. However, the stem branch of Neodermata comprised 63.6% of the total average branch length. This evolutionary period was also marked by a high rate of gene order rearrangements in the ancestral Neodermata. We discuss how this period of rapid evolution deep in the evolutionary history may have decoupled sequence evolution rates from longevity and GORR, and overestimated the explanatory power of “parasitism”. This study shows that impacts of variables often vary across lineages, and stresses the importance accounting for the episodic nature of evolutionary patterns in studies of mitogenomic evolution.