Search

This study identified potential biomarkers for feed efficiency by blood plasma proteome analysis of a tropically adapted beef cattle breed. Two experimental groups were selected based on residual feed intake (RFI). The proteome was investigated by LC-MS/MS in a data-dependent acquisition mode. After quality control, 123 differentially abundant proteins (DAPs) were identified between the two experimental groups. Among DAPs with the highest absolute log-fold change values, the PRDM2, KRT5, UGGT1, DENND5B, B2M, SLC44A2, SLC7A2, PTPRC, and FETUB were highlighted as potential biomarkers because of their functions that may contribute to RFI. Furthermore, functional enrichment analysis revealed several biological processes, molecular functions and pathways that contributes to RFI, such as cell signaling, cellular responses to stimuli, immune system, calcium, hormones, metabolism and functions of proteins, lipids and carbohydrates. Protein-protein interaction analysis identified 32 and 11 DAPs as important nodes based on their interactions in the high- and low-RFI groups, respectively. This study represents the first comprehensive profiling of the blood plasma proteome of a tropically adapted beef cattle breed and provides valuable insights into the potential roles of these DAPs in key biological processes and pathways, contributing to our understanding of the mechanisms underlying feed efficiency in tropically adapted beef cattle. SIGNIFICANCE: LC-MS/MS analysis was performed to investigate changes in the blood plasma proteome associated with residual feed intake (RFI) in a tropically adapted beef cattle breed (Bos taurus taurus). Some putative biomarkers were identified to distinguish the high-RFI to low-RFI animals, based on their log-fold change value or on their protein-protein interaction network, which provide helpful sources in developing novel selection strategies for breeding programs. Our findings also revealed valuable insights into the metabolic pathways and biological processes that contribute to RFI in beef cattle, such as those closely linked to cell signaling, cellular responses to stimuli, immune system, calcium, hormones, metabolism and functions of proteins, lipids and carbohydrates., Competing Interests: Declaration of competing interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Background: Heat stress has deleterious effects on physiological and performance traits in livestock. Within this context, using tropically adapted cattle breeds in pure herds or terminal crossbreeding schemes to explore heterosis is attractive for increasing animal production in warmer climate regions. This study aimed to identify biological processes, pathways, and potential biomarkers related to thermotolerance in Caracu, a tropically adapted beef cattle breed, by proteomic analysis of blood plasma. To achieve this goal, 61 bulls had their thermotolerance evaluated through a heat tolerance index. A subset of 14 extreme animals, including the seven most thermotolerant (HIGH group) and the seven least thermotolerant (LOW group), had their blood plasma samples used for proteomic analysis by liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). The differentially regulated proteins detected between HIGH and LOW groups were used to perform functional enrichment analysis and a protein-protein interaction network analysis., Results: A total of 217 proteins were detected only in the HIGH thermotolerant group and 51 only in the LOW thermotolerant group. In addition, 81 and 87 proteins had significantly higher and lower abundancies in the HIGH group, respectively. Regarding proteins with the highest absolute log-fold change values, we highlighted those encoded by DUSP5, IGFALS, ROCK2, RTN4, IRAG1, and NNT genes based on their functions. The functional enrichment analysis detected several biological processes, molecular functions, and pathways related to cellular responses to stress, immune system, complement system, and hemostasis in both HIGH and LOW groups, in addition to terms and pathways related to lipids and calcium only in the HIGH group. Protein-protein interaction (PPI) network revealed as important nodes many proteins with roles in response to stress, hemostasis, immune system, inflammation, and homeostasis. Additionally, proteins with high absolute log-fold change values and proteins detected as essential nodes by PPI analysis highlighted herein are potential biomarkers for thermotolerance, such as ADRA1A, APOA1, APOB, APOC3, C4BPA, CAT, CFB, CFH, CLU, CXADR, DNAJB1, DNAJC13, DUSP5, FGA, FGB, FGG, HBA, HBB, HP, HSPD1, IGFALS, IRAG1, KNG1, NNT, OSGIN1, PROC, PROS1, ROCK2, RTN4, RYR1, TGFB2, VLDLR, VTN, and VWF., Conclusions: Identifying potential biomarkers, molecular mechanisms and pathways that act in response to heat stress in tropically adapted beef cattle contributes to developing strategies to improve performance and welfare traits in livestock under tropical climates., Competing Interests: Declarations Ethics approval All the experimental procedures were approved by the Institutional Animal Care and Use Committee at the Animal Science Institute, Nova Odessa, São Paulo, Brazil (CEUA n°292 − 19, October 7th, 2019), following guidelines for animal welfare according to State Law No. 11.977 of the São Paulo state, Brazil. Consent for publication Not applicable. Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Background: Identifying molecular mechanisms responsible for the response to heat stress is essential to increase production, reproduction, health, and welfare. This study aimed to identify early biological responses and potential biomarkers involved in the response to heat stress and animal's recovery in tropically adapted beef cattle through proteomic analysis of blood plasma., Methods: Blood samples were collected from 14 Caracu males during the heat stress peak (HSP) and 16 h after it (heat stress recovery-HSR) assessed based on wet bulb globe temperature index and rectal temperature. Proteome was investigated by liquid chromatography-tandem mass spectrometry from plasma samples, and the differentially regulated proteins were evaluated by functional enrichment analysis using DAVID tool. The protein-protein interaction network was evaluated by STRING tool., Results: A total of 1,550 proteins were detected in both time points, of which 84 and 65 were downregulated and upregulated during HSR, respectively. Among the differentially regulated proteins with the highest absolute log-fold change values, those encoded by the GABBR1, EPHA2, DUSP5, MUC2, DGCR8, MAP2K7, ADRA1A, CXADR, TOPBP1, and NEB genes were highlighted as potential biomarkers because of their roles in response to heat stress. The functional enrichment analysis revealed that 65 Gene Ontology terms and 34 pathways were significant ( P < 0.05). We highlighted those that could be associated with the response to heat stress, such as those related to the immune system, complement system, hemostasis, calcium, ECM-receptor interaction, and PI3K-Akt and MAPK signaling pathways. In addition, the protein-protein interaction network analysis revealed several complement and coagulation proteins and acute-phase proteins as important nodes based on their centrality and edges., Conclusion: Identifying differentially regulated proteins and their relationship, as well as their roles in key pathways contribute to improve the knowledge of the mechanisms behind the response to heat stress in naturally adapted cattle breeds. In addition, proteins highlighted herein are potential biomarkers involved in the early response and recovery from heat stress in tropically adapted beef cattle., Competing Interests: Author RS was employed by Angus Genetics Inc. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Reolon, Abduch, Freitas, Silva, Fragomeni, Lourenco, Baldi, Paz and Stafuzza.)

Infectious hematopoietic necrosis (IHN) is a disease of salmonid fish that is caused by the IHN virus (IHNV), which can cause substantial mortality and economic losses in rainbow trout aquaculture and fisheries enhancement hatchery programs. In a previous study on a commercial rainbow trout breeding line that has undergone selection, we found that genetic resistance to IHNV is controlled by the oligogenic inheritance of several moderate and many small effect quantitative trait loci (QTL). Here we used genome wide association analyses in two different commercial aquaculture lines that were naïve to previous exposure to IHNV to determine whether QTL were shared across lines, and to investigate whether there were major effect loci that were still segregating in the naïve lines. A total of 1,859 and 1,768 offspring from two commercial aquaculture strains were phenotyped for resistance to IHNV and genotyped with the rainbow trout Axiom 57K SNP array. Moderate heritability values (0.15-0.25) were estimated. Two statistical methods were used for genome wide association analyses in the two populations. No major QTL were detected despite the naïve status of the two lines. Further, our analyses confirmed an oligogenic architecture for genetic resistance to IHNV in rainbow trout. Overall, 17 QTL with notable effect (≥1.9% of the additive genetic variance) were detected in at least one of the two rainbow trout lines with at least one of the two statistical methods. Five of those QTL were mapped to overlapping or adjacent chromosomal regions in both lines, suggesting that some loci may be shared across commercial lines. Although some of the loci detected in this GWAS merit further investigation to better understand the biological basis of IHNV disease resistance across populations, the overall genetic architecture of IHNV resistance in the two rainbow trout lines suggests that genomic selection may be a more effective strategy for genetic improvement in this trait., Competing Interests: KM who is one of the co-authors of the manuscript is an employee of Troutlodge, Inc., a daughter company of Hendrix Genetics. Rainbow trout fish from the germplasm developed by Troutlodge were used in this study and their resistance to IHNV disease is being characterized in the study including the genetic architecture of this commercial germplasm. However, there is no endorsement in the manuscript of a commercial product from Troutlodge or Hendrix Genetics. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. government. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Palti, Vallejo, Purcell, Gao, Shewbridge, Long, Setzke, Fragomeni, Cheng, Martin and Naish.)

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research., (© 2023. The Author(s).)

The accuracy of genetic selection in dairy can be increased by the adoption of new technologies, such as the inclusion of sequence data. In simulation studies, assigning different weights to causative single-nucleotide polymorphism (SNP) markers led to better predictions depending on the genomic prediction method used. However, it is still not clear how the weights should be calculated. Our objective was to evaluate the accuracy of a multi-step method (GBLUP) and single-step GBLUP with simulated data using regular SNP, causatives variants (QTN) and the combination of both. Additionally, we compared the accuracies of all previous scenarios using alternatives for SNP weighting. The data were simulated assuming a single trait with a heritability of 0.3. The effective population size (Ne) was approximately 200. The pedigree contained 440,000 animals, and approximately 16,800 individuals were genotyped. A total of 49,974 SNP markers were evenly placed throughout the genome, and 100, 1000 and 2000 causative QTN were simulated. Both GBLUP and ssGBLUP were used in this study. We evaluated quadratic and nonlinear SNP weights in addition to the unweighted G. The inclusion of QTN to panels led to significant accuracy gains. Nonlinear A was demonstrated to be superior to quadratic weighting and unweighted approaches; however, results from Nonlinear A were dependent on the equation parameters. The unweighted approach was more suitable for less polygenic scenarios. Finally, SNP weighting might help elucidate trait architecture features based on changes in the accuracy of genomic prediction., (© 2023 Wiley-VCH GmbH. Published by John Wiley & Sons Ltd.)

This study aimed to identify genomic regions, pathways, and putative candidate genes associated with resistance to gastrointestinal nematode in Santa Ines sheep. The phenotypic information comprised 5529 records from 1703 naturally infected animals. After genomic data quality control, 37,511 SNPs from 589 animals were available. The weighted single-step approach for genome-wide association study was performed to estimate the SNP effects and variances accounted by 10-SNP sliding windows. Confirming the polygenic nature of the studied traits, 20, 22, 21, and 19 genomic windows that explained more than 0.5% of the additive genetic variance were identified for fecal egg counts (FEC), Famacha© (FAM), packed cell volume (PCV), and total plasma protein (TPP), respectively. A total of 81, 122, 106, and 101 protein-coding genes were found in windows associated with FEC, FAM, PCV, and TPP, respectively. Several protein-coding genes related to the immune system and inflammatory response functions were identified within those genomic regions, such as ADCY9, ADRB2, BRAF, CADM1, CCL20, CD70, CREBBP, FNBP1, HTR4, IL16, IL22, IL26, MAPK8, NDFIP1, NLRC3, PAK5, PLCB1, PLCB4, ROCK1, TEK, TNFRSF12A, and VAV1. Functional enrichment analysis by DAVID tool also revealed many significant (P < 0.05) pathways and Gene Ontology terms that could be related to resistance to gastrointestinal nematode in Santa Ines sheep, such as chemokine signaling pathway (oas04062), cAMP signaling pathway (oas04024), cGMP-PKG signaling pathway (Oas04022), platelet activation (Oas04611), Rap1 signaling pathway (oas04015), and oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen (GO:0016705). These results contribute to improving the knowledge of the genetic architecture of resistance to gastrointestinal nematode in Santa Ines sheep., Competing Interests: Declaration of Competing Interest None., (Copyright © 2023 Elsevier B.V. All rights reserved.)

This study evaluated the influence of environmental temperature on thermoregulation, hormonal, and hematological characteristics in Caracu cattle. Blood samples, hair length, coat and muzzle colors, rectal (RT), and surface temperatures were collected from 48 males and 43 females before (morning) and after sun exposure for eight hours (afternoon). Infrared thermography (IRT) was used to identify superficial temperature that exhibits a high correlation with RT. Hematological parameters, hormone concentrations, RT, and the superficial temperature obtained by IRT that exhibited the highest correlation with RT were evaluated by variance analysis. Regarding IRT, the lower left side of the body (LS) showed the highest correlation with the RT. Interaction between period and sex was observed for LS, cortisol, and eosinophils. Cortisone, progesterone, and RT were influenced by period and sex. Neutrophils and segmented neutrophils were influenced by the period, which showed the highest concentrations after sun exposure. Platelets, leukocytes, lymphocytes, and monocytes were influenced by sex. Heat stress changes several physiological characteristics where males and females exhibited differences in their responses to heat stress. Furthermore, most characteristics evaluated remained within the regular values observed for taurine Creole breeds, showing that Caracu is adapted to tropical climates.

Feed and water efficiency are important traits to improve beef cattle production’s economic and environmental sustainability. This study evaluated residual feed intake (RFI) and residual water intake (RWI) and their relationship with performance, ingestive behavior, and carcass traits in Caracu beef cattle. The data were analyzed using a generalized linear model with least squares means. The ingestive behavior, performance, and carcass traits were influenced by sex (p < 0.05). Males showed higher dry matter intake (DMI), average daily gain (ADG), mid-test metabolic weight (BW0.75), rib eye area, and rump fat thickness than females, besides spending more time drinking and eating. Low RFI animals exhibited higher DMI than high RFI animals. Low RWI animals ingested 3.89 L/d of water further than high RWI animals. The interaction between sex and RWI influenced the DMI, BW0.75, and backfat thickness. The ingestive behavior of low and high RFI animals was similar, although high RWI animals visited a smaller number of drinkers than low RWI animals. Water intake positively affects productive efficiency, and the combined use of RWI and RFI may help improve the selection of more efficient animals contributing to reducing the costs of beef cattle production and improving environmental sustainability.

This study investigated the prophylactic and therapeutic efficacies of baicalin (BC), a plant-derived flavone glycoside, in reducing the severity of Clostridioides difficile infection (CDI) in a mouse model. In the prophylactic trial, C57BL/6 mice were provided with BC (0, 11, and 22 mg/L in drinking water) from 12 days before C. difficile challenge through the end of the experiment, whereas BC administration started day 1 post challenge in the therapeutic trial. Both challenge and control groups were infected with 10 6 CFU/mL of hypervirulent C. difficile BAA 1803 spores or sterile PBS, and the clinical and diarrheal scores were recorded for 10 days post challenge. On day 2 post challenge, fecal and tissue samples were collected from mice prophylactically administered with BC for microbiome and histopathologic analysis. Both prophylactic and therapeutic supplementation of BC significantly reduced the severity of colonic lesions and improved CDI clinical progression and outcome compared with control ( p < 0.05). Microbiome analysis revealed a significant increase in Gammaproteobacteria and reduction in the abundance of protective microbiota (Firmicutes) in antibiotic-treated and C. difficile -infected mice compared with controls ( p < 0.05). However, baicalin supplementation favorably altered the microbiome composition, as revealed by an increased abundance in beneficial bacteria, especially Lachnospiraceae and Akkermansia . Our results warrant follow-up investigations on the use of BC as an adjunct to antibiotic therapy to control gut dysbiosis and reduce C. difficile infection in humans.

Infectious hematopoietic necrosis (IHN) is an economically important disease of salmonid fish caused by the IHN virus (IHNV). Under industrial aquaculture settings, IHNV can cause substantial mortality and losses. Actually, there is no confirmed and cost-effective method for IHNV control. Clear Springs Foods, Inc. has been performing family-based selective breeding to increase genetic resistance to IHNV in their rainbow trout breeding program. In an earlier study, we used siblings cross-validation to estimate the accuracy of genomic prediction (GP) for IHNV resistance in this breeding population. In the present report, we used empirical progeny testing data to evaluate whether genomic selection (GS) can improve the accuracy of breeding value predictions over traditional pedigree-based best linear unbiased predictions (PBLUP). We found that the GP accuracy with single-step GBLUP (ssGBLUP) outperformed PBLUP by 15% (from 0.33 to 0.38). Furthermore, we found that ssGBLUP had higher GP accuracy than weighted ssGBLUP (wssGBLUP) and single-step Bayesian multiple regression (ssBMR) models with BayesB and BayesC priors which supports our previous findings that the underlying liability of genetic resistance against IHNV in this breeding population might be polygenic. Our results show that GS can be more effective than either the traditional pedigree-based PBLUP model or the marker-assisted selection approach for improving genetic resistance against IHNV in this commercial rainbow trout population., (Copyright © 2020 Vallejo, Fragomeni, Cheng, Gao, Long, Shewbridge, MacMillan, Towner and Palti.)

Causal variants inferred from sequence data analysis are expected to increase accuracy of genomic selection. In this work we evaluated the gain in reliability of genomic predictions, for stature in US Holsteins, when adding selected sequence variants to a pre-existent SNP chip. Two prediction methods were tested: de-regressed proofs assuming heterogeneous (genomic BLUP; GBLUP) residual variances and by single-step GBLUP (ssGBLUP) using actual phenotypes. Phenotypic data included 3,999,631 records for stature on 3,027,304 Holstein cows. Genotypes on 54,087 SNP markers (54k) were available for 26,877 bulls. Additionally, 16,648 selected sequence variants were combined with the 54k markers, for a total of 70,735 (70k) markers. In all methods, SNP in the genomic relationship matrix (G) were unweighted or weighted iteratively, with weights derived either by SNP effects squared or by a nonlinear method that resembles BayesA (nonlinear A). Reliability of genomic predictions were obtained by cross validation. With unweighted G derived from 54k markers, the reliabilities (× 100) were 72.4 for GBLUP and 75.3 for ssGBLUP. With unweighted G derived from 70k markers, the reliabilities were 73.4 and 76.0, respectively. Weighting by nonlinear A changed reliabilities to 73.3, and 75.9, respectively. Addition of selected sequence variants had a small effect on reliabilities. Weighting by quadratic functions reduced reliabilities. Weighting by nonlinear A increased reliabilities for GBLUP but had only a small effect in ssGBLUP. Reliabilities for direct genomic values extracted from ssGBLUP using unweighted G with 54k were higher than reliabilities by any GBLUP. Thus, ssGBLUP seems to capture more information than GBLUP and there is less room for extra reliability. Improvements in GBLUP may be because the weights in G change the covariance structure, which can explain a proportion of the variance that is accounted for when a heterogeneous residual variance is assumed by considering a different number of daughters per bull., (Copyright © 2019 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Background: Infectious hematopoietic necrosis (IHN) is a disease of salmonid fish that is caused by the IHN virus (IHNV). Under intensive aquaculture conditions, IHNV can cause significant mortality and economic losses. Currently, there is no proven and cost-effective method for IHNV control. Clear Springs Foods, Inc. has been applying selective breeding to improve genetic resistance to IHNV in their rainbow trout breeding program. The goals of this study were to elucidate the genetic architecture of IHNV resistance in this commercial population by performing genome-wide association studies (GWAS) with multiple regression single-step methods and to assess if genomic selection can improve the accuracy of genetic merit predictions over conventional pedigree-based best linear unbiased prediction (PBLUP) using cross-validation analysis., Results: Ten moderate-effect quantitative trait loci (QTL) associated with resistance to IHNV that jointly explained up to 42% of the additive genetic variance were detected in our GWAS. Only three of the 10 QTL were detected by both single-step Bayesian multiple regression (ssBMR) and weighted single-step GBLUP (wssGBLUP) methods. The accuracy of breeding value predictions with wssGBLUP (0.33-0.39) was substantially better than with PBLUP (0.13-0.24)., Conclusions: Our comprehensive genome-wide scan for QTL revealed that genetic resistance to IHNV is controlled by the oligogenic inheritance of up to 10 moderate-effect QTL and many small-effect loci in this commercial rainbow trout breeding population. Taken together, our results suggest that whole genome-enabled selection models will be more effective than the conventional pedigree-based method for breeding value estimation or the marker-assisted selection approach for improving the genetic resistance of rainbow trout to IHNV in this population.

Background: In this study we integrated the CNV (copy number variation) and WssGWAS (weighted single-step approach for genome-wide association) analyses to increase the knowledge about number of piglets born alive, an economically important reproductive trait with significant impact on production efficiency of pigs., Results: A total of 3892 samples were genotyped with the Porcine SNP80 BeadChip. After quality control, a total of 57,962 high-quality SNPs from 3520 Duroc pigs were retained. The PennCNV algorithm identified 46,118 CNVs, which were aggregated by overlapping in 425 CNV regions (CNVRs) ranging from 2.5 Kb to 9718.4 Kb and covering 197 Mb (~ 7.01%) of the pig autosomal genome. The WssGWAS identified 16 genomic regions explaining more than 1% of the additive genetic variance for number of piglets born alive. The overlap between CNVR and WssGWAS analyses identified common regions on SSC2 (4.2-5.2 Mb), SSC3 (3.9-4.9 Mb), SSC12 (56.6-57.6 Mb), and SSC17 (17.3-18.3 Mb). Those regions are known for harboring important causative variants for pig reproductive traits based on their crucial functions in fertilization, development of gametes and embryos. Functional analysis by the Panther software identified 13 gene ontology biological processes significantly represented in this study such as reproduction, developmental process, cellular component organization or biogenesis, and immune system process, which plays relevant roles in swine reproductive traits., Conclusion: Our research helps to improve the understanding of the genetic architecture of number of piglets born alive, given that the combination of GWAS and CNV analyses allows for a more efficient identification of the genomic regions and biological processes associated with this trait in Duroc pigs. Pig breeding programs could potentially benefit from a more accurate discovery of important genomic regions.

We investigated the effects of different strategies for genotyping populations on variance components and heritabilities estimated with an animal model under restricted maximum likelihood (REML), genomic REML (GREML), and single-step GREML (ssGREML). A population with 10 generations was simulated. Animals from the last one, two or three generations were genotyped with 45,116 SNP evenly distributed on 27 chromosomes. Animals to be genotyped were chosen randomly or based on EBV. Each scenario was replicated five times. A single trait was simulated with three heritability levels (low, moderate, high). Phenotypes were simulated for only females to mimic dairy sheep and also for both sexes to mimic meat sheep. Variance component estimates from genomic data and phenotypes for one or two generations were more biased than from three generations. Estimates in the scenario without selection were the most accurate across heritability levels and methods. When selection was present in the simulations, the best option was to use genotypes of randomly selected animals. For selective genotyping, heritabilities from GREML were more biased compared to those estimated by ssGREML, because ssGREML was less affected by selective or limited genotyping., (© 2018 Blackwell Verlag GmbH.)

We investigated the importance of SNP weighting in populations with 2,000 to 25,000 genotyped animals. Populations were simulated with two effective sizes (20 or 100) and three numbers of QTL (10, 50 or 500). Pedigree information was available for six generations; phenotypes were recorded for the four middle generations. Animals from the last three generations were genotyped for 45,000 SNP. Single-step genomic BLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) were used to estimate genomic EBV using a genomic relationship matrix (G). The WssGBLUP performed better in small genotyped populations; however, any advantage for WssGBLUP was reduced or eliminated when more animals were genotyped. WssGBLUP had greater resolution for genome-wide association (GWA) as did increasing the number of genotyped animals. For few QTL, accuracy was greater for WssGBLUP than ssGBLUP; however, for many QTL, accuracy was the same for both methods. The largest genotyped set was used to assess the dimensionality of genomic information (number of effective SNP). The number of effective SNP was considerably less in weighted G than in unweighted G. Once the number of independent SNP is well represented in the genotyped population, the impact of SNP weighting becomes less important., (© 2017 Blackwell Verlag GmbH.)

The Algorithm for Proven and Young (APY) enables the implementation of single-step genomic BLUP (ssGBLUP) in large, genotyped populations by separating genotyped animals into core and non-core subsets and creating a computationally efficient inverse for the genomic relationship matrix (G). As APY became the choice for large-scale genomic evaluations in BLUP-based methods, a common question is how to choose the animals in the core subset. We compared several core definitions to answer this question. Simulations comprised a moderately heritable trait for 95,010 animals and 50,000 genotypes for animals across five generations. Genotypes consisted of 25,500 SNP distributed across 15 chromosomes. Genotyping errors and missing pedigree were also mimicked. Core animals were defined based on individual generations, equal representation across generations, and at random. For a sufficiently large core size, core definitions had the same accuracies and biases, even if the core animals had imperfect genotypes. When genotyped animals had unknown parents, accuracy and bias were significantly better (p ≤ .05) for random and across generation core definitions., (© 2017 The Authors. Journal of Animal Breeding and Genetics Published by Blackwell Verlag GmbH.)

Bacterial cold water disease (BCWD) causes significant mortality and economic losses in salmonid aquaculture. In previous studies, we identified moderate-large effect quantitative trait loci (QTL) for BCWD resistance in rainbow trout ( Oncorhynchus mykiss ). However, the recent availability of a 57 K SNP array and a reference genome assembly have enabled us to conduct genome-wide association studies (GWAS) that overcome several experimental limitations from our previous work. In the current study, we conducted GWAS for BCWD resistance in two rainbow trout breeding populations using two genotyping platforms, the 57 K Affymetrix SNP array and restriction-associated DNA (RAD) sequencing. Overall, we identified 14 moderate-large effect QTL that explained up to 60.8% of the genetic variance in one of the two populations and 27.7% in the other. Four of these QTL were found in both populations explaining a substantial proportion of the variance, although major differences were also detected between the two populations. Our results confirm that BCWD resistance is controlled by the oligogenic inheritance of few moderate-large effect loci and a large-unknown number of loci each having a small effect on BCWD resistance. We detected differences in QTL number and genome location between two GWAS models (weighted single-step GBLUP and Bayes B), which highlights the utility of using different models to uncover QTL. The RAD-SNPs detected a greater number of QTL than the 57 K SNP array in one population, suggesting that the RAD-SNPs may uncover polymorphisms that are more unique and informative for the specific population in which they were discovered.

Background: Much effort is put into identifying causative quantitative trait nucleotides (QTN) in animal breeding, empowered by the availability of dense single nucleotide polymorphism (SNP) information. Genomic selection using traditional SNP information is easily implemented for any number of genotyped individuals using single-step genomic best linear unbiased predictor (ssGBLUP) with the algorithm for proven and young (APY). Our aim was to investigate whether ssGBLUP is useful for genomic prediction when some or all QTN are known., Methods: Simulations included 180,000 animals across 11 generations. Phenotypes were available for all animals in generations 6 to 10. Genotypes for 60,000 SNPs across 10 chromosomes were available for 29,000 individuals. The genetic variance was fully accounted for by 100 or 1000 biallelic QTN. Raw genomic relationship matrices (GRM) were computed from (a) unweighted SNPs, (b) unweighted SNPs and causative QTN, (c) SNPs and causative QTN weighted with results obtained with genome-wide association studies, (d) unweighted SNPs and causative QTN with simulated weights, (e) only unweighted causative QTN, (f-h) as in (b-d) but using only the top 10% causative QTN, and (i) using only causative QTN with simulated weight. Predictions were computed by pedigree-based BLUP (PBLUP) and ssGBLUP. Raw GRM were blended with 1 or 5% of the numerator relationship matrix, or 1% of the identity matrix. Inverses of GRM were obtained directly or with APY., Results: Accuracy of breeding values for 5000 genotyped animals in the last generation with PBLUP was 0.32, and for ssGBLUP it increased to 0.49 with an unweighted GRM, 0.53 after adding unweighted QTN, 0.63 when QTN weights were estimated, and 0.89 when QTN weights were based on true effects known from the simulation. When the GRM was constructed from causative QTN only, accuracy was 0.95 and 0.99 with blending at 5 and 1%, respectively. Accuracies simulating 1000 QTN were generally lower, with a similar trend. Accuracies using the APY inverse were equal or higher than those with a regular inverse., Conclusions: Single-step GBLUP can account for causative QTN via a weighted GRM. Accuracy gains are maximum when variances of causative QTN are known and blending is at 1%.

Background: Previously, we have shown that bacterial cold water disease (BCWD) resistance in rainbow trout can be improved using traditional family-based selection, but progress has been limited to exploiting only between-family genetic variation. Genomic selection (GS) is a new alternative that enables exploitation of within-family genetic variation., Methods: We compared three GS models [single-step genomic best linear unbiased prediction (ssGBLUP), weighted ssGBLUP (wssGBLUP), and BayesB] to predict genomic-enabled breeding values (GEBV) for BCWD resistance in a commercial rainbow trout population, and compared the accuracy of GEBV to traditional estimates of breeding values (EBV) from a pedigree-based BLUP (P-BLUP) model. We also assessed the impact of sampling design on the accuracy of GEBV predictions. For these comparisons, we used BCWD survival phenotypes recorded on 7893 fish from 102 families, of which 1473 fish from 50 families had genotypes [57 K single nucleotide polymorphism (SNP) array]. Naïve siblings of the training fish (n = 930 testing fish) were genotyped to predict their GEBV and mated to produce 138 progeny testing families. In the following generation, 9968 progeny were phenotyped to empirically assess the accuracy of GEBV predictions made on their non-phenotyped parents., Results: The accuracy of GEBV from all tested GS models were substantially higher than the P-BLUP model EBV. The highest increase in accuracy relative to the P-BLUP model was achieved with BayesB (97.2 to 108.8%), followed by wssGBLUP at iteration 2 (94.4 to 97.1%) and 3 (88.9 to 91.2%) and ssGBLUP (83.3 to 85.3%). Reducing the training sample size to n = ~1000 had no negative impact on the accuracy (0.67 to 0.72), but with n = ~500 the accuracy dropped to 0.53 to 0.61 if the training and testing fish were full-sibs, and even substantially lower, to 0.22 to 0.25, when they were not full-sibs., Conclusions: Using progeny performance data, we showed that the accuracy of genomic predictions is substantially higher than estimates obtained from the traditional pedigree-based BLUP model for BCWD resistance. Overall, we found that using a much smaller training sample size compared to similar studies in livestock, GS can substantially improve the selection accuracy and genetic gains for this trait in a commercial rainbow trout breeding population.

Genetically linked small and large dairy cattle populations were simulated to test the effect of different sources of information from foreign populations on the accuracy of predicting breeding values for young animals in a small population. A large dairy cattle population (P L ) with >20 generations was simulated, and a small subpopulation (P S ) with 3 generations was formed as a related population, including phenotypes and genomic information. Predicted breeding values for young animals in the small population were calculated using BLUP and single-step genomic BLUP (ssGBLUP) in 4 different scenarios: (S1) 3,166 phenotypes, 22,855 pedigree animals, and 1,000 to 6,000 genotypes for P S ; (S2) S1 plus genomic estimated breeding value (GEBV) for 4,475 sires from P L as external information; (S3) S1 plus 221,580 phenotypes, 402,829 pedigree animals, and 53,558 genotypes for P L ; and (S4) single nucleotide polymorphism (SNP) effects calculated based on P L data. The ability to predict true breeding value was assessed in the youngest third of the genotyped animals in the small population. When data only from the small population were used and 1,000 animals were genotyped, the accuracy of GEBV was only 1 point greater than the estimated breeding value accuracy (0.32 vs. 0.31). Adding external GEBV for sires from P L did not considerably increase accuracy (0.33 vs. 0.32 in S1). Combining phenotypes, pedigree, and genotypes for P S and P L was beneficial for predicting accuracy of GEBV in the small population, and the prediction accuracy of GEBV in this scenario was 0.38 compared with 0.31 from estimated breeding values. When SNP effects from P L were used to predict GEBV for young genotyped animals from P S , accuracy was greatest (0.56). With 6,000 genotyped animal in P S , accuracy was greatest (0.61) with the combined populations. In a small population with few genotypes, the highest accuracy of evaluation may be obtained by using SNP effects derived from a related large population., (Copyright © 2017 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

This paper evaluates an efficient implementation to multiply the inverse of a numerator relationship matrix for genotyped animals () by a vector (). The computation is required for solving mixed model equations in single-step genomic BLUP (ssGBLUP) with the preconditioned conjugate gradient (PCG). The inverse can be decomposed into sparse matrices that are blocks of the sparse inverse of a numerator relationship matrix () including genotyped animals and their ancestors. The elements of were rapidly calculated with the Henderson's rule and stored as sparse matrices in memory. Implementation of was by a series of sparse matrix-vector multiplications. Diagonal elements of , which were required as preconditioners in PCG, were approximated with a Monte Carlo method using 1,000 samples. The efficient implementation of was compared with explicit inversion of with 3 data sets including about 15,000, 81,000, and 570,000 genotyped animals selected from populations with 213,000, 8.2 million, and 10.7 million pedigree animals, respectively. The explicit inversion required 1.8 GB, 49 GB, and 2,415 GB (estimated) of memory, respectively, and 42 s, 56 min, and 13.5 d (estimated), respectively, for the computations. The efficient implementation required <1 MB, 2.9 GB, and 2.3 GB of memory, respectively, and <1 sec, 3 min, and 5 min, respectively, for setting up. Only <1 sec was required for the multiplication in each PCG iteration for any data sets. When the equations in ssGBLUP are solved with the PCG algorithm, is no longer a limiting factor in the computations.

The purposes of this study were to analyze the impact of seasonal losses due to heat stress in pigs from different breeds raised in different environments and to evaluate the accuracy improvement from adding genomic information to genetic evaluations. Data were available for 2 different swine populations: purebred Duroc animals raised in Texas and North Carolina and commercial crosses of Duroc and F females (Landrace × Large White) raised in Missouri and North Carolina; pedigrees provided links for animals from different states. Pedigree information was available for 553,442 animals, of which 8,232 pure breeds were genotyped. Traits were BW at 170 d for purebred animals and HCW for crossbred animals. Analyses were done with an animal model as either single- or 2-trait models using phenotypes measured in different states as separate traits. Additionally, reaction norm models were fitted for 1 or 2 traits using heat load index as a covariable. Heat load was calculated as temperature-humidity index greater than 70 and was averaged over 30 d prior to data collection. Variance components were estimated with average information REML, and EBV and genomic EBV (GEBV) with BLUP or single-step genomic BLUP (ssGBLUP). Validation was assessed for 146 genotyped sires with progeny in the last generation. Accuracy was calculated as a correlation between EBV and GEBV using reduced data (all animals, except the last generation) and using complete data. Heritability estimates for purebred animals were similar across states (varying from 0.23 to 0.26), and reaction norm models did not show evidence of a heat stress effect. Genetic correlations between states for heat loads were always strong (>0.91). For crossbred animals, no differences in heritability were found in single- or 2-trait analysis (from 0.17 to 0.18), and genetic correlations between states were moderate (0.43). In the reaction norm for crossbreeds, heritabilities ranged from 0.15 to 0.30 and genetic correlations between heat loads were as weak as 0.36, with heat load ranging from 0 to 12. Accuracies with ssGBLUP were, on average, 25% greater than with BLUP. Accuracies were greater in 2-trait reaction norm models and at extreme heat load values. Impacts of seasonality are evident only for crossbred animals. Genomic information can help producers mitigate heat stress in swine by identifying superior sires that are more resistant to heat stress.

The purpose of this study was to analyze the impact of seasonal losses due to heat stress in different environments and genetic group combinations. Data were available for 2 different swine populations: purebred Duroc animals raised in nucleus farms in Texas and North Carolina and crosses of Duroc and F females (Landrace × Large White) raised in commercial farms in Missouri and North Carolina; pedigrees provided links between animals from different states. Traits included BW at harvest age for purebred animals and HCW for crossbred animals. Weather data were collected at airports located close to the farms. Heat stress was quantified by a heat load function, defined by the units of temperature-humidity of temperature-humidity index (THI) greater than a certain threshold for 30 to 70 d before phenotype collection. Heat stress responses were quantified by a linear regression of phenotype on heat load. The greatest coefficient of determination occurred with a length of 30 d before phenotype measurements for all states and genetic groups. In the crossbreed data, THI thresholds were 67 in Missouri and 72 in North Carolina. For pure breeds, heat load had the best fit for THI thresholds greater than 70 in North Carolina, although differences in coefficient of determinations were negligible. On the other hand, no optimal THI threshold existed in Texas. In this study, heat stress had a greater impact in commercial farms than in nucleus farms and the effect of heat stress on weight varied by year and state.

This study evaluated the impact of region and season on growth in Angus seed stock. To assess geographic differences, the United States was partitioned into 9 regions based on similar climate and topography related to cow-calf production. Seasonal effects were associated with the month that animals were weighed. The American Angus Association provided growth data, and records were assigned to regions based on the owner's zip code. Most Angus cattle were in the Cornbelt, Lower Plains, Rocky Mountain, Upper Plains, and Upper South regions, with proportionally fewer Angus in Texas compared with the national cow herd. Most calves were born in the spring, especially February and March. Weaning weights (WW; = 49,886) and yearling weights (YW; = 45,168) were modeled with fixed effects of age-of-dam class (WW only), weigh month, region, month-region interaction, and linear covariate of age. Random effects included contemporary group nested within month-region combination and residual. The significant month-region interaction ( < 0.0001) was expected because of the diverse production environments across the country and cyclical fluctuations in forage availability. Additionally, significant seasonal contrasts existed for several regions. Fall-born calves were heavier ( < 0.01) than spring-born calves in the hot and humid Lower South region coinciding with fall being the primary calving season. The North and Upper Plains regions had heavier, spring-born calves ( < 0.01), more than 90% spring calving, and colder climates. Interestingly, no seasonal WW or YW differences existed between spring- and fall-born calves in the upper South region despite challenging environmental conditions. Angus seed stock producers have used calving seasons to adapt to the specific environmental conditions in their regions and to optimize growth in young animals.

The objectives were to assess the impact of heat stress and to develop a model for genetic evaluation of growth heat tolerance in Angus cattle. The American Angus Association provided weaning weight (WW) and yearling weight (YW) data, and records from the Upper South region were used because of the hot climatic conditions. Heat stress was characterized by a weaning (yearling) heat load function defined as the mean temperature-humidity index (THI) units greater than 75 (70) for 30 (150) d prior to the weigh date. Therefore, a weaning (yearling) heat load of 5 units corresponded to 80 (75) for the corresponding period prior to the weigh date. For all analyses, 82,669 WW and 69,040 YW were used with 3 ancestral generations in the pedigree. Univariate models were a proxy for the Angus growth evaluation, and reaction norms using 2 B-splines for heat load were fit separately for weaning and yearling heat loads. For both models, random effects included direct genetic, maternal genetic, maternal permanent environment (WW only), and residual. Fixed effects included a linear age covariate, age-of-dam class (WW only), and contemporary group for both models and fixed regressions on the B-splines in the reaction norm. Direct genetic correlations for WW were strong for modest heat load differences but decreased to less than 0.50 for large differences. Reranking of proven sires occurred for only WW direct effects for the reaction norms with extreme heat load differences. Conversely, YW results indicated little effect of heat stress on genetic merit. Therefore, weaning heat tolerance was a better candidate for developing selection tools. Maternal heritabilities were consistent across heat loads, and maternal genetic correlations were greater than 0.90 for nearly all heat load combinations. No evidence existed for a genotype × environment interaction for the maternal component of growth. Overall, some evidence exists for phenotypic plasticity for the direct genetic effects of WW, but traditional national cattle evaluations are likely adequately ranking sires for nonextreme environmental conditions.

Animal feeding is the most important economic component of beef production systems. Selection for feed efficiency has not been effective mainly due to difficult and high costs to obtain the phenotypes. The application of genomic selection using SNP can decrease the cost of animal evaluation as well as the generation interval. The objective of this study was to compare methods for genomic evaluation of feed efficiency traits using different cross-validation layouts in an experimental beef cattle population genotyped for a high-density SNP panel (BovineHD BeadChip assay 700k, Illumina Inc., San Diego, CA). After quality control, a total of 437,197 SNP genotypes were available for 761 Nelore animals from the Institute of Animal Science, Sertãozinho, São Paulo, Brazil. The studied traits were residual feed intake, feed conversion ratio, ADG, and DMI. Methods of analysis were traditional BLUP, single-step genomic BLUP (ssGBLUP), genomic BLUP (GBLUP), and a Bayesian regression method (BayesCπ). Direct genomic values (DGV) from the last 2 methods were compared directly or in an index that combines DGV with parent average. Three cross-validation approaches were used to validate the models: 1) YOUNG, in which the partition into training and testing sets was based on year of birth and testing animals were born after 2010; 2) UNREL, in which the data set was split into 3 less related subsets and the validation was done in each subset a time; and 3) RANDOM, in which the data set was randomly divided into 4 subsets (considering the contemporary groups) and the validation was done in each subset at a time. On average, the RANDOM design provided the most accurate predictions. Average accuracies ranged from 0.10 to 0.58 using BLUP, from 0.09 to 0.48 using GBLUP, from 0.06 to 0.49 using BayesCπ, and from 0.22 to 0.49 using ssGBLUP. The most accurate and consistent predictions were obtained using ssGBLUP for all analyzed traits. The ssGBLUP seems to be more suitable to obtain genomic predictions for feed efficiency traits on an experimental population of genotyped animals.

Bacterial cold water disease (BCWD) causes significant economic losses in salmonid aquaculture, and traditional family-based breeding programs aimed at improving BCWD resistance have been limited to exploiting only between-family variation. We used genomic selection (GS) models to predict genomic breeding values (GEBVs) for BCWD resistance in 10 families from the first generation of the NCCCWA BCWD resistance breeding line, compared the predictive ability (PA) of GEBVs to pedigree-based estimated breeding values (EBVs), and compared the impact of two SNP genotyping methods on the accuracy of GEBV predictions. The BCWD phenotypes survival days (DAYS) and survival status (STATUS) had been recorded in training fish (n = 583) subjected to experimental BCWD challenge. Training fish, and their full sibs without phenotypic data that were used as parents of the subsequent generation, were genotyped using two methods: restriction-site associated DNA (RAD) sequencing and the Rainbow Trout Axiom® 57 K SNP array (Chip). Animal-specific GEBVs were estimated using four GS models: BayesB, BayesC, single-step GBLUP (ssGBLUP), and weighted ssGBLUP (wssGBLUP). Family-specific EBVs were estimated using pedigree and phenotype data in the training fish only. The PA of EBVs and GEBVs was assessed by correlating mean progeny phenotype (MPP) with mid-parent EBV (family-specific) or GEBV (animal-specific). The best GEBV predictions were similar to EBV with PA values of 0.49 and 0.46 vs. 0.50 and 0.41 for DAYS and STATUS, respectively. Among the GEBV prediction methods, ssGBLUP consistently had the highest PA. The RAD genotyping platform had GEBVs with similar PA to those of GEBVs from the Chip platform. The PA of ssGBLUP and wssGBLUP methods was higher with the Chip, but for BayesB and BayesC methods it was higher with the RAD platform. The overall GEBV accuracy in this study was low to moderate, likely due to the small training sample used. This study explored the potential of GS for improving resistance to BCWD in rainbow trout using, for the first time, progeny testing data to assess the accuracy of GEBVs, and it provides the basis for further investigation on the implementation of GS in commercial rainbow trout populations.

The objectives of this study were to develop and evaluate an efficient implementation in the computation of the inverse of genomic relationship matrix with the recursion algorithm, called the algorithm for proven and young (APY), in single-step genomic BLUP. We validated genomic predictions for young bulls with more than 500,000 genotyped animals in final score for US Holsteins. Phenotypic data included 11,626,576 final scores on 7,093,380 US Holstein cows, and genotypes were available for 569,404 animals. Daughter deviations for young bulls with no classified daughters in 2009, but at least 30 classified daughters in 2014 were computed using all the phenotypic data. Genomic predictions for the same bulls were calculated with single-step genomic BLUP using phenotypes up to 2009. We calculated the inverse of the genomic relationship matrix GAPY(-1) based on a direct inversion of genomic relationship matrix on a small subset of genotyped animals (core animals) and extended that information to noncore animals by recursion. We tested several sets of core animals including 9,406 bulls with at least 1 classified daughter, 9,406 bulls and 1,052 classified dams of bulls, 9,406 bulls and 7,422 classified cows, and random samples of 5,000 to 30,000 animals. Validation reliability was assessed by the coefficient of determination from regression of daughter deviation on genomic predictions for the predicted young bulls. The reliabilities were 0.39 with 5,000 randomly chosen core animals, 0.45 with the 9,406 bulls, and 7,422 cows as core animals, and 0.44 with the remaining sets. With phenotypes truncated in 2009 and the preconditioned conjugate gradient to solve mixed model equations, the number of rounds to convergence for core animals defined by bulls was 1,343; defined by bulls and cows, 2,066; and defined by 10,000 random animals, at most 1,629. With complete phenotype data, the number of rounds decreased to 858, 1,299, and at most 1,092, respectively. Setting up GAPY(-1) for 569,404 genotyped animals with 10,000 core animals took 1.3h and 57 GB of memory. The validation reliability with APY reaches a plateau when the number of core animals is at least 10,000. Predictions with APY have little differences in reliability among definitions of core animals. Single-step genomic BLUP with APY is applicable to millions of genotyped animals., (Copyright © 2016 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Combining purebreed and crossbreed information is beneficial for genetic evaluation of some livestock species. Genetic evaluations can use relationships based on genomic information, relying on allele frequencies that are breed specific. Single-step genomic BLUP (ssGBLUP) does not account for different allele frequencies, which could limit the genetic gain in crossbreed evaluations. In this study, we tested the performance of different breed-specific genomic relationship matrices () in ssGBLUP for crossbreed evaluations; we also tested the importance of genotyping crossbred animals. Genotypes were available for purebreeds (AA and BB) and crossbreeds (F) in simulated and real pig populations. The number of genotyped animals was, on average, 4,315 for the simulated population and 15,798 for the real population. Cross-validation was performed on 1,200 and 3,117 F animals in the simulated and real populations, respectively. Simulated scenarios were under no artificial selection, mass selection, or BLUP selection. Two genomic relationship matrices were constructed based on breed-specific allele frequencies: 1) , a genomic relationship matrix centered by breed-specific allele frequencies, and 2) , a genomic relationship matrix centered and scaled by breed-specific allele frequencies. All (the across-breed genomic relationship matrix), , and were also tuned to account for selective genotyping. Using breed-specific allele frequencies reduced the number of negative relationships between 2 purebreeds, pulling the average closer to 0, as in the pedigree-based relationship matrix. For simulated populations that included mass selection, genomic EBV (GEBV) in F, when using and , were, on average, 10% more accurate than ; however, after tuning to account for selective genotyping, provided the same accuracy as for breed-specific genomic relationship matrices. For the real population, accuracies for litter size in F were 0.62 for , , and , and tuning had no impact on accuracy, except for , which was 1 percentage point less accurate. Accuracy of GEBV for number of stillborns in F1 was 0.5 for all tested genomic relationship matrices with no changes after tuning. We observed that genotyping F increased accuracies of GEBV for the same animals by up to 39% compared with having genotypes for only AA and BB. In crossbreed evaluations, accounting for breed-specific allele frequencies promoted changes in G that were not influential enough to improve accuracy of GEBV. Therefore, the best performance of ssGBLUP for crossbreed evaluations requires genotypes for pure- and crossbreeds and no breed-specific adjustments in the realized relationship matrix.

The purpose of this study was to examine accuracy of genomic selection via single-step genomic BLUP (ssGBLUP) when the direct inverse of the genomic relationship matrix (G) is replaced by an approximation of G(-1) based on recursions for young genotyped animals conditioned on a subset of proven animals, termed algorithm for proven and young animals (APY). With the efficient implementation, this algorithm has a cubic cost with proven animals and linear with young animals. Ten duplicate data sets mimicking a dairy cattle population were simulated. In a first scenario, genomic information for 20k genotyped bulls, divided in 7k proven and 13k young bulls, was generated for each replicate. In a second scenario, 5k genotyped cows with phenotypes were included in the analysis as young animals. Accuracies (average for the 10 replicates) in regular EBV were 0.72 and 0.34 for proven and young animals, respectively. When genomic information was included, they increased to 0.75 and 0.50. No differences between genomic EBV (GEBV) obtained with the regular G(-1) and the approximated G(-1) via the recursive method were observed. In the second scenario, accuracies in GEBV (0.76, 0.51 and 0.59 for proven bulls, young males and young females, respectively) were also higher than those in EBV (0.72, 0.35 and 0.49). Again, no differences between GEBV with regular G(-1) and with recursions were observed. With the recursive algorithm, the number of iterations to achieve convergence was reduced from 227 to 206 in the first scenario and from 232 to 209 in the second scenario. Cows can be treated as young animals in APY without reducing the accuracy. The proposed algorithm can be implemented to reduce computing costs and to overcome current limitations on the number of genotyped animals in the ssGBLUP method., (© 2015 Blackwell Verlag GmbH.)

Background: As more and more genotypes become available, accuracy of genomic evaluations can potentially increase. However, the impact of genotype data on accuracy depends on the structure of the genotyped cohort. For populations such as dairy cattle, the greatest benefit has come from genotyping sires with high accuracy, whereas the benefit due to adding genotypes from cows was smaller. In broiler chicken breeding programs, males have less progeny than dairy bulls, females have more progeny than dairy cows, and most production traits are recorded for both sexes. Consequently, genotyping both sexes in broiler chickens may be more advantageous than in dairy cattle., Methods: We studied the contribution of genotypes from males and females using a real dataset with genotypes on 15 723 broiler chickens. Genomic evaluations used three training sets that included only males (4648), only females (8100), and both sexes (12 748). Realized accuracies of genomic estimated breeding values (GEBV) were used to evaluate the benefit of including genotypes for different training populations on genomic predictions of young genotyped chickens., Results: Using genotypes on males, the average increase in accuracy of GEBV over pedigree-based EBV for males and females was 12 and 1 percentage points, respectively. Using female genotypes, this increase was 1 and 18 percentage points, respectively. Using genotypes of both sexes increased accuracies by 19 points for males and 20 points for females. For two traits with similar heritabilities and amounts of information, realized accuracies from cross-validation were lower for the trait that was under strong selection., Conclusions: Overall, genotyping males and females improves predictions of all young genotyped chickens, regardless of sex. Therefore, when males and females both contribute to genetic progress of the population, genotyping both sexes may be the best option.

The purpose of this study was to evaluate the accuracy of genomic selection in single-step genomic BLUP (ssGBLUP) when the inverse of the genomic relationship matrix (G) is derived by the "algorithm for proven and young animals" (APY). This algorithm implements genomic recursions on a subset of "proven" animals. Only a relationship matrix for animals treated as "proven" needs to be inverted, and the extra costs of adding animals treated as "young" are linear. Analyses involved 10,102,702 final scores on 6,930,618 Holstein cows. Final score, which is a composite of type traits, is popular trait in the United States and was easily available for this study. A total of 100,000 animals with genotypes were used in the analyses and included 23,000 sires (16,000 with >5 progeny), 27,000 cows, and 50,000 young animals. Genomic EBV (GEBV) were calculated with a regular inverse of G, and with the G inverse approximated by APY. Animals in the proven subset included only sires (23,000), sires+cows (50,000), only cows (27,000), or sires with >5 progeny (16,000). The correlations of GEBV with APY and regular GEBV for young genotyped animals were 0.994, 0.995, 0.992, and 0.992, respectively Later, animals in the proven subset were randomly sampled from all genotyped animals in sets of 2,000, 5,000, 10,000, 15,000, and 20,000; each sample was replicated 4 times. Respective correlations were 0.97 (5,000 sample), 0.98 (10,000 sample), and 0.99 (20,000 sample), with minimal difference between samples of the same size. Genomic EBV with APY were accurate when the number of animals used in the subset is between 10,000 and 20,000, with little difference between the ways of creating the subset. Due to the approximately linear cost of APY, ssGBLUP with APY could support any number of genotyped animals without affecting accuracy., (Copyright © 2015 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Predictive ability of genomic EBV when using single-step genomic BLUP (ssGBLUP) in Angus cattle was investigated. Over 6 million records were available on birth weight (BiW) and weaning weight (WW), almost 3.4 million on postweaning gain (PWG), and over 1.3 million on calving ease (CE). Genomic information was available on, at most, 51,883 animals, which included high and low EBV accuracy animals. Traditional EBV was computed by BLUP and genomic EBV by ssGBLUP and indirect prediction based on SNP effects was derived from ssGBLUP; SNP effects were calculated based on the following reference populations: ref&#95;2k (contains top bulls and top cows that had an EBV accuracy for BiW ≥0.85), ref&#95;8k (contains all parents that were genotyped), and ref&#95;33k (contains all genotyped animals born up to 2012). Indirect prediction was obtained as direct genomic value (DGV) or as an index of DGV and parent average (PA). Additionally, runs with ssGBLUP used the inverse of the genomic relationship matrix calculated by an algorithm for proven and young animals (APY) that uses recursions on a small subset of reference animals. An extra reference subset included 3,872 genotyped parents of genotyped animals (ref&#95;4k). Cross-validation was used to assess predictive ability on a validation population of 18,721 animals born in 2013. Computations for growth traits used multiple-trait linear model and, for CE, a bivariate CE-BiW threshold-linear model. With BLUP, predictivities were 0.29, 0.34, 0.23, and 0.12 for BiW, WW, PWG, and CE, respectively. With ssGBLUP and ref&#95;2k, predictivities were 0.34, 0.35, 0.27, and 0.13 for BiW, WW, PWG, and CE, respectively, and with ssGBLUP and ref&#95;33k, predictivities were 0.39, 0.38, 0.29, and 0.13 for BiW, WW, PWG, and CE, respectively. Low predictivity for CE was due to low incidence rate of difficult calving. Indirect predictions with ref&#95;33k were as accurate as with full ssGBLUP. Using the APY and recursions on ref&#95;4k gave 88% gains of full ssGBLUP and using the APY and recursions on ref&#95;8k gave 97% gains of full ssGBLUP. Genomic evaluation in beef cattle with ssGBLUP is feasible while keeping the models (maternal, multiple trait, and threshold) already used in regular BLUP. Gains in predictivity are dependent on the composition of the reference population. Indirect predictions via SNP effects derived from ssGBLUP allow for accurate genomic predictions on young animals, with no advantage of including PA in the index if the reference population is large. With the APY conditioning on about 10,000 reference animals, ssGBLUP is potentially applicable to a large number of genotyped animals without compromising predictive ability.

The study reported here evaluated genotype × environment interaction in individual performance and progeny tests in beef cattle. Genetic parameters for final weight (FW), ADG, and scrotal circumference (SC) of 33,013 Nellore young bulls tested on pasture or in feedlots were analyzed. The posterior means (and highest posterior density interval with 90% of samples [HPD90]) of heritability for traits measured on pasture-raised and feedlot-raised animals were 0.44 (HPD90 = 0.40 to 0.48) and 0.50 (HPD90 = 0.43 to 0.56) for FW, 0.26 (HPD90 = 0.23 to 0.29) and 0.26 (HPD90 = 0.20 to 0.32) for ADG, and 0.53 (HPD90 = 0.48 to 0.59) and 0.65 (HPD90 = 0.55 to 0.74) for SC, respectively. The posterior means (and HPD90) of genetic correlations for FW, ADG, and SC on pasture and in feedlots were 0.75 (HPD90 = 0.66 to 0.87), 0.49 (HPD90 = 0.31 to 0.66), and 0.89 (HPD90 = 0.83 to 0.97), respectively. When the selection intensity was kept the same for both the environments, the greatest direct responses for FW and ADG were exhibited by the animals reared and selected in feedlots. The correlated responses relative to production on pasture and based on selection in feedlots were similar to the direct responses, whereas the correlated responses for production in feedlots and based on selection on pasture were lower than the direct responses. When the selection intensity on pasture was higher than the selection intensity in feedlots, the responses to direct selection were similar for both the environments and correlated responses obtained in feedlots by selection on pasture were similar to the direct responses in feedlots. Analyses of few or poor indicators of genotype × environment interaction result in incorrect interpretations of its existence and implications. The present work demonstrated that traits with lower heritability are more susceptible to genotype × environment interaction and that selection intensity plays an important role in the study of genotype × environment interaction in beef cattle.

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.

Studies have shown cases of poisoning with plants from the genus Crotalaria (Leguminosae) mainly in animals. They induce damages in the central nervous system (CNS), which has been attributed to toxic effects of the pyrrolizidine alkaloid (PA) monocrotaline (MCT). Previously we demonstrated that both MCT and dehydromonocrotaline (DHMC), its main active metabolite, induce changes in the levels and patterns of expression of the main protein from astrocyte cytoskeleton, glial fibrillary acidic protein (GFAP). In this study we investigated the effect of MCT on rat cortical astrocyte/neuron primary co-cultures. Primary cultures were exposed to 10 or 100 μM MCT. The MTT test and the measurement of LDH activity on the culture medium revealed that after 24h exposure MCT was not cytotoxic to neuron/astrocyte cells. However, the cell viability after 72 h treatment decreased in 10-20%, and the LDH levels in the culture medium increased at a rate of 12% and 23%, in cultures exposed to 10 or 100 μM MCT. Rosenfeld staining showed vacuolization and increase in cell body in astrocytes after MCT exposure. Immunocytochemistry and Western blot analyses revealed changes on pattern of GFAP and βIII-tubulin expression and steady state levels after MCT treatment, with a dose and time dependent intense down regulation and depolarization of neuronal βIII-tubulin. Moreover, treatment with 100 μM MCT for 12h induced GSH depletion, which was not seen when cytochrome P450 enzyme system was inhibited indicating that it is involved in MCT induced cytotoxicity in CNS cells., (Copyright © 2011 Elsevier Inc. All rights reserved.)

In this study, we investigated the effects of the flavonoid rutin (3,3',4',5,7-pentahydroxyflavone-3-rutinoside) on glioma cells, using the highly proliferative human cell line GL-15 as a model. We observed that rutin (50-100μM) reduced proliferation and viability of GL-15 cells, leading to decreased levels of ERK1/2 phosphorylation (P-ERK1/2) and accumulation of cells in the G2 phase of the cell cycle. On the other hand, 87.4% of GL-15 cells exposed to 100μM rutin entered apoptosis, as revealed by flow cytometry after AnnexinV/PI staining. Nuclear condensation and DNA fragmentation were also observed, further confirming that apoptosis had occurred. Moreover, the remaining cells that were treated with 50μM rutin presented a morphological pattern of astroglial differentiation in culture, characterised by a condensed cell body and thin processes with overexpression of GFAP. Because of its capacity to induce differentiation and apoptosis in cultured human glioblastoma cells, rutin could be considered as a potential candidate for malignant gliomas treatment., (Copyright © 2011 Elsevier Ltd. All rights reserved.)