Abstract
The gaur is the largest extant cattle species and distributed across South and Southeast Asia. Around 85% of its current global population resides in India, however there has been a gradual decrease in the gaur population over the last two decades due to various anthropogenic activities. Mitochondrial genome is considered as an important tool for species identification and monitoring the populations of conservation concern and therefore it becomes an obligation to sequence the mitochondrial genome of Indian gaur. We report here for the first time 16,345 bp mitochondrial genome of four Indian gaur sequenced using two different approaches. Mitochondrial genome consisted of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region. Among the 37 genes, 28 were positioned on the H-strand and 9 were positioned on the L-strand. The overall base composition appeared to be 33.5% A, 27.2% T, 25.9% C and 13.4% G, which yielded a higher AT content. The phylogenetic analysis using complete mitochondrial genome sequences unambiguously suggested that gaur is the maternal ancestor of domestic mithun. Moreover, it also clearly distinguished the three sub species of B. gaurus i.e. B. gaurus gaurus, B. gaurus readei and B. gaurus hubbacki. Among the three sub species, B. gaurus gaurus was genetically closer to B. gaurus readei as compared to B. gaurus hubbacki. The findings of our study provide an insight into the genetic structure and evolutionary history of Indian gaur.
Similar content being viewed by others
Introduction
The gaur, Bos gaurus also known as “Indian bison” is the largest living wild cattle species belonging to the family Bovidae1. The historical distribution of gaur ranged throughout much of the mainland South and Southeast Asia. But, currently it occurs in a few Asian countries such as Bangladesh, Bhutan, Cambodia, China, Malaysia, Myanmar, Nepal, Thailand, and Vietnam, with about 85% of its total population surviving in India2. Gaur population has declined drastically in almost its entire geographical range primarily due to habitat loss, poaching for horn and meat, diseases and competition for food resources1,2. As a result gaur has been categorized as vulnerable species by the IUCN3 and protected under schedule I of the wild life (protection) Act 1972 in India. Therefore conservation of gaur populations is pertinent.
The B. gaurus has been classified into several subspecies by different researchers elucidating ambiguities in their taxonomy. Lydekker4,5 reported three subspecies of gaur based on morphological descriptions namely Bos gaurus gaurus, which inhabits in India, Nepal and Bhutan, Bos gaurus readei, which inhabits in Cambodia, southern China, Lao PDR, Viet Nam, Myanmar, and Thailand and Bos gaurus hubbacki, which inhabits in Malaysia. Hubback6 opined the possibility for the presence of two types of gaur in Malaysia, one with no dewlap and one possessing well developed dewlap. Recently, Groves7 & Groves and Grubb8 proposed two sub species; B. gaurus gaurus which inhabits in India and Nepal; B. gaurus laosiensis which inhabits in Cambodia, Lao PDR, west Malaysia, Myanmar, Thailand, and Vietnam, based on skull and horn size. Another subspecies B. gaurus sinhaleyus has been recorded from Sri Lanka which is now extinct9. Notably, all the above mentioned studies have classified the subspecies based on morphological/skull features of only a few specimens which confound the taxonomic status of gaur. Considering the phenotypic differences, the International Union for Conservation of Nature (IUCN) has recognized only two sub species of gaur; B. gaurus gaurus and Bos gaurus laosiensis. However, there is no adequate genetic data to confirm the validity of sub species defined by the morphological data.
It is widely believed that gaur is the ancestral species of domestic mithun (Bos frontalis)10,11. Nevertheless, a few studies have suggested that mithun is a hybrid descendant of gaur and domestic cattle12,13 while a few others have reported mithun as a descendant of an unknown wild bovine which is already extinct14,15. Therefore, the origin of domestic mithun remains unresolved. Considering the above facts, genetic characterization of gaur holds great importance. Mitochondrial genome has been widely used to study the evolutionary and phylogenetic relationship of various species due to its high mutation rate and lack of recombination16,17,18,19. In the present study for the first time, we sequenced the complete mitochondrial genome of four Indian gaur sampled from different places of India. These sequences were analysed with that of Cambodian gaur, Malayan gaur and mithun to resolve the subspecies classification of gaur and shed light on the domestication history of mithun.
Materials and methods
Sample collection
In the present study four Indian gaur samples were used. Of the four samples, two were fresh dung samples (sample ID: GR01 & GR02) of free-ranging gaur collected from unprotected forest areas in Karnataka, India. The dung samples were collected within a few minutes of defecation by watching the gaur from a distance. The remaining two were blood and muscle samples of the dead gaur obtained from Arignar Anna Zoological Park, Vandalur, Chennai, (sample ID: GR03), Tamil Nadu and Periyar National Park, Idukki, Kerala, India (sample ID: GR04) respectively. Samples were collected with the help of forest officials/veterinarians for which necessary approval was taken from the concerned state forest departments. The collected dung and tissue samples were preserved in absolute ethanol and stored at − 20 °C until DNA extraction.
Mitochondrial genome sequencing from dung samples
Genomic DNA from the dung samples (GR01 & GR02) was extracted using QIAamp DNA Stool Mini Kit (Qiagen, Germany) as per the manufacturer’s instructions with few modifications. In brief, 1.5 ml buffer ASL was added to 250 mg dung sample (which contains the mucous layer), vortexed for 1 min, and incubated at 60 °C for 3 h. Subsequently the sample was centrifuged at 14,000 rpm for 3 min and the supernatant was transferred to a new micro centrifuge tube to which one InhibitEX tablet was added and vortexed for 1 min. Following 12 min centrifugation at 14,000 rpm, the supernatant was collected in a new micro centrifuge tube containing 20 µl Proteinase K and 600 µl buffer AL was added which was then vortexed for 15 s and incubated at 70 °C for 15 min. After the incubation, 600 µl absolute ethanol was added to the lysate and mixed thoroughly by vortexing. The lysate was transferred to the QIAamp spin column and centrifuged at 14,000 rpm for 1 min. Followed by, the QIAamp spin column was washed with supplied buffers AW1 and AW2 by centrifugation at 12,000 rpm for 1 min. Then, a high spin (14,000 rpm for 1 min) was given to dry the column membrane. The purified DNA was eluted twice (50 µl per elution) in buffer AE. The quality and quantity of the eluted DNA were checked by agarose gel electrophoresis and NanoDrop One (Thermo Fisher Scientific, USA) respectively.
The complete mitochondrial genome was amplified using 23 sets of overlapping primers20 (Supplementary Information Table S1). PCR amplifications were performed in 25 µl reaction mixture which included 80 ng of genomic DNA, 12.5 µl master mix (Promega, USA) and 2 µl (10 pmol) of each primer. The following conditions were applied to the PCR: 5 min initial denaturation at 94 °C followed by 30 cycles of denaturation at 94 °C for 1 min, annealing at 46–59 °C for 1 min, extension at 72 °C for 1 min; and the final extension at 72 °C for 7 min. The amplified PCR products were purified using QIAquick PCR Purification Kit (Qiagen, Germany) as per manufacturer’s instruction and sequenced using both forward and reverse primers. Sequencing was performed in a 10 μl scale using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Thermo Fisher Scientific, USA), which included 1 μl ready reaction mix, 1.5 μl sequencing buffer, 10 pmol primer and 200 ng PCR product. The following thermal cycle was applied for the amplification: 96 °C for 1 min, followed by 25 cycles of 96 °C for 10 s, 55 °C for 5 s and 60 °C for 4 min. The excess dye labelled terminators and buffers were removed by EDTA/ethanol precipitation at room temperature. Followed by DNA was denatured by adding 10 μl formamide at 95 ºC for 4 min. Sequencing was performed on ABI 3730XL DNA analyzer (Applied Biosystems, USA). The generated sequencing data was analyzed with Sequencing Analysis Software (Applied Biosystems, USA).
Mitochondrial genome sequencing from tissue samples
Genomic DNA from tissue samples (blood-GR03 and muscle-GR04) was extracted using DNeasy Blood and Tissue Kit (Qiagen, Germany) as specified by the manufacturer. The complete mitochondrial genome was amplified by long-range PCR using two overlapping sets of primers (5′-AATATGCTCGCCATCATTCC-3′, 5′-ATTGCAGAGGGAAGTCATGG-3′) and (5′-TCACCAGCATAATTCCCACA-3′, 5′-GGCATGTCACCAAGGAGAGT-3′). PCR was performed in 50 μl reaction mixture containing 10 μl of 5XPrimeSTAR GXL Buffer (Takara, Japan), 4 μl of dNTPs (2.5 mM each), 4 μl of each primer (15 pmol), 1 μl of PrimeSTAR GXL DNA Polymerase (Takara, Japan), 23 μl of nuclease free water and 4 μl of DNA template (300 ng). The following conditions were applied for the PCR: 94 °C for 5 min; 30 cycles of 94 °C for 1 min, 68 °C for 1 min and 68 °C for 10 min; 68 °C for 10 min.
The paired end (PE) libraries were prepared using NEBNext Ultra DNA Library Prep Kit (New England BioLabs, USA) as per the manufacturer’s instructions. In brief, the two amplified DNA fragments were pooled together in equimolar concentration and sonicated to a size of 300 bp using Covaris M220 (Covaris, USA). Subsequently the DNA fragments ends were repaired, ‘A’ tailed and ligated to indexed adapters. The resulting 300 bp adaptor-ligated DNA fragments were selected using sample purification beads and enriched through PCR amplification. The purity, size and concentration of the amplified libraries were analysed by Bioanalyzer (Agilent, USA). Finally, the PE libraries (2 × 150 bp) were sequenced on Illumina HiSeq X10 platform (Illumina, USA).
Mitochondrial genome sequence analysis
A total of 324,110 and 610,358 sequence reads were generated for the samples GR03 and GR04 respectively. The adaptor sequences were removed from the raw reads by Cutadapt v1.821. Further Sickle22 and FastUniq23 were used to remove low quality and duplicate reads respectively. After quality filtering, high quality reads were assembled using de novo assembly and reference based assembly (de novo assembled GR04 sequence was used as reference for GR03) using SeqManNGen (DNASTAR, USA) assembler. Similarly, the sequences generated through Sanger method was assembled using SeqManPro (DNASTAR, USA). The assembled mitochondrial genome sequences were edited and aligned using the software MegAlign Pro (DNASTAR, USA). MITOS web server was used to annotate the mitochondrial genomes24 followed by NCBI ORFfinder (https://www.ncbi.nlm.nih.gov/orffinder) and BLAST (https://blast.ncbi.nlm.nih.gov) were used to validate the annotations. Mitochondrial genome map was constructed using OGDRAW with default parameters25. The tRNA secondary structures were analyzed by tRNAscan-SE26 and MITOS web servers24. Nucleotide composition, Relative Synonymous Codon Usage (RSCU) values and genetic divergence between species were calculated using the software MEGA27. The AT and GC skewness were calculated as follows: AT skew = (A – T)/(A + T) and GC skew = (G – C)/(G + C)28. The overlapping regions and intergenic spacers between the genes were manually calculated. To ascertain the genetic relationship of Indian gaur with other Bos species a Bayesian phylogenetic tree was constructed using the general time reversible model on MrBayes29. The MCMC chains were run for 10 × 106 cycles. A total of 20,000 trees were sampled, and a 50% majority rule consensus tree was obtained with burnin = 5,000. The Maximum parsimony (MP) tree was constructed using the software MEGA with 5,000 bootstrap value27. The analysis of molecular variance (AMOVA) was calculated with the software ARLEQUIN30. The figures were drawn/edited using Inkscape 1.0 (https://inkscape.org).
Results and discussion
Mitochondrial genome organization
We have sequenced the complete mitochondrial genome of four Indian gaur using two different approaches. The complete mitochondrial genome of Indian gaur was 16,345 bp in length, which included 22 tRNA genes, 13 protein coding genes, 2 ribosomal RNA genes, and a control region (Table 1 and Fig. 1). The genome size and gene order were similar to previously reported gaur and mithun mitochondrial genomes17,18. The heavy (H) strand encoded most of the genes (28 of the 37 genes) except NADH dehydrogenase subunit 6 (nad6) and 8 tRNA genes (trnQ, trnA, trnN, trnC, trnY, trnS2, trnE and trnP) which were encoded by the Light (L) strand. The AT and GC content of the mitochondrial genome was observed to be 60.7% and 39.3% respectively, which indicated that the nucleotide composition is overall biased towards adenine and thymine (Fig. 2). This is a common trend among bovine species including mithun16,18,31,32. Moreover, the mitochondrial genome showed positive AT (0.104), and negative GC (-0.319) skews (Fig. 2) which suggested the higher content of adenine and cytosine than their respective complementary nucleotides guanine and thymine. In total there were 72 overlapping nucleotides in the range from 1 to 40 bp, which were found at 7 distinct locations. The largest overlapping region (40 bp) was observed between the two protein coding genes ATP synthase F0 subunit 8 (atp8) and ATP synthase F0 subunit 6 (atp6). In addition to these, the intergenic spacer (IGS) was interspersed at 14 regions across the mitochondrial genome with varying range from 1 to 32 bp which summed up to a total length of 75 bp. The largest intergenic spacer was located between the two tRNA genes trnN and trnC.
Protein coding genes (PCGs)
Mitochondrial genome encoded 13 PCGs and was observed to be 11,339 bp in length which accounted for 69.37% of the mitochondrial genome. The AT and GC content was 60.4% and 39.6% respectively (Fig. 2; Supplementary Information Table S2) which revealed the nucleotide compositional biasness of the PCGs towards adenine and thymine. Moreover, the AT and GC skews of PCGs were positive (0.047) and negative (-0.336) respectively as observed in the case of whole mitochondrial genome (Fig. 2) which suggested that adenine content is relatively higher than thymine while cytosine content is higher than guanine. Also, notably, all the PCGs showed positive AT skew except for the genes cox1, cox3, and nad6, whereas GC skew was positive only in nad6 gene. As commonly observed in other bovine and vertebrate species18,19,32,33, the PCGs were subdivided into seven NADH dehydrogenase subunits (nad1, nad2, nad3, nad4, nad4L, nad5 and nad6), three cytochrome c oxidase (cox1, cox2 and cox3), two ATPase subunits (atp8 and atp6) and a cytochrome b gene (cob). The size of the PCGs varied significantly with atp8 (201 bp) being the shortest and nad5 (1821 bp) being the longest among all. Besides, there were four adjacent pairs of PCGs (atp8-atp6, atp6-cox3, nad4L-nad4 and nad5-nad6) which is a common gene positioning observed among the vertebrates18,19,33. The relative synonymous codon usage (RSCU) analysis showed the highest utilization of CUA, GGA, GUA and CGA codons among PCGs (Fig. 3A). From the RSCU pattern, it could be understood that among the synonymous alternative codons of each amino acid, the codons associated with adenine at their third codon position were more preferred. The amino acids leucine, threonine, proline and isoleucine were the most abundant among the PCGs (Fig. 3B).
Ribosomal RNA and transfer RNA
The total size of the rRNA was 2,525 bp which was formed of two subunits 12S rRNA (956 bp) and 16S rRNA (1569 bp). The nucleotide composition of 12S rRNA as well as 16S rRNA exhibited biasness towards AT content. The AT skew of 12S rRNA and 16S rRNA was positive whereas GC skew was negative, which showed the relatively higher occurrence of adenine and cytosine than thymine and guanine in the rRNAs (Fig. 2). There were 22 tRNA genes in the mitochondrial genome which varied in size from 60 bp (trnS1) to 75 bp (trnL2). Most of the tRNA genes were encoded by the H-strand (trnF, trnV, trnL2, trnI, trnM, trnW, trnD, trnK, trnG, trnR, trnH, trnS1, trnL1, trnT) while the tRNAs trnQ, trnA, trnN, trnC, trnY, trnS2, trnE and trnP were encoded by the L-strand. All the tRNA genes exhibited the cloverleaf secondary structure except trnS1 and trnK which lacked a stable dihydrouridine arm loop (Fig. 4). Such unusual tRNA structures have been commonly observed in mammals including the closely related domestic mithun18,34. The AT and GC skews were positive for the tRNAs, which revealed higher compositional count of adenine and guanine than their respective complementary nucleotides (Fig. 2).
Control region
The control region of the mitochondrial genome was 902 bp in length which was positioned between the tRNAs trnP and trnF. The AT and GC skews were positive and negative respectively (Fig. 2). The palindromic motifs ‘TACAT’ and ‘ATGTA’ that tend to form hairpin loop structures were found dispersed in the control region in multiple copies. Such identical motifs have also been observed in domestic mithun and other closely related bovine species as well as in other vertebrates18,19,35,36, and these are believed to function as the termination site for the elongation of heavy strand37. The control region terminated with a characteristic poly-C stretch (13 nucleotides) which was also observed in domestic mithun18 (Supplementary Information Table S3).
Phylogenetic analysis
Overall, there were no differences in size, and gene order and organization among the four mitochondrial genome sequences of Indian gaur. But they differed in nucleotide composition, as there were 33 substitutions for the entire 16,345 bp long mitochondrial genome. It shows extremely low levels of genetic diversity among Indian gaur even though the four samples were collected from far distant places. Our result is in line with a recent study which also observed low genetic diversity among gaur populations of Central India based on partial mtDNA D-loop sequences38. Similar low mtDNA diversity was also observed in Malayan gaur for partial D-loop sequence39. Low mtDNA diversity is likely to occur in the wild population due to the presence of small number of founder females38,39,40,41. Also, it can be taken as a sign of population decline40,42,43. Genetic variation between individuals is prerequisite for evolution and adaptive changes, and has profound implications for conservation41,44,45,46. The global gaur population is anticipated to fall by 30% within next three decades3 due to habitat loss, poaching for its meat and horns and fatal diseases. However, the decline of gaur population in India is considerably lower as compared to other countries such as Bangladesh, Cambodia, China, Laos, Malaysia, Vietnam, and Thailand. Yet the low genetic diversity estimates indicate that a substantial effort needs to be taken immediately to design and implement strategies for the conservation of threatened Indian gaur germplasm. Further analyses using SNP and microsatellite markers may reveal exactly the current status and population structure of Indian gaur.
The Indian gaur mitochondrial genomes were further compared with the mitochondrial genome of Cambodian gaur (JN632604) and Malayan gaur (MK770201). The Cambodian gaur showed 304 substitutions with Indian gaur whereas 634 substitutions were observed between Indian and Malayan gaur. The AMOVA revealed 94% variation between Indian and Cambodian gaur and that between Indian and Malayan gaur was 97%, which indicated that there is significant genetic variation among Indian, Cambodian and Malayan gaur. Thus, in order to understand the phylogenetic relationship of Indian, Cambodian and Malayan gaur, phylogenetic trees were constructed along with seven congeneric species viz B. frontalis, B. taurus, B. indicus, B. javanicus, B. mutus, B. grunniens and B. primigenius. For the congeneric species, a maximum of five sequences for each species or total available sequences were included in the analysis. The African buffalo, Syncerus caffer was used as an out-group. Bayesian phylogenetic tree showed three distinct clades and each clade was further divided by the wild and domestic individual. The MP tree was identical to the former, with strong bootstrap support. The clade one included wild and domestic cattle while the clade two was fully encompassed with wild and domestic yak. Similar to clade one and two, clade three comprised of wild gaur and domestic mithun which demonstrated the ancestral connections between gaur and domestic mithun (Fig. 5). The banteng, B. javanicus was distributed in all the clades except clade 2 along with taurine cattle, domestic mithun and gaur which shows the hybrid nature of B. javanicus as reported elsewhere47. The presence of wild gaur and domestic mithun in a single clade as observed in the case of cattle and yak (Fig. 5) unambiguously suggests that wild gaur is the maternal ancestor of domestic mithun. Similar findings have also been obtained in previous studies using 16S rRNA gene10, SNP11, cytochrome b gene48,49 and Y chromosomal DNA50 markers. Further, AMOVA revealed less (34%) variation between mithun and gaur as compared to mithun and cattle (> 97%). Genetic divergence was also much lower (0.031) between gaur and mithun than between mithun and the other two cattle species (\(\ge\) 0.052). Moreover, studies based on descriptive characteristics, protein polymorphism, karyotype, and microsatellite marker have showed significant differences between domestic mithun and cattle11,51,52,53,54,55. Therefore, based on the results of present and previous studies, we strongly suggest that, wild gaur is the maternal ancestor of domestic mithun.
On the other hand, one of the Chinese mithun clustered with B. indicus (Fig. 5), but its nuclear genome did not support this clustering13. Introgression of cattle mitochondrial DNA into mithun has been reported by Li et al.48 and Gou et al.56. These outcomes demonstrate that some mithun have mithun chromosomal genome and cattle mitochondrial genome. Hybridization has been practiced between mithun and cattle to obtain animal of higher economic value since historical times51. Studies have also reported the transportation of mithun from India to Bhutan during the ancient times, where they were extensively used to cross breed with domestic cattle57. It is therefore, possible that mithun individuals analysed in these studies were likely to be the hybrid of domestic mithun and domestic cattle which is prevalent in China16,56. Thus, in domestic animals phenotype is not always related to their mitochondrial genome58, which leads to phenotypic ambiguity among domestic species and impede conservation efforts. Identification of genetically unique population is critical for species conservation hence, an extensive study is essential to evaluate the cross bred animals on a larger scale using both chromosomal and mitochondrial DNA markers.
The Indian gaur, Cambodian gaur and Malayan gaur formed distinct clades within the third clade with strong Bayesian posterior probability and bootstrap values (Fig. 5). Our finding therefore, clearly supports the existence of three sub species of gaur i.e. B. gaurus gaurus, B. gaurus readei and B. gaurus hubbacki proposed based on the morphological features by Lydekker4,5. Further, the phylogenetic tree revealed that among three subspecies, the B. gaurus gaurus is genetically closer to B. gaurus readei as compared to B. gaurus hubbacki. Nevertheless more number of samples from each subspecies needs to be analysed using both the nuclear and mitochondrial DNA markers to confirm their identity. Also, the domestic mithun including Indian mithun samples were clustered with Cambodian gaur which explains the close genetic relationship of domestic mithun and Cambodian gaur. However, at this juncture it is premature to suggest the place of domestication of mithun while considering the sample size of the present study particularly the gaur of Northeast India which is reported to closely resemble the Southeast Asian gaur7 would have significant implications on the origin of domestic mithun. On the whole, we suggest that it is not only important but necessary to conduct a detailed study on complete mitochondrial genome of gaur from different countries particularly Northeast India, Myanmar, Malaysia and China to unveil the time and place of domestication of mithun. On the other hand, studies on whole mitochondrial genome of domestic mithun from different countries would further complement to establish the domestication history of mithun.
Conclusions
The study on complete mitochondrial genome of Indian gaur is of paramount importance primarily for two reasons (i) gaur is considered as the ancestral species of domestic mithun but still is in dispute and (ii) there is a taxonomic ambiguity with regard to the classification of gaur into sub species. This is the first study to sequence the complete mitochondrial genome of Indian gaur. The findings of our study clearly conclude that, gaur is the maternal ancestor of domestic mithun. Although, Indian gaur mitochondrial genome shared similar structural features with Cambodian and Malayan gaur, a significant genetic variation was observed between them which support the existence of three sub species of gaur. Further, the low genetic diversity of Indian gaur indicates the necessity for developing appropriate conservation strategies for gaur populations in India. The complete mitochondrial genome sequence of gaur would serve as a useful genetic resource for phylogenetic, evolutionary biology, and conservation related studies.
Data availability
The assembled mitochondrial genome sequences are available at GenBank under the following Accession Nos. MT345892, MT345893, MT360652 and MT360653. The raw sequence reads can be found at Sequence Read Archive (SRA) under the BioProject ID: PRJNA627336.
References
Ahrestani, F. S. Bos frontalis and Bos gaurus (Artiodactyla: Bovidae). Mamm. Species 50, 34–50 (2018).
Ashokkumar, M., Swaminathan, S., Nagarajan, R. & Desai, A. A. Distribution, ecology and conservation of the gaur. In Animal diversity, natural history and conversation (eds Gupta, V. K. & Verma, A. K.) 77–94 (Daya Publishing House, New Delhi, 2011).
Duckworth, J. W., Sankar, K., Williams, A. C., Samba Kumar, N. & Timmins, R. J. Bos gaurus. The IUCN Red List of Threatened Species 2016: e.T2891A46363646. https://doi.org/10.2305/IUCN.UK.2016-2.RLTS.T2891A46363646.en (2016).
Lydekker, R. The Burmese gaur or pyoung. Zoologist 7, 264–266 (1903).
Lydekker, R. The game animals of India, Burma, Malaya and Tibet (Rowland Ward, London, 1907).
Hubback, T. R. The Malayan gaur or seladang. J. Mammal. 18, 267–279 (1937).
Groves, C. P. Taxonomy of ungulates of the Indian subcontinent. J. Bombay. Nat. Hist. Soc. 100, 341–362 (2003).
Groves, C. & Grubb, P. Ungulate taxonomy (John Hopkins University Press, Baltimore, 2011).
Grubb, P. Order artiodactylain mammal species of the world: a taxonomic and geographic reference (eds. D.E. Wilson and D.M. Reeder) 637–722 (Johns Hopkins University Press, Baltimore, 2005).
Dorji, T., Mannen, H., Namikawa, T., Inamura, T. & Kawamoto, Y. Diversity and phylogeny of mitochondrial DNA isolated from mithun Bos frontalis located in Bhutan. Anim. Genet. 41, 554–556 (2010).
Mukherjee, A. et al. High-density genotyping reveals genomic characterization, population structure and genetic diversity of Indian Mithun (Bos frontalis). Sci. Rep. 8, 10316. https://doi.org/10.1038/s41598-018-28718-x (2018).
Payne, W. J. A. Cattle production in the tropics in breeds and breeding Vol. 1 (Longman, London, 1970).
Mei, C. et al. Whole-genome sequencing of the endangered bovine species Gayal (Bos frontalis) provides new insights into its genetic features. Sci. Rep. 6, 19787. https://doi.org/10.1038/srep19787 (2016).
Ma, G. et al. Phylogenetic relationships and status quo of colonies for gayal based on analysis of cytochrome b gene partial sequences. J. Genet. Genomics 34, 413–419 (2007).
Baig, M. et al. Mitochondrial DNA diversity and origin of Bos frontalis. Curr Sci. 104, 115–120 (2013).
Ren, Q. et al. Complete mitochondrial genome of bovine species Gayal (Bos frontalis). Conserv. Genet. Resour. 10, 889–891 (2018).
Rosli, N. et al. The complete mitochondrial genome of Malayan Gaur (Bos gaurus hubbacki) from Peninsular Malaysia. Mitochondrial DNA B 4, 2535–2536 (2019).
Prabhu, V. R., Arjun, M. S., Bhavana, K., Kamalakkannan, R. & Nagarajan, M. Complete mitochondrial genome of Indian mithun, Bos frontalis and its phylogenetic implications. Mol. Biol. Rep. 46, 2561–2566 (2019).
Prabhu, V. R., Singha, H. S., Kumar, R. G., Gopalakrishnan, A. & Nagarajan, M. Characterization of the complete mitochondrial genome of Barilius malabaricus and its phylogenetic implications. Genomics 112, 2154–2163 (2020).
Hassanin, A., Ropiquet, A., Couloux, A. & Cruaud, C. Evolution of the mitochondrial genome in mammals living at high altitude: new insights from a study of the tribe Caprini (Bovidae, Antilopinae). J. Mol. Evol. 68, 293–310 (2009).
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10–12 (2011).
Joshi, N. & Fass, J. Sickle: a sliding-window, adaptive, quality-based trimming tool for fastQ files (version 1.33) [software]. https://github.com/najoshi/sickle (2011).
Xu, H. et al. FastUniq: a fast de novo duplicates removal tool for paired short reads. PLoS ONE 7, e52249. https://doi.org/10.1371/journal.pone.0052249 (2012).
Bernt, M. et al. MITOS: improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol. 69, 313–319 (2013).
Lohse, M., Drechsel, O., Kahlau, S. & Bock, R. OrganellarGenomeDRAW-a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucl. Acids. Res. 41, W575–W581 (2013).
Lowe, T. M. & Chan, P. P. tRNAscan-SE on-line: integrating search and context for analysis of transfer RNA genes. Nucl. Acids. Res. 44, W54–W57 (2016).
Kumar, S., Stecher, G. & Tamura, K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874 (2016).
Perna, N. T. & Kocher, T. D. Patterns of nucleotide composition at four fold degenerate sites of animal mitochondrial genomes. J. Mol. Evol. 41, 353–358 (1995).
Ronquist, F. & Huelsenbeck, J. P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19, 1572–1574 (2003).
Excoffier, L. & Lischer, H. E. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567 (2010).
Meng, W. et al. The complete mitochondrial genome of Gazella subgutturosa yarkandens (Artiodactyla; Bovidae; Antilopinae) revealed by next-generation sequencing and its phylogenetic implications. Conserv. Genet. Resour. 10, 747–749 (2018).
Zhou, X. et al. The complete mitochondrial genome of Sunan yak (Bos grunniens) and a phylogenetic study. Mitochondrial DNA B 4, 1053–1054 (2019).
Anderson, S. et al. Sequence and organization of the human mitochondrial genome. Nature 290, 457–465 (1981).
Manee, M. M. et al. Comparative analysis of camelid mitochondrial genomes. J. Genet. 98, 88. https://doi.org/10.1007/s12041-019-1134-x (2019).
Wang, X., Wang, J., He, S. & Mayden, R. L. The complete mitochondrial genome of the Chinese hook snout carp Opsariichthys bidens (Actinopterygii: Cypriniformes) and an alternative pattern of mitogenomic evolution in vertebrate. Gene 399, 11–19 (2007).
Zhou, M., Yu, J., Li, B., Ouyang, B. & Yang, J. The complete mitochondrial genome of Budorcas taxicolor tibetana (Artiodactyla: Bovidae) and comparison with other Caprinae species: Insight into the phylogeny of the genus Budorcas. Int. J. Biol. Macromol. 121, 223–232 (2019).
Saccone, C., Pesole, G. & Sbisa, E. The main regulatory region of mammalian mitochondrial DNA: structure–function model and evolutionary pattern. J. Mol. Evol. 33, 83–91 (1991).
Atkulwar, A., Farah, S., Gadhikar, Y. & Baig, M. Mitochondrial DNA diversity in wild gaur (Bos gaurus gaurus): evidence from extant and historical samples. Mitochondrial DNA B 5, 1556–1560 (2020).
Md-Zain, B. M. et al. Sequence variation of captive Malayan gaur (Bos gaurus hubbacki) based on mitochondrial D-loop region DNA sequences. Biodiversitas 19, 1601–1606 (2018).
Wang, W., Qiao, Y., Li, S., Pan, W. & Yao, M. Low genetic diversity and strong population structure shaped by anthropogenic habitat fragmentation in a critically endangered primate Trachypithecus leucocephalus. Heredity 118, 542–553 (2017).
Wang, W., Zheng, Y., Zhao, J. & Yao, M. Low genetic diversity in a critically endangered primate: shallow evolutionary history or recent population bottleneck?. BMC Evol. Biol. 19, 134. https://doi.org/10.1186/s12862-019-1451-y (2019).
Ahmad, K. et al. Genetic diversity of the Tibetan antelope (Pantholops hodgsonii) population of Ladakh, India, its relationship with other populations and conservation implications. BMC Res. Notes. 9, 477. https://doi.org/10.1186/s13104-016-2271-4 (2016).
Khaire, D., Atkulwar, A., Farah, S. & Baig, M. Mitochondrial DNA analyses revealed low genetic diversity in the endangered Indian wild ass Equus hemionuskhur. Mitochondrial DNA A 28, 681–686 (2017).
Reed, D. H. & Frankham, R. Correlation between fitness and genetic diversity. Conserv. Biol. 17, 230–237 (2003).
Ørsted, M., Hoffmann, A. A., Sverrisdóttir, E., Nielsen, K. L. & Kristensen, T. N. Genomic variation predicts adaptive evolutionary responses better than population bottleneck history. PLoS Genet. 15, e1008205. https://doi.org/10.1371/journal.pgen.1008205 (2019).
Kleinhans, C. & Willows-Munro, S. Low genetic diversity and shallow population structure in the endangered vulture Gyps coprotheres. Sci. Rep. 9, 5536. https://doi.org/10.1038/s41598-019-41755-4 (2019).
Nijman, I. J. et al. Hybridization of banteng (Bos javanicus) and zebu (Bos indicus) revealed by mitochondrial DNA, satellite DNA, AFLP and microsatellites. Heredity 90, 10–16. https://doi.org/10.1038/sj.hdy.6800174 (2003).
Li, S. P., Chang, H., Ma, G. L. & Cheng, H. Y. Molecular phylogeny of the gayal in Yunnan China inferred from the analysis of cytochrome b gene entire sequences. Asian-Aust. J. Anim. 21, 789–793 (2008).
Tanaka, K. et al. Molecular phylogeny and diversity of Myanmar and Bhutan mithun based on mtDNA sequences. Anim. Sci. J. 82, 52–56 (2011).
Nijman, I. J. et al. Phylogeny of Y chromosomes from bovine species. Cladistics 24, 723–726 (2008).
Winter, H. et al. Karyotyping, red blood cell and haemoglobin typing of the mithun (Bos frontalis), its wild ancestor and its hybrids. Res. Vet. Sci. 36, 276–283 (1984).
Nie, L. et al. Genetic diversity of cattle in South China as revealed by blood protein electrophoresis. Biochem. Genet. 37, 257–265 (1999).
Tu, Z., Nie, L., Yu, Y., Wen, J. & Zhang, Y. Blood protein polymorphism in B. frontalis, B. grunniens, B. taurus, and B. indicus. Biochem. Genet. 38, 413–416 (2000).
Wang, M. S. et al. Draft genome of the gayal, Bos frontalis. Gigascience 6, 094. https://doi.org/10.1093/gigascience/gix094 (2017).
Ritz, L. R., Glowatzki-Mullis, M. L., MacHugh, D. E. & Gaillard, C. Phylogenetic analysis of the tribe Bovini using microsatellites. Anim. Genet. 31, 178–185 (2000).
Gou, X., Wang, Y., Yang, S., Deng, W. & Mao, H. Genetic diversity and origin of Gayal and cattle in Yunnan revealed by mtDNA control region and SRY gene sequence variation. J. Anim. Breed. Genet. 127, 154–160 (2010).
Simoons, F. J. Gayal or mithan in Evolution of domesticated animals (ed. Mason, L. L.) 34–39 (Longman, London. 1984).
Srirattana, K., McCosker, K., Schatz, T. & John, J. C. S. Cattle phenotypes can disguise their maternal ancestry. BMC Genet. 18, 59. https://doi.org/10.1186/s12863-017-0523-5 (2017).
Acknowledgements
The study was partly supported by the Science and Engineering Research Board (SERB), Department of Science and Technology, (EMR/2015/000937), Government of India, New Delhi. Karippadakam Bhavana and Vandana R. Prabhu were supported by UGC (365844) and DST-INSPIRE (IF160266) through PhD fellowship respectively. The authors are sincerely thankful to the Forest department of Kerala, Karnataka, and Tamil Nadu, for giving permission to carry out the fieldwork and their support during sample collection. We also owe thanks to the Additional Principal Chief Conservator of Forest & Director, and Veterinary Officer, Arignar Anna Zoological Park, Chennai, Tamil Nadu for his support during the sample collection.
Author information
Authors and Affiliations
Contributions
M.N. conceived the study. K.B., R.K. and D.S performed the experiments. R.K. and M.N. performed the analyses. V.R.P., R.K, H.S.S. and M.N. wrote the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kamalakkannan, R., Bhavana, K., Prabhu, V.R. et al. The complete mitochondrial genome of Indian gaur, Bos gaurus and its phylogenetic implications. Sci Rep 10, 11936 (2020). https://doi.org/10.1038/s41598-020-68724-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-68724-6
This article is cited by
-
Ancestral alleles defined for 70 million cattle variants using a population-based likelihood ratio test
Genetics Selection Evolution (2024)
-
Mitochondrial genome of the Indian spot-billed duck and its phylogenetic and conservation implications
Scientific Reports (2024)
-
Characterization of Complete Mitochondrial Genome of Badri Breed of Bos indicus (Bovidae: Bovinae): Selection Pressure and Comparative Analysis
Biochemical Genetics (2024)
-
First description of the mitochondrial genomes of the Central American brocket deer Mazama temama (Kerr, 1792) and the Yucatán Peninsula brocket deer Odocoileus pandora Merriam, 1901
Molecular Biology Reports (2023)
-
Gaur genome reveals expansion of sperm odorant receptors in domesticated cattle
BMC Genomics (2022)