Inbreeding coefficients and runs of homozygosity islands in Brazilian water buffalo.

Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.

Genomic reaction norm models exploiting genotype × environment interaction on sexual precocity indicator traits in Nellore cattle.

Brazilian beef cattle are raised predominantly on pasture in a wide range of environments. In this scenario, genotype by environment (G×E) interaction is an important source of phenotypic variation in the reproductive traits. Hence, the evaluation of G×E interactions for heifer's early pregnancy (HP) and scrotal circumference (SC) traits in Nellore cattle, belonging to three breeding programs, was carried out to determine the animal's sensitivity to the environmental conditions (EC). The dataset consisted of 85 874 records for HP and 151 553 records for SC, from which 1800 heifers and 3343 young bulls were genotyped with the BovineHD BeadChip. Genotypic information for 826 sires was also used in the analyses. EC levels were based on the contemporary group solutions for yearling body weight. Linear reaction norm models (RNM), using pedigree information (RNM_A) or pedigree and genomic information (RNM_H), were used to infer G×E interactions. Two validation schemes were used to assess the predictive ability, with the following training populations: (a) forward scheme-dataset was split based on year of birth from 2008 for HP and from 2011 for SC; and (b) environment-specific scheme-low EC (-3.0 and -1.5) and high EC (1.5 and 3.0). The inclusion of the H matrix in RNM increased the genetic variance of the intercept and slope by 18.55 and 23.00% on average respectively, and provided genetic parameter estimates that were more accurate than those considering pedigree only. The same trend was observed for heritability estimates, which were 0.28-0.56 for SC and 0.26-0.49 for HP, using RNM_H, and 0.26-0.52 for SC and 0.22-0.45 for HP, using RNM_A. The lowest correlation observed between unfavorable (-3.0) and favorable (3.0) EC levels were 0.30 for HP and -0.12 for SC, indicating the presence of G×E interaction. The G×E interaction effect implied differences in animals' genetic merit and re-ranking of animals on different environmental conditions. SNP marker-environment interaction was detected for Nellore sexual precocity indicator traits with changes in effect and variance across EC levels. The RNM_H captured G×E interaction effects better than RNM_A and improved the predictive ability by around 14.04% for SC and 21.31% for HP. Using the forward scheme increased the overall predictive ability for SC (20.55%) and HP (11.06%) compared with the environment-specific scheme. The results suggest that the inclusion of genomic information combined with the pedigree to assess the G×E interaction leads to more accurate variance components and genetic parameter estimates.

Modelagem das trajetórias médias do peso vivo e do rendimento de carcaça de tilápias do Nilo (Oreochromis niloticus) por meio de funções B-spline / Modeling the mean trajectories of live body weight and carcass yield of Nile tilapia (Oreochromis niloticus) by B-spline functions

Empregando o método dos quadrados mínimos e polinômios B-spline quadráticos, diferentes modelos estatísticos foram testados para identificar o mais apropriado para modelar as trajetórias médias do peso vivo e do rendimento de carcaça de tilápias do Nilo (Oreochromis niloticus). Dados de peso vivo (8.758) e de rendimento de carcaça (2.042) de tilápias com idades entre 106 e 245 dias foram obtidos de 72 famílias provenientes de 36 machos e 72 fêmeas. As variáveis sexo e tanque de criação foram consideradas como classificatórias, e os coeficientes dos polinômios B-spline quadráticos com dois a cinco intervalos de mesmo tamanho foram utilizados como covariáveis. Segundo a maioria dos critérios de ajuste utilizados, os modelos com polinômio B-spline quadrático com cinco intervalos de mesmo tamanho apresentaram os melhores ajustes. O aumento do número de intervalos do polinômio B-spline melhorou o ajuste dos polinômios aos dados. A inclusão dos efeitos classificatórios de sexo, tanque de criação, interação entre esses efeitos e polinômio B-spline quadrático aninhado a essa interação indicou que, com o decorrer do tempo, cada sexo, cultivado em diferente tanque, apresentou trajetória média diferente, sendo necessária a inclusão do aninhamento do tempo na interação sexo x tanque de criação para que, em programas de melhoramento genético da espécie, os valores genéticos dos candidatos à seleção não sejam sub ou superestimados.

Employing the method of least squares and quadratic B-spline polynomials, different statistical models were tested to identify the most appropriate to model the mean trajectories of live weight and carcass yield of Nile tilapia (Oreochromis niloticus). Data of live weight (8,758) and carcass yield (2,042) of tilapias with ages between 106 and 245 days were obtained from 72 families derived from 36 males and 72 females. The sex and tank variables were considered as classificatory and the coefficients of quadratic polynomials B-spline of two to five intervals of the same size were used as covariables. According to most fit criteria used, the models with quadratic B-spline polynomial with five intervals of the same size presented the best adjustments. The increase in number of intervals of B-spline polynomial improved the fit of the polynomials to the data. The inclusion of classificatory effects from sex, tank, the interaction of these effects and the quadratic polynomial B-spline nested in this interaction indicated that, over time, each sex, grown in different tank, presented different mean trajectory, necessitating the inclusion of nesting time in the interaction sex x tank in order to avoid the under or overestimation of breeding values of the selection candidates in breeding programs of this species.

Short communication: Variable number of tandem repeat polymorphisms in DGAT1 gene of buffaloes (Bubalus bubalis) is associated with milk constituents.

The diacylglycerol-O-transferase 1 gene is a positional and functional candidate for milk composition traits. The objective of this study was to evaluate the segregation of the variable number of tandem repeat polymorphisms in the regulatory region of diacylglycerol-O-transferase 1 gene in a water buffalo herd, and to assess the association of this mutation with milk production traits. For this purpose, 196 Murrah buffalo cows were genotyped by PCR. The association of the marker with total milk, fat, and protein yields at 305 d of lactation, milk fat and protein percentage, and somatic cell scores were evaluated by single-trait analyses using a generalized mixed model. Two segregating alleles were identified in the population. The allele with 2 repeats affected fat percentage favorably. The present results suggest that this polymorphism is an interesting marker to include in the genetic evaluation of buffaloes.

Polynomials to model the growth of young bulls in performance tests.

The use of polynomial functions to describe the average growth trajectory and covariance functions of Nellore and MA (21/32 Charolais+11/32 Nellore) young bulls in performance tests was studied. The average growth trajectories and additive genetic and permanent environmental covariance functions were fit with Legendre (linear through quintic) and quadratic B-spline (with two to four intervals) polynomials. In general, the Legendre and quadratic B-spline models that included more covariance parameters provided a better fit with the data. When comparing models with the same number of parameters, the quadratic B-spline provided a better fit than the Legendre polynomials. The quadratic B-spline with four intervals provided the best fit for the Nellore and MA groups. The fitting of random regression models with different types of polynomials (Legendre polynomials or B-spline) affected neither the genetic parameters estimates nor the ranking of the Nellore young bulls. However, fitting different type of polynomials affected the genetic parameters estimates and the ranking of the MA young bulls. Parsimonious Legendre or quadratic B-spline models could be used for genetic evaluation of body weight of Nellore young bulls in performance tests, whereas these parsimonious models were less efficient for animals of the MA genetic group owing to limited data at the extreme ages.