Weighted Single-step GBLUP Research Articles

Weighed ssGBLUP enhances the genomic prediction accuracy for milk citrate predicted by milk mid-infrared spectra of Holstein cows in early lactation

Previous studies have shown that milk citrate predicted by milk mid-infrared spectra (MIR) is strongly affected by a few genomic regions. This study aimed to explore the effect of the weighted single-step GBLUP on the accuracy of genomic prediction (GP) for MIR-predicted milk citrate in early lactation Holstein cows. A total of 134,517 test-day predicted milk citrate collected within the first 50 d in milk on 52,198 Holstein cows from the first 5 parities were used. There were 122,218 animals in the pedigree, of them 4,479 had genotypic data for 566,170 SNPs. Two data sets (partial and whole data sets) were used to verify whether the accuracy of GP is improved using the following different methods. The (genomic) estimated breeding values (EBV or GEBV) in the partial and whole data sets were estimated by pedigree-based BLUP (ABLUP), single-step GBLUP (ssGBLUP, pedigree-genomic combined using no weight for SNP), weighted ssGBLUP (WssGBLUP, pedigree-genomic combined using weighted SNP), respectively. The difference between the 2 data sets is that the phenotypic data from 2017 to 2019 in the partial data set were set as missing values. 181 youngest cows with genomic were selected as the validation population. A linear regression method was used to compare EBV (GEBV) predicted for partial and whole data sets. The accuracies of GP for ABLUP and ssGBLUP were 0.42 and 0.70, respectively. The accuracies of GP for WssGBLUP in the 5 iterations with different CT (constant) values (determines departure from normality for SNP effects) ranged from 0.70 to 0.86. This study showed that weighted SNP is beneficial in improving prediction accuracy for predicted milk citrate.

Genome-wide scans identify biological and metabolic pathways regulating carcass and meat quality traits in beef cattle

Genome association studies (GWAS) provides knowledge about the genetic architecture of beef-related traits that allow linking the target phenotype to genomic information aiding breeding decision. Thus, the present study aims to uncover the genetic mechanism involved in carcass (REA: rib eye area, BF: backfat thickness, and HCW: hot carcass weight) and meat quality traits (SF: shear-force, MARB: marbling score, and IMF: intramuscular fat content) in Nellore cattle. For this, 6910 young bulls with phenotypic information and 23,859 animals genotyped with 435 k markers were used to perform the weighted single-step GBLUP (WssGBLUP) approach, considering two iterations. The top 10 genomic regions explained 8.13, 11.81, and 9.58% of the additive genetic variance, harboring a total of 119, 143, and 95 positional candidate genes for REA, BF, and HCW, respectively. For meat quality traits, the top 10 windows explained a large proportion of the total genetic variance for SF (14.95%), MARB (17.56%), and IMF (21.41%) surrounding 92, 155, and 111 candidate genes, respectively. Relevant candidate genes (CAST, PLAG1, XKR4, PLAGL2, AQP3/AQP7, MYLK2, WWOX, CARTPT, and PLA2G16) are related to physiological aspects affecting growth, carcass, meat quality, feed intake, and reproductive traits by signaling pathways controlling muscle control, key signal metabolic molecules INS / IGF-1 pathway, lipid metabolism, and adipose tissue development. The GWAS results provided insights into the genetic control of the traits studied and the genes found are potential candidates to be used in the improvement of carcass and meat quality traits.

141 Exploring Alternative Approaches for the National Argentine Angus Cattle Genomic Evaluation Program

Abstract Beef cattle production is a relevant component of the Argentine economy. The Argentine Angus is the most common beef cattle breed, and this population has a large economic and cultural importance for the country. In 2019, a national genomic selection program was implemented in the Argentine Angus population. The objective of this study was to explore alternative approaches for the genomic selection program in Argentinean Angus so that it will better serve the stakeholders. For this study, a total of 9,944 phenotypes for weaning weight from Argentine Angus born in 2018 were used. All animals had pedigree information for at least 3 generations. Genomic information on 1,078 individuals for 150k SNPs was used. The genotyped individuals were connected to the population by pedigree. Weaning weight EPDs were calculated under an animal model considering contemporary group as fixed effects and additive direct and maternal genetic effects as random. Variance components were estimated using AIREMLF90, and EPDs were calculated using BLUPF90. The enhanced EPDs were the genomic estimated breeding values obtained by single-step GBLUP, by including the genomic information into the original model. Validation of the EPDs and enhanced EPDs were made by cross-validation, where a random sample of 70 genotyped individuals with more than 20 progeny was the validation population. The reference breeding values were the EPDs of the validation population with complete data. In a second analysis, the phenotype of the progeny of the validation population was removed. The correlation between the reduced and complete EPDs was the accuracy for this study, according to the LR method. The accuracy obtained by this method was overestimated due to potential relationships between animals in the training and validation population. However, the enhanced EPDs were 4% more accurate than the pedigree EPDs. The next steps of this research are exploring alternative approaches of the validation methodology, which will be accomplished by implementing a forward-in-time cross-validation. Such validation will be performed once more data is included in the model. Additionally, other models will be tested to identify the ideal combination of effects that will yield higher accuracies. Details regarding the construction of the genomic relationship matrix will be tested with the overarching goal of improving convergence without affecting accuracy. Moreover, a weighted single-step GBLUP approach will be tested to verify if prediction accuracy can improve even further. Finally, indirect prediction models will be tested so that the Argentine Angus Association can implement interim enhanced EPDs for young animals. In conclusion, including genomic data into this dataset yielded greater accuracy and such values are expected to improve even further after optimization.

Genomic association using principal components of morphometric traits in horses: identification of genes related to bone growth

The measurement of morphometric traits in horses is important for determining breed qualification and is one of the main selection criteria for the species. The development of an index (HPC) that consists of principal components weighted by additive genetic values allows to explore the most relevant relationships using a reduced number of variables that explain the greatest amount of variation in the data. Genome-wide association studies (GWAS) using HPC are a relatively new approach that permits to identify regions related to a set of traits. The aim of this study was to perform GWAS using HPC for 15 linear measurements as the explanatory variable in order to identify associated genomic regions and to elucidate the biological mechanisms linked to this index in Campolina horses. For GWAS, weighted single-step GBLUP was applied to HPC. The eight genomic windows that explained the highest proportion of additive genetic variance were identified. The sum of the additive variance explained by the eight windows was 95.89%. Genes involved in bone and cartilage development were identified (SPRY2, COL9A2, MIR30C, HEYL, BMP8B, LTBP1, FAM98A, and CRIM1). They represent potential positional candidates for the HPC of the linear measurements evaluated. The HPC is an efficient alternative to reduce the 15 usually measured traits in Campolina horses. Moreover, candidate genes inserted in region that explained high additive variance of the HPC were identified and might be fine-mapped for searching putative mutation/markers.

17 Genome-wide Association Study for Heat Stress Tolerance in Swine Using Weighted Single-step GBLUP

Abstract The purpose of this study was to perform a genome-wide association study to determine the genomic regions associated with heat stress tolerance in swine as well as analyze the accuracy of prediction. Phenotypic information on carcass weight was available for 227,043 individuals from commercial farms in North Carolina and Missouri. Individuals were a commercial cross of a Duroc sire and a dam resulting from a Landrace and Large White cross. Genotypic information was available for 8,232 animals with 33,581 SNP. The pedigree file contained 553,448 animals. A 78 on the Temperature Humidity Index (THI) was used as a threshold for heat stress. A two-trait analysis was used with the phenotypes heat stress (trait one) and non-heat stress (trait two). Variance components were calculated via AIREML and breeding values were calculated using single step GBLUP (ssGBLUP). The heritability for trait one and two were calculated at 0.25 and 0.20, respectively, and the genetic correlation was calculated as 0.63. Validation was calculated for 163 genotyped sires with progeny in the last generation. The GEBV of complete data was used as the benchmark, and the accuracy was determined as the correlation between the GEBV of the reduced and complete data for the validation sires. Weighted ssGBLUP did not increase the accuracies, both methods showed a maximum accuracy of 0.32 for trait one and 0.54 for trait two. Manhattan Plots for trait one, trait two, and the difference between the two were created from the results of the two-trait analysis. Windows explaining around 1% of the genetic variance were identified. The only difference between the two traits was a peak at chromosome 14. The genetic correlation suggests different mechanisms for growth under heat stress. The GWAS results show that both traits are highly polygenic, with few genomic regions explaining more than 1% of variance.

The accuracy of genomic predictions for bacterial cold water disease resistance remains higher than the pedigree-based model one generation after model training in a commercial rainbow trout breeding population

Bacterial cold water disease (BCWD) causes significant mortality and economic losses in salmonid aquaculture. Previously, we reported high genomic prediction accuracy of 0.72 for BCWD resistance in rainbow trout using the 57 K Axiom SNP array. Due to the high cost of phenotyping and genotyping in rainbow trout breeding programs, it is paramount to know if acceptable genomic prediction accuracy can be obtained in the subsequent generation without retraining the prediction model. In the current study, we found that the accuracy of genomic prediction without model retraining in the subsequent generation was reduced from 0.65 and 0.61 to 0.56 and 0.53 using a lower density array of 10 K SNPs and a panel of 49 QTL-linked SNPs, respectively. Although markedly lower than the genomic prediction with model retraining, the prediction accuracy without retraining was better than the pedigree-based model (PBLUP) with retraining (0.48) and substantially higher than the PBLUP model accuracy without retraining (0.22). We conclude that genomic selection provides better prediction accuracy than the traditional PBLUP model even when ‘skipping’ one generation of model retraining. The weighted single-step GBLUP (wssGBLUP) and single-step Bayesian multiple regression BayesB (ssBMR-BayesB) had higher genomic prediction accuracy than single-step GBLUP (ssGBLUP), which is consistent with the oligogenic inheritance of BCWD resistance in this population. Imputation from the 10 K array genotypes back to the 57 K array genotypes did not improve the accuracy of genomic prediction, likely due to the high linkage disequilibrium in rainbow trout aquaculture breeding populations and the oligogenic architecture of BCWD resistance.

Genomic Prediction in Local Breeds: The Rendena Cattle as a Case Study.

Simple SummaryAlthough genomic selection is being used in many livestock species, it has not yet been considered in local breeds due to the lower population size and the potential less effective impact on the genetic evaluation of these breeds. The current research aims to investigate how genomic data can impact the accuracy of genetic predictions for beef traits in Rendena, a small local cattle breed of the North-East of Italy selected for a dual purpose. Classical animal models using only phenotypic information were compared with two models that integrated genomic data with pedigree information. The genomic models presented better accuracy in estimated breeding values of the animals than the ‘classical’ animal model, especially the ‘simpler’ one assuming homogeneous variances of single nucleotide polymorphisms. Our results show that the inclusion of genomic information can be successfully applied to breeding selection scenarios even in small local cattle breeds such as Rendena.The maintenance of local cattle breeds is key to selecting for efficient food production, landscape protection, and conservation of biodiversity and local cultural heritage. Rendena is an indigenous cattle breed from the alpine North-East of Italy, selected for dual purpose, but with lesser emphasis given to beef traits. In this situation, increasing accuracy for beef traits could prevent detrimental effects due to the antagonism with milk production. Our study assessed the impact of genomic information on estimated breeding values (EBVs) in Rendena performance-tested bulls. Traits considered were average daily gain, in vivo EUROP score, and in vivo estimate of dressing percentage. The final dataset contained 1691 individuals with phenotypes and 8372 animals in pedigree, 1743 of which were genotyped. Using the cross-validation method, three models were compared: (i) Pedigree-BLUP (PBLUP); (ii) single-step GBLUP (ssGBLUP), and (iii) weighted single-step GBLUP (WssGBLUP). Models including genomic information presented higher accuracy, especially WssGBLUP. However, the model with the best overall properties was the ssGBLUP, showing higher accuracy than PBLUP and optimal values of bias and dispersion parameters. Our study demonstrated that integrating phenotypes for beef traits with genomic data can be helpful to estimate EBVs, even in a small local breed.

Genome-wide scan reveals genomic regions and candidate genes underlying direct and maternal effects of preweaning calf mortality in Nellore cattle

We conducted analysis to estimate genetic parameters and to identify genomic regions and candidate genes affecting direct and maternal effects of preweaning calf mortality (PWM) in Nellore cattle. Phenotypic records of 67,196 animals, and 8443 genotypes for 410,936 SNPs were used. Analysis were performed through the weighted single-step GBLUP approach and considering a threshold animal model via Bayesian Inference. Direct and maternal heritability estimates were of 0.2143 ± 0.0348 and 0.0137 ± 0.0066, respectively. The top 10 genomic regions accounted for 13.61 and 14.23% of the direct and maternal additive genetic variances and harbored a total of 63 and 91 positional candidate genes, respectively. Two overlapping regions on BTA2 were identified for both direct and maternal effects. Candidate genes are involved in biological mechanisms i.e. embryogenesis, immune response, feto-maternal communication, circadian rhythm, hormone alterations, myometrium adaptation, and milk secretion, which are critical for the successful calf growth and survival during preweaning period.

Genome-wide scan for common variants associated with intramuscular fat and moisture content in rainbow trout

BackgroundGenetic improvement of fillet quality attributes is a priority of the aquaculture industry. Muscle composition impacts quality attributes such as flavor, appearance, texture, and juiciness. Fat and moisture make up about ~ 80% of the tissue weight. The genetic architecture underlying the fat and moisture content of the muscle is still to be fully explored in fish. A 50 K gene transcribed SNP chip was used for genotyping 789 fish with available phenotypic data for fat and moisture content. Genotyped fish were obtained from two consecutive generations produced in the National Center for Cool and Cold Water Aquaculture (NCCCWA) growth-selective breeding program. Estimates of SNP effects from weighted single-step GBLUP (WssGBLUP) were used to perform genome-wide association (GWA) analysis to identify quantitative trait loci (QTL) associated with the studied traits.ResultsUsing genomic sliding windows of 50 adjacent SNPs, 137 and 178 SNPs were identified as associated with fat and moisture content, respectively. Chromosomes 19 and 29 harbored the highest number of SNPs explaining at least 2% of the genetic variation in fat and moisture content. A total of 61 common SNPs on chromosomes 19 and 29 affected the aforementioned traits; this association suggests common mechanisms underlying intramuscular fat and moisture content. Additionally, based on single-marker GWA analyses, 8 and 24 SNPs were identified in association with fat and moisture content, respectively.ConclusionSNP-harboring genes were primarily involved in lipid metabolism, cytoskeleton remodeling, and protein turnover. This work provides putative SNP markers that could be prioritized and used for genomic selection in breeding programs.

Genome-wide identification of loci associated with growth in rainbow trout

BackgroundGrowth is a major economic production trait in aquaculture. Improvements in growth performance will reduce time and cost for fish to reach market size. However, genes underlying growth have not been fully explored in rainbow trout.ResultsA previously developed 50 K gene-transcribed SNP chip, containing ~ 21 K SNPs showing allelic imbalances potentially associated with important aquaculture production traits including body weight, muscle yield, was used for genotyping a total of 789 fish with available phenotypic data for bodyweight gain. Genotyped fish were obtained from two consecutive generations produced in the NCCCWA growth-selection breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform a genome-wide association (GWA) analysis to identify quantitative trait loci (QTL) associated with bodyweight gain. Using genomic sliding windows of 50 adjacent SNPs, 247 SNPs associated with bodyweight gain were identified. SNP-harboring genes were involved in cell growth, cell proliferation, cell cycle, lipid metabolism, proteolytic activities, chromatin modification, and developmental processes. Chromosome 14 harbored the highest number of SNPs (n = 50). An SNP window explaining the highest additive genetic variance for bodyweight gain (~ 6.4%) included a nonsynonymous SNP in a gene encoding inositol polyphosphate 5-phosphatase OCRL-1. Additionally, based on a single-marker GWA analysis, 33 SNPs were identified in association with bodyweight gain. The highest SNP explaining variation in bodyweight gain was identified in a gene coding for thrombospondin-1 (THBS1) (R2 = 0.09).ConclusionThe majority of SNP-harboring genes, including OCRL-1 and THBS1, were involved in developmental processes. Our results suggest that development-related genes are important determinants for growth and could be prioritized and used for genomic selection in breeding programs.

Genomic regions associated with principal components for growth, visual score and reproductive traits in Nellore cattle

The search for candidate genes underlying the expression of principal components (PCs) for different traits is a useful tool that allows exploring biological mechanisms associated with the traits of interest. The aim of this study was to identify genomic regions associated with PCs for growth, visual score and reproductive traits in Nellore cattle by performing a genome-wide association study (GWAS). Phenotypic and pedigree data from 355,524 animals and genotypes from 3,519 animals, were used in this investigation. The estimated breeding values (EBV) were obtained from a multi-trait analysis using a mixed linear animal model. The eigen-decomposition of the additive genetic (co)variance matrix among traits was used to calculate the EBVs for the main PCs. The SNP effects were estimated using the weighted single-step GBLUP and the BayesC method. The top-10 ranking windows that explained the highest proportion of variance were identified for further functional analyses. The most important genomic regions were identified on BTA7 and BTA24 for PC1, BTA8 for PC2, and BTA3 and BTA10 for PC3. The functional analyses contributed to unravel biological interpretation of PCs by identifying genes potentially associated with growth, carcass traits, conformation, and fatty acid composition traits. These findings are of relevance to the biological understanding of the PCs and their related biotypes in Nellore cattle, potentially allowing for genetic selection for more specific breeding goals.

Genome-wide association analysis and accuracy of genome-enabled breeding value predictions for resistance to infectious hematopoietic necrosis virus in a commercial rainbow trout breeding population

BackgroundInfectious 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.ResultsTen 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).ConclusionsOur 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.

Genome-Wide Association Study Identifies Genomic Loci Affecting Filet Firmness and Protein Content in Rainbow Trout

Filet quality traits determine consumer satisfaction and affect profitability of the aquaculture industry. Soft flesh is a criterion for fish filet downgrades, resulting in loss of value. Filet firmness is influenced by many factors, including rate of protein turnover. A 50K transcribed gene SNP chip was used to genotype 789 rainbow trout, from two consecutive generations, produced in the USDA/NCCCWA selective breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform genome-wide association (GWA) analyses to identify quantitative trait loci affecting filet firmness and protein content. Applying genomic sliding windows of 50 adjacent SNPs, 212 and 225 SNPs were associated with genetic variation in filet shear force and protein content, respectively. Four common SNPs in the ryanodine receptor 3 gene (RYR3) affected the aforementioned filet traits; this association suggests common mechanisms underlying filet shear force and protein content. Genes harboring SNPs were mostly involved in calcium homeostasis, proteolytic activities, transcriptional regulation, chromatin remodeling, and apoptotic processes. RYR3 harbored the highest number of SNPs (n = 32) affecting genetic variation in shear force (2.29%) and protein content (4.97%). Additionally, based on single-marker analysis, a SNP in RYR3 ranked at the top of all SNPs associated with variation in shear force. Our data suggest a role for RYR3 in muscle firmness that may be considered for genomic- and marker-assisted selection in breeding programs of rainbow trout.

Similar Genetic Architecture with Shared and Unique Quantitative Trait Loci for Bacterial Cold Water Disease Resistance in Two Rainbow Trout Breeding Populations.

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.

Genome-wide association study provides strong evidence of genes affecting the reproductive performance of Nellore beef cows.

Reproductive traits are economically important for beef cattle production; however, these traits are still a bottleneck in indicine cattle since these animals typically reach puberty at older ages when compared to taurine breeds. In addition, reproductive traits are complex phenotypes, i.e., they are controlled by both the environment and many small-effect genes involved in different pathways. In this study, we conducted genome-wide association study (GWAS) and functional analyses to identify important genes and pathways associated with heifer rebreeding (HR) and with the number of calvings at 53 months of age (NC53) in Nellore cows. A total of 142,878 and 244,311 phenotypes for HR and NC53, respectively, and 2,925 animals genotyped with the Illumina Bovine HD panel (Illumina®, San Diego, CA, USA) were used in GWAS applying the weighted single-step GBLUP (WssGBLUP) method. Several genes associated with reproductive events were detected in the 20 most important 1Mb windows for both traits. Significant pathways for HR and NC53 were associated with lipid metabolism and immune processes, respectively. MHC class II genes, detected on chromosome 23 (window 25-26Mb) for NC53, were significantly associated with pregnancy success of Nellore cows. These genes have been proved previously to be associated with reproductive traits such as mate choice in other breeds and species. Our results suggest that genes associated with the reproductive traits HR and NC53 may be involved in embryo development in mammalian species. Furthermore, some genes associated with mate choice may affect pregnancy success in Nellore cattle.

Genome-Wide Association Study for Identifying Loci that Affect Fillet Yield, Carcass, and Body Weight Traits in Rainbow Trout (Oncorhynchus mykiss).

Fillet yield (FY, %) is an economically-important trait in rainbow trout aquaculture that affects production efficiency. Despite that, FY has received little attention in breeding programs because it is difficult to measure on a large number of fish and cannot be directly measured on breeding candidates. The recent development of a high-density SNP array for rainbow trout has provided the needed tool for studying the underlying genetic architecture of this trait. A genome-wide association study (GWAS) was conducted for FY, body weight at 10 (BW10) and 13 (BW13) months post-hatching, head-off carcass weight (CAR), and fillet weight (FW) in a pedigreed rainbow trout population selectively bred for improved growth performance. The GWAS analysis was performed using the weighted single-step GBLUP method (wssGWAS). Phenotypic records of 1447 fish (1.5 kg at harvest) from 299 full-sib families in three successive generations, of which 875 fish from 196 full-sib families were genotyped, were used in the GWAS analysis. A total of 38,107 polymorphic SNPs were analyzed in a univariate model with hatch year and harvest group as fixed effects, harvest weight as a continuous covariate, and animal and common environment as random effects. A new linkage map was developed to create windows of 20 adjacent SNPs for use in the GWAS. The two windows with largest effect for FY and FW were located on chromosome Omy9 and explained only 1.0–1.5% of genetic variance, thus suggesting a polygenic architecture affected by multiple loci with small effects in this population. One window on Omy5 explained 1.4 and 1.0% of the genetic variance for BW10 and BW13, respectively. Three windows located on Omy27, Omy17, and Omy9 (same window detected for FY) explained 1.7, 1.7, and 1.0%, respectively, of genetic variance for CAR. Among the detected 100 SNPs, 55% were located directly in genes (intron and exons). Nucleotide sequences of intragenic SNPs were blasted to the Mus musculus genome to create a putative gene network. The network suggests that differences in the ability to maintain a proliferative and renewable population of myogenic precursor cells may affect variation in growth and fillet yield in rainbow trout.

Assessing the value of phenotypic information from non-genotyped animals for QTL mapping of complex traits in real and simulated populations.

BackgroundQTL mapping through genome-wide association studies (GWAS) is challenging, especially in the case of low heritability complex traits and when few animals possess genotypic and phenotypic information. When most of the phenotypic information is from non-genotyped animals, GWAS can be performed using the weighted single-step GBLUP (WssGBLUP) method, which permits to combine all available information, even that of non-genotyped animals. However, it is not clear to what extent phenotypic information from non-genotyped animals increases the power of QTL detection, and whether factors such as the extent of linkage disequilibrium (LD) in the population and weighting SNPs in WssGBLUP affect the importance of using information from non-genotyped animals in GWAS. These questions were investigated in this study using real and simulated data.ResultsAnalysis of real data showed that the use of phenotypes of non-genotyped animals affected SNP effect estimates and, consequently, QTL mapping. Despite some coincidence, the most important genomic regions identified by the analyses, either using or ignoring phenotypes of non-genotyped animals, were not the same. The simulation results indicated that the inclusion of all available phenotypic information, even that of non-genotyped animals, tends to improve QTL detection for low heritability complex traits. For populations with low levels of LD, this trend of improvement was less pronounced. Stronger shrinkage on SNPs explaining lower variance was not necessarily associated with better QTL mapping.ConclusionsThe use of phenotypic information from non-genotyped animals in GWAS may improve the ability to detect QTL for low heritability complex traits, especially in populations in which the level of LD is high.