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.     

Genome-wide linkage disequilibrium analysis in bread wheat and durum wheat.

Bread wheat and durum wheat were examined for linkage disequilibrium (LD) using microsatellite markers distributed across the genome. The allele database consisted of 189 bread wheat accessions genotyped at 370 loci and 93 durum wheat accessions genotyped at 245 loci. A significance level of p < 0.001 was set for all comparisons. The bread and durum wheat collections showed that 47.9% and 14.0% of all locus pairs were in LD, respectively. LD was more prevalent between loci on the same chromosome compared with loci on independent chromosomes and was highest between adjacent loci. Only a small fraction (bread wheat, 0.9%; durum wheat, 3.2%) of the locus pairs in LD showed R2 values > 0.2. The LD between adjacent locus pairs extended (R2 > 0.2) approximately 2-3 cM, on average, but some regions of the bread and durum wheat genomes showed high levels of LD (R2 = 0.7 and 1.0, respectively) extending 41.2 and 25.5 cM, respectively. The wheat collections were clustered by similarity into subpopulations using unlinked microsatellite data and the software Structure. Analysis within subpopulations showed 14- to 16-fold fewer locus pairs in LD, higher R2 values for those pairs in LD, and LD extending further along the chromosome. The data suggest that LD mapping of wheat can be performed with simple sequence repeats to a resolution of <5 cM.     

The distribution of long range admixture linkage disequilibrium in an African-American population

OBJECTIVES: To better understand the effect of admixture on long range linkage disequilibrium (LD), we characterized extended LD in gene-rich regions of an African-American population. METHODS: Approximately 290 cM of chromosomes 1, 3, 6, 11-17, 20 and 22 were scanned using 109 polymorphic microsatellite markers spaced an average of 3 cM apart. Disequilibrium between loci (D') was based on maximum-likelihood estimates of haplotype frequencies computed for 200 unrelated African Americans. RESULTS: Mean D' values were highest on chromosomes 6p23-p21.3 (D' = 0.33) and 15p22.2-p25.3 (D' = 0.34), and lowest on chromosome 12p11.2-q14 (D' = 0.21). Overall, the variance in LD among chromosomes accounted for approximately two-thirds of the total LD variance. Of the 434 locus pairs spaced between 0.3 and 38.7 cM apart, there was no detectable correlation between LD and recombination distance and a weak negative correlation between LD and physical distance (r(s) = -0.12; p = 0.031). For the 192 intrachromosomal locus pairs where allele frequency data were available from the Centre d'Etude du Polymorphisme humain (CEPH), we found a statistically significant positive correlation between LD and the allelic frequency differences (delta) between the African-American study population and Caucasian reference CEPH population (r(s) = 0.53; p < 0.0001). The correlation between LD and both recombination and physical distance was markedly increased for locus pairs with high delta levels. CONCLUSIONS: Our results suggest that recent Caucasian admixture maintains a high level of long range LD in African Americans on a genomic scale, and selected markers with large African American/Caucasian delta levels may be useful in association studies.     

Genetic diversity and linkage disequilibrium in Chinese bread wheat (Triticum aestivum L.) revealed by SSR markers

Two hundred and fifty bread wheat lines, mainly Chinese mini core accessions, were assayed for polymorphism and linkage disequilibrium (LD) based on 512 whole-genome microsatellite loci representing a mean marker density of 5.1 cM. A total of 6,724 alleles ranging from 1 to 49 per locus were identified in all collections. The mean PIC value was 0.650, ranging from 0 to 0.965. Population structure and principal coordinate analysis revealed that landraces and modern varieties were two relatively independent genetic sub-groups. Landraces had a higher allelic diversity than modern varieties with respect to both genomes and chromosomes in terms of total number of alleles and allelic richness. 3,833 (57.0%) and 2,788 (41.5%) rare alleles with frequencies of <5% were found in the landrace and modern variety gene pools, respectively, indicating greater numbers of rare variants, or likely new alleles, in landraces. Analysis of molecular variance (AMOVA) showed that A genome had the largest genetic differentiation and D genome the lowest. In contrast to genetic diversity, modern varieties displayed a wider average LD decay across the whole genome for locus pairs with r(2)>0.05 (P<0.001) than the landraces. Mean LD decay distance for the landraces at the whole genome level was <5 cM, while a higher LD decay distance of 5-10 cM in modern varieties. LD decay distances were also somewhat different for each of the 21 chromosomes, being higher for most of the chromosomes in modern varieties (<5 ～ 25 cM) compared to landraces (<5 ～ 15 cM), presumably indicating the influences of domestication and breeding. This study facilitates predicting the marker density required to effectively associate genotypes with traits in Chinese wheat genetic resources.     

The extent and distribution of linkage disequilibrium in a multi-hierarchic outbred canine pedigree

A canine integrated linkage-radiation map has been recently constructed by using microsatellite markers. This map, with a good coverage of the canine genome, allows for a genome-wide search for the extent and distribution of linkage disequilibrium derived from linkage and evolutionary forces. In this study, we genotyped an outbred pedigree between Labrador retriever and Greyhound breeds with a set of microsatellite markers (240) from the canine linkage map. Linkage disequilibrium was measured between all syntenic and nonsyntenic marker pairs. Analysis of syntenic pairs revealed a significant correlation (–0.229, P < 0.001) between linkage disequilibrium and genetic distance (log transformed). Significant linkage disequilibria were observed more frequently between syntenic pairs spaced <40 cM than those paced >40 cM. There is a clear trend for linkage disequilibrium to decline with marker distance. From our results, a genome-wide screen with markers at low to moderate density (1–2 per 10 cM) should take full advantage of linkage disequilibrium for quantitative trait locus mapping in dogs. This study supports the appropriateness of linkage disequilibrium analysis to detect and map quantitative trait loci underlying complex traits in dogs.     

Levels of linkage disequilibrium in a wild bird population

Population-based mapping approaches are attractive for tracing the genetic background to phenotypic traits in wild species, given that it is often difficult to gather extensive and well-defined pedigrees needed for quantitative trait locus analysis. However, the feasibility of association or hitch-hiking mapping is dependent on the degree of linkage disequilibrium (LD) in the population, on which there is yet limited information for wild species. Here we use single nucleotide polymorphism (SNP) markers from 23 genes in a recently established linkage map of the Z chromosome of the collared flycatcher, to study the extent of LD in a natural bird population. In most but not all cases we find SNPs within the same intron (less than 500 bp) to be in perfect LD. However, LD then decays to background level at a distance 1cM or 400-500 kb. Although LD seems more extensive than in other species, if the observed pattern is representative for other regions of the genome and turns out to be a general feature of natural bird populations, dense marker maps might be needed for genome scans aimed at identifying association between marker and trait loci.     

Fine mapping of the chromosome 11q22-23 region using PFGE, linkage and haplotype analysis; localization of the gene for ataxia telangiectasia to a 5cM region flanked by NCAM/DRD2 and STMY/CJ52.75, phi 2.22.

Using pulsed-field gel electrophoresis, and a range of different enzyme digests, we have established that both markers of each of the pairs CJ52.208/YNB3.12, NCAM/DRD2, and STMY/CJ52.75, on chromosome 11q22-23, show physical linkage on a single DNA fragment. We have also shown, using genetic linkage and haplotype analyses, that these markers lie within a region of approximately 18cM, which, it has been shown previously, is likely to contain the A-T gene. The relative positions of these marker loci, and the distance between them was determined in order to construct a detailed map which has allowed a more precise localization of the A-T gene. We have shown that in pairwise linkage analysis the strongest support for linkage to the A-T gene was with the STMY/CJ52.75 locus (Z = 5.59, theta = 0.0). A three-point analysis using the results from STMY/CJ52.75 and the closely linked marker phi 2.22 gave Z = 5.55, theta = 0.03. Despite persisting evidence of some linkage to Thy-1 our results are consistent with the existence of a single A-T locus on chromosome 11q22-23 and our best estimate of the position of this locus places it between NCAM/DRD2 and (STMY/CJ52.75, F2.22) (Z = 6.74), a region of approximately 5cM in males.     

Linkage mapping of a polymorphic plumage locus associated with intermorph incompatibility in the Gouldian finch (Erythrura gouldiae)

Colour polymorphism is known to facilitate speciation but the genetic basis of animal pigmentation and how colour polymorphisms contribute to speciation is poorly understood. Restricted recombination may promote linkage disequilibrium between the colour locus and incompatibility genes. Genomic rearrangement and the position of relevant loci within a chromosome are important factors that influence the frequency of recombination. Therefore, it is important to know the position of the colour locus, gene order and recombination landscape of the chromosome to understand the mechanism that generates incompatibilities between morphs. Recent studies showed remarkable pre- and postzygotic incompatibilities between sympatric colour morphs of the Gouldian finch (Erythrura gouldiae), in which head feather colour is genetically determined by a single sex-linked locus, Red. We constructed a genetic map for the Z chromosome of the Gouldian finch (male-specific map distance=131 cM), using 618 captive-bred birds and 34 microsatellite markers, to investigate the extent of inter- and intraspecific genomic rearrangements and variation in recombination rate within the Z chromosome. We refined the location of the Red locus to a ~7.2-cM interval in a region with a moderate recombination rate but outside the least-recombining, putative centromeric region. There was no evidence of chromosome-wide genomic rearrangements between the chromosomes carrying the red or black alleles with the current marker resolution. This work will contribute to identifying the causal gene, which will in turn enable alternative explanations for the association between incompatibility and colouration, such as fine-scale linkage disequilibrium, genomic rearrangements and pleiotropy, to be tested.     

Linkage disequilibrium in related breeding lines of chickens

High-density genotyping of single-nucleotide polymorphisms (SNPs) enables detection of quantitative trait loci (QTL) by linkage disequilibrium (LD) mapping using LD between markers and QTL and the subsequent use of this information for marker-assisted selection (MAS). The success of LD mapping and MAS depends on the extent of LD in the populations of interest and the use of associations across populations requires LD between loci to be consistent across populations. To assess the extent and consistency of LD in commercial broiler breeding populations, we used genotype data for 959 and 398 SNPs on chromosomes 1 and 4 on 179-244 individuals from each of nine commercial broiler chicken breeding lines. Results show that LD measured by r(2) extends over shorter distances than reported previously in other livestock breeding populations. The LD at short distance (within 1 cM) tended to be consistent across related populations; correlations of LD measured by r for pairs of lines ranged from 0.17 to 0.94 and closely matched the line relationships based on marker allele frequencies. In conclusion, LD-based correlations are good estimates of line relationships and the relationship between a pair of lines a good predictor of LD consistency between the lines.     

Detecting population structure using STRUCTURE software: effect of background linkage disequilibrium

STRUCTURE is the most widely used clustering software to detect population genetic structure. The last version of this software (STRUCTURE 2.1) has been enhanced recently to take into account the occurrence of linkage disequilibrium (LD) caused by admixture between populations. This last version, however, still does not consider the effects of strong background LD caused by genetic drift, and which may cause spurious results. STRUCTURE authors have, therefore, suggested a rough threshold value of the distance (1.0 cM) between two loci below which the pair of loci should not be used. Because of the sensitiveness of LD to demographic events, the distance between loci is not always a good indicator of the strength of LD. In this study, we examine the link between genomic distance and the strength of the correlation between loci (r(LD)) in a free-ranging population of mouflon (Ovis aries), and we present an empirical test of effect of r(LD) on the clustering results provided by the linkage model in STRUCTURE. We showed that a high r(LD) value increases the probability of detecting spurious clustering. We propose to use r(LD) as an index to base a decision on whether or not to use a pair of loci in a clustering analysis.     

An AFLP linkage map for Douglas-fir based upon multiple full-sib families

An amplified fragment length polymorphism (AFLP) linkage map for coastal Douglas-fir (Pseudotsuga menziesii) was constructed from eight full-sib families each consisting of 40 progeny. These families were part of the British Columbia Ministry of Forests second-generation progeny test program and represent typical family sizes used in progeny trials. For map construction, ten primer pairs using EcoRI+3 and MseI+4 were employed to identify and assay AFLP loci that segregated in backcross configurations. A new technique was used to obtain a single recombination rate for each pair of marker loci: for each locus pair, a recombination rate and log-odd value were estimated across all segregating families using a joint maximum likelihood function that considered the full dataset of segregating genotypes. The resulting matrix of recombination rates between all pairs of loci was used to construct an integrated linkage map using JoinMap. The final map consisted of 19 linkage groups spanning 938.6 cM at an average distance of 9.3 cM between markers. The simultaneous integration of data from multiple families may provide an effective way to construct a linkage map, using the genetic resources inherent in most tree improvement programs, where progeny tests of small size are conducted. The statistical property of number of families used is briefly discussed. For our data, at least three to four families greatly increased the chance of obtaining an informative locus in at least one family. Families as small as ten are adequate for closely linked loci (<10 cM), while the size used in our study (40) is adequate for loci within 30 cM.     

Absence of population genetic differentiation in the New Zealand greenshell mussel Perna canaliculus (Gmelin 1791) as assessed by allozyme variation

Genetic variation in the endemic New Zealand greenshell mussel, Perna canaliculus (Gmelin 1791), was examined using starch-gel electrophoresis at seven protein-coding loci (Idh; Acon-1; Acon-2; Gpd; Pgi; Pgm; Pgd) in 35 populations (N=1038 mussels). For all loci and all populations, Fisher's exact tests indicated highly significant departures from Hardy-Weinberg equilibrium (HWE), but this overall result was caused by significant heterozygote deficiencies at only two loci (Pgm and Pgd), and in only three northern populations (Kuaotunu, Te Haumi and Days Bay). Allelic and genotypic differentiation between population pairs at individual loci and across all loci were nonsignificant, and genotypic disequilibrium at each locus pair was also nonsignificant for all populations. Genetic variation in all populations was high (mean heterozygosity, 0.210+/-0.113), while Nei's D among populations was very low (0.002+/-0.002). Low population subdivision (θ=-0.001+/-0.002) and high levels of gene flow (Nm(p)=10.18; Nm(θ)=infinity) also indicated that the single panmictic unit model best explains population genetic homogeneity in P. canaliculus over a north-south distance >2000 km. Lack of genetic subdivision in this species is discussed in light of two previous allozyme studies, with differing results: one suggested that a north-south division exists between greenshell mussel stocks, and the other suggested that population structure in this species can be explained through isolation by distance model modified by local hydrology.     

Genetic mapping of DNA segments relative to the locus for the fragile-X syndrome at Xq27.3.

We have tested linkage between the locus for the fragile-X [fra(X)] syndrome at Xq27.3 and five polymorphic restriction sites identified by four DNA probes mapping distal to Xq26.1. A maximum distance of approximately 15 centimorgans (cM) between Xq27.3 and the marker loci mapping to this region was predicted based on the physical chromosome length. Close linkage between the disease and marker loci was excluded for probes DXS19 and DXS37 (theta = .05, Z = -2.94 and Z = -4.17, respectively). These marker loci were estimated to be less than five cM apart but approximately 40 cM proximal to the fragile site, indicating that there is a significantly greater frequency of recombination in this region of the X chromosome than expected from the physical length. Linkage results for the other marker loci and the fra(X) syndrome were inconclusive. However, the pX45d probe locus appears very closely linked to the factor IX locus (Z = 1.94 at theta = 0) and is approximately 20 cM proximal to Xq27.3. A relative map of the polymorphic restriction sites, fra(X) syndrome locus, and factor IX locus was constructed by maximizing lod scores over the Xq26.1----q27.3 region.     

Demarcation of informative chromosomes in tropical sweet corn inbred lines using microsatellite DNA markers

A study of genetic variation among 10 pairs of chromosomes extracted from 13 tropical sweet corn inbred lines, using 99 microsatellite markers, revealed a wide range of genetic diversity. Allelic richness and the number of effective alleles per chromosome ranged from 2.78 to 4.33 and 1.96 to 3.47, respectively, with respective mean values of 3.62 and 2.73. According to the Shannon’s information index (I) and Nei’s gene diversity coefficient (Nei), Chromosome 10 was the most informative chromosome (I = 1.311 and Nei = 0.703), while Chromosome 2 possessed the least (I = 0.762 and Nei = 0.456). Based on linkage disequilibrium (LD) measurements for loci less than 50 cM apart on the same chromosome, all loci on Chromosomes 1, 6 and 7 were in equilibrium. Even so, there was a high proportion of genetic variation in Chromosomes 4, 5, 8, 9 and 10, thereby revealing their appropriateness for use in the genetic diversity investigations among tropical sweet corn lines. Chromosome 4, with the highest number of loci in linkage disequilibrium, was considered the best for marker-phenotype association and QTL mapping, followed by Chromosomes 5, 8, 9 and 10.     

Ten novel microsatellite loci characterized for a remarkably widespread fish: Galaxias maculatus (Galaxiidae)

Ten polymorphic microsatellite markers (five tetra‐, one compound tetra‐, one octa‐ and three dinucleotides) were isolated and characterized for Galaxias maculatus, a fish species widely distributed in the Southern Hemisphere. Markers were tested in 89 individual samples from a single location and the number of alleles ranged between 2 and 28. Observed and expected heterozygosities ranged from 0.103 to 0.910 and 0.098 to 0.935 respectively. No evidence was detected for either linkage disequilibrium (P‐values > 0.05 for each locus pair) or deviations from HWE (P‐values > 0.05 for every loci).     

Population structure and long-range linkage disequilibrium in a durum wheat elite collection

A collection of 134 durum wheat accessions, mainly including cultivars (cvs.) representative of the major gene pools, was assembled and characterized with 70 SSRs for genetic diversity and level of long-range linkage disequilibrium (LD). Results of both a distance-based and a model-based (Bayesian) cluster analysis evidenced the presence of a structured diversity. In the model-based analysis, six to eight main distinct subpopulations were identified based on the molecular data. Only a relatively small portion (20%) of the molecular variation was accounted for by the geographical origin of the accessions. Major differences were detected between the North American and the Mediterranean cvs., while a considerable overlap characterized the cvs. from CIMMYT-ICARDA and Italy. The North American cvs. showed the highest within group mean genetic similarity (GS_m = 0.68). French cvs. revealed sizeable similarities with both the North American as well as the Italian and CIMMYT-ICARDA pools. Considering the germplasm as a whole, high levels of LD were found both at locus pairs with an intrachromosomal distance <50 cM as well as at those with distances more than 50 cM and independent (86, 52 and 54% of SSR pairs at p < 0.01, respectively). After re-evaluating LD within each of the three main subgroups identified through the analysis of the germplasm structure, the LD level remained high for tightly to moderately linked locus pairs (<20 cM apart), but was greatly reduced in the loosely linked (more than 50 cM apart) and independent locus pairs. The implications of these findings as to the possibility of using association mapping for gene/QTL discovery in durum wheat are discussed.     

Evidence of a non-MHC susceptibility locus in type I diabetes linked to HLA on chromosome 6.

Linkage studies have led to the identification of several chromosome regions that may contain susceptibility loci to type I diabetes (IDDM), in addition to the HLA and INS loci. These include two on chromosome 6q, denoted IDDM5 and IDDM8, that are not linked to HLA. In a previous study, we noticed that the evidence for linkage to IDDM susceptibility around the HLA locus extended over a total distance of 100 cM, which suggested to us that another susceptibility locus could reside near HLA. We developed a statistical method to test this hypothesis in a panel of 523 multiplex families from France, the United States, and Denmark (a total of 667 affected sib pairs, 536 with both parents genotyped), and here present evidence (P = .00003) of a susceptibility locus for IDDM located 32 cM from HLA in males but not linked to HLA in females and distinct from IDDM5 and IDDM8. A new statistical method to test for the presence of a second susceptibility locus linked to a known first susceptibility locus (here HLA) is presented. In addition, we analyzed our current family panel with markers for IDDM5 and IDDM8 on chromosome 6 and found suggestions of linkage for both of these loci (P = .002 and .004, respectively, on the complete family panel). When cumulated with previously published results, with overlapping families removed, the affected-sib-pair tests had a significance of P = .0001 for IDDM5 and P = .00004 for IDDM8.     

Osteoarthritis-susceptibility locus on chromosome 11q, detected by linkage.

We present a two-stage genomewide scan for osteoarthritis-susceptibility loci, using 481 families that each contain at least one affected sibling pair. The first stage, with 272 microsatellite markers and 297 families, involved a sparse map covering 23 chromosomes at intervals of approximately 15 cM. Sixteen markers that showed evidence of linkage at nominal P相似文献   

Linkage disequilibrium and persistence of phase in Holstein-Friesian, Jersey and Angus cattle

When a genetic marker and a quantitative trait locus (QTL) are in linkage disequilibrium (LD) in one population, they may not be in LD in another population or their LD phase may be reversed. The objectives of this study were to compare the extent of LD and the persistence of LD phase across multiple cattle populations. LD measures r and r(2) were calculated for syntenic marker pairs using genomewide single-nucleotide polymorphisms (SNP) that were genotyped in Dutch and Australian Holstein-Friesian (HF) bulls, Australian Angus cattle, and New Zealand Friesian and Jersey cows. Average r(2) was approximately 0.35, 0.25, 0.22, 0.14, and 0.06 at marker distances 10, 20, 40, 100, and 1000 kb, respectively, which indicates that genomic selection within cattle breeds with r(2) >or= 0.20 between adjacent markers would require approximately 50,000 SNPs. The correlation of r values between populations for the same marker pairs was close to 1 for pairs of very close markers (<10 kb) and decreased with increasing marker distance and the extent of divergence between the populations. To find markers that are in LD with QTL across diverged breeds, such as HF, Jersey, and Angus, would require approximately 300,000 markers.