The domestic water buffalo is native to the Asian continent but through historical migrations and recent importations, nowadays has a worldwide distribution. The two types of water buffalo, i.e., river and swamp, display distinct morphological and behavioral traits, different karyotypes and also have different purposes and geographical distributions. River buffaloes from Pakistan, Iran, Turkey, Egypt, Romania, Bulgaria, Italy, Mozambique, Brazil and Colombia, and swamp buffaloes from China, Thailand, Philippines, Indonesia and Brazil were genotyped with a species-specific medium-density 90K SNP panel. We estimated the levels of molecular diversity and described population structure, which revealed historical relationships between populations and migration events. Three distinct gene pools were identified in pure river as well as in pure swamp buffalo populations. Genomic admixture was seen in the Philippines and in Brazil, resulting from importations of animals for breed improvement. Our results were largely consistent with previous archeological, historical and molecular-based evidence for two independent domestication events for river- and swamp-type buffaloes, which occurred in the Indo-Pakistani region and close to the China/Indochina border, respectively. Based on a geographical analysis of the distribution of diversity, our evidence also indicated that the water buffalo spread out of the domestication centers followed two major divergent migration directions: river buffaloes migrated west from the Indian sub-continent while swamp buffaloes migrated from northern Indochina via an east-south-eastern route. These data suggest that the current distribution of water buffalo diversity has been shaped by the combined effects of multiple migration events occurred at different stages of the post-domestication history of the species.

BackgroundThe availability of the bovine genome sequence and SNP panels has improved various genomic analyses, from exploring genetic diversity to aiding genetic selection. However, few of the SNP on the bovine chips are polymorphic in buffalo, therefore a panel of single nucleotide DNA markers exclusive for buffalo was necessary for molecular genetic analyses and to develop genomic selection approaches for water buffalo. The creation of a 90K SNP panel for river buffalo and testing in a genome wide association study for milk production is described here.MethodsThe genomes of 73 buffaloes of 4 different breeds were sequenced and aligned against the bovine genome, which facilitated the identification of 22 million of sequence variants among the buffalo genomes. Based on frequencies of variants within and among buffalo breeds, and their distribution across the genome, inferred from the bovine genome sequence, 90,000 putative single nucleotide polymorphisms were selected to create an Axiom® Buffalo Genotyping Array 90K.ResultsThis 90K “SNP-Chip” was tested in several river buffalo populations and found to have ∼70% high quality and polymorphic SNPs. Of the 90K SNPs about 24K were also found to be polymorphic in swamp buffalo. The SNP chip was used to investigate the structure of buffalo populations, and could distinguish buffalo from different farms. A Genome Wide Association Study identified genomic regions on 5 chromosomes putatively involved in milk production.ConclusionThe 90K buffalo SNP chip described here is suitable for the analysis of the genomes of river buffalo breeds, and could be used for genetic diversity studies and potentially as a starting point for genome-assisted selection programmes. This SNP Chip could also be used to analyse swamp buffalo, but many loci are not informative and creation of a revised SNP set specific for swamp buffalo would be advised.

No abstract

Water buffalo is a globally important species for agriculture and local economies. A de novo assembled, well-annotated reference sequence for the water buffalo is an important prerequisite for studying the biology of this species, and is necessary to manage genetic diversity and to use modern breeding and genomic selection techniques. However, no such genome assembly has been previously reported. There are 2 species of domestic water buffalo, the river (2n = 50) and the swamp (2n = 48) buffalo. Here we describe a draft quality reference sequence for the river buffalo created from Illumina GA and Roche 454 short read sequences using the MaSuRCA assembler. The assembled sequence is 2.83 Gb, consisting of 366 983 scaffolds with a scaffold N50 of 1.41 Mb and contig N50 of 21 398 bp. Annotation of the genome was supported by transcriptome data from 30 tissues and identified 21 711 predicted protein coding genes. Searches for complete mammalian BUSCO gene groups found 98.6% of curated single copy orthologs present among predicted genes, which suggests a high level of completeness of the genome. The annotated sequence is available from NCBI at accession GCA_000471725.1.

https://github.com/nicolazzie/AffyPipe.git.

The domestic water buffalo ( Bubalus bubalis ) makes a major contribution to the global agricultural economy in the form of milk, meat, hides, and draught power. The global water buffalo population is predominantly found in Asia, and per head of population more people depend upon the buffalo than on any other livestock species. Despite its agricultural importance, there are comparatively fewer genomic and transcriptomic resources available for buffalo than for other livestock species. We have generated a large-scale gene expression atlas covering multiple tissue and cell types from all major organ systems collected from three breeds of riverine water buffalo (Mediterranean, Pandharpuri and Bhadawari) and used the network analysis tool Graphia Professional to identify clusters of genes with similar expression profiles. Alongside similar data, we and others have generated for ruminants as part of the Functional Annotation of Animal Genomes Consortium; this comprehensive transcriptome supports functional annotation and comparative analysis of the water buffalo genome.

Paratuberculosis disease is a chronic bacterial disease infection of ruminants of global relevance, caused by MAP (Mycobacterium avium subsp. paratuberculosis). The present study was conducted on the Garfagnina goat breed that is an Italian native goat population registered on the Tuscan regional repertory of genetic resources at risk of extinction. Forty-eight adult goats (27 serologically positive to MAP-positive and 21 serologically negative to MAP-negative) belonging to a single flock that had experienced annual mortalities due to MAP infection were identified and genotyped with the Illumina GoatSNP60 BeadChip. Diagnosis was achieved by serological tests, as well as post-mortem examination of affected animals. A genome-wide scan was then performed on the individual marker genotypes, in an attempt to identify genomic regions associated with MAP infection disease. Nine significant markers were highlighted and they were located within, or nearby, annotated genes. Two genes found in this study encode are linked to protein kinases that are among the most important enzymes involved in the immune response to Johne's disease, and four genes are involved in the functions of the Golgi complex.

Two methods have been developed for the assessment of conservation priorities on the basis of molecular markers. According to the Weitzman approach, contributions to genetic diversity are derived from genetic distances between populations. Alternatively, diversity within and across populations is optimized by minimizing marker-estimated kinships. We have applied, for the first time, both methods to a comprehensive data set of 69 European cattle breeds, including all cosmopolitan breeds and several local breeds, for which genotypes of 30 microsatellite markers in 25-50 animals per breed have been obtained. Both methods were used to calculate the gain in diversity if a breed was added to a set of nine non-endangered breeds. Weitzman-derived diversities were confounded by genetic drift in isolated populations, which dominates the genetic distances but does not necessarily increase the conservation value of a breed. Marker-estimated kinships across populations were less disturbed by genetic drift than the Weitzman diversities and assigned high conservation values to Mediterranean breeds, which indeed have genetic histories that differ from the non-endangered breeds. Prospects and limitations of marker-assisted decisions on conservation priorities are discussed.