Genes of the RNASE5 pathway contain SNP associated with milk production traits in dairy cattle

Background

Identification of the processes and mutations responsible for the large genetic variation in milk production among dairy cattle has proved challenging. One approach is to identify a biological process potentially involved in milk production and to determine the genetic influence of all the genes included in the process or pathway. Angiogenin encoded by angiogenin, ribonuclease, RNase A family 5 (RNASE5) is relatively abundant in milk, and has been shown to regulate protein synthesis and act as a growth factor in epithelial cells in vitro. However, little is known about the role of angiogenin in the mammary gland or if the polymorphisms present in the bovine RNASE5 gene are associated with lactation and milk production traits in dairy cattle. Given the high economic value of increased protein in milk, we have tested the hypothesis that RNASE5 or genes in the RNASE5 pathway are associated with milk production traits. First, we constructed a “RNASE5 pathway” based on upstream and downstream interacting genes reported in the literature. We then tested SNP in close proximity to the genes of this pathway for association with milk production traits in a large dairy cattle dataset.

Results

The constructed RNASE5 pathway consisted of 11 genes. Association analysis between SNP in 1 Mb regions surrounding these genes and milk production traits revealed that more SNP than expected by chance were associated with milk protein percent (P < 0.05 significance). There was no significant association with other traits such as milk fat content or fertility.

Conclusions

These results support a role for the RNASE5 pathway in milk production, specifically milk protein percent, and indicate that polymorphisms in or near these genes explain a proportion of the variation for this trait. This method provides a novel way of understanding the underlying biology of lactation with implications for milk production and can be applied to any pathway or gene set to test whether they are responsible for the variation of complex traits.     

Partitioning of genomic variance reveals biological pathways associated with udder health and milk production traits in dairy cattle

Background

We have used a linear mixed model (LMM) approach to examine the joint contribution of genetic markers associated with a biological pathway. However, with these markers being scattered throughout the genome, we are faced with the challenge of modelling the contribution from several, sometimes even all, chromosomes at once. Due to linkage disequilibrium (LD), all markers may be assumed to account for some genomic variance; but the question is whether random sets of markers account for the same genomic variance as markers associated with a biological pathway?

Results

We applied the LMM approach to identify biological pathways associated with udder health and milk production traits in dairy cattle. A random gene sampling procedure was applied to assess the biological pathways in a dataset that has an inherently complex genetic correlation pattern due to the population structure of dairy cattle, and to linkage disequilibrium within the bovine genome and within the genes associated to the biological pathway.

Conclusions

Several biological pathways that were significantly associated with health and production traits were identified in dairy cattle; i.e. the markers linked to these pathways explained more of the genomic variance and provided a better model fit than 95 % of the randomly sampled gene groups. Our results show that immune related pathways are associated with production traits, and that pathways that include a causal marker for production traits are identified with our procedure.We are confident that the LMM approach provides a general framework to exploit and integrate prior biological information and could potentially lead to improved understanding of the genetic architecture of complex traits and diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0132-6) contains supplementary material, which is available to authorized users.     

Characterization of bioactive recombinant human lysozyme expressed in milk of cloned transgenic cattle

Background

There is great potential for using transgenic technology to improve the quality of cow milk and to produce biopharmaceuticals within the mammary gland. Lysozyme, a bactericidal protein that protects human infants from microbial infections, is highly expressed in human milk but is found in only trace amounts in cow milk.

Methodology/Principal Findings

We have produced 17 healthy cloned cattle expressing recombinant human lysozyme using somatic cell nuclear transfer. In this study, we just focus on four transgenic cattle which were natural lactation. The expression level of the recombinant lysozyme was up to 25.96 mg/L, as measured by radioimmunoassay. Purified recombinant human lysozyme showed the same physicochemical properties, such as molecular mass and bacterial lysis, as its natural counterpart. Moreover, both recombinant and natural lysozyme had similar conditions for reactivity as well as for pH and temperature stability during in vitro simulations. The gross composition of transgenic and non-transgenic milk, including levels of lactose, total protein, total fat, and total solids were not found significant differences.

Conclusions/Significance

Thus, our study not only describes transgenic cattle whose milk offers the similar nutritional benefits as human milk but also reports techniques that could be further refined for production of active human lysozyme on a large scale.     

Targeted resequencing of GWAS loci reveals novel genetic variants for milk production traits

Background

Genome wide association study (GWAS) has been proven to be a powerful tool for detecting genomic variants associated with complex traits. However, the specific genes and causal variants underlying these traits remain unclear.

Results

Here, we used target-enrichment strategy coupled with next generation sequencing technique to study target regions which were found to be associated with milk production traits in dairy cattle in our previous GWAS. Among the large amount of novel variants detected by targeted resequencing, we selected 200 SNPs for further association study in a population consisting of 2634 cows. Sixty six SNPs distributed in 53 genes were identified to be associated significantly with on milk production traits. Of the 53 genes, 26 were consistent with our previous GWAS results. We further chose 20 significant genes to analyze their mRNA expression in different tissues of lactating cows, of which 15 were specificly highly expressed in mammary gland.

Conclusions

Our study illustrates the potential for identifying causal mutations for milk production traits using target-enrichment resequencing and extends the results of GWAS by discovering new and potentially functional mutations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1105) contains supplementary material, which is available to authorized users.     

The bovine lactation genome: insights into the evolution of mammalian milk

Background

The newly assembled Bos taurus genome sequence enables the linkage of bovine milk and lactation data with other mammalian genomes.

Results

Using publicly available milk proteome data and mammary expressed sequence tags, 197 milk protein genes and over 6,000 mammary genes were identified in the bovine genome. Intersection of these genes with 238 milk production quantitative trait loci curated from the literature decreased the search space for milk trait effectors by more than an order of magnitude. Genome location analysis revealed a tendency for milk protein genes to be clustered with other mammary genes. Using the genomes of a monotreme (platypus), a marsupial (opossum), and five placental mammals (bovine, human, dog, mice, rat), gene loss and duplication, phylogeny, sequence conservation, and evolution were examined. Compared with other genes in the bovine genome, milk and mammary genes are: more likely to be present in all mammals; more likely to be duplicated in therians; more highly conserved across Mammalia; and evolving more slowly along the bovine lineage. The most divergent proteins in milk were associated with nutritional and immunological components of milk, whereas highly conserved proteins were associated with secretory processes.

Conclusions

Although both copy number and sequence variation contribute to the diversity of milk protein composition across species, our results suggest that this diversity is primarily due to other mechanisms. Our findings support the essentiality of milk to the survival of mammalian neonates and the establishment of milk secretory mechanisms more than 160 million years ago.     

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle

Background

Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits.

Results

Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits.

Conclusions

The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1595-0) contains supplementary material, which is available to authorized users.     

Whole genome association study identifies regions of the bovine genome and biological pathways involved in carcass trait performance in Holstein-Friesian cattle

Background

Four traits related to carcass performance have been identified as economically important in beef production: carcass weight, carcass fat, carcass conformation of progeny and cull cow carcass weight. Although Holstein-Friesian cattle are primarily utilized for milk production, they are also an important source of meat for beef production and export. Because of this, there is great interest in understanding the underlying genomic structure influencing these traits. Several genome-wide association studies have identified regions of the bovine genome associated with growth or carcass traits, however, little is known about the mechanisms or underlying biological pathways involved. This study aims to detect regions of the bovine genome associated with carcass performance traits (employing a panel of 54,001 SNPs) using measures of genetic merit (as predicted transmitting abilities) for 5,705 Irish Holstein-Friesian animals. Candidate genes and biological pathways were then identified for each trait under investigation.

Results

Following adjustment for false discovery (q-value < 0.05), 479 quantitative trait loci (QTL) were associated with at least one of the four carcass traits using a single SNP regression approach. Using a Bayesian approach, 46 QTL were associated (posterior probability > 0.5) with at least one of the four traits. In total, 557 unique bovine genes, which mapped to 426 human orthologs, were within 500kbs of QTL found associated with a trait using the Bayesian approach. Using this information, 24 significantly over-represented pathways were identified across all traits. The most significantly over-represented biological pathway was the peroxisome proliferator-activated receptor (PPAR) signaling pathway.

Conclusions

A large number of genomic regions putatively associated with bovine carcass traits were detected using two different statistical approaches. Notably, several significant associations were detected in close proximity to genes with a known role in animal growth such as glucagon and leptin. Several biological pathways, including PPAR signaling, were shown to be involved in various aspects of bovine carcass performance. These core genes and biological processes may form the foundation for further investigation to identify causative mutations involved in each trait. Results reported here support previous findings suggesting conservation of key biological processes involved in growth and metabolism.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-837) contains supplementary material, which is available to authorized users.     

Differential temporal expression of milk miRNA during the lactation cycle of the marsupial tammar wallaby (Macropus eugenii)

Background

Lactation is a key aspect of mammalian evolution for adaptation of various reproductive strategies along different mammalian lineages. Marsupials, such as tammar wallaby, adopted a short gestation and a relatively long lactation cycle, the newborn is immature at birth and significant development occurs postnatally during lactation. Continuous changes of tammar milk composition may contribute to development and immune protection of pouch young. Here, in order to address the putative contribution of newly identified secretory milk miRNA in these processes, high throughput sequencing of miRNAs collected from tammar milk at different time points of lactation was conducted. A comparative analysis was performed to find distribution of miRNA in milk and blood serum of lactating wallaby.

Results

Results showed that high levels of miRNA secreted in milk and allowed the identification of differentially expressed milk miRNAs during the lactation cycle as putative markers of mammary gland activity and functional candidate signals to assist growth and timed development of the young. Comparative analysis of miRNA distribution in milk and blood serum suggests that milk miRNAs are primarily expressed from mammary gland rather than transferred from maternal circulating blood, likely through a new putative exosomal secretory pathway. In contrast, highly expressed milk miRNAs could be detected at significantly higher levels in neonate blood serum in comparison to adult blood, suggesting milk miRNAs may be absorbed through the gut of the young.

Conclusion

The function of miRNA in mammary gland development and secretory activity has been proposed, but results from the current study also support a differential role of milk miRNA in regulation of development in the pouch young, revealing a new potential molecular communication between mother and young during mammalian lactation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1012) contains supplementary material, which is available to authorized users.     

Genomic analysis of dominance effects on milk production and conformation traits in Fleckvieh cattle

Background

Estimates of dominance variance in dairy cattle based on pedigree data vary considerably across traits and amount to up to 50% of the total genetic variance for conformation traits and up to 43% for milk production traits. Using bovine SNP (single nucleotide polymorphism) genotypes, dominance variance can be estimated both at the marker level and at the animal level using genomic dominance effect relationship matrices. Yield deviations of high-density genotyped Fleckvieh cows were used to assess cross-validation accuracy of genomic predictions with additive and dominance models. The potential use of dominance variance in planned matings was also investigated.

Results

Variance components of nine milk production and conformation traits were estimated with additive and dominance models using yield deviations of 1996 Fleckvieh cows and ranged from 3.3% to 50.5% of the total genetic variance. REML and Gibbs sampling estimates showed good concordance. Although standard errors of estimates of dominance variance were rather large, estimates of dominance variance for milk, fat and protein yields, somatic cell score and milkability were significantly different from 0. Cross-validation accuracy of predicted breeding values was higher with genomic models than with the pedigree model. Inclusion of dominance effects did not increase the accuracy of the predicted breeding and total genetic values. Additive and dominance SNP effects for milk yield and protein yield were estimated with a BLUP (best linear unbiased prediction) model and used to calculate expectations of breeding values and total genetic values for putative offspring. Selection on total genetic value instead of breeding value would result in a larger expected total genetic superiority in progeny, i.e. 14.8% for milk yield and 27.8% for protein yield and reduce the expected additive genetic gain only by 4.5% for milk yield and 2.6% for protein yield.

Conclusions

Estimated dominance variance was substantial for most of the analyzed traits. Due to small dominance effect relationships between cows, predictions of individual dominance deviations were very inaccurate and including dominance in the model did not improve prediction accuracy in the cross-validation study. Exploitation of dominance variance in assortative matings was promising and did not appear to severely compromise additive genetic gain.     

Causal relationships between milk quality and coagulation properties in Italian Holstein-Friesian dairy cattle

Background

Recently, selection for milk technological traits was initiated in the Italian dairy cattle industry based on direct measures of milk coagulation properties (MCP) such as rennet coagulation time (RCT) and curd firmness 30 min after rennet addition (a₃₀) and on some traditional milk quality traits that are used as predictors, such as somatic cell score (SCS) and casein percentage (CAS). The aim of this study was to shed light on the causal relationships between traditional milk quality traits and MCP. Different structural equation models that included causal effects of SCS and CAS on RCT and a₃₀ and of RCT on a₃₀ were implemented in a Bayesian framework.

Results

Our results indicate a non-zero magnitude of the causal relationships between the traits studied. Causal effects of SCS and CAS on RCT and a₃₀ were observed, which suggests that the relationship between milk coagulation ability and traditional milk quality traits depends more on phenotypic causal pathways than directly on common genetic influence. While RCT does not seem to be largely controlled by SCS and CAS, some of the variation in a₃₀ depends on the phenotypes of these traits. However, a₃₀ depends heavily on coagulation time. Our results also indicate that, when direct effects of SCS, CAS and RCT are considered simultaneously, most of the overall genetic variability of a₃₀ is mediated by other traits.

Conclusions

This study suggests that selection for RCT and a₃₀ should not be performed on correlated traits such as SCS or CAS but on direct measures because the ability of milk to coagulate is improved through the causal effect that the former play on the latter, rather than from a common source of genetic variation. Breaking the causal link (e.g. standardizing SCS or CAS before the milk is processed into cheese) would reduce the impact of the improvement due to selective breeding. Since a₃₀ depends heavily on RCT, the relative emphasis that is put on this trait should be reconsidered and weighted for the fact that the pure measure of a₃₀ almost double-counts RCT.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0123-7) contains supplementary material, which is available to authorized users.     

Over-expression of human lipoprotein lipase in mouse mammary glands leads to reduction of milk triglyceride and delayed growth of suckling pups

Background

The mammary gland is a conserved site of lipoprotein lipase expression across species and lipoprotein lipase attachment to the luminal surface of mammary gland vascular endothelial cells has been implicated in the direction of circulating triglycerides into milk synthesis during lactation.

Principal Findings

Here we report generation of transgenic mice harboring a human lipoprotein lipase gene driven by a mammary gland-specific promoter. Lipoprotein lipase levels in transgenic milk was raised to 0.16 mg/ml, corresponding to an activity of 8772.95 mU/ml. High lipoprotein lipase activity led to a significant reduction of triglyceride concentration in milk, but other components were largely unchanged. Normal pups fed with transgenic milk showed inferior growth performances compared to those fed with normal milk.

Conclusion

Our study suggests a possibility to reduce the triglyceride content of cow milk using transgenic technology.     

Effectiveness of genomic prediction on milk flow traits in dairy cattle

Background

Milkability, primarily evaluated by measurements of milking speed and time, has an economic impact in milk production of dairy cattle. Recently the Italian Brown Swiss Breeders Association has included milking speed in genetic evaluations. The main objective of this study was to investigate the possibility of implementing genomic selection for milk flow traits in the Italian Brown Swiss population and thereby evaluate the potential of genomic selection for novel traits in medium-sized populations. Predicted breeding values and reliabilities based on genomic information were compared with those obtained from traditional breeding values.

Methods

Milk flow measures for total milking time, ascending time, time of plateau, descending time, average milk flow and maximum milk flow were collected on 37 213 Italian Brown Swiss cows. Breeding values for genotyped sires (n = 1351) were obtained from standard BLUP and genome-enhanced breeding value techniques utilizing two-stage and single-step methods. Reliabilities from a validation dataset were estimated and used to compare accuracies obtained from parental averages with genome-enhanced predictions.

Results

Genome-enhanced breeding values evaluated using two-stage methods had similar reliabilities with values ranging from 0.34 to 0.49 for the different traits. Across two-stage methods, the average increase in reliability from parental average was approximately 0.17 for all traits, with the exception of descending time, for which reliability increased to 0.11. Combining genomic and pedigree information in a single-step produced the largest increases in reliability over parent averages: 0.20, 0.24, 0.21, 0.14, 0.20 and 0.21 for total milking time, ascending time, time of plateau, descending time, average milk flow and maximum milk flow, respectively.

Conclusions

Using genomic models increased the accuracy of prediction compared to traditional BLUP methods. Our results show that, among the methods used to predict genome-enhanced breeding values, the single-step method was the most successful at increasing the reliability for most traits. The single-step method takes advantage of all the data available, including phenotypes from non-genotyped animals, and can easily be incorporated into current breeding evaluations.     

Overexpression of miR-30b in the Developing Mouse Mammary Gland Causes a Lactation Defect and Delays Involution

Background

MicroRNA (miRNA) are negative regulators of gene expression, capable of exerting pronounced influences upon the translation and stability of mRNA. They are potential regulators of normal mammary gland development and of the maintenance of mammary epithelial progenitor cells. This study was undertaken to determine the role of miR-30b on the establishment of a functional mouse mammary gland. miR-30b is a member of the miR-30 family, composed of 6 miRNA that are highly conserved in vertebrates. It has been suggested to play a role in the differentiation of several cell types.

Methodology/Principal Findings

The expression of miR-30b was found to be regulated during mammary gland development. Transgenic mice overexpressing miR-30b in mammary epithelial cells were used to investigate its role. During lactation, mammary histological analysis of the transgenic mice showed a reduction in the size of alveolar lumen, a defect of the lipid droplets and a growth defect of pups fed by transgenic females. Moreover some mammary epithelial differentiated structures persisted during involution, suggesting a delay in the process. The genes whose expression was affected by the overexpression of miR-30b were characterized by microarray analysis.

Conclusion/Significance

Our data suggests that miR-30b is important for the biology of the mammary gland and demonstrates that the deregulation of only one miRNA could affect lactation and involution.     

Database of cattle candidate genes and genetic markers for milk production and mastitis

A cattle database of candidate genes and genetic markers for milk production and mastitis has been developed to provide an integrated research tool incorporating different types of information supporting a genomic approach to study lactation, udder development and health. The database contains 943 genes and genetic markers involved in mammary gland development and function, representing candidates for further functional studies. The candidate loci were drawn on a genetic map to reveal positional overlaps. For identification of candidate loci, data from seven different research approaches were exploited: (i) gene knockouts or transgenes in mice that result in specific phenotypes associated with mammary gland (143 loci); (ii) cattle QTL for milk production (344) and mastitis related traits (71); (iii) loci with sequence variations that show specific allele-phenotype interactions associated with milk production (24) or mastitis (10) in cattle; (iv) genes with expression profiles associated with milk production (207) or mastitis (107) in cattle or mouse; (v) cattle milk protein genes that exist in different genetic variants (9); (vi) miRNAs expressed in bovine mammary gland (32) and (vii) epigenetically regulated cattle genes associated with mammary gland function (1). Fourty-four genes found by multiple independent analyses were suggested as the most promising candidates and were further in silico analysed for expression levels in lactating mammary gland, genetic variability and top biological functions in functional networks. A miRNA target search for mammary gland expressed miRNAs identified 359 putative binding sites in 3'UTRs of candidate genes.     

Mapping quantitative trait loci (QTL) in sheep. IV. Analysis of lactation persistency and extended lactation traits in sheep

Background

In sheep dairy production, total lactation performance, and length of lactation of lactation are of economic significance. A more persistent lactation has been associated with improved udder health. An extended lactation is defined by a longer period of milkability. This study is the first investigation to examine the presence of quantitative trait loci (QTL) for extended lactation and lactation persistency in sheep.

Methods

An (Awassi × Merino) × Merino single-sire backcross family with 172 ewes was used to map QTL for lactation persistency and extended lactation traits on a framework map of 189 loci across all autosomes. The Wood model was fitted to data from multiple lactations to estimate parameters of ovine lactation curves, and these estimates were used to derive measures of lactation persistency and extended lactation traits of milk, protein, fat, lactose, useful yield, and somatic cell score. These derived traits were subjected to QTL analyses using maximum likelihood estimation and regression analysis.

Results

Overall, one highly significant (LOD > 3.0), four significant (2.0 < LOD < 3.0) and five suggestive (1.7 < LOD < 2.0) QTL were detected across all traits in common by both mapping methods. One additional suggestive QTL was identified using maximum likelihood estimation, and four suggestive (0.01 < P < 0.05) and two significant (P < 0.01) QTL using the regression approach only. All detected QTL had effect sizes in the range of 0.48 to 0.64 SD, corresponding to QTL heritabilities of 3.1 to 8.9%. The comparison of the detected QTL with results in cattle showed conserved linkage regions. Most of the QTL identified for lactation persistency and extended lactation did not coincide. This suggests that persistency and extended lactation for the same as well as different milk yield and component traits are not controlled by the same genes.

Conclusion

This study identified ten novel QTL for lactation persistency and extended lactation in sheep, but results suggest that lactation persistency and extended lactation do not have a major gene in common. These results provide a basis for further validation in extended families and other breeds as well as targeting regions for genome-wide association mapping using high-density SNP arrays.     

The SNPs in the ACACA gene are effective on fatty acid composition in holstein milk

Fatty acid composition is an important economic trait for both dairy and beef cattle and controlled by genetic factors. Candidate genes controlling fatty acid composition may be found in fat synthesis and metabolism pathways. Acetyl-CoA carboxylase is the flux-determining enzyme in the regulation of fatty acid synthesis in animal tissues. One of two isozymes of this enzyme, acetyl-CoA carboxylase-α (ACACA), catalyses the first committed step of fatty acid synthesis in mammalian cytosol, leading to the biosynthesis of long-chain fatty acids. In the current study, the sequence comparison of the coding sequence (CDS) and two promoter regions (PIA and PIII) in bovine ACACA gene was performed between Japanese Black and Holstein cattle to detect nucleotide polymorphisms influencing fatty acid composition in milk and beef. Five single nucleotide polymorphisms (SNPs) were identified in the CDS region, 28 SNPs in the PIA region and three SNPs in the PIII region. Association study revealed that CCT/CCT type of PIII_#1, #2/PIA_#26 indicated a higher percentage of C14:0 in the milk of the Holstein cattle than CCT/GTC type (p = 0.050) and that a difference of the percentage of C16:0 was observed between CCT/CCT and GTC/GTC type (p = 0.023). CDS_#2 T/T type indicated a higher percentage of C18:0 than T/C type (p = 0.008). In addition, the Japanese Black cattle with CC/GT type of PIII_#1, #2 showed a higher percentage of C18:2 in the meat than those with GT/GT type (p = 0.025). Since PIII is the promoter specific to mammary gland during lactation, the altered expression of the ACACA gene owing to the SNPs in the PIII region may influence the fatty acid composition in the milk.     

Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins

Background

Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits.

Results

In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs.

Conclusion

We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-683) contains supplementary material, which is available to authorized users.