RESUMEN
Feed efficiency plays a major role in the overall profitability and sustainability of the beef cattle industry, as it is directly related to the reduction of the animal demand for input and methane emissions. Traditionally, the average daily feed intake and weight gain are used to calculate feed efficiency traits. However, feed efficiency traits can be analysed longitudinally using random regression models (RRMs), which allow fitting random genetic and environmental effects over time by considering the covariance pattern between the daily records. Therefore, the objectives of this study were to: (1) propose genomic evaluations for dry matter intake (DMI), body weight gain (BWG), residual feed intake (RFI) and residual weight gain (RWG) data collected during an 84-day feedlot test period via RRMs; (2) compare the goodness-of-fit of RRM using Legendre polynomials (LP) and B-spline functions; (3) evaluate the genetic parameters behaviour for feed efficiency traits and their implication for new selection strategies. The datasets were provided by the EMBRAPA-GENEPLUS beef cattle breeding program and included 2920 records for DMI, 2696 records for BWG and 4675 genotyped animals. Genetic parameters and genomic breeding values (GEBVs) were estimated by RRMs under ssGBLUP for Nellore cattle using orthogonal LPs and B-spline. Models were compared based on the deviance information criterion (DIC). The ranking of the average GEBV of each test week and the overall GEBV average were compared by the percentage of individuals in common and the Spearman correlation coefficient (top 1%, 5%, 10% and 100%). The highest goodness-of-fit was obtained with linear B-Spline function considering heterogeneous residual variance. The heritability estimates across the test period for DMI, BWG, RFI and RWG ranged from 0.06 to 0.21, 0.11 to 0.30, 0.03 to 0.26 and 0.07 to 0.27, respectively. DMI and RFI presented within-trait genetic correlations ranging from low to high magnitude across different performance test-day. In contrast, BWG and RWG presented negative genetic correlations between the first 3 weeks and the other days of performance tests. DMI and RFI presented a high-ranking similarity between the GEBV average of week eight and the overall GEBV average, with Spearman correlations and percentages of individuals selected in common ranging from 0.95 to 1.00 and 93 to 100, respectively. Week 11 presented the highest Spearman correlations (ranging from 0.94 to 0.98) and percentages of individuals selected in common (ranging from 85 to 94) of BWG and RWG with the average GEBV of the entire period of the test. In conclusion, the RRM using linear B-splines is a feasible alternative for the genomic evaluation of feed efficiency. Heritability estimates of DMI, RFI, BWG and RWG indicate enough additive genetic variance to achieve a moderate response to selection. A new selection strategy can be adopted by reducing the performance test to 56 days for DMI and RFI selection and 77 days for BWG and RWG selection.
Asunto(s)
Genoma , Genómica , Humanos , Bovinos/genética , Animales , Fenotipo , Aumento de Peso/genética , Genotipo , Ingestión de Alimentos/genética , Alimentación AnimalRESUMEN
The aims of this study were to: (1) estimate genetic correlation for milk production traits (milk, fat and protein yields and fat and protein contents) and fatty acids (FA: C16:0, C18:1 cis-9, LCFA, SFA, and UFA) over days in milk, (2) investigate the performance of genomic predictions using single-step GBLUP (ssGBLUP) based on random regression models (RRM), and (3) identify the optimal scaling and weighting factors to be used in the construction of the H matrix. A total of 302 684 test-day records of 63.875 first lactation Walloon Holstein cows were used. Positive genetic correlations were found between milk yield and fat and protein yield (rg from 0.46 to 0.85) and between fat yield and milk FA (rg from 0.17 to 0.47). On the other hand, negative correlations were estimated between fat and protein contents (rg from -0.22 to -0.59), between milk yield and milk FA (rg from -0.22 to -0.62), and between protein yield and milk FA (rg from -0.11 to -0.19). The selection for high fat content increases milk FA throughout lactation (rg from 0.61 to 0.98). The test-day ssGBLUP approach showed considerably higher prediction reliability than the parent average for all milk production and FA traits, even when no scaling and weighting factors were used in the H matrix. The highest validation reliabilities (r2 from 0.09 to 0.38) and less biased predictions (b1 from 0.76 to 0.92) were obtained using the optimal parameters (i.e., ω = 0.7 and α = 0.6) for the genomic evaluation of milk production traits. For milk FA, the optimal parameters were ω = 0.6 and α = 0.6. However, biased predictions were still observed (b1 from 0.32 to 0.81). The findings suggest that using ssGBLUP based on RRM is feasible for the genomic prediction of daily milk production and FA traits in Walloon Holstein dairy cattle.
RESUMEN
Digital image analysis is a practical, non-invasive, and relatively low-cost tool that may assist in the evaluation of body traits in Nile tilapia, being particularly useful for assessing difficult-to-measure variables, such as body areas. In this study, we aimed to estimate variance components and genetic parameters for body areas of Nile tilapia obtained by digital images. The data set comprised body weight (BW) records of 1,917 pond-reared fish at 366 days of age. Of this total, 656 animals were photographed and subjected to image analysis of trunk area (TA), head area (HA), caudal fin area (CFA) and fillet area (FA). Heritabilities and genetic correlations were estimated through multiple-trait models based on Bayesian inference. Heritability estimates for BW, TA, HA, CFA and FA were 0.25, 0.23, 0.26, 0.21 and 0.25, respectively. Genetic correlations between the traits were high and positive, ranging from 0.70 to 0.98. We highlight the genetic correlation between BW and TA (rG = 0.98) and FA (rG = 0.97). In view of the observed results, it can be concluded that trunk and fillet areas obtained by digital image analysis can lead to indirect genetic gains in weight and other body areas. In addition, the areas studied have potential as a selection criterion and may assist in studies on changes in the body shape in Nile tilapia.
Asunto(s)
Cíclidos , Animales , Teorema de Bayes , Peso Corporal , Cíclidos/genética , FenotipoRESUMEN
Nellore is the main cattle breed used in Brazil, being the largest commercial herd in the world. Beyond the importance of male reproductive efficiency for farm profit, the use of reproductive techniques, mainly artificial insemination, turns the evaluation of male reproductive traits even more important. Estimation of genetic parameters increases the knowledge on traits variances and allows envisaging the possibility of the inclusion of new traits as selection criterion. Genetic parameters for fifteen traits that can be classified as testicular biometry or physical and morphological semen traits were estimated for a Nellore bull population ranging from 18 to 36 months. Single-trait and bi-trait animal models were used for (co)variance components estimation. The contemporary group was considered as fixed effect and age at measurement as covariable. Scrotal circumference presented heritability of 0.47 ± 0.12. This value is similar to the heritabilities found for all testicular biometry traits (0.34-0.48). Sperm progressive motility, which has a direct effect on bull fertility, presented low heritability (0.07 ± 0.08). Major and total sperm defects presented moderate to high heritabilities (0.49 ± 0.18 and 0.39 ± 0.15, respectively), indicating that great genetic gain can be obtained through selection against sperm defects. High and positive genetic correlations were observed among testicular biometry traits, which also presented favourable genetic correlations with physical and morphological traits of the semen with magnitude ranging from high to low. Scrotal circumference presented moderate to high and favourable genetic correlations with sperm progressive motility, sperm turbulence, major sperm defects and total sperm defects. Thus, the selection for scrotal circumference results in favourable correlated genetic response for semen quality. The results show that the use of scrotal circumference as reference trait for bull fertility is appropriate, since it presents high heritability and favourable genetic correlation with semen quality.
Asunto(s)
Bovinos/genética , Fertilidad/genética , Testículo/anatomía & histología , Animales , Cruzamiento , Bovinos/anatomía & histología , Masculino , Carácter Cuantitativo Heredable , Escroto/anatomía & histología , Análisis de Semen/veterinaria , Motilidad Espermática/genética , Espermatozoides/anomalíasRESUMEN
The objectives of this study were to assess linkage disequilibrium (LD) and selection-induced changes in single nucleotide polymorphism (SNP) frequency, and to perform association mapping in popcorn chromosome regions containing quantitative trait loci (QTLs) for quality traits. Seven tropical and two temperate popcorn populations were genotyped for 96 SNPs chosen in chromosome regions containing QTLs for quality traits. The populations were phenotyped for expansion volume, 100-kernel weight, kernel sphericity, and kernel density. The LD statistics were the difference between the observed and expected haplotype frequencies (D), the proportion of D relative to the expected maximum value in the population, and the square of the correlation between the values of alleles at two loci. Association mapping was based on least squares and Bayesian approaches. In the tropical populations, D-values greater than 0.10 were observed for SNPs separated by 100-150 Mb, while most of the D-values in the temperate populations were less than 0.05. Selection for expansion volume indirectly led to increase in LD values, population differentiation, and significant changes in SNP frequency. Some associations were observed for expansion volume and the other quality traits. The candidate genes are involved with starch, storage protein, lipid, and cell wall polysaccharides synthesis.
RESUMEN
This study used qRT-PCR to examine variation in the expression of 13 myogenes during muscle development in four prenatal periods (21, 40, 70 and 90 days post-insemination) in commercial (the three-way Duroc, Landrace and Large-White cross) and local Piau pig breeds that differ in muscle mass. There was no variation in the expression of the CHD8, EID2B, HIF1AN, IKBKB, RSPO3, SOX7 and SUFU genes at the various prenatal ages or between breeds. The MAP2K1 and RBM24 genes showed similar expression between commercial and Piau pigs but greater expression (p < 0.05) in at least one prenatal period. Pair-wise comparisons of prenatal periods in each breed showed that only the CSRP3, LEF1, MRAS and MYOG genes had higher expression (p < 0.05) in at least one prenatal period in commercial and Piau pigs. Overall, these results identified the LEF1 gene as a primary candidate to account for differences in muscle mass between the pig breeds since activation of this gene may lead to greater myoblast fusion in the commercial breed compared to Piau pigs. Such fusion could explain the different muscularity between breeds in the postnatal periods.
RESUMEN
In the current post-genomic era, the genetic basis of pig growth can be understood by assessing SNP marker effects and genomic breeding values (GEBV) based on estimates of these growth curve parameters as phenotypes. Although various statistical methods, such as random regression (RR-BLUP) and Bayesian LASSO (BL), have been applied to genomic selection (GS), none of these has yet been used in a growth curve approach. In this work, we compared the accuracies of RR-BLUP and BL using empirical weight-age data from an outbred F2 (Brazilian Piau X commercial) population. The phenotypes were determined by parameter estimates using a nonlinear logistic regression model and the halothane gene was considered as a marker for evaluating the assumptions of the GS methods in relation to the genetic variation explained by each locus. BL yielded more accurate values for all of the phenotypes evaluated and was used to estimate SNP effects and GEBV vectors. The latter allowed the construction of genomic growth curves, which showed substantial genetic discrimination among animals in the final growth phase. The SNP effect estimates allowed identification of the most relevant markers for each phenotype, the positions of which were coincident with reported QTL regions for growth traits.
RESUMEN
Efficient computing techniques allow the estimation of variance components for virtually any traditional dataset. When genomic information is available, variance components can be estimated using genomic REML (GREML). If only a portion of the animals have genotypes, single-step GREML (ssGREML) is the method of choice. The genomic relationship matrix (G) used in both cases is dense, limiting computations depending on the number of genotyped animals. The algorithm for proven and young (APY) can be used to create a sparse inverse of G (GAPY~-1) with close to linear memory and computing requirements. In ssGREML, the inverse of the realized relationship matrix (H-1) also includes the inverse of the pedigree relationship matrix, which can be dense with a long pedigree, but sparser with short. The main purpose of this study was to investigate whether costs of ssGREML can be reduced using APY with truncated pedigree and phenotypes. We also investigated the impact of truncation on variance components estimation when different numbers of core animals are used in APY. Simulations included 150K animals from 10 generations, with selection. Phenotypes (h2 = 0.3) were available for all animals in generations 1-9. A total of 30K animals in generations 8 and 9, and 15K validation animals in generation 10 were genotyped for 52,890 SNP. Average information REML and ssGREML with G-1 and GAPY~-1 using 1K, 5K, 9K, and 14K core animals were compared. Variance components are impacted when the core group in APY represents the number of eigenvalues explaining a small fraction of the total variation in G. The most time-consuming operation was the inversion of G, with more than 50% of the total time. Next, numerical factorization consumed nearly 30% of the total computing time. On average, a 7% decrease in the computing time for ordering was observed by removing each generation of data. APY can be successfully applied to create the inverse of the genomic relationship matrix used in ssGREML for estimating variance components. To ensure reliable variance component estimation, it is important to use a core size that corresponds to the number of largest eigenvalues explaining around 98% of total variation in G. When APY is used, pedigrees can be truncated to increase the sparsity of H and slightly reduce computing time for ordering and symbolic factorization, with no impact on the estimates.
The estimation of variance components is computationally expensive under large-scale genetic evaluations due to several inversions of the coefficient matrix. Variance components are used as parameters for estimating breeding values in mixed model equations (MME). However, resulting breeding values are not Best Linear Unbiased Predictions (BLUP) unless the variance components approach the true parameters. The increasing availability of genomic data requires the development of new methods for improving the efficiency of variance component estimations. Therefore, this study aimed to reduce the costs of single-step genomic REML (ssGREML) with the Algorithm for Proven and Young (APY) for estimating variance components with truncated pedigree and phenotypes using simulated data. In addition, we investigated the influence of truncation on variance components and genetic parameter estimates. Under APY, the size of the core group influences the similarity of breeding values and their reliability compared to the full genomic matrix. In this study, we found that to ensure reliable variance component estimation, it is required to consider a core size that corresponds to the number of largest eigenvalues explaining around 98% of the total variation in G to avoid biased parameters. In terms of costs, the use of APY slightly decreased the time for ordering and symbolic factorization with no impact on estimations.
Asunto(s)
Genoma , Modelos Genéticos , Algoritmos , Animales , Genómica/métodos , Genotipo , Linaje , FenotipoRESUMEN
OBJECTIVE: The aim of this study was to estimate genetic parameters for 305-day cumulative milk yield and components, growth, and reproductive traits in Guzerá cattle. METHODS: The evaluated traits were 305-day first-lactation cumulative yields (kg) of milk (MY305), fat (FY305), protein (PY305), lactose (LY305), and total solids (SY305); age at first calving (AFC) in days; adjusted scrotal perimeter (cm) at the ages of 365 (SP365) and 450 (SP450) days; and adjusted body weight (kg) at the ages of 210 (W210), 365 (W365), and 450 (W450) days. The (co)variance components were estimated using the restricted maximum likelihood method for single-trait, bi-trait and tri-trait analyses. Contemporary groups and additive genetic effects were included in the general mixed model. Maternal genetic and permanent environmental effects were also included for W210. RESULTS: The direct heritability estimates ranged from 0.16 (W210) to 0.32 (MY305). The maternal heritability estimate for W210 was 0.03. Genetic correlation estimates among milk production traits and growth traits ranged from 0.92 to 0.99 and from 0.92 to 0.99, respectively. For milk production and growth traits, the genetic correlations ranged from 0.33 to 0.56. The genetic correlations among AFC and all other traits were negative (-0.43 to -0.27). Scrotal perimeter traits and body weights showed genetic correlations ranging from 0.41 to 0.46, and scrotal perimeter and milk production traits showed genetic correlations ranging from 0.11 to 0.30. The phenotypic correlations were similar in direction (same sign) and lower than the corresponding genetic correlations. CONCLUSION: These results suggest the viability and potential of joint selection for dairy and beef traits in Guzerá cattle, taking into account reproductive traits.
RESUMEN
This study aimed to evaluate the effect of parity order on milk yield (MY) and composition over time of grazing beef cows and to evaluate non-linear models to describe the lactation curve. Thirty-six pregnant Nellore cows (12 nulliparous, 2 years; 12 primiparous, 3 years; and 12 multiparous, 4-6 years) were included in the study. With calving day assigned as day 0, milking was performed using a milking machine to estimate MY on days 7, 14, 21, 42, 63, 91, 119, 154, and 203. Dummy variable analyses were applied to estimate its effects on MY, composition (kg and percentage), afternoon/morning, and afternoon/total proportions. Since multiparous cows had higher MY than nulliparous and primiparous cows, two different groups were used for lactation curve analysis: Mult (multiparous) and Null/Prim (nulliparous and primiparous). The MY estimated by the last edition of BR-Corte (Nutrient Requirements of Zebu and Crossbred Cattle) equation was compared with the observed values from this study. Five nonlinear models proposed by Wood (WD), Jenkins & Ferrell (JF), Wilmink (WK), Henriques (HR) and Cobby & Le Du (CL) were evaluated. Models were validated using an independent dataset of multiparous and primiparous cows. The estimates for parameters a, b, and c of the CL equation were compared between groups, and the BR-Corte equation used the model identity methodology. Nulliparous and primiparous cows displayed similar MY (P > 0.05); however, multiparous cows had an average MY that is 0.70 kg/day greater than that of nulliparous and primiparous cows (P < 0.05). Milk protein and total solids were higher for multiparous cows (P < 0.05). Effect of days in milking was found for milk fat, protein, and total solids (P < 0.05). The yield of all milk components was higher for multiparous cows than for nulliparous and primiparous cows. The afternoon/morning and afternoon/total proportions of milk production were not affected by parities and days in milking (P > 0.05), with an average of 0.76 and 0.42, respectively. The BR-Corte equation did not correctly estimate the MY (P < 0.05). The equations of WD, WK, and CL had the best estimate of MY for both Mult and Null/Prim datasets. The equations had a very similar Akaike's information criterion with correction and mean square error of prediction.
RESUMEN
To evaluate the effect of an Escherichia coli lipopolysaccharide (LPS) challenge on the digestible lysine (Lys) requirement for growing pigs, a nitrogen (N) balance assay was performed. Seventy-two castrated male pigs (19 ± 1.49 kg body weight [BW]) were allocated in a 2 × 6 factorial design composed of two immune activation states (control and LPS-challenged) and six dietary treatments with N levels of 0.94, 1.69, 2.09, 3.04, 3.23, and 3.97% N, as fed, where Lys was limiting, with six replicates and one pig per unit. The challenge consisted of an initial LPS dose of 30 µg/kg BW via intramuscular (IM) injection and a subsequent dose of 33.6 µg/kg BW after 48 h. The experimental period lasted 11 d and was composed of a 7-d adaptation and a subsequent 4-d sampling period in which N intake (NI), N excretion (NEX), and N deposition (ND) were evaluated. Inflammatory mediators and rectal temperature were assessed during the 4-d collection period. A three-way interaction (N levels × LPS challenge × time, P < 0.05) for IgG was observed. Additionally, two-way interactions (challenge × time, P < 0.05) were verified for IgA, ceruloplasmin, transferrin, haptoglobin, α-1-acid glycoprotein, total protein, and rectal temperature; and (N levels × time, P < 0.05) for transferrin, albumin, haptoglobin, total protein, and rectal temperature. LPS-challenged pigs showed lower (P < 0.05) feed intake. A two-way interaction (N levels × LPS challenge, P < 0.05) was observed for NI, NEX, and ND, with a clear dose-response (P < 0.05). LPS-challenged pigs showed lower NI and ND at 2.09% N and 1.69 to 3.97% N (P < 0.05), respectively, and higher NEX at 3.23% N (P < 0.05). The parameters obtained by a nonlinear model (N maintenance requirement, NMR and theoretical maximum N deposition, NDmaxT) were 152.9 and 197.1 mg/BWkg0.75/d for NMR, and 3,524.7 and 2,077.8 mg/BWkg0.75/d for NDmaxT, for control and LPS-challenged pigs, respectively. The estimated digestible Lys requirements were 1,994.83 and 949.16 mg/BWkg0.75/d for control and LPS-challenged pigs, respectively. The daily digestible Lys intakes required to achieve 0.68 and 0.54 times the NRmaxT value were 18.12 and 8.62 g/d, respectively, and the optimal dietary digestible Lys concentration may change depending on the feed intake levels. Based on the derived model parameters obtained in the N balance trial with lower cost and time, it was possible to differentiate the digestible Lys requirement for swine under challenging conditions.
Asunto(s)
Alimentación Animal , Lisina , Alimentación Animal/análisis , Fenómenos Fisiológicos Nutricionales de los Animales , Animales , Dieta/veterinaria , Ingestión de Alimentos , Lipopolisacáridos , Masculino , PorcinosRESUMEN
Pig survival is an economically important trait with relevant social welfare implications, thus standing out as an important selection criterion for the current pig farming system. We aimed to estimate (co)variance components for survival in different production phases in a crossbred pig population as well as to investigate the benefit of including genomic information through single-step genomic best linear unbiased prediction (ssGBLUP) on the prediction accuracy of survival traits compared with results from traditional BLUP. Individual survival records on, at most, 64,894 crossbred piglets were evaluated under two multi-trait threshold models. The first model included farrowing, lactation, and combined postweaning survival, whereas the second model included nursery and finishing survival. Direct and maternal breeding values were estimated using BLUP and ssGBLUP methods. Furthermore, prediction accuracy, bias, and dispersion were accessed using the linear regression validation method. Direct heritability estimates for survival in all studied phases were low (from 0.02 to 0.08). Survival in preweaning phases (farrowing and lactation) was controlled by the dam and piglet additive genetic effects, although the maternal side was more important. Postweaning phases (nursery, finishing, and the combination of both) showed the same or higher direct heritabilities compared with preweaning phases. The genetic correlations between survival traits within preweaning and postweaning phases were favorable and strong, but correlations between preweaning and postweaning phases were moderate. The prediction accuracy of survival traits was low, although it increased by including genomic information through ssGBLUP compared with the prediction accuracy from BLUP. Direct and maternal breeding values were similarly accurate with BLUP, but direct breeding values benefited more from genomic information. Overall, a slight increase in bias was observed when genomic information was included, whereas dispersion of breeding values was greatly reduced. Combined postweaning survival presented higher direct heritability than in the preweaning phases and the highest prediction accuracy among all evaluated production phases, therefore standing out as a candidate trait for improving survival. Survival is a complex trait with low heritability; however, important genetic gains can still be obtained, especially under a genomic prediction framework.
Asunto(s)
Genoma , Modelos Genéticos , Animales , Femenino , Genómica , Genotipo , Linaje , Fenotipo , Porcinos/genéticaRESUMEN
Pedigree information is incomplete by nature and commonly not well-established because many of the genetic ties are not known a priori or can be wrong. The genomic era brought new opportunities to assess relationships between individuals. However, when pedigree and genomic information are used simultaneously, which is the case of single-step genomic BLUP (ssGBLUP), defining the genetic base is still a challenge. One alternative to overcome this challenge is to use metafounders, which are pseudo-individuals that describe the genetic relationship between the base population individuals. The purpose of this study was to evaluate the impact of metafounders on the estimation of breeding values for tick resistance under ssGBLUP for a multibreed population composed by Hereford, Braford, and Zebu animals. Three different scenarios were studied: pedigree-based model (BLUP), ssGBLUP, and ssGBLUP with metafounders (ssGBLUPm). In ssGBLUPm, a total of four different metafounders based on breed of origin (i.e., Hereford, Braford, Zebu, and unknown) were included for the animals with missing parents. The relationship coefficient between metafounders was in average 0.54 (ranging from 0.34 to 0.96) suggesting an overlap between ancestor populations. The estimates of metafounder relationships indicate that Hereford and Zebu breeds have a possible common ancestral relationship. Inbreeding coefficients calculated following the metafounder approach had less negative values, suggesting that ancestral populations were large enough and that gametes inherited from the historical population were not identical. Variance components were estimated based on ssGBLUPm, ssGBLUP, and BLUP, but the values from ssGBLUPm were scaled to provide a fair comparison with estimates from the other two models. In general, additive, residual, and phenotypic variance components in the Hereford population were smaller than in Braford across different models. The addition of genomic information increased heritability for Hereford, possibly because of improved genetic relationships. As expected, genomic models had greater predictive ability, with an additional gain for ssGBLUPm over ssGBLUP. The increase in predictive ability was greater for Herefords. Our results show the potential of using metafounders to increase accuracy of GEBV, and therefore, the rate of genetic gain in beef cattle populations with partial levels of missing pedigree information.
RESUMEN
Identifying maize inbred lines that are more efficient in nitrogen (N) use is an important strategy and a necessity in the context of environmental and economic impacts attributed to the excessive N fertilization. N-uptake efficiency (NUpE) and N-utilization efficiency (NUtE) are components of N-use efficiency (NUE). Despite the most maize breeding data have a multi-trait structure, they are often analyzed under a single-trait framework. We aimed to estimate the genetic parameters for NUpE and NUtE in contrasting N levels, in order to identify superior maize inbred lines, and to propose a Bayesian multi-trait multi-environment (MTME) model. Sixty-four tropical maize inbred lines were evaluated in two experiments: at high (HN) and low N (LN) levels. The MTME model was compared to single-trait multi-environment (STME) models. Based on deviance information criteria (DIC), both multi- and single-trait models revealed genotypes x environments (G x E) interaction. In the MTME model, NUpE was found to be weakly heritable with posterior modes of heritability of 0.016 and 0.023 under HN and LN, respectively. NUtE at HN was found to be highly heritable (0.490), whereas under LN condition it was moderately heritable (0.215). We adopted the MTME model, since combined analysis often presents more accurate breeding values than single models. Superior inbred lines for NUpE and NUtE were identified and this information can be used to plan crosses to obtain maize hybrids that have superior nitrogen use efficiency.
Asunto(s)
Teorema de Bayes , Ambiente , Nitrógeno/metabolismo , Carácter Cuantitativo Heredable , Zea mays/genética , Zea mays/metabolismo , Algoritmos , Modelos Estadísticos , Clima TropicalRESUMEN
Nonlinear mixed models were used to describe longitudinal scrotal circumference (SC) measurements of Nellore bulls. Models comparisons were based on Akaike's information criterion, Bayesian information criterion, error sum of squares, adjusted R2 and percentage of convergence. Sequentially, the best model was used to compare the SC growth curve in bulls divergently classified according to SC at 18-21 months of age. For this, bulls were classified into five groups: SC < 28cm; 28cm ≤ SC < 30cm, 30cm ≤ SC < 32cm, 32cm ≤ SC < 34cm and SC ≥ 34cm. Michaelis-Menten model showed the best fit according to the mentioned criteria. In this model, ß1 is the asymptotic SC value and ß2 represents the time to half-final growth and may be related to sexual precocity. Parameters of the individual estimated growth curves were used to create a new dataset to evaluate the effect of the classification, farms, and year of birth on ß1 and ß2 parameters. Bulls of the largest SC group presented a larger predicted SC along all analyzed periods; nevertheless, smaller SC group showed predicted SC similar to intermediate SC groups (28cm ≤ SC < 32cm), around 1200 days of age. In this context, bulls classified as improper for reproduction at 18-21 months old can reach a similar condition to those considered as good condition. In terms of classification at 18-21 months, asymptotic SC was similar among groups, farms and years; however, ß2 differed among groups indicating that differences in growth curves are related to sexual precocity. In summary, it seems that selection based on SC at too early ages may lead to discard bulls with suitable reproductive potential.
Asunto(s)
Bovinos/anatomía & histología , Escroto/anatomía & histología , Animales , Teorema de Bayes , Cruzamiento/métodos , Bovinos/fisiología , Masculino , Dinámicas no Lineales , Reproducción , Escroto/fisiología , Maduración SexualRESUMEN
This study aimed to evaluate the correlation between testicular biometry and semen variables, as well as, to relate testicular variables to the probability of selecting Nellore bulls with desirable sperm morphology when conducting breeding soundness evaluations (BSE). A total of 2055 BSEs from 506 bulls comprised the dataset. Biometric variables evaluated were: scrotal circumference, testicular volume, width, length, ratio and eccentricity; and semen variables were sperm motility, major sperm defects, minor sperm defects and normal sperm. Data of testicular biometry were correlated with data for semen variables using the Pearson's correlation assessment. Effects of testicular variables in selecting for sperm morphology of bulls in the BSE were evaluated by logistic regression. Scrotal circumference, testicular volume, length and width were positively correlated to sperm motility (0.18 to 0.19) and normal sperm (0.24 to 0.27) and negatively correlated with values for major defects (-0.24 to -0.27), but for testicular ratio and eccentricity there were coefficients near zero for all semen traits. Testicular ratio and eccentricity were not suitable for predicting the probability of selecting a bull based on semen variables using the BSE, but scrotal circumference, testicular volume, length and width were highly significant (P < 0.0001) with moderate values of area under ROC (Receiver Operating Characteristics) curve (0.608 to 0.620).
Asunto(s)
Cruzamiento , Bovinos/fisiología , Semen/fisiología , Testículo/fisiología , Animales , Biometría , Masculino , Escroto , Motilidad Espermática , EspermatozoidesRESUMEN
Flowering is an important agronomic trait. Quantile regression (QR) can be used to fit models for all portions of a probability distribution. In Genome-wide association studies (GWAS), QR can estimate SNP (Single Nucleotide Polymorphism) effects on each quantile of interest. The objectives of this study were to estimate genetic parameters and to use QR to identify genomic regions for phenological traits (Days to first flower-DFF; Days for flowering-DTF; Days to end of flowering-DEF) in common bean. A total of 80 genotypes of common beans, with 3 replicates were raised at 4 locations and seasons. Plants were genotyped for 384 SNPs. Traditional single-SNP and 9 QR models, ranging from equally spaced quantiles (τ) 0.1 to 0.9, were used to associate SNPs to phenotype. Heritabilities were moderate high, ranging from 0.32 to 0.58. Genetic and phenotypic correlations were all high, averaging 0.66 and 0.98, respectively. Traditional single-SNP GWAS model was not able to find any SNP-trait association. On the other hand, when using QR methodology considering one extreme quantile (τ = 0.1) we found, respectively 1 and 7, significant SNPs associated for DFF and DTF. Significant SNPs were found on Pv01, Pv02, Pv03, Pv07, Pv10 and Pv11 chromosomes. We investigated potential candidate genes in the region around these significant SNPs. Three genes involved in the flowering pathways were identified, including Phvul.001G214500, Phvul.007G229300 and Phvul.010G142900.1 on Pv01, Pv07 and Pv10, respectively. These results indicate that GWAS-based QR was able to enhance the understanding on genetic architecture of phenological traits (DFF and DTF) in common bean.
Asunto(s)
Fabaceae/genética , Flores/genética , Genes de Plantas , Estudio de Asociación del Genoma Completo , Polimorfismo de Nucleótido SimpleRESUMEN
Gene expression time series (GETS) analysis aims to characterize sets of genes according to their longitudinal patterns of expression. Due to the large number of genes evaluated in GETS analysis, an useful strategy to summarize biological functional processes and regulatory mechanisms is through clustering of genes that present similar expression pattern over time. Traditional cluster methods usually ignore the challenges in GETS, such as the lack of data normality and small number of temporal observations. Independent Component Analysis (ICA) is a statistical procedure that uses a transformation to convert raw time series data into sets of values of independent variables, which can be used for cluster analysis to identify sets of genes with similar temporal expression patterns. ICA allows clustering small series of distribution-free data while accounting for the dependence between subsequent time-points. Using temporal simulated and real (four libraries of two pig breeds at 21, 40, 70 and 90 days of gestation) RNA-seq data set we present a methodology (ICAclust) that jointly considers independent components analysis (ICA) and a hierarchical method for clustering GETS. We compare ICAclust results with those obtained for K-means clustering. ICAclust presented, on average, an absolute gain of 5.15% over the best K-means scenario. Considering the worst scenario for K-means, the gain was of 84.85%, when compared with the best ICAclust result. For the real data set, genes were grouped into six distinct clusters with 89, 51, 153, 67, 40, and 58 genes each, respectively. In general, it can be observed that the 6 clusters presented very distinct expression patterns. Overall, the proposed two-step clustering method (ICAclust) performed well compared to K-means, a traditional method used for cluster analysis of temporal gene expression data. In ICAclust, genes with similar expression pattern over time were clustered together.
Asunto(s)
Algoritmos , Perfilación de la Expresión Génica , Análisis de Secuencia por Matrices de Oligonucleótidos , ARN , Animales , Análisis por Conglomerados , Simulación por Computador , Perfilación de la Expresión Génica/métodos , Regulación del Desarrollo de la Expresión Génica , Modelos Moleculares , Músculo Esquelético/embriología , Músculo Esquelético/metabolismo , Análisis de Secuencia por Matrices de Oligonucleótidos/métodos , ARN/metabolismo , Porcinos , Factores de TiempoRESUMEN
ABSTRACT: The objective of this study was to adjust nonlinear quantile regression models for the study of dry matter accumulation in garlic plants over time, and to compare them to models fitted by the ordinary least squares method. The total dry matter of nine garlic accessions belonging to the Vegetable Germplasm Bank of Universidade Federal de Viçosa (BGH/UFV) was measured in four stages (60, 90, 120 and 150 days after planting), and those values were used for the nonlinear regression models fitting. For each accession, there was an adjustment of one model of quantile regression (τ=0.5) and one based on the least squares method. The nonlinear regression model fitted was the Logistic. The Akaike Information Criterion was used to evaluate the goodness of fit of the models. Accessions were grouped using the UPGMA algorithm, with the estimates of the parameters with biological interpretation as variables. The nonlinear quantile regression is efficient for the adjustment of models for dry matter accumulation in garlic plants over time. The estimated parameters are more uniform and robust in the presence of asymmetry in the distribution of the data, heterogeneous variances, and outliers.
RESUMO: Este trabalho teve como objetivo ajustar modelos de regressão quantílica não linear para o estudo do acúmulo de matéria seca total em plantas de alho ao longo do tempo, e compará-los com modelos ajustados pelo método dos mínimos quadrados. A matéria seca total de nove acessos de alho pertencentes ao Banco de Germoplasma de Hortaliças da Universidade Federal de Viçosa (BGH/UFV) foi avaliada em quatro períodos (60, 90, 120 e 150 dias após plantio), e estes valores foram utilizados para o ajuste de modelos de regressão - não linear - logística. Para cada acesso, foram ajustados um modelo de regressão quantílica (τ=0,5) e um modelo pela metodologia dos mínimos quadrados. Para avaliar a qualidade de ajuste dos modelos foi utilizado o Critério de Informação de Akaike. Os acessos foram agrupados pelo algoritmo UPGMA, utilizando as estimativas dos parâmetros com interpretação biológica como variáveis. A regressão quantílica não linear foi eficiente no ajuste de modelos para descrição do acúmulo de matéria seca ao longo do tempo. As estimativas de parâmetros foram mais uniformes e robustas na presença de dados assimétricos, variâncias heterogêneas e de valores discrepantes.
RESUMEN
Canine soft tissue sarcomas (STS) comprise a heterogeneous group of malignancies that share similar histopathological features, a low to moderate recurrence rate and low metastatic potential. In human medicine, the expression of estrogen receptors (ER) and progesterone receptors (PR) in sarcomas has been studied to search for prognostic factors and new treatment targets. Similar studies have yet to be conducted in veterinary medicine. The objective of this study was therefore to investigate by immunohistochemistry (IHC) the ER and PR expression in a series of 80 cutaneous and subcutaneous sarcomas in dogs with histopathological features of peripheral nerve sheath tumor (PNST) and perivascular wall tumor (PWT). All cases were positive for PR and negative for ER. Tumors of high malignancy grade (grade III) exhibited higher PR expression than low-grade tumors (grade I). Tumors with mitotic activity greater than 9 mitotic figures/10 high power fields also exhibited higher PR expression. In addition, there was a positive correlation between cell proliferation (Ki67) and PR expression. Therefore, it is possible that progesterone plays a greater role than estrogen in the pathogenesis of these tumors. Future studies should explore the potential for selective progesterone receptor modulators as therapeutic agents in canine STS, as well as evaluating PR expression as a predictor of prognosis.(AU)
Sarcomas de tecidos moles (STM) caninos compreendem um grupo heterogêneo de neoplasias malignas, que apresentam alterações histopatológicas similares, baixa a moderada taxa de recorrência e baixo potencial metastático. Em medicina humana, a expressão de receptor para estrógeno (RE) e receptor para progesterona (RP) nos sarcomas tem sido estudada, visando a busca por fatores prognósticos e novos alvos para tratamentos. Na medicina veterinária, ainda não foram realizados estudos similares. O objetivo deste trabalho foi investigar por imuno-histoquímica a expressão de RE e RP em uma série de 80 sarcomas cutâneos e subcutâneos de cães, com características histopatológicas de tumor de bainha de nervo periférico e tumor de parede perivascular. Todos os casos foram positivos para RP e negativos para RE. Tumores de alto grau de malignidade (grau III) exibiram maior expressão deste receptor que os tumores de baixo grau (grau I). Tumores com atividade mitótica maior que 9 figuras mitóticas/10 campos de grande aumento também exibiram maior expressão do RP. Em adição, houve correlação positiva entre o índice de proliferação celular (Ki67) e a expressão de RP. Assim, é possível que a progesterona desempenhe maior papel que o estrógeno na patogênese desses tumores. Futuros trabalhos poderão explorar o potencial dos moduladores seletivos de RP como agente terapêutico em STM caninos, bem como avaliar a expressão de RP como preditiva de prognóstico.(AU)