Pleistocene climatic cycling and diversification of the Andean treefrog, Hypsiboas andinus

Our understanding of the causes of diversification of Neotropical organisms lags behind that of Northern Hemisphere biota, especially for montane and temperate regions of southern South America. We investigated the mitochondrial DNA genealogical patterns in 262 individuals of the frog Hypsiboas andinus from 26 sites across the eastern ranges of the Andes Mountains in Argentina and Bolivia. Our phylogenetic analyses indicate at least three distinct lineages: one representing H. andinus from Northwestern Argentina and southern Bolivia, at least one H. andinus lineage from northern Bolivia, and one clade containing both H. andinus (from the southern portion of the species range) and its putative sister taxon Hypsiboas riojanus. Hypsiboas andinus samples from northern Bolivia are well differentiated and may represent distinct species. The northern Argentine H. andinus lineage and southern H. andinus/H. riojanus lineage likely diverged between 2 and 6 million years ago; their current sympatry may be the result of secondary contact due to range expansion after isolation during Andean uplift or may reflect cryptic species. Within the geographically extensive northern H. andinus clade, we found significant geographical structuring consistent with historical fragmentation and subsequent range expansion. The timing of this fragmentation and range expansion coincide with the Pleistocene, a time of extensive climatic cycling and vegetational shifts. Average divergence among clades is lower than those found for other Neotropical taxa, highlighting the potential importance of recent climatic history in diversification in the southern Andes.     

The last diadectomorph sheds light on Late Palaeozoic tetrapod biogeography

Diadectomorpha is a clade of Late Palaeozoic vertebrates widely recognized as the sister group of crown-group Amniota and the first tetrapod lineage to evolve high-fibre herbivory. Despite their evolutionary importance, diadectomorphs are restricted stratigraphically and geographically, with all records being from the Upper Carboniferous and Lower Permian of North America and Germany. We describe a new diadectomorph, Alveusdectes fenestralis, based on a partial skull from the Upper Permian of China. The new species exhibits the derived mechanism for herbivory and is recovered phylogenetically as a deeply nested diadectid. Approximately 16 Myr younger than any other diadectomorph, Alveusdectes is the product of at least a 46 Myr ghost lineage. How much of this time was probably spent in Russia and/or central Asia will remain unclear until a specimen is described that subdivides this cryptic history, but the lineage assuredly crossed this region before entering the relatively isolated continent of North China. The discovery of Alveusdectes raises important questions regarding diadectomorph extinction dynamics including what, if any, ecological factors limited the diversity of this group in eastern Pangea. It also suggests that increased sampling in Asia will likely significantly affect our views of clade and faunal insularity leading up to the Permo-Triassic extinction.     

Cladistic relationships of extant and fossil hominoids

There is general agreement that the hominoid primates form a monophyletic group, that the extant great apes and humans form a second clade within that group with the gibbons as the sister group, and that the African apes and humans form a third clade. Although it has recently been proposed that humans and orang utans are sister taxa and also that the great apes form a clade to the exclusion of humans, our analysis, particularly of the molecular evidence, supports the existence of an African ape and human clade. The major problem in hominoid phylogeny at present is the relationships of the species within this clade: morphological data generally support the existence of an African ape clade which is the sister group to humans; some molecular data also support this conclusion, but most molecular evidence indicates the existence of a chimpanzee/human clade. We have cladistically re-analysed the DNA and protein sequence data for which apomorphic character states can be assessed. It is clear that there is a high degree of homoplasy whichever branching pattern is produced, with some characters supporting the existence of a chimpanzee/human clade and others supporting an African ape clade. When the cladistic analyses of morphological and molecular data are combined we believe that the most parsimonious interpretation of the data is that the African apes form a clade which is the sister taxon of the human (i.e., Australopithecus, Homo and Paranthropus) clade.This paper is not intended as a survey of all hominoid fossils but as a study of branching points in hominoid evolution and fossils are included which are relevant to this branching pattern. The analysis of fossil taxa in this study leads us to conclude that Proconsul is the sister taxon to the later Hominoidea. A number of middle Miocene forms such as Dryopithecus, Kenyapithecus, Heliopithecus and Afropithecus are shown to share derived characters with great apes and humans and provide evidence for the divergence of that clade from the gibbon lineage prior to 18 Ma. The position that Sivapithecus represents the sister group of the orang utan clade is supported here and shows that the orang utan lineage had diverged from the African ape and human lineage prior to 11·5 Ma. There is unfortunately no definitive fossil cvidence on branching sequences within the African ape and human clade, although a new specimen from Samburu, Kenya may be related to the gorilla.     

Discovery of a relict lineage and monotypic family of passerine birds

Analysis of one of the most comprehensive datasets to date of the largest passerine bird clade, Passerida, identified 10 primary well-supported lineages corresponding to Sylvioidea, Muscicapoidea, Certhioidea, Passeroidea, the ‘bombycillids’ (here proposed to be recognized as Bombycilloidea), Paridae/Remizidae (proposed to be recognized as Paroidea), Stenostiridae, Hyliotidae, Regulidae (proposed to be recognized as Reguloidea) and spotted wren-babbler Spelaeornis formosus. The latter was found on a single branch in a strongly supported clade with Muscicapoidea, Certhioidea and Bombycilloidea, although the relationships among these were unresolved. We conclude that the spotted wren-babbler represents a relict basal lineage within Passerida with no close extant relatives, and we support the already used name Elachura formosa and propose the new family name Elachuridae for this single species.     

Phylogenetics of the Antopocerus-Modified Tarsus Clade of Hawaiian Drosophila: Diversification across the Hawaiian Islands

The Hawaiian Drosophilidae radiation is an ecologically and morphologically diverse clade of almost 700 described species. A phylogenetic approach is key to understanding the evolutionary forces that have given rise to this diverse lineage. Here we infer the phylogeny for the antopocerus, modified tarsus and ciliated tarsus (AMC) clade, a lineage comprising 16% (91 of 687 species) of the described Hawaiian Drosophilidae. To improve on previous analyses we constructed the largest dataset to date for the AMC, including a matrix of 15 genes for 68 species. Results strongly support most of the morphologically defined species groups as monophyletic. We explore the correlation of increased diversity in biogeography, sexual selection and ecology on the present day diversity seen in this lineage using a combination of dating methods, rearing records, and distributional data. Molecular dating analyses indicate that AMC lineage started diversifying about 4.4 million years ago, culminating in the present day AMC diversity. We do not find evidence that ecological speciation or sexual selection played a part in generating this diversity, but given the limited number of described larval substrates and secondary sexual characters analyzed we can not rule these factors out entirely. An increased rate of diversification in the AMC is found to overlap with the emergence of multiple islands in the current chain of high islands, specifically Oahu and Kauai.     

Cryptic diversity within the invasive virile crayfish <Emphasis Type="Italic">Orconectes virilis</Emphasis> (Hagen, 1870) species complex: new lineages recorded in both native and introduced ranges

The virile crayfish (Orconectes virilis) represents a cryptic species complex with several lineages known in the USA, and a wide introduced range. In Europe, O. virilis is an emerging invader, established during the last decade in at least two areas—one in the United Kingdom and another in the Netherlands. We assessed the position of both known European populations within the species complex by sequencing part of the mitochondrial gene for cytochrome c oxidase subunit I. Tested UK and Dutch individuals did not belong to any mitochondrial lineage recorded in North America so far but formed a separate clade, the original distribution area of which is unknown. Additionally sequenced virile crayfish from Iowa (USA) also represented a new clade, suggesting that undiscovered lineage variation within O. virilis remains high. This exemplifies that genetic analyses of invading populations may provide new insights into diversity of a taxon in its original range.     

Virus-Like Attachment Sites and Plastic CpG Islands: Landmarks of Diversity in Plant Del Retrotransposons

Full-length Del elements from ten angiosperm genomes, 5 monocot and 5 dicot, were retrieved and putative attachment (att) sites were identified. In the 2432 Del elements, two types of U5 att sites and a single conserved type of U3 att site were identified. Retroviral att sites confer specificity to the integration process, different att sites types therefore implies lineage specificity. While some features are common to all Del elements, CpG island patterns within the LTRs were particular to lineage specific clusters. All eudicot copies grouped into one single clade while the monocots harbour a more diverse collection of elements. Furthermore, full-length Del elements and truncated copies were unevenly distributed amongst chromosomes. Elements of Del lineage are organized in plants into three clusters and each cluster is composed of elements with distinct LTR features. Our results suggest that the Del lineage efficiently amplified in the monocots and that one branch is probably a newly emerging sub-lineage. Finally, sequences in all groups are under purifying selection. These results show the LTR region is dynamic and important in the evolution of LTR-retrotransposons, we speculate that it is a trigger for retrotransposon diversification.     

Chlorophyll-binding proteins revisited - a multigenic family of light-harvesting and stress proteins from a brown algal perspective

Background

Mitochondrial introgression may result in the mitochondrial genome of one species being replaced by that of another species without leaving any trace of past hybridization in its nuclear genome. Such introgression can confuse the species genealogy estimates and lead to absurd inferences of species history. We used a phylogenetic approach to explore the potential mitochondrial genome introgression event(s) between two closely related green pond frog species, Pelophylax nigromaculatus and P. plancyi.

Results

DNA sequence data of one mitochondrial and two nuclear genes from an extensive sampling of the two species were collected, and the genealogies of the three genes were constructed and compared. While the two nuclear genes congruently showed mutual reciprocal monophyly of both species, the mitochondrial phylogeny separated a Korean P. nigromaculatus clade, a paraphyletic central China P. plancyi assemblage, and a large well-supported introgression clade. Within the introgression clade, the mitochondrial haplotypes of the two species were mixed together. This reticulated pattern can be most parsimoniously explained by an ancient mitochondrial introgression event from P. plancyi to P. nigromaculatus that occurred at least 1.36 MYA, followed by multiple recent introgression events from P. nigromaculatus back to P. plancyi within the last 0.63 MY. The re-constitution of previously co-adapted genomes in P. plancyi may be responsible for the recent rampant introgression events. The Korean P. nigromaculatus clade likely represents the only surviving "true" mitochondrial lineage of P. nigromaculatus, and the central China P. plancyi assemblage likely represents the "original" P. plancyi mitochondrial lineage. Refugia in the Korean Peninsula and central China may have played a significant role in preserving these ancient lineages.

Conclusions

The majority of individuals in the two species have either introgressed (P. nigromaculatus) or reclaimed (P. plancyi) mitochondrial genomes while no trace of past hybridization in their nuclear genomes was detected. Asymmetrical reproductive ability of hybrids and continuous backcrossing are likely responsible for the observed mitochondrial introgression. This case is unique in that it includes an ancient "forward" introgression and many recent "backward" introgressions, which re-constitutes the original nuclear and mitochondrial genomes of P. plancyi. This hybrid system provides an excellent opportunity to study cyto-nuclear interaction and co-adaptation.     

Extrachromosomal, extraordinary and essential—the plasmids of the Roseobacter clade

The alphaproteobacterial Roseobacter clade (Rhodobacterales) is one of the most important global players in carbon and sulfur cycles of marine ecosystems. The remarkable metabolic versatility of this bacterial lineage provides access to diverse habitats and correlates with a multitude of extrachromosomal elements. Four non-homologous replication systems and additional subsets of individual compatibility groups ensure the stable maintenance of up to a dozen replicons representing up to one third of the bacterial genome. This complexity presents the challenge of successful partitioning of all low copy number replicons. Based on the phenomenon of plasmid incompatibility, we developed molecular tools for target-oriented plasmid curing and could generate customized mutants lacking hundreds of genes. This approach allows one to analyze the relevance of specific replicons including so-called chromids that are known as lifestyle determinants of bacteria. Chromids are extrachromosomal elements with a chromosome-like genetic imprint (codon usage, GC content) that are essential for competitive survival in the natural habitat, whereas classical dispensable plasmids exhibit a deviating codon usage and typically contain type IV secretion systems for conjugation. The impact of horizontal plasmid transfer is exemplified by the scattered occurrence of the characteristic aerobic anoxygenic photosynthesis among the Roseobacter clade and the recently reported transfer of the 45-kb photosynthesis gene cluster to extrachromosomal elements. Conjugative transmission may be the crucial driving force for rapid adaptations and hence the ecological prosperousness of this lineage of pink bacteria.     

Global Distribution of Bartonella Infections in Domestic Bovine and Characterization of Bartonella bovis Strains Using Multi-Locus Sequence Typing

Bartonella bovis is commonly detected in cattle. One B. bovis strain was recently isolated from a cow with endocarditis in the USA, suggesting its role as an animal pathogen. In the present study, we investigated bartonella infections in 893 cattle from five countries (Kenya, Thailand, Japan, Georgia, and Guatemala) and 103 water buffaloes from Thailand to compare the prevalence of the infection among different regions and different bovid hosts. We developed a multi-locus sequence typing (MLST) scheme based on nine loci (16S rRNA, gltA, ftsZ, groEL, nuoG, ribC, rpoB, ssrA, and ITS) to compare genetic divergence of B. bovis strains, including 26 representatives from the present study and two previously described reference strains (one from French cows and another from a cow with endocarditis in the USA). Bartonella bacteria were cultured in 6.8% (7/103) of water buffaloes from Thailand; all were B. bovis. The prevalence of bartonella infections in cattle varied tremendously across the investigated regions. In Japan, Kenya, and the Mestia district of Georgia, cattle were free from the infection; in Thailand, Guatemala, and the Dusheti and Marneuli districts of Georgia, cattle were infected with prevalences of 10–90%. The Bartonella isolates from cattle belonged to three species: B. bovis (n=165), B. chomelii (n=9), and B. schoenbuchensis (n=1), with the latter two species found in Georgia only. MLST analysis suggested genetic variations among the 28 analyzed B. bovis strains, which fall into 3 lineages (I, II, and III). Lineages I and II were found in cattle while lineage III was restricted to water buffaloes. The majority of strains (17/28), together with the strain causing endocarditis in a cow in the USA, belonged to lineage I. Further investigations are needed to determine whether B. bovis causes disease in bovids.     

Target enrichment improves phylogenetic resolution in the genus Zanthoxylum (Rutaceae) and indicates both incomplete lineage sorting and hybridization events

Background and Aims Zanthoxylum is the only pantropical genus within Rutaceae, with a few species native to temperate eastern Asia and North America. Efforts using Sanger sequencing failed to resolve the backbone phylogeny of Zanthoxylum. In this study, we employed target-enrichment high-throughput sequencing to improve resolution. Gene trees were examined for concordance and sectional classifications of Zanthoxylum were evaluated. Off-target reads were investigated to identify putative single-copy markers for bait refinement, and low-copy markers for evidence of putative hybridization events.MethodsA custom bait set targeting 354 genes, with a median of 321 bp, was designed for Zanthoxylum and applied to 44 Zanthoxylum species and one Tetradium species as the outgroup. Illumina reads were processed via the HybPhyloMaker pipeline. Phylogenetic inferences were conducted using coalescent and maximum likelihood methods based on concatenated datasets. Concordance was assessed using quartet sampling. Additional phylogenetic analyses were performed on putative single and low-copy genes extracted from off-target reads.Key ResultsFour major clades are supported within Zanthoxylum: the African clade, the Z. asiaticum clade, the Asian–Pacific–Australian clade and the American–eastern Asian clade. While overall support has improved, regions of conflict are similar to those previously observed. Gene tree discordances indicate a hybridization event in the ancestor of the Hawaiian lineage, and incomplete lineage sorting in the American backbone. Off-target putative single-copy genes largely confirm on-target results, and putative low-copy genes provide additional evidence for hybridization in the Hawaiian lineage. Only two of the five sections of Zanthoxylum are resolved as monophyletic.ConclusionsTarget enrichment is suitable for assessing phylogenetic relationships in Zanthoxylum. Our phylogenetic analyses reveal that current sectional classifications need revision. Quartet tree concordance indicates several instances of reticulate evolution. Off-target reads are proven useful to identify additional phylogenetically informative regions for bait refinement or gene tree based approaches.     

Diversification of Fungal Chitinases and Their Functional Differentiation in Histoplasma capsulatum

Chitinases enzymatically hydrolyze chitin, a highly abundant and utilized polymer of N-acetyl-glucosamine. Fungi are a rich source of chitinases; however, the phylogenetic and functional diversity of fungal chitinases are not well understood. We surveyed fungal chitinases from 373 publicly available genomes, characterized domain architecture, and conducted phylogenetic analyses of the glycoside hydrolase (GH18) domain. This large-scale analysis does not support the previous division of fungal chitinases into three major clades (A, B, C) as chitinases previously assigned to the “C” clade are not resolved as distinct from the “A” clade. Fungal chitinase diversity was partly shaped by horizontal gene transfer, and at least one clade of bacterial origin occurs among chitinases previously assigned to the “B” clade. Furthermore, chitin-binding domains (including the LysM domain) do not define specific clades, but instead are found more broadly across clades of chitinases. To gain insight into biological function diversity, we characterized all eight chitinases (Cts) from the thermally dimorphic fungus, Histoplasma capsulatum: six A clade, one B clade, and one formerly classified C clade chitinases. Expression analyses showed variable induction of chitinase genes in the presence of chitin but preferential expression of CTS3 in the mycelial stage. Activity assays demonstrated that Cts1 (B-I), Cts2 (A-V), Cts3 (A-V), Cts4 (A-V) have endochitinase activities with varying degrees of chitobiosidase function. Cts6 (C-I) has activity consistent with N-acetyl-glucosaminidase exochitinase function and Cts8 (A-II) has chitobiase activity. These results suggest chitinase activity is variable even within subclades and that predictions of functionality require more sophisticated models.     

Multiple causation of phylogeographical pattern as revealed by nested clade analysis of the bamboo viper (Trimeresurus stejnegeri) within Taiwan

In order to assess the utility of nested clade analysis, both standard phylogenetic algorithms and nested clade analysis were performed on a geographically widespread survey of mitochondrial DNA haplotypes of the bamboo viper, Trimeresurus stejnegeri, within Taiwan. Gross tree topologies were congruent for all analyses and indicated the presence of two geographically overlapping clades within Taiwan. The smaller lineage was restricted to the north and east coasts, whereas the larger lineage occupied all but the northern range of the species within Taiwan including the Pacific offshore populations of Green and Orchid Islands. The phylogeographical pattern supports the existence of at least one colonization event from the continent since the initial isolation of Taiwan from the mainland in the Pliocene. However, determining the exact number of colonization events was not possible due to the simultaneous vicariant forces of hypothesized continental landbridge connections and the occurrence of dramatic in situ orogenesis throughout the Pleistocene. Nested clade analysis provided multiple temporal and spatial population historical inferences that are not possible with standard analyses and therefore should become widely applied to future phylogeographical studies.     

Complex Population Structure and Virulence Differences among Serotype 2 Streptococcus suis Strains Belonging to Sequence Type 28

Streptococcus suis is a major swine pathogen and a zoonotic agent. Serotype 2 strains are the most frequently associated with disease. However, not all serotype 2 lineages are considered virulent. Indeed, sequence type (ST) 28 serotype 2 S. suis strains have been described as a homogeneous group of low virulence. However, ST28 strains are often isolated from diseased swine in some countries, and at least four human ST28 cases have been reported. Here, we used whole-genome sequencing and animal infection models to test the hypothesis that the ST28 lineage comprises strains of different genetic backgrounds and different virulence. We used 50 S. suis ST28 strains isolated in Canada, the United States and Japan from diseased pigs, and one ST28 strain from a human case isolated in Thailand. We report a complex population structure among the 51 ST28 strains. Diversity resulted from variable gene content, recombination events and numerous genome-wide polymorphisms not attributable to recombination. Phylogenetic analysis using core genome single-nucleotide polymorphisms revealed four discrete clades with strong geographic structure, and a fifth clade formed by US, Thai and Japanese strains. When tested in experimental animal models, strains from this latter clade were significantly more virulent than a Canadian ST28 reference strain, and a closely related Canadian strain. Our results highlight the limitations of MLST for both phylogenetic analysis and virulence prediction and raise concerns about the possible emergence of ST28 strains in human clinical cases.     

Ancient wolf lineages in India

All previously obtained wolf (Canis lupus) and dog (Canis familiaris) mitochondrial (mt) DNA sequences fall within an intertwined and shallow clade (the 'wolf-dog' clade). We sequenced mtDNA of recent and historical samples from 45 wolves from throughout lowland peninsular India and 23 wolves from the Himalayas and Tibetan Plateau and compared these sequences with all available wolf and dog sequences. All 45 lowland Indian wolves have one of four closely related haplotypes that form a well-supported, divergent sister lineage to the wolf-dog clade. This unique lineage may have been independent for more than 400,000 years. Although seven Himalayan wolves from western and central Kashmir fall within the widespread wolf-dog clade, one from Ladakh in eastern Kashmir, nine from Himachal Pradesh, four from Nepal and two from Tibet form a very different basal clade. This lineage contains five related haplotypes that probably diverged from other canids more than 800,000 years ago, but we find no evidence of current barriers to admixture. Thus, the Indian subcontinent has three divergent, ancient and apparently parapatric mtDNA lineages within the morphologically delineated wolf. No haplotypes of either novel lineage are found within a sample of 37 Indian (or other) dogs. Thus, we find no evidence that these two taxa played a part in the domestication of canids.     

Kretzoiarctos gen. nov., the Oldest Member of the Giant Panda Clade

The phylogenetic position of the giant panda, Ailuropoda melanoleuca (Carnivora: Ursidae: Ailuropodinae), has been one of the most hotly debated topics by mammalian biologists and paleontologists during the last century. Based on molecular data, it is currently recognized as a true ursid, sister-taxon of the remaining extant bears, from which it would have diverged by the Early Miocene. However, from a paleobiogeographic and chronological perspective, the origin of the giant panda lineage has remained elusive due to the scarcity of the available Miocene fossil record. Until recently, the genus Ailurarctos from the Late Miocene of China (ca. 8–7 mya) was recognized as the oldest undoubted member of the Ailuropodinae, suggesting that the panda lineage might have originated from an Ursavus ancestor. The role of the purported ailuropodine Agriarctos, from the Miocene of Europe, in the origins of this clade has been generally dismissed due to the paucity of the available material. Here, we describe a new ailuropodine genus, Kretzoiarctos gen. nov., based on remains from two Middle Miocene (ca. 12–11 Ma) Spanish localities. A cladistic analysis of fossil and extant members of the Ursoidea confirms the inclusion of the new genus into the Ailuropodinae. Moreover, Kretzoiarctos precedes in time the previously-known, Late Miocene members of the giant panda clade from Eurasia (Agriarctos and Ailurarctos). The former can be therefore considered the oldest recorded member of the giant panda lineage, which has significant implications for understanding the origins of this clade from a paleobiogeographic viewpoint.     

Origin and Biology of Simian Immunodeficiency Virus in Wild-Living Western Gorillas

Western lowland gorillas (Gorilla gorilla gorilla) are infected with a simian immunodeficiency virus (SIVgor) that is closely related to chimpanzee and human immunodeficiency viruses (SIVcpz and HIV-1, respectively) in west central Africa. Although existing data suggest a chimpanzee origin for SIVgor, a paucity of available sequences has precluded definitive conclusions. Here, we report the molecular characterization of one partial (BQ664) and three full-length (CP684, CP2135, and CP2139) SIVgor genomes amplified from fecal RNAs of wild-living gorillas at two field sites in Cameroon. Phylogenetic analyses showed that all SIVgor strains clustered together, forming a monophyletic lineage throughout their genomes. Interestingly, the closest relatives of SIVgor were not SIVcpzPtt strains from west central African chimpanzees (Pan troglodytes troglodytes) but human viruses belonging to HIV-1 group O. In trees derived from most genomic regions, SIVgor and HIV-1 group O formed a sister clade to the SIVcpzPtt lineage. However, in a tree derived from 5′ pol sequences (~900 bp), SIVgor and HIV-1 group O fell within the SIVcpzPtt radiation. The latter was due to two SIVcpzPtt strains that contained mosaic pol sequences, pointing to the existence of a divergent SIVcpzPtt lineage that gave rise to SIVgor and HIV-1 group O. Gorillas appear to have acquired this lineage at least 100 to 200 years ago. To examine the biological properties of SIVgor, we synthesized a full-length provirus from fecal consensus sequences. Transfection of the resulting clone (CP2139.287) into 293T cells yielded infectious virus that replicated efficiently in both human and chimpanzee CD4⁺ T cells and used CCR5 as the coreceptor for viral entry. Together, these results provide strong evidence that P. t. troglodytes apes were the source of SIVgor. These same apes may also have spawned the group O epidemic; however, the possibility that gorillas served as an intermediary host cannot be excluded.     

The molecular phylogenetic signature of Bali cattle revealed by maternal and paternal markers

Bali cattle is a domestic cattle breed that can be found in Malaysia. It is a domestic cattle that was purely derived from a domestication event in Banteng (Bos javanicus) around 3,500 BC in Indonesia. This research was conducted to portray the phylogenetic relationships of the Bali cattle with other cattle species in Malaysia based on maternal and paternal lineage. We analyzed the cytochrome c oxidase I (COI) mitochondrial gene and SRY of Y chromosome obtained from five species of the Bos genus (B. javanicus, Bos gaurus, Bos indicus, Bos taurus, and Bos grunniens). The water buffalo (Bubalus bubalis) was used as an outgroup. The phylogenetic relationships were observed by employing several algorithms: Neighbor-Joining (PAUP version 4.0), Maximum parsimony (PAUP version 4.0) and Bayesian inference (MrBayes 3.1). Results from the maternal data showed that the Bali cattle formed a monophyletic clade, and together with the B. gaurus clade formed a wild cattle clade. Results were supported by high bootstrap and posterior probability values together with genetic distance data. For the paternal lineage, the sequence variation is low (with parsimony informative characters: 2/660) resulting an unresolved Neighbor-Joining tree. However, Bali cattle and other domestic cattle appear in two monophyletic clades distinct from yak, gaur and selembu. This study expresses the potential of the COI gene in portraying the phylogenetic relationships between several Bos species which is important for conservation efforts especially in decision making since cattle is highly bred and hybrid breeds are often formed. Genetic conservation for this high quality beef cattle breed is important by maintaining its genetic characters to prevent extinction or even decreased the genetic quality.