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 effects of heat stress on milk fatty acids in Brazilian Holstein cattle.

The present study aimed to estimate covariance components of milk fatty acids (FA) and to compare the genomic estimated breeding values under general and heat-stress effects. Data consisted of 38,762 test-day records from 6,344 Holstein cows obtained from May 2012 through January 2018 on 4 dairy herds from Brazil. Single-trait repeatability test-day models with random regressions as a function of temperature-humidity index values were used for genetic analyses. The models included contemporary groups, parity order (1-6), and days in milk classes as fixed effects, and general and thermotolerance additive genetic and permanent environmental as random effects. Notably, differences in heritability estimates between environments (general and heat stress) increased (0.03 to 0.06) for unsaturated FA traits, such as unsaturated, monounsaturated, and polyunsaturated, at higher heat-stress levels. In contrast, heritability estimated between environments for saturated FA traits, including saturated FA, palmitic acid (C16:0), and stearic acid (C18:0) did not observe significant differences between environments. In addition, our study revealed negative genetic correlations between general and heat-stress additive genetic effects (antagonistic effect) for the saturated FA, C16:0, C18:0, and C18:1, which ranged from -0.007 to -0.32. Spearman's ranking correlation between genomic estimated breeding values ranged from -0.27 to 0.99. Results indicated a moderate to strong interaction of genotype by the environment for most FA traits comparing a heat-stress environment with thermoneutral conditions. Our findings point out novel opportunities to explore the use of FA milk profile and heat-stress models.

Genome-wide association study reveals genes associated with the absence of intermuscular bones in tambaqui (Colossoma macropomum).

The presence of intermuscular bones in fisheries products limits the consumption and commercialization potential of many fish species, including tambaqui (Colossoma macropomum). These bones have caused medical emergencies and are an undesirable characteristic for fish farming because their removal is labor-intensive during fish processing. Despite the difficulty in identifying genes related to the lack of intermuscular bone in diverse species of fish, the discovery of individuals lacking intermuscular bones in a Neotropical freshwater characiform fish has provided a unique opportunity to delve into the genetic mechanisms underlying the pathways of intermuscular bone formation. In this study, we carried out a GWAS among boneless and wt tambaqui populations to identify markers associated with a lack of intermuscular bone. After analyzing 11 416 SNPs in 360 individuals (12 boneless and 348 bony), we report 675 significant (Padj  < 0.003) associations for this trait. Of those, 13 associations were located near candidate genes related to the reduction of bone mass, promotion of bone formation, inhibition of bone resorption, central control of bone remodeling, bone mineralization and other related functions. To the best of our knowledge, for the first time, we have successfully identified genes related to a lack of intermuscular bones using GWAS in a non-model species.

Genome-wide association study for milk production traits in a Brazilian Holstein population.

Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.

Estimates of genomic heritability and genome-wide association study for fatty acids profile in Santa Inês sheep.

BACKGROUND: Despite the health concerns and nutritional importance of fatty acids, there is a relative paucity of studies in the literature that report genetic or genomic parameters, especially in the case of sheep populations. To investigate the genetic architecture of fatty acid composition of sheep, we conducted genome-wide association studies (GWAS) and estimated genomic heritabilities for fatty acid profile in Longissimus dorsi muscle of 216 male sheep. RESULTS: Genomic heritability estimates for fatty acid content ranged from 0.25 to 0.46, indicating that substantial genetic variation exists for the evaluated traits. Therefore, it is possible to alter fatty acid profiles through selection. Twenty-seven genomic regions of 10 adjacent SNPs associated with fatty acids composition were identified on chromosomes 1, 2, 3, 5, 8, 12, 14, 15, 16, 17, and 18, each explaining ≥0.30% of the additive genetic variance. Twenty-three genes supporting the understanding of genetic mechanisms of fat composition in sheep were identified in these regions, such as DGAT2, TRHDE, TPH2, ME1, C6, C7, UBE3D, PARP14, and MRPS30. CONCLUSIONS: Estimates of genomic heritabilities and elucidating important genomic regions can contribute to a better understanding of the genetic control of fatty acid deposition and improve the selection strategies to enhance meat quality and health attributes.

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.

Variation in myogenic differentiation 1 mRNA abundance is associated with beef tenderness in Nelore cattle.

The myogenic differentiation 1 gene (MYOD1) has a key role in skeletal muscle differentiation and composition through its regulation of the expression of several muscle-specific genes. We first used a general linear mixed model approach to evaluate the association of MYOD1 expression levels on individual beef tenderness phenotypes. MYOD1mRNA levels measured by quantitative polymerase chain reactions in 136 Nelore steers were significantly associated (P ≤ 0.01) with Warner-Bratzler shear force, measured on the longissimus dorsi muscle after 7 and 14 days of beef aging. Transcript abundance for the muscle regulatory gene MYOD1 was lower in animals with more tender beef. We also performed a co-expression network analysis using whole transcriptome sequence data generated from 30 samples of longissimus muscle tissue to identify genes that are potentially regulated by MYOD1. The effect of MYOD1 gene expression on beef tenderness may emerge from its function as an activator of muscle-specific gene transcription such as for the serum response factor (C-fos serum response element-binding transcription factor) gene (SRF), which determines muscle tissue development, composition, growth and maturation.

A single nucleotide polymorphism in NEUROD1 is associated with production traits in Nelore beef cattle.

Feed efficiency and carcass characteristics are late-measured traits. The detection of molecular markers associated with them can help breeding programs to select animals early in life, and to predict breeding values with high accuracy. The objective of this study was to identify polymorphisms in the functional and positional candidate gene NEUROD1 (neurogenic differentiation 1), and investigate their associations with production traits in reference families of Nelore cattle. A total of 585 steers were used, from 34 sires chosen to represent the variability of this breed. By sequencing 14 animals with extreme residual feed intake (RFI) values, seven single nucleotide polymorphisms (SNPs) in NEUROD1 were identified. The investigation of marker effects on the target traits RFI, backfat thickness (BFT), ribeye area (REA), average body weight (ABW), and metabolic body weight (MBW) was performed with a mixed model using the restricted maximum likelihood method. SNP1062, which changes cytosine for guanine, had no significant association with RFI or REA. However, we found an additive effect on ABW (P ≤ 0.05) and MBW (P ≤ 0.05), with an estimated allele substitution effect of -1.59 and -0.93 kg0.75, respectively. A dominant effect of this SNP for BFT was also found (P ≤ 0.010). Our results are the first that identify NEUROD1 as a candidate that affects BFT, ABW, and MBW. Once confirmed, the inclusion of this SNP in dense panels may improve the accuracy of genomic selection for these traits in Nelore beef cattle as this SNP is not currently represented on SNP chips.

Allele- and parent-of-origin-specific effects on expression of the KCNJ11 gene: A candidate for meat tenderness in cattle.

In contrast to the Mendelian inheritance model, parental alleles can contribute unequally to gene expression, which may result in phenotypic variance among individuals and bias in the predicted additive effect of molecular markers associated with production traits. Given the need to understand the effects of allelic variation and parent-of-origin effects on the expression of genes with a commercial interest in cattle, we analyzed the expression of KCNJ11 (potassium inwardly rectifying channel, subfamily J, member 11), which was previously described as a functional candidate gene for meat tenderness. Allele-specific and parent-of-origin-dependent expression of this gene were assessed in bovine muscle using the rs379610823 single nucleotide polymorphism as a reference. Biallelic expression was observed; however, the T allele was expressed at significantly higher levels than the C allele. Furthermore, increased expression of KCNJ11 was found in animals harboring the maternal T allele. This study is the first to describe the differential allelic expression of bovine KCNJ11. Our findings are important for understanding the mechanisms that underlie the pattern of KCNJ11 expression and its potential impact on the phenotypic variation of meat tenderness in Nelore beef cattle. This reinforces the need for further investigation of allelic- and parent-of-origin expression deviation in genetic markers eligible for the selection of target traits.

Genetic parameters for milk fatty acids, milk yield and quality traits of a Holstein cattle population reared under tropical conditions.

Information about genetic parameters is essential for selection decisions and genetic evaluation. These estimates are population specific; however, there are few studies with dairy cattle populations reared under tropical and sub-tropical conditions. Thus, the aim was to obtain estimates of heritability and genetic correlations for milk yield and quality traits using pedigree and genomic information from a Holstein population maintained in a tropical environment. Phenotypic records (n = 36 457) of 4203 cows as well as the genotypes for 57 368 single nucleotide polymorphisms from 755 of these cows were used. Covariance components were estimated using the restricted maximum likelihood method under a mixed animal model, considering a pedigree-based relationship matrix or a combined pedigree-genomic matrix. High heritabilities (around 0.30) were estimated for lactose and protein content in milk whereas moderate values (between 0.19 and 0.26) were obtained for percentages of fat, saturated fatty acids and palmitic acid in milk. Genetic correlations ranging from -0.38 to -0.13 were determined between milk yield and composition traits. The smaller estimates compared to other similar studies can be due to poor environmental conditions, which may reduce genetic variability. These results highlight the importance in using genetic parameters estimated in the population under evaluation for selection decisions.

The Receptive Endometrial Transcriptomic Signature Indicates an Earlier Shift from Proliferation to Metabolism at Early Diestrus in the Cow.

This study aimed to characterize the endometrial transcriptome and functional pathways overrepresented in the endometrium of cows treated to ovulate larger (≥13 mm) versus smaller (≤12 mm) follicles. Nelore cows were presynchronized prior to receiving cloprostenol (large follicle [LF] group) or not (small follicle [SF] group), along with a progesterone (P4) device on Day (D) -10. Devices were withdrawn and cloprostenol administered 42-60 h (LF) or 30-36 h (SF) before GnRH agonist treatment (D0). Tissues were collected on D4 (experiment [Exp.] 1; n = 24) or D7 (Exp. 2; n = 60). Endometrial transcriptome was obtained by RNA-Seq, whereas proliferation and apoptosis were assessed by immunohistochemistry. Overall, LF cows developed larger follicles and corpora lutea, and produced greater amounts of estradiol (D-1, Exp. 1, SF: 0.7 ± 0.2; LF: 2.4 ± 0.2 pg/ml; D-1, Exp. 2, SF: 0.5 ± 0.1; LF: 2.3 ± 0.6 pg/ml) and P4 (D4, Exp. 1, SF: 0.8 ± 0.1; LF: 1.4 ± 0.2 ng/ml; D7, Exp. 2, SF: 2.5 ± 0.4; LF: 3.7 ± 0.4 ng/ml). Functional enrichment indicated that biosynthetic and metabolic processes were enriched in LF endometrium, whereas SF endometrium transcriptome was biased toward cell proliferation. Data also suggested reorganization of the extracellular matrix toward a proliferation-permissive phenotype in SF endometrium. LF endometrium showed an earlier onset of proliferative activity, whereas SF endometrium expressed a delayed increase in glandular epithelium proliferation. In conclusion, the periovulatory endocrine milieu regulates bovine endometrial transcriptome and seems to determine the transition from a proliferation-permissive to a biosynthetic and metabolically active endometrial phenotype, which may be associated with the preparation of an optimally receptive uterine environment.

SNP and INDEL detection in a QTL region on chicken chromosome 2 associated with muscle deposition.

Genetic improvement is important for the poultry industry, contributing to increased efficiency of meat production and quality. Because breast muscle is the most valuable part of the chicken carcass, knowledge of polymorphisms influencing this trait can help breeding programs. Therefore, the complete genome of 18 chickens from two different experimental lines (broiler and layer) from EMBRAPA was sequenced, and SNPs and INDELs were detected in a QTL region for breast muscle deposition on chicken chromosome 2 between microsatellite markers MCW0185 and MCW0264 (105,849-112,649 kb). Initially, 94,674 unique SNPs and 10,448 unique INDELs were identified in the target region. After quality filtration, 77% of the SNPs (85,765) and 60% of the INDELs (7828) were retained. The studied region contains 66 genes, and functional annotation of the filtered variants identified 517 SNPs and three INDELs in exonic regions. Of these, 357 SNPs were classified as synonymous, 153 as non-synonymous, three as stopgain, four INDELs as frameshift and three INDELs as non-frameshift. These exonic mutations were identified in 37 of the 66 genes from the target region, three of which are related to muscle development (DTNA, RB1CC1 and MOS). Fifteen non-tolerated SNPs were detected in several genes (MEP1B, PRKDC, NSMAF, TRAPPC8, SDR16C5, CHD7, ST18 and RB1CC1). These loss-of-function and exonic variants present in genes related to muscle development can be considered candidate variants for further studies in chickens. Further association studies should be performed with these candidate mutations as should validation in commercial populations to allow a better explanation of QTL effects.

Variant discovery in a QTL region on chromosome 3 associated with fatness in chickens.

Abdominal fat content is an economically important trait in commercially bred chickens. Although many quantitative trait loci (QTL) related to fat deposition have been detected, the resolution for these regions is low and functional variants are still unknown. The current study was conducted aiming at increasing resolution for a region previously shown to have a QTL associated with fat deposition, to detect novel variants from this region and to annotate those variants to delineate potentially functional ones as candidates for future studies. To achieve this, 18 chickens from a parental generation used in a reciprocal cross between broiler and layer lines were sequenced using the Illumina next-generation platform with an initial coverage of 18X/chicken. The discovery of genetic variants was performed in a QTL region located on chromosome 3 between microsatellite markers LEI0161 and ADL0371 (33,595,706-42,632,651 bp). A total of 136,054 unique SNPs and 15,496 unique INDELs were detected in this region, and after quality filtering, 123,985 SNPs and 11,298 INDELs were retained. Of these variants, 386 SNPs and 15 INDELs were located in coding regions of genes related to important metabolic pathways. Loss-of-function variants were identified in several genes, and six of those, namely LOC771163, EGLN1, GNPAT, FAM120B, THBS2 and GGPS1, were related to fat deposition. Therefore, these loss-of-function variants are candidate mutations for conducting further studies on this important trait in chickens.

Identification of polymorphisms associated with production traits on chicken (Gallus gallus) chromosome 4.

Genetic selection for production traits has resulted in a rapid improvement in animal performance and development. Previous studies have mapped quantitative trait loci for body weight at 35 and 41 days, and drum and thigh yield, onto chicken chromosome 4. We investigated this region for single nucleotide polymorphisms and their associations with important economic traits. Three positional candidate genes were studied: KLF3 (Krüeppel-like factor 3), SLIT2 (Slit homolog 2), and PPARGC1A (peroxisome proliferator-activated receptor gamma, coactivator 1 alpha). Fragment sequencing of these genes was conducted in 11 F1 animals, and one polymorphism in each gene was selected and genotyped in an F2 population (N = 276) and a paternal broiler line TT (N = 840). Associations were identified with growth, carcass, and fat traits in the F2 and the paternal line (P < 0.05). Using single markers in both the F2 and the TT line, KLF3 was associated with weight gain (P < 0.05), PPPARGC1A was associated with liver and wing-parts weights and yields (P < 0.05), and SLIT2 was associated with back yield (P < 0.05) and fat traits (P < 0.05). Using multiple markers, KLF3 lost its significance in both populations, and SLIT2 was associated with feed conversion only in the TT population (P < 0.05). The QTLs mapped in the F2 population could be partly explained by PPARGC1A and SLIT2, which were associated with body weight at 35 and 41 days, respectively, and with drum and thigh yield in the same population. The results of this study indicate the importance of these genes for production traits.

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.

A genome-wide scan for selection signatures in Nellore cattle.

Brazilian Nellore cattle (Bos indicus) have been selected for growth traits for over more than four decades. In recent years, reproductive and meat quality traits have become more important because of increasing consumption, exports and consumer demand. The identification of genome regions altered by artificial selection can potentially permit a better understanding of the biology of specific phenotypes that are useful for the development of tools designed to increase selection efficiency. Therefore, the aims of this study were to detect evidence of recent selection signatures in Nellore cattle using extended haplotype homozygosity methodology and BovineHD marker genotypes (>777,000 single nucleotide polymorphisms) as well as to identify corresponding genes underlying these signals. Thirty-one significant regions (P < 0.0001) of possible recent selection signatures were detected, and 19 of these overlapped quantitative trait loci related to reproductive traits, growth, feed efficiency, meat quality, fatty acid profiles and immunity. In addition, 545 genes were identified in regions harboring selection signatures. Within this group, 58 genes were associated with growth, muscle and adipose tissue metabolism, reproductive traits or the immune system. Using relative extended haplotype homozygosity to analyze high-density single nucleotide polymorphism marker data allowed for the identification of regions potentially under artificial selection pressure in the Nellore genome, which might be used to better understand autozygosity and the effects of selection on the Nellore genome.

Response of resistant and susceptible Brazilian Somalis crossbreed sheep naturally infected by Haemonchus contortus.

This study was carried out to evaluate the performance of Brazilian Somalis sheep to natural infections by gastrointestinal nematodes. During 98 days, 75 weaned sheep, initially 3-4 months old, were kept on the same pasture and evaluated. Fecal and blood samples were collected for parasitological and hematological exams. After this period, the eight most resistant and the eight most susceptible animals were selected based on their individual averages of nematode fecal egg counts and were slaughtered for worm burden determination and nematodes identification. Abomasum and abomasum lymph nodes were also recovered for gene expression analysis. The animals selected as resistant had lower fecal egg counts during experimental period and smaller worm burdens than the susceptible ones (P < 0.05). The genus Haemonchus, followed by Trischostrongylus and Oesophagostomum, were identified in composite cultures. Haemonchus contortus was the specie identified in the abomasum. Packed cell volume and total plasma protein means were higher in the resistant group (27.2% and 6.1 g/dL) than in the susceptible one (22.5% and 5.3 g/dL), respectively. Regarding cytokine gene expression, IL-4 (P < 0.05) was up-regulated in the abomasum of resistant animals and TNF-α (P < 0.03) and IFN-γ (P < 0.03) in susceptible ones. In abomasum lymph nodes, IL-4 (P < 0.04) and IL-13 (P < 0.05) were up-regulated in the resistant animals and IFN-γ in the susceptible one (P < 0.01). This work provides further evidence that, within a given animal breed, individuals have different responses when infected by gastrointestinal nematodes. Resistant animals who responded more quickly and efficiently to these infections activated a TH2-type response.

Identification and characterization of microRNAs expressed in chicken skeletal muscle.

MicroRNAs (miRNAs, miRs) encompass a class of small non-coding RNAs that often negatively regulate gene expression. miRNAs play an essential role in skeletal muscle, determining the proper development and maintenance of this tissue. In comparison to other organs and tissues, the full set of muscle miRNAs and its expression patterns are still poorly understood. In this report, a chicken skeletal muscle miRNA library was constructed, and the expression of selected miRNAs was further characterized during muscle development in chicken lines with distinct muscling phenotypes. Clone library sequence analysis revealed 40 small RNAs with similarities to previously described chicken miRNAs, seven miRNAs that were never identified before in chicken, and some sequence clusters representing other possible novel miRNAs. Temporal expression profiles of three miRNAs associated with cell proliferation and differentiation (miR-125b, miR-221, and miR-206) in two chicken lines (broiler and layer) revealed the differential steady-state levels of these miRs during skeletal muscle growth and suggests that miR-206 is involved in the muscling phenotype that is observed in growth-selected chicken lines.

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.