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.

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.

The aim of this study was to estimate effects of CSN1S1-CSN3 (α(S1)-κ-casein) composite genotypes on milk production traits and milk coagulation properties (MCP) in Mediterranean water buffalo. Genotypes at CSN1S1 and CSN3 and coagulation properties [rennet clotting time (RCT), curd firming time (K₂₀), and curd firmness (A₃₀)] were assessed by reversed-phase HPLC and computerized renneting meter analysis, respectively, using single test-day milk samples of 536 animals. Alternative protein variants of α(S1)-CN and κ-CN were detected by HPLC, and identification of the corresponding genetic variants was carried out by DNA analysis. Two genetic variants were detected at CSN1S1 (A and B variants) and 2 at CSN3 (X1 and X2 variants). Statistical inference was based on a linear model including the CSN1S1-CSN3 composite genotype effect (7 genotypes), the effects of herd-test-day (8 levels), and a combined days in milk (DIM)-parity class. Composite genotype AB-X2X2 was associated with decreased test-day milk yield [-0.21 standard deviation (SD) units of the trait] relative to genotype BB-X2X2. Genotypes did not affect milk protein content, but genotype AB-X1X1 was associated with increased fat content compared with genotype BB-X2X2 (+0.28 SD units of the trait) and AB-X1X1 (+0.43 SD units of the trait). For RCT, the largest difference (+1.91 min; i.e., 0.61 SD units of the trait) was observed between genotype AA-X1X2 and AB-X1X1. Direction of genotype effects on K(20) was consistent with that for RCT. The maximum variation in K₂₀ due to genotype effects (between AA-X1X2 and AB-X1X1 genotypes) was almost 0.9 SD units of the trait. Magnitude of genotype effects was smaller for A₃₀ than for RCT and K₂₀, with a maximum difference of 0.5 SD units of the trait between genotype AA-X1X2 and AA-X1X1. The B allele at CSN1S1 was associated with increased RCT and K₂₀ and with weaker curds compared with allele A. Allele X2 at CSN3 exerted opposite effects on MCP relative to CSN1S1 B. Because of linkage disequilibrium, allele B at CSN1S1 and allele X2 at CSN3 tend to be associated and this likely makes their effects cancel each other. This study indicates a role for casein genes in variation of MCP of buffalo milk. Further studies are necessary to estimate the effects of casein genetic variants on variation of cheese yield.

An association study between the milk yield trait and the stearoyl-CoA desaturase (SCD) polymorphism (g.133A > C) in Italian Mediterranean river buffalo was carried out. A full characterization of the river buffalo SCD promoter region was presented. Genotyping information was provided and a quick method for allelic discrimination was developed. The frequency of the C allele was 0·16. Test-day (TD) records (43 510) of milk production belonging to 226 lactations of 169 buffalo cows were analysed with a mixed linear model in order to estimate the effect of g.133A > C genotype, as well as the effect of parity and calving season. The SCD genotype was significantly associated with milk yield (P = 0·02). The genotype AC showed an over-dominance effect with an average daily milk yield approximately 2 kg/d higher than CC buffaloes. Such a difference represents about 28% more milk/d. The effect of the genotype was constant across lactation stages. The contribution of SCD genotype (r(2)SCD) to the total phenotypic variance in milk yield was equal to 0·12. This report is among the first indications of genetic association between a trait of economic importance in river buffalo. Although such results need to be confirmed with large-scale studies in the same and other buffalo populations, they might offer useful indications for the application of MAS programmes in river buffalo and in the future they might be of great economic interest for the river buffalo dairy industry.

No abstract