Linkage disequilibrium and signatures of selection on chromosomes 19 and 29 in beef and dairy cattle

The objective of this study was to quantify the extent of linkage disequilibrium (LD) on bovine chromosomes 19 and 29 and to study the pattern of selection signatures in beef and dairy breeds (Angus and Holstein) of Bos taurus. The extent of LD was estimated for 370 and 186 single nucleotide polymorphism markers on BTA19 and 29 respectively using the square of the correlation coefficient (r(2)) among alleles at pairs of loci. A comparison of the extent of LD found that the decline of LD followed a similar pattern in both breeds. We observed long-range LD and found that LD dissipates to background levels at a locus separation of about 20 Mb on both chromosomes. Along each chromosome, patterns of LD were variable in both breeds. We find that a minimum of 30 000 informative and evenly spaced markers would be required for whole-genome association studies in cattle. In addition, we have identified chromosomal regions that show some evidence of selection for economically important traits in Angus and Holstein cattle. The results of this study are of importance for the design and application of association studies.     

Linkage disequilibrium structures in cattle and their application to breed identification testing

We examined the extent of linkage disequilibrium (LD) block lengths in four breed populations: Japanese Black, Angus, Hereford and Holstein. Three chromosomal regions in which QTL were previously mapped in Japanese Black populations were scanned with 84 microsatellite markers. The estimated LD lengths in these four purebred populations varied from 535 to 683 kb, which is much shorter than the values reported previously. Our findings suggest that QTL can be mapped in sub-centimorgan regions in these populations using an LD-mapping method. We also developed breed identification methods to distinguish Japanese Black from Angus, Hereford, Holstein and F₁ animals (Japanese Black × Holstein) respectively using the haplotypic frequencies of a pair of markers in the breed populations. After assessing the distributions of posterior probabilities to be Japanese Black, we obtained several pairs of markers that completely distinguished Japanese Black from the other breeds. We also obtained several combinations of six markers that completely distinguished Japanese Black animals from F₁ animals.     

Dissecting linkage disequilibrium in African-American genomes: roles of markers and individuals

Substantial increases of linkage disequilibrium (LD) both in magnitude and in range have been observed in recently admixed populations such as African-American (AfA). On the other hand, it has also been shown that LD in AfAs was very similar to that of African. In this study, we attempted to resolve these contradicting observations by conducting a systematic examination of the LD structure in AfAs by genotyping a sample of AfA individuals at 24,341 single nucleotide polymorphisms (SNPs) spanning almost the entire chromosome 21, with an average density of 1.5 kb/SNP. The overall LD in AfAs is similar to that in African populations and much less than that in European populations. Even when the ancestry-informative markers (AIMs) were used, extended LD in AfA was found to be limited to certain magnitude range (0.2 < or = r(2) < or = 0.8) and certain distance range, that is, between-marker distance more than 200 kb. Furthermore, the inclusion of AfA individuals with predominant African ancestry was found to reduce the overall magnitude of LD. Elevation of LD in the AfA population, compared with its parental populations, can only be observed at the markers with large allele frequency differences between 2 parental populations at limited scenario. AfA individuals of wholly African ancestry contribute little to the extended LD in the AfA population, and further genotyping or association analysis conducted using only admixed individuals may lead to higher statistical power and possibly reduced cost.     

Genome-wide scan for bovine twinning rate QTL using linkage disequilibrium

Twinning is a complex trait with negative impacts on health and reproduction, which cause economic loss in dairy production. Several twinning rate quantitative trait loci (QTL) have been detected in previous studies, but confidence intervals for QTL location are broad and many QTL are unreplicated. To identify genomic regions or genes associated with twinning rate, QTL analysis based on linkage combined with linkage disequilibrium (LLD) and individual marker associations was conducted across the genome using high-throughput single nucleotide polymorphism (SNP) genotypes. A total of 9919 SNP markers were genotyped with 200 sires and sons in 19 half-sib North American Holstein dairy cattle families. After SNPs were genotyped, informative markers were selected for genome-wide association tests and QTL searches. Evidence for twinning rate QTL was found throughout the genome. Thirteen markers significantly associated with twinning rate were detected on chromosomes 2, 5 and 14 ( P  < 2.3 × 10⁻⁵). Twenty-six regions on fourteen chromosomes were identified by LLD analysis at P  < 0.0007. Seven previously reported ovulation or twinning rate QTL were supported by results of single marker association or LLD analyses. Single marker association analysis and LLD mapping were complementary tools for the identification of putative QTL in this genome scan.     

利用全基因组连锁不平衡估计中国荷斯坦牛有效群体大小

有效群体大小是群体遗传学研究的一个重要内容,有助于我们更清楚地了解群体的遗传变异、进化和复杂性状的遗传机制等。随着高密度SNP标记的出现,越来越多的研究利用SNP标记间连锁不平衡估计有效群体大小。文章采集北京地区中国荷斯坦牛2 093个样本,并利用牛SNP芯片(Illumina BovineSNP50,含5 4001 SNPs)进行基因型测定,估计不同世代中国荷斯坦牛的有效群体大小。质量控制标准设定为SNP检出率0.95,最小等位基因频率>0.05,样本检出率0.95,哈代温伯格平衡检验显著性水平P<0.0001。经过质量控制,共1 968个样本和38 796个SNPs用于连锁不平衡分析。文章选取SNP间距0.1、0.2、0.5、1、2、5、10、15(Mb),估计中国荷斯坦牛在4世代之前有效群体大小。结果表明,中国荷斯坦牛的有效群体呈逐代下降趋势,至4世代前,中国荷斯坦牛平均有效群体为45头左右。     

单核苷酸多态性基因分型技术原理与进展

在基因组规模了解遗传变异与生物功能之间的关系可望为生物学带来全新的深入认识。本从等位基因分型机理、反应形式和检测方法等三个方面讨论SNP分型方法的现状,并简要介绍了目前应用的一些分型方法。     

Genome-wide linkage disequilibrium in two Japanese beef cattle breeds

There is little knowledge about the degree of linkage disequilibrium (LD) in beef cattle. This study aims to perform a genome-wide search for LD in Japanese Black and Japanese Brown beef cattle and to compare the level of LD between these two breeds. Parameter D' (the LD coefficient) was used as a measure of LD, and LD was tested for significance of allelic associations between syntenic and between non-syntenic marker pairs. Effects of breed, chromosome, genetic map distance and their interactions with D' were tested based on least squares analyses. Both breeds showed high levels of LD, which ranged over several tens of cM and declined as the marker distance increased for syntenic marker pairs. A rapid decline of the D' value was observed between markers that were spaced 5 and 20 cM apart. LD was significant in most cases for marker pairs <40 cM apart but was not significant between non-syntenic loci. The pattern of LD found in these two breeds was similar to that previously published for dairy cattle. The D' value between breeds was not significantly different (P > 0.05), but the interaction between breed and chromosome was highly significant (P < 0.001). Genetic selection seems to have caused the heterogeneity of the D' values among chromosomes within breed. These results indicate that LD mapping is a useful tool for fine-mapping quantitative trait loci of economically important traits in Japanese beef cattle.     

Genomic hybridization of bovine class II major histocompatibility genes: 2. Polymorphism of DR genes and linkage disequilibrium in the DQ-DR region

Summary. Class II genes of the bovine major histocompatibility complex have been investigated by Southern blot analysis using human cDNA probes for DQα, DQβ, DRα and DRβ. In this report restriction fragment length polymorphisms of DR α and DR β are described. The polymorphisms were interpreted genetically by analysing five paternal half-sib families of the Swedish Red and White Breed, comprising altogether 28 offspring. Using the restriction enzymes Bam HI, Eco RI and Pvu II, three DRα and three DRβ allelic fragment patterns were resolved. The DR α and DR β genes thus appear to be much less polymorphic than the previously described DQ α and DQ β genes. Also, the observed linkage disequilibrium between DR genes was less pronounced than that between DQ genes, whereas the association between DR and DQ haplotypes was very strong. The family data available indicated strongly that the DQ α, DQ β, DR α and DR β genes are all closely linked.     

Linkage disequilibrium for different scales and applications

Assessing the patterns of linkage disequilibrium (LD) has become an important issue in both evolutionary biology and medical genetics since the rapid accumulation of densely spaced DNA sequence variation data in several organisms. LD deals with the correlation of genetic variation at two or more loci or sites in the genome within a given population. There are a variety of LD measures which range from traditional pairwise LD measures such as D' or r2 to entropy-based multi-locus measures or haplotype-specific approaches. Understanding the evolutionary forces (in particular recombination) that generate the observed variation of LD patterns across genomic regions is addressed by model-based LD analysis. Marker type and its allelic composition also influence the observed LD pattern, microsatellites having a greater power to detect LD in population isolates than SNPs. This review aims to explain basic LD measures and their application properties.     

The extent of linkage disequilibrium in a large cattle population of western Africa and its consequences for association studies

Several previous studies concluded that linkage disequilibrium (LD) in livestock populations from developed countries originated from the impact of strong selection. Here, we assessed the extent of LD in a cattle population from western Africa that was bred in an extensive farming system. The analyses were performed on 363 individuals in a Bos indicus x Bos taurus population using 42 microsatellite markers on BTA04, BTA07 and BTA13. A high level of expected heterozygosity (0.71), a high mean number of alleles per locus (9.7) and a mild shift in Hardy-Weinberg equilibrium were found. Linkage disequilibrium extended over shorter distances than what has been observed in cattle from developed countries. Effective population size was assessed using two methods; both methods produced large values: 1388 when considering heterozygosity (assuming a mutation rate of 10(-3)) and 2344 when considering LD on whole linkage groups (assuming a constant population size over generations). However, analysing the decay of LD as a function of marker spacing indicated a decreasing trend in effective population size over generations. This decrease could be explained by increasing selective pressure and/or by an admixture process. Finally, LD extended over small distances, which suggested that whole-genome scans will require a large number of markers. However, association studies using such populations will be effective.     

Linkage disequilibrium fine mapping of quantitative trait loci: A simulation study

Recently, the use of linkage disequilibrium (LD) to locate genes which affect quantitative traits (QTL) has received an increasing interest, but the plausibility of fine mapping using linkage disequilibrium techniques for QTL has not been well studied. The main objectives of this work were to (1) measure the extent and pattern of LD between a putative QTL and nearby markers in finite populations and (2) investigate the usefulness of LD in fine mapping QTL in simulated populations using a dense map of multiallelic or biallelic marker loci. The test of association between a marker and QTL and the power of the test were calculated based on single-marker regression analysis. The results show the presence of substantial linkage disequilibrium with closely linked marker loci after 100 to 200 generations of random mating. Although the power to test the association with a frequent QTL of large effect was satisfactory, the power was low for the QTL with a small effect and/or low frequency. More powerful, multi-locus methods may be required to map low frequent QTL with small genetic effects, as well as combining both linkage and linkage disequilibrium information. The results also showed that multiallelic markers are more useful than biallelic markers to detect linkage disequilibrium and association at an equal distance.     

Analysis of population differentiation in North Eurasian cattle (Bos taurus) using single nucleotide polymorphisms in three genes associated with production traits

Single nucleotide polymorphisms (SNPs) in growth hormone 1 (GH1), insulin-like growth factor 1 (IGF1) and leptin (LEP), all candidates for production traits in cattle, were characterized in North Eurasian cattle breeds. Allele frequencies of IGF1 exhibited significant (P < 0.05) deviation from neutral expectation and therefore, might be associated with divergence in North Eurasian cattle because of genetic selection. Allele frequencies and lower heterozygosity of LEP may indicate a recent introduction of an alternative allele in this geographic region. Locus F(ST) estimates were highest for IGF1 (0.151, sigma = 0.042) and lowest for GH (0.062, sigma = 0.020). Our results suggest a slightly higher population differentiation across the candidate genes (FST = 0.108) than across microsatellites (FST = 0.095), possibly because of selection and stochastic effects.     

Nucleotide variability and linkage disequilibrium patterns at the porcine FABP5 gene

Fatty acid binding protein 5 (FABP5) is a major positional and physiological candidate gene for the porcine FAT1 QTL on SSC4. Here we characterize the nucleotide polymorphism and haplotype variability of FABP5 and we compare it with that of FABP4, given their close physical location and similar metabolic roles. DNA resequencing of the FABP5 gene region in 29 pigs from 14 breeds and in European and Japanese wild boars revealed 36 polymorphisms in 5.2 kb, and a nucleotide diversity of 0.19%, comparable to values reported in other domestic species but sixfold lower than that previously found for FABP4. Remarkably, both the nucleotide variability and the haplotype structure of FABP5 and FABP4 were dramatically different, and the Hudson-Kreitman-Aguadé test was highly significant. Nevertheless, both genes also had similarities. The neighbour-joining trees of their haplotypes did not show a geographical arrangement for any of the genes. Besides, both genes presented a similar extent and pattern of linkage disequilibrium. Haplotype blocks did not extend for large stretches ( approximately 1 kb in both genes), and the number of tag SNPs required to capture all variability was higher than previously expected. Our findings indicate that FABP4 and FABP5 have undergone different selective or evolutive processes. The fact that haplotype blocks were so small may require us to increase the number of SNPs in prospective whole-genome association studies in the pig.     

Linkage mapping of gene-associated SNPs to pig chromosome 11

Single nucleotide polymorphisms (SNPs) were discovered in porcine expressed sequence tags (ESTs) orthologous to genes from human chromosome 13 (HSA13) and predicted to be located on pig chromosome 11 (SSC11). The SNPs were identified as sequence variants in clusters of EST sequences from pig cDNA libraries constructed in the Sino-Danish pig genome project. In total, 312 human gene sequences from HSA13 were used for similarity searches in our pig EST database. Pig ESTs showing significant similarity with HSA13 genes were clustered and candidate SNPs were identified. Allele frequencies for 26 SNPs were estimated in a group of 80 unrelated pigs from Danish commercial pig breeds: Duroc, Hampshire, Landrace and Large White. Eighteen of the 26 SNPs genotyped in the PiGMaP Reference Families were mapped by linkage analysis to SSC11. The EST-based SNPs published here are new genetic markers useful for linkage and association studies in commercial and experimental pig populations. This study represents the first gene-associated SNP linkage map of pig chromosome 11 and adds new comparative mapping information between SSC11 and HSA13. Furthermore, our data facilitate future studies aimed at the identification of interesting regions on pig chromosome 11, positional cloning and fine mapping of quantitative trait loci in pig.     

Linkage disequilibrium in Bougainville Island

Linkage disequilibria are estimated for three 2-locus systems in 18 samples from Bougainville Island, Solomon Islands. The systems are haptoglobin, acid phosphatase and MN blood group. The disequilibria are estimated two ways: by maximum likelihood (ML) and by the covariance between the non-alleles. Though seven of the 52 ML estimates are statistically different than zero, none of the covariance estimates are significant. We conclude that because linkage disequilibrium for loosely linked loci is a small quantity and because the sample sizes for most populations studied by anthropologists are small, linkage disequilibrium will not be a useful parameter for the study of natural selection in these populations.     

Genome-Wide Association Studies Using Single Nucleotide Polymorphism Markers Developed by Re-Sequencing of the Genomes of Cultivated Tomato

With the aim of understanding relationship between genetic and phenotypic variations in cultivated tomato, single nucleotide polymorphism (SNP) markers covering the whole genome of cultivated tomato were developed and genome-wide association studies (GWAS) were performed. The whole genomes of six tomato lines were sequenced with the ABI-5500xl SOLiD sequencer. Sequence reads covering ∼13.7× of the genome for each line were obtained, and mapped onto tomato reference genomes (SL2.40) to detect ∼1.5 million SNP candidates. Of the identified SNPs, 1.5% were considered to confer gene functions. In the subsequent Illumina GoldenGate assay for 1536 SNPs, 1293 SNPs were successfully genotyped, and 1248 showed polymorphisms among 663 tomato accessions. The whole-genome linkage disequilibrium (LD) analysis detected highly biased LD decays between euchromatic (58 kb) and heterochromatic regions (13.8 Mb). Subsequent GWAS identified SNPs that were significantly associated with agronomical traits, with SNP loci located near genes that were previously reported as candidates for these traits. This study demonstrates that attractive loci can be identified by performing GWAS with a large number of SNPs obtained from re-sequencing analysis.     

单核苷酸多态性在作物遗传及改良中的应用

单核苷酸多态性（single nucleotide polymorphism,SNP）是等位基因间序列差异最为普遍的类型,可作为一种高通量的遗传标记。已建立了PCR扩增目标序列及其产物测序和电子SNP（eSNP）等多种发现和检测SNP的方法。玉米和大豆等作物也已开展了SNP分析。一些栽培作物种质的多样性不断减少,其结果使连锁不平衡（linkage disequilibrium,LD）增加,这有利于目的基因座上SNP单元型（haplotype）与表型的相关性分析。SNP已在作物基因作图及其整合、分子标记辅助育种和功能基因组学等领域展示了广泛的应用价值。Abstract:Single nucleotide polymorphism(SNP) is the most common type of sequence difference between alleles,which can be used as a kind of high-throughput genetic marker.Several different routes have been developed to discover and identify SNP.These include the direct sequencing of PCR amplicons,electronic SNP(eSNP) and so on.SNP assays have been made in many crop species such as maize and soybean.The elite germplasm of some crops have been narrowed in genetic diversity,increasing the amount of linkage disequilibrium(LD) present and facilitating the association of SNP haplotypes at candidate gene loci with phenotypes.SNP analysis has been broadly used in the field of plant gene mapping,integration of genetic and physical maps,DNA marker-assisted breeding and functional genomics.