Unknown parent groups and truncated pedigree in single-step genomic evaluations of Murrah buffaloes.

Missing pedigree may produce bias in genomic evaluations. Thus, strategies to deal with this problem have been proposed as using unknown parent groups (UPG) or truncated pedigrees. The aim of this study was to investigate the impact of modeling missing pedigree under ssGBLUP evaluations for productive and reproductive traits in dairy buffalos using different approaches: 1) traditional BLUP without UPG (BLUP), 2) traditional BLUP including UPG (BLUP/UPG), 3) ssGBLUP without UPG (ssGBLUP), 4) ssGBLUP including UPG in the A and A22 matrices (ssGBLUP/A_UPG), 5) ssGBLUP including UPG in all elements of the H matrix (ssGBLUP/H_UPG), 6) BLUP with pedigree truncation for the last 3 generations (BLUP/truncated), and 7) ssGBLUP with pedigree truncation for the last 3 generations (ssGBLUP/ truncated). UPGs were not used in the scenarios with truncated pedigree. A total of 3,717, 4,126 and 3,823 records of the first lactation for accumulated 305 d milk yield (MY), age at first calving (AFC) and lactation length (LL), respectively were used. Accuracies ranged from 0.27 for LL (BLUP) to 0.46 for MY (BLUP), bias ranged from -0.62 for MY (ssGBLUP) to 0.0002 for AFC (BLUP/truncated), and dispersion ranged from 0.88 for MY (BLUP/ A_UPG) to 1.13 for LL (BLUP). Genetic trend showed genetic gains for all traits across 20 years of selection and the impact of including either genomic information, UPG or pedigree truncation under GEBV accuracies ranged among the evaluated traits. Overall, methods using UPGs, truncation pedigree and genomic information exhibited potential to improve GEBV accuracies, bias and dispersion for all traits compared with other methods. Truncated scenarios promoted high genetic gains. In small populations with few genotyped animals, combining truncated pedigree or UPG with genomic information is a feasible approach to deal with missing pedigrees.

Application of proteomic to investigate the different degrees of meat tenderness in Nellore breed.

This study describes the association between meat tenderness and abundance of soluble muscle proteins in Nellore bulls (Bos indicus) using a proteomic approach. We evaluated shear force (SF) of Longissimus thoracis muscle 24 h after slaughter and selected three experimental groups of animals with moderately tender (TE; SF = 3.9 ± 0.7 kg), moderately tough (TO; SF = 5.6 ± 0.7 kg) and very tough meat (TO+; SF = 7.9 ± 1.4 kg). Proteome was investigated by two-dimensional electrophoresis (2D-PAGE) in combination with electrospray ionization-tandem mass spectrometry (ESI-MS/MS). The metabolic proteins triosephosphate isomerase (TPI1) and phosphoglucomutase 1 (PGM1), the structural protein profilin 1 (PFN1), and cytosol aminopeptidase (LAP3) were up-regulated (P < 0.05) in the TE meat group when compared to the TO and TO+ groups. Actin structural proteins (ACTA1, ACTB, and ACTG1), the oxidative stress protein peroxiredoxin (PRDX6, PRDX2, PRDX1, and PARK7), heat shock protein isoforms, and co-chaperones (CDC37 and STIP1) were up-regulated (P < 0.05) in the TO and TO+ meat groups. In addition, we also identified proteins PFN1, LAP3, PRDX1, PRDX2, HSPD1, and ARHGDIA to be associated with beef tenderness. The results reported herein demonstrated that meat tenderness in Nellore cattle depends on the modulation and expression of a set of proteins involved in different biological pathways. SIGNIFICANCE: The manuscript entitled "Application of proteomic to investigate the different degrees of meat tenderness in Nellore breed" describes a classical proteomics work using two-dimensional gel electrophoresis (2D-PAGE), followed by mass spectrometry coupled to electrospray ionization ion trap (ESI-MS/MS) in order to understand the biochemical engineering involved in the process of meat tenderness. We evaluated shear force (SF) of Longissimus thoracis muscle samples of Nellore cattle (n = 90) and select three experimental groups of animals with moderately tender (TE; SF = 3.9 ± 0.7), moderately tough (TO; SF = 5.6 ± 0.7) and very tough meat (TO+; SF = 7.9 ± 1.4). The proteomic approach allowed observing that meat tenderness is influenced by structural proteins (ACTA1, ACTG1, ACTB, MYL1 and PFN1), co-chaperones (CDC37 and STIP1), heat shock proteins (HSP90AA1, HSP90AB1, HSPD1, HSPA1L, HSPA1A and HSPB1), regulatory protein (ARHGDIA), metabolic proteins (TPI1 and PGM1) and oxidative stress proteins (PRDX1, PRDX2, PRDX6, PARK7). Our results suggest that meat tenderness in Nellore depends on the modulation and expression of a set of proteins involved in different biological pathways.

Genome-wide association study for growth traits in Nelore cattle.

The objective of this study was to investigate the association of single nucleotide polymorphisms (SNPs) with birth weight, weight gain from birth to weaning and from weaning to yearling, yearling height and cow weight in Nelore cattle. Data from 5064 animals participating in the DeltaGen and PAINT breeding programs were used. The animals were genotyped with a panel of 777 962 SNPs (Illumina BovineHD BeadChip) and 412 993 SNPs remained after quality control analysis of the genomic data. A genome-wide association study was performed using a single-step methodology. The analyses were processed with the BLUPF90 family of programs. When applied to a genome-wide association studies, the single-step GBLUP methodology is an iterative process that estimates weights for the SNPs. The weights of SNPs were included in all analyses by iteratively applying the single-step GBLUP methodology and repeated twice so that the effect of the SNP and the effect of the animal were recalculated in order to increase the weight of SNPs with large effects and to reduce the weight of those with small effects. The genome-wide association results are reported based on the proportion of variance explained by windows of 50 adjacent SNPs. Considering the two iterations, only windows with an additive genetic variance >1.5% were presented in the results. Associations were observed with birth weight on BTA 14, with weight gain from birth to weaning on BTA 5 and 29, with weight gain from weaning to yearling on BTA 11, and with yearling height on BTA 8, showing the genes TMEM68 (transmembrane protein 8B) associated with birth weight and yearling height, XKR4 (XK, Kell blood group complex subunit-related family, member 4) associated with birth weight, NPR2 (natriuretic peptide receptor B) associated with yearling height, and REG3G (regenerating islet-derived 3-gamma) associated with weight gain from weaning to yearling. These genes play an important role in feed intake, weight gain and the regulation of skeletal growth.

SNP detection using RNA-sequences of candidate genes associated with puberty in cattle.

Fertility traits, such as heifer pregnancy, are economically important in cattle production systems, and are therefore, used in genetic selection programs. The aim of this study was to identify single nucleotide polymorphisms (SNPs) using RNA-sequencing (RNA-Seq) data from ovary, uterus, endometrium, pituitary gland, hypothalamus, liver, longissimus dorsi muscle, and adipose tissue in 62 candidate genes associated with heifer puberty in cattle. RNA-Seq reads were assembled to the bovine reference genome (UMD 3.1.1) and analyzed in five cattle breeds; Brangus, Brahman, Nellore, Angus, and Holstein. Two approaches used the Brangus data for SNP discovery 1) pooling all samples, and 2) within each individual sample. These approaches revealed 1157 SNPs. These were compared with those identified in the pooled samples of the other breeds. Overall, 172 SNPs within 13 genes (CPNE5, FAM19A4, FOXN4, KLF1, LOC777593, MGC157266, NEBL, NRXN3, PEPT-1, PPP3CA, SCG5, TSG101, and TSHR) were concordant in the five breeds. Using Ensembl's Variant Effector Predictor, we determined that 12% of SNPs were in exons (71% synonymous, 29% nonsynonymous), 1% were in untranslated regions (UTRs), 86% were in introns, and 1% were in intergenic regions. Since these SNPs were discovered in RNA, the variants were predicted to be within exons or UTRs. Overall, 160 novel transcripts in 42 candidate genes and five novel genes overlapping five candidate genes were observed. In conclusion, 1157 SNPs were identified in 62 candidate genes associated with puberty in Brangus cattle, of which, 172 were concordant in the five cattle breeds. Novel transcripts and genes were also identified.

Physical characteristics and fertility of fractionated donkey semen cooled at 5°C / Características físicas e fertilidade do sêmen asinino fracionado e resfriado a 5ºC

The aim of this study was to evaluate the effect of two different extenders (Skimmed Milk Glucose - SMG or Lactose - Egg Yolk - LEY) on physical characteristics and fertility of fractionated donkey semen cooled at 5°C. For this, four Pêga donkeys were used as semen donors. The sperm rich fraction of the ejaculate was diluted preparing insemination doses containing 400 x 106 motile spermatozoa in a volume of 22 mL, cooled to 5°C and stored up to 48 hours in a container proposed by Palhares (1997). Sperm motility and vigor were assessed in fresh semen, after first semen dilution, before insemination, at 24 and 48 hours after storage. For the fertility evaluation, 44 mares were inseminated with semen stored for a period between 12 and 24 hours. The mares were inseminated on fixed days (Mondays, Wednesdays, and Fridays) after the detection of a follicle greater than a 30mm diameter in one of the ovaries through ovulation. Pregnancy diagnosis was performed on day 12 post-ovulation, using transrectal ultrasonography. Semen diluted in SMG showed superior sperm motility than LEY, at the Pre-AI evaluation (P<0.05). At 48 hours of storage, all donkeys had motility values between 45 and 53% for semen diluted in SMG, while only one donkey showed motility greater than 30% in the LEY treatment. The pregnancy rate/cycle for mares inseminated with semen diluted in SMG was superior than that obtained using LEY (56.52% vs 4.76%, respectively).(AU)

Objetivou-se com o presente experimento avaliar o efeito de dois diferentes diluidores (leite em pó desnatado glicose - SMG ou lactose gema de ovo - LEY) sobre as características físicas e a fertilidade do sêmen asinino coletado de forma fracionada e resfriado a 5ºC. Para isso, quatro jumentos da raça Pêga foram utilizados como doadores de sêmen. A fração espermática rica do ejaculado foi diluída preparando-se doses inseminantes contendo 400 x 106 espermatozoides móveis em um volume de 22 mL, resfriadas a 5ºC e armazenadas por até 48 horas em contêiner proposto por Palhares (1997). A motilidade e o vigor espermáticos foram avaliados no sêmen fresco, após a pré-diluição, antes das inseminações, às 24 e 48 horas de armazenamento. Para avaliação de fertilidade, 44 éguas foram inseminadas com sêmen armazenado por um período entre 12 e 24 horas, em dias fixos (segundas, quartas e sextas-feiras), após a detecção de um folículo de diâmetro maior ou igual a 30mm em um dos ovários, até a ovulação. O diagnóstico de gestação foi realizado a partir de 12 dias após a ovulação, por meio de ultrassonografia transretal. O sêmen diluído em SMG apresentou motilidade espermática superior à do LEY, já a partir do tempo pré-IA. Às 48 horas de armazenamento, todos os jumentos apresentaram valores de motilidade entre 45% e 53%, quando o sêmen foi diluído em SMG, enquanto apenas um jumento apresentou motilidade superior a 30% no tratamento utilizando LEY. A taxa de concepção/ciclo das éguas inseminadas também foi superior para o sêmen diluído em SMG em relação ao diluído em LEY (56,52% versus 4,76%, respectivamente).(AU)

Genome scan for postmortem carcass traits in Nellore cattle.

Carcass traits measured after slaughter are economically relevant traits in beef cattle. In general, the slaughter house payment system is based on HCW. Ribeye area (REA) is associated with the amount of the meat in the carcass, and a minimum of backfat thickness (BFT) is necessary to protect the carcass during cooling. The aim of this study was to identify potential genomic regions harboring candidate genes affecting those traits in Nellore cattle. The data set used in the present study consisted of 1,756 Nellore males with phenotype records. A subset of 1,604 animals had both genotypic and phenotypic information. Genotypes were generated based on a panel with 777,962 SNPs from the Illumina Bovine HD chip. The SNP effects were calculated based on the genomic breeding values obtained by using the single-step GBLUP approach and a genomic matrix re-weighting procedure. The proportion of the variance explained by moving windows of 100 consecutive SNPs was used to assess potential genomic regions harboring genes with major effects on each trait. The top 10 non-overlapping SNP-windows explained 8.72%, 11.38%, and 9.31% of the genetic variance for REA, BFT, and HCW, respectively. These windows are located on chromosomes 5, 7, 8, 10, 12, 20, and 29 for REA; chromosomes 6, 8, 10, 13, 16, 17, 18, and 24 for BFT; and chromosomes 4, 6, 7, 8, 14, 16, 17, and 21 for HCW. For REA, there were identified genes ( and ) involved in the cell cycle biological process which affects many aspects of animal growth and development. The and genes, both from AA transporter family, was also associated with REA. The AA transporters are essential for cell growth and proliferation, acting as carriers of tissue nutrient supplies. Various genes identified for BFT (, , , , , and ) have been associated with lipid metabolism in different mammal species. One of the most promising genes identified for HCW was the . There is evidence, in the literature, that this gene is located in putative QTL affecting carcass weight in beef cattle. Our results showed several genomic regions containing plausible candidate genes that may be associated with carcass traits in Nellore cattle. Besides contributing to a better understanding of the genetic control of carcass traits, the identified genes can also be helpful for further functional genomic studies.

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.

Reproductive biology of Pêga Jackasses: characteristics of sperm-rich fraction collected using the split-ejaculate method / Biologia reprodutiva de jumentos da raça Pêga: características da fração rica em espermatozoides coletada de forma fracionada

The aim of this study was to evaluate the physical and morphological characteristics of the sperm-rich fraction of jackass semen. To this end, 130 ejaculates from five Pêga jackasses were collected using an open model artificial vagina. The sperm-rich fraction was collected using the split-ejaculate method and assessed for the number of mounts/ejaculate, for physical and morphological characteristics of the semen and number of doses produced/ejaculate. It was observed that all characteristics evaluated differed among the five jackasses, except for the head defect rates. The mean values obtained for the jackass sperm-rich fraction collected were: number of mounts/ejaculate - 1.27; semen volume - 20.21mL; motility - 84.53%; vigor - 4.46; motility after dilution - 80.10%; sperm concentration/mL - 894.38 x 10 6; total sperm/ejaculate - 16.14 x 10 9; number of insemination doses/ejaculate (400x10 6 motile sptz) - 33.39; number of insemination doses/ejaculate (800 x 10 6 motile sptz) - 16.69; and percentage of normal sperm - 90.46%. Thus, in the present experiment the split-ejaculate method using an open artificial vagina worked well with the jackasses, and the sperm-rich fraction of the ejaculate of Pêga jackasses had high quality and sperm concentration, allowing its use for semen processing without reducing the number of insemination doses produced per ejaculate.(AU)

Foram coletados 130 ejaculados de cinco jumentos da raça Pêga, utilizando-se vagina artificial modelo aberta, com o objetivo de se avaliar as características físicas e morfológicas da fração rica em espermatozoides do ejaculado de asininos. Para tal, a fração rica em espermatozoides, composta pelos três primeiros jatos ejaculados, foi coletada e avaliada quanto ao número de montas/ejaculado, quanto às características físicas e morfológicas do sêmen, bem como quanto ao número de doses inseminantes produzidas/ejaculado. Observou-se que todas as características avaliadas diferiram entre os cinco reprodutores avaliados, com exceção do percentual de defeitos de cabeça. Os valores médios obtidos da coleta da fração rica do ejaculado de jumentos foram: número de montas/ejaculado - 1,27; volume de sêmen - 20,21mL; motilidade - 84,53%; vigor - 4,46; motilidade pós-diluição - 80,10%; espermatozoides/mL - 894,38 x 10 6; espermatozoides/ejaculado - 16,14 x 10 9; número de doses inseminantes/ejaculado (400x10 6 sptz móveis) - 33,39; doses inseminantes/ejaculado (800 x 10 6 sptz/móveis) - 16,69; e percentual de espermatozoides normais - 90,46%. Assim, no presente experimento, observou-se boa aceitação dos reprodutores à coleta fracionada utilizando-se a vagina artificial aberta, sendo a fração rica do ejaculado de jumentos da raça Pêga caracterizada por alta qualidade e concentração espermática, o que viabilizou sua utilização para o processamento do sêmen, sem prejuízos quanto ao número de doses inseminantes produzidas/ejaculado.(AU)

Physical and morphological characteristics of the first three jets of Pêga jackasses sperm-rich fraction / Características físicas e morfológicas dos três primeiros jatos da fração rica do ejaculado de jumentos da raça Pêga

The first three jets of the sperm-rich fraction of Pêga jackasses were collected and assessed separately. Five fertile Pêga jackasses were used as semen donors and underwent fractionated semen collection, using an open model artificial vagina. The first three jets of the semen were collected separately and assessed for volume, sperm motility, vigor, concentration/mL of semen, and sperm morphology. These characteristics were compared between first, second and third jets and between jackasses. It was observed that the jet volume differed (P<0.05) between jackasses, although it was similar (P>0.05) between first, second and third jets. Sperm motility did not differ (P>0.05) between jets and jackasses. Vigor was similar (P>0.05) between jets of the same jackass, and only the first jet differed (P<0.05) between jackasses. The first, second and third jets of the sperm-rich fraction had decreased sperm concentrations (P<0.05) of 955.56, 725.56 and 280.56x 106 sperm/mL of semen, respectively. Sperm morphology differed between the first three jets only for the incidence of mid-piece defect, higher in the third one (4.26%), compared to the first (3.36%) and second (3.38%) ones. When comparing the morphological characteristics of the sperm-rich fraction between five jackasses, regardless of the jet, there were differences in the percentage of normal sperm, proximal cytoplasmic droplet, mid-piece and head defects.(AU)

Objetivou-se, no presente experimento, caracterizar os três primeiros jatos da fração rica do ejaculado de jumentos da raça Pêga. Para tal, cinco jumentos foram submetidos à coleta fracionada do sêmen, utilizando-se vagina artificial modelo aberta. Os três primeiros jatos da fração rica do ejaculado foram coletados separadamente e avaliados quanto ao volume, à motilidade e ao vigor espermáticos, à concentração espermática/mL de sêmen e à morfologia espermática. Comparações foram realizadas entre jatos e entre jumentos. Observou-se que o volume do jato diferiu (P<0,05) entre os jumentos, embora fosse similar (P>0,05) entre os jatos. A motilidade espermática não diferiu (P>0,05) entre jatos nem entre jumentos. O vigor espermático foi similar (P>0,05) entre os jatos de um mesmo jumento, e apenas o vigor do jato 1 diferiu (P<0,05) entre os jumentos. Independentemente do jumento, a fração rica do ejaculado foi composta por três jatos apresentando concentrações espermáticas decrescentes (P<0,05), com 955,56; 725,56 e 280,56 x 106 espermatozoides/mL de sêmen. A morfologia espermática diferiu entre os três jatos avaliados apenas para a incidência de defeitos de peça intermediária, sendo maior no jato três (4,26%), em relação aos jatos um (3,36%) e dois (3,38%). Comparando-se as características morfológicas do sêmen entre os cinco jumentos avaliados, independentemente do jato, observaram-se diferenças entre os reprodutores quanto ao percentual de espermatozoides normais, com gota citoplasmática proximal, com defeitos de peça intermediária e de cabeça.(AU)

Genetic parameter estimates for carcass traits and visual scores including or not genomic information.

The objective of this study was to determine whether visual scores used as selection criteria in Nellore breeding programs are effective indicators of carcass traits measured after slaughter. Additionally, this study evaluated the effect of different structures of the relationship matrix ( and ) on the estimation of genetic parameters and on the prediction accuracy of breeding values. There were 13,524 animals for visual scores of conformation (CS), finishing precocity (FP), and muscling (MS) and 1,753, 1,747, and 1,564 for LM area (LMA), backfat thickness (BF), and HCW, respectively. Of these, 1,566 animals were genotyped using a high-density panel containing 777,962 SNP. Six analyses were performed using multitrait animal models, each including the 3 visual scores and 1 carcass trait. For the visual scores, the model included direct additive genetic and residual random effects and the fixed effects of contemporary group (defined by year of birth, management group at yearling, and farm) and the linear effect of age of animal at yearling. The same model was used for the carcass traits, replacing the effect of age of animal at yearling with the linear effect of age of animal at slaughter. The variance and covariance components were estimated by the REML method in analyses using the numerator relationship matrix () or combining the genomic and the numerator relationship matrices (). The heritability estimates for the visual scores obtained with the 2 methods were similar and of moderate magnitude (0.23-0.34), indicating that these traits should response to direct selection. The heritabilities for LMA, BF, and HCW were 0.13, 0.07, and 0.17, respectively, using matrix and 0.29, 0.16, and 0.23, respectively, using matrix . The genetic correlations between the visual scores and carcass traits were positive, and higher correlations were generally obtained when matrix was used. Considering the difficulties and cost of measuring carcass traits postmortem, visual scores of CS, FP, and MS could be used as selection criteria to improve HCW, BF, and LMA. The use of genomic information permitted the detection of greater additive genetic variability for LMA and BF. For HCW, the high magnitude of the genetic correlations with visual scores was probably sufficient to recover genetic variability. The methods provided similar breeding value accuracies, especially for the visual scores.

Babesia bovis and Babesia bigemina infection levels estimated by qPCR in Angus cattle from an endemic area of São Paulo state, Brazil.

The levels of infection by Babesia bovis and Babesia bigemina were estimated by absolute quantification through the quantitative PCR technique (qPCR). Fifty-one contemporaneous Angus cattle were evaluated on two occasions. The number of standard female Rhipicephalus microplus ticks present on the left side of the body was counted and blood samples were drawn from the tail vein into tubes containing the anticoagulant EDTA. The blood samples were submitted to DNA extraction and used to quantify the number of copies (NC) of DNA from B. bovis and B. bigemina by qPCR. The data on tick count and number of DNA copies were transformed for normalization and analyzed by a mixed model method. A multivariate model with repeated measures of the same animal, including the effects of collection, parasite species and their interaction, was used. The repeatability values were obtained from the matrix of (co)variances and were expressed for each species. The correlations between the counts of different species on the same animal, in the same collection or different collections, were also estimated. The results showed the qPCR could distinguish the two between infection by the two Babesia species. Infection levels by B. bovis and B. bigemina were detected in 100% and 98% of the animals, respectively. Significant differences were found (P<0.05) between the NC of the two Babesia species, B. bovis 1.49±0.07 vs. B. bigemina 0.82±0.06. Low repeatabilities were found for the counts of R. microplus and NC of B. bovis and B. bigemina: 0.05, 0.10 and 0.02, respectively. The correlations between R. microplus count and NC of B. bovis and B. bigemina were both very near zero. However, an association was observed between the NC of the two species, with a correlation coefficient of 0.30 for measures from the same collection. The absence of associations between the quantity of DNA from B. bovis and B. bigemina and the tick counts suggests that the variation of parasitemia by the hemoparasites did not depend on the tick infestation levels at the moment of each collection. The repeatability values estimated indicate that under the study conditions, the variations in the tick infestation levels and of parasitemia by B. bovis and B. bigemina depend more on factors related to each collection than on intrinsic factors of the animal.

Genetic association of productive and reproductive traits with stayability in Nellore cattle: analysis using Bayesian models.

We evaluated the genetic association of growth traits [weight adjusted to 205 days of age (W205), 365 days of age (W365), and 550 days of age (W550); weight gain between 205 days of age and 365 days of age (WG1) and between 365 days of age and 550 days of age (WG2)] and reproductive traits [age at first calving (AFC); first calving interval (FCI)] with stayability in the herd (STAY), using Bayesian inference in linear and threshold models. We defined STAY as the probability of a cow calving three or more times before the age of 76 months, given that she had calved at least once. We assigned binary codes (0, failure; 1, success) to each female. We used a sire model for analysis and formed different contemporary groups for the investigated traits. We analyzed the results by applying a two-trait sire model that included STAY (threshold trait) and linear traits (W205, W365, W550, WG1, WG2, AFC, and FCI). We used Gibbs sampling to estimate variance components and heritabilities. In all the analyses, we found that the mean heritability estimates for STAY were of moderate magnitude (0.20-0.25). The mean heritabilities for W205, W365, W550, WG1, WG2, AFC, and FCI were 0.20, 0.23, 0.39, 0.08, 0.14, 0.12, and 0.11, respectively. We observed wide variation in the posterior distributions of genetic correlations; however, with the exception of those obtained for the reproductive traits, the mean estimates were of low magnitude. Selection for WG2 can results in favorable correlated response in STAY.

Associations between single nucleotide polymorphisms and carcass traits in Nellore cattle using high-density panels.

The objective of this study was to evaluate associations between single nucleotide polymorphism (SNP) markers and carcass traits measured postmortem in Nellore cattle. Records of loin eye area (LEA) and backfat thickness (BF) from 740 males and records of hot carcass weight (HCW) from 726 males were analyzed. All of the animals were genotyped using the BovineHD BeadChip. Association analyses were performed by the restricted maximum likelihood method that considered one SNP at a time. Significant SNPs were identified on chromosomes 2 and 6 for LEA and on chromosomes 7, 1, and 2 for BF. For HCW, associations with SNPs were found on chromosomes 13, 14, and 28, in addition to genome regions that were directly related to this trait, such as the EFCAB8 and VSTM2L genes, and to bone development (RHOU). Some SNPs were located in very close proximity to genes involved in basal metabolism (BLCAP, NNAT, CTNNBL1, TGM2, and LOC100296770) and the immune system (BPI).

Reaction norms for the study of genotype-environment interaction for growth and indicator traits of sexual precocity in Nellore cattle.

The objective of this study was to quantify the magnitude of genotype-environment interaction (GxE) effects on age at first calving (AFC), scrotal circumference (SC), and yearling weight (YW) in Nellore cattle using reaction norms. For the study, 89,152 weight records of female and male Nellore animals obtained at yearling age were used. Genetic parameters were estimated with a single-trait random-regression model using Legendre polynomials as base functions. The heritability estimates were of low to medium magnitude for AFC (0.05 to 0.47) and of medium to high magnitude for SC (0.32 to 0.51) and YW (0.13 to 0.72), and increased as the environmental gradient became more favorable. The genetic correlation estimates ranged from 0.25 to 1.0 for AFC, from 0.71 to 1.0 for SC, and from 0.42 to 1.0 for YW. High Spearman correlation coefficients were obtained for the three traits, ranging from 0.97 to 0.99. The reaction norms along the environmental gradient of 10 sires each with the highest or lowest breeding value for YW predicted by single-trait analysis demonstrated more plastic phenotypes for YW and more robust phenotypes for SC. The effect of GxE was most important for YW and AFC with respect to SC. When animals are selected for higher SC or YW or lower AFC, considering or not the GxE effect, it is expected that the same animals will be selected. The reaction norms obtained based on sire breeding values along the environmental gradient showed that animals with extreme breeding values respond differently as environmental conditions improve.

Genotype × environment interaction for age at first calving, scrotal circumference, and yearling weight in Nellore cattle using reaction norms in multitrait random regression models.

The aim of this study was to evaluate the effect of genotype × environment interaction (G×E) on age at first calving (AFC), scrotal circumference (SC), and yearling weight (YW) and to estimate genetic correlations between these traits in Nellore cattle using reaction norms in multitrait random regression models. In this study, 28,871, 41,386, and 89,152 records of Nellore cattle for AFC, SC, and YW, respectively, were used. The data were obtained from farms located in the north, northeast, midwest, and southeast regions of Brazil that participate in the DeltaGen Breeding Program. Environmental levels were defined as a function of contemporary groups, that is, animals born in the same herd and year, from the same management group (from birth to yearling), and of the same sex. Postweaning weight gain was used as a criterion to evaluate the environmental conditions for all traits. For reaction norm analyses, residual variances were modeled with homogeneous and heterogeneous classes. The model for SC and YW included the fixed effects of contemporary group and age of the animal as a covariate as well as random direct additive genetic and residual effects. The same model, excluding the covariate age of the animal, was used for AFC. The heritability estimates were low to high for AFC (0.09 to 0.50), high for SC (0.51 to 0.67), and moderate to high for YW (0.33 to 0.71). The genetic correlations (within each trait) along the environmental levels varied from -0.27 to 1.0 for AFC, from 0.73 to 1.0 for SC, and from 0.26 to 1.0 for YW. The genetic correlations between different traits in different environments varied from -0.14 to -0.60 between AFC and SC, from -0.05 to -0.32 between AFC and YW, and from -0.05 to 0.72 between YW and SC. The genetic correlations have had different magnitudes for AFC, SC, and YW, which could indicate the presence of G×E. The present results should support researchers and farmers in defining selection criteria to improve growth traits and sexual precocity. Our results suggest that animals for breeding have to be selected in the same environment and management conditions as their progeny will be reared.

Study of lipid metabolism-related genes as candidate genes of sexual precocity in Nellore cattle.

The objective of this study was to identify associations between known polymorphisms in genes related to adipose tissue and sexual precocity in Nellore cattle. A total of 1689 precocious and non-precocious heifers belonging to farms participating in Conexão Delta G breeding program were studied. SNPs from the Illumina High-Density Bovine SNP BeadChip were used. This chip contains 777,000 SNPs located within the region of the candidate genes at a distance of up to 5 kb, considering that linkage disequilibrium (LD) exists at this distance. Linear models were used for statistical analysis. The fastPHASE and GenomeStudio programs were used for haplotype reconstruction and LD analysis based on r2 statistics. Fifty-seven candidate genes and 443 SNPs were analyzed: among the latter, 370 SNPs formed 83 haplotypes, while the remaining SNPs were studied separately. Statistical analysis showed that only three haplotypes, one haplotype consisting of two SNPs located in the FABP4 gene and two haplotypes consisting of four and two SNPs located in the PPP3CA gene, had a significant effect on sexual precocity at P < 0.05. It can be concluded that the FABP4 and PPP3CA genes influence sexual precocity and may therefore be used in selection programs designed to improve sexual precocity in Nellore cattle.

Polymorphisms in the GHRL gene and their associations with traits of economic interest in beef cattle.

The hormone ghrelin is produced in the stomach wall, has an orexigenic function, stimulates growth hormone secretion, and affects the energy balance of the animal. Therefore, the ghrelin gene (GHRL) is considered to be a good candidate marker for the identification of traits of great economic importance in cattle, such as those associated with feed intake, growth, and carcass quality. The use of molecular genetic markers associated with such traits permits the earlier and more accurate identification of superior animals, thus reducing the interval between generations, and increasing the genetic gain. Six SNPs were found in the GHRL gene, located in intron 3, intron 4, and exon 5. The positions of the SNPs on the gene and the substitutions were: g.2184A>G, g.2347T>C, g.4469T>C, g.4548A>G, g.4663T>C, and g.4729T>C (GenBank accession No. JX565585). After analysis of linkage disequilibrium, association tests were performed between four SNPs with the traits year weight for males, yearling weight for females, dry matter intake, loin eye area, and rump fat thickness (P ≤ 0.05). Therefore, GHRL is an important candidate gene that may be used to identify genetic variations that influence traits of economic importance in beef cattle.

Genetic association of growth traits with carcass and meat traits in Nellore cattle.

The aim of this study was to estimate genetic and phenotypic associations of growth traits with carcass and meat traits in Nellore cattle. Data from male and female animals were used for weaning weight (WW; N = 241,416), yearling weight (YW, N = 126,596), weight gain from weaning to yearling (GWY, N = 78,687), and yearling hip height (YHH, N = 90,720), respectively; 877 male animals were used for hot carcass weight (HCW) and 884 for longissimus muscle area (LMA), backfat thickness (BT), marbling score (MS), and shear force (SF). The variance components were estimated by the restricted maximum likelihood method using three-trait animal models that included WW. The model for WW included direct and maternal additive genetic, maternal permanent environmental, and residual effects as random effects; contemporary group as fixed effects; and age of dam at calving and age of animal as covariates (linear and quadratic effects). For the other traits, maternal effects and the effect of age of dam at calving were excluded from the model. Heritability ranged from 0.10 ± 0.12 (LMA) to 0.44 ± 0.007 (YW). Genetic correlations ranged from -0.40 ± 0.38 (WW x LMA) to 0.55 ± 0.10 (HCW x YW). Growth, carcass, and meat traits have sufficient genetic variability to be included as selection criteria in animal breeding programs.

Prospection of genomic regions divergently selected in racing line of Quarter Horses in relation to cutting line.

Selection of Quarter Horses for different purposes has led to the formation of lines, including racing and cutting horses. The objective of this study was to identify genomic regions divergently selected in racing line of Quarter Horses in relation to cutting line applying relative extended haplotype homozygosity (REHH) analysis, an extension of extended haplotype homozygosity (EHH) analysis, and the fixation index (F ST) statistic. A total of 188 horses of both sexes, born between 1985 and 2009 and registered at the Brazilian Association of Quarter Horse Breeders, including 120 of the racing line and 68 of the cutting line, were genotyped using single nucleotide polymorphism arrays. On the basis of 27 genomic regions identified as selection signatures by REHH and F ST statistics, functional annotations of genes were made in order to identify those that could have been important during formation of the racing line and that could be used subsequently for the development of selection tools. Genes involved in muscle growth (n=8), skeletal growth (n=10), muscle energy metabolism (n=15), cardiovascular system (n=14) and nervous system (n=23) were identified, including the FKTN, INSR, GYS1, CLCN1, MYLK, SYK, ANG, CNTFR and HTR2B.

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.