Linkage Disequilibrium-Based Inference of Genome Homology and Chromosomal Rearrangements Between Species.

The aim of this study was to analyze the genomic homology between cattle (Bos taurus) and buffaloes (Bubalus bubalis) and to propose a rearrangement of the buffalo genome through linkage disequilibrium analyses of buffalo SNP markers referenced in the cattle genome assembly and also compare it to the buffalo genome assembly. A panel of bovine SNPs (single nucleotide polymorphisms) was used for hierarchical, non-hierarchical and admixture cluster analyses. Thus, the linkage disequilibrium information between markers of a specific panel of buffalo was used to infer chromosomal rearrangement. Haplotype diversity and imputation accuracy of the submetacentric chromosomes were also analyzed. The genomic homology between the species enabled us to use the bovine genome assembly to recreate a buffalo genomic reference by rearranging the submetacentric chromosomes. The centromere of the submetacentric chromosomes exhibited high linkage disequilibrium and low haplotype diversity. It allowed hypothesizing about chromosome evolution. It indicated that buffalo submetacentric chromosomes are a centric fusion of ancestral acrocentric chromosomes. The chronology of fusions was also suggested. Moreover, a linear regression between buffalo and cattle rearranged assembly and the imputation accuracy indicated that the rearrangement of the chromosomes was adequate. When using the bovine reference genome assembly, the rearrangement of the buffalo submetacentric chromosomes could be done by SNP BTA (chromosome of Bos taurus) calculations: shorter BTA (shorter arm of buffalo chromosome) was given as [(shorter BTA length - SNP position in shorter BTA)] and larger BTA length as [shorter BTA length + (larger BTA length - SNP position in larger BTA)]. Finally, the proposed linkage disequilibrium-based method can be applied to elucidate other chromosomal rearrangement events in other species with the possibility of better understanding the evolutionary relationship between their genomes.

Genome-wide association study applied to type traits related to milk yield in water buffaloes (Bubalus bubalis).

This research aimed to estimate genetic parameters for milk yield and type traits [withers height (WH), croup height (CH), body length (BL), croup length (CL), iliac width (ILW), ischial width (ISW), and thoracic circumference] in Murrah buffaloes and to identify genomic regions related to type traits by applying a single-step genome-wide association study. Data used to estimate the genetic parameters consisted of 601 records of milk yield in the first lactation and the aforementioned type traits. For the single-step genome-wide association study, 322 samples genotyped with a 90K Axiom Buffalo Genotyping array (Thermo Fisher Scientific, Santa Clara, CA) were used. Bivariate analysis revealed that heritability for milk yield (kg) at 305 d was 0.31 ± 0.11, whereas it ranged from 0.22 ± 0.07 to 0.34 ± 0.09 for the studied conformation traits. Based on the percentages of genetic variance explained by windows of 10 markers, there were 16 genomic regions explaining more than 0.5% of the variance for WH, CH, BL, CL, ILW, ISW, and thoracic circumference. Between those regions, 4 were associated with more than 1 trait, suggesting pleiotropic roles for some genes of Bos taurus autosome (BTA) 12 on CL and WH, BTA13 on ISW and ILW, BTA23 on CH and BL, and BTA28 on ISW and BL. Most of these regions coincide with known quantitative trait loci for milk traits. Thus, further studies based on sequence data will help to validate the association of this region with type traits and likely identify the causal mutations.

Genomic analysis of stayability in Nellore cattle.

Stayability, which can be defined as the probability of a cow calving at a certain age when given the opportunity, is an important reproductive trait in beef cattle because it is directly related to herd profitability. The objective of this study was to estimate genetic parameters and to identify possible genomic regions associated with the phenotypic expression of stayability in Nellore cows. The variance components were estimated by Bayesian inference using a threshold animal model that included the systematic effects of contemporary group and sexual precocity and the random effects of animal and residual. The SNP effects were estimated by the single-step genomic BLUP method using information of 2,838 animals (2,020 females and 930 sires) genotyped with the Illumina High-Density BeadChip Array (San Diego, CA, USA). The variance explained by windows formed by 200 consecutive SNPs was used to identify genomic regions of largest effect on the expression of stayability. The heritability was 0.11 ± 0.01 when A matrix (pedigree) was used and 0.14 ± 0.01 when H matrix (relationship matrix that combines pedigree information and SNP data) was used. A total of 147 candidate genes for stayability were identified on chromosomes 1, 2, 5, 6, 9 and 20 and on the X chromosome. New candidate regions for stayability were detected, most of them related to reproductive, immunological and central nervous system functions.

Genome-Wide Association Study for Carcass Traits in an Experimental Nelore Cattle Population.

The purpose of this study was to identify genomic regions associated with carcass traits in an experimental Nelore cattle population. The studied data set contained 2,306 ultrasound records for longissimus muscle area (LMA), 1,832 for backfat thickness (BF), and 1,830 for rump fat thickness (RF). A high-density SNP panel (BovineHD BeadChip assay 700k, Illumina Inc., San Diego, CA) was used for genotyping. After genomic data quality control, 437,197 SNPs from 761 animals were available, of which 721 had phenotypes for LMA, 669 for BF, and 718 for RF. The SNP solutions were estimated using a single-step genomic BLUP approach (ssGWAS), which calculated the variance for windows of 50 consecutive SNPs and the regions that accounted for more than 0.5% of the additive genetic variance were used to search for candidate genes. The results indicated that 12, 18, and 15 different windows were associated to LMA, BF, and RF, respectively. Confirming the polygenic nature of the studied traits, 43, 65, and 53 genes were found in those associated windows, respectively for LMA, BF, and RF. Among the candidate genes, some of them, which already had their functions associated with the expression of energy metabolism, were found associated with fat deposition in this study. In addition, ALKBH3 and HSD17B12 genes, which are related in fibroblast death and metabolism of steroids, were found associated with LMA. The results presented here should help to better understand the genetic and physiologic mechanism regulating the muscle tissue deposition and subcutaneous fat cover expression of Zebu animals. The identification of candidate genes should contribute for Zebu breeding programs in order to consider carcass traits as selection criteria in their genetic evaluation.