Production of recombinant buffalo (Bubalus bubalis) and goat (Capra hircus) growth hormones from genetically modified E. coli strains

The growth hormone cDNAs of Indian reverine buffalo (Bubalus bubalis) and beetal goat (Capra hircus) were cloned in Escherichia coli through RT-PCR technique. Nucleotide sequencing revealed several silent mutations in both cDNAs and only one amino acid change in the case of goat when compared to reported bovine (Bos taurus) sequence. The high level expression of both the polypeptide hormones was achieved in E. coli (> or =30% of soluble intracellular proteins) through the construction of two-cistronic gene expression system. The solubilisation of recombinant growth hormones from inclusion bodies and subsequent oxidation to correctly folded monomeric form was also carried out. A combination of reverse-phase HPLC and non-reducing SDS-PAGE was successfully applied to distinguish between reduced and oxidised forms of growth hormones. A moderate yield ( approximately 40% of starting material, with potential for upscaling), two-step purification process comprising of hydrophobic interaction and ion-exchange chromatographies was developed. The process eliminates the need for costly, laborious and time-consuming steps of ultrafiltration and dialysis, as reported earlier for the purification of many recombinant animal growth hormones. The biophysical, biochemical and functional analyses of purified refolded polypeptides showed that the hormones produced in this study were identical to natural pituitary bovine growth hormone.     

An insight into the possible mechanism of working of two-cistronic gene expression systems and rational designing of newer systems

The initial attempts at hyper-expressing buffalo/goat growth hormone (GH)-ORFs inEscherichia coli directly under various strong promoters were not successful despite the presence of a functional gene. High level expression of GH was achieved as a fusion protein with glutathione-S-transferase (GST). To produce native GH in an unfused state, we adapted an established strategy of two-cistronic approach in our system. In this strategy, utilizing one of the highly efficient reported sequences as the first cistron led to a nearly 1000-fold enhancement in the level of expression under anE. coli promoter (trc). In search of a newer first-cistron sequence as well as to see the generality of the two-cistronic approach, we explored the ability of different lengths of a highly expressing natural gene to act as an efficient first cistron. Surprisingly,GST, which is naturally highly expressible inE. coli, could not be fitted into a successful two-cistronic construct. In addition, placement of the entire two-cistronic expression cassette (which had earlier given high-level GH expression undertrc promoter) under theT7 promoter inE. coli failed to hyper-express GH. These results suggest that the successful exploitation of the two-cistron arrangement for hyper-expression of eukaryotic ORFs in bacteria is not as straightforward as was previously thought. It appears probable that factors such as the sequence context, together with the length and codons used in the first cistron are important as well.     

Molecular cloning,sequence characterization and recombinant expression of Nanog gene in goat fibroblast cells using lentiviral based expression system

Nanog is a homeodomain containing protein which plays important roles in regulation of signaling pathways for maintenance and induction of pluripotency in stem cells. Because of its unique expression in stem cells it is also regarded as pluripotency marker. In this study goat Nanog (gNanog) gene has been amplified, cloned and characterized at sequence level with successful over-expression in CHO-K1 cell line using a lentiviral based system. gNanog ORF is 903 bp long which codes for Nanog protein of size 300 amino acids (aas). Complete nucleotide sequence shows some evolutionary mutation in goat in comparision to other species. Protein sequence of goat is highly similar to other species. Overall, gNanog nucleotide sequence and predicted protein sequence showed high similarity and minimum divergence with cattle (96 % identity/4 % divergence) and buffalo (94/5 %) while low similarity and high divergence with pig (84/15 %), human (81/23 %) and mouse (69/40 %) indicating evolutionary closeness of gNanog to cattle and buffalo. gNanog lentiviral expression construct was prepared for over-expression of Nanog gene in adult goat fibroblast cells. Lentiviral expression construct of Nanog enabled continuous protein expression for induction and maintenance of pluripotency. Western blotting revealed the expression of Nanog gene at protein level which supported that the lentiviral expression system is highly promising for Nanog protein expression in differentiated goat cell.     

Sequence analysis of UTR and coding region of kappa-casein gene of Indian riverine buffalo (Bubalus bubalis).

In this study, complete nucleotide as well as derived amino acid sequence characterization of water buffalo (Bubalus bubalis) kappa-casein gene has been presented. Kappa-casein cDNA clones were identified and isolated from a buffalo lactating mammary gland cDNA library. Sequence analysis of kappa-casein cDNA revealed 850 nucleotides with an open reading frame (ORF) of 573 nucleotides, encoding mature peptide of 169 amino acids. The 5' untranslated region (UTR) comprised 71 nucleotides, while 3' UTR was of 206 nucleotides. A total of 11 nucleotide and seven amino acid changes were observed in, buffalo (Bubalus bubalis) as compared to cattle (Bos taurus), sheep (Ovis aries) and goat (Capra hircus). Among these nucleotide changes, eight were unique in buffalo as they were fully conserved in cattle, sheep and goat. Majority of the nucleotide changes and all the amino acid changes; 14 (Asp-Glu), 19(Asp/Ser-Asn), 96(Ala-Thr), 126(Ala-Val), 128(Ala/Gly-Val), 156(Ala/Pro-Val) and 168(Ala/Glu-Val) were limited to exon IV. Three glycosylation sites, Thr 131, Thr 133 and Thr 142 reported in cattle and goat kappa-casein gene were also conserved in buffalo, however, in sheep Thr 142 was replaced by Ala. Chymosin hydrolysis site, between amino acids Phe 105 and Met 106, important for rennet coagulation process, were found to be conserved across four bovid species. Buffalo kappa-casein with the presence of amino acids Thr 136 and Ala 148 seems to be an intermediate of "A" and "B" variants of cattle. Comparison with other livestock species revealed buffalo kappa-casein sharing maximum nucleotide (95.5%) and amino acid (92.6%) similarity with cattle, whereas with pig it showed least sequence similarity of 76.0% and 53.2%, respectively. Phylogenetic analysis based on both nucleotide and amino acid sequence indicated buffalo kappa-casein grouping with cattle, while sheep and goat forming a separate cluster close to them. The non-ruminant species viz. camel, horse and pig were distantly placed, in separate lineages.     

Molecular Cloning,Sequence, and Phylogenetic Analysis of T Helper 1 Cytokines of Pashmina Goats

Cytokines play an important role in regulation of immune responses either in health or disease. In the present study, the cDNAs encoding mature Interleukin (IL)-2, interferon gamma (IFN-γ), and IL-12 p35 and p40 of Pashmina goat were cloned and sequenced. The amino acid sequence was deduced from nucleotide sequence and compared with those available in GeneBank. Mature forms of goat IL-2, IFN-γ, IL-12 p35, and IL-12 p40 composed of 135, 143, 196, and 305 amino acid residues, respectively. Comparison of amino acid sequence of goat IL-2 with sheep, buffalo, cattle, pig, camel, cat, and human sequences showed homology percentages of 100, 97.8, 96.3, 72.4, 72.4, 67.2, and 64.7, respectively. Amino acid sequence of goat IFN-γ showed 98.6, 95.8, 81.1, 81.8, 80.4, and 62.9 percent homology with sheep, bovine, pig, horse, dog, and humans, respectively. Homology ranging from 81.6 to 99% for IL-12 p35 sequences and 85.6 to 100% for IL-12 p40 sequences at amino acid level were observed across these species. Multiple sequence alignment and phylogenetic analysis of goat cytokines revealed close relationship with sheep sequence.     

Pstl repeat: a family of short interspersed nucleotide element (SINE)-like sequences in the genomes of cattle, goat, and buffalo.

The PstI family of elements are short, highly repetitive DNA sequences interspersed throughout the genome of the Bovidae. We have cloned and sequenced some members of the PstI family from cattle, goat, and buffalo. These elements are approximately 500 bp, have a copy number of 2 x 10(5) - 4 x 10(5), and comprise about 4% of the haploid genome. Studies of nucleotide sequence homology indicate that the buffalo and goat PstI repeats (type II) are similar types of short interspersed nucleotide element (SINE) sequences, but the cattle PstI repeat (type I) is considerably more divergent. Additionally, the goat PstI sequence showed significant sequence homology with bovine serine tRNA, and is therefore likely derived from serine tRNA. Interestingly, Southern hybridization suggests that both types of SINEs (I and II) are present in all the species of Bovidae. Dendrogram analysis indicates that cattle PstI SINE is similar to bovine Alu-like SINEs. Goat and buffalo SINEs formed a separate cluster, suggesting that these two types of SINEs evolved separately in the genome of the Bovidae.     

Detection of polymorphism and sequence characterization of toll-like receptor 7 gene of Indian goat revealing close relationship between ruminant species

In this study, approximately 3.4?kb nucleotide sequence of caprine TLR7 (Toll-like receptor 7) gene was generated from twelve different Indian goat breeds belonging to different geographical regions. Goat TLR7 gene ORF (Open Reading Frame) was found to be 3141 nucleotides long coding for 1046?amino acids similar to sheep. The sequence analysis at nucleotide level revealed goat TLR7 having 99.5% homology with sheep, followed by other livestock species. Simple Modular Architecture Research Tool (SMART) was used for the structural analysis of goat TLR7 that showed the presence of 22 leucine rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domains. TIR domain, when compared, was found to be similar in ruminant species, goat, sheep, cattle, and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to sheep followed by cattle and buffalo. A total of 22 polymorphic sites were observed in TLR7 gene of 24 goats representing 12 different breeds, out of which 19 were present within the coding region and three in 3'UTR. Out of the seven nonsynonymous SNPs, two were in ectodomains and one in TIR domain. Overall our results indicate substantial variation within goat TLR7 gene, which could be exploited for association with disease susceptibility.     

Detection of Polymorphism and Sequence Characterization of Toll-Like Receptor 7 Gene of Indian Goat Revealing Close Relationship Between Ruminant Species

In this study, approximately 3.4 kb nucleotide sequence of caprine TLR7 (Toll-like receptor 7) gene was generated from twelve different Indian goat breeds belonging to different geographical regions. Goat TLR7 gene ORF (Open Reading Frame) was found to be 3141 nucleotides long coding for 1046 amino acids similar to sheep. The sequence analysis at nucleotide level revealed goat TLR7 having 99.5% homology with sheep, followed by other livestock species. Simple Modular Architecture Research Tool (SMART) was used for the structural analysis of goat TLR7 that showed the presence of 22 leucine rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domains. TIR domain, when compared, was found to be similar in ruminant species, goat, sheep, cattle, and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to sheep followed by cattle and buffalo. A total of 22 polymorphic sites were observed in TLR7 gene of 24 goats representing 12 different breeds, out of which 19 were present within the coding region and three in 3'UTR. Out of the seven nonsynonymous SNPs, two were in ectodomains and one in TIR domain. Overall our results indicate substantial variation within goat TLR7 gene, which could be exploited for association with disease susceptibility.     

Isolation and sequence characterization of mammary derived growth inhibitor gene of riverine buffalo (Bubalus bubalis)

In this study, attempts have been made to identify and characterize water buffalo (Bubalus bubalis) mammary derived growth inhibitor (MDGI) gene, isolated from a mammary gland cDNA library of lactating buffalo. The complete MDGI cDNA was of 698 nucleotides, consisting 61 nucleotides in 5' UTR, coding region of 402 nucleotides, and 235 nucleotides representing the 3' UTR. Comparison of nucleotide and deduced amino acid sequence data with that of MDGI/fatty acid binding protein (FABP) of other species shows three buffalo specific nucleotide changes while seven nucleotide changes were common to cattle and buffalo. Buffalo and cattle MDGI had 100% amino acid sequence similarity, which also shared three amino acid changes: 34 (Ala-Gly), 109 (Leu-Met), and 132 (Glu-Gln) as compared to other species. Comparison with FABPs reported from other cattle tissues revealed highest amino acid sequence similarity with FABP-heart (100%) and least with FABP-liver (20.5%). Phylogenetic analysis revealed cattle MDGI to be closest to buffalo, while mouse MDGI was distantly placed, whereas different tissue derived FABPs of cattle showed FABP-heart closest and FABP-epidermis most distantly placed from buffalo MDGI. This report also differs from the earlier findings that MDGI is intermediate of FABP-heart and adipose.     

Isolation and Sequence Characterization of Mammary Derived Growth Inhibitor Gene of Riverine Buffalo (Bubalus Bubalis)

In this study, attempts have been made to identify and characterize water buffalo (Bubalus bubalis) mammary derived growth inhibitor (MDGI) gene, isolated from a mammary gland cDNA library of lactating buffalo. The complete MDGI cDNA was of 698 nucleotides, consisting 61 nucleotides in 5′ UTR, coding region of 402 nucleotides, and 235 nucleotides representing the 3′ UTR. Comparison of nucleotide and deduced amino acid sequence data with that of MDGI//fatty acid binding protein (FABP) of other species shows three buffalo specific nucleotide changes while seven nucleotide changes were common to cattle and buffalo. Buffalo and cattle MDGI had 100% amino acid sequence similarity, which also shared three amino acid changes: 34 (Ala-Gly), 109 (Leu-Met), and 132 (Glu-Gln) as compared to other species. Comparison with FABPs reported from other cattle tissues revealed highest amino acid sequence similarity with FABP-heart (100%) and least with FABP-liver (20.5%). Phylogenetic analysis revealed cattle MDGI to be closest to buffalo, while mouse MDGI was distantly placed, whereas different tissue derived FABPs of cattle showed FABP-heart closest and FABP-epidermis most distantly placed from buffalo MDGI. This report also differs from the earlier findings that MDGI is intermediate of FABP-heart and adipose.     

Analysis of conserved microsatellite sequences suggests closer relationship between water buffalo Bubalus bubalis and sheep Ovis aries.

The distribution and evolutionary pattern of the conserved microsatellite repeat sequences (CA)n, (TGG)6, and (GGAT)4 were studied to determine the divergence time and phylogenetic position of the water buffalo, Bubalus bubalis. The mean allelic frequencies of these repeat loci showed a high level of heterozygosity among the euartiodactyls (buffalo, cattle, sheep, and goat). Genetic distances calculated from the allelic frequencies of these microsatellites were used to position Bubalus bubalis in the phylogenetic tree. The tree topology revealed a closer proximity of the Bubalus bubalis to the Ovis aries (sheep) genome than to other domestic species. The estimated time of divergence of the water buffalo genome relative to cattle, goat, sheep, pig, rabbit, and horse was found to be 21, 0.5, 0.7, 94, 20.3, and 408 million years (Myr), respectively. Although water buffaloes share morphological and biochemical similarities with cattle, our study using the microsatellite sequences places the bubaline species in an entirely new phylogenetic position. Our results also suggest that with respect to these repeat loci, the water buffalo genome shares a common ancestry with sheep and goat after the divergence of subfamily Bovinae (Bos taurus) from the family Bovidae.     

Identification of Cu/Zn superoxide dismutase in cattle and river buffaloes

All air-living organisms produce superoxide dismutase (SOD) and several antioxidant enzymes that limit oxidative stress by detoxifying reactive oxygen species. SOD1 gene has been investigated in Egyptian native cattle and buffalos at the level of genomic DNA and cDNA that were extracted from leucocytes. An unexpected band at approximately 370 bp was obtained in cattle genomic DNA and cDNA as well as in buffalo cDNA. SOD1 amplified sequence of native cattle genomic DNA and cDNA showed a 93% alignment. Native cattle genomic DNA SOD1 amplicon shares sequence homology with mRNAs of Bos taurus “similar to superoxide dismutase” (SOD1) sequence of the GeneBank database. It also shares sequence homology with “similar to superoxide dismutase” on B. taurus chromosome BTA13. The results indicate that the genomic DNA of Egyptian native cattle contains SOD1 processed pseudo gene. SOD1 primers amplified three fragments in buffalo genomic DNA which indicates that buffalo genome has different copies of SOD1 due to alternative splicing. It failed to produce the 370 bp fragments found in cattle DNA. The protein analysis revealed no differences between Egyptian native cattle and B. taurus SOD1 mRNA. However, one amino acid, aspartic acid (Asp), in Egyptian native cattle and B. taurus SOD1, is substituted with asparagine (Asn) (D26N) in buffaloes. This amino acid substitution may be due to non-synonymous single nucleotide polymorphisms (nsSNPs). The nsSNPs detected in buffaloes may affect the function of the encoded protein. This study is the first investigation reporting that the resistance of the buffalo to diseases and parasites that afflict cattle may not be acquired but may have a genetic basis.     

牛科动物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)更为合理。     

Molecular cloning and characterization of beta-defensin cDNA expressed in distal ileum of buffalo (Bubalus bubalis).

Defensins play a prominent role in protection of various epithelial surfaces. In this study, we have cloned and characterized the mRNA from the distal ileum of Bubalus bubalis. Total RNA after isolation from ileal epithelium was reverse transcribed to synthesize cDNA using primers designed by taking conserved region of cattle enteric beta-defensin (EBD) mRNA, goat beta-defensin 2 (BD 2) and cattle lingual antimicrobial peptide (LAP) mRNA sequences. The PCR amplified cDNA of 254 bp was ligated to pDrive cloning vector and transformed into XL-blue strain of E coli. The sequence analysis indicated 29 nucleotide substitutions with reported cattle EBD mRNA sequence sharing 86.2% homology, 92.1% with cattle LAP, 81.6% with cattle tracheal antimicrobial peptide and 84.6% with goat BD 2. The deduced amino acid sequence encodes for a 64 amino acid precursor peptide. Both nucleotide and amino acid sequence homology shows that the cloned sequence is closer to cattle LAP.     

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.     

内蒙古白绒山羊4E-BP1基因cDNA克隆及组织表达特异性分析

旨在克隆内蒙古白绒山羊4E-BP1(真核细胞翻译起始因子4E结合蛋白1)基因并进行生物信息学及表达模式分析。根据已报道物种4E-BP1基因cDNA序列,用primer premier5软件设计引物,通过RT-PCR从绒山羊胎儿成纤维细胞总RNA中扩增出4E-BP1基因编码区cDNA序列,对目的片段进行测序及表达模式分析。克隆到的内蒙古白绒山羊4E-BP1基因cDNA全长357 bp,包含了完整的的ORF,编码118个氨基酸残基。核酸序列与牛、马、人、大鼠及小鼠的同源性分别为98%、90%、90%、88%和87%。4E-BP1基因在绒山羊脑、心脏、睾丸及胰腺组织中均有表达。     

Polymorphism and gene organization of water buffalo MHC-DQB genes show homology to the BoLA DQB region

In cattle (Bos taurus), there is evidence of more than 50 alleles of BoLA-DQB (bovine lymphocyte antigen DQB) that are distributed across at least five DQB loci, making this region one of the most complex in the BoLA gene family. In this study, DQB alleles were analysed for the water buffalo (Bubalus bubalis), another economically important bovine species. Twelve alleles for Bubu-DQB (Bubalis bubalis DQB) were determined by nucleotide sequence analysis. A phylogenetic analysis revealed numerous trans-species polymorphisms, with alleles from water buffalo assigned to at least three different loci (BoLA-DQB1, BoLA-DQB3 and BoLA-DQB4) that are also found in cattle. These presumptive loci were analysed for patterns of synonymous (d(S)) and non-synonymous (d(N)) substitution. Like BoLA-DQB1, Bubu-DQB1 was observed to be under strong positive selection for polymorphism. We conclude that water buffalo and cattle share the current arrangement of their DQB region because of their common ancestry.     

Sequence of specific mitochondrial 16S rRNA gene fragment from Egyptian buffalo is used as a pattern for discrimination between river buffaloes,cattle, sheep and goats

Characterization of molecular markers and the development of better assays for precise and rapid detection of domestic species are always in demand. This is particularly due to recent food scares and the crisis of biodiversity resulting from the huge ongoing illegal traffic of endangered species. The aim of this study was to develop a new and easy method for domestic species identification (river buffalo, cattle, sheep and goat) based on the analysis of a specific mitochondrial nucleotide sequence. For this reason, a specific fragment of Egyptian buffalo mitochondrial 16S rRNA gene (422 bp) was amplified by PCR using two universal primers. The sequence of this specific fragment is completely conserved between all tested Egyptian buffaloes and other river buffaloes in different places in the world. Also, the lengths of the homologous fragments were less by one nucleotide (421 bp) in case of goats and two nucleotides (420 bp) in case of both cattle and sheep. The detection of specific variable sites between investigated species within this fragment was sufficient to identify the biological origin of the samples. This was achieved by alignment between the unknown homologous sequence and the reference sequences deposited in GenBank database (accession numbers, FJ748599–FJ748607). Considering multiple alignment results between 16S rRNA homologous sequences obtained from GenBank database with the reference sequence, it was shown that definite nucleotides are specific for each of the four studied species of the family Bovidae. In addition, other nucleotides are detected which can allow discrimination between two groups of animals belonging to two subfamilies of family Bovidae, Group one (closely related species like cattle and buffalo, Subfamily Bovinae) and Group two (closely related species like sheep and goat, Subfamily Caprinae). This 16S DNA barcode character-based approach could be used to complement cytochrome c oxidase I (COI) in DNA barcoding. Also, it is a good tool for identification of unknown sample belonging to one of the four domestic animal species of family Bovidae quickly and easily.