new - generation sequencing technologies have substantially enhanced the development of genomic tools to assist in breeding decisions . among many techniques ,
genotyping technologies have enabled breeders or researchers to identify dna regions or quantitative trait loci ( qtl ) associated with a particular phenotypic trait of domesticated animals .
especially , several single - nucleotide polymorphisms ( snps ) and qtl associated with holstein complex traits , including economic traits related to milk production , were investigated in the approximate 10-year wave of genome - wide association studies ( gwass ) . in the near past
, researchers discovered genetic markers and created marker panels for use in performing marker - assisted selection ( mas ) .
although mas provided a first genomic approach to achieve breeding goals to animal breeders , the domesticated animal genome in mas was taken lightly with narrow perspectives .
mas is restricted when it comes to predicting animal capacity , as it depends on a small number of qtls that are tagged by markers associated with each trait .
additionally , the gap between the proportion of phenotypic variance accounted for by the top snps that reach genome - wide significance level in a gwas and the heritability estimated from pedigree analyses remained unexplained for most complex traits of human .
this was called the missing heritability problem , explanations to which have been debated in the field .
recent studies using whole - genome estimation approaches demonstrated that a large proportion of heritability for complex traits of humans can be captured by all common snps on current genotyping arrays .
it implies that there are a large number of variants with an effect that is too small to pass the stringent genome - wide significance level . in animal breeding ,
genomics selection ( gs ) has been developed to overcome restriction factors of mas ; it has covered a large number of variants containing small effects . the fundamental difference between mas and
gs as important animal selection tools is that scale gs embraces a much larger number of markers than mas in predicting animal capacity .
gs uses a dense set of markers from across the entire genome to use all markers containing snps with small effects .
for predictions based on genomic data for humans , a previous study has already estimated the proportion of phenotypic variance .
it was elucidated by the common snps all together on a genotyping array for a range of quantitative traits in a large homogenous human sample , using the whole - genome estimation and partitioning approaches .
the novelty of our research lies in animal capacity prediction based on genomic data due to the fact that animal evaluation is a core content in the breeding industry .
finally , to conduct our animal breeding research on polygenic inheritance , which was likely to be ubiquitous for complex holstein traits associated with milk production traits , the genomic selection technique was performed along with various different analyses .
this study used 462 holsteins from korea , and they had information on three phenotypes related to milk production for parity 1 ( fig .
1 ) . three traits used in this study were milk yield , milk fat , and milk protein ; 339 of 462 holsteins in records of parity 1 had records of parity 2 .
all candidates were measured for a range of quantitative traits through public surveys for livestock improvement .
the genomic dnas were isolated from snivel by using the nasal collection kit and were genotyped with 54,609 snps on the illumina bovinesnp50 beadchip .
we excluded the snps with a missingness rate of 0.05 , minor allele frequency ( maf ) 0.01 , and hardy - weinberg equilibrium test p - value 10 using plink .
then , we removed snps on the sex chromosome and retained 41,099 autosomal snps for further analysis .
we estimated the genetic relationship matrix ( grm ) between all pairs of individuals from all genotyped snps for each trait , we then estimated the variance that can be captured by all snps using the restricted maximum likelihood approach in a mixed linear model : y = xb + gg + e where y is a vector of phenotypes , b is a vector of fixed effects with its incidence matrix x , and gg is a vector of aggregate effects of all snps .
var(gg ) = agg ag is the snp - derived grm , and g is the additive genetic variance .
the proportion of variance explained by all snps is defined as hg = g p p is the phenotypic variance .
when we estimated genetic variance using snp data , we have to consider seasonal effects ( winter , december  february ; spring , march  may ; summer , june  august ; and autumn , september 
november based on birth date ) as environmental effects , because milk production traits are closely related to season .
then , we performed anova test to investigate the relationship between each trait ( parity 1 and 2 ) and seasonal effect .
milk fat had only a close relation to season in the anova test of both parity 1 ( p = 0.0267 ) and 2 ( p = 0.000432 ) .
in addition , using the same method as above but allowing multiple genetic components to be fit simultaneously in the model , we partitioned hg into the contributions of genic ( hgg ) and intergenic ( hgi ) regions of the whole genome across all traits .
the genic regions were defined as 0 kb of the 3 ' and 5 ' untranslated regions ( utrs ) .
a total of 13,297 snps were located within the boundaries of 7001 protein - coding genes for this definition ( 0 kb ) .
the total length of the 7001 protein - coding genes was approximately 585 mb , which covered 21.91% of the genome .
this study used 462 holsteins from korea , and they had information on three phenotypes related to milk production for parity 1 ( fig .
1 ) . three traits used in this study were milk yield , milk fat , and milk protein ; 339 of 462 holsteins in records of parity 1 had records of parity 2 .
all candidates were measured for a range of quantitative traits through public surveys for livestock improvement .
the genomic dnas were isolated from snivel by using the nasal collection kit and were genotyped with 54,609 snps on the illumina bovinesnp50 beadchip .
we excluded the snps with a missingness rate of 0.05 , minor allele frequency ( maf ) 0.01 , and hardy - weinberg equilibrium test p - value 10 using plink .
then , we removed snps on the sex chromosome and retained 41,099 autosomal snps for further analysis .
we estimated the genetic relationship matrix ( grm ) between all pairs of individuals from all genotyped snps for each trait , we then estimated the variance that can be captured by all snps using the restricted maximum likelihood approach in a mixed linear model : y = xb + gg + e where y is a vector of phenotypes , b is a vector of fixed effects with its incidence matrix x , and gg is a vector of aggregate effects of all snps .
var(gg ) = agg ag is the snp - derived grm , and g is the additive genetic variance .
the proportion of variance explained by all snps is defined as hg = g p p is the phenotypic variance .
when we estimated genetic variance using snp data , we have to consider seasonal effects ( winter , december  february ; spring , march  may ; summer , june  august ; and autumn , september 
november based on birth date ) as environmental effects , because milk production traits are closely related to season .
then , we performed anova test to investigate the relationship between each trait ( parity 1 and 2 ) and seasonal effect .
milk fat had only a close relation to season in the anova test of both parity 1 ( p = 0.0267 ) and 2 ( p = 0.000432 ) .
in addition , using the same method as above but allowing multiple genetic components to be fit simultaneously in the model , we partitioned hg into the contributions of genic ( hgg ) and intergenic ( hgi ) regions of the whole genome across all traits .
the genic regions were defined as 0 kb of the 3 ' and 5 ' untranslated regions ( utrs ) .
a total of 13,297 snps were located within the boundaries of 7001 protein - coding genes for this definition ( 0 kb ) .
the total length of the 7001 protein - coding genes was approximately 585 mb , which covered 21.91% of the genome .
these data were collected from 462 cows recruited from 63 farms in south korea , genotyped at 54,609 snps on the illumina bovinesnp50 beadchip .
there were 462 individuals and 41,099 autosomal snps after quality control ( methods section ) .
all individuals were measured for three traits ( milk yield , milk fat , and milk protein ) , which were related to milk production .
we estimated the proportion of variance explained by fitting all snps in a mixed linear model for each of the three traits .
in general , there was a substantial amount of variance explained by all autosomal snps after quality control on the illumina bovinesnp50 beadchip ( g ) for three traits of the records of two parities , with a mean of 43.51% ( a range from 36.3 to 47% ) across all three traits of the records of two parities ( table 1 ) . in the records ,
the estimate of hg was non - zero and reached the nominal significance level ( likelihood ratio test , p = 0.05 ) .
we compared the estimates of hg with the narrow - sense heritability h , estimated from pedigree analyses in the literature based on canada holstein .
the value of hg was similar to h from pedigree analysis in the literature , and all common snps explained most ( average , 43.51% ; range from 36.3% to 47% ) of the narrow - sense heritability ( average , 41.38% ; range from 36.3% to 47% ) , despite the estimates of h being from different country populations ( table 1 ) .
in contrast , when we performed a genome - wide association analysis in the same sample , we tried to identify genome - wide - significant ( bonfferoni correction p = 0.05 ) snps for three traits of parity 1 .
three snps and one snp were significant in milk yield and fat , respectively , and there were no significant snps for milk protein ( supplementary fig .
we hypothesized that there was a major reason for only a few significant large - effect snps associated with holstein traits related to milk production .
first of all , many snps with small effects affected the three traits related to milk production .
second , according to previous studies , there are many common variants associated with the traits in humans at a nominally significant level , but their effect sizes are too small to be genome - wide - significant . using the same method as above but allowing multiple genetic components to be fit simultaneously in the model , we then partitioned hg into the contributions of individual chromosomes for each of the three traits of two parities and plotted the estimate of variance explained by each chromosome ( hc ) against the chromosome length for each trait .
in a previous study , no linear correlation between hc and chromosome length for any particular traits was observed for most traits in humans as strongly as shown in the other studies .
then , the average value of estimates of hg over all of tens traits to reduce the sampling error variance , and the average estimated hc was strongly correlated with chromosome length . however , we dealt with only three traits related to milk production , and we aimed to identify a linear correlation between hc and chromosome length for each trait .
the squared correlation between hc and chromosome length using trait information of parity 1 was 0.46 , 0.532 , and 0.509 for milk yield , fat , and protein , respectively ( fig .
the regression slope of the proportion of the genetic variance attributed to each chromosome on the proportion of the genome represented by each chromosome was significant ( p<0.1 ) , suggesting a proportional relationship between genome length and genetic variance in holstein .
consistent results were achieved from myriad genetic variants , each possessing a small effect throughout the widespread whole genome in holstein .
in addition , we partitioned hg into the contributions of genic ( hgg ) and intergenic ( hgi ) regions of the whole genome ( methods section ) .
we estimated the variance explained by the genic ( hcg ) and intergenic ( hci ) regions of each chromosome .
the numbers of genic and intergenic snps on each chromosome are presented in supplementary table 2 .
we showed that the variance explained by the genic ( intergenic ) regions on each chromosome is also proportional to the total length of the genic ( intergenic ) regions ( fig .
previous studies using the whole - genome estimation approach have shown that common snps explicate a large proportion of heritability for traits in human . the reason why gwass in this field
have not yet identified all common snps that contain information on the amount of variation is mainly that many peculiar effects are too small to pass the stringent genome - wide significance level .
therefore , we hypothesized that inheritance of korean holstein milk production traits would consist of a lot of genetic variants with small effects .
finally , we estimated and partitioned the genetic variance that met certain conditions , tagged by all common snps for three complex traits related to milk production and showing ubiquitous polygenic inheritance of korean holstein .
the estimates of hg for three traits were different from 0 at the significance level ( p<0.1 ) .
the estimate of hg for milk yield of parity 1 was 43.0% ( se = 10.3% ) , which was smaller than the estimate from a study in canada ( hg = 52% ) .
there could be three possible reasons : ( 1 ) there is a difference in the reference population size between korea and canada .
in fact , koreans import a lot of canada holstein semen , but the present reference population is smaller than canada . (
2 ) the korean environment is different from canada . ( 3 ) we can consider only one environmental effect ( season ) due to data size in this study .
the estimate for milk fat ( hg = 47.0% , se = 10.1% ) and protein ( hg = 44.3% , se = 10.2% ) in korea was larger than in canada ( fat hg = 36.9% , protein hg = 42.3% ) .
it is demonstrated by the genome partitioning analysis that there was a significant linear relationship between the estimates of variance explained by individual chromosomes and chromosome length ( fig .
2 ) . if we define the genic region as 0 kb of the utrs , the correlation between variance explained and chromosome length was stronger in the intergenic regions than in the genic regions ( fig .
3 ) . a previous study showed by a number of analyses that the result in this was driven neither by the difference between the number of snps in genic regions and in intergenic regions nor by the difference in maf distribution between genic and intergenic snps .
if trait - associated genetic variants are enriched in functional elements , such as introns and utrs , and diluted in exons , the relationship between the explained variance and dna length will be attenuated in the genic region , as in human .
however , it could possibly be just a sampling problem . because our data size was smaller than our expectancy , some chromosomes
did not reach the significance level ( p = 0.01 ) in the genome partitioning analysis . in the genome partitioning analysis for milk yield using snps in the intergenic region and milk protein using snps in the genic region , the linear regression results were also insignificant , and we could identify some outliers . moreover , when these outliers were removed , the linear regression results were significant ( e in fig .
we hypothesize that some chromosomes could have a slightly larger effect than our expectation , based on the fitted line in special cases ( chromosome 17 or 4 in milk protein ) .
we showed by whole - genome estimation and partitioning analyses that holstein milk production traits appeared to be highly polygenic , which means that there were a large number of genetic variants segregating in the population with a small effect widely distributed across the whole genome .
all common snps expressed most of the heritability on average over all 3 traits analyzed in this study .
the slight , remaining unexplained heritability could be due to causal variants , including common and rare ones that are not well tagged by snps on this illumina snp chip , or possibly because the heritability was over - estimated in the pedigree study . in conclusion , heritability consists of many variants with small effects for holstein milk production traits in korea .
overall , this research will become a cornerstone for genomic selection applications in animal breeding .
supplementary data , including two figures and two tables , can be found with this article online at http://www / genominfo.org / src / sm / gni-13 - 146-s001.pdf
. numbers of genic and intergenic snps on each chromosome variance explained by all the genic and intergenic snps on individual chromosomes for the 3 traits all pairwise phenotypic correlation of the three traits related to milk production .
manhattan plot of genome - wide association studies results for three traits of parity 1 ( milk yield , milk fat , and milk protein , respectively ) .
 previous studies in holstein have shown 35% to 51.8% heritability in milk production traits , such as milk yield , fat , and protein , using pedigree data . 
 other studies in complex human traits could be captured by common single - nucleotide polymorphisms ( snps ) , and their genetic variations , attributed to chromosomes , are in proportion to their length . using genome - wide estimation and partitioning approaches , we analyzed three quantitative holstein traits relevant to milk production in korean holstein data harvested from 462 individuals genotyped for 54,609 snps . 
 for all three traits ( milk yield , fat , and protein ) , we estimated a nominally significant ( p = 0.1 ) proportion of variance explained by all snps on the illumina bovinesnp50 beadchip ( h2 g ) . 
 these common snps explained approximately most of the narrow - sense heritability . 
 longer genomic regions tended to provide more phenotypic variation information , with a correlation of 0.46~0.53 between the estimate of variance explained by individual chromosomes and their physical length . 
 these results suggested that polygenicity was ubiquitous for holstein milk production traits . 
 these results will expand our knowledge on recent animal breeding , such as genomic selection in holstein . 
Introduction
Methods
Korean Holstein data
Genotyped and imputed data
Estimating and partitioning genetic variance using SNP data
Results
Discussion
Supplementary materials