Identification of genomic regions associated with feed efficiency in Nelore cattle

Detection of genomic variations and selection signatures in Wagyu using whole-genome sequencing data

Wagyu is recognized for producing marbled beef with high nutritional value and flavor. Reportedly, Wagyu has been widely used to improve the meat quality of local breeds around the world. However, studies on the genetic mechanism of meat quality in Wagyu at the whole-genome level are rarely reported. Here, whole-genome sequencing data of 11 Wagyu and 115 other individuals were used to explore the genomic variations and genes under selection pressure in Wagyu. A total of 31 349 non-synonymous variants and 53 102 synonymous variants were identified in Wagyu. The population structure analysis showed that Wagyu had the closest genetic relationship with Mishima-Ushi cattle and was apparently separated from other cattle breeds. Then, composite likelihood ratio (CLR), integrated haplotype score, fixation index and cross-population composite likelihood ratio (XP-CLR) tests were performed to identify the candidate genes under positive selection in Wagyu. In total, 770 regions containing 312 genes were identified by at least three methods. Among them, 97 regions containing 27 genes were detected by all four methods. We specifically illustrate a list of interesting genes, including LRP2BP, GAA, CACNG6, CXADR, GPCPD1, KLF2, KLF13, SOX5, MYBPC1, SLC25A10, ATP8A1 and MYH15, which are associated with lipid metabolism, fat deposition, muscle development, bone development, feed intake and growth traits in Wagyu. This is the first study to explore the genomic variations and selection signatures of Wagyu at the whole-genome level. These results will provide significant help to beef cattle improvement and breeding.

Intergenerational impact of dietary protein restriction in dairy ewes on epigenetic marks in the perirenal fat of their suckling lambs

In sheep, nutrition during the prepubertal stage is essential for growth performance and mammary gland development. However, the potential effects of nutrient restriction in a prepuberal stage over the progeny still need to be better understood. Here, the intergenerational effect of maternal protein restriction at prepubertal age (2 months of age) on methylation patterns was evaluated in the perirenal fat of Assaf suckling lambs. In total, 17 lambs from ewes subjected to dietary protein restriction (NPR group, 44% less protein) and 17 lambs from control ewes (C group) were analyzed. These lambs were ranked based on their carcass proportion of perirenal and cavitary fat and classified into HighPCF and LowPCF groups. The perirenal tissue from 4 NPR-LowPCF, 4 NPR-HighPCF, 4 C-LowPCF, and 4 C-HighPCF lambs was subjected to whole-genome bisulfite sequencing and differentially methylated regions (DMRs) were identified. Among other relevant processes, these DMRs were mapped in genes responsible for regulating the transition of brown to white adipose tissue and nonshivering thermoregulation, which might be associated with better adaptation/survival of lambs in the perinatal stage. The current study provides important biological insights about the intergenerational effect on the methylation pattern of an NPR in replacement ewes.

Bivariate GWAS reveals pleiotropic regions among feed efficiency and beef quality-related traits in Nelore cattle

Feed-efficient cattle selection is among the most leading solutions to reduce cost for beef cattle production. However, technical difficulties in measuring feed efficiency traits had limited the application in livestock. Here, we performed a Bivariate Genome-Wide Association Study (Bi-GWAS) and presented candidate biological mechanisms underlying the association between feed efficiency and meat quality traits in a half-sibling design with 353 Nelore steers derived from 34 unrelated sires. A total of 13 Quantitative Trait Loci (QTL) were found explaining part of the phenotypic variations. An important transcription factor of adipogenesis in cattle, the TAL1 (rs133408775) gene located on BTA3 was associated with intramuscular fat and average daily gain (IMF-ADG), and a region located on BTA20, close to CD180 and MAST4 genes, both related to fat accumulation. We observed a low positive genetic correlation between IMF-ADG (r = 0.30 ± 0.0686), indicating that it may respond to selection in the same direction. Our findings contributed to clarifying the pleiotropic modulation of the complex traits, indicating new QTLs for bovine genetic improvement.

Genomic and transcriptomic analyses enable the identification of important genes associated with subcutaneous fat deposition in Holstein cows

Subcutaneous fat deposition has many important roles in dairy cattle, including immunological defense and mechanical protection. The main objectives of this study are to identify key candidate genes regulating subcutaneous fat deposition in high-producing dairy cows by integrating genomic and transcriptomic datasets. A total of 1,654 genotyped Holstein cows are used to perform a genome-wide association study (GWAS) aiming to identify genes associated with subcutaneous fat deposition. Subsequently, weighted gene co-expression network analyses (WGCNA) are conducted based on RNA-sequencing data of 34 cows and de-regressed estimated breeding values of subcutaneous fat deposition. Lastly, differentially expressed (DE) mRNA, lncRNA, and differentially alternative splicing genes are obtained for 12 Holstein cows with extreme and divergent phenotypes for subcutaneous fat deposition. Forty-six protein-coding genes are identified as candidate genes regulating subcutaneous fat deposition in Holstein cattle based on the GWAS. Eleven overlapping genes are identified based on the analyses of DE genes and WGCNA. Furthermore, the candidate genes identified based on the GWAS, WGCNA, and analyses of DE genes are significantly enriched for pathways involved in metabolism, oxidative phosphorylation, thermogenesis, fatty acid degradation, and glycolysis/gluconeogenesis pathways. Integrating all findings, the NID2, STARD3, UFC1, DEDD, PPP1R1B, and USP21 genes are considered to be the most important candidate genes influencing subcutaneous fat deposition traits in Holstein cows. This study provides novel insights into the regulation mechanism underlying fat deposition in high-producing dairy cows, which will be useful when designing management and breeding strategies.

Genomic prediction ability and genetic parameters for residual feed intake calculated using different approaches and their associations with growth, reproductive, and carcass traits in Nellore cattle

This study aimed to estimate prediction ability and genetic parameters for residual feed intake (RFI) calculated using a regression equation for each test (RFI_test) and for the whole population (RFI_pop) in Nellore beef cattle. It also aimed to evaluate the correlations between RFI_pop and RFI_test with growth, reproductive, and carcass traits. Genotypic and phenotypic records from 8354 animals were used. An analysis of variance (ANOVA) was performed to verify the adequacy of the regression equations applied to estimate the RFI_test and RFI_pop. The (co)variance components were obtained using the single-step genomic best linear unbiased prediction under single and two-trait animal model analyses. The genetic and phenotypic correlations between RFI_test and RFI_pop with dry matter intake, frame, growth, reproduction, and carcass-related traits were evaluated. The prediction ability and bias were estimated to compare the RFI_test and RFI_pop genomic breeding values (GEBV). The RFI_pop ANOVA showed a higher significance level (p < 0.0001) than did the RFI_test for the fixed effects. The RFI_pop displayed higher additive genetic variance estimated than the RFI_test, although the RFI_pop and RFI_test displayed similar heritabilities. Overall, the RFI_test showed higher residual correlations with growth, reproductive, and carcass traits, while the RFI_pop displayed higher genetic correlations with such traits. The GEBV for the RFI_test was slightly biased than GEBV RFI_pop. The approach to calculate the RFI influenced the decomposition and estimation of variance components and genomic prediction for RFI. The application of RFI_pop would be more appropriate for genetic evaluation purpose to adjust or correct for non-genetic effects and to decrease the prediction bias for RFI.

Transcriptome Profiling of the Liver in Nellore Cattle Phenotypically Divergent for RFI in Two Genetic Groups

The identification and selection of genetically superior animals for residual feed intake (RFI) could enhance productivity and minimize environmental impacts. The aim of this study was to use RNA-seq data to identify the differentially expressed genes (DEGs), known non-coding RNAs (ncRNAs), specific biomarkers and enriched biological processes associated with RFI of the liver in Nellore cattle in two genetic groups. In genetic group 1 (G1), 24 extreme RFI animals (12 low RFI (LRFI) versus 12 high RFI (HRFI)) were selected from a population of 60 Nellore bulls. The RNA-seq of the samples from their liver tissues was performed using an Illumina HiSeq 2000. In genetic group 2 (G2), 20 samples of liver tissue of Nellore bulls divergent for RFI (LRFI, n = 10 versus HRFI, n = 10) were selected from 83 animals. The raw data of the G2 were chosen from the ENA repository. A total of 1811 DEGs were found for the G1 and 2054 for the G2 (p-value ≤ 0.05). We detected 88 common genes in both genetic groups, of which 33 were involved in the immune response and in blocking oxidative stress. In addition, seven (B2M, ADSS, SNX2, TUBA4A, ARHGAP18, MECR, and ABCF3) possible gene biomarkers were identified through a receiver operating characteristic analysis (ROC) considering an AUC > 0.70. The B2M gene was overexpressed in the LRFI group. This gene regulates the lipid metabolism protein turnover and inhibits cell death. We also found non-coding RNAs in both groups. MIR25 was up-regulated and SNORD16 was down-regulated in the LRFI for G1. For G2, up-regulated RNase_MRP and SCARNA10 were found. We highlight MIR25 as being able to act by blocking cytotoxicity and oxidative stress and RMRP as a blocker of mitochondrial damage. The biological pathways associated with RFI of the liver in Nellore cattle in the two genetic groups were for energy metabolism, protein turnover, redox homeostasis and the immune response. The common transcripts, biomarkers and metabolic pathways found in the two genetic groups make this unprecedented work even more relevant, since the results are valid for different herds raised in different ways. The results reinforce the biological importance of these known processes but also reveal new insights into the complexity of the liver tissue transcriptome of Nellore cattle.

Variance component estimation and genome-wide association of predicted methane production in crossbred beef steers

Enteric methane is a potent greenhouse gas and represents an escape of energy from the ruminant digestive system. Additive genetic variation in methane production suggests that genetic selection offers an opportunity to diminish enteric methane emissions. Logistic and monetary difficulties in directly measuring methane emissions can make genetic evaluation on an indicator trait such as predicted methane production a more appealing option, and inclusion of genotyping data can result in greater genetic progress. Three predicted methane production traits were calculated for 830 crossbred steers fed in seven groups. The methane prediction equations used included mathematical models from Ellis et al. (2007), Mills et al. (2003), and IPCC (2019). Pearson correlations between the traits were all greater than 0.99, indicating that each prediction equation behaved similarly. Further, the Spearman correlations between the estimated breeding values for each trait were also 0.99, which suggests any of the predicted methane models could be used without substantially changing the ranking of the selection candidates. The heritabilities of Ellis, Mills, and IPCC predicted methane production were 0.60, 0.62, and 0.59, respectively. A genome-wide association study identified one single nucleotide polymorphism (SNP) that reached the threshold for significance for all of the traits on chromosome 7 related to oxidoreductase activity. Additionally, the SNP slightly below the significance threshold indicate genes related to collagen, intracellular microtubules, and DNA transcription may play a role in predicted methane production or its component traits.

Untargeted metabolomics confirms the association between plasma branched chain amino acids and residual feed intake in beef heifers

This study explored plasma biomarkers and metabolic pathways underlying feed efficiency measured as residual feed intake (RFI) in Charolais heifers. A total of 48 RFI extreme individuals (High-RFI, n = 24; Low-RFI, n = 24) were selected from a population of 142 heifers for classical plasma metabolite and hormone quantification and plasma metabolomic profiling through untargeted LC-MS. Most efficient heifers (Low-RFI) had greater (P = 0.03) plasma concentrations of IGF-1 and tended to have (P = 0.06) a lower back fat depth compared to least efficient heifers. However, no changes were noted (P ≥ 0.10) for plasma concentrations of glucose, insulin, non-esterified fatty acids, β-hydroxybutyrate and urea. The plasma metabolomic dataset comprised 3,457 ions with none significantly differing between RFI classes after false discovery rate correction (FDR > 0.10). Among the 101 ions having a raw P < 0.05 for the RFI effect, 13 were putatively annotated by using internal databases and 6 compounds were further confirmed with standards. Metabolic pathway analysis from these 6 confirmed compounds revealed that the branched chain amino acid metabolism was significantly (FDR < 0.05) impacted by the RFI classes. Our results confirmed for the first time in beef heifers previous findings obtained in male beef cattle and pointing to changes in branched-chain amino acids metabolism along with that of body composition as biological mechanisms related to RFI. Further studies are warranted to ascertain whether there is a cause-and-effect relationship between these mechanisms and RFI.

Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers

Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.

EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle

Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new insights into the impacts of single-nucleotide polymorphisms (SNPs) on transcription factors (TFs) and transcription cofactors (TcoF) coding genes, we carried out a multi-omic analysis to identify cis-regulatory effects of SNPs on these genes' expression in muscle and describe their association with feed efficiency-related traits in Nelore cattle. As a result, we identified one SNP, the rs137256008C > T, predicted to impact the EEF1A1 gene expression (β = 3.02; P-value = 3.51E-03) and the residual feed intake trait (β = - 3.47; P-value = 0.02). This SNP was predicted to modify transcription factor sites and overlaps with several QTL for feed efficiency traits. In addition, co-expression network analyses showed that animals containing the T allele of the rs137256008 SNP may be triggering changes in the gene network. Therefore, our analyses reinforce and contribute to a better understanding of the biological mechanisms underlying gene expression control of feed efficiency traits in bovines. The cis-regulatory SNP can be used as biomarker for feed efficiency in Nelore cattle.

Genome-Wide Selection Signatures and Human-Mediated Introgression Events in Bos taurus indicus-influenced Composite Beef Cattle

Genetic introgression from interbreeding hybridization of European Bos taurus taurus (EBT) and Indian Bos taurus indicus (IBI) cattle breeds have been widely used to combine the climatic resilience of the IBI cattle and the higher productivity of EBT when forming new composite beef cattle (CB) populations. The subsequent breeding strategies have shifted their initial genomic compositions. To uncover population structure, signatures of selection, and potential introgression events in CB populations, high-density genotypes [containing 492,954 single nucleotide polymorphisms (SNPs) after the quality control] of 486 individuals from 15 cattle breeds, including EBT, IBI, and CB populations, along with two Bos grunniens genotypes as outgroup were used in this study. Then, in-depth population genetics analyses were performed for three CB breeds of Beefmaster, Brangus, and Santa Gertrudis. Neighbor-joining, principal components, and admixture analyses confirmed the historical introgression of EBT and IBI haplotypes into CB breeds. The f_dM statistics revealed that only 12.9% of CB populations' genetic components are of IBI origin. The results of signatures of selection analysis indicated different patterns of selection signals in the three CB breeds with primary pressure on pathways involved in protein processing and stress response in Beefmaster, cell proliferation regulation and immune response in Brangus, and amino acids and glucose metabolisms in Santa Gertrudis. An average of >90% of genomic regions underlying selection signatures were of EBT origin in the studied CB populations. Investigating the CB breeds' genome allows the estimation of EBT and IBI ancestral proportions and the locations within the genome where either taurine or indicine origin alleles are under selective pressure. Such findings highlight various opportunities to control the selection process more efficiently and explore complementarity at the genomic level in CB populations.

Meta-analysis across Nellore cattle populations identifies common metabolic mechanisms that regulate feed efficiency-related traits

Background

Feed efficiency (FE) related traits play a key role in the economy and sustainability of beef cattle production systems. The accurate knowledge of the physiologic background for FE-related traits can help the development of more efficient selection strategies for them. Hence, multi-trait weighted GWAS (MTwGWAS) and meta-analyze were used to find genomic regions associated with average daily gain (ADG), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), and residual feed intake (RFI). The FE-related traits and genomic information belong to two breeding programs that perform the FE test at different ages: post-weaning (1,024 animals IZ population) and post-yearling (918 animals for the QLT population).

Results

The meta-analyze MTwGWAS identified 14 genomic regions (-log10(p -value) > 5) regions mapped on BTA 1, 2, 3, 4, 7, 8, 11, 14, 15, 18, 21, and 29. These regions explained a large proportion of the total genetic variance for FE-related traits across-population ranging from 20% (FCR) to 36% (DMI) in the IZ population and from 22% (RFI) to 28% (ADG) in the QLT population. Relevant candidate genes within these regions (LIPE, LPL, IGF1R, IGF1, IGFBP5, IGF2, INS, INSR, LEPR, LEPROT, POMC, NPY, AGRP, TGFB1, GHSR, JAK1, LYN, MOS, PLAG1, CHCD7, LCAT, and PLA2G15) highlighted that the physiological mechanisms related to neuropeptides and the metabolic signals controlling the body's energy balance are responsible for leading to greater feed efficiency. Integrated meta-analysis results and functional pathway enrichment analysis highlighted the major effect of biological functions linked to energy, lipid metabolism, and hormone signaling that mediates the effects of peptide signals in the hypothalamus and whole-body energy homeostasis affecting the genetic control of FE-related traits in Nellore cattle.

Conclusions

Genes and pathways associated with common signals for feed efficiency-related traits provide better knowledge about regions with biological relevance in physiological mechanisms associated with differences in energy metabolism and hypothalamus signaling. These pleiotropic regions would support the selection for feed efficiency-related traits, incorporating and pondering causal variations assigning prior weights in genomic selection approaches.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08671-w.

Host Genome-Metagenome Analyses Using Combinatorial Network Methods Reveal Key Metagenomic and Host Genetic Features for Methane Emission and Feed Efficiency in Cattle

Cattle production is one of the key contributors to global warming due to methane emission, which is a by-product of converting feed stuff into milk and meat for human consumption. Rumen hosts numerous microbial communities that are involved in the digestive process, leading to notable amounts of methane emission. The key factors underlying differences in methane emission between individual animals are due to, among other factors, both specific enrichments of certain microbial communities and host genetic factors that influence the microbial abundances. The detection of such factors involves various biostatistical and bioinformatics methods. In this study, our main objective was to reanalyze a publicly available data set using our proprietary Synomics Insights platform that is based on novel combinatorial network and machine learning methods to detect key metagenomic and host genetic features for methane emission and residual feed intake (RFI) in dairy cattle. The other objective was to compare the results with publicly available standard tools, such as those found in the microbiome bioinformatics platform QIIME2 and classic GWAS analysis. The data set used was publicly available and comprised 1,016 dairy cows with 16S short read sequencing data from two dairy cow breeds: Holstein and Nordic Reds. Host genomic data consisted of both 50 k and 150 k SNP arrays. Although several traits were analyzed by the original authors, here, we considered only methane emission as key phenotype for associating microbial communities and host genetic factors. The Synomics Insights platform is based on combinatorial methods that can identify taxa that are differentially abundant between animals showing high or low methane emission or RFI. Focusing exclusively on enriched taxa, for methane emission, the study identified 26 order-level taxa that combinatorial networks reported as significantly enriched either in high or low emitters. Additionally, a Z-test on proportions found 21/26 (81%) of these taxa were differentially enriched between high and low emitters (p value <.05). In particular, the phylum of Proteobacteria and the order Desulfovibrionales were found enriched in high emitters while the order Veillonellales was found to be more abundant in low emitters as previously reported for cattle (Wallace et al., 2015). In comparison, using the publicly available tool ANCOM only the order Methanosarcinales could be identified as differentially abundant between the two groups. We also investigated a link between host genome and rumen microbiome by applying our Synomics Insights platform and comparing it with an industry standard GWAS method. This resulted in the identification of genetic determinants in cows that are associated with changes in heritable components of the rumen microbiome. Only four key SNPs were found by both our platform and GWAS, whereas the Synomics Insights platform identified 1,290 significant SNPs that were not found by GWAS. Gene Ontology (GO) analysis found transcription factor as the dominant biological function. We estimated heritability of a core 73 taxa from the original set of 150 core order-level taxonomies and showed that some species are medium to highly heritable (0.25–0.62), paving the way for selective breeding of animals with desirable core microbiome characteristics. We identified a set of 113 key SNPs associated with >90% of these core heritable taxonomies. Finally, we have characterized a small set (<10) of SNPs strongly associated with key heritable bacterial orders with known role in methanogenesis, such as Desulfobacterales and Methanobacteriales.

Association between host genetics of sheep and the rumen microbial composition

A synergy between the rumen microbiota and the host genetics has created a symbiotic relationship, beneficial to the host's health. In this study, the association between the host genetics and rumen microbiome of Damara and Meatmaster sheep was investigated. The composition of rumen microbiota was estimated through the analysis of the V3-V4 region of the 16S rRNA gene, while the sheep blood DNA was genotyped with Illumina OvineSNP50 BeadChip and the genome-wide association (GWA) was analyzed. Sixty significant SNPs dispersed in 21 regions across the Ovis aries genome were found to be associated with the relative abundance of seven genera: Acinetobacter, Bacillus, Clostridium, Flavobacterium, Prevotella, Pseudomonas, and Streptobacillus. A total of eighty-four candidate genes were identified, and their functional annotations were mainly associated with immunity responses and function, metabolism, and signal transduction. Our results propose that those candidate genes identified in the study may be modulating the composition of rumen microbiota and further indicating the significance of comprehending the interactions between the host and rumen microbiota to gain better insight into the health of sheep.

Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.

The Quest for Genes Involved in Adaptation to Climate Change in Ruminant Livestock

Livestock radiated out from domestication centres to most regions of the world, gradually adapting to diverse environments, from very hot to sub-zero temperatures and from wet and humid conditions to deserts. The climate is changing; generally global temperature is increasing, although there are also more extreme cold periods, storms, and higher solar radiation. These changes impact livestock welfare and productivity. This review describes advances in the methodology for studying livestock genomes and the impact of the environment on animal production, giving examples of discoveries made. Sequencing livestock genomes has facilitated genome-wide association studies to localize genes controlling many traits, and population genetics has identified genomic regions under selection or introgressed from one breed into another to improve production or facilitate adaptation. Landscape genomics, which combines global positioning and genomics, has identified genomic features that enable animals to adapt to local environments. Combining the advances in genomics and methods for predicting changes in climate is generating an explosion of data which calls for innovations in the way big data sets are treated. Artificial intelligence and machine learning are now being used to study the interactions between the genome and the environment to identify historic effects on the genome and to model future scenarios.

Expression of candidate genes for residual feed intake in tropically adapted Bos taurus and Bos indicus bulls under thermoneutral and heat stress environmental conditions

The objectives of this study were to measure the relative expression of the ATP1A1, NR3C1, POMC, NPY, and LEP genes in Caracu (Bos taurus) and Nelore (Bos indicus) bulls submitted to feed efficiency tests at high environmental temperatures, and to evaluate differences in adaptability to tropical conditions between breeds. Thirty-five Caracu and 30 Nelore bulls were submitted to a feed efficiency test using automated feeding stations. At the end of the test, the animals were subjected to thermoneutral (TN) and heat stress (HS) conditions. Blood samples were collected after the exposure to the TN and HS conditions and the relative expression of genes was measured by qPCR. The bulls exhibited lower expression of ATP1A1 in the HS condition than in the TN condition (1.98 ± 0.27 and 2.86 ± 0.26, P = 0.02), while the relative expression of NR3C1, POMC, and LEP did not differ (P > 0.05) between climatic conditions. The breed and feed intake influenced NPY and LEP expression levels (P < 0.05). Different climate conditions associated with residual feed intake can modify the gene expression patterns of ATP1A1 and NPY. The association observed among all genes studied shows that they are involved in appetite control. Bos taurus and Bos indicus bulls exhibited similar adaptability to tropical climate conditions.

Detecting the footprint of selection on the genomes of Murciano-Granadina goats

Artificial selection is one of the major forces modifying the genetic composition of livestock populations. Identifying genes under selection could be useful to elucidate their impact on phenotypic variation. We aimed to identify genomic regions targeted by selection for dairy and pigmentation traits in Murciano-Granadina goats. Performance of a selection scan based on the integrated haplotype score test in a population of 1183 Murciano-Granadina goats resulted in the identification of 77 candidate genomic regions/SNPs. The most significant selective sweeps mapped to chromosomes 1 (69.86 Mb), 4 (41.80-49.95 Mb), 11 (65.74 Mb), 12 (31.24 and 52.51 Mb), 17 (34.76-37.67 Mb), 22 (31.75 Mb), and 26 (26.69-31.05 Mb). By using previously generated RNA-Seq data, we built a catalogue of 6414 genes that are differentially expressed across goat lactation (i.e. 78 days post-partum, early lactation; 216 days post-partum, late lactation; 285 days post-partum, dry period). Interestingly, 183 of these genes mapped to selective sweeps and several of them display functions related with lipid, protein, and carbohydrate metabolism, insulin signaling, cell proliferation, as well as mammary development and involution. Of particular interest are the CSN3 and CSN1S2 genes, which encode two major milk proteins. Additionally, we found three pigmentation genes (GLI3, MC1R, and MITF) co-localizing with selective sweeps. Performance of a genome-wide association study and Sanger sequencing and TaqMan genotyping experiments revealed that the c.801C>G (p.Cys267Trp) polymorphism in the melanocortin 1 receptor (MC1R) gene is the main determinant of the black (GG or GC genotypes) and brown (CC genotypes) colorations of Murciano-Granadina goats.

Muscle allele-specific expression QTLs may affect meat quality traits in Bos indicus

Single nucleotide polymorphisms (SNPs) located in transcript sequences showing allele-specific expression (ASE SNPs) were previously identified in the Longissimus thoracis muscle of a Nelore (Bos indicus) population consisting of 190 steers. Given that the allele-specific expression pattern may result from cis-regulatory SNPs, called allele-specific expression quantitative trait loci (aseQTLs), in this study, we searched for aseQTLs in a window of 1 Mb upstream and downstream from each ASE SNP. After this initial analysis, aiming to investigate variants with a potential regulatory role, we further screened our aseQTL data for sequence similarity with transcription factor binding sites and microRNA (miRNA) binding sites. These aseQTLs were overlapped with methylation data from reduced representation bisulfite sequencing (RRBS) obtained from 12 animals of the same population. We identified 1134 aseQTLs associated with 126 different ASE SNPs. For 215 aseQTLs, one allele potentially affected the affinity of a muscle-expressed transcription factor to its binding site. 162 aseQTLs were predicted to affect 149 miRNA binding sites, from which 114 miRNAs were expressed in muscle. Also, 16 aseQTLs were methylated in our population. Integration of aseQTL with GWAS data revealed enrichment for traits such as meat tenderness, ribeye area, and intramuscular fat . To our knowledge, this is the first report of aseQTLs identification in bovine muscle. Our findings indicate that various cis-regulatory and epigenetic mechanisms can affect multiple variants to modulate the allelic expression. Some of the potential regulatory variants described here were associated with the expression pattern of genes related to interesting phenotypes for livestock. Thus, these variants might be useful for the comprehension of the genetic control of these phenotypes.

Exploring the genetic architecture of feed efficiency traits in chickens

Chicken feed efficiency (FE) traits are the most important economic traits in broiler production. Several studies evaluating genetic factors affecting food consumption in chickens are available. However, most of these studies identified genomic regions containing putative quantitative trait loci for each trait separately. It is still a challenge to find common gene networks related to these traits. Therefore, here, a genome-wide association study (GWAS) was conducted to explore candidate genomic regions responsible for Feed Intake (FI), Body Weight Gain (BWG) and Feed Conversion Ratio (FCR) traits and their gene networks. A total of 1430 broilers from an experimental population was genotyped with the high density Affymetrix 600K SNP array. A total of 119 associated SNPs located in 20 chromosomes were identified, where some of them were common in more than one FE trait. In addition, novel genomic regions were prospected considering the SNPs dominance effects and sex interaction, identifying putative candidate genes only when these effects were fit in the model. Relevant candidate genes such as ATRNL1, PIK3C2A, PTPRN2, SORCS3 and gga-mir-1759 were highlighted in this study helping to elucidate the genomic architecture of feed efficiency traits. These results provide new insights on the mechanisms underlying the consumption and utilization of food in chickens.

Selection criteria for feed efficiency-related traits and their association with growth, reproductive and carcass traits in Nelore cattle

Context Livestock feed costs have a higher impact on the profitability of beef production systems and are directly related to feed efficiency. However, these traits are hard and have high costs to measure, reducing the availability of phenotypic records and reliability of genetic evaluations. Thus, the use of genomic information can increase the robustness of genetic studies that address them. Aims The aim of the present study was to estimate genetic parameters for feed efficiency, growth, reproductive and carcass traits in Nelore cattle and the correlated response among them, using genomic information. Methods Residual feed intake (RFI), dry-matter intake, feed conversion ratio, feed efficiency (FE), residual average daily gain (RG), residual feed intake and average daily gain (RIG), birthweight, weight at 120, 240, 365 and 450 days of age, scrotal circumference at 365 and 450 days of age, rib-eye area, backfat thickness and rump fat thickness were evaluated. The genetic parameters were estimated using the single-step genomic best linear unbiased prediction approach. Key results The FE-related traits showed low to moderate heritability ranging from 0.07 to 0.23. Feed efficiency-related traits showed low genetic correlations with reproductive (–0.24 to 0.27), carcass (–0.17 to 0.27) and growth (–0.19 to 0.24) traits, except for growth with dry-matter intake (0.32–0.56) and weight at 365 days of age with FE (–0.40). Conclusions The selection to improve growth, reproductive and carcass traits would not change RFI, RG and RIG. The choice of the most adequate selection criterion depends on the production system, that is, RFI might be used for low-input beef cattle systems, and RIG would be used for more intensive and without-any-dietary-restrictions beef cattle systems. Implications The estimates of heritability and genetic correlations suggest that genetic selection for feed efficiency using RFI, RG and RIG in Nellore cattle leads to higher genetic gain than does that using FE and feed conversion ratio without affecting other profitability traits.

Genomic prediction ability for feed efficiency traits using different models and pseudo-phenotypes under several validation strategies in Nelore cattle

There is a growing interest to improve feed efficiency (FE) traits in cattle. The genomic selection was proposed to improve these traits since they are difficult and expensive to measure. Up to date, there are scarce studies about the implementation of genomic selection for FE traits in indicine cattle under different scenarios of pseudo-phenotypes, models, and validation strategies on a commercial large scale. Thus, the aim was to evaluate the feasibility of genomic selection implementation for FE traits in Nelore cattle applying different models and pseudo-phenotypes under validation strategies. Phenotypic and genotypic information from 4 329 and 3 467 animals were used, respectively, which were tested for residual feed intake, DM intake, feed efficiency, feed conversion ratio, residual BW gain, and residual intake and BW gain. Six prediction methods were used: single-step genomic best linear unbiased prediction, Bayes A, Bayes B, Bayes Cπ, Bayesian least absolute shrinkage and selection operator (BLASSO), and Bayes R. Phenotypes adjusted for fixed effects (Y*), estimated breeding value (EBV), and EBV deregressed (DEBV) were used as pseudo-phenotypes. The validation approaches used were: (1) random: the data was randomly divided into ten subsets and the validation was done in each subset at a time; (2) age: the partition into training and testing sets was based on year of birth and testing animals were born after 2016; and (3) EBV accuracy: the data was split into two groups, being animals with accuracy above 0.45 the training set; and below 0.45 the validation set. In the analyses that used the Y* as pseudo-phenotype, prediction ability (PA) was obtained by dividing the correlation between pseudo-phenotype and genomic EBV (GEBV) by the square root of the heritability of the trait. When EBV and DEBV were used as the pseudo-phenotype, the simple correlation of this quantity with the GEBV was considered as PA. The prediction methods show similar results for PA and bias. The random cross-validation presented higher PA (0.17) than EBV accuracy (0.14) and age (0.13). The PA was higher for Y* than for EBV and DEBV (30.0 and 34.3%, respectively). Random validation presented the highest PA, being indicated for use in populations composed mainly of young animals and traits with few generations of data recording. For high heritability traits, the validation can be done by age, enabling the prediction of the next-generation genetic merit. These results would support breeders to identify genomic approaches that are more viable for genomic prediction for FE-related traits.

Gene networks for three feed efficiency criteria reveal shared and specific biological processes

Background

French beef producers suffer from the decrease in profitability of their farms mainly because of the continuous increase in feed costs. Selection for feed efficiency in beef cattle represents a relevant solution to face this problem. However, feed efficiency is a complex trait that can be assessed by three major criteria: residual feed intake (RFI), residual gain (RG) and feed efficiency ratio (FE), which involve different genetic determinisms. An analysis that combines phenotype and whole-genome sequence data provides a unique framework for genomic studies. The aim of our study was to identify the gene networks and the biological processes that are responsible for the genetic determinism that is shared between these three feed efficiency criteria.

Results

A population of 1477 French Charolais young bulls was phenotyped for feed intake (FI), average daily gain (ADG) and final weight (FW) to estimate RFI, RG and FE. A subset of 789 young bulls was genotyped on the BovineSNP50 single nucleotide polymorphism (SNP) array and imputed at the sequence level using RUN6 of the 1000 Bull Genomes Project. We conducted a genome-wide association study (GWAS) to estimate the individual effect of 8.5 million SNPs and applied an association weight matrix (AWM) approach to analyse the results, one for each feed efficiency criterion. The results highlighted co-association networks including 626 genes for RFI, 426 for RG and 564 for FE. Enrichment assessment revealed the biological processes that show the strongest association with RFI, RG and FE, i.e. digestive tract (salivary, gastric and mucin secretion) and metabolic processes (cellular and cardiovascular). Energetic functions were more associated with RFI and FE and cardio-vascular and cellular processes with RG. Several hormones such as apelin, glucagon, insulin, aldosterone, the gonadotrophin releasing hormone and the thyroid hormone were also identified, and these should be tested in future studies as candidate biomarkers for feed efficiency.

Conclusions

The combination of network and pathway analyses at the sequence level led to the identification of both common and specific mechanisms that are involved in RFI, RG and FE, and to a better understanding of the genetic determinism underlying these three criteria. The effects of the genes involved in each of the identified processes need to be tested in genomic evaluations to confirm the potential gain in reliability of using functional variants to select animals for feed efficiency.

Bovine hepatic miRNAome profiling and differential miRNA expression analyses between beef steers with divergent feed efficiency phenotypes

MicroRNAs (miRNAs) are small RNA molecules involved in regulation of multiple biological processes through modulating expression of their target genes. Here we employed RNAseq to profile liver tissue miRNAome of 60 steers from Angus, Charolais, and Kinsella Composite (KC) populations. Of these animals, 36 animals (n = 12 for each breed) were utilized to identify differentially expressed (DE) miRNAs between animals with high (n = 6) or low (n = 6) phenotypic values of residual feed intake (RFI), a common measurement of feed efficiency. At a threshold of fold-change > 1.5 and P-value < 0.05, we detected 12 (7 up- and 5 downregulated in low-RFI animals), 18 (12 up- and 6 downregulated), and 13 (8 up- and 5 downregulated) DE miRNAs for Angus, Charolais, and KC steers, respectively. Most of the DE miRNAs were breed specific, with bta-miR-449a and bta-miR-AB-2 being differentially expressed in all three breeds. The predicted target genes of the identified DE miRNA are mainly involved in cell cycle, cell death and survival, cell signaling, cellular growth and proliferation, protein trafficking, cell morphology, cell-to-cell signaling and interaction, cellular development, molecular transport, post-translational modification, as well as nutrient metabolism (lipids, carbohydrates, protein and amino acid). Our results provide insights into the bovine hepatic miRNAome and their potential roles in molecular regulation of RFI in beef cattle.

Weighted single-step genome-wide association study and pathway analyses for feed efficiency traits in Nellore cattle

The aim was to conduct a weighted single-step genome-wide association study to detect genomic regions and putative candidate genes related to residual feed intake, dry matter intake, feed efficiency (FE), feed conversion ratio, residual body weight gain, residual intake and weight gain in Nellore cattle. Several protein-coding genes were identified within the genomic regions that explain more than 0.5% of the additive genetic variance for these traits. These genes were associated with insulin, leptin, glucose, protein and lipid metabolisms; energy balance; heat and oxidative stress; bile secretion; satiety; feed behaviour; salivation; digestion; and nutrient absorption. Enrichment analysis revealed functional pathways (p-value < .05) such as neuropeptide signalling (GO:0007218), negative regulation of canonical Wingless/Int-1 (Wnt) signalling (GO:0090090), bitter taste receptor activity (GO:0033038), neuropeptide hormone activity (GO:0005184), bile secretion (bta04976), taste transduction (bta0742) and glucagon signalling pathway (bta04922). The identification of these genes, pathways and their respective functions should contribute to a better understanding of the genetic and physiological mechanisms regulating Nellore FE-related traits.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Genetic regulators of mineral amount in Nelore cattle muscle predicted by a new co-expression and regulatory impact factor approach

Mineral contents in bovine muscle can affect meat quality, growth, health, and reproductive traits. To better understand the genetic basis of this phenotype in Nelore (Bos indicus) cattle, we analysed genome-wide mRNA and miRNA expression data from 114 muscle samples. The analysis implemented a new application for two complementary algorithms: the partial correlation and information theory (PCIT) and the regulatory impact factor (RIF), in which we included the estimated genomic breeding values (GEBVs) for the phenotypes additionally to the expression levels, originally proposed for these methods. We used PCIT to determine putative regulatory relationships based on significant associations between gene expression and GEBVs for each mineral amount. Then, RIF was adopted to determine the regulatory impact of genes and miRNAs expression over the GEBVs for the mineral amounts. We also investigated over-represented pathways, as well as pieces of evidences from previous studies carried in the same population and in the literature, to determine regulatory genes for the mineral amounts. For example, NOX1 expression level was positively correlated to Zinc and has been described as Zinc-regulated in humans. Based on our approach, we were able to identify genes, miRNAs and pathways not yet described as underlying mineral amount. The results support the hypothesis that extracellular matrix interactions are the core regulator of mineral amount in muscle cells. Putative regulators described here add information to this hypothesis, expanding the knowledge on molecular relationships between gene expression and minerals.

Identification and Annotation of Potential Function of Regulatory Antisense Long Non-Coding RNAs Related to Feed Efficiency in Bos taurus Bulls

Long non-coding RNAs (lncRNAs) can influence transcriptional and translational processes in mammalian cells and are associated with various developmental, physiological and phenotypic conditions. However, they remain poorly understood and annotated in livestock species. We combined phenotypic, metabolomics and liver transcriptomic data of bulls divergent for residual feed intake (RFI) and fat accretion. Based on a project-specific transcriptome annotation for the bovine reference genome ARS-UCD.1.2 and multiple-tissue total RNA sequencing data, we predicted 3590 loci to be lncRNAs. To identify lncRNAs with potential regulatory influence on phenotype and gene expression, we applied the regulatory impact factor algorithm on a functionally prioritized set of loci (n = 4666). Applying the algorithm of partial correlation and information theory, significant and independent pairwise correlations were calculated and co-expression networks were established, including plasma metabolites correlated with lncRNAs. The network hub lncRNAs were assessed for potential cis-actions and subjected to biological pathway enrichment analyses. Our results reveal a prevalence of antisense lncRNAs positively correlated with adjacent protein-coding genes and suggest their participation in mitochondrial function, acute phase response signalling, TCA-cycle, fatty acid β-oxidation and presumably gluconeogenesis. These antisense lncRNAs indicate a stabilizing function for their cis-correlated genes and a putative regulatory role in gene expression.

Potential Biomarkers for Feed Efficiency-Related Traits in Nelore Cattle Identified by Co-expression Network and Integrative Genomics Analyses

Feed efficiency helps to reduce environmental impacts from livestock production, improving beef cattle profitability. We identified potential biomarkers (hub genes) for feed efficiency, by applying co-expression analysis in Longissimus thoracis RNA-Seq data from 180 Nelore steers. Six co-expression modules were associated with six feed efficiency-related traits (p-value ≤ 0.05). Within these modules, 391 hub genes were enriched for pathways as protein synthesis, muscle growth, and immune response. Trait-associated transcription factors (TFs) ELF1, ELK3, ETS1, FLI1, and TCF4, were identified with binding sites in at least one hub gene. Gene expression of CCDC80, FBLN5, SERPINF1, and OGN was associated with multiple feed efficiency-related traits (FDR ≤ 0.05) and were previously related to glucose homeostasis, oxidative stress, fat mass, and osteoblastogenesis, respectively. Potential regulatory elements were identified, integrating the hub genes with previous studies from our research group, such as the putative cis-regulatory elements (eQTLs) inferred as affecting the PCDH18 and SPARCL1 hub genes related to immune system and adipogenesis, respectively. Therefore, our analyses contribute to a better understanding of the biological mechanisms underlying feed efficiency in bovine and the hub genes disclosed can be used as biomarkers for feed efficiency-related traits in Nelore cattle.

Genetic parameters and genome-wide association study regarding feed efficiency and slaughter traits in Charolais cows

Residual energy intake (REI) on two successive diets (hay and maize based) and slaughter traits, including visceral organs, were phenotyped in 584 adult purebred Charolais cows. To investigate the relationships between these traits and their genetic determinism, we first estimated the genetic parameters, including correlations, using REML modeling under WOMBAT software. The animals were then genotyped on the BovineSNP50 SNPchip before being imputed to the 600K density and genome wide association study was performed with GCTA software. We found low heritability for REI (h2 = 0.12 in each of the diet phases). Although the phenotypic correlation between the two diet phases was moderate (0.36), the genetic correlation was high (0.83), indicating a common genetic determinism for feed efficiency regardless of the diet. Correlations between REI and slaughter traits were negative regarding muscle-related traits and positive for fat-related traits, indicating that efficient animals generally had a more muscular carcass. It was also seen that feed efficiency was genetically and phenotypically correlated with smaller organs when expressed as a proportion of their empty body weight. From the GWAS analysis, seven QTLs were found to be associated with a trait at the genome-wide level of significance and 18 others at the chromosome-wide level. One important QTL was detected in BTA 2, reflecting the essential effect of the myostatin gene on both carcass composition and relative organ weight. Three QTLs were detected for REI during the maize diet phase on BTA 13, 19, and 28, the latter being significant at the genome-wide level. The QTLs on BTA 19 mapped into the TANC2 gene and the QTLs on BTA 28 into the KIF1BP gene, which are both known to interact with the same protein (KIF1A). However, no obvious functional link between these genes and feed efficiency could be made. Among the other QTLs detected, one association on BTA 4 with liver proportion mapped to the candidate gene WASL, which has previously been shown to be differentially expressed in liver cells and linked to feed restriction or cancer development. No QTLs were found to be common between feed efficiency and any slaughter traits.

Genome-Wide Association Studies for Methane Production in Dairy Cattle

Genomic selection has been proposed for the mitigation of methane (CH₄) emissions by cattle because there is considerable variability in CH₄ emissions between individuals fed on the same diet. The genome-wide association study (GWAS) represents an important tool for the detection of candidate genes, haplotypes or single nucleotide polymorphisms (SNP) markers related to characteristics of economic interest. The present study included information for 280 cows in three dairy production systems in Mexico: 1) Dual Purpose (n = 100), 2) Specialized Tropical Dairy (n = 76), 3) Familiar Production System (n = 104). Concentrations of CH₄ in a breath of individual cows at the time of milking (MEIm) were estimated through a system of infrared sensors. After quality control analyses, 21,958 SNPs were included. Associations of markers were made using a linear regression model, corrected with principal component analyses. In total, 46 SNPs were identified as significant for CH₄ production. Several SNPs associated with CH₄ production were found at regions previously described for quantitative trait loci of composition characteristics of meat, milk fatty acids and characteristics related to feed intake. It was concluded that the SNPs identified could be used in genomic selection programs in developing countries and combined with other datasets for global selection.

Gender status effect on carcass and meat quality traits of feedlot Angus × Nellore cattle

The study evaluated the effect of gender status on carcass and meat quality of feedlot Angus × Nellore cattle. A total of 176 cattle, 20 months old, were confined for 190-days and assigned to four treatments: bulls, immunocastrated, steers, and heifers. Bulls had greater rib eye area and HCW (p = 0.0001). Heifers had increased fat thickness (p = 0.0001). Steers and heifers had higher marbling scores (p = 0.0001). There was interaction between gender and aging time for Warner-Bratzler Shear Force (p = 0.0002), L* (p = 0.0118), and b* (p = 0.0113) values of beef. The sensory panel results showed that beef from bulls had the lowest consumer overall acceptance (p = 0.0278). Especially, regardless tenderness, steers and immunocastrated beef were considered tender, independent of aging time. Beef produced by heifers, steers, and immunocastrated is considered to be of higher quality than bulls. Thus, it is may be an interesting alternative to produce high-quality beef than bulls, to attend the consumer demand for high-quality products. Additionally, the low fatty acids n6 levels and low n6:n3 ratio, high levels of CLA, MUFAs, and oleic acid suggests that the heifer meat is favorable for human health.

Differentially expressed mRNAs, proteins and miRNAs associated to energy metabolism in skeletal muscle of beef cattle identified for low and high residual feed intake

Background

Feed efficiency is one of the most important parameters that affect beef production costs. The energy metabolism of skeletal muscle greatly contributes to variations in feed efficiency. However, information regarding differences in proteins involved in the energy metabolism of the skeletal muscle in beef cattle divergently identified for feed efficiency is scarce. In this study, we aimed to investigate energy metabolism of skeletal muscle of Nellore beef cattle, identified for low and high residual feed intake using a proteomics approach. We further assessed the expression of candidate microRNAs as a one of the possible mechanisms controlling the biosynthesis of the proteins involved in energy metabolism that were differentially abundant between high and low residual feed intake animals.

Results

A greater abundance of 14–3-3 protein epsilon (P = 0.01) was observed in skeletal muscle of residual feed intake (RFI) high animals (RFI-High). Conversely, a greater abundance of Heat Shock Protein Beta 1 (P < 0.01) was observed in the skeletal muscle of RFI-Low cattle. A greater mRNA expression of YWHAE, which encodes the 14–3-3 protein epsilon, was also observed in the skeletal muscle of RFI-High animals (P = 0.01). A lower mRNA expression of HSPB1, which encodes the Heat Shock Protein Beta 1, was observed in the skeletal muscle of RFI-High animals (P = 0.01). The miR-665 was identified as a potential regulator of the 14–3-3 protein epsilon, and its expression was greater in RFI-Low animals (P < .001). A greater expression of miR-34a (P = 0.01) and miR-2899 (P < .001) was observed in the skeletal muscle of RFI-High animals, as both miRNAs were identified as potential regulators of HSPB1 expression.

Conclusion

Our results show that Nellore cattle divergently identified for feed efficiency by RFI present changes in the abundance of proteins involved in energy expenditure in skeletal muscle. Moreover, our data point towards that miR-665, miR34a and miR-2899 are likely involved in controlling both 14-3-3 epsilon and HSPB1 proteins identified as differentially abundant in the skeletal muscle of RFI-High and RFI-Low Nellore cattle.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5890-z) contains supplementary material, which is available to authorized users.

Host genetics influence the rumen microbiota and heritable rumen microbial features associate with feed efficiency in cattle

BACKGROUND

The symbiotic rumen microbiota is essential for the digestion of plant fibers and contributes to the variation of production and health traits in ruminants. However, to date, the heritability of rumen microbial features and host genetic components associated with the rumen microbiota, as well as whether such genetic components are animal performance relevant, are largely unknown.

RESULTS

In the present study, we assessed rumen microbiota from a cohort of 709 beef cattle and showed that multiple factors including breed, sex, and diet drove the variation of rumen microbiota among animals. The diversity indices, the relative abundance of ~ 34% of microbial taxa (59 out of 174), and the copy number of total bacteria had a heritability estimate (h²) ≥ 0.15, suggesting that they are heritable elements affected by host additive genetics. These moderately heritable rumen microbial features were also found to be associated with host feed efficiency traits and rumen metabolic measures (volatile fatty acids). Moreover, 19 single nucleotide polymorphisms (SNPs) located on 12 bovine chromosomes were found to be associated with 14 (12 of them had h² ≥ 0.15) rumen microbial taxa, and five of these SNPs were known quantitative trait loci for feed efficiency in cattle.

CONCLUSIONS

These findings suggest that some rumen microbial features are heritable and could be influenced by host genetics, highlighting a potential to manipulate and obtain a desirable and efficient rumen microbiota using genetic selection and breeding. It could be a useful strategy to further improve feed efficiency and optimize rumen fermentation through targeting both cattle and their rumen microbiota.

Genome-Wide Scan Identifies Selection Signatures in Chinese Wagyu Cattle Using a High-Density SNP Array

Selective breeding can lead to genetic diversity and diverse phenotypes in farm animals. Analysis of the genomic regions under selection can provide important insights into the genetic basis of complex traits. In this study, a high-density SNP array was used for analysis of genome selection signatures in Chinese Wagyu cattle. In total, we obtained 478,903 SNPs and 24,820 no-overlap regions for |iHS| (integrated haplotype score) estimations. Under the threshold of the top 1%, 239 regions were finally identified as candidate selected regions and 162 candidate genes were found based on the UMD3.1 genome assembly. These genes were reported to be associated with fatty acids, such as Bos taurus nitric oxide synthase 1 adaptor protein (NOS1AP), Bos taurus hydroxysteroid 17-beta dehydrogenase 7 (HSD17B7), Bos taurus WD repeat domain 7 (WDR7), Bos taurus ELOVL fatty acid elongase 2 (ELOVL2), Bos taurus calpain 1 (CAPN1), Bos taurus parkin RBR E3 ubiquitin protein ligase (PRKN, also known as PARK2), Bos taurus mitogen-activated protein kinase kinase 6 (MAP2K6), meat quality, including Bos taurus ADAM metallopeptidase domain 12 (ADAM12), Bos taurus 5'-aminolevulinate synthase 1 (ALAS1), Bos taurus small integral membrane protein 13 (SMIM13) and Bos taurus potassium two pore domain channel subfamily K member 2 (KCNK2), growth, and developmental traits, such as Bos taurus insulin like growth factor 2 receptor (IGF2R), Bos taurus RAR related orphan receptor A (RORA), Bos taurus fibroblast growth factor 14 (FGF14), Bos taurus paired box 6 (PAX6) and Bos taurus LIM homeobox 6 (LHX6). In addition, we identified several genes that are associated with body size and weight, including Bos taurus sorting nexin 29 (SNX29), Bos taurus zinc finger imprinted 2 (ZIM2), Bos taurus family with sequence similarity 110 member A (FAM110A), immune system, including Bos taurus toll like receptor 9 (TLR9), Bos taurus TAFA chemokine like family member 1 (TAFA1), Bos taurus glutathione peroxidase 8 (putative) (GPX8), Bos taurus interleukin 5 (IL5), Bos taurus PR domain containing 9 (PRDM9), Bos taurus glutamate ionotropic receptor kainate type subunit 2 (GRIK2) and feed intake efficiency, Bos taurus sodium voltage-gated channel alpha subunit 9 (SCN9A), Bos taurus relaxin family peptide/INSL5 receptor 4 (RXFP4), Bos taurus RNA polymerase II associated protein 3 (RPAP3). Moreover, four GO terms of biological regulation (GO:0009987, GO:0008152) and metabolic process (GO:0003824, GO:0005488) were found based on these genes. In addition, we found that 232 candidate regions (~18 Mb) overlapped with the Quantitative trait loci (QTL)regions extracted from cattle QTLdb. Our findings imply that many genes were selected for important traits in Chinese Wagyu cattle. Moreover, these results can contribute to the understanding of the genetic basis of the studied traits during the formation of this population.

Liver proteomics unravel the metabolic pathways related to Feed Efficiency in beef cattle

Improving nutrient utilization efficiency is essential for livestock, given the current scenario of increasing demand for animal protein and sustainable resource use. In this context, understanding the biology of feed efficiency (FE) in beef cattle allows the development of markers for identification and selection of best animals for animal production. Thus, 98 young Nellore bulls were evaluated for FE and at the end of the experiment liver samples from six High Feed Efficient (HFE) and six Low Feed Efficient (LFE) animals were collected for protein extraction, digestion and analysis by HPLC-MS/MS. Data were analyzed for differential abundant proteins (DAPs), protein networks, and functional enrichment. Serum endotoxin was also quantified. We found 42 DAPs and 3 protein networks significantly related to FE. The main pathways associated with FE were: microbial metabolism; biosynthesis of fatty acids, amino acids and vitamins; glycolysis/gluconeogenesis; xenobiotic metabolism and; antigen processing and presentation. Serum endotoxins were significantly higher in LFE animals supporting the results. Therefore, the findings presented here confirmed the altered hepatic metabolism and pronounced hepatic inflammation in LFE animals supporting that the increased bacterial load is at least in part responsible for the hepatic lesions and inflammation in LFE animals.

Systems Biology Reveals NR2F6 and TGFB1 as Key Regulators of Feed Efficiency in Beef Cattle

Systems biology approaches are used as strategy to uncover tissue-specific perturbations and regulatory genes related to complex phenotypes. We applied this approach to study feed efficiency (FE) in beef cattle, an important trait both economically and environmentally. Poly-A selected RNA of five tissues (adrenal gland, hypothalamus, liver, skeletal muscle and pituitary) of eighteen young bulls, selected for high and low FE, were sequenced (Illumina HiSeq 2500, 100 bp, pared-end). From the 17,354 expressed genes considering all tissues, 1,335 were prioritized by five selection categories (differentially expressed, harboring SNPs associated with FE, tissue-specific, secreted in plasma and key regulators) and used for network construction. NR2F6 and TGFB1 were identified and validated by motif discovery as key regulators of hepatic inflammatory response and muscle tissue development, respectively, two biological processes demonstrated to be associated with FE. Moreover, we indicated potential biomarkers of FE, which are related to hormonal control of metabolism and sexual maturity. By using robust methodologies and validation strategies, we confirmed the main biological processes related to FE in Bos indicus and indicated candidate genes as regulators or biomarkers of superior animals.

An integrative transcriptome analysis indicates regulatory mRNA-miRNA networks for residual feed intake in Nelore cattle

Residual Feed Intake (RFI) is an economically relevant trait in beef cattle. Among the molecular regulatory mechanisms, microRNAs (miRNAs) are an important dimension in post-transcriptional regulation and have been associated with different biological pathways. Here, we performed differential miRNAs expression and weighted gene co-expression network analyses (WGCNA) to better understand the complex interactions between miRNAs and mRNAs expressed in bovine skeletal muscle and liver. MiRNA and mRNA expression data were obtained from Nelore steers that were genetically divergent for RFI (N = 10 [low RFI or feed efficient]; N = 10 [high RFI or feed inefficient]). Differentially expressed and hub miRNAs such as bta-miR-486, bta-miR-7, bta-miR15a, bta-miR-21, bta-miR 29, bta- miR-30b, bta-miR-106b, bta-miR-199a-3p, bta-miR-204, and bta-miR 296 may have a potential role in variation of RFI. Functional enrichment analysis of differentially expressed (DE) miRNA's target genes and miRNA-mRNA correlated modules revealed that insulin, lipid, immune system, oxidative stress and muscle development signaling pathways might potentially be involved in RFI in this population. Our study identified DE miRNAs, miRNA - mRNA regulatory networks and hub miRNAs related to RFI. These findings suggest a possible role of miRNAs in regulation of RFI, providing new insights into the potential molecular mechanisms that control feed efficiency in Nelore cattle.

Development and validation of a small SNP panel for feed efficiency in beef cattle

The objective of this study was to develop and validate a customized cost-effective single nucleotide polymorphism (SNP) panel for genetic improvement of feed efficiency in beef cattle. The SNPs identified in previous association studies and through extensive analysis of candidate genomic regions and genes, were screened for their functional impact and allele frequency in Angus and Hereford breeds used as validation candidates for the panel. Association analyses were performed on genotypes of 159 SNPs from new samples of Angus (n = 160), Hereford (n = 329), and Angus-Hereford crossbred (n = 382) cattle using allele substitution and genotypic models in ASReml. Genomic heritabilities were estimated for feed efficiency traits using the full set of SNPs, SNPs associated with at least one of the traits (at P ≤ 0.05 and P < 0.10), as well as the Illumina bovine 50K representing a widely used commercial genotyping panel. A total of 63 SNPs within 43 genes showed association (P ≤ 0.05) with at least one trait. The minor alleles of SNPs located in the GHR and CAST genes were associated with decreasing effects on residual feed intake (RFI) and/or RFI adjusted for backfat (RFIf), whereas minor alleles of SNPs within MKI67 gene were associated with increasing effects on RFI and RFIf. Additionally, the minor allele of rs137400016 SNP within CNTFR was associated with increasing average daily gain (ADG). The SNPs genotypes within UMPS, SMARCAL, CCSER1, and LMCD1 genes showed significant over-dominance effects whereas other SNPs located in SMARCAL1, ANXA2, CACNA1G, and PHYHIPL genes showed additive effects on RFI and RFIf. Gene enrichment analysis indicated that gland development, as well as ion and cation transport are important physiological mechanisms contributing to variation in feed efficiency traits. The study revealed the effect of the Jak-STAT signaling pathway on feed efficiency through the CNTFR, OSMR, and GHR genes. Genomic heritability using the 63 significant (P ≤ 0.05) SNPs was 0.09, 0.09, 0.13, 0.05, 0.05, and 0.07 for ADG, dry matter intake, midpoint metabolic weight, RFI, RFIf, and backfat, respectively. These SNPs contributed to genetic variation in the studied traits and thus can potentially be used or tested to generate cost-effective molecular breeding values for feed efficiency in beef cattle.

Transcriptome analyses reveal reduced hepatic lipid synthesis and accumulation in more feed efficient beef cattle

The genetic mechanisms controlling residual feed intake (RFI) in beef cattle are still largely unknown. Here we performed whole transcriptome analyses to identify differentially expressed (DE) genes and their functional roles in liver tissues between six extreme high and six extreme low RFI steers from three beef breed populations including Angus, Charolais, and Kinsella Composite (KC). On average, the next generation sequencing yielded 34 million single-end reads per sample, of which 87% were uniquely mapped to the bovine reference genome. At false discovery rate (FDR) < 0.05 and fold change (FC) > 2, 72, 41, and 175 DE genes were identified in Angus, Charolais, and KC, respectively. Most of the DE genes were breed-specific, while five genes including TP53INP1, LURAP1L, SCD, LPIN1, and ENSBTAG00000047029 were common across the three breeds, with TP53INP1, LURAP1L, SCD, and LPIN1 being downregulated in low RFI steers of all three breeds. The DE genes are mainly involved in lipid, amino acid and carbohydrate metabolism, energy production, molecular transport, small molecule biochemistry, cellular development, and cell death and survival. Furthermore, our differential gene expression results suggest reduced hepatic lipid synthesis and accumulation processes in more feed efficient beef cattle of all three studied breeds.

Genome-wide scan reveals population stratification and footprints of recent selection in Nelore cattle

Background

This study aimed at (1) assessing the genomic stratification of experimental lines of Nelore cattle that have experienced different selection regimes for growth traits, and (2) identifying genomic regions that have undergone recent selection. We used a sample of 763 animals genotyped with the Illumina BovineHD BeadChip, among which 674 animals originated from two lines that are maintained under directional selection for increased yearling body weight and 89 animals from a control line that is maintained under stabilizing selection.

Results

Multidimensional analysis of the genomic dissimilarity matrix and admixture analysis revealed a substantial level of population stratification between the directional selection lines and the stabilizing selection control line. Two of the three tests used to detect selection signatures (F_ST, XP-EHH and iHS) revealed six candidate regions with indications of selection, which strongly indicates truly positive signals. The set of identified candidate genes included several genes with roles that are functionally related to growth metabolism, such as COL14A1, CPT1C, CRH, TBC1D1, and XKR4.

Conclusions

The current study identified genetic stratification that resulted from almost four decades of divergent selection in an experimental Nelore population, and highlighted autosomal genomic regions that present patterns of recent selection. Our findings provide a basis for a better understanding of the metabolic mechanism that underlies the growth traits, which are modified by selection for yearling body weight.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0381-2) contains supplementary material, which is available to authorized users.

Genome-wide association study of Stayability and Heifer Pregnancy in Red Angus cattle

Reproductive performance is the most important component of cattle production from the standpoint of economic sustainability of commercial beef enterprises. Heifer Pregnancy (HPG) and Stayability (STAY) genetic predictions are 2 selection tools published by the Red Angus Association of America (RAAA) to assist with improvements in reproductive performance. Given the importance of HPG and STAY to the profitability of commercial beef enterprises, the objective of this study was to identify QTL associated with both HPG and STAY in Red Angus cattle. A genome-wide association study (GWAS) was performed using deregressed HPG and STAY EBV, calculated using a single-trait animal model and a 3-generation pedigree with data from the Spring 2015 RAAA National Cattle Evaluation. Each individual animal possessed 74,659 SNP genotypes. Individual animals with a deregressed EBV reliability > 0.05 were merged with the genotype file and marker quality control was performed. Criteria for sifting genotypes consisted of removing those markers where any of the following were found: average call rate less than 0.85, minor allele frequency < 0.01, lack of Hardy-Weinberg equilibrium (P < 0.0001), or extreme linkage disequilibrium (r2 > 0.99). These criteria resulted in 2,664 animals with 62,807 SNP available for GWAS. Association studies were performed using a Bayes Cπ model in the BOLT software package. Marker significance was calculated as the posterior probability of inclusion (PPI), or the number of instances a specific marker was sampled divided by the total number of samples retained from the Markov chain Monte Carlo chains. Nine markers, with a PPI ≥ 3% were identified as QTL associated with HPG on BTA 1, 11, 13, 23, and 29. Twelve markers, with a PPI ≥ 75% were identified as QTL associated with STAY on BTA 6, 8, 9, 12, 15, 18, 22, and 23.

Genome-wide association analyses based on a multiple-trait approach for modeling feed efficiency

Genome-wide association (GWA) of feed efficiency (FE) could help target important genomic regions influencing FE. Data provided by an international dairy FE research consortium consisted of phenotypic records on dry matter intakes (DMI), milk energy (MILKE), and metabolic body weight (MBW) on 6,937 cows from 16 stations in 4 counties. Of these cows, 4,916 had genotypes on 57,347 single nucleotide polymorphism (SNP) markers. We compared a GWA analysis based on the more classical residual feed intake (RFI) model with one based on a previously proposed multiple trait (MT) approach for modeling FE using an alternative measure (DMI|MILKE,MBW). Both models were based on a single-step genomic BLUP procedure that allowed the use of phenotypes from both genotyped and nongenotyped cows. Estimated effects for single SNP markers were small and not statistically important but virtually identical for either FE measure (RFI vs. DMI|MILKE,MBW). However, upon further refining this analysis to develop joint tests within nonoverlapping 1-Mb windows, significant associations were detected between either measure of FE with a window on each of Bos taurus autosomes BTA12 and BTA26. There was, as expected, no overlap between detected genomic regions for DMI|MILKE,MBW and genomic regions influencing the energy sink traits (i.e., MILKE and MBW) because of orthogonal relationships clearly defined between the various traits. Conversely, GWA inferences on DMI can be demonstrated to be partly driven by genetic associations between DMI with these same energy sink traits, thereby having clear implications when comparing GWA studies on DMI to GWA studies on FE-like measures such as RFI.

Widespread modulation of gene expression by copy number variation in skeletal muscle

Copy number variation (CNV) is a frequently observed deviation from the diploid state due to duplication or deletion of genomic regions. Although intensively analyzed for association with diseases and production traits, the specific mechanisms and extent by which such variations affect the phenotype are incompletely understood. We present an integrative study on CNV and genome-wide gene expression in Brazilian Bos indicus cattle. We analyzed CNVs inferred from SNP-chip data for effects on gene expression measured with RNA-seq in skeletal muscle samples of 183 steers. Local effects, where expression changes coincided with CNVs in the respective genes, were restricted to immune genes. Distal effects were attributable to several high-impact CNVs that modulated remote expression in an orchestrated and intertwined fashion. These CNVs were located in the vicinity of major skeletal muscle pathway regulators and associated genes were enriched for proteolysis, autophagy, and muscle structure development. From association analysis between CNVs and several meat quality and production traits, we found CNV-associated expression effects to also manifest at the phenotype level. Based on genome sequences of the population founders, we further demonstrate that CNVs with impact on expression and phenotype are passed on from one generation to another.

Association of residual feed intake with growth and slaughtering performance, blood metabolism, and body composition in growing lambs

The aim of this study was to determine the association of residual feed intake (RFI) with growth performance, blood metabolic parameters, and body composition factors in growing lambs. Individual body weight (BW) and dry matter intake (DMI) were determined in 137 male Hu lambs that were given a pellet feed four times a day for 50 d. RFI did not show a correlation with metabolic BW (MBW) or average daily gain (ADG), but it showed a positive correlation with DMI and feed conversation ratio (FCR). Organ weight and intestine length had a large influence on RFI in lambs. The low-RFI lambs have smaller rumen and longer duodenum indicating the less feed intake and more sufficient absorption rate of low-RFI lambs. The smaller organs like liver, lung and kidney in low-RFI lambs may be related to lower energy consumption and slower metabolic rate. The observed bigger testis was in low-RFI lambs was another cause of the improved feed efficiency. Finally, the plasma concentrations of thyroxine (T4) and adrenocorticotropic hormone (ACTH) were lower in the ELow-RFI group than in the EHigh-RFI group. This study provides new insight into the biological processes underlying variations in feed efficiency in growing lambs.

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.

Genome-wide association analysis and functional annotation of positional candidate genes for feed conversion efficiency and growth rate in pigs

Feed conversion efficiency is a measure of how well an animal converts feed into live weight and it is typically expressed as feed conversion ratio (FCR). FCR and related traits like growth rate (e.g. days to 110 kg—D110) are of high interest for animal breeders, farmers and society due to implications on animal performance, feeding costs and environmental sustainability. The objective of this study was to identify genomic regions associated with FCR and D110 in pigs. A total of 952 terminal line boars, showing an individual variation in FCR, were genotyped using 60K SNP-Chips. Markers were tested for associations with estimated breeding values (EBV) for FCR and D110. For FCR, the largest number of associated SNPs was located on chromosomes 4 (30 SNPs), 1 (25 SNPs), X (15 SNPs) and 6 (12 SNPs). The most prominent genomic regions for D110 were identified on chromosomes 15 (10 SNPs), 1 and 4 (both 9 SNPs). The most significantly associated SNPs for FCR and D110 mapped 129.8 Kb from METTL11B (chromosome 4) and 32Kb from MBD5 (chromosome 15), respectively. A list of positional genes, closest to significantly associated SNPs, was used to identify enriched pathways and biological functions related to the QTL for both traits. A number of candidate genes were significantly overrepresented in pathways of immune cell trafficking, lymphoid tissue structure, organ morphology, endocrine system function, lipid metabolism, and energy production. After resequencing the coding region of selected positional and functional candidate genes, six SNPs were genotyped in a subset of boars. SNPs in PRKDC, SELL, NR2E1 and AKRIC3 showed significant associations with EBVs for FCR/D110. The study revealed a number of chromosomal regions and candidate genes affecting FCR/D110 and pointed to corresponding biological pathways related to lipid metabolism, olfactory reception, and also immunological status.

Genomewide association of piglet responses to infection with one of two porcine reproductive and respiratory syndrome virus isolates

Porcine reproductive and respiratory syndrome (PRRS) is a devastating disease in the swine industry. Identification of host genetic factors that enable selection for improved performance during PRRS virus (PRRSV) infection would reduce the impact of this disease on animal welfare and production efficiency. We conducted genomewide association study (GWAS) analyses of data from 13 trials of approximately 200 commercial crossbred nursery-age piglets that were experimentally infected with 1 of 2 type 2 isolates of PRRSV (NVSL 97-7985 [NVSL] and KS2006-72109 [KS06]). Phenotypes analyzed were viral load (VL) in blood during the first 21 d after infection (dpi) and weight gain (WG) from 0 to 42 dpi. We accounted for the previously identified QTL in the region on SSC4 in our models to increase power to identify additional regions. Many regions identified by single-SNP analyses were not identified using Bayes-B, but both analyses identified the same regions on SSC3 and SSC5 to be associated with VL in the KS06 trials and on SSC6 in the NVSL trials ( < 5 × 10); for WG, regions on SSC5 and SSC17 were associated in the NVSL trials ( < 3 × 10). No regions were identified with either method for WG in the KS06 trials. Except for the region on SSC4, which was associated with VL for both isolates (but only with WG for NVSL), identified regions did not overlap between the 2 PRRSV isolate data sets, despite high estimates of the genetic correlation between isolates for traits based on these data. We also identified genomic regions whose associations with VL or WG interacted with either PRRSV isolate or with genotype at the SSC4 QTL. Gene ontology (GO) annotation terms for genes located near moderately associated SNP ( < 0.003) were enriched for multiple immunologically (VL) and metabolism- (WG) related GO terms. The biological relevance of these regions suggests that, although it may increase the number of false positives, the use of single-SNP analyses and a relaxed threshold also increased the identification of true positives. In conclusion, although only the SSC4 QTL was associated with response to both PRRSV isolates, genes near associated SNP were enriched for the same GO terms across PRRSV isolates, suggesting that host responses to these 2 isolates are affected by the actions of many genes that function together in similar biological processes.