A Gene Expression Atlas of the Domestic Water Buffalo (Bubalus bubalis).

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.

Computational detection and experimental validation of segmental duplications and associated copy number variations in water buffalo ( Bubalus bubalis ).

Duplicated sequences are an important source of gene evolution and structural variation within mammalian genomes. Using a read depth approach based on next-generation sequencing, we performed a genome-wide analysis of segmental duplications (SDs) and associated copy number variations (CNVs) in the water buffalo (Bubalus bubalis). By aligning short reads of Olimpia (the reference water buffalo) to the UMD3.1 cattle genome, we identified 1,038 segmental duplications comprising 44.6 Mb (equivalent to ~1.73% of the cattle genome) of the autosomal and X chromosomal sequence in the buffalo genome. We experimentally validated 70.3% (71/101) of these duplications using fluorescent in situ hybridization. We also detected a total of 1,344 CNV regions across 14 additional water buffaloes, amounting to 59.8 Mb of variable sequence or the equivalent of 2.2% of the cattle genome. The CNV regions overlap 1,245 genes that are significantly enriched for specific biological functions including immune response, oxygen transport, sensory system and signal transduction. Additionally, we performed array Comparative Genomic Hybridization (aCGH) experiments using the 14 water buffaloes as test samples and Olimpia as the reference. Using a linear regression model, a high Pearson correlation (r = 0.781) was observed between the log2 ratios between copy number estimates and the log2 ratios of aCGH probes. We further designed Quantitative PCR assays to confirm CNV regions within or near annotated genes and found 74.2% agreement with our CNV predictions. These results confirm sub-chromosome-scale structural rearrangements present in the cattle and water buffalo. The information on genome variation that will be of value for evolutionary and phenotypic studies, and may be useful for selective breeding of both species.

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle.

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.

Genomic structural differences between cattle and River Buffalo identified through comparative genomic and transcriptomic analysis.

Water buffalo (Bubalus bubalis L.) is an important livestock species worldwide. Like many other livestock species, water buffalo lacks high quality and continuous reference genome assembly, required for fine-scale comparative genomics studies. In this work, we present a dataset, which characterizes genomic differences between water buffalo genome and the extensively studied cattle (Bos taurus Taurus) reference genome. This data set is obtained after alignment of 14 river buffalo whole genome sequencing datasets to the cattle reference. This data set consisted of 13,444 deletion CNV regions, and 11,050 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Gene expression data from cattle and buffalo were also presented for genes impacted by these regions. Public assessment of this dataset will allow for further analyses and functional annotation of genes that are potentially associated with phenotypic difference between cattle and water buffalo.

New Insights on Water Buffalo Genomic Diversity and Post-Domestication Migration Routes From Medium Density SNP Chip Data.

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.

Genome assembly and transcriptome resource for river buffalo, Bubalus bubalis (2n = 50).

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 (2 n = 50) and the swamp (2 n = 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.

Design and validation of a 90K SNP genotyping assay for the water buffalo (Bubalus bubalis).

BACKGROUND: The 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. METHODS: The 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. RESULTS: This 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. CONCLUSION: The 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.

Identification of candidate genes for paratuberculosis resistance in the native Italian Garfagnina goat breed.

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.

Evaluation of genetic variability within PrP genotyped sheep of endangered Italian Altamurana breed.

In the last few decades, there has been increased awareness of preservation and exploitation of endangered breeds for the maintenance of biodiversity, as well as the concern for diseases in sheep breeding. This study was carried out in native endangered Altamurana dairy sheep breed from Southern Italy. The Altamurana breed was considered as two populations (Alt-Cav and Alt-Cra-Zoe), based on presumed cross-breed and remarkable differences in the PrP genotypes frequencies. The genetic diversity between the two Altamurana populations (Alt-Cav and Alt-Cra-Zoe) was evaluated in comparison to three Italian dairy breeds through fourteen microsatellite markers. Both measures of genetic distance and the population structure analysis highlighted that the Alt-Cav and Alt-Cra-Zoe sheep have a particular genetic background. The estimated fixation index (FST) and the genetic Nei's distances among the populations showed a higher level of genetic differentiation for Alt-Cav than Alt-Cra-Zoe. The Bayesian clustering analysis implemented in the STRUCTURE software showed clear and distinct clusters for the two Altamurana populations, confirming the hypothesis of Alt-Cav as a genetic group well differentiated from Alt-Cra-Zoe. Alt-Cav likely can be considered as belonging to the original strain of the Altamurana breed. This findings may be used to assist the programme for conservation and selection of scrapie resistance genotypes in endangered Altamurana sheep breed.

PrP genotype frequencies and risk evaluation for scrapie in dairy sheep breeds from southern Italy.

Concerns regarding scrapie in sheep breeding have increased in the last few decades. The present study was carried out in dairy sheep breeds from southern Italy. In order to find breeding animals resistant to scrapie, the PrP genes of 1,205 animals from entire flocks of dairy native Apulian Leccese and Altamurana breeds, and Sicilian Comisana breed, were analysed for polymorphisms at codons 136, 154, and 171 related to scrapie resistance/susceptibility. The Altamurana breed was considered as two populations (Alt-Cav and Alt-Cra-Zoe), based on presumed cross-breeding. A total of five alleles and ten different genotypes were found. The ARQ allele was predominant for all breeds followed by ARR, the most resistant allele to scrapie, which was highly prevalent in Comisana (50%) and in native Alt-Cav (42.4%). The VRQ allele, associated with the highest susceptibility to scrapie, was detected at not negligeable levels in allocthonous Comisana (3.5%), at a low frequency (0.2%) in native Leccese and Alt-Cra-Zoe, while it was absent in Alt-Cav. The frequencies of PrP genotypes with a very low susceptibility risk to scrapie (R1) was higher in Comisana and Alt-Cav. The most susceptible genotype, ARQ/VRQ, was found only in Comisana. Within the Altamurana breed, there were notable differences between Alt-Cav and Alt-Cra-Zoe sheep. The Alt-Cav was characterised by the absence of VRQ and AHQ alleles and by the higher frequency of the ARR/ARR genotype (18.7%). Breeding programs, mainly in endangered breeds such as Altamurana, should be conducted gradually, combining resistance to scrapie, maintenance of genetic variability, and production.

AffyPipe: an open-source pipeline for Affymetrix Axiom genotyping workflow.

UNLABELLED: The Affymetrix Axiom genotyping standard and 'best practice' workflow for Linux and Mac users consists of three stand-alone executable programs (Affymetrix Power Tools) and an R package (SNPolisher). Currently, SNP analysis has to be performed in a step-by-step procedure. Manual intervention and/or programming skills by the user is required at each intermediate point, as Affymetrix Power Tools programs do not produce input files for the program next-in-line. An additional problem is that the output format of genotypes is not compatible with most analysis software currently available. AffyPipe solves all the above problems, by automating both standard and 'best practice' workflows for any species genotyped with the Axiom technology. AffyPipe does not require programming skills and performs all the steps necessary to obtain a final genotype file. Furthermore, users can directly edit SNP probes and export genotypes in PLINK format. AVAILABILITY AND IMPLEMENTATION: https://github.com/nicolazzie/AffyPipe.git.