Microsatellite markers of water buffalo, Bubalus bubalis--development, characterisation and linkage disequilibrium studies

Deciphering genetic diversity in conserved cattle bulls to achieve sustainable development goals

The primary objective of Sustainable Development Goal target 2.5 established by the United Nations is to ensure the preservation of genetic diversity in domesticated animals. The ICAR-National Bureau of Animal Genetic Resources in India has been actively engaged in the conservation of cattle and buffalo bull semen for long-term storage. This present study aimed to assess the genetic diversity present in the conserved cattle bull semen, which would aid in determining the most suitable strategy for future conservation management. A total of 192 bull semen belonging to 19 cattle breeds were selected to evaluate genetic diversity using 17 pairs of FAO recommended microsatellite primers. Total 267 alleles were detected across all the samples which indicates substantial amount of allelic variation is being maintained in conserved bulls. Further, all cattle bulls semen conserved showed higher observed heterozygosity than expected heterozygosity which indicates excess genetic diversity in all the populations. The FST, F IT and FIS value across the loci and population is 0.146 ± 0.009, 0.054 ± 0.038, and - 0.105 ± 0.035, respectively, which suggests lack of inbreeding in conserved cattle bull semen. This study has established genetic diversity in conserved cattle semen samples to achieve sustainable development goals. In addition, it provides compelling evidence that the current approach for conserving cattle bull semen is heading in the correct direction.

Microsatellite Diversity and Phylogenetic Relationships among East Eurasian Bos taurus Breeds with an Emphasis on Rare and Ancient Local Cattle

We report the genetic analysis of 18 population samples of animals, which were taken from cattle (Bos taurus) breeds of European and Asian origins. The main strength of our study is the use of rare and ancient native cattle breeds: the Altai, Ukrainian Grey, Tagil, and Buryat ones. The cattle samples studied have different production purposes, belong to various eco-geographic regions, and consequently have distinct farming conditions. In order to clarify the genetic diversity, phylogenetic relationships and historical origin of the studied breeds, we carried out an analysis of the genetic variation of 14 high-variability microsatellite loci at 1168 genotyped animals. High levels of heterozygosity and allelic richness were identified in four of the ancient local breeds, namely the Kalmyk, Tagil, Kyrgyz native, and Buryat breeds. The greatest phylogenetic distances from a common ancestor were observed for the Yakut and Ukrainian Grey breeds, while the Tagil breed showed the smallest difference. By using clustering approaches, we found that the Altai cattle is genetically close to the Kyrgyz one. Moreover, both the Altai and Kyrgyz breeds exposed genetic divergences from other representatives of the Turano-Mongolian type and genetic relationships with the Brown Swiss and Kostroma breeds. This phenomenon can be explained by the extensive use of the Brown Swiss and Kostroma breeds in the breeding and improvement processes for the Kyrgyz breeds, which have been involved in the process of keeping the Altai cattle. Our results can be valuable for conservation and management purposes.

Study of whole genome linkage disequilibrium patterns of Iranian water buffalo breeds using the Axiom Buffalo Genotyping 90K Array

Accuracy of genome-wide association studies, and the successful implementation of genomic selection depends on the level of linkage disequilibrium (LD) across the genome and also the persistence of LD phase between populations. In the present study LD between adjacent SNPs and LD decay between SNPs was calculated in three Iranian water buffalo populations. Persistence of LD phase was evaluated across these populations and effective population size (Ne) was estimated from corrected r² information. A set of 404 individuals from three Iranian buffalo populations were genotyped with the Axiom Buffalo Genotyping 90K Array. Average r² and |D'| between adjacent SNP pairs across all chromosomes was 0.27 and 0.66 for AZI, 0.29 and 0.68 for KHU, and 0.32 and 0.72 for MAZ. The LD between the SNPs decreased with increasing physical distance from 100Kb to 1Mb between markers, from 0.234 to 0.018 for AZI, 0.254 to 0.034 for KHU, and 0.297 to 0.119 for MAZ, respectively. These results indicate that a density of 90K SNP is sufficient for genomic analyses relying on long range LD (e.g. GWAS and genomic selection). The persistence of LD phase decreased with increasing marker distances across all the populations, but remained above 0.8 for AZI and KHU for marker distances up to 100Kb. For multi-breed genomic evaluation, the 90K SNP panel is suitable for AZI and KHU buffalo breeds. Estimated effective population sizes for AZI, KHU and MAZ were 477, 212 and 32, respectively, for recent generations. The estimated effective population sizes indicate that the MAZ is at risk and requires careful management.

Analysis of microsatellite markers in a Cuban water buffalo breed

The aim of this Regional Research Communication was to validate a panel of 30 microsatellite markers recommended by FAO/ISAG for studies of biodiversity in cattle to improve the characterisation of Cuban buffalo populations. The water buffalo (Bubalus bubalis) is an economically important livestock species. Therefore, research focused on the study of the genetic relationships among water buffalo populations is useful to support conservation decisions and to design breeding schemes. Twenty-eight of the 30 tested regions were amplified, one of which (ETH10) turned out to be monomorphic. A total of 143 alleles were observed in the Cuban water buffalo population. The average number of alleles per locus was 5·04. The number of alleles per polymorphic locus ranged from two (INRA 63 and MM12) to nine (ETH185). The observed and expected heterozygosity ranged from 0·108 (HAUT24) to 0·851 (CSSM66) and 0·104 (MM12) to 0·829(INRA32), respectively. The polymorphic information content (PIC) ranged from 0·097 (MM12) to 0·806 (INRA32), and the overall value for these markers was 0·482. Within the population, inbreeding estimates (FIS) was positive in 14 of the 30 loci analysed. This study thus highlights the usefulness of heterologous bovine microsatellite markers to assess the genetic variability in Cuban water buffalo breeds. Furthermore, the results can be utilised for future breeding strategies and conservation.

In silico mining of putative microsatellite markers from whole genome sequence of water buffalo (Bubalus bubalis) and development of first BuffSatDB

BACKGROUND

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.

DESCRIPTION

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.

CONCLUSION

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.

Power of exclusion of 19 microsatellite markers for parentage testing in river buffalo (Bubalus bubalis)

In the present study, 19 microsatellite markers were assessed for their power of exclusion to test parentage in river buffalo. Microsatellite genotypes of 216 unrelated buffaloes belonging to five different breeds were utilized for the study. The probabilities of exclusion were calculated for three hypothetical situations viz. paternity testing (PE1), one parental genotype unavailable (PE2) and exclusion of both parents i.e. substituted offspring (PE3). The mean probability of exclusion across 19 investigated markers in buffalo was 0.578 (PE1), 0.405 (PE2) and 0.764 (PE3) respectively. The probability of exclusion for paternity (PE1) ranged between 0.297 and 0.814 across different markers. The exclusion probability for the cases one parent unavailable (PE2) and substituted offspring (PE3) varied from 0.143 to 0.688 and 0.465 to 0.946 respectively. Polymorphism information content and expected heterozygosity were found to have significantly high correlation with probability of exclusion of microsatellite markers. The cumulative PE1 of nine marker loci was estimated to be 0.9999 while in case of absence of one of the parental genotypes, a minimum of 11 markers were required to achieve a cumulative PE2 of 0.999. In conclusion, the present study proposes two multiplex sets with four and five markers respectively for routine parentage testing in buffalo and an additional set of four markers for doubtful cases of paternity.