Phylogenetic relationship among all living species of the genus Bubalus based on DNA sequences of the cytochrome b gene

The cytochrome b genes of all living species of Bubalus, including the river type and the swamp type of domestic buffaloes (Bubalus bubalis), were sequenced to clarify their phylogenetic relationships. These sequences were compared together with the African buffalo (Syncerus caffer) and banteng (Bos javanicus) sequences as an outgroup. Phylogenetic trees of Bubalus species based on the DNA sequences of the cytochrome b gene demonstrated that the tamarraw (Bubalus mindorensis), endemic to the Philippines, could be classified into the subgenus Bubalus, not the subgenus Anoa. The divergence time between the lowland anoa (B. depressicornis) and the mountain anoa (B. quarlesi) was estimated at approximately 2.0 million yeats (Myr), which is almost the same as the coalescence time for the Bubalus sequences. This large genetic distance supports the idea that the lowland anoa and the mountain anoa are different species. An unexpectedly large genetic distance between the river and the swamp type of domestic buffaloes suggests a divergence time of about 1.7 Myr, while the swamp type was noticed to have the closest relationship with the tamarraw (1.5 Myr). This result implies that the two types of domestic buffaloes have differentiated at the full species level.