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.

Genetic evaluation for latent variables derived from factor analysis in broilers.

1. The aim of this study was to investigate the associations between several carcass, performance and meat quality traits in broilers through factor analysis and use the latent variables (i.e. factors) as pseudo-phenotypes in genetic evaluations.2. Factors were extracted using the principal components method and varimax rotation algorithm. Genetic parameters were estimated via Bayesian inference under a multiple-trait animal model.3. All factors taken together explained 71% of the original variance of the data. The first factor, denominated as 'weight', was associated with carcass and body weight traits; and the second factor, defined as 'tenderness', represented traits related to water-holding capacity and shear force. The third factor, 'colour', was associated with traits related to meat colour, whereas the fourth, referenced as 'viscera', was related to heart, liver and abdominal fat.4. The four biological factors presented moderate to high heritability (ranging from 0.35 to 0.75), which may confer genetic gains in this population.5. In conclusion, it seems possible to reduce the number of traits in the genetic evaluation of broilers using latent variables derived from factor analysis.

Genotypic and allelic frequencies of gene polymorphisms associated with meat tenderness in Nellore beef cattle.

The objectives of this study were to characterize the allelic and genotypic frequencies of polymorphisms in the µ-calpain and calpastatin genes, and to assess their association with meat tenderness and animal growth in Nellore cattle. We evaluated 605 Nellore animals at 24 months of age, on average, at slaughter. The polymorphisms were determined for the molecular markers CAPN316, CAPN530, CAPN4751, CAPN4753, and UOGACAST1. Analyses of meat tenderness at 7, 14, and 21 days of maturation were performed in samples of longissimus thoracis obtained between the 12th and 13th rib and sheared using a Warner Bratzler Shear Force. Significant effects were observed for meat tenderness at days 7, 14, and 21 of maturation for the marker CAPN4751, at day 21 for the marker CAPN4753, and at days 14 and 21 for the marker UOGCAST1. For genotypic combinations of markers, the results were significant for the combination CAPN4751/UOGCAST1 in the three maturation periods and CAPN4753/UOGCAST1 at days 14 and 21 of maturation.

Implications of SNP weighting on single-step genomic predictions for different reference population sizes.

We investigated the importance of SNP weighting in populations with 2,000 to 25,000 genotyped animals. Populations were simulated with two effective sizes (20 or 100) and three numbers of QTL (10, 50 or 500). Pedigree information was available for six generations; phenotypes were recorded for the four middle generations. Animals from the last three generations were genotyped for 45,000 SNP. Single-step genomic BLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) were used to estimate genomic EBV using a genomic relationship matrix (G). The WssGBLUP performed better in small genotyped populations; however, any advantage for WssGBLUP was reduced or eliminated when more animals were genotyped. WssGBLUP had greater resolution for genome-wide association (GWA) as did increasing the number of genotyped animals. For few QTL, accuracy was greater for WssGBLUP than ssGBLUP; however, for many QTL, accuracy was the same for both methods. The largest genotyped set was used to assess the dimensionality of genomic information (number of effective SNP). The number of effective SNP was considerably less in weighted G than in unweighted G. Once the number of independent SNP is well represented in the genotyped population, the impact of SNP weighting becomes less important.

Multibreed genetic evaluation in bovines using simulated data employing a composite population.

Genetic evaluations in Brazil are performed within each animal breed; however, with the wide range of extant genetic groups in the country and the increased use of genetic crossing as a form of rapid meat production, more elaborate programs that can jointly evaluate animals of different genetic groups are needed. Genetic evaluation of a composite breed is difficult because of the variation in the genetic composition of a given herd, as well as the inclusion of non-additive genetic effects among breeds that can be important for selecting traits in certain breed combinations. Newer models include additive and non-additive effects; however, few studies have investigated these aspects in tropical breeds. The aim of this study was to simulate genetic values to compare different models. Non-inclusion of maternal effects in models leads to overestimation of variance and direct heritability. Estimates of the biological additive effects are influenced by the maternal effects; however, estimates of the non-additive effects are minimally influenced by the maternal effects and are well estimated in any situation. The studied models effectively predict the direct genetic values.

Parameters of a dynamic mechanistic model of cattle growth retain enough biological interpretation for genotype-to-phenotype mapping.

This study aimed to investigate the predictability of a phenotype when using a dynamic model of cattle growth. Genotypic and phenotypic information on Nellore (Bos indicus) cattle were used in a genome-wide association analysis designed to contrast the biological interpretation of core parameters [conversion efficiency of metabolizable energy to net energy for gain (kg) and adjusted final shrunk body weight (AFSBW)] to their associated genomic regions and nearby quantitative trait loci (QTLs). Single nucleotide polymorphisms (SNPs) were used to develop prediction equations for kg and AFSBW, which enter the model for simulative prediction purposes. QTLs and genes, one related to mature body weight and another to growth efficiency, are consistent with the model equations. Significantly associated SNPs were used to compute parameters, which yielded reasonable model outcomes when compared with regular parameter computations. Our results provide evidence of the biological validity of using such parameters as component traits of higher phenotypes and the possibility of using genomic data for genotype-to-parameter mapping.

Inclusion of cytoplasmic lineage effect and direct-maternal genetic covariance for genetic evaluation of growth traits in Nellore cattle.

We evaluated the impact of cytoplasmic lineage effects (Lc) for growth traits on genetic evaluation, including the genetic covariance between direct and maternal effects (σam). Pedigree data from 496,190 Nellore animals and observations on birth weight (BW, N = 243,391), weaning weight (WW, N = 431,681), and post-weaning weight gain adjusted to 345 days (PWG, N = 172,131) were analyzed. Four univariate models were used to obtain estimates of (co)variance components using the restricted maximum likelihood method in the BLUPF90 program. Model 1 included Lc and σam. Model 2 included Lc and σam was set to zero. Model 3 did not include Lc. Model 4 did not include Lc and σam was set to zero. These models considered the effects of the Lc as random. Phenotypic variance obtained through cytoplasmic lineage effects was determined for all traits, ranging from 0.07 to 0.15, 0.15 to 0.03, and 0.05 to 0.03% for BW, WW, and PWG, respectively, for models 1 and 2. Correlations between direct and maternal genetic components were positive for WW and negative for BW and PWG. No differences were observed for genetic parameter estimates or animal ranking with the inclusion of σam. For BW, the likelihood ratio suggested that model 1 best fits the data, while model 4 was the most appropriate for WW and PWG. Thus, these models are recommended for genetic evaluations. Despite the low magnitude of cytoplasmic lineages, this effect could predict breeding value and improve the selection of animals for BW in this Nellore population.

Systems genetics and genome-wide association approaches for analysis of feed intake, feed efficiency, and performance in beef cattle.

Feed intake, feed efficiency, and weight gain are important economic traits of beef cattle in feedlots. In the present study, we investigated the physiological processes underlying such traits from the point of view of systems genetics. Firstly, using data from 1334 Nellore (Bos indicus) cattle and 943,577 single nucleotide polymorphisms (SNPs), a genome-wide association analysis was performed for dry matter intake, average daily weight gain, feed conversion ratio, and residual feed intake with a Bayesian Lasso procedure. Genes within 50-kb SNPs, most relevant for explaining the genomic variance, were annotated and the biological processes underlying the traits were inferred from Database for Annotation, Visualization and Integrated Discovery (DAVID) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Our results indicated several putative genomic regions associated with the target phenotypes and showed that almost all genomic variances were in the SNPs located in the intergenic and intronic regions. We further identified five main metabolic pathways related to ion transport, body composition, and feed intake control, which influenced the four phenotypes simultaneously. The systems genetics approach used in this study revealed novel pathways related to feed efficiency traits in beef cattle.

Genetic variation of the weaning weight of beef cattle as a function of accumulated heat stress.

The objective of this study was to identify the genetic variation in the weaning weight (WW) of beef cattle as a function of heat stress. The WWs were recorded at approximately 205 days of age in three Brazilian beef cattle populations: Nelore (93,616), Brangus (18,906) and Tropical Composite (62,679). In view of the cumulative nature of WW, the effect of heat stress was considered as the accumulation of temperature and humidity index units (ACTHI) from the animal's birth to weaning. A reaction norm model was used to estimate the (co)variance components of WW across the ACTHI scale. The accumulation of THI units from birth to weaning negatively affected the WW. The definition of accumulated THI units as an environmental descriptor permitted to identify important genetic variation in the WW as a function of heat stress. As evidence of genotype by environment interaction, substantial heterogeneity was observed in the (co)variance components for WW across the environmental gradient. In this respect, the best animals in less stressful environments are not necessarily the best animals in more stressful environments. Furthermore, the response to selection for WW is expected to be lower in more stressful environments.

Genome-wide association with residual body weight gain in Bos indicus cattle.

Weight gain is a key performance trait for beef cat-tle; however, attention should be given to the production costs for better profitability. Therefore, a feed efficiency trait based on per-formance can be an interesting approach to improve performance without increasing food costs. To identify candidate genes and ge-nomic regions associated with residual body weight gain (RWG), we conducted a genome-wide association study (GWAS) with 720 Nellore cattle using the GRAMMAR-Gamma association test. We identified 30 significant single nucleotide polymorphisms (SNPs), especially on chromosomes 2, 8, 12, and 17. Several genes and quantitative train loci (QTLs) present in the regions identified were appointed; we highlight DMRT2 (doublesex and mab-3 related tran-scription factor 2), IFFO2 (intermediate filament family orphan 2), LNX2 (ligand of numb-protein X 2), MTIF3 (mitochondrial transla-tional initiation factor 3), and TRNAG-CCC (transfer RNA glycine anticodon CCC). The metabolic pathways that can explain part of the phenotypic variation in RWG are related to oxidative stress and muscle control.

A genomewide association mapping study using ultrasound-scanned information identifies potential genomic regions and candidate genes affecting carcass traits in Nellore cattle.

The aim of this study was to identify candidate genes and genomic regions associated with ultrasound-derived measurements of the rib-eye area (REA), backfat thickness (BFT) and rumpfat thickness (RFT) in Nellore cattle. Data from 640 Nellore steers and young bulls with genotypes for 290 863 single nucleotide polymorphisms (SNPs) were used for genomewide association mapping. Significant SNP associations were explored to find possible candidate genes related to physiological processes. Several of the significant markers detected were mapped onto functional candidate genes including ARFGAP3, CLSTN2 and DPYD for REA; OSBPL3 and SUDS3 for BFT; and RARRES1 and VEPH1 for RFT. The physiological pathway related to lipid metabolism (CLSTN2, OSBPL3, RARRES1 and VEPH1) was identified. The significant markers within previously reported QTLs reinforce the importance of the genomic regions, and the other loci offer candidate genes that have not been related to carcass traits in previous investigations.

Comparison of bivariate and multivariate joint analyses on the selection loss of beef cattle.

For genetic evaluation of beef cattle, univariate or bivariate analyses are often performed as an alternative to decrease the complexity of matrices and mathematical models compared to multivariate analysis, which considers a larger number of joint traits. The use of bivariate methods to calculate genetic predictors may cause bias in the estimation of breeding values and, as a consequence, reclassification of the rank of top-selected sires, resulting in a loss of genetic gain in future generations. The objective of this study was to compare the bivariate and multivariate joint methods of genetic evaluation, verifying the selection loss, and reclassification of the ranking of the best animals with different selection intensities. Records of 431,224 Nellore breed animals were evaluated for birth weight, weaning weight, post-weaning gain, muscle score, scrotal circumference, and selection index. The pedigree file consisted of 505,848 animals, including 218,727 males and 287,121 females. The predicted breeding values were obtained using the program PEST 2, and the complete pedigree analysis was performed by the PopReport software. The results showed that, for the four different selection intensities considered (TOP 10 and 1, 10, and 30%), selection loss and reclassification of animals in ranking, were detected for all traits evaluated when the two methods of analysis were compared.

Use of weaning management group as a random effect for a more robust estimation of genetic parameters for post-weaning traits in Nellore cattle.

Data from 69,525 animals were used to compare two types of analyses, one of them having the weaning management group (WEMANG) included as an effect in the contemporary group (F_WEMANG) and the other considering the weaning management group as a random effect, not related to the mathematical model (R_WEMANG) for post-weaning traits. The components of (co)variance were estimated for pre-weaning traits (birth weight and weaning weight) and for post-weaning traits [scrotal circumference (SC), weight gain from weaning to 18 months of age (WG) and muscle score (MUSC)] in Nellore cattle, based on a complete animal model. Heritability of SC, WG and MUSC for the F_WEMANG model was equal to 0.46 ± 0.02, 0.38 ± 0.03 and 0.26 ± 0.01, and for the R_WEMANG model it was 0.45 ± 0.02, 0.31 ± 0.03 and 0.25 ± 0.01, respectively. Genetic correlations between all the studied traits varied between 0.07 ± 0.01 and 0.77 ± 0.03 in F_WEMANG and between 0.02 ± 0.01 and 0.76 ± 0.04 in R_WEMANG. The R_ WEMANG model allowed a decrease in the number of contemporary groups as well as an increase in the number of observations per group without significant alterations in heritability coefficients, for the post-weaning traits. Consequently, the analysis became more robust and avoided having contemporary groups with low variability.

Genotype by environment interaction for post-weaning weight gain, scrotal circumference, and muscling score of composite beef cattle in different regions of Brazil.

The objectives of this study were to characterize and define homogenous production environments of composite beef cattle in Brazil in terms of climatic and geographic variables by using multivariate exploratory techniques; to evaluate the presence of genotype by environment interaction (GxE) for post-weaning weight gain (PWG), yearling scrotal circumference (SC), and yearling muscling (MUS). Hierarchical and nonhierarchical cluster analysis was used to group farms located in regions with similar environmental variables into clusters. Six clusters of farms were formed. The effect of sire-cluster interaction was tested by single-trait analysis. Genetic parameters were estimated by multi-trait analysis considering the same trait to be different in each cluster. The effect of sire-cluster interaction was significant (P <0.01) for PWG and MUS. Estimates of genetic correlations among clusters ranged from 0.31 to 0.93 for PWG, 0.64 to 0.89 for SC, and 0.18 to 0.80 for MUS. These results indicate the need for a genetic analysis on a regional basis or inclusion of the GxE effect in the statistical model to permit appropriate evaluation of the animals.

Sire effects on carcass and meat quality traits of young Nellore bulls.

Meat quality is being increasingly demanded by consumers in recent years. Several factors can affect meat quality, ranging from animal traits such as breed and genetic heritage to pre- and post-slaughter processes. This study investigated the influence of Nellore bulls on carcass and meat quality traits. We used 475 young uncastrated males, the progeny of 54 bulls, to evaluate characteristics of the following carcass traits: hot carcass weight, rib-eye area, and fat thickness. We also evaluated the following beef quality traits: marbling, color, drip loss, cooking loss, and shear force at 0, 7, and 14 days of aging. Bulls had a significant influence (P≤0.05) on rib-eye area, fat thickness, marbling, drip loss at 14 days of aging and color at all aging periods. Based on these results, the use of bulls with high breeding values for these traits can provide important advances in carcass traits and meat quality in breeding programs of Nellore cattle that are raised in tropical conditions.

Quantitative genetic study of age at subsequent rebreeding in Nellore cattle by using survival analysis.

The continuous trait age at subsequent rebreeding (ASR) was evaluated using survival analysis in Nellore breed cows that conceived for the first time at approximately 14 months of age. This methodology was chosen because the restricted breeding season produces censored data. The dataset contained 2885 records of ASR (in days). Records of females that did not produce calves in the following year after being exposed to a sire were considered censored (48.3% of the total). The statistical model used was a Weibull mixed survival model, which included fixed effects of contemporary groups (CG) and period and a random effect of individual animal. The effect of contemporary groups on ASR was significant (P < 0.01). Heritabilities obtained for ASR were 0.03 and 0.04 in logarithmic and original scales, respectively. These results indicate that the genetic selection response for subsequent reproduction of 2-year-old Nellore breed females is not expected to be effective based on survival analysis. Furthermore, these results suggest that environmental improvement is fundamental to this important trait. It should be highlighted that an increase in the average date of birth can produce an adverse effect in the future, since this cannot be compensated by genetic improvement.

Genome-wide association study for feedlot average daily gain in Nellore cattle (Bos indicus).

The genome-wide association study (GWAS) results are presented for average daily gain (ADG) in Nellore cattle. Phenotype of 720 male Bos indicus animals with information of ADG in feedlots and 354,147 single-nucleotide polymorphisms (SNPs) obtained from a database added by information from Illumina Bovine HD (777,962 SNPs) and Illumina BovineSNP50 (54,609) by imputation were used. After quality control and imputation, 290,620 SNPs remained in the association analysis, using R package Genome-wide Rapid Association using Mixed Model and Regression method GRAMMAR-Gamma. A genomic region with six significant SNPs, at Bonferroni-corrected significance, was found on chromosome 3. The most significant SNP (rs42518459, BTA3: 85849977, p = 9.49 × 10(-8)) explained 5.62% of the phenotypic variance and had the allele substitution effect of -0.269 kg/day. Important genes such as PDE4B, LEPR, CYP2J2 and FGGY are located near this region, which is overlapped by 12 quantitative trait locus (QTLs) described for several production traits. Other regions with markers with suggestive effects were identified in BTA6 and BTA10. This study showed regions with major effects on ADG in Bos indicus in feedlots. This information may be useful to increase the efficiency of selecting this trait and to understand the physiological processes involved in its regulation.

Polymorphisms in FGFBP1 and FGFBP2 genes associated with carcass and meat quality traits in chickens.

In the past, the focus of broiler breeding programs on yield and carcass traits improvement led to problems related to meat quality. Awareness of public concern for quality resulted in inclusion of meat quality traits in the evaluation process. Nevertheless, few genes associated with meat quality attributes are known. Previous studies mapped quantitative trait loci for weight at 35 and 42 days in a region of GGA4 flanked by the microsatellite markers, MCW0240 and LEI0063. In this region, there are 2 fibroblast growth factor binding protein (FGFBP) genes that play an important role in embryogenesis, cellular differentiation, and proliferation in chickens. The objective of this study was to identify and associate single nucleotide polymorphisms (SNPs) in FGFBP1 and FGFBP2 with performance, carcass, and meat quality in experimental and commercial chicken populations. In the commercial population, SNP g.2014G>A in FGFBP1 was associated with decreased carcass weight (P < 0.05), and SNP g.651G>A in FGFBP2 was associated with thawing loss and meat redness content (P < 0.05). Four haplotypes were constructed based on 2 SNPs and were associated with breast weight, thawing loss, and meat redness content. The diplotypes were associated with thawing loss, lightness, and redness content. The SNPs evaluated in the present study may be used as markers in poultry breeding programs to aid in improving growth and meat quality traits.

Identification and association of polymorphisms in CAPN1 and CAPN3 candidate genes related to performance and meat quality traits in chickens.

Meat quality is an important feature for the poultry industry and is associated with consumer satisfaction. The calpain 1 (CAPN1) gene is related to the tenderness process of meat post- mortem, and the calpain 3 (CAPN3) gene plays an important role in myofibrillar organization and growth. The objective of the present study was to identify polymorphisms in these genes and to determine the association between these polymorphisms and traits of economic interest in poultry. Eleven animals (F1) from an experimental poultry population at Embrapa Swine and Poultry were used to identify the polymorphisms. Four single nucleotide polymorphisms (SNPs) were found in the CAPN1 gene, and one SNP was found in the CAPN3 gene. A polymorphism from each gene was selected for genotyping in 152 chickens from the Embrapa F2 experimental population and 311 chickens from a commercial population. Polymorphism g.2554T>C (CAPN1) was associated with body weight at 35 to 42 days, thigh weight, breast weight, carcass weight, and meat lightness content. SNP g.15486C>T (CAPN3) was associated with thigh yield, thawing-cooking loss, and shear force. Results suggest the possibility of using molecular markers in CAPN1 and CAPN3 genes as a tool for performance and meat quality traits in poultry breeding programs.

Protein synthesis and degradation gene SNPs related to feed intake, feed efficiency, growth, and ultrasound carcass traits in Nellore cattle.

We looked for possible associations of SNPs in genes related to protein turnover, with growth, feed efficiency and carcass traits in feedlot Nellore cattle. Purebred Nellore bulls and steers (N = 290; 378 ± 42 kg body weight, 23 months ± 42 days old) were evaluated for daily feed intake, body weight gain (BWG), gross feed efficiency, feed conversion ratio, partial efficiency of growth, residual feed intake (RFI), ultrasound backfat, rump fat, and ribeye area. Genotypes were obtained for SNPs in the growth hormone receptor (GHR-1 and GHR-2); calpain (CAPN4751); calpastatin (UoGCAST); ubiquitin-conjugating enzyme 2I (UBE2I-1 and UBE2I-2); R3H domain containing 1 (R3HDM1-1, -2, -3, and -4), ring finger protein 19 (RNF19); proteasome 26S subunit, non-ATPase, 13 (PSMD13); ribosomal protein, large, P2 (RPLP2); and isoleucine-tRNA synthetase 2, mitochondrial (IARS2) genes. Allelic substitution, additive and dominant effects were tested and molecular breeding values were computed. CAPN4751, GHR-1 and -2, IARS2, R3HDM1-4, and UoGCAST were found to be normally segregating polymorphisms. Additive and dominance effects were observed on BWG, feed efficiency and carcass traits, although dominant effects predominated. Significant allelic substitution effects were observed for CAPN4751, GHR-1 and -2, and UoGCAST on BWG, gross feed efficiency, RFI, and carcass traits, under single- or multiple-marker analyses. Correlations between molecular breeding values and phenotypes were low, excepted for RFI, based on allelic substitution estimates obtained by stepwise linear regression. We conclude that SNPs in genes related to protein turnover are related to economically important traits in Nellore cattle.