Performance of using opposing homozygotes for paternity testing in Japanese Black cattle

Genome-wide single nucleotide polymorphism (SNP) markers in Japanese Black cattle enable genomic prediction and verifying parent–offspring relationships. We assessed the performance of opposing homozygotes (OH) for paternity testing in Japanese Black cattle, using SNP genotype information of 50 sires and 3,420 fattened animals, 1,945 of which were fathered by the 50 genotyped sires. The number of OH was counted for each sire–progeny pair in 28,764 SNPs with minor allele frequencies of ≥0.05 in this population. Across all pairs of animals, the number of OH tended to increase as the pedigree-based coefficient of relationship decreased. With a threshold of 288 (1% of SNPs) for paternity testing, most sire–progeny pairs were detected as true relationships. The frequency of Mendelian inconsistencies was 2.4%, reflecting the high accuracy of pedigree information in Japanese Black cattle population. The results indicate the utility of OH for paternity testing in Japanese Black cattle.     

Genetic diversity and genomic inbreeding in Japanese Black cows in the islands of Okinawa Prefecture evaluated using single-nucleotide polymorphism array

Maintaining genetic diversity and inbreeding control are important in Japanese Black cattle production, especially in remote areas such as the islands of Okinawa Prefecture. Using a single-nucleotide polymorphism (SNP) array, we evaluated the genetic diversity and genomic inbreeding in Japanese Black cows from the islands of Okinawa Prefecture and compared them to those from other locations across Japan. Linkage disequilibrium decay was slower in cows in the islands of Okinawa Prefecture. The estimated effective population size declined over time in both populations. The genomic inbreeding coefficient (F_ROH) was estimated using long stretches of consecutive homozygous SNPs (runs of homozygosity; ROH). F_ROH was higher in the cows on the islands of Okinawa Prefecture than on other locations. In total, 818 ROH fragments, including those containing NCAPG and PLAG1, which are major quantitative trait loci for carcass weight in Japanese Black cattle, were present at significantly higher frequencies in cows in the islands of Okinawa Prefecture. This suggests that the ROH fragments are under strong selection and that cows in the islands of Okinawa Prefecture have low genetic diversity and high genomic inbreeding relative to those at other locations. SNP arrays are useful tools for evaluating genetic diversity and genomic inbreeding in cattle.     

Estimation of variance and genomic prediction using genotypes imputed from low‐density marker subsets for carcass traits in Japanese black cattle

The influence of genotype imputation using low‐density single nucleotide polymorphism (SNP) marker subsets on the genomic relationship matrix (G matrix), genetic variance explained, and genomic prediction (GP) was investigated for carcass weight and marbling score in Japanese Black fattened steers, using genotype data of approximately 40,000 SNPs. Genotypes were imputed using equally spaced SNP subsets of different densities. Two different linear models were used. The first (model 1) incorporated one G matrix, while the second (model 2) used two different G matrices constructed using the selected and remaining SNPs. When using model 1, the estimated additive genetic variance was always larger when using all SNPs obtained via genotype imputation than when using only equally spaced SNP subsets. The correlations between the genomic estimated breeding values obtained using genotype imputation with at least 3,000 SNPs and those using all available SNPs without imputation were higher than 0.99 for both traits. While additive genetic variance was likely to be partitioned with model 2, it did not enhance the accuracy of GP compared with model 1. These results indicate that genotype imputation using an equally spaced low‐density panel of an appropriate size can be used to produce a cost‐effective, valid GP.     

Estimation of variance components for carcass traits in Japanese Black cattle using 50K SNP genotype data

Genomic selection using high‐density single nucleotide polymorphism (SNP) genotype data may accelerate genetic improvements in livestock animals. In this study, we attempted to estimate the variance components of six carcass traits in fattened Japanese Black steers using SNP genotype data. Six hundred and seventy‐three steers were genotyped using an Illumina Bovine SNP50 BeadChip and phenotyped for cold carcass weight, ribeye area, rib thickness, subcutaneous fat thickness, estimated yield percent and marbling score. Additive polygenic variance and the variance attributable to a set of SNPs that had statistically significant effects on the trait were estimated via Gibbs sampling with two models: (i) a model with the chosen SNPs and the additive polygenic effects; and (ii) a model with the polygenic effects alone. The proportion of the estimated variance attributable to the SNPs became higher as the number of SNP effects that fit increased. High correlations between breeding values estimated with the model containing the polygenic effect alone and those estimated by chosen SNPs were obtained. No fraction of the total genetic variance was explained by SNPs associated with the trait at P ≥ 0.1. Our results suggest that for the carcass traits of Japanese Black cattle, a maximum of half of the total additive genetic variance may be explained by SNPs between 100 several tens to several 100s.     

Accuracy of imputation of single nucleotide polymorphism marker genotypes from low‐density panels in Japanese Black cattle

Using target and reference fattened steer populations, the performance of genotype imputation using lower‐density marker panels in Japanese Black cattle was evaluated. Population imputation was performed using BEAGLE software. Genotype information for approximately 40 000 single nucleotide polymorphism (SNP) markers by Illumina BovineSNP50 BeadChip was available, and imputation accuracy was assessed based on the average concordance rates of the genotypes, varying equally spaced SNP densities, and the number of individuals in the reference population. Two additional statistics were also calculated as indicators of imputation performance. The concordance rates tended to be lower for SNPs with greater minor allele frequencies, or those located near the ends of the chromosomes. Longer autosomes yielded greater imputation accuracies than shorter ones. When SNPs were selected based on linkage disequilibrium information, relative imputation accuracy was slightly improved. When 3000 and 10 000 equally spaced SNPs were used, the imputation accuracies were greater than 90% and approximately 97%, respectively. These results indicate that combining genotyping using a lower‐density SNP chip with genotype imputation based on a population of individuals genotyped using a higher‐density SNP chip is a cost‐effective and valid approach for genomic prediction.     

A genome‐wide association study for fat‐related traits computed by image analysis in Japanese Black cattle

The objective of this study was to identify genomic regions associated with fat‐related traits using a Japanese Black cattle population in Hyogo. From 1836 animals, those with high or low values were selected on the basis of corrected phenotype and then pooled into high and low groups (n = 100 each), respectively. DNA pool‐based genome‐wide association study (GWAS) was performed using Illumina BovineSNP50 BeadChip v2 with three replicate assays for each pooled sample. GWAS detected that two single nucleotide polymorphisms (SNPs) on BTA7 (ARS‐BFGL‐NGS‐35463 and Hapmap23838‐BTA‐163815) and one SNP on BTA12 (ARS‐BFGL‐NGS‐2915) significantly affected fat percentage (FAR). The significance of ARS‐BFGL‐NGS‐35463 on BTA7 was confirmed by individual genotyping in all pooled samples. Moreover, association analysis between SNP and FAR in 803 Japanese Black cattle revealed a significant effect of SNP on FAR. Thus, further investigation of these regions is required to identify FAR‐associated genes and mutations, which can lead to the development of DNA markers for marker‐assisted selection for the genetic improvement of beef quality.     

Genomic prediction for meat and carcass traits in Nellore cattle using a Markov blanket algorithm

This study was carried out to evaluate the advantage of preselecting SNP markers using Markov blanket algorithm regarding the accuracy of genomic prediction for carcass and meat quality traits in Nellore cattle. This study considered 3675, 3680, 3660 and 524 records of rib eye area (REA), back fat thickness (BF), rump fat (RF), and Warner–Bratzler shear force (WBSF), respectively, from the Nellore Brazil Breeding Program. The animals have been genotyped using low-density SNP panel (30 k), and subsequently imputed for arrays with 777 k SNPs. Four Bayesian specifications of genomic regression models, namely Bayes A, Bayes B, Bayes Cπ and Bayesian Ridge Regression methods were compared in terms of prediction accuracy using a five folds cross-validation. Prediction accuracy for REA, BF and RF was all similar using the Bayesian Alphabet models, ranging from 0.75 to 0.95. For WBSF, the predictive ability was higher using Bayes B (0.47) than other methods (0.39 to 0.42). Although the prediction accuracies using Markov blanket of SNP markers were lower than those using all SNPs, for WBSF the relative gain was lower than 13%. With a subset of informative SNPs markers, identified using Markov blanket, probably, is possible to capture a large proportion of the genetic variance for WBSF. The development of low-density and customized arrays using Markov blanket might be cost-effective to perform a genomic selection for this trait, increasing the number of evaluated animals, improving the management decisions based on genomic information and applying genomic selection on a large scale.     

Genome‐wide association studies identified variants for taurine concentration in Japanese Black beef

We performed genome‐wide association studies (GWAS) using the BovineSNP50 array to detect significant single nucleotide polymorphisms (SNPs) that may affect the concentration of 22 free amino acids and three peptides in Japanese Black beef cattle. A total of 574 Japanese Black cattle and 40,657 SNPs from the array were used for this study. Genome‐wide significant SNPs were detected for β‐alanine (three SNPs on chromosomes 22 and 29) and taurine (26 SNPs on chromosome 22). Importantly, the top two SNPs for taurine were highly significant (p = 6.2 × 10^?21), and the frequency of the increase‐concentration allele (Q) for taurine was found to be 0.73. The Q allele frequency of this population was similar to that of the other unrelated Japanese Black cattle, but different from that of the other breeds. In addition, the significant SNPs were not associated with carcass traits or fatty acid compositions. Interestingly, the top three of the four most significant SNPs for taurine were located near solute carrier family 6, member 6 (SLC6A6), which is a membrane transporter for taurine. We also found two associated variants in the 5′‐upstream region of SLC6A6; however, they were less significantly associated than the SNPs from the BovineSNP50 array.     

Novel SNP in 5' flanking region of EDG1 associated with marbling in Japanese Black beef cattle

Marbling, defined by the amount and distribution of intramuscular fat, is an economically important trait of beef cattle in Japan. The endothelial differentiation , sphingolipid G-protein-coupled receptor , 1 ( EDG1 ) gene has been considered as a positional functional candidate for the gene responsible for marbling. We have recently reported that 2 single nucleotide polymorphisms (SNPs), c.-312A > G in the 5' untranslated region (UTR) and c.*446G > A in the 3' UTR in EDG1 were associated with marbling in Japanese Black beef cattle, but this was not functional and a causal mutation for marbling. In the present study, we detected 2 novel SNPs, referred to as g.1475435G > A and g.1471620G > T , in the 5' flanking region of the EDG1 between low-marbled and high-marbled steer groups, which were previously shown to have EDG1 expression differences in musculus longissimus muscle. The g.1475435G > A SNP seemed not to segregate in Japanese Black beef cattle. The g.1471620G > T SNP was associated with the predicted breeding value for beef marbling standard number by the analyses using Japanese Black beef cattle population. Based on these findings, we hypothesized that the g.1471620G > T SNP might have an impact on EDG1 expression and also marbling.     

The T allele at the g.1471620G>T in the EDG1 gene associated with high marbling in Japanese Black cattle is at a low frequency in breeds not selected for marbling

Our previous study detected a single nucleotide polymorphism (SNP), g.1471620G > T , in the 5' flanking region of the endothelial differentiation sphingolipid G-protein-coupled receptor 1 ( EDG1 ) gene, which has been considered as a positional functional candidate for the gene responsible for marbling, and showed association of the g.1471620G > T SNP with marbling in Japanese Black beef cattle. In the present study, we investigated the allele frequency distribution of the g.1471620G > T SNP among the 5 cattle breeds, Japanese Black, Japanese Brown, Japanese Short Horn, Holstein, and Brown Swiss breeds. The T allele at the g.1471620G > T SNP associated with high marbling was found at high frequency in Japanese Black breed that has been subjected to a strong selection for high marbling, while the allele was absent or at very low frequencies in the other breeds that have not been strongly selected for high marbling. Based on this finding, we hypothesized that the pressure of the strong selection for high marbling in Japanese Black breed has increased the frequency of the T allele at the g.1471620G > T SNP in the EDG1 .     

Genome‐wide association study for fatty acid composition in Japanese Black cattle

Fatty acid composition is one of the important traits in beef. The aim of this study was to identify candidate genomic regions for fatty acid composition by genome‐wide association study with 50 K single nucleotide polymorphism (SNP) array in Japanese Black cattle. A total of 461 individuals and 40 657 SNPs were used in this study. We applied genome‐wide rapid association using mixed model and regression (GRAMMAR) and genomic control approaches to estimate the associations between genotypes and fatty acid composition. In addition, two SNPs in fatty acid synthase (FASN) (T1952A) and stearoyl‐CoA desaturase (SCD) (V293A) genes were also genotyped. Association analysis revealed that 30 significant SNPs for several fatty acids (C14:0, C14:1, C16:1 and C18:1) were located in the BTA19 FASN gene located within this region but the FASN mutation had no significant effect on any traits. We also detected one significant SNP for C18:1 on BTA23 and two SNPs for C16:0 on BTA25. The region around 17 Mb on BTA26 harbored two significant SNPs for C14:1 and SNP in SCD in this region showed the strongest association with C14:1. This study demonstrated novel candidate regions in BTA19, 23 and 25 for fatty acid composition.     

布莱凯特黑牛心脏脂肪酸结合蛋白基因单核苷酸多态性的研究

利用PCR-SSCP方法测定96头布莱凯特黑牛H-FABP基因外显子2和内含子2的部分序列多态性,并对含有多态序列的片段进行了序列分析。结果表明,引物1、3的扩增序列不存在多态性,引物2、4的扩增序列存在多态性;引物2在内含子第323位点处存在G→A突变,在群体中产生了AA、AB、BB 3种基因型,基因型频率分别为0.448、0.146、0.406,多态信息含量(PIC)为0.375,属于中度多态(0.25χ²检验该位点基因型频率不处于Hardy-Weinberg平衡状态;引物4在内含子2第872位点处存在C→G的突变,产生了HH、Hh、hh3种基因型,基因型频率分别为0.104、0.417、0.479,多态信息含量(PIC)为0.337,为中度多态(0.25χ²检验该位点基因型频率达到Hardy-Weinberg平衡状态。     

Development of SNP markers for individual identification and parentage test in a Japanese Black cattle population

This study describes the development of efficient single nucleotide polymorphism (SNP) markers for individual identification and parentage tests in a Japanese Black cattle population. An amplified fragment length polymorphism method was employed to detect informative candidate markers, and yielded 44 SNP markers from 220 primer combinations. 29 unlinked SNPs were finally selected as diagnostic markers. The allelic frequencies for each marker were estimated by using PCR‐RFLP in the Japanese Black population. Based on the frequency data, the estimated identity power of these markers was 2.73 × 10^?12. Parentage exclusion probabilities, when both suspected parents' genotypes were known and when only one suspected parent was genotyped, were estimated as 0.96929 and 0.99693, respectively. This panel of SNP markers is theoretically sufficient for individual identification, and would also be a powerful tool for a parentage test in Japanese Black cattle. The markers could contribute to the management of the beef industry in Japan.     

SNPpath: Characterizing cattle SNPs by enriched pathway terms

High‐density single nucleotide polymorphism (SNP) microarrays have made large‐scale genome‐wide association studies (GWAS) and genomic selection (GS) feasible. Valuable insight into the genetic basis underlying complex polygenic traits will likely be gained by considering functionally related sets of genes simultaneously. SNPpath, a suite of computer‐generated imagery‐based web servers has been developed to automatically annotate and characterize cattle SNPs by enriched KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway terms. The SNPpath allows users to navigate and analysis large SNP sets and is the only web server currently providing pathway annotations of cattle SNPs in National Center for Biotechnology Information's dbSNP database and three commercial platforms. Hence, we describe SNPpath and provide details of the query options, as well as biological examples of use. The SNPpath may be favorable for the analysis of combining SNP association analysis with pathway‐driven gene set enrichment analysis and is freely available at http://klab.sjtu.edu.cn/SNPpath .     

Genomic relationships computed from either next‐generation sequence or array SNP data

The use of sequence data in genomic prediction models is a topic of high interest, given the decreasing prices of current ‘next’‐generation sequencing technologies (NGS) and the theoretical possibility of directly interrogating the genomes for all causal mutations. Here, we compare by simulation how well genetic relationships (G) could be estimated using either NGS or ascertained SNP arrays. DNA sequences were simulated using the coalescence according to two scenarios: a ‘cattle’ scenario that consisted of a bottleneck followed by a split in two breeds without migration, and a ‘pig’ model where Chinese introgression into international pig breeds was simulated. We found that introgression results in a large amount of variability across the genome and between individuals, both in differentiation and in diversity. In general, NGS data allowed the most accurate estimates of G, provided enough sequencing depth was available, because shallow NGS (4×) may result in highly distorted estimates of G elements, especially if not standardized by allele frequency. However, high‐density genotyping can also result in accurate estimates of G . Given that genotyping is much less noisy than NGS data, it is suggested that specific high‐density arrays (~3M SNPs) that minimize the effects of ascertainment could be developed in the population of interest by sequencing the most influential animals and rely on those arrays for implementing genomic selection.     

KIT基因多态性与荷斯坦牛着色性状的关联分析

牛的被毛颜色是重要的外貌性状。本研究采集北京地区401头荷斯坦母牛血样及21头公牛冻精样品．提取基因组DNA,选取牛6号染色体上KIT基因的4个SNPs位点,利用飞行时间质谱技术检测基因型。通过数码相机照相并进行图像分析计算牛的被毛黑白比例,记录乳头颜色。将基因型和表型数据进行关联分析,结果显示KIT基因的G72792875T住点对被毛比例有显著影响（P〈0．05）,推测KIT基因可能与荷斯坦牛的着色性状的形成相关。     

Genetic and genomic analyses for predicted methane‐related traits in Japanese Black steers

The objectives of this study were to estimate genetic parameters and to perform a genome‐wide association study (GWAS) for predicted methane‐related traits in Japanese Black steers. The methane production and yield traits were predicted using on‐farm measurable traits, such as dry matter intake and average daily gain. A total of 4,578 Japanese Black steers, which were progenies of 362 sires genotyped with imputed 551,995 single nucleotide polymorphisms (SNPs), had phenotypes of predicted methane‐related traits during the total fattening period (52 weeks). For the estimation of genetic parameters, the estimated heritabilities were moderate (ranged from 0.57 to 0.60). In addition, the estimated genetic correlations of methane production traits with most of carcass traits and feed‐efficiency traits were unfavorable, but those of methane yield traits were favorable or low. For the GWAS, no genome‐wide significant SNP was detected, but a total of four quantitative trait locus (QTL) regions that explained more than 5.0% of genetic variance were localized on the genome, and some candidate genes associated with growth and feed‐efficiency traits were located on the regions. Our results suggest that the predicted methane‐related traits are heritable and some QTL regions for the traits are localized on the genome in Japanese Black steers.     

Comparison between haplotype-based and individual snp-based genomic predictions for beef fatty acid profile in Nelore cattle

The aim of this study was to evaluate the genomic predictions using the single-step genomic best linear unbiased predictor (ssGBLUP) method based on SNPs and haplotype markers associated with beef fatty acids (FAs) profile in Nelore cattle. The data set contained records from 963 Nelore bulls finished in feedlot (±90 days) and slaughtered with approximately 24 months of age. Meat samples from the Longissimus dorsi muscle were taken for FAs profile measurement. FAs were quantified by gas chromatography using a SP-2560 capillary column. Animals were genotyped with the high-density SNP panel (BovineHD BeadChip assay) containing 777,962 markers. SNPs with a minor allele frequency and a call rate lower than 0.05 and 0.90, respectively, monomorphic, located on sex chromosomes, and with unknown position were removed from the data set. After genomic quality control, a total of 469,981 SNPs and 892 samples were available for subsequent analyses. Missing genotypes were imputed and phased using the FImpute software. Haplotype blocks were defined based on linkage disequilibrium using the Haploview software. The model to estimate variance components and genetic parameters and to predict the genomic values included the random genetic additive effects, fixed effects of the contemporary group and the age at slaughter as a linear covariate. Accuracies using the haplotype-based approach ranged from 0.07 to 0.31, and those SNP-based ranged from 0.06 to 0.33. Regression coefficients ranged from 0.07 to 0.74 and from 0.08 to 1.45 using the haplotype- and SNP-based approaches, respectively. Despite the low to moderate accuracies for the genomic values, it is possible to obtain genetic progress trough selection using genomic information based either on SNPs or haplotype markers. The SNP-based approach allows less biased genomic evaluations, and it is more feasible when taking into account the computational and operational cost underlying the haplotypes inference.     

Identification of the putative ancestral allele of bovine single‐nucleotide polymorphisms

Identifying the action of natural selection from patterns of standing genetic variation has long been of interest to the population genetic community. Thanks to the availability of large single‐nucleotide polymorphism (SNP) data sets for many species and of high‐throughput SNP genotyping methods, whole‐genomic surveys to detect selective sweeps are now possible. Knowing the ancestral allele increases the power to detect selection. We present here a comparative genomic approach to determine the putative ancestral allele of bovine SNPs deposited in public databases. We analysed 19 551 488 SNPs and identified the putative ancestral allele for 14 339 107 SNPs. Our predicted ancestral alleles were in agreement with ancestral alleles detected by genotyping outgroup species for 97% SNPs from the BovineSNP50 BeadChip. This comparison indicates that our comparative genomic‐based approach to identify putative ancestral alleles is reliable.     

利用50K芯片解析和田羊和策勒黑羊的遗传规律

旨在探究和田羊和策勒黑羊遗传规律的不同,筛选优异基因及其分子标记,推动新疆南疆绵羊品种遗传资源的开发和利用.本研究分别选取健康、成年和田母羊84只、策勒黑羊41只,采集耳组织,放入75％的酒精.基于酚氯仿法提取DNA,经电泳和核酸仪检测,制备Illumina Ovine SNP50芯片.用Plink1.07对基因组数据...