F. Schiavini

Quantitative trait loci affecting milk yield and protein percentage in a three-country Brown Swiss population.

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.

DGAT1 p.K232A polymorphism in dairy and dual purpose Italian cattle breeds

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 (<1%). In Italian Brown and Modenese, this allele was not detected. These results indicated that this polymorphic site can be considered for association studies only in Italian Holstein and Reggiana breeds. Deviation from Hardy- Weinberg equilibrium was observed in the Reggiana breed (P<0.01) in which there was an excess of heterozygous sires and absence of animals with the p.232KK genotype. This result should be further evaluated because the analysed sires represented almost all bulls available for artificial insemination in this breed. Comparison of allele frequencies at the DGAT1 locus with several other Holstein populations showed a wide range of variability, probably due to different selection strategies adopted.

Extensive Long-Range and Nonsyntenic Linkage Disequilibrium in Livestock Populations: Deconstruction of a Conundrum

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.

Genome-wide association study for somatic cell score in Valdostana Red Pied cattle breed using pooled DNA

BackgroundMastitis 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.ResultsA 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.ConclusionsThe 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.

A whole genome scan for QTL affecting milk protein percentage in Italian Holstein cattle, applying selective milk DNA pooling and multiple marker mapping in a daughter design

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.

Genomic and genetic variability of six chicken populations using single nucleotide polymorphism and copy number variants as markers

Genomic and genetic variation among six Italian chicken native breeds (Livornese, Mericanel della Brianza, Milanino, Bionda Piemontese, Bianca di Saluzzo and Siciliana) were studied using single nucleotide polymorphism (SNP) and copy number variants (CNV) as markers. A total of 94 DNA samples genotyped with Axiom® Genome-Wide Chicken Genotyping Array (Affymetrix) were used in the analyses. The results showed the genetic and genomic variability occurring among the six Italian chicken breeds. The genetic relationship among animals was established with a principal component analysis. The genetic diversity within breeds was calculated using heterozygosity values (expected and observed) and with Wrights F-statistics. The individual-based CNV calling, based on log R ratio and B-allele frequency values, was done by the Hidden-Markov Model (HMM) of PennCNV software on autosomes. A hierarchical agglomerative clustering was applied in each population according to the absence or presence of definite CNV regions (CNV were grouped by overlapping of at least 1 bp). The CNV map was built on a total of 1003 CNV found in individual samples, after grouping by overlaps, resulting in 564 unique CNV regions (344 gains, 213 losses and 7 complex), for a total of 9.43 Mb of sequence and 1.03% of the chicken assembly autosome. All the approaches using SNP data showed that the Siciliana breed clearly differentiate from other populations, the Livornese breed separates into two distinct groups according to the feather colour (i.e. white and black) and the Bionda Piemontese and Bianca di Saluzzo breeds are closely related. The genetic variability found using SNP is comparable with that found by other authors in the same breeds using microsatellite markers. The CNV markers analysis clearly confirmed the SNP results.

Quantitative trait loci mapping for conjugated linoleic acid, vaccenic acid and ∆9-desaturase in Italian Brown Swiss dairy cattle using selective DNA pooling

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.

Identification and validation of copy number variants in Italian Brown Swiss dairy cattle using Illumina Bovine SNP50 Beadchip

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.

Detection of QTL for milk protein percentage in Italian Friesian cattle by AFLP markers and selective genotyping

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 six Bos 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.

The Use of Kosher Phenotyping for Mapping QTL Affecting Susceptibility to Bovine Respiratory Disease

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.