Seventy-two pigs of three genetic groups (Brazilian indigenous breed Piau, Commercial line and Crossbred) of both sexes were slaughtered at four live weights (30, 60, 90 and 120 kg). Intramuscular fat (IMF) content in Longissimus dorsi muscle of each animal was extracted and correlated with candidate gene mRNA expression (ATN1, EEF1A2, FABP3, LDLR, MGP, OBSCN, PDHB, TRDN and RYR1). Within slaughter weight of 120 kg, Piau and Crossbred pigs showed higher IMF content (p < 0.05) than commercial animals, with 2.48, 2.08 and 1.00% respectively. Barrows presented higher values of IMF (p < 0.05) than gilts (1.54 and 1.30% respectively). Gene expression of EEF1A2, FABP3, LDLR, OBSCN, PDHB, TRDN and RYR1 were correlated with IMF (p < 0.05) using the whole dataset. For Piau data only, expression of FABP3, LDLR, MGP, OBSCN, PDHB, TRDN and RYR1 showed correlation with IMF (p < 0.05). Genes that have important roles in lipid transportation inside the cell (FABP3) and tissues (LDLR) showed correlation with IMF of, respectively, 0.68 and 0.63 using the whole data set, and 0.90 and 0.91 using data from Piau animals. The highly positive correlation of the LDLR and FAPB3 expression with IMF content may confirm that these genes are important for fat deposition in the porcine L. dorsi muscle.

BackgroundReproductive traits such as number of stillborn piglets (SB) and number of teats (NT) have been evaluated in many genome-wide association studies (GWAS). Most of these GWAS were performed under the assumption that these traits were normally distributed. However, both SB and NT are discrete (e.g. count) variables. Therefore, it is necessary to test for better fit of other appropriate statistical models based on discrete distributions. In addition, although many GWAS have been performed, the biological meaning of the identified candidate genes, as well as their functional relationships still need to be better understood. Here, we performed and tested a Bayesian treatment of a GWAS model assuming a Poisson distribution for SB and NT in a commercial pig line. To explore the biological role of the genes that underlie SB and NT and identify the most likely candidate genes, we used the most significant single nucleotide polymorphisms (SNPs), to collect related genes and generated gene-transcription factor (TF) networks.ResultsComparisons of the Poisson and Gaussian distributions showed that the Poisson model was appropriate for SB, while the Gaussian was appropriate for NT. The fitted GWAS models indicated 18 and 65 significant SNPs with one and nine quantitative trait locus (QTL) regions within which 18 and 57 related genes were identified for SB and NT, respectively. Based on the related TF, we selected the most representative TF for each trait and constructed a gene-TF network of gene-gene interactions and identified new candidate genes.ConclusionsOur comparative analyses showed that the Poisson model presented the best fit for SB. Thus, to increase the accuracy of GWAS, counting models should be considered for this kind of trait. We identified multiple candidate genes (e.g. PTP4A2, NPHP1, and CYP24A1 for SB and YLPM1, SYNDIG1L, TGFB3, and VRTN for NT) and TF (e.g. NF-κB and KLF4 for SB and SOX9 and ELF5 for NT), which were consistent with known newborn survival traits (e.g. congenital heart disease in fetuses and kidney diseases and diabetes in the mother) and mammary gland biology (e.g. mammary gland development and body length).Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-016-0189-x) contains supplementary material, which is available to authorized users.

We performed a genome-wide mapping for the age at first calving (AFC) with the goal of annotating candidate genes that regulate fertility in Nellore cattle. Phenotypic data from 762 cows and 777k SNP genotypes from 2,992 bulls and cows were used. Single nucleotide polymorphism (SNP) effects based on the single-step GBLUP methodology were blocked into adjacent windows of 1 Megabase (Mb) to explain the genetic variance. SNP windows explaining more than 0.40% of the AFC genetic variance were identified on chromosomes 2, 8, 9, 14, 16 and 17. From these windows, we identified 123 coding protein genes that were used to build gene networks. From the association study and derived gene networks, putative candidate genes (e.g., PAPPA, PREP, FER1L6, TPR, NMNAT1, ACAD10, PCMTD1, CRH, OPKR1, NPBWR1 and NCOA2) and transcription factors (TF) (STAT1, STAT3, RELA, E2F1 and EGR1) were strongly associated with female fertility (e.g., negative regulation of luteinizing hormone secretion, folliculogenesis and establishment of uterine receptivity). Evidence suggests that AFC inheritance is complex and controlled by multiple loci across the genome. As several windows explaining higher proportion of the genetic variance were identified on chromosome 14, further studies investigating the interaction across haplotypes to better understand the molecular architecture behind AFC in Nellore cattle should be undertaken.

The genetic improvement of reproductive traits such as the number of teats is essential to the success of the pig industry. As opposite to most SNP association studies that consider continuous phenotypes under Gaussian assumptions, this trait is characterized as a discrete variable, which could potentially follow other distributions, such as the Poisson. Therefore, in order to access the complexity of a counting random regression considering all SNPs simultaneously as covariate under a GWAS modeling, the Bayesian inference tools become necessary. Currently, another point that deserves to be highlighted in GWAS is the genetic dissection of complex phenotypes through candidate genes network derived from significant SNPs. We present a full Bayesian treatment of SNP association analysis for number of teats assuming alternatively Gaussian and Poisson distributions for this trait. Under this framework, significant SNP effects were identified by hypothesis tests using 95% highest posterior density intervals. These SNPs were used to construct associated candidate genes network aiming to explain the genetic mechanism behind this reproductive trait. The Bayesian model comparisons based on deviance posterior distribution indicated the superiority of Gaussian model. In general, our results suggest the presence of 19 significant SNPs, which mapped 13 genes. Besides, we predicted gene interactions through networks that are consistent with the mammals known breast biology (e.g., development of prolactin receptor signaling, and cell proliferation), captured known regulation binding sites, and provided candidate genes for that trait (e.g., TINAGL1 and ICK).

Heat stress is an important issue in the global dairy industry. In tropical areas, an alternative to overcome heat stress is the use of crossbred animals or synthetic breeds, such as the Girolando. In this study, we performed a genome-wide association study (GWAS) and post-GWAS analyses for heat stress in an experimental Gir × Holstein F 2 population. Rectal temperature (RT) was measured in heat-stressed F 2 animals, and the variation between 2 consecutive RT measurements (ΔRT) was used as the dependent variable. Illumina BovineSNP50v1 BeadChip (Illumina Inc., San Diego, CA) and single-SNP approach were used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene-transcription factor (TF) networks, generated from enriched TF. The breed origin of marker alleles in the F 2 population was assigned using the breed of origin of alleles (BOA) approach. Heritability and repeatability estimates (± standard error) for ΔRT were 0.13 ± 0.08 and 0.29 ± 0.06, respectively. Association analysis revealed 6 SNP significantly associated with ΔRT. Genes involved with biological processes in response to heat stress effects (LIF, OSM, TXNRD2, and DGCR8) were identified as putative candidate genes. After performing the BOA approach, the 10% of F 2 animals with the lowest breeding values for ΔRT were classified as low-ΔRT, and the 10% with the highest breeding values for ΔRT were classified as high-ΔRT. On average, 49.4% of low-ΔRT animals had 2 alleles from the Holstein breed (HH), and 39% had both alleles from the Gir breed (GG). In high-ΔRT animals, the average proportion of animals for HH and GG were 1.4 and 50.2%, respectively. This study allowed the identification of candidate genes for ΔRT in Gir × Holstein crossbred animals. According to the BOA approach, Holstein breed alleles could be associated with better response to heat stress effects, which could be explained by the fact that Holstein animals are more affected by heat stress than Gir animals and thus require a genetic architecture to defend the body from the deleterious effects of heat stress. Future studies can provide further knowledge to uncover the genetic architecture underlying heat stress in crossbred cattle.

Genetic diversity is one of the most important issues in studies on conservation of cattle breeds and endangered species. The objective of this study was to estimate the levels of genetic differentiation between locally adapted taurine (Bos taurus taurus) and zebu (Bos taurus indicus) breeds in Brazil, which were genotyped for more than 777,000 SNPs. The fixation index (F ), principal component analysis (PCA), and Bayesian clustering were estimated. The F highlighted genetic differentiation between taurine and zebu breeds. The taurine lines, Caracu and Caracu Caldeano, had significant genetic differentiation (F close to 5%) despite their recent selection for different uses (meat and milk). This genetic variability can be used for conservation of locally adapted animals, as well as for breeding programs on zebu breeds. Introgression of zebu in locally adapted breeds was identified, especially in Curraleiro Pé-Duro breed. The Gyr breed, however, had low breed purity at genomic level due to its very heterogeneous mixing pattern.

Fertility traits are economically important in cattle breeding programs. Scrotal circumference (SC) measures are repeatable, easily obtained, highly heritable, and positively correlated with female fertility traits and sperm quality traits in males. A useful approach to summarize SC measures over time is using nonlinear models, which summarize specific measures of SC in a few parameters with biological interpretation. This approach facilitates the selection of bulls with larger SC and maturity index (K), that is, early maturing animals. Because SC is a sex-limited trait, identifying the underlying genomics of growth curve parameters will allow selection across both males and females. We reported the first multitrait genomewide association study (GWAS) of estimated growth curve parameters for SC data in Brahman cattle. Five widely used nonlinear models were tested to fit a total of 3,612 SC records, measured at 6, 12, 18, and 24 mo of age. The von Bertalanffy model, individually fitted for each animal, best fit this SC data. Parameter estimates SC at maturity (A) and K as well as SC at all ages were jointly analyzed in a GWAS to identify 1-Mb regions most strongly associated with each trait. Heritabilities were 0.25 for K and 0.32 for A and ranged from 0.51 to 0.72 for SC at 6 (SC6), 12 (SC12), 18 (SC18), and 24 mo of age (SC24). An overlapping window on chromosome 14 explaining around 0.8% of genetic variance for K, SC12, SC18, and SC24 was observed. The major positional candidate genes within 1 Mb upstream and downstream of this overlapping window were , , , and . Windows of 1 Mb explaining more than 0.4% of each trait on chromosomes 1, 3, 6, 7, 14, 17, 18, 24, 25, and 26 were identified. Pathways and net-work analyses were indicated through transcription factors playing a role on fertility traits: , , , , , , and . Further validation studies on larger populations or other breeds are required to validate these findings and to improve our understanding of the biology and complex genetic architecture of traits associated with scrotal growth and male fertility in cattle.

Rhipicephalus (Boophilus) microplus is the main cattle ectoparasite in tropical areas. Gir × Holstein crossbred cows are well adapted to different production systems in Brazil. In this context, we performed genome-wide association study (GWAS) and post-GWAS analyses for R. microplus resistance in an experimental Gir × Holstein F population. Single nucleotide polymorphisms (SNP) identified in GWAS were used to build gene networks and to investigate the breed of origin for its alleles. Tick artificial infestations were performed during the dry and rainy seasons. Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) and single-step BLUP procedure was used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene transcription factors networks, generated from enriched transcription factors, identified from the promoter sequences of selected gene sets. The genetic origin of marker alleles in the F population was assigned using the breed of origin of alleles approach. Heritability estimates for tick counts were 0.40 ± 0.11 in the rainy season and 0.54 ± 0.11 in the dry season. The top ten 0.5-Mbp windows with the highest percentage of genetic variance explained by SNP markers were found in chromosomes 10 and 23 for both the dry and rainy seasons. Gene network analyses allowed the identification of genes involved with biological processes relevant to immune system functions (TREM1, TREM2, and CD83). Gene-transcription factors network allowed the identification of genes involved with immune functions (MYO5A, TREML1, and PRSS16). In resistant animals, the average proportion of animals showing significant SNPs with paternal and maternal alleles originated from Gir breed was 44.8% whereas the proportion of animals with both paternal and maternal alleles originated from Holstein breed was 11.3%. Susceptible animals showing both paternal and maternal alleles originated from Holstein breed represented 44.6% on average, whereas both paternal and maternal alleles originated from Gir breed animals represented 9.3%. This study allowed us to identify candidate genes for tick resistance in Gir × Holstein crossbreds in both rainy and dry seasons. According to the origin of alleles analysis, we found that most animals classified as resistant showed 2 alleles from Gir breed, while the susceptible ones showed alleles from Holstein. Based on these results, the identified genes may be thoroughly investigated in additional experiments aiming to validate their effects on tick resistance phenotype in cattle.