Comparative Mitogenomic Analysis Reveals Gene and Intron Dynamics in Rubiaceae and Intra-Specific Diversification in Damnacanthus indicus

The dynamic evolution of mitochondrial gene and intron content has been reported across the angiosperms. However, a reference mitochondrial genome (mitogenome) is not available in Rubiaceae. The phylogenetic utility of mitogenome data at a species level is rarely assessed. Here, we assembled mitogenomes of six Damnacanthus indicus (Rubiaceae, Rubioideae) representing two varieties (var. indicus and var. microphyllus). The gene and intron content of D. indicus was compared with mitogenomes from representative angiosperm species and mitochondrial contigs from the other Rubiaceae species. Mitogenome structural rearrangement and sequence divergence in D. indicus were analyzed in six individuals. The size of the mitogenome in D. indicus varied from 417,661 to 419,435 bp. Comparing the number of intact mitochondrial protein-coding genes in other Gentianales taxa (38), D. indicus included 32 genes representing several losses. The intron analysis revealed a shift from cis to trans splicing of a nad1 intron (nad1i728) in D. indicus and it is a shared character with the other four Rubioideae taxa. Two distinct mitogenome structures (type A and B) were identified. Two-step direct repeat-mediated recombination was proposed to explain structural changes between type A and B mitogenomes. The five individuals from two varieties in D. indicus diverged well in the whole mitogenome-level comparison with one exception. Collectively, our study elucidated the mitogenome evolution in Rubiaceae along with D. indicus and showed the reliable phylogenetic utility of the whole mitogenome data at a species-level evolution.