Genomic selection of milk fatty acid composition in Sarda dairy sheep: Effect of different phenotypes and relationship matrices on heritability and breeding value accuracy.

Fatty acid (FA) composition is one of the most important aspects of milk nutritional quality. However, the inclusion of this trait as a breeding goal for dairy species is hampered by the logistics and high costs of phenotype recording. Fourier-transform infrared spectroscopy (FTIR) is a valid and cheap alternative to laboratory gas chromatography (GC) for predicting milk FA composition. Moreover, as for other novel phenotypes, the efficiency of selection for these traits can be enhanced by using genomic data. The objective of this research was to compare traditional versus genomic selection approaches for estimating genetic parameters and breeding values of milk fatty acid composition in dairy sheep using either GC-measured or FTIR-predicted FA as phenotypes. Milk FA profiles were available for a total of 923 Sarda breed ewes. The youngest 100 had their own phenotype masked to mimic selection candidates. Pedigree relationship information and genotypes were available for 923 and 769 ewes, respectively. Three statistical approaches were used: the classical-pedigree-based BLUP, the genomic BLUP that considers the genomic relationship matrix G, and the single-step genomic BLUP (ssGBLUP) where pedigree and genomic relationship matrices are blended into a single H matrix. Heritability estimates using pedigree were lower than ssGBLUP, and very similar between GC and FTIR regarding the statistical approach used. For some FA, mostly associated with animal diet (i.e., C18:2n-6, C18:3n-3), random effect of combination of flock and test date explained a relevant quota of total variance, reducing the heritability estimates accordingly. Genomic approaches (genomic BLUP and ssGBLUP) outperformed the traditional pedigree method both for GC and FTIR FA. Prediction accuracies in the older cohort were larger than the young cohort. Genomic prediction accuracies (obtained using either G or H relationship matrix) in the young cohort of animals, where their own phenotypes were masked, were similar for GC and FTIR. Multiple-trait analysis slightly affected genomic breeding value accuracies. These results suggest that FTIR-predicted milk FA composition could represent a valid option for inclusion in breeding programs.

Genetic parameters of milk fatty acid profile in sheep: comparison between gas chromatographic measurements and Fourier-transform IR spectroscopy predictions.

Fatty acid (FA) composition is a key component of sheep milk nutritional quality. However, breeding for FA composition in dairy sheep is hampered by the logistic and phenotyping costs. This study was aimed at estimating genetic parameters for sheep milk FA and to test the feasibility of their routine measurement by using Fourier-transform IR (FTIR) spectroscopy. Milk FA composition of 989 Sarda ewes farmed in 48 flocks was measured by gas chromatography (FAGC). Moreover, FTIR spectrum was collected for each sample, and it was used to predict FA composition (FAFTIR) by partial least squares regression: 700 ewes were used for estimating model parameters, whereas the remaining 289 animals were used to validate the model. One hundred replicates were performed by randomly assigning animals to estimation and validation data set, respectively. Variance components for both measured and predicted traits were estimated with an animal model that included the fixed effects of parity, days in milking interval, lambing month, province, altitude of flock location, the random effects of flock-test-date and animal genetic additive. Genetic correlations among FAGC, and between corresponding FAGC and FAFTIR were estimated by bivariate analysis. Coefficients of determination between FAGC and FAFTIR ranged from moderate (about 0.50 for odd- and branched-chain FA) to high (about 0.90 for de novo FA) values. Low-to-moderate heritabilities were observed for individual FA (ranging from 0.01 to 0.47). The largest value was observed for GC measured C16:0. Low-to-moderate heritabilities were estimated for FA groups. In most of cases, heritabilites were slightly larger for FAGC than FAFTIR. Estimates of genetic correlations among FAGC showed a large variability in magnitude and sign. The genetic correlation between FAFTIR and FAGC was higher than 60% for all investigated traits. Results of the present study confirm the existence of genetic variability of the FA composition in sheep and suggest the feasibility of using FAFTIR as proxies for these traits in large scale breeding programs.

Genome wide association study on beef production traits in Marchigiana cattle breed.

A genomewide association study was carried out on a sample of Marchigiana breed cattle to detect markers significantly associated with carcass and meat traits. Four hundred and nine young bulls from 117 commercial herds were genotyped by Illumina 50K BeadChip assay. Eight growth and carcass traits (average daily gain, carcass weight, dressing percentage, body weight, skin weight, shank circumference, head weight and carcass conformation) and two meat quality traits (pH at slaughter and pH 24 h after slaughter) were measured. Data were analysed with a linear mixed model that included fixed effects of herd, slaughter date, fixed covariables of age at slaughter and SNP genotype, and random effects of herd and animal. A permutation test was performed to correct SNP genotype significance level for multiple testing. A total of 96 SNPs were significantly associated at genomewide level with one or more of the considered traits. Gene search was performed on genomic regions identified on the basis of significant SNP position and level of linkage disequilibrium. Interesting loci affecting lipid metabolism (SOAT1), bone (BMP4) and muscle (MYOF) biology were highlighted. These results may be useful to better understand the genetic architecture of growth and body composition in cattle.

Genome-wide association analysis in Italian Simmental cows for lactation curve traits using a low-density (7K) SNP panel.

High-throughput cow genotyping has opened new perspectives for genome-wide association studies (GWAS). Directly recorded phenotypes and several records per animal could be used. In this study, a GWAS on lactation curve traits of 337 Italian Simmental cows genotyped with the Illumina (San Diego, CA) low-density BeadChip (7K) was carried out. Scores of the first 2 principal components extracted from test-day records (7 for each lactation) for milk yield, fat and protein percentages, and somatic cell score were used as phenotypes. The first component described the average level of the lactation curve, whereas the second summarized its shape. Data were analyzed with a mixed linear model that included fixed effects of herd, calving month, calving year, parity, SNP genotype, and random effects of animal and permanent environment. All statistically significant markers (Bonferroni corrected) were associated with the average level component (2 for milk yield, 9 for fat percentage, 6 for protein percentages, and 1 for somatic cell score). No markers were found to be associated with the lactation curve shape. Gene discovery was performed using windows of variable size, according to the linkage disequilibrium level of the specific genomic region. Several suggestive candidate genes were identified, some of which already reported to be associated with dairy traits, such as DGAT1. Others were involved in lipid metabolism, in protein synthesis, in the immune response, in cellular processes, and in early development. The large number of genes flagged in the present study suggests interesting perspectives for the use of low-density genotyped females for GWAS, also for novel phenotypes that are not currently considered as breeding goals.

Dissection of genomic correlation matrices of US Holsteins using multivariate factor analysis.

The aim of this study was to compare correlation matrices between direct genomic predictions for 31 traits at the genomic and chromosomal levels in US Holstein bulls. Multivariate factor analysis carried out at the genome level identified seven factors associated with conformation, longevity, yield, feet and legs, fat and protein content traits. Some differences were found at the chromosome level; variations in covariance structure on BTA 6, 14, 18 and 20 were interpreted as evidence of segregating QTL for different groups of traits. For example, milk yield and composition tended to join in a single factor on BTA 14, which is known to harbour the DGAT1 locus that affects these traits. Another example was on BTA 18, where a factor strongly correlated with sire calving ease and conformation traits was identified. It is known that in US Holstein, there is a segregating QTL on BTA18 influencing these traits. Moreover, a possible candidate gene for daughter pregnancy rate was suggested for BTA28. The methodology proposed in this study could be used to identify individual chromosomes, which have covariance structures that differ from the overall (whole genome) covariance structure. Such differences can be difficult to detect when a large number of traits are evaluated, and covariances may be affected by QTL that do not have large allele substitution effects.

Multiple-breed genomic evaluation by principal component analysis in small size populations.

In this study, the effects of breed composition and predictor dimensionality on the accuracy of direct genomic values (DGV) in a multiple breed (MB) cattle population were investigated. A total of 3559 bulls of three breeds were genotyped at 54 001 single nucleotide polymorphisms: 2093 Holstein (H), 749 Brown Swiss (B) and 717 Simmental (S). DGV were calculated using a principal component (PC) approach for either single (SB) or MB scenarios. Moreover, DGV were computed using all SNP genotypes simultaneously with SNPBLUP model as comparison. A total of seven data sets were used: three with a SB each, three with different pairs of breeds (HB, HS and BS), and one with all the three breeds together (HBS), respectively. Editing was performed separately for each scenario. Reference populations differed in breed composition, whereas the validation bulls were the same for all scenarios. The number of SNPs retained after data editing ranged from 36 521 to 41 360. PCs were extracted from actual genotypes. The total number of retained PCs ranged from 4029 to 7284 in Brown Swiss and HBS respectively, reducing the number of predictors by about 85% (from 82% to 89%). In all, three traits were considered: milk, fat and protein yield. Correlations between deregressed proofs and DGV were used to assess prediction accuracy in validation animals. In the SB scenarios, average DGV accuracy did not substantially change when either SNPBLUP or PC were used. Improvement of DGV accuracy were observed for some traits in Brown Swiss, only when MB reference populations and PC approach were used instead of SB-SNPBLUP (+10% HBS, +16%HB for milk yield and +3% HBS and +7% HB for protein yield, respectively). With the exclusion of the abovementioned cases, similar accuracies were observed using MB reference population, under the PC or SNPBLUP models. Random variation owing to sampling effect or size and composition of the reference population may explain the difficulty in finding a defined pattern in the results.

Use of the canonical discriminant analysis to select SNP markers for bovine breed assignment and traceability purposes.

Several market research studies have shown that consumers are primarily concerned with the provenance of the food they eat. Among the available identification methods, only DNA-based techniques appear able to completely prevent frauds. In this study, a new method to discriminate among different bovine breeds and assign new individuals to groups was developed. Bulls of three cattle breeds farmed in Italy - Holstein, Brown, and Simmental - were genotyped using the 50K SNP Illumina BeadChip. Multivariate canonical discriminant analysis was used to discriminate among breeds, and discriminant analysis (DA) was used to assign new observations. This method was able to completely identify the three groups at chromosome level. Moreover, a genome-wide analysis developed using 340 linearly independent SNPs yielded a significant separation among groups. Using the reduced set of markers, the DA was able to assign 30 independent individuals to the proper breed. Finally, a set of 48 high discriminant SNPs was selected and used to develop a new run of the analysis. Again, the procedure was able to significantly identify the three breeds and to correctly assign new observations. These results suggest that an assay with the selected 48 SNP could be used to routinely track monobreed products.

The impact of the rank of marker variance-covariance matrix in principal component evaluation for genomic selection applications.

In genomic selection (GS) programmes, direct genomic values (DGV) are evaluated using information provided by high-density SNP chip. Being DGV accuracy strictly dependent on SNP density, it is likely that an increase in the number of markers per chip will result in severe computational consequences. Aim of present work was to test the effectiveness of principal component analysis (PCA) carried out by chromosome in reducing the marker dimensionality for GS purposes. A simulated data set of 5700 individuals with an equal number of SNP distributed over six chromosomes was used. PCs were extracted both genome-wide (ALL) and separately by chromosome (CHR) and used to predict DGVs. In the ALL scenario, the SNP variance-covariance matrix (S) was singular, positive semi-definite and contained null information which introduces 'spuriousness' in the derived results. On the contrary, the S matrix for each chromosome (CHR scenario) had a full rank. Obtained DGV accuracies were always better for CHR than ALL. Moreover, in the latter scenario, DGV accuracies became soon unsettled as the number of animals decreases, whereas in CHR, they remain stable till 900-1000 individuals. In real applications where a 54k SNP chip is used, the largest number of markers per chromosome is approximately 2500. Thus, a number of around 3000 genotyped animals could lead to reliable results when the original SNP variables are replaced by a reduced number of PCs.