Compression distance can discriminate animals by genetic profile, build relationship matrices and estimate breeding values

A conserved haplotype in Wagyu cattle contains RAB4A whose encoded protein regulates glucose trafficking in muscle and fat cells

The Wagyu breed of taurine cattle possess favourable genetics for intramuscular fat (IMF) but genomic loci associated with the trait remain under characterised. Here, we report the identification of a previously unidentified genomic region possessing a particular haplotype structure in Wagyu. Through deployment of a genome-wide haplotype detection analysis that captures regions conserved in a target population but not other populations we screened 100 individual Wagyu and contrasted them with 100 individuals from two independent comparison breeds, Charolais and Angus, using high-density SNPs. An extreme level of Wagyu conservation was assigned to a single genomic window (spanning genomic coordinates BTA28:41 088-300 265 bp). In fact, a five-SNP region spanning 27 096 bp is almost perfectly conserved among the 100 Wagyu individuals assayed and partially overlaps RAB4A. Focussing in, two consecutive SNPs (genomic coordinates 236 949 and 239 950) are apparently fixed within the Wagyu (BB and AA respectively), but at mixed frequencies in the other two breeds. These SNPs are located in the two introns straddling exon 7. In a separate analysis using the 1000 Bulls database, we found that, coincident with exon 7 of RAB4A first allele frequencies were highest in the high IMF Japanese Native (Wagyu) breeds (0.78) and lowest in the low IMF indicine breeds (Nelore and Brahman), with intermediate marbling breeds (Angus and Charolais) assigned intermediate rankings (0.42). RAB4A is known to encode a protein that regulates intracellular trafficking of the insulin-regulated glucose transporter GLUT4. RAB4A can be considered an attractive new positional candidate for IMF development.

RAPID COMMUNICATION: A haplotype information theory method reveals genes of evolutionary interest in European vs. Asian pigs

Asian and European wild boars were independently domesticated ca. 10,000 yr ago. Since the 17th century, Chinese breeds have been imported to Europe to improve the genetics of European animals by introgression of favorable alleles, resulting in a complex mosaic of haplotypes. To interrogate the structure of these haplotypes further, we have run a new haplotype segregation analysis based on information theory, namely compression efficiency (CE). We applied the approach to sequence data from individuals from each phylogeographic region (n = 23 from Asia and Europe) including a number of major pig breeds. Our genome-wide CE is able to discriminate the breeds in a manner reflecting phylogeography. Furthermore, 24,956 nonoverlapping sliding windows (each comprising 1,000 consecutive SNP) were quantified for extent of haplotype sharing within and between Asia and Europe. The genome-wide distribution of extent of haplotype sharing was quite different between groups. Unlike European pigs, Asian pigs haplotype sharing approximates a normal distribution. In line with this, we found the European breeds possessed a number of genomic windows of dramatically higher haplotype sharing than the Asian breeds. Our CE analysis of sliding windows captures some of the genomic regions reported to contain signatures of selection in domestic pigs. Prominent among these regions, we highlight the role of a gene encoding the mitochondrial enzyme LACTB which has been associated with obesity, and the gene encoding MYOG a fundamental transcriptional regulator of myogenesis. The origin of these regions likely reflects either a population bottleneck in European animals, or selective targets on commercial phenotypes reducing allelic diversity in particular genes and/or regulatory regions.

Data compression can discriminate broilers by selection line, detect haplotypes, and estimate genetic potential for complex phenotypes

Accurately establishing the relationships among individuals lays the foundation for genetic analyses such as genome-wide association studies and identification of selection signatures. Of particular interest to the poultry industry are estimates of genetic merit based on molecular data. These estimates can be commercially exploited in marker-assisted breeding programs to accelerate genetic improvement. Here, we test the utility of a new method we have recently developed to estimate animal relatedness and applied it to genetic parameter estimation in commercial broilers. Our approach is based on the concept of data compression from information theory. Using the real-world compressor gzip to estimate normalized compression distance (NCD) we have built compression-based relationship matrices (CRM) for 988 chickens from 4 commercial broiler lines-2 male and 2 female lines. For all pairs of individuals, we found a strong negative relationship between the commonly used genomic relationship matrix (GRM) and NCD. This reflects the fact that "similarity" is the inverse of "distance." The CRM explained more genetic variation than the corresponding GRM in 2 of 3 phenotypes, with corresponding improvements in accuracy of genomic-enabled predictions of breeding value. A sliding-window version of the analysis highlighted haplotype regions of the genome apparently under selection in a line-specific manner. In the male lines, we retrieved high population-specific scores for IGF-1 and a cognate receptor, INSR. For the female lines, we detected an extreme score for a region containing a reproductive hormone receptor (GNRHR). We conclude that our compression-based method is a valid approach to established relationships and identify regions under selective pressure in commercial lines of broiler chickens.

The Bos taurus-Bos indicus balance in fertility and milk related genes

Numerical approaches to high-density single nucleotide polymorphism (SNP) data are often employed independently to address individual questions. We linked independent approaches in a bioinformatics pipeline for further insight. The pipeline driven by heterozygosity and Hardy-Weinberg equilibrium (HWE) analyses was applied to characterize Bos taurus and Bos indicus ancestry. We infer a gene co-heterozygosity network that regulates bovine fertility, from data on 18,363 cattle with genotypes for 729,068 SNP. Hierarchical clustering separated populations according to Bos taurus and Bos indicus ancestry. The weights of the first principal component were subjected to Normal mixture modelling allowing the estimation of a gene’s contribution to the Bos taurus-Bos indicus axis. We used deviation from HWE, contribution to Bos indicus content and association to fertility traits to select 1,284 genes. With this set, we developed a co-heterozygosity network where the group of genes annotated as fertility-related had significantly higher Bos indicus content compared to other functional classes of genes, while the group of genes associated with milk production had significantly higher Bos taurus content. The network analysis resulted in capturing novel gene associations of relevance to bovine domestication events. We report transcription factors that are likely to regulate genes associated with cattle domestication and tropical adaptation. Our pipeline can be generalized to any scenarios where population structure requires scrutiny at the molecular level, particularly in the presence of a priori set of genes known to impact a phenotype of evolutionary interest such as fertility.