Search

Detecting genetic variation such as Copy Number Variants (CNVs) in cattle provides the opportunity to study their association with quantitative traits. CNVs are DNA sequences of 50 bp up to several Mb long, which can vary in copy number in comparison with a reference genome. The aim of this study was to investigate CNVs in 1410 samples of the Brown Swiss cattle breed using Illumina Bovine HD SNP chip information, which includes 777,962 SNPs. After stringent quality control, CNVs were called with the Golden Helix SVS 8.3.1 (SVS) and PennCNV software and were summarized to CNV regions (CNVRs) at a population level (i.e. overlapping CNVs), using BEDTools. Additionally, common CNVRs between the two software were set as consensus regions. Genes within consensus CNVRs were annotated with a GO analysis using the DAVID Bioinformatics Resources 6.7. In order to validate these results, quantitative PCRs were executed on 15 selected CNVRs. The SVS software identified 25,030 CNVs summarized to 398 CNVRs, which comprised 30 gains, 344 losses and 24 complex CNVRs (i.e. containing both losses and gains), covering 3.92% of the bovine genome. The PennCNV software identified 62,341 CNVs summarized to 5578 CNVRs, which comprised 2638 gains, 2404 losses and 537 complex CNVRs, covering 7.68% of the bovine genome. The length of these CNVRs ranged from 1244 bp to 1,381,355 bp. A total of 563 consensus CNVRs were found covering 2.29% of the UMD 3.1 bovine genome assembly. Of these, 24 were gains, 300 were losses and 239 were complex CNVRs. A total of 775 official gene IDs were annotated in the consensus CNVRs. Among the 537 genes with functional information, the GO and pathway analysis was reported for those who clustered with a p-value

We targeted quantitative trait loci (QTL) for milk protein percentage (P%) in two Italian Holstein granddaughter design families using selective genotyping in combination with high throughput amplified fragment length polymorphism (AFLP) markers. A total of 64 extreme high and low sires in respect to estimated breeding value (EBV) for P% (EBVP%) were genotyped with 25 AFLP primer combinations that revealed 305 and 291 polymorphisms in the two families. Association between markers and EBVP% was investigated by a linear model only on bands having paternal origin (105 and 96 AFLP bands in family D and S, respectively). Although no marker was significantly associated with the target trait after correction for multiple comparisons, 17 AFLP markers, significant without correction for multiple tests, were considered suggestive of the presence of a QTL. Eleven of these were successfully located on sixBos taurus(BTA) chromosomes by radiation hybrid or in-silico mapping. Ten of these mapped in the immediate neighbourhood (less than 10 cM) of already described QTL for P%. Suggestive association was verified in four regions by microsatellites analysis: one on BTA 10; one on BTA 28; and two on BTA 18. Microsatellites identified significant effects by single marker and interval mapping analyses on BTA 10 and BTA 28, while they were only suggestive of the presence of QTL on BTA 18. In summary, our results firstly indicate that AFLP markers may be used to seek QTL exploiting a selective genotyping approach in GDD, a wide used experimental design in cattle; secondly, propose two approaches for AFLP mapping, namely in-silico mapping exploiting most updated release from the bovine whole genome sequencing project, and physical mapping exploiting a panel of Bovine/Hamster Radiation Hybrids; and thirdly, provide new information on QTLs for an economic important trait in a never investigated Holstein cattle population. AFLP in combination with selective genotyping can be a useful strategy for QTL searching in minor livestock species, sometimes having large economic impact in marginal areas, where more informative markers are still poorly developed.

Quantitative trait loci (QTL) mapping projects have been implemented mainly in the Holstein dairy cattle breed for several traits. The aim of this study is to map QTL for milk yield (MY) and milk protein percent (PP) in the Brown Swiss cattle populations of Austria, Germany, and Italy, considered in this study as a single population. A selective DNA pooling approach using milk samples was applied to map QTL in 10 paternal half-sib daughter families with offspring spanning from 1,000 to 3,600 individuals per family. Three families were sampled in Germany, 3 in Italy, 1 in Austria and 3 jointly in Austria and Italy. The pools comprised the 200 highest and 200 lowest performing daughters, ranked by dam-corrected estimated breeding value for each sire-trait combination. For each tail, 2 independent pools, each of 100 randomly chosen daughters, were constructed. Sire marker allele frequencies were obtained by densitometry and shadow correction analyses of 172 genome-wide allocated autosomal markers. Particular emphasis was placed on Bos taurus chromosomes 3, 6, 14, and 20. Marker association for MY and PP with a 10% false discovery rate resulted in nominal P-values of 0.071 and 0.073 for MY and PP, respectively. Sire marker association tested at a 20% false discovery rate (within significant markers) yielded nominal P-values of 0.031 and 0.036 for MY and PP, respectively. There were a total of 36 significant markers for MY, 33 for PP, and 24 for both traits; 75 markers were not significant for any of the traits. Of the 43 QTL regions found in the present study, 10 affected PP only, 8 affected MY only, and 25 affected MY and PP. Remarkably, all 8 QTL regions that affected only MY in the Brown Swiss, also affected MY in research reported in 3 Web-based QTL maps used for comparison with the findings of this study (http://www.vetsci.usyd.edu.au/reprogen/QTL_Map/; http://www.animalgenome.org/QTLdb/cattle.html; http://bovineqtl.tamu.edu/). Similarly, all 10 QTL regions in the Brown Swiss that affected PP only, affected only PP in the databases. Thus, many QTL appear to be common to Brown Swiss and other breeds in the databases (mainly Holstein), and an appreciable fraction of QTL appears to affect MY or PP primarily or exclusively, with little or no effect on the other trait. Although QTL information available today in the Brown Swiss population can be utilized only in a within family marker-assisted selection approach, knowledge of QTL segregating in the whole population should boost gene identification and ultimately the implementation and efficiency of an individual genomic program.

The aim of this study was to evaluate the variations of protein, casein, saturated (SFA), unsaturated (UFA), monounsaturated (MUFA), polyunsaturated (PUFA) fatty acids contents and cheese yield in the milk of two groups of Italian Brown cows conventionally reared in indoor period of housing or consuming pasture during the summer months in 2008 and 2013. Milk components were obtained from samples collected during the national routine (conventionally reared) and ‘extraordinary’ (pasture period) milk recording scheme in herds located near Sondrio (Lombardia, Italy). Milk samples were processed with the MilkoScanTM FT6000 for the identification of milk casein, SFA, UFA, MUFA and PUFA composition. The groups were analysed separately per year and the environmental factors affecting milk protein, casein, and fatty acids contents (pasture/indoor, parity, data of sampling, days in milk, days from collection to analysis) were included in the MIXED procedure of SAS 9.3. A total of 778 milk samples were available, including 234 records from indoor and 544 observations from pasture feeding. Pasture intake affected the content of casein (%) and the proportion of fat in milk (g/100 g), enhancing milk casein levels (from 2·90 to 3) and reducing the concentration of milk SFA in milk from grazing cows (from 2·29 to 1·92). Additionally, the cheese yield was calculated as ‘kg of cheese per 100 kg of milk’ and resulted to be 10·4 and 12 in 2008 from milk of cows reared indoor and with pasture based diet, respectively. The dairy industry should take advantage of the milk production during grazing periods from which high quality products may be obtained.

The determination of copy number variation (CNV) is very important for the evaluation of genomic traits in several species because they are a major source for the genetic variation, influencing gene expression, phenotypic variation, adaptation and the development of diseases. The aim of this study was to obtain a CNV genome map using the Illumina Bovine SNP50 BeadChip data of 651 bulls of the Italian Brown Swiss breed. PennCNV and SVS7 (Golden Helix) software were used for the detection of the CNVs and Copy Number Variation Regions (CNVRs). A total of 5,099 and 1,289 CNVs were identified with PennCNV and SVS7 software, respectively. These were grouped at the population level into 1101 (220 losses, 774 gains, 107 complex) and 277 (185 losses, 56 gains and 36 complex) CNVR. Ten of the selected CNVR were experimentally validated with a qPCR experiment. The GO and pathway analyses were conducted and they identified genes (false discovery rate corrected) in the CNVR related to biological processes cellular component, molecular function and metabolic pathways. Among those, we found the FCGR2B, PPARα, KATNAL1, DNAJC15, PTK2, TG, STAT family, NPM1, GATA2, LMF1, ECHS1 genes, already known in literature because of their association with various traits in cattle. Although there is variability in the CNVRs detection across methods and platforms, this study allowed the identification of CNVRs in Italian Brown Swiss, overlapping those already detected in other breeds and finding additional ones, thus producing new knowledge for association studies with traits of interest in cattle.

A selective DNA pooling approach was applied to identify QTL for conjugated linoleic acid (CLA), vaccenic acid (VA) and Δ(9) -desaturase (D9D) milk content in Italian Brown Swiss dairy cattle. Milk samples from 60 animals with higher values (after correction for environmental factors) and 60 animals with lower values for each of these traits from each of five half-sib families were pooled separately. The pools were genotyped using the Illumina BovineSNP50 BeadChip. Sire allele frequencies were compared between high and low tails at the sire and marker level for SNPs for which the sires were heterozygous. An r procedure was implemented to perform data analysis in a selective DNA pooling design. A correction for multiple tests was applied using the proportion of false positives among all test results. BTA 19 showed the largest number of markers in association with CLA. Associations between SNPs and the VA and Δ(9) -desaturase traits were found on several chromosomes. A bioinformatics survey identified genes with an important role in pathways for milk fat and fatty acids metabolism within 1 Mb of SNP markers associated with fatty acids contents.

The aim of this study is to evaluate the effects of interaction between sire of cow and service sire on the success/unsuccess of inseminations. Data from insemination events of Italian Brown Swiss cows collected from January 1993 through August 2007 were restricted to repeat breeder cows. A cluster analysis was carried out to group herds with very few observations in clusters with at least 15 observations. The edited data set included 102,710 services of 10,708 cows, daughters of 1,716 sires and mated to 3,108 service sires. The success or unsuccess at each insemination was evaluated by a linear mixed model including the fixed effects of herd-year interaction, month of insemination, age, and the random effects of sire service-sire of cow interaction and residual. The distribution of bull combination estimates was bimodal. When the tails of distribution (best and worst 5% of estimates) were considered, 271 service sires were included in both tails. Results suggest that major gene can affect the survival of embryos and that positive or negative interactions between paternal and maternal genotype can affect this reproductive trait.

Background Mastitis is a major disease of dairy cattle occurring in response to environmental exposure to infective agents with a great economic impact on dairy industry. Somatic cell count (SCC) and its log transformation in somatic cell score (SCS) are traits that have been used as indirect measures of resistance to mastitis for decades in selective breeding. A selective DNA pooling (SDP) approach was applied to identify Quantitative Trait Loci (QTL) for SCS in Valdostana Red Pied cattle using the Illumina Bovine HD BeadChip. Results A total of 171 SNPs reached the genome-wide significance for association with SCS. Fifty-two SNPs were annotated within genes, some of those involved in the immune response to mastitis. On BTAs 1, 2, 3, 4, 9, 13, 15, 17, 21 and 22 the largest number of markers in association to the trait was found. These regions identified novel genomic regions related to mastitis (1-Mb SNP windows) and confirmed those already mapped. The largest number of significant SNPs exceeding the threshold for genome-wide significant signal was found on BTA 15, located at 50.43-51.63 Mb. Conclusions The genomic regions identified in this study contribute to a better understanding of the genetic control of the mastitis immune response in cattle and may allow the inclusion of more detailed QTL information in selection programs. Electronic supplementary material The online version of this article (doi:10.1186/s12863-014-0106-7) contains supplementary material, which is available to authorized users.

Bovine respiratory disease (BRD) is the leading cause of morbidity and mortality in feedlot cattle, caused by multiple pathogens that become more virulent in response to stress. As clinical signs often go undetected and various preventive strategies failed, identification of genes affecting BRD is essential for selection for resistance. Selective DNA pooling (SDP) was applied in a genome wide association study (GWAS) to map BRD QTLs in Israeli Holstein male calves. Kosher scoring of lung adhesions was used to allocate 122 and 62 animals to High (Glatt Kosher) and Low (Non-Kosher) resistant groups, respectively. Genotyping was performed using the Illumina BovineHD BeadChip according to the Infinium protocol. Moving average of -logP was used to map QTLs and Log drop was used to define their boundaries (QTLRs). The combined procedure was efficient for high resolution mapping. Nineteen QTLRs distributed over 13 autosomes were found, some overlapping previous studies. The QTLRs contain polymorphic functional and expression candidate genes to affect kosher status, with putative immunological and wound healing activities. Kosher phenotyping was shown to be a reliable means to map QTLs affecting BRD morbidity.

Summary We report on a complete genome scan for quantitative trait loci (QTL) affecting milk protein percentage (PP) in the Italian Holstein-Friesian cattle population, applying a selective DNA pooling strategy in a daughter design. Ten Holstein-Friesian sires were chosen, and for each sire, about 200 daughters, each from the high and low tails of estimated breeding value for PP, were used to construct milk DNA pools. Sires and pools were genotyped for 181 dinucleotide microsatellites covering all cattle autosomes. Sire marker allele frequencies in the pools were obtained by shadow correction of peak height in the electropherograms. After quality control, pool data from eight sires were used for all subsequent analyses. The QTL heterozygosity estimate was lower than that of similar studies in other cattle populations. Multiple marker mapping identified 19 QTL located on 14 chromosomes (BTA1, 2, 3, 4, 5, 6, 8, 9, 12, 14, 17, 20, 23 and 27). The sires were also genotyped for seven polymorphic sites in six candidate genes (ABCG2, SPP1, casein kappa, DGAT1, GHR and PRLR) located within QTL regions of BTA6, 14 and 20 found in this study. The results confirmed or excluded the involvement of some of the analysed markers as the causative polymorphic sites of the identified QTL. The QTL identified, combined with genotype data of these candidate genes, will help to identify other quantitative trait genes and clarify the complex QTL patterns observed for a few chromosomes. Overall, the results are consistent with the Italian Holstein population having been under long-term selection for high PP.

The aim of this study was to evaluate the allele frequency distribution at the DGAT1 p.K232A polymorphic site in seven Italian dairy and dual purpose cattle breeds. On the whole, 651 animals belonging to Italian Holstein (116), Italian Brown (105), Italian Simmental (95), Valdostana Red Pied (95), Rendena (62), Reggiana (128) and Modenese (50) breeds were genotyped by PCR-RFLP. Sequencing was carried out to confirm results of the genotyping protocol. The DGAT1 p.232K allele was identified in Italian Holstein (25.4%), Reggiana (17.2%), and with very low frequency in Italian Simmental, Valdostana Red Pied and Rendena (

Great interest was aroused by reports, based on microsatellite markers, of high levels of statistically significant long-range and nonsyntenic linkage disequilibrium (LD) in livestock. Simulation studies showed that this could result from population family structure. In contrast, recent SNP-based studies of livestock populations report much lower levels of LD. In this study we show, on the basis of microsatellite data from four cattle populations, that high levels of long-range LD are indeed obtained when using the multi-allelic D′ measure of LD. Long-range and nonsyntenic LD are exceedingly low, however, when evaluated by the standardized chi-square measure of LD, which stands in relation to the predictive ability of LD. Furthermore, specially constructed study populations provided no evidence for appreciable LD resulting from family structure at the grandparent level. We propose that the high statistical significance and family structure effects observed in the earlier studies are due to the use of large sample sizes, which accord high statistical significance to even slight deviations from asymptotic expectations under the null hypothesis. Nevertheless, even after taking sample size into account, our results indicate that microsatellites testify to the presence of usable LD at considerably wider separation distances than SNPs, suggesting that use of SNP haplotypes may considerably increase the usefulness of a given fixed SNP array.

Summary The bovine oxidized LDL receptor 1 (OLR1) was chosen as a candidate gene for association tests with milk composition traits. Genotyping of 773 Italian Brown Swiss for a SNP at position 8232 in OLR1 (NW_215807:g.8232C>A) revealed a frequency of 0.95 for the g.8232C allele. The University of Wisconsin Holstein resource population was genotyped for the OLR1 NW_215807:g.[7160C>T; 7161A>G; 8232C>A] SNPs, and four haplotypes were inferred based on the genotypes of sires and their daughters. Oxidized LDL receptor 1 haplotypes were significantly associated with fat percentage (P = 0.0015). Haplotype [C; A; C] was associated with a significant increase in fat percentage when compared with the other haplotypes.

Numerous studies have found a large number of QTL associated with productive and functional traits in the cattle genome. Several Countries have already established research programs aiming at identification and exploitation of QTLs in dairy cattle (Bovenhuis and Schrooten 2002). However, the published results are not directly exploitable for MAS, as commercial interests limit the information in the literature, and because marker-QTL phase is specific to each family.

A selective DNA pooling approach using milk samples was employed to map QTL affecting milk somatic cells count (MSCC) in the Italian Brown Swiss dairy cattle population. The mapping population consisted of five half sib daughter families of Brown Swiss bulls, sires of 1000 to 3600 daughters. Two hundred highest and 200 lowest daughters, ranked by dam-corrected EBV, were selected from the high and the low tail. Four independent replicate pools, each made of 50 randomly chosen daughters, were prepared for each sire-tail combination. Dinucleotide microsatellite markers were used to scan the genome. Sire marker allele frequencies were estimated by densitometry and shadow correction analysis. Significance threshold of 10% aFDR was used at the marker level, and resulted in a critical CWER P-value of 0.054.A threshold of 20% aFDR within the significant markers was used at the sire-marker level and resulted in a critical P-values of 0.058. Out of 145 markers, 41 were significant. Out of 122 sire-marker tests, at the significant markers, 58 resulted significant. QTL regions will be selected for further intensive study. This is the first complete genome scan for MSCC in the Brown Swiss breed.

Milk contains a number of micro-components (“functional foods”) having specific health promoting activities. Important among these compounds are the conjugated linoleic acids (CLAs) consistingof a mixture of isomers of C18:2 fatty acid with conjugated doublebonds. The predominant form in milk fat is the cis-9,trans-11 isomer, accounting for 75 to 90% of the total CLA content. Biomedical studies with animal models have demonstrated a variety of preventive health effects from CLA including anticarcinogenic, antiatherogenic, antiobesity, immune system enhancement and antidiabetic benefits. Linolenic acid does not form CLA as an intermediate in rumen biohydrogenation, but it does form vaccenic acid (C18:1,t11) that enters the bloodstream from the rumen, through the abomasum and intestine, and is transported to the mammary gland. Several investigations using direct and indirect approaches established that endogenous synthesis of CLA via ∆9-desaturase action on vaccenic acid (VA, C18:1,t11) is the predo...