A first-generation metric linkage disequilibrium map of bovine chromosome 6

We constructed a metric linkage disequilibrium (LD) map of bovine chromosome 6 (BTA6) on the basis of data from 220 SNPs genotyped on 433 Australian dairy bulls. This metric LD map has distances in LD units (LDUs) that are analogous to centimorgans in linkage maps. The LD map of BTA6 has a total length of 8.9 LDUs. Within the LD map, regions of high LD (represented as blocks) and regions of low LD (steps) are observed, when plotted against the integrated map in kilobases. At the most stringent block definition, namely a set of loci with zero LDU increase over the span of these markers, BTA6 comprises 40 blocks, accounting for 41% of the chromosome. At a slightly lower stringency of block definition (a set of loci covering a maximum of 0.2 LDUs on the LD map), up to 81% of BTA6 is spanned by 46 blocks and with 13 steps that are likely to reflect recombination hot spots. The mean swept radius (the distance over which LD is likely to be useful for mapping) is 13.3 Mb, confirming extensive LD in Holstein-Friesian dairy cattle, which makes such populations ideal for whole-genome association studies.     

Genome association studies of complex diseases by case-control designs

One way to perform linkage-disequilibrium (LD) mapping of genetic traits is to use single markers. Since dense marker maps-such as single-nucleotide polymorphism and high-resolution microsatellite maps-are available, it is natural and practical to generalize single-marker LD mapping to high-resolution haplotype or multiple-marker LD mapping. This article investigates high-resolution LD-mapping methods, for complex diseases, based on haplotype maps or microsatellite marker maps. The objective is to explore test statistics that combine information from haplotype blocks or multiple markers. Based on two coding methods, genotype coding and haplotype coding, Hotelling's T2 statistics TG and TH are proposed to test the association between a disease locus and two haplotype blocks or two markers. The validity of the two T2 statistics is proved by theoretical calculations. A statistic TC, an extension of the traditional chi2 method of comparing haplotype frequencies, is introduced by simply adding the chi2 test statistics of the two haplotype blocks together. The merit of the three methods is explored by calculation and comparison of power and of type I errors. In the presence of LD between the two blocks, the type I error of TC is higher than that of TH and TG, since TC ignores the correlation between the two blocks. For each of the three statistics, the power of using two haplotype blocks is higher than that of using only one haplotype block. By power comparison, we notice that TC has higher power than that of TH, and TH has higher power than that of TG. In the absence of LD between the two blocks, the power of TC is similar to that of TH and higher than that of TG. Hence, we advocate use of TH in the data analysis. In the presence of LD between the two blocks, TH takes into account the correlation between the two haplotype blocks and has a lower type I error and higher power than TG. Besides, the feasibility of the methods is shown by sample-size calculation.     

An integrated comparative map of the porcine X chromosome

The objectives of this study were to assign both microsatellite and gene-based markers on porcine chromosome X to two radiation hybrid (RH) panels and to develop a more extensive integrated map of SSC-X. Thirty-five microsatellite and 20 gene-based markers were assigned to T43RH, and 16 previously unreported microsatellite and 15 gene-based markers were added to IMpRH map. Of these, 30 microsatellite and 12 gene-based markers were common to both RH maps. Twenty-two gene-based markers were submitted to BLASTN analysis for identification of orthologues of genes on HSA-X. Single nucleotide polymorphisms (SNPs) were detected for 12 gene-based markers, and nine of these were placed on the genetic map. A total of 92 known loci are present on at least one porcine chromosome X map. Thirty-seven loci are present on all three maps; 31 loci are found on only one map. Location of 33 gene-based markers on the comprehensive map translates into an integrated comparative map that supports conservation of gene order between SSC-X and HSA-X. This integrated map will be valuable for selection of candidate genes for porcine quantitative trait loci (QTLs) that map to SSC-X.     

An exploration of sex-specific linkage disequilibrium on chromosome X in Caucasians from the COGA study

This paper explores the decay of linkage disequilibrium (LD) on the autosomes and chromosome X. The extent of marker-marker LD is important for both linkage and association studies. The analysis of the Caucasian sample from the Collaborative Study on the Genetics of Alcoholism study revealed the expected negative relationship between the magnitude of the marker-marker LD and distance (cM), with the male and female subgroups exhibiting similar patterns of LD. The observed extent of LD in females was less across the pseudoautosomal markers relative to the heterosomal region of chromosome X. Marked differences in LD patterns were also observed between chromosomes X and the 22 autosomes in both males and females.     

A collection of ordered tetranucleotide-repeat markers from the human genome. The Utah Marker Development Group.

A collection of 1,069 human PCR-based genetic markers has been developed, and their distribution over the 22 autosomes and the X chromosome has been determined. Each marker was developed around a short-tandem-repeat DNA sequence. The majority (85%) of the markers described here were selected to contain tetranucleotide repeats, because these repeats show better stability during PCR than do dinucleotide repeats. Linkage maps constructed from genotypes collected with these markers in four CEPH pedigrees (1331, 1332, 1362, and 884) covered 3,417 cM of the human genome. More than 600 of the loci revealed heterozygosities > .70. Overall, 444 loci were ordered, with odds > 100:1 against inversion of adjacent loci. The average distance between markers was 7.4 cM on the autosomes and 24.8 cM on the X chromosome. Likely locations (100:1 odds intervals) were assigned for the remaining 621 short-tandem-repeat polymorphisms, as well as for 160 other markers that are present on the framework maps published by the Cooperative Human Linkage Center. Four markers specific to the Y chromosome are also reported here. From our maps, 347 markers were chosen to define "index" maps for each of the 22 autosomes. The index markers detect loci with an average heterozygosity of .85 and cover 3,169 cM of the autosomes, with an average distance between markers of 9.2 cM. These polymorphic short tandem repeats will be highly useful as reagents for the ongoing genetic and physical mapping of the human genome and for characterization of genetic changes in cancer.     

Analysis of genome-wide linkage disequilibrium in the highly polyploid sugarcane

Linkage disequilibrium (LD) in crops, established by domestication and early breeding, can be a valuable basis for mapping the genome. We undertook an assessment of LD in sugarcane (Saccharum spp), characterized by one of the most complex crop genomes, with its high ploidy level (>or=8) and chromosome number (>100) as well as its interspecific origin. Using AFLP markers, we surveyed 1,537 polymorphisms among 72 modern sugarcane cultivars. We exploited information from available genetic maps to determine a relevant statistical threshold that discriminates marker associations due to linkage from other associations. LD is very common among closely linked markers and steadily decreases within a 0-30 cM window. Many instances of linked markers cannot be recognized due to the confounding effect of polyploidy. However, LD within a sample of cultivars appears as efficient as linkage analysis within a controlled progeny in terms of assigning markers to cosegregation groups. Saturating the genome coverage remains a challenge, but applying LD-based mapping within breeding programs will considerably speed up the localization of genes controlling important traits by making use of phenotypic information produced in the course of selection.     

The extent of linkage disequilibrium in beef cattle breeds using high-density SNP genotypes

Background

The extent of linkage disequilibrium (LD) between molecular markers impacts genome-wide association studies and implementation of genomic selection. The availability of high-density single nucleotide polymorphism (SNP) genotyping platforms makes it possible to investigate LD at an unprecedented resolution. In this work, we characterised LD decay in breeds of beef cattle of taurine, indicine and composite origins and explored its variation across autosomes and the X chromosome.

Findings

In each breed, LD decayed rapidly and r² was less than 0.2 for marker pairs separated by 50 kb. The LD decay curves clustered into three groups of similar LD decay that distinguished the three main cattle types. At short distances between markers (< 10 kb), taurine breeds showed higher LD (r² = 0.45) than their indicine (r² = 0.25) and composite (r² = 0.32) counterparts. This higher LD in taurine breeds was attributed to a smaller effective population size and a stronger bottleneck during breed formation. Using all SNPs on only the X chromosome, the three cattle types could still be distinguished. However for taurine breeds, the LD decay on the X chromosome was much faster and the background level much lower than for indicine breeds and composite populations. When using only SNPs that were polymorphic in all breeds, the analysis of the X chromosome mimicked that of the autosomes.

Conclusions

The pattern of LD mirrored some aspects of the history of breed populations and showed a sharp decay with increasing physical distance between markers. We conclude that the availability of the HD chip can be used to detect association signals that remained hidden when using lower density genotyping platforms, since LD dropped below 0.2 at distances of 50 kb.     

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.     

Populationsgenetik des humanen X-Chromosoms

Mitochondrial DNA and the Y chromosome (ChrY) are both highly informative regarding human evolution, demographic history, and the genetic relationships between extant populations. The major reason for this is that both genomic compartments do not recombine, except for the pseudo-autosomal regions of ChrY, and that typing of haploid markers automatically allows the identification of haplotypes. In terms of its recombination behaviour, the X chromosome (ChrX) falls between autosomes and ChrY. The significance of ChrX in terms of population genetics is partially based on the fact that its haplotypes are easier to determine than those of autosomes. While ChrY and mtDNA each represent a single locus only, with a common evolutionary history of all their constituents, ChrX comprises several regions that may each reflect its own history. Therefore, ChrX studies seem most suitable for distinguishing between subpopulations or for research on regional ethnic structures. The analysis of linkage disequilibrium (LD) is one of the key aspects of population genetics studies of ChrX markers because LD may be an indicator of genetic isolation or of the emergence from a small founder population. Populations with high LD play an important role in the identification of genes involved in the aetiology of multifactorial diseases.     

Recombination is suppressed over a large region of the rainbow trout Y chromosome

The previous genetic mapping data have suggested that most of the rainbow trout sex chromosome pair is pseudoautosomal, with very small X-specific and Y-specific regions. We have prepared an updated genetic and cytogenetic map of the male rainbow trout sex linkage group. Selected sex-linked markers spanning the X chromosome of the female genetic map have been mapped cytogenetically in normal males and genetically in crosses between the OSU female clonal line and four different male clonal lines as well as in outcrosses involving outbred OSU and hybrids between the OSU line and the male clonal lines. The cytogenetic maps of the X and Y chromosomes were very similar to the female genetic map for the X chromosome. Five markers on the male maps are genetically very close to the sex determination locus ( SEX ), but more widely spaced on the female genetic map and on the cytogenetic map, indicating a large region of suppressed recombination on the Y chromosome surrounding the SEX locus. The male map is greatly extended at the telomere. A BAC clone containing the SCAR (sequence characterized amplified region) Omy - 163 marker, which maps close to SEX , was subjected to shotgun sequencing. Two carbonyl reductase genes and a gene homologous to the vertebrate skeletal ryanodine receptor were identified. Carbonyl reductase is a key enzyme involved in production of trout ovarian maturation hormone. This brings the number of type I genes mapped to the sex chromosome to six and has allowed us to identify a region on zebrafish chromosome 10 and medaka chromosome 13 which may be homologous to the distal portion of the long arm of the rainbow trout Y chromosome.     

Positional cloning by linkage disequilibrium

Recently, metric linkage disequilibrium (LD) maps that assign an LD unit (LDU) location for each marker have been developed (Maniatis et al. 2002). Here we present a multiple pairwise method for positional cloning by LD within a composite likelihood framework and investigate the operating characteristics of maps in physical units (kb) and LDU for two bodies of data (Daly et al. 2001; Jeffreys et al. 2001) on which current ideas of blocks are based. False-negative indications of a disease locus (type II error) were examined by selecting one single-nucleotide polymorphism (SNP) at a time as causal and taking its allelic count (0, 1, or 2, for the three genotypes) as a pseudophenotype, Y. By use of regression and correlation, association between every pseudophenotype and the allelic count of each SNP locus (X) was based on an adaptation of the Malecot model, which includes a parameter for location of the putative gene. By expressing locations in kb or LDU, greater power for localization was observed when the LDU map was fitted. The efficiency of the kb map, relative to the LDU map, to describe LD varied from a maximum of 0.87 to a minimum of 0.36, with a mean of 0.62. False-positive indications of a disease locus (type I error) were examined by simulating an unlinked causal SNP and the allele count was used as a pseudophenotype. The type I error was in good agreement with Wald's likelihood theorem for both metrics and all models that were tested. Unlike tests that select only the most significant marker, haplotype, or haploset, these methods are robust to large numbers of markers in a candidate region. Contrary to predictions from tagging SNPs that retain haplotype diversity, the sample with smaller size but greater SNP density gave less error. The locations of causal SNPs were estimated with the same precision in blocks and steps, suggesting that block definition may be less useful than anticipated for mapping a causal SNP. These results provide a guide to efficient positional cloning by SNPs and a benchmark against which the power of positional cloning by haplotype-based alternatives may be measured.     

The extent of linkage disequilibrium in four populations with distinct demographic histories

The design and feasibility of whole-genome-association studies are critically dependent on the extent of linkage disequilibrium (LD) between markers. Although there has been extensive theoretical discussion of this, few empirical data exist. The authors have determined the extent of LD among 38 biallelic markers with minor allele frequencies >.1, since these are most comparable to the common disease-susceptibility polymorphisms that association studies aim to detect. The markers come from three chromosomal regions-1,335 kb on chromosome 13q12-13, 380 kb on chromosome 19q13.2, and 120 kb on chromosome 22q13.3-which have been extensively mapped. These markers were examined in approximately 1,600 individuals from four populations, all of European origin but with different demographic histories; Afrikaners, Ashkenazim, Finns, and East Anglian British. There are few differences, either in allele frequencies or in LD, among the populations studied. A similar inverse relationship was found between LD and distance in each genomic region and in each population. Mean D' is.68 for marker pairs <5 kb apart and is.24 for pairs separated by 10-20 kb, and the level of LD is not different from that seen in unlinked marker pairs separated by >500 kb. However, only 50% of marker pairs at distances <5 kb display sufficient LD (delta>.3) to be useful in association studies. Results of the present study, if representative of the whole genome, suggest that a whole-genome scan searching for common disease-susceptibility alleles would require markers spaced < or = 5 kb apart.     

Integration of the PiGMaP and USD A maps for porcine chromosome 14

In order to align two previously published genetic linkage maps, a set of four of the United States Department of Agriculture (USDA) microsatellite linkage markers was mapped in the International Pig Gene Mapping Project (PiGMaP) reference families. Two-point linkage analysis was used between these USDA markers and the set of genes and markers previously mapped on the PiGMaP chromosome 14 map-Markers with threshold lod scores of three or greater were used for multipoint map construction. The USDA and PigGMaP linkage maps of chromosome 14 were aligned using the four USDA microsatellite markers along with three markers that are common to both maps. The PiGMaP genetic linkage map order for chromosome 14 was confirmed and the map was expanded to 193 cM with addition of the new markers.     

Genetic localization and sequential electrophoresis of glucose-6-phosphate dehydrogenase in Drosophila melanogaster

Studies were undertaken to investigate two critical aspects of the glucose-6-phosphate dehydrogenase polymorphism in Drosophila melanogaster. The first investigation unequivocally maps the genetic site of the G6PD locus to the X chromosome. The second study subjects a set of isochromosomal lines to sequential electrophoresis in an attempt to uncover common molecular heterogeneity within the global polymorphism, assuming that this variation may have gone undetected under conventional electrophoretic conditions. The genetic site was mapped following the segregation of the two common electrophoretic alleles, a so-called null allele, and two rare electrophoretic variants. From the pooled results, the Zw locus mapped to 62.9 on the X chromosome relative to the flanking markers car (at 62.5) and sw (at 64.7). A set of 126 iso-X chromosomal lines of diverse geographic origin was subjected to sequential electrophoresis under three different acrylamide conditions in addition to the conventional starch electrophoretic system. No additional variation beyond the common diallele polymorphism was seen.     

Mapping of 262 DNA markers into 24 intervals on human chromosome 11.

We have extended our mapping effort on human chromosome 11 to encompass a total of 262 DNA markers, which have been mapped into 24 intervals on chromosome 11; 123 of the markers reveal RFLPs. These clones are scattered throughout the chromosome, although some clustering occurs in R-positive bands (p15.1, p11.2, q13, and q23.3). Fifty-two of the markers were found to contain DNA sequences conserved in Chinese hamster, and some of these 52 also cross-hybridized with DNA from other mammals and/or chicken. As the length of chromosome 11 is estimated at nearly 130 cM, the average distance between RFLP markers is roughly 1 cM. The large panel of DNA markers on our map should contribute to investigations of hereditary diseases on this chromosome, and it will also provide reagents for constructing either fine-scale linkage and physical maps or contig maps of cosmids or yeast artificial chromosomes.     

Comparative mapping of human chromosome 10 to pig chromosomes 10 and 14

Identification of predictive markers in QTL regions that impact production traits in commercial populations of swine is dependent on construction of dense comparative maps with human and mouse genomes. Chromosomal painting in swine suggests that large genomic blocks are conserved between pig and human, while mapping of individual genes reveals that gene order can be quite divergent. High-resolution comparative maps in regions affecting traits of interest are necessary for selection of positional candidate genes to evaluate nucleotide variation causing phenotypic differences. The objective of this study was to construct an ordered comparative map of human chromosome 10 and pig chromosomes 10 and 14. As a large portion of both pig chromosomes are represented by HSA10, genes at regularly spaced intervals along this chromosome were targeted for placement in the porcine genome. A total of 29 genes from human chromosome 10 were mapped to porcine chromosomes 10 (SSC10) and 14 (SSC14) averaging about 5 Mb distance of human DNA per marker. Eighteen genes were assigned by linkage in the MARC mapping population, five genes were physically assigned with the IMpRH mapping panel and seven genes were assigned on both maps. Seventeen genes from human 10p mapped to SSC10, and 12 genes from human 10q mapped to SSC14. Comparative maps of mammalian species indicate that chromosomal segments are conserved across several species and represent syntenic blocks with distinct breakpoints. Development of comparative maps containing several species should reveal conserved syntenic blocks that will allow us to better define QTL regions in livestock.     

Mapping genes for common diseases: the case for genetic (LD) maps

We examine the current effort to develop a haplotype map of the human genome and suggest an alternative approach which represents linkage disequilibrium patterns in the form of a metric LD map. LD maps have some of the useful properties of genetic linkage maps but have a much higher resolution which is optimal for SNP-based association mapping of common diseases. The studies that have been undertaken to date suggest that LD and recombination maps show some close similarities because of abundant, narrow, recombination hot spots. These hot spots are co-localised in all populations but, unlike linkage maps, LD maps differ in scale for different populations because of differences in population history. The prospects for developing optimized panels of SNPs and the use of linkage disequilibrium maps in disease gene localisation are assessed in the light of recent evidence.     

A comparative radiation hybrid map of sheep chromosome 10

Comparative radiation hybrid (RH) maps of individual ovine chromosomes are essential to identify genes governing traits of economic importance in sheep, a livestock species for which whole genome sequence data are not yet available. The USUoRH5000 radiation hybrid panel was used to generate a RH map of sheep chromosome 10 (OAR10) with 59 markers that span 1,422 cR over an estimated 92 Mb of the chromosome, thus providing markers every 2 Mb (equivalent to every 24 cR). The markers were derived from 46 BAC end sequences (BESs), a single EST, and 12 microsatellites. Comparative analysis showed that OAR10 shares remarkable conservation of gene order along the entire length of cattle chromosome 12 and that OAR10 contains four major homologous synteny blocks, each related to segments of the homologous human chromosome 13. Extending the comparison to the horse, dog, mouse, and chicken genome showed that these blocks share conserved synteny across species.     

The variable number of tandem repeats element in DAT1 regulates in vitro dopamine transporter density

Background

The selection of markers in association studies can be informed through the use of haplotype blocks. Recent reports have determined the genomic architecture of chromosomal segments through different haplotype block definitions based on linkage disequilibrium (LD) measures or haplotype diversity criteria. The relative applicability of distinct block definitions to association studies, however, remains unclear. We compared different block definitions in 6.1 Mb of chromosome 17q in 189 unrelated healthy individuals. Using 137 single nucleotide polymorphisms (SNPs), at a median spacing of 15.5 kb, we constructed haplotype block maps using published methods and additional methods we have developed. Haplotype tagging SNPs (htSNPs) were identified for each map.

Results

Blocks were found to be shorter and coverage of the region limited with methods based on LD measures, compared to the method based on haplotype diversity. Although the distribution of blocks was highly variable, the number of SNPs that needed to be typed in order to capture the maximum number of haplotypes was consistent.

Conclusion

For the marker spacing used in this study, choice of block definition is not important when used as an initial screen of the region to identify htSNPs. However, choice of block definition has consequences for the downstream interpretation of association study results.     

Comparison of type I error for multiple test corrections in large single-nucleotide polymorphism studies using principal components versus haplotype blocking algorithms

Although permutation testing has been the gold standard for assessing significance levels in studies using multiple markers, it is time-consuming. A Bonferroni correction to the nominal p-value that uses the underlying pair-wise linkage disequilibrium (LD) structure among the markers to determine the number of effectively independent tests has recently been proposed. We propose using the number of independent LD blocks plus the number of independent single-nucleotide polymorphisms for correction. Using the Collaborative Study on the Genetics of Alcoholism LD data for chromosome 21, we simulated 1,000 replicates of parent-child trio data under the null hypothesis with two levels of LD: moderate and high. Assuming haplotype blocks were independent, we calculated the number of independent statistical tests using 3 haplotype blocking algorithms. We then compared the type I error rates using a principal components-based method, the three blocking methods, a traditional Bonferroni correction, and the unadjusted p-values obtained from FBAT. Under high LD conditions, the PC method and one of the blocking methods were slightly conservative, whereas the 2 other blocking methods exceeded the target type I error rate. Under conditions of moderate LD, we show that the blocking algorithm corrections are closest to the desired type I error, although still slightly conservative, with the principal components-based method being almost as conservative as the traditional Bonferroni correction.