Resolving a zoological mystery: the kouprey is a real species

The kouprey is a rare and enigmatic forest ox discovered by scientists in Cambodia only in 1937. Numerous morphological hypotheses have been proposed for the origin of the kouprey: that it is a species closely related to banteng and gaur, two other wild oxen of southeast Asia; a morphologically divergent species placed in a separate genus, named Novibos; a wild species linked to aurochs and domestic cattle; a vicariant population of banteng; a feral cattle; or a hybrid of banteng with either zebu cattle, gaur or water buffalo. In a recent paper, which gained a lot of media coverage, Galbreath et al. analysed mitochondrial DNA sequences and concluded that the kouprey never existed as a wild, natural species, and that it was a feral hybrid between banteng and zebu cattle. Here we analyse eight DNA markers-three mitochondrial regions and five nuclear fragments-representing an alignment of 4582 nucleotides for the holotype of the kouprey and all related species. Our results demonstrate that the kouprey is a real and naturally occurring species, and show that Cambodian populations of banteng acquired a mitochondrial genome of kouprey by natural introgressive hybridization during the Pleistocene epoch.     

Maternal and paternal lineages in cross-breeding bovine species. Has wisent a hybrid origin?

The tribe Bovini comprises cattle and cattle-like species. Reconstructions of their phylogeny have so far been incomplete and have yielded conflicting conclusions about the relationship of American bison and wisent (European bison). We have compared the sequences of three mitochondrial and two Y-chromosomal DNA segments. Mitochondrial DNA indicates that four distinct maternal lineages diverged after an early split-off of the buffalo species, leading to (1) taurine cattle and zebu, (2) wisent, (3) American bison and yak, and (4) banteng, gaur, and gayal, respectively. At a higher level, lineages (1) and (2) and lineages (3) and (4) are probably associated. In contrast, Y-chromosomal sequences indicate a close association of American and European bison, which is in agreement with their morphological similarity, complete fertility of hybrid offspring, and amplified fragment length polymorphism (AFLP) fingerprints of nuclear DNA. One explanation for the anomalous divergence of the mitochondrial DNA from the two bison species is lineage sorting, which implies that two distinct mitochondrial lineages coexisted in the bison-yak branch until the recent divergence of American bison and wisent. Alternatively, the wisent may have emerged by species hybridization initiated by introgression of bison bulls in another ancestral species. This "transpatric" mode of species formation would be consistent with the recent appearance of the wisent in the fossil record without clearly identifiable ancestors.     

What is the taxonomic status of the Cambodian banteng and does it have close genetic links with the kouprey?

In a recent paper, Galbreath, Mordacq & Weiler (2006) concluded that the kouprey Bos sauveli was not a natural species, but rather a feral animal derived from hybridization between banteng Bos javanicus and zebu Bos taurus indicus . Here, we analyze two mitochondrial genes (cytochrome b and subunit II of the cytochrome c oxidase) for all the seven species of the subtribe Bovina, including new sequences for several specimens of banteng, zebu and gaur of Cambodia. Our analyses indicate that mitochondrial sequences of Cambodian banteng are divergent from those of Javan banteng (mean difference =4.27%), but similar to those of kouprey (1.25%). We propose two conflicting hypotheses to interpret these results: (1) the Cambodian and Indonesian banteng belong to two distinct species, and the kouprey derived from Cambodian banteng; (2) all subspecies of banteng belong to Bo. javanicus , but the mitochondrial genome of kouprey was transferred by natural hybridization into the ancestor of Cambodian banteng. Morphological, ecological and ethological characteristics of banteng and kouprey rather support the second hypothesis. However, we need to sequence nuclear markers, and to analyze banteng from Lao, Myanmar, Thailand and Vietnam, to give a definitive conclusion on the taxonomic status of banteng and kouprey.     

Molecular phylogeny of the tribe Bovini (Bovidae, Bovinae) and the taxonomic status of the Kouprey, Bos sauveli Urbain 1937

The kouprey is a very rare bovid species of the Indochinese peninsula, and no living specimen has been described for a long time, suggesting that it is possibly extinct. Its systematic position within the tribe Bovini remains confused since the analyses of morphological characters have led to several conflicting hypotheses. Some authors have also suggested that it could be a hybrid species produced by the crossing of the banteng with gaur, zebu, or water buffalo. Here we performed a molecular phylogeny of the tribe Bovini to determine the taxonomic status of the kouprey. DNA was extracted from the holotype specimen preserved in the MNHN collections. Phylogenetic analyses were carried out on a matrix including all the taxonomic diversity described in the tribe Bovini, and 2065 nucleotide characters, representing three different markers, i.e., the promotor of the lactoferrin and two mitochondrial genes (cytochrome b and subunit II of the cytochrome c oxidase). The results show that the kouprey belongs to the subtribe Bovina, and that three different clades can be evidenced into this group: the first includes the domestic ox, zebu, and European bison; the second incorporates the yak and American bison; and the third contains the kouprey, banteng and gaur. All hypotheses involving hybridization for the origin of the kouprey can be rejected, confirming that it is a real wild species. Molecular datings and biogeographic inferences suggest that the kouprey diverged from banteng and gaur during the Plio-Pleistocene of Asia. In addition, several molecular signatures were detected in the cytochrome b gene, permitting a molecular identification of the kouprey. We propose a conservation project based on a molecular taxonomy approach for tracking the kouprey in Indochina in order to determine whether some populations still survive in the wild.     

Paternally inherited markers in bovine hybrid populations

The genetic integrity of crossfertile bovine- or cattle-like species may be endangered by species hybridization. Previously, amplified fragment length polymorphism, satellite fragment length polymorphism and microsatellite assays have been used to analyze the species composition of nuclear DNA in taurine cattle, zebu, banteng and bison populations, while mitochondrial DNA reveals the origin of the maternal lineages. Here, we describe species-specific markers of the paternally transmitted Y-chromosome for the direct detection of male-mediated introgression. Convenient PCR-restriction fragment length polymorphism and competitive PCR assays are shown to differentiate the Y-chromosomes of taurine cattle, American bison and European bison, and to detect the banteng origin of Indonesian Madura and Bali cattle bulls.     

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species – taurine cattle, zebu and yak – and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.     

On the Origin of Indonesian Cattle

Background

Two bovine species contribute to the Indonesian livestock, zebu (Bos indicus) and banteng (Bos javanicus), respectively. Although male hybrid offspring of these species is not fertile, Indonesian cattle breeds are supposed to be of mixed species origin. However, this has not been documented and is so far only supported by preliminary molecular analysis.

Methods and Findings

Analysis of mitochondrial, Y-chromosomal and microsatellite DNA showed a banteng introgression of 10–16% in Indonesian zebu breeds. East-Javanese Madura and Galekan cattle have higher levels of autosomal banteng introgression (20–30%) and combine a zebu paternal lineage with a predominant (Madura) or even complete (Galekan) maternal banteng origin. Two Madura bulls carried taurine Y-chromosomal haplotypes, presumably of French Limousin origin. In contrast, we did not find evidence for zebu introgression in five populations of the Bali cattle, a domestic form of the banteng.

Conclusions

Because of their unique species composition Indonesian cattle represent a valuable genetic resource, which potentially may also be exploited in other tropical regions.     

Levels of conservation and variation of heterochromatin and nucleolus organizers in the Bovidae

Chromomycin A3 banding of the mitotic sets of 10 species of Bovidae (cattle, wisent, yak, banteng, gaur, red buffalo, swamp buffalo, sheep, mufflon, and goat) serves to demarcate both centromeric constitutive heterochromatin and R-banding patterns capable of identifying all the chromosomes within a given complement. In all species significant amounts of chromomycin-bright heterochromatin are present at the centromeres of all autosomes, though there was a high degree of intra- and inter-individual variation in the size of the heterochromatic blocks. Marked interspecies differences in the centromeric patterns were evident. The X chromosomes contained appreciable amounts of centromeric heterochromatin only in the two buffaloes. All the animals studied lacked distamycin A - diamidinophenylindole type heterochromatin. AgNO3 staining was applied sequentially to detect the location of active nucleolus organizer regions (NORs). The distribution of NORs was reasonably conservative in most of the species. An exceptional situation was found in the two buffaloes, where only one NOR pair matched with the standard karyotype of the Bovidae.     

Diversity and phylogeny of mitochondrial DNA isolated from mithun Bos frontalis located in Bhutan

We sequenced the 16S rRNA gene in mitochondrial DNA to characterize mithun located in Bhutan and to increase our understanding of its origin. We compared mithun with yak, European cattle, Bhutanese zebu and Indian zebu. Sequencing revealed low nucleotide diversity within the mithun population and their phylogenetic proximity to gaur. A close relationship between Bhutanese mithun and gaur was confirmed by an additional comparison with wild gaur specimens from three locations in Bhutan. Direct domestication of mithun from gaur was supported, while maternal contribution from the cattle lineage during domestication was not supported.     

Genetically solving a zoological mystery: was the kouprey (Bos sauveli) a feral hybrid?

A famous zoological discovery of the 20th century was that of the kouprey Bos sauveli , a medium-sized ox inhabiting Cambodian forests. The kouprey was suspiciously intermediate between banteng oxen and domestic zebu cattle in its structure. Mitochondrial DNA sequences of mainland banteng are compared here with a published kouprey sequence, and the comparison demonstrates a close relationship. Either the kouprey derives partly from banteng or (less likely) these particular banteng acquired kouprey DNA via recent genetic introgression. The kouprey may have been a feral hybrid form, a descendant of domestic oxen, rather than a natural species.     

A bubaline-derived satellite DNA probe uncovers generic affinities of gaur with other bovids

DNA typing using genome derived cloned probes may be conducted for ascertaining genetic affinities of closely related species. We analysed gaurBos gaurus, cattleBos indicus, buffaloBubalus bubalis, sheepOvis aries and goatCapra hircus DNA using buffalo derived cloned probe pDS5 carrying an array ofBamHI satellite fraction of 1378 base residues to uncover its genomic organization. Zoo-blot analysis showed that pDS5 does not cross hybridize with non-bovid animals and surprisingly with female gaur genomic DNA. The presence of pDS5 sequences in the gaur males suggests their possible location on the Y chromosome. Genotyping of pDS5 withBamHI enzyme detected mostly monomorphic bands in the bubaline samples and polymorphic ones in cattle and gaur giving rise to clad specific pattern. Similar typing withRsaI enzyme also revealed clad specific band pattern detecting more number of bands in buffalo and fewer in sheep, goat and gaur samples. Copy number variation was found to be prominent in cattle and gaur withRsaI typing. Our data based on matched band profiles (MBP) suggest that gaur is genetically closer to cattle than buffalo contradicting the age-old notion held by some that gaur is a wild buffalo. The pDS5 clone has a potential for estimating the generic and genetic relationship amongst closely related bovid species.     

Weights of gaur (Bos gaurus) and banteng (Bos javanicus) killed by tigers in Thailand

The primary prey of tigers across much of South‐East Asia has been depleted, reducing the ability of already limited habitat to support tigers. To better understand the extent to which two of the largest prey species, gaur (Bos gaurus) and banteng (Bos javanicus), contribute to the tiger's diet, we estimated the average size of these species killed by tigers. This information is needed to more accurately calculate biomass of these species in the tiger's diet and to devise strategies to increase tiger carrying capacity where habitat is fragmented and limited in west‐central Thailand. We used temporally clumped locations of 24 satellite radio‐collared tigers to identify their kill sites and obtained mandibles from 82 gaur and 79 banteng. Kills were aged by teeth eruption sequence, sectioning the M1 molar and counting cementum annuli. Of all gaur killed, 45.2% were adults; of all banteng killed, 55.7% were adults. The average weight of banteng killed was 423.9 kg, which was similar to the 397.9 kg average weight for gaur. The mean weight of both prey species is 3.5–4.5 times greater than the predicted 1:1 preferred prey to predator ratio. In the absence of medium‐sized prey, killing these larger animals may be especially critical for female tigers provisioning nearly independent young when male offspring are already larger than the mother. This is the first study to present data on the average weights of gaur and banteng killed in South‐East Asia, and these results suggest that these are key prey species to target in tiger prey recovery efforts.     

DNA fingerprinting of individual species and intergeneric and interspecific hybrids of genera Bos and Bison,subfamily Bovinae

Genome fingerprinting with a hypervariable minisatellite sequence of phage M13 DNA was used to study the genetic variation in individual species of the genera Bos and Bison (subfamily Bovinae) and in their interspecific and intergeneric hybrids. DNA fingerprints were obtained for domestic cow Bos taurus primigenius, vatussy Bos taurus macroceros, banteng Bos javanicus, gaur Bos gaurus, wisent Bison bonasus, bison Bison bison, and for the interspecific and intergeneric hybrids. Compared with the original species, most hybrids showed a greater variation in number and size of hybridization fragments. An association was revealed between the number of hybridization fragments and blood composition of interspecific hybrids resulting from unique crossing of domestic cow and banteng. Pairwise similarity coefficients were calculated to construct a dendrogram of genetic similarity, which reflected the relationships between the parental species and hybrids varying in blood composition. The applicability of the method for identifying interspecific and intergeneric hybrids and for studying the consequences of distant hybridization in the subfamily Bovinae is discussed.     

牛科动物HSL基因序列分析及其分子进化研究

在对牛科中4种动物即牦牛、瘤牛、普通牛和水牛HSL基因外显子Ⅰ部分核苷酸序列进行测定的基础上,与Gen-Bank中其他物种相应基因核苷酸序列、氨基酸序列进行了比对分析,并构建了牦牛与其他物种间分子系统进化树。结果表明:牦牛与普通牛、瘤牛、水牛、猪、人、小鼠、大鼠7个物种HSL基因外显子Ⅰ部分核苷酸序列间保守性较高,同源性大小依次为99.8%、99.6%、97.4%、90.6%、88.4%、83.5%、82.3%。相应氨基酸序列间保守性更高,同源性分别为100%、100%、98.2%、94.0%、92.2%、89.8%、89.8%。牦牛与各物种该基因部分核苷酸序列间碱基变异类型主要表现为碱基转换和颠换,无碱基插入和缺失发生,碱基转换的频率高于颠换的频率;在核苷酸水平上的多数碱基替换都是同义替换;序列间单碱基变异位点大多出现在同一位点,多发生在密码子第3位,其次是第1位,最少发生在第2位,符合分子进化的中性学说。HSL基因外显子Ⅰ部分核苷酸序列进行多序列对位排列构建的各物种间分子系统进化树结果表明,普通牛和瘤牛首先聚为一类,再分别与牦牛、水牛、猪、人聚类,最后与大鼠、小鼠聚为一类。该聚类结果与动物学上的分类结果一致,表明HSL基因外显子Ⅰ部分核苷酸序列适合于构建物种间分子系统进化树。研究表明,牦牛、普通牛和瘤牛3个物种间的遗传距离大小相近,牦牛和水牛间的遗传距离与普通牛、瘤牛和水牛间的遗传距离大小相当。牦牛、普通牛和瘤牛3个物种间的遗传距离远小于它们各自与水牛这一物种的遗传距离,它们三者之间的亲缘关系也相对于它们各自与水牛间的亲缘关系都较近,故将牦牛、普通牛和瘤牛划分在同一个属——牛属(Bos)更为合理。     

Phylogenetical and ontogenetical studies on the molecular weight heterogeneity of bovine serum transferrin

Antitransferrin (Tf) rabbit serum was highly specific: it reacted with Tfs of ruminants, such as European breeds and Zebu breeds of cattle, Bali cattle, banteng, swamp and river types of water buffalo, anoa, goat, sheep, deer, antelope, camel, and giraffe, but did not react with serum of other non-ruminant species, such as pig, wild boar, hippopotamus, horse, rabbit, rat, chicken, etc. Electrophoresis of Tf and immunoglobulin G (IgG) complexes was carried out using sodium dodecyl sulfate--polyacrylamide gel electrophoresis (SDS-PAGE). Within ruminants, the following species showed two Tf molecules on SDS-PAGE; European and Zebu cattle, Bali cattle, banteng, two types of water buffalo, and two species of anoa. Other ruminants, sheep, goat, deer, antelope, camel, and giraffe, etc., showed only one Tf molecule. The Tf heterogeneity in molecular weight was, thus, restricted to Bos, Bubalus, and Anoa. The molecular weight of Tf of water buffalo was slightly larger than that of cattle on the gel. The peptide pattern from cyanogen bromide cleavage of Tf of the water buffalo differed clearly from that of cattle. Fetal Tf showed only one molecule during development, but a newborn calf has two Tf molecules, (one large and one small) within 18 hr after birth. We suggest, therefore, that the small molecules formed during the last month of gestation. The peptide patterns of adult and fetal Tfs cleaved by cyanogen bromide differed with regard to the two large peptides; fetal Tf, lacking the second-largest peptide, had twice the amount of the largest peptide compared with adult Tf. From these results, we suggest that a change in peptide sequence occurs from the last month of gestation, when the largest peptide is degraded to the second largest. However, a Tf-like protein detected in the liver microsomal fraction has only one molecular size, both in adult and in fetal livers.     

DNA Fingerprinting of Individual Species and Intergeneric and Interspecific Hybrids of the Genera Bos and Bison, Subfamily Bovinae

Genome fingerprinting with a hypervariable minisatellite sequence of phage M13 DNA was used to study the genetic variation in individual species of the genera Bosand Bison(subfamily Bovinae) and in their interspecific and intergeneric hybrids. DNA fingerprints were obtained for domestic cow Bos taurus primigenius, vatussy Bos taurus macroceros, banteng Bos javanicus, gaur Bos gaurus, wisent Bison bonasus, bison Bison bison, and for the interspecific and intergeneric hybrids. Compared with the original species, most hybrids showed a greater variation in number and size of hybridization fragments. An association was revealed between the number of hybridization fragments and degree of consanguinity of interspecific hybrids resulting from unique crossing of domestic cow and banteng. Pairwise similarity coefficients were calculated to construct a dendrogram of genetic similarity, which reflected the relationships between the parental species and hybrids varying in degree of consanguinity. The applicability of the method for identifying interspecific and intergeneric hybrids and for studying the consequences of hybridization in the subfamily Bovinae is discussed.     

Y‐chromosomal variation confirms independent domestications of swamp and river buffalo

Y‐chromosomal variation in the water buffalo was analysed by sequencing of DBY, ZFY and SRY gene segments. A clear separation of the paternal lineages of the river and swamp types parallels the differences between their maternal lineages and nuclear DNA. Sequence divergence was found to be comparable to the divergence of taurine cattle and zebu, and this divergence predated domestication, confirming that river and swamp buffalo originated from different wild populations. Within a sample of 23 Thai swamp buffaloes, we identified four haplotypes with different geographical distributions, two of which were shared by Thai wild buffaloes.     

Using carbon isotope analysis of the diet of two introduced Australian megaherbivores to understand Pleistocene megafaunal extinctions

Aim Australia lost a diverse assemblage of large marsupial herbivores in the late Pleistocene, with suggestions that the extinctions were biased towards browsers. In modern times two bovines, the Asian water buffalo (Bubalus bubalis) and banteng (Bos javanicus), have established feral populations in the Northern Territory, Australia. Buffalo have aggressively expanded throughout the savanna landscape, yet banteng remain near their point of introduction on the Cobourg Peninsula. We hypothesized that this difference is related to feeding ecology, possibly reflecting a legacy of the Pleistocene extinctions. Location Western Arnhem Land, Northern Territory, Australia. Methods Analysing a previously published dataset of body mass and feeding ecology of extinct and extant marsupial herbivores, we evaluated whether browsers were at greater risk of extinction than grazers. We compared the carbon isotope composition and nitrogen content of banteng and buffalo dung in order to evaluate the hypotheses that the differences in invasion success are related to feeding ecology, and that seasonal variation in browse consumption is linked to changing nutritional quality of grass. Results Controlling for body mass, the Pleistocene extinctions were clearly biased towards browsers. Introduced banteng appear to be primarily browsers, with their diets comprising 40% grass in the wet season and 15% in the late dry season. Buffalo have a more variable diet, with an increasing proportion of browse from the wet (30%) to the late dry season (75%), and can therefore be described as switching from grazer to browser. The decline of grass in the diet of both species appears to reflect the decline in the nutritional value of grass through the dry season, an inference supported by the negative relationship between δ¹³C values and the nitrogen content of dung. Main conclusions Banteng and buffalo are much larger than extant native herbivores, of which browsers are restricted to isolated rocky habitats. This suggests that banteng and buffalo have filled niches made vacant following the Pleistocene extinctions. The success of buffalo appears to be related to their greater dietary breadth, which enables them to graze and browse in eucalypt savannas, whilst the browsing banteng remain tethered to a mosaic of rain forest patches. The restriction of browsers may be a long‐range consequence of habitat transformations associated with Aboriginal landscape burning.