Genome-wide association study of reproductive traits in Nellore heifers using Bayesian inference

An important goal of Zebu breeding programs is to improve reproductive performance. A major problem faced with the genetic improvement of reproductive traits is that recording the time for an animal to reach sexual maturity is costly. Another issue is that accurate estimates of breeding values are obtained only a long time after the young bulls have gone through selection. An alternative to overcome these problems is to use traits that are indicators of the reproductive efficiency of the herd and are easier to measure, such as age at first calving. Another problem is that heifers that have conceived once may fail to conceive in the next breeding season, which increases production costs. Thus, increasing heifer's rebreeding rates should improve the economic efficiency of the herd. Response to selection for these traits tends to be slow, since they have a low heritability and phenotypic information is provided only later in the life of the animal. Genome-wide association studies (GWAS) are useful to investigate the genetic mechanisms that underlie these traits by identifying the genes and metabolic pathways involved. Data from 1853 females belonging to the Agricultural Jacarezinho LTDA were used. Genotyping was performed using the BovineHD BeadChip (777 962 single nucleotide polymorphisms (SNPs)) according to the protocol of Illumina - Infinium Assay II ® Multi-Sample HiScan with the unit SQ ™ System. After quality control, 305 348 SNPs were used for GWAS. Forty-two and 19 SNPs had a Bayes factor greater than 150 for heifer rebreeding and age at first calving, respectively. All significant SNPs for age at first calving were significant for heifer rebreeding. These 42 SNPs were next or within 35 genes that were distributed over 18 chromosomes and comprised 27 protein-encoding genes, six pseudogenes and two miscellaneous noncoding RNAs. The use of Bayes factor to determine the significance of SNPs allowed us to identify two sets of 42 and 19 significant SNPs for heifer rebreeding and age at first calving, respectively, which explain 11.35 % and 6.42 % of their phenotypic variance, respectively. These SNPs provide relevant information to help elucidate which genes affect these traits

Association of ADIPOQ, OLR1 and PPARGC1A gene polymorphisms with growth and carcass traits in Nelore cattle

In beef cattle farming, growth and carcass traits are important for genetic breeding programs. Molecular markers can be used to assist selection and increase genetic gain. The ADIPOQ, OLR1 and PPARGC1A genes are involved in lipid synthesis and fat accumulation in adipose tissue. The objective of this study was to identify polymorphisms in these genes and to assess the association with growth and carcass traits in Nelore cattle. A total of 639 animals were genotyped by PCR-RFLP for rs208549452, rs109019599 and rs109163366 in ADIPOQ, OLR1 and PPARGC1A gene, respectively. We analyzed the association of SNPs identified with birth weight, weaning weight, female yearling weight, female hip height, male yearling weight, male hip height, loin eye area, rump fat thickness, and backfat thickness. The OLR1 marker was associated with rump fat thickness and weaning weight (P < 0.05) and the PPARGC1 marker was associated with female yearling weight

Study of whole genome linkage disequilibrium in Nellore cattle

Background: Knowledge of the linkage disequilibrium (LD) between markers is important to establish the number of markers necessary for association studies and genomic selection. The objective of this study was to evaluate the extent of LD in Nellore cattle using a high density SNP panel and 795 genotyped steers.Results: After data editing, 446,986 SNPs were used for the estimation of LD, comprising 2508.4 Mb of the genome. The mean distance between adjacent markers was 4.90 ± 2.89 kb. The minor allele frequency (MAF) was less than 0.20 in a considerable proportion of SNPs. The overall mean LD between marker pairs measured by r2 and |D'| was 0.17 and 0.52, respectively. The LD (r2) decreased with increasing physical distance between markers from 0.34 (1 kb) to 0.11 (100 kb). In contrast to this clear decrease of LD measured by r2, the changes in |D'| indicated a less pronounced decline of LD. Chromosomes BTA1, BTA27, BTA28 and BTA29 showed lower levels of LD at any distance between markers. Except for these four chromosomes, the level of LD (r2) was higher than 0.20 for markers separated by less than 20 kb. At distances < 3 kb, the level of LD was higher than 0.30. The LD (r2) between markers was higher when the MAF threshold was high (0.15), especially when the distance between markers was short.Conclusions: The level of LD estimated for markers separated by less than 30 kb indicates that the High Density Bovine SNP BeadChip will likely be a suitable tool for prediction of genomic breeding values in Nellore cattle. © 2013 Espigolan et al.; licensee BioMed Central Ltd

Genomic association for sexual precocity in beef heifers using pre-selection of genes and haplotype reconstruction

<div><p>Reproductive traits are of the utmost importance for any livestock farming, but are difficult to measure and to interpret since they are influenced by various factors. The objective of this study was to detect associations between known polymorphisms in candidate genes related to sexual precocity in Nellore heifers, which could be used in breeding programs. Records of 1,689 precocious and non-precocious heifers from farms participating in the Conexão Delta G breeding program were analyzed. A subset of single nucleotide polymorphisms (SNP) located in the region of the candidate genes at a distance of up to 5 kb from the boundaries of each gene, were selected from the panel of 777,000 SNPs of the High-Density Bovine SNP BeadChip. Linear mixed models were used for statistical analysis of early heifer pregnancy, relating the trait with isolated SNPs or with haplotype groups. The model included the contemporary group (year and month of birth) as fixed effect and parent of the animal (sire effect) as random effect. The fastPHASE® and GenomeStudio® were used for reconstruction of the haplotypes and for analysis of linkage disequilibrium based on r² statistics. A total of 125 candidate genes and 2,024 SNPs forming haplotypes were analyzed. Statistical analysis after Bonferroni correction showed that nine haplotypes exerted a significant effect (p<0.05) on sexual precocity. Four of these haplotypes were located in the Pregnancy-associated plasma protein-A2 gene (<i>PAPP-A2)</i>, two in the Estrogen-related receptor gamma gene <i>(ESRRG)</i>, and one each in the Pregnancy-associated plasma protein-A gene <i>(PAPP-A)</i>, Kell blood group complex subunit-related family <i>(XKR4)</i> and <i>mannose-binding lectin genes (MBL-1)</i> genes. Although the present results indicate that the <i>PAPP-A2</i>, <i>PAPP-A</i>, <i>XKR4</i>, <i>MBL-1</i> and <i>ESRRG</i> genes influence sexual precocity in Nellore heifers, further studies are needed to evaluate their possible use in breeding programs.</p></div

Frequency, average effect, and respective standard error of haplotype alleles located in the <i>XKR4</i> gene.

<p>Frequency, average effect, and respective standard error of haplotype alleles located in the <i>XKR4</i> gene.</p

Frequency, average effect, and respective standard error of the haplotype alleles located in the <i>MBL-1</i> gene.

<p>Frequency, average effect, and respective standard error of the haplotype alleles located in the <i>MBL-1</i> gene.</p