Estimating genetic parameters of reproductive, carcass, and meat quality traits in Polled Nellore cattle.

Considering the economic and commercial efficiency of the beef production chain, the yield and quality of the meat produced must also be included in breeding programs. For the Nellore breed, including the polled herd, these aspects have not been much studied. The aim of this study was to estimate genetic parameters for scrotal circumference adjusted to 365 (SC365) and 450 (SC450) days of age, age at first calving (AFC), accumulated productivity (AP), stayability (STAY), longissimus muscle area (LMA), thickness of subcutaneous fat over the 12th-13th ribs (BF), thickness of subcutaneous fat over the rump (RF), and shear force measured by Warner-Bratzler shear force (WBSF) of polled Nellore cattle. Bayesian analyses were performed by adopting a linear animal model, whereas STAY analyses used the linear threshold model. Heritability estimates were 0.31 (SC365), 0.37 (SC450), 0.16 (AFC), 0.25 (AP), 0.16 (STAY), 0.30 (LMA), 0.13 (BF), 0.24 (RF), and 0.15 (WBSF), indicating moderate response to selection. Genetic and residual correlations between SC365 and SC450 were high (0.91 and 0.74, respectively), as well as the genetic correlations of AP with SC365, SC450, AFC, and STAY (0.61, 0.62, - 0.69, and 0.83, respectively). Genetic and residual correlations of WBSF with reproductive and carcass characteristics exhibited high standard deviations, however favorable. Based on the results, it is expected that in the medium term, animals with greater sexual precocity will also have greater accumulated productivity and longer permanence of females in the herd, along with superior carcass traits. However, due to the low heritabilities and small genetic associations with reproductive traits, fat thickness characteristics (BF and RF) will still require direct selection.

Norovirus contamination in different shellfish species harvested in the same production areas.

AIMS: To investigate Norovirus (NoV) contamination of mussels, clams and oysters harvested in two class B harvesting areas of the delta of the Po river, to choose a species as an indicator. METHODS AND RESULTS: Environmental parameters (temperature and salinity) and hydrometric levels of the tributary river were measured. Seventy shellfish samples (35 samples per area) were examined for Escherichia coli and NoV (GI and GII). NoV contamination was found in 51.4% of samples, of which, 2.9% contained only NoV GI, 14.3% only NoV GII, while the majority of the samples (34.3%) contained both genogroups. Most of the positive results (90.0%) were obtained in the period between November 2008 and April 2009. CONCLUSIONS: No significant differences were found between the results from the two harvesting areas and the three shellfish species. However, on the basis of the average C(t) values, the recovery rate (from 0.46 to 1.15%) and the distribution of positive results in the samplings, mussels seem to be a suitable indicator species to monitor viral contamination in these areas. SIGNIFICANCE AND IMPACT OF THE STUDY: The data allow the optimization of monitoring plans to improve the prevention strategies in terms of money and time, by the intensification of controls in the cold season and the use of one species as indicator.

Effect of high hydrostatic pressure on murine norovirus in Manila clams.

AIMS: Eating raw or insufficiently cooked bivalve molluscs contaminated with human noroviruses (NVs) can result in acute cases of gastroenteritis in humans. Manila clams (Ruditapes philippinarum) are particularly prone to exposure to NVs due to the brackish environment in which they are farmed which is known to be susceptible to human faecal contamination. High hydrostatic pressure processing (HHP) is a food treatment technique that has been shown to inactivate NV. METHODS AND RESULTS: In this study we investigated the ability of HHP to inactivate murine norovirus (MNV-1), a recognised surrogate for NV, in experimentally contaminated manila clams. Pools of contaminated live clams were subjected to hydrostatic pressure ranging from 300 to 500 MPa for different time intervals of between one and 10 min. The trial was repeated three times, at monthly intervals. CONCLUSIONS: Virus vitality post-treatment was assessed and the data obtained indicates that the use of high hydrostatic pressures of at least 500 MPa for 1 min was effective in inactivating MNV-1. SIGNIFICANCE AND IMPACT OF THE STUDY: HHP results to be an effective technique that could be applied to industrial process to obtain safe Manila clams ready to eat.

Genome-enabled prediction of meat and carcass traits using Bayesian regression, single-step genomic best linear unbiased prediction and blending methods in Nelore cattle.

Several methods have been used for genome-enabled prediction (or genomic selection) of complex traits, for example, multiple regression models describing a target trait with a linear function of a set of genetic markers. Genomic selection studies have been focused mostly on single-trait analyses. However, most profitability traits are genetically correlated, and an increase in prediction accuracy of genomic breeding values for genetically correlated traits is expected when using multiple-trait models. Thus, this study was carried out to assess the accuracy of genomic prediction for carcass and meat quality traits in Nelore cattle, using single- and multiple-trait approaches. The study considered 15 780, 15 784, 15 742 and 526 records of rib eye area (REA, cm2), back fat thickness (BF, mm), rump fat (RF, mm) and Warner-Bratzler shear force (WBSF, kg), respectively, in Nelore cattle, from the Nelore Brazil Breeding Program. Animals were genotyped with a low-density single nucleotide polymorphism (SNP) panel and subsequently imputed to arrays with 54 and 777 k SNPs. Four Bayesian specifications of genomic regression models, namely, Bayes A, Bayes B, Bayes Cπ and Bayesian Ridge Regression; blending methods, BLUP; and single-step genomic best linear unbiased prediction (ssGBLUP) methods were compared in terms of prediction accuracy using a fivefold cross-validation. Estimates of heritability ranged from 0.20 to 0.35 and from 0.21 to 0.46 for RF and WBSF on single- and multiple-trait analyses, respectively. Prediction accuracies for REA, BF, RF and WBSF were all similar using the different specifications of regression models. In addition, this study has shown the impact of genomic information upon genetic evaluations in beef cattle using the multiple-trait model, which was also advantageous compared to the single-trait model because it accounted for the selection process using multiple traits at the same time. The advantage of multi-trait analyses is attributed to the consideration of correlations and genetic influences between the traits, in addition to the non-random association of alleles.

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

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

Genetic correlations between categorical morphological traits in Nelore cattle by applying Bayesian analysis under a threshold animal model.

In this study, Bayesian analysis under a threshold animal model was used to estimate genetic correlations between morphological traits (body structure, finishing precocity and muscling) in Nelore cattle evaluated at weaning and yearling. Visual scores obtained from 7651 Nelore cattle at weaning and from 4155 animals at yearling, belonging to the Brazilian Nelore Program, were used. Genetic parameters for the morphological traits were estimated by two-trait Bayesian analysis under a threshold animal model. The genetic correlations between the morphological traits evaluated at two ages of the animal (weaning and yearling) were positive and high for body structure (0.91), finishing precocity (0.96) and muscling (0.94). These results indicate that the traits are mainly determined by the same set of genes of additive action and that direct selection at weaning will also result in genetic progress for the same traits at yearling. Thus, selection of the best genotypes during only one phase of life of the animal is suggested. However, genetic differences between morphological traits were better detected during the growth phase to yearling. Direct selection for body structure, finishing precocity and muscling at only one age, preferentially at yearling, is recommended as genetic differences between traits can be detected at this age.

New ecosystems in the deep subsurface follow the flow of water driven by geological activity.

Eukarya have been discovered in the deep subsurface at several locations in South Africa, but how organisms reach the subsurface remains unknown. We studied river-subsurface fissure water systems and identified Eukarya from a river that are genetically identical for 18S rDNA. To further confirm that these are identical species one metazoan species recovered from the overlying river interbred successfully with specimen recovered from an underlying mine at -1.4 km. In situ seismic simulation experiments were carried out and show seismic activity to be a major force increasing the hydraulic conductivity in faults allowing organisms to create ecosystems in the deep subsurface. As seismic activity is a non-selective force we recovered specimen of algae and Insecta that defy any obvious other explanation at a depth of -3.4 km. Our results show there is a steady flow of surface organisms to the deep subsurface where some survive and adapt and others perish. As seismic activity is also present on other planets and moons in our solar system the mechanism elucidated here may be relevant for future search and selection of landing sites in planetary exploration.

Dating phototrophic microbial lineages with reticulate gene histories.

Phototrophic bacteria are among the most biogeochemically significant organisms on Earth and are physiologically related through the use of reaction centers to collect photons for energy metabolism. However, the major phototrophic lineages are not closely related to one another in bacterial phylogeny, and the origins of their respective photosynthetic machinery remain obscured by time and low sequence similarity. To better understand the co-evolution of Cyanobacteria and other ancient anoxygenic phototrophic lineages with respect to geologic time, we designed and implemented a variety of molecular clocks that use horizontal gene transfer (HGT) as additional, relative constraints. These HGT constraints improve the precision of phototroph divergence date estimates and indicate that stem green non-sulfur bacteria are likely the oldest phototrophic lineage. Concurrently, crown Cyanobacteria age estimates ranged from 2.2 Ga to 2.7 Ga, with stem Cyanobacteria diverging ~2.8 Ga. These estimates provide a several hundred Ma window for oxygenic photosynthesis to evolve prior to the Great Oxidation Event (GOE) ~2.3 Ga. In all models, crown green sulfur bacteria diversify after the loss of the banded iron formations from the sedimentary record (~1.8 Ga) and may indicate the expansion of the lineage into a new ecological niche following the GOE. Our date estimates also provide a timeline to investigate the temporal feasibility of different photosystem HGT events between phototrophic lineages. Using this approach, we infer that stem Cyanobacteria are unlikely to be the recipient of an HGT of photosystem I proteins from green sulfur bacteria but could still have been either the HGT donor or the recipient of photosystem II proteins with green non-sulfur bacteria, prior to the GOE. Together, these results indicate that HGT-constrained molecular clocks are useful tools for the evaluation of various geological and evolutionary hypotheses, using the evolutionary histories of both genes and organismal lineages.

Fluctuations in populations of subsurface methane oxidizers in coordination with changes in electron acceptor availability.

The concentrations of electron donors and acceptors in the terrestrial subsurface biosphere fluctuate due to migration and mixing of subsurface fluids, but the mechanisms and rates at which microbial communities respond to these changes are largely unknown. Subsurface microbial communities exhibit long cellular turnover times and are often considered relatively static-generating just enough ATP for cellular maintenance. Here, we investigated how subsurface populations of CH4 oxidizers respond to changes in electron acceptor availability by monitoring the biological and geochemical composition in a 1339 m-below-land-surface (mbls) fluid-filled fracture over the course of both longer (2.5 year) and shorter (2-week) time scales. Using a combination of metagenomic, metatranscriptomic, and metaproteomic analyses, we observe that the CH4 oxidizers within the subsurface microbial community change in coordination with electron acceptor availability over time. We then validate these findings through a series of 13C-CH4 laboratory incubation experiments, highlighting a connection between composition of subsurface CH4 oxidizing communities and electron acceptor availability.

Genomewide association mapping and pathway analysis of meat tenderness in Polled Nellore cattle.

Brazil is one of the world's largest beef exporters, although the product has a low price due to quality issues. The meat exported by Brazil is considered medium and low quality by international buyers, mainly due to lack of tenderness. The predominant Zebu breeds (80% Nellore) are known for producing tougher beef than taurine breeds. Nonetheless, some studies have shown that there is substantial genetic variability for tenderness within the Nellore breed, although it is a difficult trait to improve by conventional selection methods. Therefore, the aim of this study was to perform a genomewide association study (GWAS) and a gene set enrichment analysis to identify genomic regions and biologically relevant pathways associated with meat tenderness in Polled Nellore cattle. Data consisted of Warner-Bratzler shear force values of LM from 427 Polled Nellore animals divided into 3 experimental slaughters (years 2005, 2008, and 2010). The animals were genotyped with either the Illumina BovineHD BeadChip (777k, on 61 samples) or the GGP Indicus HD chip (77k, on 366 samples). Single nucleotide polymorphisms were excluded when the call rate was <90%, the Hardy-Weinberg proportions -value was <1% (Fisher exact test, Bonferroni adjusted), and the minor allele frequency was <1%. Imputation from the GGP Indicus HD chip to the Illumina BovineHD BeadChip was performed using the FImput program. Genomewide association analysis was performed using the Efficient Mixed Model Association eXpedited (EMMAx) and the population parameters previously determined (P3D) methods. The GWAS was complemented with a gene set enrichment analysis performed using the FatiGO procedure. Significant markers ( < 0.0001) explaining a larger proportion of variation than other significant SNPs were located on chromosomes 3, 13, 17, 20, 21, and 25, indicating QTL associated with meat tenderness throughout the genome. Additionally, gene set analysis identified 22 Gene Ontology functional terms and 2 InterPro entries that showed significant enrichment of genes associated with tenderness. The functional categories included protein tyrosine and serine/threonine kinase activity, calcium ion binding, lipid metabolic process, and growth factors, among others. These results help to elucidate the genetic architecture and metabolic pathways underlying this trait, which is of extreme economic and social importance to Brazil, because Nellore is the dominant beef cattle breed in the country.

Accuracy of genomic breeding values for meat tenderness in Polled Nellore cattle.

Zebu () cattle, mostly of the Nellore breed, comprise more than 80% of the beef cattle in Brazil, given their tolerance of the tropical climate and high resistance to ectoparasites. Despite their advantages for production in tropical environments, zebu cattle tend to produce tougher meat than Bos taurus breeds. Traditional genetic selection to improve meat tenderness is constrained by the difficulty and cost of phenotypic evaluation for meat quality. Therefore, genomic selection may be the best strategy to improve meat quality traits. This study was performed to compare the accuracies of different Bayesian regression models in predicting molecular breeding values for meat tenderness in Polled Nellore cattle. The data set was composed of Warner-Bratzler shear force (WBSF) of longissimus muscle from 205, 141, and 81 animals slaughtered in 2005, 2010, and 2012, respectively, which were selected and mated so as to create extreme segregation for WBSF. The animals were genotyped with either the Illumina BovineHD (HD; 777,000 from 90 samples) chip or the GeneSeek Genomic Profiler (GGP Indicus HD; 77,000 from 337 samples). The quality controls of SNP were Hard-Weinberg Proportion -value ≥ 0.1%, minor allele frequency > 1%, and call rate > 90%. The FImpute program was used for imputation from the GGP Indicus HD chip to the HD chip. The effect of each SNP was estimated using ridge regression, least absolute shrinkage and selection operator (LASSO), Bayes A, Bayes B, and Bayes Cπ methods. Different numbers of SNP were used, with 1, 2, 3, 4, 5, 7, 10, 20, 40, 60, 80, or 100% of the markers preselected based on their significance test (-value from genomewide association studies [GWAS]) or randomly sampled. The prediction accuracy was assessed by the correlation between genomic breeding value and the observed WBSF phenotype, using a leave-one-out cross-validation methodology. The prediction accuracies using all markers were all very similar for all models, ranging from 0.22 (Bayes Cπ) to 0.25 (Bayes B). When preselecting SNP based on GWAS results, the highest correlation (0.27) between WBSF and the genomic breeding value was achieved using the Bayesian LASSO model with 15,030 (3%) markers. Although this study used relatively few animals, the design of the segregating population ensured wide genetic variability for meat tenderness, which was important to achieve acceptable accuracy of genomic prediction. Although all models showed similar levels of prediction accuracy, some small advantages were observed with the Bayes B approach when higher numbers of markers were preselected based on their -values resulting from a GWAS analysis.

Survey of Anisakis sp. and Hysterothylacium sp. in sardines and anchovies from the North Adriatic Sea.

The occurrence of larval Anisakidae and Raphidascarididae in anchovies and sardines from the North Adriatic Sea has been estimated. Anisakis pegreffii and Hysterothylacium aduncum were reported, with low prevalence values. In brief, a total amount of 7650 fish specimens collected between September 2011 and 2012 were analysed using three different inspection analyses: a visual inspection of the coelomic cavity, an examination of the viscera exploiting the positive hydro-tropism of the larvae (modified Baermann technique) and enzymatic digestion of muscular tissue pools. Low level of infestation was reported for Anisakis sp. in both in anchovies and sardines, while higher values were reported for Hysterothylacium sp. Subsamples of nematodes collected were characterized at species level using the molecular diagnostic key based on ITS nuclear ribosomal region, and A. pegreffii and H. aduncum were identified. The low prevalence of Anisakis sp. in sardines and anchovies from the North Adriatic Sea could be related to the peculiar distribution of cetaceans and carnivorous zooplankton in the investigated region and could be used as a potential tag to define oily fishes from this specific fishing area as at low-risk for anisakiasis.

Prevalence of Toxoplasma gondii antibodies in sera of domestic cats from Guarulhos and São Paulo, Brazil.

Antibodies to Toxoplasma gondii were determined in serum samples of 502 domestic cats from Brazil by the modified agglutination test (MAT), using formalin-fixed whole tachyzoites and mercaptoethanol. Antibodies (MAT > or = 1:20) were found in 132 (26.3%) of 502 cats. With respect to origin, antibodies were found in 26.7% of 430 stray cats from São Paulo, 10% of 40 stray cats from Guarulhos, and 40.6% of 32 cats from a cat breeder in São Paulo. Antibody titers were: 1:20 in 10 cats, 1:25 in 40 cats, 1:50 in 73 cats, and > or = 1:500 in 9 cats. Exposure rates of T. gondii in cats from São Paulo, Brazil are similar to that in domestic cats in North America.

Does aspartic acid racemization constrain the depth limit of the subsurface biosphere?

Previous studies of the subsurface biosphere have deduced average cellular doubling times of hundreds to thousands of years based upon geochemical models. We have directly constrained the in situ average cellular protein turnover or doubling times for metabolically active micro-organisms based on cellular amino acid abundances, D/L values of cellular aspartic acid, and the in vivo aspartic acid racemization rate. Application of this method to planktonic microbial communities collected from deep fractures in South Africa yielded maximum cellular amino acid turnover times of ~89 years for 1 km depth and 27 °C and 1-2 years for 3 km depth and 54 °C. The latter turnover times are much shorter than previously estimated cellular turnover times based upon geochemical arguments. The aspartic acid racemization rate at higher temperatures yields cellular protein doubling times that are consistent with the survival times of hyperthermophilic strains and predicts that at temperatures of 85 °C, cells must replace proteins every couple of days to maintain enzymatic activity. Such a high maintenance requirement may be the principal limit on the abundance of living micro-organisms in the deep, hot subsurface biosphere, as well as a potential limit on their activity. The measurement of the D/L of aspartic acid in biological samples is a potentially powerful tool for deep, fractured continental and oceanic crustal settings where geochemical models of carbon turnover times are poorly constrained. Experimental observations on the racemization rates of aspartic acid in living thermophiles and hyperthermophiles could test this hypothesis. The development of corrections for cell wall peptides and spores will be required, however, to improve the accuracy of these estimates for environmental samples.

Características reprodutivas e suas associações com outras características de importância econômica na raça Nelore / Reproductive traits and their associations with other economically important traits in Nelore

Objetivou-se com este trabalho estimar as herdabilidades (h²) e as correlações genéticas (r g) entre idade ao primeiro parto (IPP) e primeiro intervalo de partos (PIEP) e outras características como peso (PS) ao ano (A) e ao sobreano (S), altura do posterior (ALT) e perímetro escrotal (PE450) em animais da raça Nelore. Os parâmetros genéticos foram estimados em uma análise multicaracterística por modelo animal, utilizando-se a inferência bayesiana via algoritmo de "Gibbs Sampling". Os parâmetros genéticos estimados sugerem a existência de variabilidade genética para IPP (h² = 0,26), sendo que a seleção para a diminuição da IPP de fêmeas Nelore deve responder à seleção individual, sem causar antagonismo do valor genético dos animais para PS (r g = -0,22 (A) e -0,44 (S)) e PE450 (r g = 0,02). A seleção para a diminuição da IPP, no longo prazo, pode levar a um aumento da ALT dos animais, embora essa associação seja relativamente baixa (-0,35). A estimativa de herdabilidade a posteriori para a característica PIEP foi baixa, 0,11±0,03. As r g entre PIEP e as demais características estudadas indicam que a seleção para essas características de crescimento não afetará o PIEP.

Heritability (h²) and genetic correlations (r g) were estimated between reproductive traits such as age at first calving (AFC), first calving interval (FCI) and other economically relevant traits, i.e., weight (W) at year (Y) and at 18 months of age (S), scrotal circumference (SC), and hip height (HH) in Nelore cattle. The genetic parameters were estimated in a multiple-trait analysis, with animal models using the Bayesian inference by Gibbs Sampling algorithm. The genetic parameters estimated in this work suggest the existence of genetic variability for AFC (h² = 0.26), where the selection for the reduction of Nelore females AFC should respond to mass selection, without causing genetic antagonism in the selection of W (r g = -0,22 (Y) and -0,44 (S)), and SC (r g = 0,02). The selection for the AFC in the long term could lead to an increase in the animal's frame, although this association is relatively low (-0.35). The posteriori heritability estimate for FCI was low, 0.11±0.03. The r g between FCI and the other traits studied indicate that selection for these growth traits will not affect the FCI.

Análise bayesiana na estimação de correlações genéticas entre escores visuais e características reprodutivas de bovinos Nelore utilizando modelos linear-limiar / Bayesian analysis in the estimation of genetic correlations between visual scores and reproductive traits in Nelore cattle using linear-threshold models

Estimaram-se as correlações genéticas entre os escores visuais e as características reprodutivas, utilizando a estatística bayesiana sob modelo animal linear-limiar, em bovinos da raça Nelore. Foram estudadas características categóricas morfológicas, avaliadas visualmente aos oito, 15 e 22 meses de idade; e características contínuas de perímetro escrotal padronizado aos 365 e 450 dias de idade, além da idade ao primeiro parto. As estimativas de correlações genéticas foram de sentido favorável à seleção, apresentando magnitudes moderadas, sugerindo que a seleção de animais para um biótipo desejável pode levar a animais com maior fertilidade e precocidade sexual. As estimativas de correlação genética para o perímetro escrotal padronizado aos 450 dias e a idade ao primeiro parto com as características morfológicas avaliadas aos 22 meses de idade foram maiores do que as obtidas entre as características de escores visuais avaliadas aos oito e 15 meses de idade. A utilização de escores visuais como critério de seleção trará progresso genético também para as características reprodutivas.

The genetic correlations between visual scores and reproductive traits, using the bayesian statistic under the linear-threshold animal model, in bovines of the Nelore breed were estimated. Categorical morphological traits were also evaluated at eight, 15, and 22 month-old, concerning musculature, physical structure, conformation and sacrum, the escrotal circumferences measured at the ages of 365 and 450-day-old, and the age at first calving. The estimates of genetic correlations were moderates, suggesting that the selection of animals with better morphological conformation may lead to animals more precocious and with greater fertility. The estimates of genetic correlation concerning the escrotal circumference at 450 day-old and the age at first calving with the evaluated visual scores at 22 month-old were higher than the ones obtained between visual scores at eight and 15 month-old. The use of these visual scores as a criterion of selection may also bring genetic progress to the reproductive traits.