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 Δ⁹‐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 Δ⁹‐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.     

Contributions of FASN and SCD gene polymorphisms on fatty acid composition in muscle from Japanese Black cattle

The fatty acid synthase (FASN) and stearoyl‐CoA desaturase (delta‐9‐desaturase) (SCD) genes affect fatty acid composition. This study evaluated the contributions of polymorphisms of these genes on fatty acid composition in muscle in two different populations: 1189 and 1058 Japanese Black cattle from the Miyagi and the Yamagata populations respectively. We sampled intramuscular fat from the longissimus thoracis muscle in the Miyagi population and from the trapezius muscle in the Yamagata population. The collective contributions of FASN and SCD polymorphisms to total additive genetic variance for oleic acid were 13.46% in the Miyagi population and 16.29% in the Yamagata population and to phenotypic variance were 5.45% and 6.54% respectively. Although the individual effects of FASN and SCD polymorphisms on fatty acid composition were small, overall gene substitution may effectively improve fatty acid composition. In addition, we found that gene polymorphism contributions of fatty acids varied by population even in the same breed.     

Fatty acid synthase effects on bovine adipose fat and milk fat

A quantitative trait locus (QTL) was identified by linkage analysis on bovine Chromosome 19 that affects the fatty acid, myristic acid (C14:0), in subcutaneous adipose tissue of pasture-fed beef cattle (99% level: experiment-wise significance). The QTL was also shown to have significant effects on ten fatty acids in the milk fat of pasture-fed dairy cattle. A positional candidate gene for this QTL was identified as fatty acid synthase (FASN), which is a multifunctional enzyme with a central role in the metabolism of lipids. Five single nucleotide polymorphisms (SNPs) were identified in the bovine FASN gene, and animals were genotyped for FASN SNPs in three different cattle resource populations. Linkage and association mapping results using these SNPs were consistent with FASN being the gene underlying the QTL. SNP substitution effects for C14:0 percentage were found to have an effect in the opposite direction in adipose fat to that in milk fat. It is concluded that SNPs in the bovine FASN gene are associated with variation in the fatty acid composition of adipose fat and milk fat.     

A Global View of 54,001 Single Nucleotide Polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and Their Transferability to Water Buffalo

The Illumina BovineSNP50 BeadChip features 54,001 informative single nucleotide polymorphisms (SNPs) that uniformly span the entire bovine genome. Among them, 52,255 SNPs have locations assigned in the current genome assembly (Btau_4.0), including 19,294 (37%) intragenic SNPs (i.e., located within genes) and 32,961 (63%) intergenic SNPs (i.e., located between genes). While the SNPs represented on the Illumina Bovine50K BeadChip are evenly distributed along each bovine chromosome, there are over 14,000 genes that have no SNPs placed on the current BeadChip. Kernel density estimation, a non-parametric method, was used in the present study to identify SNP-poor and SNP-rich regions on each bovine chromosome. With bandwidth = 0.05 Mb, we observed that most regions have SNP densities within 2 standard deviations of the chromosome SNP density mean. The SNP density on chromosome X was the most dynamic, with more than 30 SNP-rich regions and at least 20 regions with no SNPs. Genotyping ten water buffalo using the Illumina BovineSNP50 BeadChip revealed that 41,870 of the 54,001 SNPs are fully scored on all ten water buffalo, but 6,771 SNPs are partially scored on one to nine animals. Both fully scored and partially/no scored SNPs are clearly clustered with various sizes on each chromosome. However, among 43,687 bovine SNPs that were successfully genotyped on nine and ten water buffalo, only 1,159 were polymorphic in the species. These results indicate that the SNPs sites, but not the polymorphisms, are conserved between two species. Overall, our present study provides a solid foundation to further characterize the SNP evolutionary process, thus improving understanding of within- and between-species biodiversity, phylogenetics and adaption to environmental changes.     

A global view of 54,001 single nucleotide polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and their transferability to water buffalo

The Illumina BovineSNP50 BeadChip features 54,001 informative single nucleotide polymorphisms (SNPs) that uniformly span the entire bovine genome. Among them, 52,255 SNPs have locations assigned in the current genome assembly (Btau_4.0), including 19,294 (37%) intragenic SNPs (i.e., located within genes) and 32,961 (63%) intergenic SNPs (i.e., located between genes). While the SNPs represented on the Illumina Bovine50K BeadChip are evenly distributed along each bovine chromosome, there are over 14,000 genes that have no SNPs placed on the current BeadChip. Kernel density estimation, a non-parametric method, was used in the present study to identify SNP-poor and SNP-rich regions on each bovine chromosome. With bandwidth = 0.05 Mb, we observed that most regions have SNP densities within 2 standard deviations of the chromosome SNP density mean. The SNP density on chromosome X was the most dynamic, with more than 30 SNP-rich regions and at least 20 regions with no SNPs. Genotyping ten water buffalo using the Illumina BovineSNP50 BeadChip revealed that 41,870 of the 54,001 SNPs are fully scored on all ten water buffalo, but 6,771 SNPs are partially scored on one to nine animals. Both fully scored and partially/no scored SNPs are clearly clustered with various sizes on each chromosome. However, among 43,687 bovine SNPs that were successfully genotyped on nine and ten water buffalo, only 1,159 were polymorphic in the species. These results indicate that the SNPs sites, but not the polymorphisms, are conserved between two species. Overall, our present study provides a solid foundation to further characterize the SNP evolutionary process, thus improving understanding of within- and between-species biodiversity, phylogenetics and adaption to environmental changes.     

A region on BTA14 that includes the positional candidate genes LYPLA1, XKR4 and TMEM68 is associated with feed intake and growth phenotypes in cattle(1)

Feed cost for beef cattle is the largest expense incurred by cattle producers. The development of genetic markers to enhance selection of more efficient animals that require less feed while still achieving acceptable levels of production has the potential to substantially reduce production costs. A genome‐wide marker association approach based on the Illumina BovineSNP50 BeadChip? was used to identify genomic regions affecting average daily feed intake (ADFI), average daily gain (ADG) and residual feed intake traits in a population of 1159 crossbred steers. This approach identified a region on BTA14 from 22.02 to 23.92 Mb containing several single‐nucleotide polymorphisms (SNPs) that have significant association with at least one of the traits. Two genes in this region, lysophospholipase 1 (LYPLA1) and transmembrane protein 68 (TMEM68), appeared to be logical positional and functional candidate genes. LYPLA1 deacylates ghrelin, a hormone involved in the regulation of appetite in the rat stomach, while TMEM68 is expressed in bovine rumen, abomasum, intestine and adipose tissue in cattle, and likely affects lipid biosynthetic processes. SNPs lying in or near these two genes were identified by sequencing a subset of animals with extreme phenotypes. A total of 55 SNPs were genotyped and tested for association with the same population of steers. After correction for multiple testing, five markers within 22.79–22.84 Mb, located downstream of TMEM68, and between TMEM68 and the neighbouring gene XKR4, were significant for both ADFI and ADG. Genetic markers predictive of feed intake and weight gain phenotypes in this population of cattle may be useful for the identification and selection of animals that consume less feed, although further evaluation of these markers for effects on other production traits and validation in additional populations will be required.     

Genome Wide Association Study Identifies 20 Novel Promising Genes Associated with Milk Fatty Acid Traits in Chinese Holstein

Detecting genes associated with milk fat composition could provide valuable insights into the complex genetic networks of genes underling variation in fatty acids synthesis and point towards opportunities for changing milk fat composition via selective breeding. In this study, we conducted a genome-wide association study (GWAS) for 22 milk fatty acids in 784 Chinese Holstein cows with the PLINK software. Genotypes were obtained with the Illumina BovineSNP50 Bead chip and a total of 40,604 informative, high-quality single nucleotide polymorphisms (SNPs) were used. Totally, 83 genome-wide significant SNPs and 314 suggestive significant SNPs associated with 18 milk fatty acid traits were detected. Chromosome regions that affect milk fatty acid traits were mainly observed on BTA1, 2, 5, 6, 7, 9, 13, 14, 18, 19, 20, 21, 23, 26 and 27. Of these, 146 SNPs were associated with more than one milk fatty acid trait; most of studied fatty acid traits were significant associated with multiple SNPs, especially C18:0 (105 SNPs), C18 index (93 SNPs), and C14 index (84 SNPs); Several SNPs are close to or within the DGAT1, SCD1 and FASN genes which are well-known to affect milk composition traits of dairy cattle. Combined with the previously reported QTL regions and the biological functions of the genes, 20 novel promising candidates for C10:0, C12:0, C14:0, C14:1, C14 index, C18:0, C18:1n9c, C18 index, SFA, UFA and SFA/UFA were found, which composed of HTR1B, CPM, PRKG1, MINPP1, LIPJ, LIPK, EHHADH, MOGAT1, ECHS1, STAT1, SORBS1, NFKB2, AGPAT3, CHUK, OSBPL8, PRLR, IGF1R, ACSL3, GHR and OXCT1. Our findings provide a groundwork for unraveling the key genes and causal mutations affecting milk fatty acid traits in dairy cattle.     

Genome‐wide association study for semen traits of the bulls in Chinese Holstein

A genome‐wide association study (GWAS) was performed to identify markers and candidate genes for five semen traits in the Holstein bull population in China. The analyzed dataset consisted of records from 692 bulls from eight bull stations; each bull was genotyped using the Illumina BovineSNP50 BeadChip. Association tests between each trait and the 41 188 informative high‐quality SNPs were achieved with gapit software. In total, 19 suggestive significant SNPs, partly located within the reported QTL regions or within or close to the reported candidate genes, associated with five semen traits were detected. By combining our GWAS results with the biological functions of these genes, eight novel promising candidate genes, including ETNK1, PDE3A, PDGFRB, CSF1R, WT1, DSCAML1, SOD1 and RUNX2, were identified that potentially relate to semen traits. Our findings may provide a basis for further research on the genetic mechanism of semen traits and marker‐assisted selection of such traits in Holstein bulls.     

Design of a bovine low-density SNP array optimized for imputation

The Illumina BovineLD BeadChip was designed to support imputation to higher density genotypes in dairy and beef breeds by including single-nucleotide polymorphisms (SNPs) that had a high minor allele frequency as well as uniform spacing across the genome except at the ends of the chromosome where densities were increased. The chip also includes SNPs on the Y chromosome and mitochondrial DNA loci that are useful for determining subspecies classification and certain paternal and maternal breed lineages. The total number of SNPs was 6,909. Accuracy of imputation to Illumina BovineSNP50 genotypes using the BovineLD chip was over 97% for most dairy and beef populations. The BovineLD imputations were about 3 percentage points more accurate than those from the Illumina GoldenGate Bovine3K BeadChip across multiple populations. The improvement was greatest when neither parent was genotyped. The minor allele frequencies were similar across taurine beef and dairy breeds as was the proportion of SNPs that were polymorphic. The new BovineLD chip should facilitate low-cost genomic selection in taurine beef and dairy cattle.     

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.     

A genome‐wide association study suggests several novel candidate genes for carcass traits in Chinese Simmental beef cattle

A genome‐wide association study (GWAS) was conducted for two carcass traits in Chinese Simmental beef cattle. The experimental population consisted of 1301 individuals genotyped with the Illumina BovineHD SNP BeadChip (770K). After quality control, 671 990 SNPs and 1217 individuals were retained for the GWAS. The phenotypic traits included carcass weight and bone weight, which were measured after the cattle were slaughtered at 16 to 18 months of age. Three statistical models—a fixed polygene model, a random polygene model and a composite interval mapping polygene model—were used for the GWAS. The genome‐wide significance threshold after Bonferroni correction was 7.44E‐08 (= 0.05/671 990). In this study, we detected eight and seven SNPs significantly associated with carcass weight and bone weight respectively. In total, 11 candidate genes were identified within or close to these significant SNPs. Of these, we found several novel candidate genes, including PBX1, GCNT4, ALDH1A2, LCORL and WDFY3, to be associated with carcass weight and bone weight in Chinese Simmental beef cattle, and their functional roles need to be verified in further studies.     

Whole‐genome prediction of fatty acid composition in meat of Japanese Black cattle

Because fatty acid composition influences the flavor and texture of meat, controlling it is particularly important for cattle breeds such as the Japanese Black, characterized by high meat quality. We evaluated the predictive ability of single‐step genomic best linear unbiased prediction (ssGBLUP) in fatty acid composition of Japanese Black cattle by assessing the composition of seven fatty acids in 3088 cattle, of which 952 had genome‐wide marker genotypes. All sires of the genotyped animals were genotyped, but their dams were not. Cross‐validation was conducted for the 952 animals. The prediction accuracy was higher with ssGBLUP than with best linear unbiased prediction (BLUP) for all traits, and in an empirical investigation, the gain in accuracy of using ssGBLUP over BLUP increased as the deviations in phenotypic values of the animals increased. In addition, the superior accuracy of ssGBLUP tended to be more evident in animals whose maternal grandsire was genotyped than in other animals, although the effect was small.     

Association analyses of single nucleotide polymorphisms in bovine stearoyl-CoA desaturase and fatty acid synthase genes with fatty acid composition in commercial cross-bred beef steers

Two previously reported non‐synonymous coding single nucleotide polymorphisms (SNPs) of bovine stearoyl‐CoA desaturase (delta‐9‐desaturase) (SCD) (c.878C>T) and fatty acid synthase (FASN) (g:17924A>G) were assessed for their associations with 72 individual and 12 groups of fatty acids in brisket adipose tissue of 223 Canadian commercial cross‐bred beef steers. It was found that the ‘CC’ genotype of the SCD SNP was significantly associated with lower concentrations of saturated fatty acids (SFA) including 10:0, 14:0 and 20:0, higher concentrations of monounsaturated fatty acids including 9c‐14:1, 12c‐16:1 and 13c‐18:1, higher concentrations of polyunsaturated fatty acids (PUFA) including 9c,15c‐18:2, 10c,12c‐18:2, 11c,13t‐18:2 and 12c,14t‐18:2, but lower concentrations of other PUFA of 9c,13t/8t,12c and 20:2n‐6 (P < 0.05). The ‘AA’ genotype of the FASN SNP was significantly associated with higher concentrations of SFAs of 10:0, 12:0, 13:0, 14:0 and 15:0, lower concentrations of unsaturated fatty acids of 9c‐18:1 and 20:3n‐6, and higher concentrations of unsaturated fatty acids of 9c‐14:1 and 12c‐16:1 (P < 0.05). Significant epistatic effects between the SCD and FASN SNP genotypes were also found for several fatty acids including 10:0, 23:0, 6t/7t/8t‐18:1, 12t‐18:1, 13t/14t‐18:1, 16t‐18:1, total trans18:1 and 9c,13t/8t,12c‐18:2 (P < 0.05). These results further suggest that SCD and FASN are strong candidate genes influencing fatty acid composition in beef cattle.     

基于高密度SNP标记的肉牛人工选择痕迹筛查

近年来随着遗传改良工作的实施,人工选择大大提高了肉牛的生产性能并使其遗传基础发生巨大改变。文章基于Illumina BovineSNP50(54K)和BovineHD(770K)两款芯片数据,采用FST检验方法分析牛群的遗传分化,并筛查人工选择在牛的基因组留下的印记。通过全基因组范围内的扫描,共发现47 104个"离群"位点和3064个群体特异的人工选择"候选基因",如CLIC5、TG、CACNA2D1、FSHR等。通过基因注释对基因的生物学过程和分子功能进行富集分析。文章构建了我国肉牛的全基因组的选择信号图谱,为深入研究人工选择和理解生物进化提供线索,且研究结果也显示人工选择对基因组的影响在牛品种遗传改良中发挥了重要作用。     

Mapping of Calf Death in Japanese Black Cattle

Weak calf syndrome (WCS) is a major cause of calf death in Japanese Black cattle. Among IARS disorders, the isoleucyl-tRNA synthetase c.235G>C mutation has been identified as one of the causes of WCS. However, calf deaths differing from those attributed to IARS disorder has been occurring. To identify other genes potentially responsible for these calf deaths, we constructed three populations of three bulls (Bull-1, -2 and -3) that did not carry the IARS mutation, and dead calves (18, 28, and 31 calves) and healthy cattle (18, 15, and 10 cattle) sired by these bulls. The populations were genotyped using the BovineSNP50 BeadChip, but homozygosity mapping did not detect any associated genomic regions with calf death. Linkage analysis performed using each population as a paternal half-sib family of Bull-1, Bull-2, and Bull-3 revealed that, in the Bull-1 population, calf death was mapped to the 8.94?Mb–14.53?Mb and 29.82?Mb–33.77?Mb regions of BTA29. The findings suggested that the incidence of calf death in calves sired by Bull-1 was a hereditary disease exhibiting a dominant, not recessive, inheritance pattern.     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

Accuracy of imputation to whole-genome sequence data in Holstein Friesian cattle

Background

The use of whole-genome sequence data can lead to higher accuracy in genome-wide association studies and genomic predictions. However, to benefit from whole-genome sequence data, a large dataset of sequenced individuals is needed. Imputation from SNP panels, such as the Illumina BovineSNP50 BeadChip and Illumina BovineHD BeadChip, to whole-genome sequence data is an attractive and less expensive approach to obtain whole-genome sequence genotypes for a large number of individuals than sequencing all individuals. Our objective was to investigate accuracy of imputation from lower density SNP panels to whole-genome sequence data in a typical dataset for cattle.

Methods

Whole-genome sequence data of chromosome 1 (1737 471 SNPs) for 114 Holstein Friesian bulls were used. Beagle software was used for imputation from the BovineSNP50 (3132 SNPs) and BovineHD (40 492 SNPs) beadchips. Accuracy was calculated as the correlation between observed and imputed genotypes and assessed by five-fold cross-validation. Three scenarios S40, S60 and S80 with respectively 40%, 60%, and 80% of the individuals as reference individuals were investigated.

Results

Mean accuracies of imputation per SNP from the BovineHD panel to sequence data and from the BovineSNP50 panel to sequence data for scenarios S40 and S80 ranged from 0.77 to 0.83 and from 0.37 to 0.46, respectively. Stepwise imputation from the BovineSNP50 to BovineHD panel and then to sequence data for scenario S40 improved accuracy per SNP to 0.65 but it varied considerably between SNPs.

Conclusions

Accuracy of imputation to whole-genome sequence data was generally high for imputation from the BovineHD beadchip, but was low from the BovineSNP50 beadchip. Stepwise imputation from the BovineSNP50 to the BovineHD beadchip and then to sequence data substantially improved accuracy of imputation. SNPs with a low minor allele frequency were more difficult to impute correctly and the reliability of imputation varied more. Linkage disequilibrium between an imputed SNP and the SNP on the lower density panel, minor allele frequency of the imputed SNP and size of the reference group affected imputation reliability.     

Investigation of transferability of BovineSNP50 BeadChip from cattle to water buffalo for genome wide association study

Cattle and water buffalo belong to the same subfamily Bovinae and share chromosome banding and gene order homology. In this study, we used genome-wide Illumina BovineSNP50 BeadChip to analyze 91 DNA samples from three breeds of water buffalo (Nili-Ravi, Murrah and their crossbred with local GuangXi buffalos in China), to demonstrate the genetic divergence between cattle and water buffalo through a large single nucleotide polymorphism (SNP) transferability study at the whole genome level, and performed association analysis of functional traits in water buffalo as well. A total of 40,766 (75.5 %) bovine SNPs were found in the water buffalo genome, but 49,936 (92.5 %) were with only one allele, and finally 935 were identified to be polymorphic and useful for association analysis in water buffalo. Therefore, the genome sequences of water buffalo and cattle shared a high level of homology but the polymorphic status of the bovine SNPs varied between these two species. The different patterns of mutations between species may associate with their phenotypic divergence due to genome evolution. Among 935 bovine SNPs, we identified a total of 9 and 7 SNPs significantly associated to fertility and milk production traits in water buffalo, respectively. However, more works in larger sample size are needed in future to verify these candidate SNPs for water buffalo.     

Genome‐wide association study for pigmentation traits in Chinese Holstein population

With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.     

Fine Mapping for Weaver Syndrome in Brown Swiss Cattle and the Identification of 41 Concordant Mutations across NRCAM,PNPLA8 and CTTNBP2

Bovine Progressive Degenerative Myeloencephalopathy (Weaver Syndrome) is a recessive neurological disease that has been observed in the Brown Swiss cattle breed since the 1970’s in North America and Europe. Bilateral hind leg weakness and ataxia appear in afflicted animals at 6 to 18 months of age, and slowly progresses to total loss of hind limb control by 3 to 4 years of age. While Weaver has previously been mapped to Bos taurus autosome (BTA) 4∶46–56 Mb and a diagnostic test based on the 6 microsatellite (MS) markers is commercially available, neither the causative gene nor mutation has been identified; therefore misdiagnosis can occur due to recombination between the diagnostic MS markers and the causative mutation. Analysis of 34,980 BTA 4 SNPs genotypes derived from the Illumina BovineHD assay for 20 Brown Swiss Weaver carriers and 49 homozygous normal bulls refined the Weaver locus to 48–53 Mb. Genotyping of 153 SNPs, identified from whole genome sequencing of 10 normal and 10 carrier animals, across a validation set of 841 animals resulted in the identification of 41 diagnostic SNPs that were concordant with the disease. Except for one intergenic SNP all are associated with genes expressed in nervous tissues: 37 distal to NRCAM, one non-synonymous (serine to asparagine) in PNPLA8, one synonymous and one non-synonymous (lysine to glutamic acid) in CTTNBP2. Haplotype and imputation analyses of 7,458 Brown Swiss animals with Illumina BovineSNP50 data and the 41 diagnostic SNPs resulted in the identification of only one haplotype concordant with the Weaver phenotype. Use of this haplotype and the diagnostic SNPs more accurately identifies Weaver carriers in both Brown Swiss purebred and influenced herds.