Parameters affecting genome simulation for evaluating genomic selection method

The present study investigated the parameter settings for obtaining a simulated genome at steady state of allele frequency (mutation–drift equilibrium) and linkage disequilibrium (LD), and evaluated the impact of whether or not the simulated genome reached steady state of allele frequency and LD on the accuracy of genomic estimated breeding values (GEBVs). After 500 to 50 000 historical generations, the base population and subsequent seven generations were generated as recent populations. The allele frequency distribution of the last generations of the historical population and LD in the base population were calculated when varying the values of five parameters: initial minor allele frequency, mutation rate, effective population size, number of markers and chromosome length. The accuracies of GEBVs in the last generation of the recent population were calculated by genomic best linear unbiased prediction. The number of historical generations required to reach mutation–drift equilibrium depended on the initial allele frequency and mutation rate. Regardless of the parameters, LD reached a steady state before allele frequency distribution reached mutation–drift equilibrium. The accuracies of GEBVs largely reflect the extent of linkage disequilibrium with the exception of varying chromosome length, although there were no associations between the accuracies of GEBVs and allele frequency distribution.     

Power of a genome scan to detect and locate quantitative trait loci in cattle using dense single nucleotide polymorphisms

There is increasing use of dense single nucleotide polymorphisms (SNPs) for whole‐genome association studies (WGAS) in livestock to map and identify quantitative trait loci (QTL). These studies rely on linkage disequilibrium (LD) to detect an association between SNP genotypes and phenotypes. The power and precision of these WGAS are unknown, and will depend on the extent of LD in the experimental population. One complication for WGAS in livestock populations is that they typically consist of many paternal half‐sib families, and in some cases full‐sib families; unless this subtle population stratification is accounted for, many spurious associations may be reported. Our aim was to investigate the power, precision and false discovery rates of WGAS for QTL discovery, with a commercial SNP array, given existing patterns of LD in cattle. We also tested the efficiency of selective genotyping animals. A total of 365 cattle were genotyped for 9232 SNPs. We simulated a QTL effect as well as polygenic and environmental effects for all animals. One QTL was simulated on a randomly chosen SNP and accounted for 5%, 10% or 18% of the total variance. The power to detect a moderate‐sized additive QTL (5% of the phenotypic variance) with 365 animals genotyped was 37% (p < 0.001). Most importantly, if pedigree structure was not accounted for, the number of false positives significantly increased above those expected by chance alone. Selective genotyping also resulted in a significant increase in false positives, even when pedigree structure was accounted for.     

A whole genome scan to map QTL for milk production traits and somatic cell score in Canadian Holstein bulls

The detection and mapping of genetic markers linked to quantitative trait loci (QTL) can be utilized to enhance genetic improvement of livestock populations. With the completion of the bovine genome sequence assembly, single nucleotide polymorphisms (SNP) assays spanning the whole bovine genome and research work on large scale identification, validation and analysis of genotypic variation in cattle has become possible. The objective of the present study was to perform a whole genome scan to identify and map QTL affecting milk production traits and somatic cell scores using linkage disequilibrium (LD) regression and 1536 SNP markers. Three and 18 SNP were found to be associated with only milk yield (MY) at a genome and chromosome wise significance (p < 0.05) level respectively. Among the 21 significant SNP, 16 were in a region reported to have QTL for MY in other dairy cattle populations and while the rest five were new QTL finding. Four SNP out of 21 are significant for the milk production traits (MY, fat yield, protein yield (PY), and milk contents) in the present study. Six and nine SNP were associated with PY at a genome and chromosome wise significant (p < 0.05) level respectively. Three and 17 SNP were found to be associated with FY at a genome and chromosome wise significant (p < 0.05) level. Five and seven SNP were mapped with somatic cell score at a genome and chromosome wise significant (p < 0.05) level respectively. The results of this study have revealed QTL for MY, PY, protein percentage, FY, fat percentage, somatic cell score and persistency of milk in the Canadian dairy cattle population. The chromosome regions identified in this study should be further investigated to potentially identify the causative mutations underlying the QTL.     

Persistency of accuracy of genomic breeding values for different simulated pig breeding programs in developing countries

Genetic improvement of pigs in tropical developing countries has focused on imported exotic populations which have been subjected to intensive selection with attendant high population‐wide linkage disequilibrium (LD). Presently, indigenous pig population with limited selection and low LD are being considered for improvement. Given that the infrastructure for genetic improvement using the conventional BLUP selection methods are lacking, a genome‐wide selection (GS) program was proposed for developing countries. A simulation study was conducted to evaluate the option of using 60 K SNP panel and observed amount of LD in the exotic and indigenous pig populations. Several scenarios were evaluated including different size and structure of training and validation populations, different selection methods and long‐term accuracy of GS in different population/breeding structures and traits. The training set included previously selected exotic population, unselected indigenous population and their crossbreds. Traits studied included number born alive (NBA), average daily gain (ADG) and back fat thickness (BFT). The ridge regression method was used to train the prediction model. The results showed that accuracies of genomic breeding values (GBVs) in the range of 0.30 (NBA) to 0.86 (BFT) in the validation population are expected if high density marker panels are utilized. The GS method improved accuracy of breeding values better than pedigree‐based approach for traits with low heritability and in young animals with no performance data. Crossbred training population performed better than purebreds when validation was in populations with similar or a different structure as in the training set. Genome‐wide selection holds promise for genetic improvement of pigs in the tropics.     

主要家畜种基因组连锁图现状和展望

对主要家畜种基因组研究中用多态遗传标记,如微卫星标记、RFLP等,制作遗传连锁图现状进行了综述。在过去近半个世纪中,家畜基因组遗传连锁图的研究随着遗传标记技术,特别是DNA 微卫星等多态标记的迅速发展,取得了很大进展。其中主要介绍了主要家畜种牛、猪、绵羊、马和家禽基因组遗传连锁图发展的现状以及将来的发展趋势。     

Genomic best linear unbiased prediction method reflecting the degree of linkage disequilibrium

The degree of linkage disequilibrium (LD) between markers differs depending on the location of the genome; this difference biases genetic evaluation by genomic best linear unbiased prediction (GBLUP). To correct this bias, we used three GBLUP methods reflecting the degree of LD (GBLUP‐LD). In the three GBLUP‐LD methods, genomic relationship matrices were conducted from single nucleotide polymorphism markers weighted according to local LD levels. The predictive abilities of GBLUP‐LD were investigated by estimating variance components and assessing the accuracies of estimated breeding values using simulation data. When quantitative trait loci (QTL) were located at weak LD regions, the predictive abilities of the three GBLUP‐LD methods were superior to those of GBLUP and Bayesian lasso except when the number of QTL was small. In particular, the superiority of GBLUP‐LD increased with decreasing trait heritability. The rates of QTL at weak LD regions would increase when selection by GBLUP continues; this consequently decreases the predictive ability of GBLUP. Thus, the GBLUP‐LD could be applicable for populations selected by GBLUP for a long time. However, if QTL were located at strong LD regions, the accuracies of three GBLUP‐LD methods were lower than GBLUP and Bayesian lasso.     

Haplotype phasing after joint estimation of recombination and linkage disequilibrium in breeding populations

A novel method for haplotype phasing in families after joint estimation of recombination fraction and linkage disequilibrium is developed. Results from Monte Carlo computer simulations show that the newly developed E.M. algorithm is accurate if true recombination fraction is 0 even for single families of relatively small sizes. Estimates of recombination fraction and linkage disequilibrium were 0.00 （SD 0.00） and 0.19 （SD 0.03） for simulated recombination fraction and linkage disequilibrium of 0.00 and 0.20, respectively. A genome fragmentation phasing strategy was developed and used for phasing haplotypes in a sire and 36 progeny using the 50 k Illumina BeadChip by： a） estimation of the recombination fraction and LD in consecutive SNPs using family information, b） linkage analyses between fragments, c） phasing of haplotypes in parents and progeny and in following generations. Homozygous SNPs in progeny allowed determination of paternal fragment inheritance, and deduction of SNP sequence information of haplotypes from dams. The strategy also allowed detection of genotyping errors. A total of 613 recombination events were detected after linkage analysis was carried out between fragments. Hot and cold spots were identified at the individual （sire level）. SNPs for which the sire and calf were heterozygotes became informative （over 90%） after the phasing of haplotypes. Average of regions of identity between half-sibs when comparing its maternal inherited haplotypes （with at least 20 SNP） in common was 0.11 with a maximum of 0.29 and a minimum of 0.05. A Monte-Carlo simulation of BTA1 with the same linkage disequilibrium structure and genetic linkage as the cattle family yielded a 99.98 and 99.94% of correct phases for informative SNPs in sire and calves, respectively.     

Haplotype phasing after joint estimation of recombination and linkage disequilibrium in breeding populations

A novel method for haplotype phasing in families after joint estimation of recombination fraction and linkage disequilibrium is developed. Results from Monte Carlo computer simulations show that the newly developed E.M. algorithm is accurate if true recombination fraction is 0 even for single families of relatively small sizes. Estimates of recombination fraction and linkage disequilibrium were 0.00 (SD 0.00) and 0.19 (SD 0.03) for simulated recombination fraction and linkage disequilibrium of 0.00 and 0.20, respectively. A genome fragmentation phasing strategy was developed and used for phasing haplotypes in a sire and 36 progeny using the 50 k Illumina BeadChip by: a) estimation of the recombination fraction and LD in consecutive SNPs using family information, b) linkage analyses between fragments, c) phasing of haplotypes in parents and progeny and in following generations. Homozygous SNPs in progeny allowed determination of paternal fragment inheritance, and deduction of SNP sequence information of haplotypes from dams. The strategy also allowed detection of genotyping errors. A total of 613 recombination events were detected after linkage analysis was carried out between fragments. Hot and cold spots were identified at the individual (sire level). SNPs for which the sire and calf were heterozygotes became informative (over 90%) after the phasing of haplotypes. Average of regions of identity between half-sibs when comparing its maternal inherited haplotypes (with at least 20 SNP) in common was 0.11 with a maximum of 0.29 and a minimum of 0.05. A Monte-Carlo simulation of BTA1 with the same linkage disequilibrium structure and genetic linkage as the cattle family yielded a 99.98 and 99.94% of correct phases for informative SNPs in sire and calves, respectively.     

Extent of linkage disequilibrium and effective population size in Finnish Landrace and Finnish Yorkshire pig breeds

The extent of linkage disequilibrium (LD) and effective population size in Finnish Landrace and Finnish Yorkshire pig populations were studied using a whole genome SNP panel (Illumina PorcineSNP60 BeadChip) and pedigree data. Genotypic data included 86 Finnish Landrace and 32 Finnish Yorkshire boars. Pedigree data included 608,138 Finnish Landrace 554,237 and Finnish Yorkshire pigs, and on average 15 ancestral generations were known for the reference animals, born in 2005 to 2009. The breeding animals of the 2 populations have been kept separate in the breeding programs. Based on the pedigree data, the current effective population size for Finnish Landrace is 91 and for Finnish Yorkshire 61. Linkage disequilibrium measures (D' and r(2)) were estimated for over 1.5 million pairs of SNP. Average r(2) for SNP 30 kb apart was 0.47 and 0.49 and for SNP 5 Mb apart 0.09 and 0.12 for Finnish Landrace and Finnish Yorkshire, respectively. Average LD (r(2)) between adjacent SNP in the Illumina PorcineSNP60 BeadChip was 0.43 (57% of the adjacent SNP pairs had r(2) > 0.2) for Finnish Landrace and 0.46 (60% of the adjacent SNP pairs had r(2) > 0.2) for Finnish Yorkshire, and average r(2) > 0.2 extended to 1.0 and 1.5 Mb for Finnish Landrace and Finnish Yorkshire, respectively. Effective population size estimates based on the decay of r(2) with distance were similar to those based on the pedigree data: 80 and 55 for Finnish Landrace and Finnish Yorkshire, respectively. Thus, the results indicate that the effective population size of Finnish Yorkshire is smaller than of Finnish Landrace and has a clear effect on the extent of LD. The current effective population size of both breeds is above the recommended minimum of 50 but may get smaller than that in the near future, if no action is taken to balance the inbreeding rate and selection response. Because a moderate level of LD extends over a long distance, selection based on whole genome SNP markers (genomic selection) is expected to be efficient for both breeds.     

基于单倍型肉牛屠宰性状全基因组关联分析研究

单倍型标记与数量性状基因座(quantitative trait loci，QTL)之间具有较强的连锁不平衡(linkage disequilibrium，LD)关系，在基因定位和因果突变鉴定方面具有较高的应用价值。为了评估单倍型标记在基因组研究中的作用，本研究在华西牛资源群体中，选取该群体于2008—2021年间屠宰的共计1 478头平均月龄为24个月的个体进行研究，其中公牛1 333头，母牛145头。利用770K高密度芯片数据，基于LD阈值(r²>0.3)及固定单核苷酸多态(single nucleotide polymorphism,SNP)个数(5个连续SNP)两种方法进行单倍型构建，分别采用单位点SNP标记和两种单倍型标记共3种标记，基于GCTA的混合线性模型(mixed linear model,MLM)，开展宰前活重(LW)和屠宰率(DP)等屠宰性状的全基因组关联分析(genome-wide association study，GWAS)，定位影响屠宰性状的显著(P<0.05) SNPs、单倍型块和候选基因，同时比较3种标记的GWAS结果，评估3种标记的优劣。结果显示，3种标记在全基因组范围内共找到16个的显著SNPs及单倍型区域，主要分布于1、5、6、14、16、17和28号染色体上，同时鉴定到FAM184B、PPM1K、LCORL、RIMS2等10个与屠宰性状相关的候选基因，其中，基于SNP标记方法鉴定到的3个候选基因，在利用基于单倍型标记的方法中也鉴定到，且单倍型鉴定到的显著性位点或区域大多位于基因内部。在两种单倍型构建方法中，与基于固定SNP个数构建单倍型进行GWAS相比，基于LD阈值的构建方法鉴定到了更多候选基因。本研究结果表明，以单倍型开展GWAS可以综合考虑SNP位点间连锁关系，能较好地揭示复杂性状的遗传结构。     

Estimation of linkage disequilibrium in a sample of the United Kingdom dairy cattle population using unphased genotypes

The association between genetic marker alleles was estimated for two regions of the bovine genome from a random sample of 50 young dairy bulls born in the United Kingdom between 1988 and 1995. Microsatellite marker genotypes were obtained for six markers on chromosome 2 and seven markers on chromosome 6, spanning 38 and 20 cM, respectively. Two different methods, which do not require family information, were used to estimate population haplotype frequencies. Haplotype frequencies were estimated for pairs of loci using the expectation-maximization algorithm and for all linked loci using a Bayesian approach via a Markov chain-Monte Carlo algorithm. Significant (P = 0.0007) linkage disequilibrium was detected between pairs of loci in syntenic groups (that is, loci in the same linkage group), extending to about 10 cM. No significant linkage disequilibrium was detected between markers in nonsyntenic regions. Given the observed level of linkage disequilibrium, mapping methods based on population-wide association might provide a better resolution than traditional quantitative trait loci mapping methods in the U.K. dairy cattle population and may reduce the required sample sizes of the experiments.     

Genome‐wide linkage disequilibrium in the Blonde d'Aquitaine cattle breed

We present here the first genome‐wide characterization of linkage disequilibrium (LD) in the French Blonde d'Aquitaine (BLA) breed, a well‐muscled breed renowned for producing high‐yielding beef carcasses. To assess the pattern and extent of LD, we used a sample of 30 unrelated bulls and 36 923 single nucleotide polymorphisms (SNPs) covering all cattle autosomes. The squared correlation of the alleles at two loci (r²) was used as a measure of LD. The analysis of adjacent marker pairs revealed that the level of LD decreases rapidly with physical distance between SNPs. Overall mean r² was 0.205 (±0.262). Strong LD (r² > 0.8) and useful LD (measured as r² > 0.2) were observed within genomic regions of up to 720 and 724 kb, respectively. We analysed the genetic structure of the BLA population and found stratification. The observed genetic sub‐structuring is consistent with the known recent demographic history that occurred during BLA breed formation. Our results indicate that LD mapping of phenotypic traits in the BLA population is feasible; however, because of this sub‐structuring, special care is needed to reduce the likelihood of false‐positive associations between marker loci and traits of interest.     

Effective population size of an indigenous Swiss cattle breed estimated from linkage disequilibrium

Effective population size is an important parameter for the assessment of genetic diversity within a livestock population and its development over time. If pedigree information is not available, linkage disequilibrium (LD) analysis might offer an alternative perspective for the estimation of effective population size. In this study, 128 individuals of the Swiss Eringer breed were genotyped using the Illumina BovineSNP50 beadchip. We set bin size at 50 kb for LD analysis, assuming that LD for proximal single nucleotide polymorphism (SNP)‐pairs reflects distant breeding history while LD from distal SNP‐pairs would reflect near history. Recombination rates varied among different regions of the genome. The use of physical distances as an approximation of genetic distances (e.g. setting 1 Mb = 0.01 Morgan) led to an upward bias in LD‐based estimates of effective population size for generations beyond 50, while estimates for recent history were unaffected. Correction for restricted sample size did not substantially affect these results. LD‐based actual effective population size was estimated in the range of 87–149, whereas pedigree‐based effective population size resulted in 321 individuals. For conservation purposes, requiring knowledge of recent history (<50 generations), approximation assuming constant recombination rate seemed adequate.     

Refined mapping of quantitative trait loci for somatic cell score on BTA02 in the German Holstein

The aim of this study was to more precisely map a previously reported quantitative trait locus (QTL) affecting somatic cell score on Bos taurus autosome 2 by increasing the number of markers fourfold, analysing more families and exploiting within‐population linkage disequilibrium (LD). A granddaughter design of 10 German Holstein grandsire families with 1121 progeny tested sons was used. Twenty‐six markers with an average marker spacing of 3.14 cM were genotyped along 81.6 cM. Linkage analysis (LA) was performed using variance‐component methodology. The incorporation of LD was first done using variance‐component methods followed by regression on marker alleles. LA revealed genome‐wide significance (LOD > 3) at 15 contiguous marker‐intervals, with the maximum test‐statistic between DIK2862 and BMS778 and a 1‐lod drop‐off interval of 38 cM. While the variance‐component methods could not detect any LD, two individual markers with a significant effect (ILSTS098, p < 0.05; BMS778, p < 0.01) were found by regression analysis. Compared with previous results QTL‐localisation was substantially narrowed; further fine‐mapping should focus on the close vicinity of BMS778.     

Linkage of dermatomyositis in the Shetland Sheepdog to chromosome 35

Dermatomyositis is an inflammatory disease of the skin and muscle and is most commonly found in the Shetland sheepdog. Both the clinical presentation and the age of onset of dermatomyositis vary widely, and the inability to diagnose dermatomyositis before clinical symptoms ensue has made control of the disease difficult. Identification of a genetic marker that cosegregates with dermatomyositis would facilitate the development of a DNA-based test for the early detection of affected dogs. We report the use of linkage disequilibrium (LD) mapping to identify linkage to phenotypic dermatomyositis in the Shetland sheepdog. One marker, microsatellite marker FH3570 on canine chromosome 35, had evidence of LD (P=0.00002). Further studies are necessary to narrow the region harbouring the dermatomyositis locus, identify candidate genes and determine mode of inheritance.     

Genes controlling vaccine responses and disease resistance to respiratory viral pathogens in cattle

Farm animals remain at risk of endemic, exotic and newly emerging viruses. Vaccination is often promoted as the best possible solution, and yet for many pathogens, either there are no appropriate vaccines or those that are available are far from ideal. A complementary approach to disease control may be to identify genes and chromosomal regions that underlie genetic variation in disease resistance and response to vaccination. However, identification of the causal polymorphisms is not straightforward as it generally requires large numbers of animals with linked phenotypes and genotypes. Investigation of genes underlying complex traits such as resistance or response to viral pathogens requires several genetic approaches including candidate genes deduced from knowledge about the cellular pathways leading to protection or pathology, or unbiased whole genome scans using markers spread across the genome. Evidence for host genetic variation exists for a number of viral diseases in cattle including bovine respiratory disease and anecdotally, foot and mouth disease virus (FMDV). We immunised and vaccinated a cattle cross herd with a 40-mer peptide derived from FMDV and a vaccine against bovine respiratory syncytial virus (BRSV). Genetic variation has been quantified. A candidate gene approach has grouped high and low antibody and T cell responders by common motifs in the peptide binding pockets of the bovine major histocompatibility complex (BoLA) DRB3 gene. This suggests that vaccines with a minimal number of epitopes that are recognised by most cattle could be designed. Whole genome scans using microsatellite and single nucleotide polymorphism (SNP) markers has revealed many novel quantitative trait loci (QTL) and SNP markers controlling both humoral and cell-mediated immunity, some of which are in genes of known immunological relevance including the toll-like receptors (TLRs). The sequencing, assembly and annotation of livestock genomes and is continuing apace. In addition, provision of high-density SNP chips should make it possible to link phenotypes with genotypes in field populations without the need for structured populations or pedigree information. This will hopefully enable fine mapping of QTL and ultimate identification of the causal gene(s). The research could lead to selection of animals that are more resistant to disease and new ways to improve vaccine efficacy.     

Board-invited review: Applications of genomic information in livestock

The availability of whole genome sequences for individual species will change the landscape for livestock genomic research. Animal scientists will have access to whole-genome sequence-based technologies such as high-throughput SNP genotyping assays, gene expression profiling, methylation profiling, RNA interference, and genome resequencing that will revolutionize the scale upon which research will be conducted. These technologies will also alter the ways we think about addressing industry and scientific problems. In this review, we discuss the scientific bases for these emerging technologies and present recent highlights of their application in human, model species, and livestock as well as their potential for future applications in livestock. Additionally, we discuss strategies for their use in the genetic improvement and management of livestock. In particular, we present a strategy for the simultaneous identification of causal mutations underlying phenotypic traits in livestock and discuss issues that will arise in the application of whole genome selection for the prediction of genetic merit in livestock. We also point out that the statistical analysis that underlies the whole genome selection methodology is a sophisticated enhancement of single marker association mapping analysis to allow the entire genome to be simultaneously analyzed.     

绵羊基因组研究进展

过去几年中,家畜基因组计划取得了巨大进展。已经构建了猪、鸡、牛、绵羊和马的遗传图谱,其遗传标记间距在5-20cM。这些图谱对于家畜中与重要经济性状相关的基因或遗传标记的鉴定非常重要。文章从绵羊的遗传图谱、物理图谱、重要经济性状基因及QTL定位等方面对绵羊基因组的研究进展作了简要阐述。