The Mitochondrial Genome of Melipona fasciculata (Apidae, Meliponini): Genome Organization and Comparative Analyses, Phylogenetic Implications and Divergence Time Estimations

The native stingless bee Melipona fasciculata is economically and ecologically important to the Brazilian Northeast, providing a sustainable source of income to family farmers and being considered an effective pollinator in most ecosystems and crops. This study describes, for the first time, the mitogenome of the species and its phylogenetic position. The mitochondrial genome was sequenced using a MiSeq Sequencer (Illumina Inc.) and compared with other GenBank bee mitogenomes. The length of the mitochondrial DNA, excluding most of the control region, is 14,753 bp, and contains 13 protein-coding genes (PCGs), 21 transfer RNAs, 2 ribosomal RNAs (16S and 12S), and 1 AT-rich region. The GC-content of the M. fasciculata mitogenome was 13.4%. Of the 36 coding regions, 12 tRNAs and 9 PCGs were encoded in the heavy strand, and 9 tRNAs, 4 PCGs and 2 rRNAs were encoded in the light strand. The relative orientation and gene order was the same as other stingless bee mitogenomes. Phylogenetic inference produced well-resolved relationships with high statistical support for concordant branch topologies, under different optimization schemes and model parameters, within and among Melipona, Bombus, Apis, and related clades of Hymenoptera. In general, our divergence time estimates, which were based on the concatenated gene sequences (PCGs + rRNAs) from various groups, overlapped estimations captured by Bayesian analysis from different studies. The divergence time among Melipona species was estimated to occur during the Oligocene, approximately 24 Mya (95% HPD 14-36 Mya). Our results represent a valuable addition to help understanding not only the taxonomy and evolution of Brazilian stingless bee species, but also to uncover historical dispersal and isolation patterns in Meliponinae.