The apicoplast (ap) is a unique, non-photosynthetic organelle found in most apicomplexan parasites. Due to the essential roles that this organelle has, it has been widely considered as target for drugs against diseases caused by apicomplexans. Exploring the ap genomes of such parasites would provide a better understanding of their systematics and their basic molecular biology for therapeutics. However, there is limited information available on the ap genomes of apicomplexan parasites. In the present study, the ap genomes of two operational taxonomic units of Babesia (known as Babesia sp. Lintan [Bl] and Babesia sp. Xinjiang [Bx]) from sheep were sequenced, assembled and annotated using a massive parallel sequencing based approach. Then, the gene content and gene order in these ap genomes (similar to 30.7 kb in size) were defined and compared, and the genetic differences were assessed. In addition, a phylogenetic analysis of ap genomic data sets was carried out to assess the relationships of these taxonomic units with other apicomplexan parasites for which complete ap genomic data sets were publicly available. The results showed that the ap genomes of Bl and Bx encode 59 and 57 genes, respectively, including 2 ribosomal RNA genes, 25 transfer RNA genes and 30-32 protein-encoding genes, being similar in content to those of Babesia bovis and B. orientalis. Ap gene regions that might serve as markers for future epidemiological and population genetic studies of Babesia species were identified. Using sequence data for a subset of six protein-encoding genes, a close relationship of Bl and Bx with Babesia bovis from cattle and B. orientalis from water buffalo was inferred. Although the focus of the present study was on Babesia, we propose that the present sequencing-bioinformatic approach should be applicable to organellar genomes of a wide range of apicomplexans of veterinary importance. (C) 2016 Elsevier B.V. All rights reserved.