In silico mining of putative microsatellite markers from whole genome sequence of water buffalo (Bubalus bubalis) and development of first BuffSatDB
© Sarika et al.; licensee BioMed Central Ltd. 2013
Received: 16 June 2012
Accepted: 9 January 2013
Published: 19 January 2013
Though India has sequenced water buffalo genome but its draft assembly is based on cattle genome BTau 4.0, thus de novo chromosome wise assembly is a major pending issue for global community. The existing radiation hybrid of buffalo and these reported STR can be used further in final gap plugging and “finishing” expected in de novo genome assembly. QTL and gene mapping needs mining of putative STR from buffalo genome at equal interval on each and every chromosome. Such markers have potential role in improvement of desirable characteristics, such as high milk yields, resistance to diseases, high growth rate. The STR mining from whole genome and development of user friendly database is yet to be done to reap the benefit of whole genome sequence.
By in silico microsatellite mining of whole genome, we have developed first STR database of water buffalo, BuffSatDb (Buffalo MicroSatellite Database (http://cabindb.iasri.res.in/buffsatdb/) which is a web based relational database of 910529 microsatellite markers, developed using PHP and MySQL database. Microsatellite markers have been generated using MIcroSAtellite tool. It is simple and systematic web based search for customised retrieval of chromosome wise and genome-wide microsatellites. Search has been enabled based on chromosomes, motif type (mono-hexa), repeat motif and repeat kind (simple and composite). The search may be customised by limiting location of STR on chromosome as well as number of markers in that range. This is a novel approach and not been implemented in any of the existing marker database. This database has been further appended with Primer3 for primer designing of the selected markers enabling researcher to select markers of choice at desired interval over the chromosome. The unique add-on of degenerate bases further helps in resolving presence of degenerate bases in current buffalo assembly.
Being first buffalo STR database in the world , this would not only pave the way in resolving current assembly problem but shall be of immense use for global community in QTL/gene mapping critically required to increase knowledge in the endeavour to increase buffalo productivity, especially for third world country where rural economy is significantly dependent on buffalo productivity.
Keywordsde novo Microsatellites Primers Radiation hybrid Water buffalo
Water buffalo (Bubalus bubalis) contributes immensely to the agricultural economy of Indian subcontinent, South East Asian countries through milk, meat, hides, fertilizer, fuel and draught animal power. A large part of human population depends on this species than any other livestock species in the world . There is 188.3 million buffalo population in the world which contributes around 55 – 60% of total milk production . Asia has nearly 97% of buffaloes and is an integral part of agriculture in India, China, Pakistan, Nepal, Bangladesh, Thailand, Myanmar and Malaysia. The productivity of buffaloes in these regions is higher as compared to cattle .
Molecular markers can play a significant role for livestock improvement through conventional breeding strategies. Scientific resources are limited in many of the countries where buffaloes are economically important livestock and as a consequence, genome research has not been supported at the level of some of the other species . Limited number of researches has been conducted globally exploring the genetic diversity on molecular genetic basis in buffalo in comparison with other farm animal genetic resources. This depends, in part on the knowledge of their genetic structure based on molecular markers like microsatellites .
Microsatellites are sequences made up of a simple sequence motif, not more than six bases long, that is tandemly repeated and arranged head to tail without interruption by any other base or motif. Simple, tandemly repeated di- and tri- nucleotide sequences have been demonstrated to be polymorphic in length in a number of eukaryotic genome . The frequency with which they occur (once every 50,000–60,000 bp), the high degree of polymorphism displayed, and their random distribution across the genome  make them potentially very useful as DNA markers in gene mapping studies. Furthermore, two or more microsatellites may be analyzed simultaneously [8, 9], opening new opportunity for genetic analysis of large number of samples.
To cater the need of microsatellite especially for biodiversity analysis, cattle microsatellite markers have been used in heterologous mode in buffalo and up to 56% of them have been found polymorphic . Cattle microsatellite markers have many disadvantages in such diversity analysis like low polymorphism and loss of amplification due to null alleles, size biasness, hitch hiking and potential exclusion of abundant STR in gene pool . Even there is limited work of STR mining using partial enriched genomic library .
There is no thorough in silico STR marker mining from buffalo genome to represent more holistic and cumulative variability of genome to be used in gene pool or biodiversity analysis and gene/QTL mapping. Though India has sequenced water buffalo genome but its draft assembly is based on cattle genome BTau 4.0, thus de novo chromosome wise assembly is a major pending issue for global community . The existing radiation hybrid of buffalo by Amaral et al. and these reported STR can be used further in final gap plugging and “finishing” expected in de novo genome assembly. Such work needs extensive STR mining from buffalo genome at equal interval on each and every chromosome. In order to cater this urgent need in resolving assembly, mapping issues and biodiversity analysis, we have developed first STR database of water buffalo, BuffSatDb (Buffalo MicroSatellite Database) which is a web based relational database of microsatellites.
Construction and content
Data collection and architecture
BuffSatDb has eight tabs (Home, About, Database, Analysis, Tutorial, Links, Contact, Team). General information of the developed microsatellite database, information about Water buffalo, microsatellite markers, comparative analysis of the buffalo genome has been discussed. The tutorial of this database contains the guidelines for users and terminologies used in the database contents. BuffSatDb is appended with other useful links, the team and contact persons.
In silico mining of microsatellite from whole genome of water buffalo
The user can query for microsatellites, chromosome wise (1–24, X, M and U), where more than one chromosome may be selected at a time from water buffalo genome. These searches may further be customised based on microsatellite characteristics like motif type (mono, di, tri, tetra, penta, hexa), repeat motif and repeat kind (simple and composite). The user may further go for advance search like limiting the location on chromosome as well as the number of markers in that range. This is a novel approach and to the best of our knowledge, it has not been implemented in any of the existing marker database which may be useful for the researchers. Identification of QTL and fine mapping of economically important genes based on LOD (Logarithm of the Odds) score also needs STR preferably at equal interval. Also other parameters like GC content, range of STR location and copy number may be customised for the above selection according to the requirement of researchers. The results are then displayed in tabular format, giving chromosome number, motif type, motif, copy number, basepair, start and end position along with the GC content.
BuffSatDb is further appended with Primer3 tool . The STRs traced by the query, may be selected with the help of radiobutton for generation of primers. Primer for selected STR locus may be designed with a template of approximately 1000 base pairs by selecting upto 500 base pairs of both flanking regions. These flexibilities would enable researchers to select markers of choice at desired interval over the chromosomes. Further one can use each individual STR of a targeted region over chromosome to narrow down location of gene of interest or linked QTL. A novel add-on for degenerate bases has been incorporated in this database search, where the users are given flexibility to replace degenerate bases with any of the alternative bases (A,T,G,C). This feature has been added to resolve the issue of some of the degenerate bases present in current buffalo genome assembly making the primer designing very difficult otherwise.
Frequencies of STRs based on their sizes
Size of STRs
Number of STRs
Contribution in percentage
Chromosome wise distribution of STRs
STRs validation result of homologous and heterologous primer pairs of water buffalo
ISAG–FAO recommended STRs from cattle
ISAG–FAO recommended STRs from buffalo
Nagarajan et al, monomorphic STRs from buffalo
Nagarajan et al, polymorphic STR loci from buffalo
Total no. of primer pairs reported
No. of positive primers (Forward)
No. of positive primers (Reverse)
No. of positive primers (common to both forward and reverse)
Discussion and utility
A total of 910529 microsatellite markers have been searched by in silico mining. Simple STR were found to be most abundant (91.16%). Microsatellite density has been found positively correlated with genome size [18–20]. Among fully sequenced eukaryotic genomes, microsatellite density is highest in mammals. However in case of plant, microsatellite frequency is negatively correlated with genome size .
In the present study of water buffalo, mono- motif was found to be most abundant. Relative distributions of different microsatellite motif length classes in genomes differ considerably from species to species .
In case of water buffalo, it was found that longer repeats are less in abundance which is expected as reported and described in various studies [23, 24]. It was also observed that microsatellite size range is increasing from 10 up to 14–16, however beyond this size range, it again starts decreasing. This is due to cyclical nature of microsatellite marker per say in its course of evolution. The birth of microsatellite starts with, out of register loop in event of DNA replication with a threshold size of 8 repeat unit or more, in the form of simple repeat. Gradually due to background mutation simple repeat gets converted in compound repeat. At the stage of simple repeat, the rate of mutation is high and predominantly it is addition of repeat unit and hence size increases. But once background mutation converts simple repeat into compound interrupted repeat, the smaller size simple repeat of less than 8 unit gets pinched off in subsequent replications. This maintains the size of microsatellite as evolutionary constraints otherwise microsatellite marker would have been always increasing in length during course of evolution. Thus individual microsatellites arrays have a “life cycle” of sorts, they are born, they grow and ultimately they perish. These events may stretch over tens or even hundreds of millions of years [25, 26].
Water buffalo microsatellite profile exhibits the similar pattern. The relative abundance of repeat motif were in order of mono, di, tri, penta, tetra and hexa (Table 2). Though di-nucleotide repeats are most abundant in eukaryotic genome [27, 28] but we found most abundance of mononucleotide repeats across all chromosome. This relatively higher abundance of mono over di nucleotide repeat type might be due to inherent limitation of the NGS technology which adds more mono nucleotide causing sequencing error . The longer the chromosome proportionately higher the total repeat content as expected in ubiquitously distributed STR markers .
In order to validate the previously reported STR markers, two sets viz. heterologous (cattle original species and buffalo focal species), homologous STR (developed from buffalo and validated in buffalo) were considered. The heterologous markers recommended by FAO-IASG  and homologous marker  were used. It was observed that both subsets of heterologous ISAG-FAO recommended primer for cattle and buffalo diversity analysis gave less validation results i.e. 10% and 13.33% respectively. Cross species amplifiability is due to conservation of cattle STR and its flanking regions in other species . Though some of the primers showed validation up to 36.67% (Table 3). In the cross species amplifiability of bovidae species, such data are usually expected due to null alleles and genomic changes during speciation . In validation of homologous STR, it was found that both subsets reported higher percentage of monomorphic (28.57%) and polymorphic (24.30%) loci. The validation results are limited as the first draft genome assembly of buffalo is based on cattle and it is not completely finished.
The findings of this study has limitations which need to be addressed. As genome of water buffalo is just draft assembly based on cow assembly Btau 4.0, thus de novo assembly is needed to have the buffalo specific chromosome wise microsatellite profile. The current database is based on chromosome number of cattle which is certainly not the same in case of buffalo. For example cattle chromosome 4 is actually buffalo chromosome 8. In fact only chromosome number common between cattle and buffalo are just 5 viz 1, 2, 17, 18, and X . The splitting and translocation has rendered syntenic relationship between these two species which are well documented. Nevertheless the microsatellites in our database with option of primer designing at desired place over “chromosome” will be of immense use especially over radiation hybrid of buffalo to resolve the problem and current issue of de novo assembly. Besides this, these markers can be further used for QTL, gene mapping as well as biodiversity analysis in setting the conservation priorities. The markers present in our database need further wet lab validation. Being first database of water buffalo microsatellite especially at juncture where de novo genome assembly is yet to be done, the use of these markers are highly warranted in order to “finishing” of water buffalo genome assembly. This will further lead to next version of buffalo microsatellite database base with proper buffalo specific chromosome wise data which is hitherto missing but critically needed. Such endeavour will fetch not only increase in buffalo productivity but also greater food security especially in third and new world countries.
Being first buffalo STR database in the world, this would not only pave the way in resolving current water buffalo genome assembly problem but shall be of immense use for global community in QTL/gene mapping critically required to increase knowledge in the endeavour to increase buffalo productivity, especially for third world country where rural economy is significantly dependent on buffalo productivity.
Availability and requirement
BuffSatDb, the buffalo microsatellite marker database is freely accessible for research purposes for non-profit and academic organizations at http://cabindb.iasri.res.in/buffsatdb/.
Financial assistance received by National Agricultural Innovation Project, Indian Council of Agricultural Research, New Delhi entitled “Establishment of National Agriculture Bioinformatics Grid in ICAR” is gratefully acknowledged. The technical assistance of Jai Bhagwan in maintaining the web server and AR Paul in designing the logo of BuffSatDb are thankfully acknowledged. Authors acknowledge the critical input of all the three anonymous reviewers and the editor in improvement of the manuscript.
- Scherf BD: World Watch List. From World watch list for domestic animal diversity. 2000, Rome, Italy: FAOGoogle Scholar
- FAO: FAO Statistical Yearbook 2009. 2009, Rome, (http://www.fao.org/economic/the-statistics-division-ess/publications-studies/statistical-yearbook/fao-statistical-yearbook-2009/en/)Google Scholar
- Tantia MS, Vijh RK, Bhasin V, Sikka P, Vij PK, Kataria RS, Mishra BP, Yadav SP, Pandey AK, Sethi RK, Joshi BK, Gupta SC, Pathak KML: Whole-genome sequence assembly of the water buffalo (Bubalus bubalis). Indian J Anim Sci. 2011, 81 (5): 38-Google Scholar
- Womack JE: Advances in livestock genomics: opening the barn door. Genome Res. 2005, 15: 1699-1705. 10.1101/gr.3809105.View ArticlePubMedGoogle Scholar
- Gonçalves EC, Silva A, Barbosa MSR, Schneider MPC: Isolation and characterization of microsatellite loci in Amazonian red-handed howlers Alouatta belzebul (Primates, Plathyrrini). Mol Ecol. 2004, 4: 406-408. 10.1111/j.1471-8286.2004.00667.x.View ArticleGoogle Scholar
- Litt M, Luty JA: A hypervariable microsatellite revealed by in vitro amplification of a dinucleotide repeat within the cardiac muscle actin gene. Am J Hum Genet. 1989, 44: 397-401.PubMed CentralGoogle Scholar
- Luty JA, Guo Z, Willard HF, Ledbetter DH, Ledbetter S, Litt M: Five polymorphic microsatellite VNTRs on the human X chromosome. Am J Hum Genet. 1990, 46: 776-783.PubMed CentralPubMedGoogle Scholar
- Weber JL, May PE: Abundant class of human DNA polymorphisms which can be typed using the polymerase chain reaction. Am J Hum Genet. 1989, 44: 388-396.PubMed CentralPubMedGoogle Scholar
- Georges M, Mishra A, Sargeant L, Steele M, Zhao X: Progress towards a primary DNA marker map in cattle. 4th World Congress Genetics Applied Livestock Production. 1990, 13: 107-112.Google Scholar
- Navani N, Jain PK, Gupta S, Sisodia BS, Kumar S: A set of cattle microsatellite DNA markers for genome analysis of riverine buffalo (Bubalus bubalis). Anim Genet. 2002, 33: 149-154. 10.1046/j.1365-2052.2002.00823.x.View ArticlePubMedGoogle Scholar
- Ellegren H, Primmer CR, Sheldon BC: Microsatellite evolution-directionality or bias. Nat Genet. 1995, 11: 360-362. 10.1038/ng1295-360.View ArticlePubMedGoogle Scholar
- Venancio LPR, Amaral MEJ, Schneider MPC, Silva AL: Microsatellite loci isolation from river buffalo using enriched partial genomic libraries. Ital J Anim Sci. 2007, 6 (2): 338-341.Google Scholar
- Amaral MEJ, Grant JR, Riggs PK, Stafuzza NB, Filho EAR, Goldammer T, Weikard R, Brunner RM, Kochan KJ, Greco AJ, Jeong J, Cai Z, Lin G, Prasad A, Kumar S, Saradhi GP, Mathew B, Kumar MA, Miziara MN, Mariani P, Caetano AR, Galvão SR, Tantia MS, Vijh RK, Mishra B, Kumar STB, Pelai VAP, Santana AM, Fornitano LC, Jones BC, Tonhati H, et al: A first generation whole genome RH map of the river buffalo with comparison to domestic cattle. BMC Genomics. 2008, 9: 631-10.1186/1471-2164-9-631.PubMed CentralView ArticlePubMedGoogle Scholar
- Thiel T, Michalek W, Varshney RK, Graner A: Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003, 106: 411-422.PubMedGoogle Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, NJ: Humana Press, 365-386.Google Scholar
- FAO: Molecular genetic characterization of animal genetic resources. 2011, Rome: FAO Animal Production and Health Guidelines. No. 9Google Scholar
- Nagarajan MN, Nishanth G, Haribaskar R, Paranthaman K, Gupta J, Mishra M, Vaidhegi R, Kumar S, Ranjan AK, Kumar S: Microsatellite markers of water buffalo, Bubalus bubalis -development, characterisation and linkage disequilibrium studies. BMC Genet. 2009, 10: 68-74.PubMed CentralView ArticlePubMedGoogle Scholar
- Hancock JM: Simple sequences and the expanding genome. BioEssays. 1996, 18: 421-425. 10.1002/bies.950180512.View ArticlePubMedGoogle Scholar
- Toth G, Gáspári Z, Jurka J: Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res. 2000, 10 (7): 967-981. 10.1101/gr.10.7.967.PubMed CentralView ArticlePubMedGoogle Scholar
- Katti MV, Ranjekar PK, Gupta VS: Differential distribution of simple sequence repeats in eukaryotic genome sequences. Mol Biol Evol. 2001, 18: 1161-1167. 10.1093/oxfordjournals.molbev.a003903.View ArticleGoogle Scholar
- Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nature Genet. 2002, 30: 194-200. 10.1038/ng822.View ArticlePubMedGoogle Scholar
- Sharma PC, Grover A, Kahl G: Mining microsatellites in eukaryotic genomes. Trends Biotechnol. 2007, 25 (11): 490-498. 10.1016/j.tibtech.2007.07.013.View ArticlePubMedGoogle Scholar
- Temnykh S, Declercke G, Lukashova A, Lipovich L, Cartinhour S, McCouch SR: Computational and experimental analysis of microsatellite in rice (O. sativa L) frequency, length variation, transposon associations and genetic marker potential. Genome Res. 2001, 11: 1441-1452. 10.1101/gr.184001.PubMed CentralView ArticlePubMedGoogle Scholar
- Grover A, Aishwarya V, Sharma PC: Biased distribution of microsatellite in rice genome. Mol Genet Genomics. 2007, 277: 469-480. 10.1007/s00438-006-0204-y.View ArticlePubMedGoogle Scholar
- Messier W, Li SH, Stewart CB: The birth of microsatellites. Nature. 1996, 381: 483-View ArticlePubMedGoogle Scholar
- Primmer CR, Ellegren H: Patterns of molecular evolution in avian microsatellites. Mol Biol Evol. 1998, 15: 997-1008. 10.1093/oxfordjournals.molbev.a026015.View ArticlePubMedGoogle Scholar
- Kariin S, Burge C: Dinucleotide relative abundance extremes: a genomic signature. Trends Genet. 1995, 11 (7): 283-290. 10.1016/S0168-9525(00)89076-9.View ArticleGoogle Scholar
- Shioiri C, Takahata N: Skew of mononucleotide frequencies, relative abundance of dinucleotides and DNA strand asymmetry. J Mol Evol. 2001, 53: 364-376. 10.1007/s002390010226.View ArticlePubMedGoogle Scholar
- Haseneyer G, Schmutzer T, Seidel M: From RNA-seq to large-scale genotyping-genomics resources for rye (Secale cereale L.). BMC Plant Biol. 2011, 11: 131-143. 10.1186/1471-2229-11-131.PubMed CentralView ArticlePubMedGoogle Scholar
- Zietkiewicz E, Rafalski A, Labuda D: Genome finger priting by Simple Sequence Repeat (SSR)-Anchored Polymerase Chain Reaction Amplification. Genomics. 1994, 20: 176-183. 10.1006/geno.1994.1151.View ArticlePubMedGoogle Scholar
- Nguyen TT, Genini S, Bui LC, Voegeli P, Stranzinger G, Renard JP, Maillard JC, Nguyen BX: Genomic conservation of cattle microsatellite loci in wild gaur (Bos gaurus) and current genetic status of this species in Vietnam. BMC Genet. 2007, 8: 77-84.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim KS, Min MS, An JH, Lee H: Cross-species amplification of bovidae microsatellites and Low diversity of the endangered Korean goral. J Hered. 2004, 95 (6): 521-525. 10.1093/jhered/esh082.View ArticlePubMedGoogle Scholar