Genetic factors underlying host resistance to Rhipicephalus microplus tick infestation in Braford cattle: a systems biology perspective.

Approximately 80% of the world's cattle are raised in regions with a high risk of tick-borne diseases, resulting in significant economic losses due to parasitism by Rhipicephalus (Boophilus) microplus. However, the lack of a systemic biology approach hampers a comprehensive understanding of tick-host interactions that mediate tick resistance phenotypes. Here, we conducted a genome-wide association study (GWAS) of 2933 Braford cattle and found 340 single-nucleotide polymorphisms (SNPs) associated with tick counts. Gene expression analyses were performed on skin samples obtained from previously tick-exposed heifers with extremely high or low estimated breeding values for R. microplus counts. Evaluations were performed both before and after artificial infestation with ticks. Differentially expressed genes were found within 1-Mb windows centered at significant SNPs from GWAS. A total of 330 genes were related to the breakdown of homeostasis that was induced by larval attachment to bovine skin. Enrichment analysis pointed to a key role of proteolysis and signal transduction via JAK/STAT, NFKB and WNT/beta catenin signaling pathways. Integrative analysis on matrixEQTL revealed two cis-eQTLs and four significant SNPs in the genes peptidyl arginine deiminase type IV (PADI4) and LOC11449251. The integration of genomic data from QTL maps and transcriptome analyses has identified a set of twelve key genes that show significant associations with tick loads. These genes could be key candidates to improve the accuracy of genomic predictions for tick resistance in Braford cattle.

Tag SNP selection for prediction of adaptation traits in Braford and Hereford cattle using Bayesian methods.

This study utilized Bayesian inference in a genome-wide association study (GWAS) to identify genetic markers associated with traits relevant to the adaptation of Hereford and Braford cattle breeds. We focused on eye pigmentation (EP), weaning hair coat (WHC), yearling hair coat (YHC), and breeding standard (BS). Our dataset comprised 126,290 animals in the pedigree. Out of these, 233 sires were genotyped using high-density (HD) chips, and 3750 animals with medium-density (50 K) single-nucleotide polymorphism (SNP) chips. Employing the Bayes B method with a prior probability of π = 0.99, we identified and tagged single nucleotide polymorphisms (Tag SNPs), ranging from 18 to 117 SNPs depending on the trait. These Tag SNPs facilitated the construction of reduced SNP panels. We then evaluated the predictive accuracy of these panels in comparison to traditional medium-density SNP chips. The accuracy of genomic predictions using these reduced panels varied significantly depending on the clustering method, ranging from 0.13 to 0.65. Additionally, we conducted functional enrichment analysis that found genes associated with the most informative SNP markers in the current study, thereby providing biological insights into the genomic basis of these traits.

Breeding objectives of Brangus cattle in Brazil.

Our goal was to define a breeding objective for Brangus beef cattle in Brazil. Bioeconomic models were produced and used to estimate economic values (EVs). The scenarios simulated were typical full-cycle beef production systems that are used in tropical and subtropical regions. The breeding objective contained pregnancy rate (PR), warm carcass weight (WCW), mature cow weight (MCW), number of nematode eggs per gram of faeces (EPG) and tick count (TICK). Two models were used in series to estimate the EV. A deterministic model was used to simulate effects of PR, WCW and MCW on profitability with a constant parasite load. Subsequently, stochastic models were used to estimate economic values for TICK and EPG as consequences of their environmental effects on weight gains, mortality and health costs. The EV of PR, WCW, MCW, EPG and TICK, was US$1.59, US$2.11, -US$0.24, -US$5.35 and -US$20.88, respectively. Results indicate positive emphasis should be placed on PR (12.49%) and WCW (65.07%) with negative emphasis on MCW (13.92%), EPG (2.77%) and TICK (5.75%). In comparison with the indexes usually used, these results suggest a reformulation in the selection indexes of the beef production system in tropical and subtropical regions in order to obtain greater profitability.

Network analysis uncovers putative genes affecting resistance to tick infestation in Braford cattle skin.

BACKGROUND: Genetic resistance in cattle is considered a suitable way to control tick burden and its consequent losses for livestock production. Exploring tick-resistant (R) and tick-susceptible (S) hosts, we investigated the genetic mechanisms underlying the variation of Braford resistance to tick infestation. Skin biopsies from four-times-artificially infested R (n = 20) and S (n = 19) hosts, obtained before the first and 24 h after the fourth tick infestation were submitted to RNA-Sequencing. Differential gene expression, functional enrichment, and network analysis were performed to identify genetic pathways and transcription factors (TFs) affecting host resistance. RESULTS: Intergroup comparisons of hosts before (Rpre vs. Spre) and after (Rpost vs. Spost) tick infestation found 51 differentially expressed genes (DEGs), of which almost all presented high variation (TopDEGs), and 38 were redundant genes. Gene expression was consistently different between R and S hosts, suggesting the existence of specific anti-tick mechanisms. In the intragroup comparisons, Rpost vs. Rpre and Spost vs. Spre, we found more than two thousand DEGs in response to tick infestation in both resistance groups. Redundant and non-redundant TopDEGs with potential anti-tick functions suggested a role in the development of different levels of resistance within the same breed. Leukocyte chemotaxis was over-represented in both hosts, whereas skin degradation and remodeling were only found in TopDEGs from R hosts. Also, these genes indicated the participation of cytokines, such as IL6 and IL22, and the activation of Wingless (WNT)-signaling pathway. A central gene of this pathway, WNT7A, was consistently modulated when hosts were compared. Moreover, the findings based on a genome-wide association study (GWAS) corroborate the prediction of the WNT-signaling pathway as a candidate mechanism of resistance. The regulation of immune response was the most relevant pathway predicted for S hosts. Members of Ap1 and NF-kB families were the most relevant TFs predicted for R and S, respectively. CONCLUSION: This work provides indications of genetic mechanisms presented by Braford cattle with different levels of resistance in response to tick infestation, contributing to the search of candidate genes for tick resistance in bovine.

Genomic predictions for economically important traits in Brazilian Braford and Hereford beef cattle using true and imputed genotypes.

BACKGROUND: Genomic selection (GS) has played an important role in cattle breeding programs. However, genotyping prices are still a challenge for implementation of GS in beef cattle and there is still a lack of information about the use of low-density Single Nucleotide Polymorphisms (SNP) chip panels for genomic predictions in breeds such as Brazilian Braford and Hereford. Therefore, this study investigated the effect of using imputed genotypes in the accuracy of genomic predictions for twenty economically important traits in Brazilian Braford and Hereford beef cattle. Various scenarios composed by different percentages of animals with imputed genotypes and different sizes of the training population were compared. De-regressed EBVs (estimated breeding values) were used as pseudo-phenotypes in a Genomic Best Linear Unbiased Prediction (GBLUP) model using two different mimicked panels derived from the 50 K (8 K and 15 K SNP panels), which were subsequently imputed to the 50 K panel. In addition, genomic prediction accuracies generated from a 777 K SNP (imputed from the 50 K SNP) were presented as another alternate scenario. RESULTS: The accuracy of genomic breeding values averaged over the twenty traits ranged from 0.38 to 0.40 across the different scenarios. The average losses in expected genomic estimated breeding values (GEBV) accuracy (accuracy obtained from the inverse of the mixed model equations) relative to the true 50 K genotypes ranged from -0.0007 to -0.0012 and from -0.0002 to -0.0005 when using the 50 K imputed from the 8 K or 15 K, respectively. When using the imputed 777 K panel the average losses in expected GEBV accuracy was -0.0021. The average gain in expected EBVs accuracy by including genomic information when compared to simple BLUP was between 0.02 and 0.03 across scenarios and traits. CONCLUSIONS: The percentage of animals with imputed genotypes in the training population did not significantly influence the validation accuracy. However, the size of the training population played a major role in the accuracies of genomic predictions in this population. The losses in the expected accuracies of GEBV due to imputation of genotypes were lower when using the 50 K SNP chip panel imputed from the 15 K compared to the one imputed from the 8 K SNP chip panel.

Tag SNP selection for prediction of tick resistance in Brazilian Braford and Hereford cattle breeds using Bayesian methods.

BACKGROUND: Cattle resistance to ticks is known to be under genetic control with a complex biological mechanism within and among breeds. Our aim was to identify genomic segments and tag single nucleotide polymorphisms (SNPs) associated with tick-resistance in Hereford and Braford cattle. The predictive performance of a very low-density tag SNP panel was estimated and compared with results obtained with a 50 K SNP dataset. RESULTS: BayesB (π = 0.99) was initially applied in a genome-wide association study (GWAS) for this complex trait by using deregressed estimated breeding values for tick counts and 41,045 SNP genotypes from 3455 animals raised in southern Brazil. To estimate the combined effect of a genomic region that is potentially associated with quantitative trait loci (QTL), 2519 non-overlapping 1-Mb windows that varied in SNP number were defined, with the top 48 windows including 914 SNPs and explaining more than 20% of the estimated genetic variance for tick resistance. Subsequently, the most informative SNPs were selected based on Bayesian parameters (model frequency and t-like statistics), linkage disequilibrium and minor allele frequency to propose a very low-density 58-SNP panel. Some of these tag SNPs mapped close to or within genes and pseudogenes that are functionally related to tick resistance. Prediction ability of this SNP panel was investigated by cross-validation using K-means and random clustering and a BayesA model to predict direct genomic values. Accuracies from these cross-validations were 0.27 ± 0.09 and 0.30 ± 0.09 for the K-means and random clustering groups, respectively, compared to respective values of 0.37 ± 0.08 and 0.43 ± 0.08 when using all 41,045 SNPs and BayesB with π = 0.99, or of 0.28 ± 0.07 and 0.40 ± 0.08 with π = 0.999. CONCLUSIONS: Bayesian GWAS model parameters can be used to select tag SNPs for a very low-density panel, which will include SNPs that are potentially linked to functional genes. It can be useful for cost-effective genomic selection tools, when one or a few key complex traits are of interest.

Linkage disequilibrium, persistence of phase and effective population size estimates in Hereford and Braford cattle.

BACKGROUND: The existence of moderate to high levels of linkage disequilibrium (LD) between genetic markers and quantitative trait loci (QTL) affecting traits of interest is fundamental for the success of genome-wide association (GWAS) and genomic selection (GS) studies. Knowledge about the extent and the pattern of LD in livestock populations is essential to determine the density of single nucleotide polymorphisms (SNP) required for accurate GWAS and GS. Moreover, observed LD is related to historical effective population sizes (Ne), and can provide insights into the genetic diversity history of populations. Estimates of the consistency of linkage phase across breeds (R H,B ) can be used to determine if there is sufficient relationship to use pooled reference populations in multi-breed GS programs. The objective of this study was to estimate LD levels, persistence of phase and effective population size in Hereford and Braford cattle populations sampled in Brazil. RESULTS: Mean LD estimates, measured using the squared correlation of alleles at two loci (r (2)), obtained between adjacent SNP across all chromosomes were 0.21 ± 0.27 for Herefords (391 samples with 41,241 SNP) and 0.16 ± 0.22 for Brafords (2044 samples and 41,207 SNP). Estimated r(2) was > 0.2 and 0.3, respectively, for 34 and 25% of adjacent markers in Herefords, and 26 and 17% in Brafords. Estimated N e for Brafords and Herefords at the current generation was 220 and 153 individuals, respectively. The two breeds demonstrated moderate to strong persistence of phase at all distances (R H,B = 0.53 to 0.97). The largest phase correlations were found in the 0 to 50 Kb bins (R H,B = 0.92 to 0.97). Estimated LD decreased rapidly with increasing distance between SNP, however, useful linkage for GWAS and GS (r(2) > 0.2) was found spanning to ~50 Kb. CONCLUSIONS: Panels containing about 50,000 and 150,000 SNP markers are necessary to detect minimal levels of LD between adjacent markers that would be useful for GWAS and GS studies to Hereford and Braford breeds, respectively. Markers are expected to be linked to the same QTL alleles in distances < 50 Kb in both populations due to observed high persistence of phase levels.

Genotype by environment interaction for tick resistance of Hereford and Braford beef cattle using reaction norm models.

BACKGROUND: The cattle tick is a parasite that adversely affects livestock performance in tropical areas. Although countries such as Australia and Brazil have developed genetic evaluations for tick resistance, these evaluations have not considered genotype by environment (G*E) interactions. Genetic gains could be adversely affected, since breedstock comparisons are environmentally dependent on the presence of G*E interactions, particularly if residual variability is also heterogeneous across environments. The objective of this study was to infer upon the existence of G*E interactions for tick resistance of cattle based on various models with different assumptions of genetic and residual variability. METHODS: Data were collected by the Delta G Connection Improvement program and included 10,673 records of tick counts on 4363 animals. Twelve models, including three traditional animal models (AM) and nine different hierarchical Bayesian reaction norm models (HBRNM), were investigated. One-step models that jointly estimate environmental covariates and reaction norms and two-step models based on previously estimated environmental covariates were used to infer upon G*E interactions. Model choice was based on the deviance criterion information. RESULTS: The best-fitting model specified heterogeneous residual variances across 10 subclasses that were bounded by every decile of the contemporary group (CG) estimates of tick count effects. One-step models generally had the highest estimated genetic variances. Heritability estimates were normally higher for HBRNM than for AM. One-step models based on heterogeneous residual variances also usually led to higher heritability estimates. Estimates of repeatability varied along the environmental gradient (ranging from 0.18 to 0.45), which implies that the relative importance of additive and permanent environmental effects for tick resistance is influenced by the environment. Estimated genetic correlations decreased as the tick infestation level increased, with negative correlations between extreme environmental levels, i.e., between more favorable (low infestation) and harsh environments (high infestation). CONCLUSIONS: HBRNM can be used to describe the presence of G*E interactions for tick resistance in Hereford and Braford beef cattle. The preferred model for the genetic evaluation of this population for tick counts in Brazilian climates was a one-step model that considered heteroscedastic residual variance. Reaction norm models are a powerful tool to identify and quantify G*E interactions and represent a promising alternative for genetic evaluation of tick resistance, since they are expected to lead to greater selection efficiency and genetic progress.

Genetic diversity analysis of two commercial breeds of pigs using genomic and pedigree data.

BACKGROUND: Genetic improvement in livestock populations can be achieved without significantly affecting genetic diversity if mating systems and selection decisions take genetic relationships among individuals into consideration. The objective of this study was to examine the genetic diversity of two commercial breeds of pigs. Genotypes from 1168 Landrace (LA) and 1094 Large White (LW) animals from a commercial breeding program in Brazil were obtained using the Illumina PorcineSNP60 Beadchip. Inbreeding estimates based on pedigree (F x) and genomic information using runs of homozygosity (F ROH) and the single nucleotide polymorphisms (SNP) by SNP inbreeding coefficient (F SNP) were obtained. Linkage disequilibrium (LD), correlation of linkage phase (r) and effective population size (N e ) were also estimated. RESULTS: Estimates of inbreeding obtained with pedigree information were lower than those obtained with genomic data in both breeds. We observed that the extent of LD was slightly larger at shorter distances between SNPs in the LW population than in the LA population, which indicates that the LW population was derived from a smaller N e . Estimates of N e based on genomic data were equal to 53 and 40 for the current populations of LA and LW, respectively. The correlation of linkage phase between the two breeds was equal to 0.77 at distances up to 50 kb, which suggests that genome-wide association and selection should be performed within breed. Although selection intensities have been stronger in the LA breed than in the LW breed, levels of genomic and pedigree inbreeding were lower for the LA than for the LW breed. CONCLUSIONS: The use of genomic data to evaluate population diversity in livestock animals can provide new and more precise insights about the effects of intense selection for production traits. Resulting information and knowledge can be used to effectively increase response to selection by appropriately managing the rate of inbreeding, minimizing negative effects of inbreeding depression and therefore maintaining desirable levels of genetic diversity.

Accuracy of genome-wide imputation in Braford and Hereford beef cattle.

BACKGROUND: Strategies for imputing genotypes from the Illumina-Bovine3K, Illumina-BovineLD (6K), BeefLD-GGP (8K), a non-commercial-15K and IndicusLD-GGP (20K) to either Illumina-BovineSNP50 (50K) or to Illumina-BovineHD (777K) SNP panel, as well as for imputing from 50K, GGP-IndicusHD (90iK) and GGP-BeefHD (90tK) to 777K were investigated. Imputation of low density (<50K) genotypes to 777K was carried out in either one or two steps. Imputation of ungenotyped parents (n = 37 sires) with four or more offspring to the 50K panel was also assessed. There were 2,946 Braford, 664 Hereford and 88 Nellore animals, from which 71, 59 and 88 were genotyped with the 777K panel, while all others had 50K genotypes. The reference population was comprised of 2,735 animals and 175 bulls for 50K and 777K, respectively. The low density panels were simulated by masking genotypes in the 50K or 777K panel for animals born in 2011. Analyses were performed using both Beagle and FImpute software. Genotype imputation accuracy was measured by concordance rate and allelic R(2) between true and imputed genotypes. RESULTS: The average concordance rate using FImpute was 0.943 and 0.921 averaged across all simulated low density panels to 50K or to 777K, respectively, in comparison with 0.927 and 0.895 using Beagle. The allelic R(2) was 0.912 and 0.866 for imputation to 50K or to 777K using FImpute, respectively, and 0.890 and 0.826 using Beagle. One and two steps imputation to 777K produced averaged concordance rates of 0.806 and 0.892 and allelic R(2) of 0.674 and 0.819, respectively. Imputation of low density panels to 50K, with the exception of 3K, had overall concordance rates greater than 0.940 and allelic R(2) greater than 0.919. Ungenotyped animals were imputed to 50K panel with an average concordance rate of 0.950 by FImpute. CONCLUSION: FImpute accuracy outperformed Beagle on both imputation to 50K and to 777K. Two-step outperformed one-step imputation for imputing to 777K. Ungenotyped animals that have four or more offspring can have their 50K genotypes accurately inferred using FImpute. All low density panels, except the 3K, can be used to impute to the 50K using FImpute or Beagle with high concordance rate and allelic R(2).

Genome-wide association and genomic prediction for scrotal circumference in Hereford and Braford bulls.

Scrotal circumference (SC) is widely used as a selection criterion for bulls in breeding programs, since it is easily assessed and correlated with several desirable reproductive traits. The objectives of this study were: to perform a genome-wide association study (GWAS) to identify genomic regions associated with SC adjusted for age (SCa) and for both age and weight (SCaw); to select Tag SNPs from GWAS to construct low-density panel for genomic prediction; and to compare the prediction accuracy of the SC through different methods for Braford and Hereford bulls from the same genetic breeding program. Data of SC from 18,172 bulls (30.4 ± 3.7 cm) and of genotypes from 131 sires and 3,545 animals were used. From GWAS, the top 1% of 1-Mb windows were observed on chromosome (BTA) 2, 20, 7, 8, 15, 3, 16, 27, 6 and 8 for SCa and on BTA 8, 15, 16, 21, 19, 2, 6, 5 and 10 for SCaw, representing 17.4% and 18.8% of the additive genetic variance of SCa and SCaw, respectively. The MeSH analysis was able to translate genomic information providing biological meanings of more specific gene functions related to the SCa and SCaw. The genomic enhancement methods, especially single step GBLUP, that combined phenotype and pedigree data with direct genomic values generated gains in accuracy in relation to pedigree BLUP, suggesting that genomic predictions should be applied to improve genetic gain and to narrow the generation interval compared to traditional methods. The proposed Tag-SNP panels may be useful for lower-cost commercial genomic prediction applications in the future, when the number of bulls in the reference population increases for SC in Hereford and Braford breeds.

Selective transcriptional profiling and data analysis strategies for expression quantitative trait loci mapping in outbred F2 populations.

Genetic analysis of transcriptional profiling experiments is emerging as a promising approach for unraveling genes and pathways that underlie variation of complex biological traits. However, these genetical genomics approaches are currently limited by the high cost of microarrays. We studied five different strategies to optimally select subsets of individuals for transcriptional profiling, including (1) maximizing genetic dissimilarity between selected individuals, (2) maximizing the number of recombination events in selected individuals, (3) selecting phenotypic extremes within inferred genotypes of a previously identified quantitative trait locus (QTL), (4) purely random selection, and (5) profiling animals with the highest and lowest phenotypic values within each family-gender subclass. A simulation study was conducted on the basis of a linkage map and marker genotypes were derived from data on chromosome 6 for 510 F2 animals from an existing pig resource population and on a simulated biallelic QTL with pleiotropic effects on performance and gene expression traits. Bivariate analyses were conducted for selected subset sample sizes of 80, 160, and 240 individuals under three different correlation scenarios between the two traits. The genetic dissimilarity and phenotypic extremes within genotype methods had the smallest mean square error on QTL effects and maximum sensitivity on QTL detection, thereby outperforming all other selection strategies, particularly at the smallest proportion of samples selected for gene expression profiling (80/510).

Crossbreeding effects on growth and efficiency in beef cow-calf systems: evaluation of Angus, Caracu, Hereford and Nelore breed direct, maternal and heterosis effects.

The objective of this study was to determine breed additive and heterosis effects on growth curves, total milk yield (TMY), calf weaning weight (WW), predicted energy intake (EI), and cow efficiency (CE) of purebred and crossbred beef cows raised in Southern Brazil. The data were from 175 purebred and crossbred cows representing eight genetic groups: Angus (A), Hereford (H), Nelore (N), A × H (AH), H × A (HA), A × N (AN), N × A (NA), and Caracu (C) × A (CA). Growth of the cows was modeled using the nonlinear Brody function and machine milking was used to assess TMY. WW was linearly adjusted to 210 d. EI was predicted with an equation in which the independent variables were estimates of parameters of the Brody function and TMY. The ratio of WW to EI estimated CE. Taurine-indicine heterosis effects were significant for all traits, and greater than those for taurine breed crosses. In general, crossbred cows were heavier at maturity, matured earlier, produced more milk, weaned heavier calves, and were predicted to consume more energy. Thus, they were more efficient than purebred cows, despite their greater predicted feed intake. Among the purebreds, A cows matured most rapidly, weighed the least at maturity, produced the most milk, weaned the heaviest calves, were predicted to consume the least energy; and were therefore most efficient among the breeds that were evaluated. These results are useful as inputs to bioeconomic models that can be used to predict productive and economic outcomes from crossbreeding and to facilitate recommendations for beef producers of southern Brazil and other similar subtropical climatic regions.

Causes of bruising in carcasses of beef cattle during farm, transport, and slaughterhouse handling in Brazil.

This study aimed to assess preslaughter factors responsible for the occurrence of bruises to front, rib, loin, hip, and round sites of cattle carcasses. Data from the slaughter of 148 loads of cattle were assessed, totaling 4,611 carcasses. Evaluated bruising factor variables were animal gender, journey distance (km), vehicle type (m²), journey time (hours), load density (kg/m2 ), handling and facilities conditions (optimal, acceptable, and unacceptable), and unloading waiting time (min). Gender had the strongest influence on carcass bruising, being the first variable to enter in the majority of models related to the carcass sites. Female carcasses displayed higher bruising in all sites. Densities above 401 kg/m² caused more bruising. With the exception of the loin, unacceptable, and acceptable conditions assigned to farm facilities and handling resulted in higher average bruising. Articulated vehicles with a greater load capacity were more likely to cause bruising. The effects of variables causing bruises act differently in the various carcass sites. In this sense, the results of this study suggested possible management practices for specific problems during the preslaughter period to minimize losses by bruising to different carcass sites.

Direct and maternal breed additive and heterosis effects on growth traits of beef cattle raised in southern Brazil.

The objective of this study was to compare growth from birth to slaughter of different breed groups that were raised in Rio Grande do Sul, Brazil and estimate the consequent breed additive and heterosis effects. Caracu (C), Hereford (H), and Nelore (N) sires were mated with Angus (A) dams, and A sires were mated with H and N dams to produce a first generation of crossbred progeny that was contemporary with purebred A, H, and N calves. Heifers from this first generation (G1) were mated with Brangus (BN) and Braford (BO) sires to produce a second generation (G2) of progeny. Data were analyzed to estimate breed group means, individual and maternal breed additive effects, and heterosis effects on birth weight, weaning weight, preweaning average daily gain, yearling weight, postweaning average daily gain, fattening phase initial weight (around 19 mo), final weight (around 24 mo), average daily gain in the fattening phase, and age at slaughter. In general, crossbred calves outperformed purebred calves. Angus-N and CA crossbred cows weaned heavier calves. Individual taurine-indicine heterosis (Z) significantly increased weaning weight. The AN, NA, and CA steers were heaviest at yearling, whereas NA, CA, AN, and HA had the greatest final weights. However, AH steers were 1 mo older at slaughter than NA contemporaries. Taurine breed effects on postweaning traits and final weight were greater than for N. Maternal breed effects on birth weight and average daily gain in the fattening phase were greater for A and H than for N. In conclusion, heterosis effects were sufficiently large for use of N to be recommended as a component of such systems, despite their relatively low-breed additive effects compared with taurine breeds. Moreover, germplasm from the tropically adapted Bos taurus C may be particularly useful when increased milk production is desired. With the breed and heterosis effects derived in the present study, it is possible to predict the performance and infer which breed and breed crosses will perform better in crossbreeding systems designed for the subtropical conditions of southern Brazil and similar regions.

Genome-wide linkage analysis of global gene expression in loin muscle tissue identifies candidate genes in pigs.

BACKGROUND: Nearly 6,000 QTL have been reported for 588 different traits in pigs, more than in any other livestock species. However, this effort has translated into only a few confirmed causative variants. A powerful strategy for revealing candidate genes involves expression QTL (eQTL) mapping, where the mRNA abundance of a set of transcripts is used as the response variable for a QTL scan. METHODOLOGY/PRINCIPAL FINDINGS: We utilized a whole genome expression microarray and an F(2) pig resource population to conduct a global eQTL analysis in loin muscle tissue, and compared results to previously inferred phenotypic QTL (pQTL) from the same experimental cross. We found 62 unique eQTL (FDR <10%) and identified 3 gene networks enriched with genes subject to genetic control involved in lipid metabolism, DNA replication, and cell cycle regulation. We observed strong evidence of local regulation (40 out of 59 eQTL with known genomic position) and compared these eQTL to pQTL to help identify potential candidate genes. Among the interesting associations, we found aldo-keto reductase 7A2 (AKR7A2) and thioredoxin domain containing 12 (TXNDC12) eQTL that are part of a network associated with lipid metabolism and in turn overlap with pQTL regions for marbling, % intramuscular fat (% fat) and loin muscle area on Sus scrofa (SSC) chromosome 6. Additionally, we report 13 genomic regions with overlapping eQTL and pQTL involving 14 local eQTL. CONCLUSIONS/SIGNIFICANCE: Results of this analysis provide novel candidate genes for important complex pig phenotypes.

Bayesian inference on genetic merit under uncertain paternity.

A hierarchical animal model was developed for inference on genetic merit of livestock with uncertain paternity. Fully conditional posterior distributions for fixed and genetic effects, variance components, sire assignments and their probabilities are derived to facilitate a Bayesian inference strategy using MCMC methods. We compared this model to a model based on the Henderson average numerator relationship (ANRM) in a simulation study with 10 replicated datasets generated for each of two traits. Trait 1 had a medium heritability (h2) for each of direct and maternal genetic effects whereas Trait 2 had a high h2 attributable only to direct effects. The average posterior probabilities inferred on the true sire were between 1 and 10% larger than the corresponding priors (the inverse of the number of candidate sires in a mating pasture) for Trait 1 and between 4 and 13% larger than the corresponding priors for Trait 2. The predicted additive and maternal genetic effects were very similar using both models; however, model choice criteria (Pseudo Bayes Factor and Deviance Information Criterion) decisively favored the proposed hierarchical model over the ANRM model.