Juxtaposed regions of extensive and minimal linkage disequilibrium in human Xq25 and Xq28

Sequence analysis of the mannose-binding lectin (MBL2) gene reveals a high degree of heterozygosity with evidence of selection

Human mannose-binding protein (MBL) is a component of innate immunity. To capture the common genetic variants of MBL2, we resequenced a 10.0 kb region that includes MBL2 in 102 individuals representing four major US ethnic groups. In all, 87 polymorphic sites were observed, indicating a high level of heterozygosity (total pi=18.3 x 10(-4)). Estimates of linkage disequilibrium across MBL2 indicate that it is divided into two blocks, with a probable recombination hot spot in the 3' end. Three non-synonymous SNPs in exon 1 of the encoding MBL2 gene and three upstream SNPs form common 'secretor haplotypes' that can predict circulating levels. Common variants have been associated with increased susceptibility to infection and autoimmune diseases. The high frequencies of B, C and D alleles in certain populations suggest a possible selective advantage for heterozygosity. There is limited diversity of haplotype structure; the 'secretor haplotypes' lie on a restricted number of extended haplotypes, which could include additional linked SNPs, which might also have possible functional implications. There is evidence for gene conversion in the region between the two blocks, in the last exon. Our data should form the basis for conducting MBL2 candidate gene association studies using a locus-wide approach.

Evaluation of microsatellite markers in association studies: a search for an immune-related susceptibility gene in sarcoidosis

Association studies using linkage disequilibrium (LD) between candidate loci and nearby markers have been proposed to identify susceptibility genes for complex diseases. We analyzed polymorphisms of microsatellites (MSs) and LD patterns of the regions in which candidate genes related to the Th1 immune response have been annotated and attempted to identify a susceptibility gene for sarcoidosis in a marker-based association study. Nineteen MSs were identified in six Th1-related genes (IFNGR1, IFNGR2, IL12RB1, IL12RB2, STAT1 and STAT4) and then eight were further characterized as useful polymorphic markers. Most of these MSs showed LD with single nucleotide polymorphisms (SNPs) on both 5' and 3' ends of these candidate genes, in which r(2) values between at least one of the MS marker alleles and the SNPs were higher than 0.1. A significant association with one MS allele near STAT4 was shown and a cluster of SNPs in LD with the MS marker was associated with sarcoidosis. These results suggest that association studies using not only SNPs but also multi-allelic MS within or near candidate loci would be useful markers to search for a disease susceptibility gene, especially in populations with unknown LD structure.

Association between evolutionary history of angiotensinogen haplotypes and plasma levels

Over the last decade, considerable effort has been invested in studying the associations between angiotensinogen (AGT) variants, AGT plasma levels and high blood pressure. Evidence accumulated to date consistently supports the relationship between the AGT locus and the protein level, while an influence on blood pressure has been difficult to establish; in both instances the predisposing molecular variants are not fully defined. An evolutionary approach, taking into account the phylogenetic relationship between all the polymorphisms at this locus, may improve our understanding of the genetic nature of these quantitative phenotypes. Accordingly we sequenced a 6.8 kb region of the AGT gene in 57 Nigerian individuals (29 with high AGT plasma levels and 28 with low AGT plasma levels). Haplotypes were grouped into seven major haplogroups and their phylogenetic relationship was established. The association between haplogroups and AGT plasma levels was investigated. A significant linear correlation was detected between haplogroup genetic distance and AGT levels, suggesting a nonrandom accumulation of risk-associated mutations during the evolutionary history of the AGT gene.

Distribution of genome-wide linkage disequilibrium based on microsatellite loci in the Samoan population

Whole genome-wide scanning for susceptibility loci based on linkage disequilibrium (LD) has been proposed as a powerful strategy for mapping common complex diseases, especially in isolated populations. We recruited 389 individuals from 175 families in the US territory of American Samoa, and 96 unrelated individuals from American Samoa and the independent country of Samoa in order to examine background LD by using a 10 centimorgan (cM) map containing 381 autosomal and 18 X-linked microsatellite markers. We tested the relationship between LD and recombination fraction by fitting a regression model. We estimated a slope of -0.021 (SE 0.00354; p < 0.0001). Based on our results, LD in the Samoan population decays steadily as the recombination fraction between autosomal markers increases. The patterns of LD observed in the Samoan population are quite similar to those previously observed in Palau but markedly contrast with those observed in a non-isolated Caucasian sample, where there is essentially no marker-to-marker LD. Our analyses support the hypothesis of a recent bottleneck, which is consistent with the known demographic history of the Samoan population. Furthermore, population substructure tests support the hypothesis that self-identified Samoans represent one homogenous genetic population.

Neither single-marker nor haplotype analyses support an association between genetic variation near NOTCH4 and bipolar disorder

Markers near the NOTCH4 locus on chromosome 6p21.3 have been reported to be associated with schizophrenia in some studies. Since schizophrenia and bipolar affective disorder (BPAD) may share genetic determinants, we tested markers in and near NOTCH4 in a sample of 153 parent-offspring triads ascertained through a sibling pair with BPAD for evidence of association. This sample would have 80% power to detect an association at or above a genotype relative risk of 2.4 at the 10(-7) level of significance. In addition to the two markers previously showing the most significant association with schizophrenia, three additional nearby markers were studied. The five markers were genotyped using validated methods. Both single-marker and 3-marker haplotype data was analyzed using family-based association methods. No genome-wide significant association was detected between any of the five SNP-markers and BPAD in this sample. One marker showed nominal evidence of association (P = 0.049), but this evidence was not supported by haplotype analyses including nearby flanking markers or by case-control analysis using 93 Caucasian controls. These results do not support an association between genetic variation near NOTCH4 and BPAD in this sample.

Genetic studies of neuropsychiatric disorders in Costa Rica: a model for the use of isolated populations

The importance of genetics in understanding the etiology of mental illness has become increasingly clear in recent years, as more evidence has mounted that almost all neuropsychiatric disorders have a genetic component. It has also become clear, however, that these disorders are etiologically complex, and multiple genetic and environmental factors contribute to their makeup. So far, traditional linkage mapping studies have not definitively identified specific disease genes for neuropsychiatric disorders, although some potential candidates have been identified via these methods (e.g. the dysbindin gene in schizophrenia; Straub et al., 2002; Schwab et al., 2003). For this reason, alternative approaches are being attempted, including studies in genetically isolated populations. Because isolated populations have a high degree of genetic homogeneity, their use may simplify the process of identifying disease genes in disorders where multiple genes may play a role. Several areas of Latin America contain genetically isolated populations that are well suited for the study of neuropsychiatric disorders. Genetic studies of several major psychiatric illnesses, including bipolar disorder, major depression, schizophrenia, Tourette Syndrome, alcohol dependence, attention deficit hyperactivity disorder, and obsessive-compulsive disorder, are currently underway in these regions. In this paper we highlight the studies currently being conducted by our groups in the Central Valley of Costa Rica to illustrate the potential advantages of this population for genetic studies.

Genetic isolates in Corsica (France): linkage disequilibrium extension analysis on the Xq13 region

Genetic isolates with a history of a small founder population, long-lasting isolation and population bottlenecks represent exceptional resources in the identification of genes involved in the pathogenesis of multifactorial diseases. In these populations, the disease allele reveals linkage disequilibrium (LD) with markers over significant genetic intervals, therefore facilitating disease locus identification. This study has been designed to examine the background LD extension in some subpopulations of Corsica. Our interest in the island of Corsica is due to its geographical and genetic proximity to the other Mediterranean island of Sardinia. Sardinian isolates in which the extension of the background LD is particularly high have been recently identified and are now the object of studies aimed at the mapping of genes involved in complex diseases. Recent evidence has highlighted that the genetic proximity between the populations of Corsica and Sardinia is particularly true for the internal conservative populations. Given these considerations, Sardinia and Corsica may represent a unique system to carry out parallel association studies whose results could be validated by comparison. In the present study, we have analyzed the LD extension on the Xq13 genomic region in three subpopulations of Corsica: Corte, Niolo and Bozio, all located in the mountainous north-center of the island. Our results show a strong degree of LD over long distance for the population of Bozio and to a less extent for the population of Niolo. Their LD extent is comparable to or higher than that reported for other isolates.

Comparative study of the haplotype structure and linkage disequilibrium of chromosome 1p36.2 region in the Korean and Japanese populations

The patterns of linkage disequilibrium (LD) in the human genome provide important information for disease gene mapping. LDs may vary depending on chromosomal regions and populations. We have compared LD and haplotypes defined by SNPs in the chromosome 1p36.2 region of the Korean and Japanese populations. Fifty-eight SNPs in about 418 kb ranging from tumor necrosis factor receptor 2 (TNFR2:TNFRSF1B) to procollagen-lysine, 2-oxoglutarate 5-dioxygenase (PLOD) gene were examined in 96 healthy Koreans and Japanese each by direct sequencing and fluorescence correlation spectroscopy combined with the PCR-sequence specific primer method (PCR-SSP-FCS), respectively. Upon pair-wise LD analysis, a total of 25 and 16 out of 58 SNPs greater than MAF 10% were included in LD blocks, encompassing almost 81 kb and 55 kb in total, in Koreans and Japanese, respectively. Both similarities and differences were observed in LD strength and haplotype frequencies between the populations. Considerable similarities were observed in the telomeric region where a long-range block of approximately 80 kb including three genes was found to have strong LDs in both Koreans and Japanese. Significant difference in LD strength was present near the TNFR2 region between the Japanese and Korean populations.

Association of multiple nucleotide variations in the pituitary glutaminyl cyclase gene (QPCT) with low radial BMD in adult women

Correlation between 13 genetic variations of the glutaminyl-peptide cyclotransferase gene and adjusted aBMD was tested among 384 adult women. Among 13 variations with strong linkage disequilibrium, R54W showed a prominent association (p = 0.0003), which was more striking when examined among 309 elder subjects (> or =50 years; p = 0.0001). Contribution for postmenopausal bone loss was suggested.

INTRODUCTION

Alterations in homeostatic regulation of estrogen through the hypothalamus-pituitary-gonadal axis (HPG axis) importantly affect the pathogenesis of osteoporosis. Osteoporosis-susceptibility genes have been proposed in this hormonal axis, such as estrogen receptor genes and the gonadotropin-releasing hormone gene (GnRH). Here we report another example of genes: glutaminyl-peptide cyclotransferase gene (QPCT), an essential modifier of pituitary peptide hormones, including GnRH.

MATERIALS AND METHODS

Analyses of association of 13 single nucleotide polymorphisms (SNPs) at the QPCT locus with adjusted areal BMD (adj-aBMD) were carried out among 384 adult women. Linkage disequilibrium (LD) was analyzed by haplotype estimation and calculation of D' and r2. Multiple regression analysis was applied for evaluating the combined effects of the variations.

RESULTS AND CONCLUSIONS

LD analysis indicated strong linkage disequilibrium within the entire 30-kb region of the QPCT gene. Significant correlations were observed between the genotypes of the six SNPs and the radial adj-aBMD, among which R54W (nt + 160C>T) presented the most prominent association (p = 0.0003). Striking association was observed for these SNPs among the 309 subjects >50 years of age (R54W, p = 0.0001; -1095T>C, p = 0.0002; -1844C>T, p = 0.0002). Multiple regression analyses indicated that multiple SNPs in the gene might act in combination to determine the radial adj-aBMD. These results indicate that genetic variations in QPCT are the important factors affecting the BMD of adult women that contribute to susceptibility for osteoporosis. The data should provide new insight into the etiology of the disease and may suggest a new target to be considered during treatment.

Global genetic analysis

The introduction of molecular markers in genetic analysis has revolutionized medicine. These molecular markers are genetic variations associated with a predisposition to common diseases and individual variations in drug responses. Identification and genotyping a vast number of genetic polymorphisms in large populations are increasingly important for disease gene identification, pharmacogenetics and population-based studies. Among variations being analyzed, single nucleotide polymorphisms seem to be most useful in large-scale genetic analysis. This review discusses approaches for genetic analysis, use of different markers, and emerging technologies for large-scale genetic analysis where millions of genotyping need to be performed.

DXS10011: studies on structure, allele distribution in three populations and genetic linkage to further q-telomeric chromosome X markers

The hypervariable tetranucleotide STR polymorphism DXS10011 is a powerful marker for forensic purposes. Investigation of this STR led to an allele nomenclature which is in consensus with the ISFG recommendations. DXS10011 is located at Xq28 and genetically closely linked to DXS7423 and DXS8377 but is unlinked to HPRTB and more distant X-chromosomal STRs. DXS10011 is a very complex marker exhibiting some structural variants within alleles of identical length. Two types of repeat structure (regular and inter-alleles) are known and described as types A and B. Two SNPs which are in strong linkage disequilibrium to the different sequence types were found in the repeat flanking region. The type A sequence consists of a long stretch of uninterrupted homogenous repeats which is highly susceptible to slippage mutation during male meiosis.

Haplotypic relationship between SNP and microsatellite markers at the NOS2A locus in two populations

The density of genetic markers required for successful association mapping of complex diseases depends on linkage disequilibrium (LD) between non-functional markers and functional variants. The haplotypic relationship between stable markers and potentially unstable but highly informative markers (e.g. microsatellites) indicates that LD might be maintained over considerable genetic distance in non-African populations, supporting the use of such 'mixed marker haplotypes' in LD-based mapping, and allowing inferences to be drawn about human origins. We investigated sequence variation in the proximal 2.6 kb of the inducible nitric oxide synthase (NOS2A) promoter and the relationship between SNP haplotypes and a pentanucleotide microsatellite (the 'NOS2A(-2.6) microsatellite') in Gambians and UK Caucasians. UK Caucasians exhibited a subset of sequence diversity observed in Gambians, sharing four of 11 SNPs and a similar haplotypic structure. Five SNPs were found in the sequence of interspersed repetitive DNA elements. In both populations, there was dramatic loss of LD between SNP haplotypes and microsatellite alleles across a very short physical distance, suggesting a high intrinsic mutation rate of the NOS2A(-2.6) microsatellite, the SNP haplotypes are relatively ancient, or that this was a region of frequent recombination. Understanding locus- and population-specific LD is essential when designing and interpreting genetic association studies.

Positive associations of polymorphisms in the metabotropic glutamate receptor type 8 gene (GRM8) with schizophrenia

The glutamatergic dysfunction has been implicated in pathophysiology of schizophrenia. The Group III metabotropic glutamate receptor 4 (mGluR4), 6, 7, and 8 are thought to modulate glutamatergic transmission in the brain by inhibiting glutamate release at the synapse. We tested association of schizophrenia with GRM8 using 22 single nucleofide polymorphisms (SNPs) with the average intervals of 40.3 kb in the GRM8 region in 100 case-control pairs for the SNPs. Although we observed significant associations of schizophrenia with two SNPs, SNP18 (rs2237748, allele: P = 0.0279; genotype: P = 0.0124) and SNP19 (rs2299472, allele: P = 0.0302; genotype: P = 0.0127), none of two SNPs showed significant association with disease after Bonferroni correction. Both SNP18 and SNP19 were included in a large region (>330 kb) in which SNPs are in linkage disequilibrium (LD) at the 3' region of GRM8. We also tested haplotype association of schizophrenia with constructed haplotypes of the SNPs in LD. Significant associations were detected for the combinations of SNP5-SNP6 (chi(2) = 18.12, df = 3, P = 0.0004, P corr = 0.0924 with Bonferroni correction), SNP4-SNP5-SNP6 (chi(2) = 27.50, df = 7, P = 0.0075, P corr = 0.015 with Bonferroni correction), and SNP5-SNP6-SNP7 (chi(2) = 23.92, df = 7, P = 0.0011, P corr = 0.0022 with Bonferroni correction). Thus, we conclude that at least one susceptibility locus for schizophrenia is located within the GRM8 region in Japanese.

Patterns of linkage disequilibrium and haplotype distribution in disease candidate genes

Background

The adequacy of association studies for complex diseases depends critically on the existence of linkage disequilibrium (LD) between functional alleles and surrounding SNP markers.

Results

We examined the patterns of LD and haplotype distribution in eight candidate genes for osteoporosis and/or obesity using 31 SNPs in 1,873 subjects. These eight genes are apolipoprotein E (APOE), type I collagen α1 (COL1A1), estrogen receptor-α (ER-α), leptin receptor (LEPR), parathyroid hormone (PTH)/PTH-related peptide receptor type 1 (PTHR1), transforming growth factor-β1 (TGF-β1), uncoupling protein 3 (UCP3), and vitamin D (1,25-dihydroxyvitamin D₃) receptor (VDR). Yin yang haplotypes, two high-frequency haplotypes composed of completely mismatching SNP alleles, were examined. To quantify LD patterns, two common measures of LD, D' and r², were calculated for the SNPs within the genes. The haplotype distribution varied in the different genes. Yin yang haplotypes were observed only in PTHR1 and UCP3. D' ranged from 0.020 to 1.000 with the average of 0.475, whereas the average r² was 0.158 (ranging from 0.000 to 0.883). A decay of LD was observed as the intermarker distance increased, however, there was a great difference in LD characteristics of different genes or even in different regions within gene.

Conclusion

The differences in haplotype distributions and LD patterns among the genes underscore the importance of characterizing genomic regions of interest prior to association studies.

Haplotype block partitioning and tag SNP selection using genotype data and their applications to association studies

Recent studies have revealed that linkage disequilibrium (LD) patterns vary across the human genome with some regions of high LD interspersed by regions of low LD. A small fraction of SNPs (tag SNPs) is sufficient to capture most of the haplotype structure of the human genome. In this paper, we develop a method to partition haplotypes into blocks and to identify tag SNPs based on genotype data by combining a dynamic programming algorithm for haplotype block partitioning and tag SNP selection based on haplotype data with a variation of the expectation maximization (EM) algorithm for haplotype inference. We assess the effects of using either haplotype or genotype data in haplotype block identification and tag SNP selection as a function of several factors, including sample size, density or number of SNPs studied, allele frequencies, fraction of missing data, and genotyping error rate, using extensive simulations. We find that a modest number of haplotype or genotype samples will result in consistent block partitions and tag SNP selection. The power of association studies based on tag SNPs using genotype data is similar to that using haplotype data.

Quantitative Analysis of Nucleic Acids - the Last Few Years of Progress

DNA and RNA quantifications are widely used in biological and biomedical research. In the last ten years, many technologies have been developed to enable automated and high-throughput analyses. In this review, we first give a brief overview of how DNA and RNA quantifications are carried out. Then, five technologies (microarrays, SAGE, differential display, real time PCR and real competitive PCR) are introduced, with an emphasis on how these technologies can be applied and what their limitations are. The technologies are also evaluated in terms of a few key aspects of nucleic acids quantification such as accuracy, sensitivity, specificity, cost and throughput.

The International HapMap Project

The goal of the International HapMap Project is to determine the common patterns of DNA sequence variation in the human genome and to make this information freely available in the public domain. An international consortium is developing a map of these patterns across the genome by determining the genotypes of one million or more sequence variants, their frequencies and the degree of association between them, in DNA samples from populations with ancestry from parts of Africa, Asia and Europe. The HapMap will allow the discovery of sequence variants that affect common disease, will facilitate development of diagnostic tools, and will enhance our ability to choose targets for therapeutic intervention.

Defining haplotype blocks and tag single-nucleotide polymorphisms in the human genome

Recent studies suggest that the genome is organized into blocks of haplotypes, and efforts to create a genome-wide haplotype map of single-nucleotide polymorphisms (SNPs) are already underway. Haplotype blocks are defined algorithmically and to date several algorithms have been proposed. However, little is known about their relative performance in real data or about the impact of allele frequencies and parameter choices on the detection of haplotype blocks and the markers that tag them. Here we present a formal comparison of two major algorithms, a linkage disequilibrium (LD)-based method and a dynamic programming algorithm (DPA), in three chromosomal regions differing in gene content and recombination rate. The two methods produced strikingly different results. DPA identified fewer and larger haplotype blocks as well as a smaller set of tag SNPs than the LD method. For both methods, the results were strongly dependent on the allele frequency. Decreasing the minor allele frequency led to an up to 3.7-fold increase in the number of haplotype blocks and tag SNPs. Definition of haploytpe blocks and tag SNPs was also sensitive to parameter changes, but the results could not be reconciled simply by parameter adjustment. These results show that two major methods for detecting haplotype blocks and tag SNPs can produce different results in the same data and that these results are sensitive to marker allele frequencies and parameter choices. More information is needed to guide the choice of method, marker allele frequencies, and parameters in the development of a haplotype map.

Striking differentiation of sub-populations within a genetically homogeneous isolate (Ogliastra) in Sardinia as revealed by mtDNA analysis

Since the reduced genetic diversity found in isolates should simplify the study of complex traits, analyses of patterns of homogeneity within populations are of particular interest. We analysed the mtDNA haplogroups and hypervariable segment I (HVS-I) sequences of 475 individuals from a geographically restricted and isolated area (Ogliastra) within Sardinia, comprehending 175 random samples from 20 out of 23 villages. The remaining 300 subjects were chosen from the other three villages, Talana, Urzulei and Perdasdefogu, by sampling all maternal lineages. A comparison with other European populations reveals that Ogliastra ranks among the most genetically homogenous population and that it has been small and isolated throughout its history. The lack of variation and the high genetic homogeneity indicate that an important founder event and a demographic expansion took place during the Neolithic (approximately 7700 years before present) in Ogliastra's mtDNA gene pool. We present highly resolved phylogenetic networks for Ogliastra and for the three sub-isolates. MtDNA differentiation in the sub-populations versus Ogliastra is revealed by a strong demarcation in their genetic pools due to distinctive founder effects and genetic drift. We found that genetic homogeneity strictly depends on a scale factor in population size and on sampling methodology. The outstanding homogeneity and the reduced female gene pool observed in Ogliastra, in the European context, hide an extremely marked differentiation in sub-isolates originated from the same archaic population. Although Ogliastra can be considered a genetically homogeneous isolate, small villages' divergent genetic histories underline the importance of more systematic analysis of DNA variation between and within populations.

Protein kinase C/ζ (PRKCZ) Gene is associated with type 2 diabetes in Han population of North China and analysis of its haplotypes

AIM: To identify the susceptible gene (s) for type 2 diabetes in the prevousely mapped region, 1p36.33-p36.23, in Han population of North China using single nucleotide polymorphisms (SNPs) and to analyze the haplotypes of the gene (s) related to type 2 diabetes.

METHODS: Twenty three SNPs located in 10 candidate genes in the mapped region were chosen from public SNP domains with bioinformatic methods, and the single base extension (SBE) method was used to genotype the loci for 192 sporadic type 2 diabetes patients and 172 normal individuals, all with Han ethical origin, to perform this case-control study. The haplotypes with significant difference in the gene (s) were further analyzed.

RESULTS: Among the 23 SNPs, 8 were found to be common in Chinese Han population. Allele frequency of one SNP, rs436045 in the protein kinase C/ζgene (PRKCZ) was statistically different between the case and control groups (P < 0.05). Furthermore, haplotypes at five SNP sites of PRKCZ gene were identified.

CONCLUSION: PRKCZ gene may be associated with type 2 diabetes in Han population in North China. The haplotypes at five SNP sites in this gene may be responsible for this association.

Assessing the performance of the haplotype block model of linkage disequilibrium

Several recent studies have suggested that linkage disequilibrium (LD) in the human genome has a fundamentally "blocklike" structure. However, thus far there has been little formal assessment of how well the haplotype block model captures the underlying structure of LD. Here we propose quantitative criteria for assessing how blocklike LD is and apply these criteria to both real and simulated data. Analyses of several large data sets indicate that real data show a partial fit to the haplotype block model; some regions conform quite well, whereas others do not. Some improvement could be obtained by genotyping higher marker densities but not by increasing the number of samples. Nonetheless, although the real data are only moderately blocklike, our simulations indicate that, under a model of uniform recombination, the structure of LD would actually fit the block model much less well. Simulations of a model in which much of the recombination occurs in narrow hotspots provide a much better fit to the observed patterns of LD, suggesting that there is extensive fine-scale variation in recombination rates across the human genome.

Haplotypes and the systematic analysis of genetic variation in genes and genomes

Haplotypes have been used in various fields of genetics for a long time, in a variety of contexts, and for different purposes. Now, haplotype-based approaches to the analysis of candidate genes and genome-wide linkage disequilibrium (LD) mapping have gained center stage. It is time to explicitly distinguish the different concepts implied in the present haplotype approaches: haplotypes are not haplotypes, after all. The distinction of three different categories, ancestral, common haplotypes or haplotype blocks, gene-based haplotypes as complex genetic markers and gene-based functional haplotypes, is proposed. These categories serve as framework to review and analyze in particular the recent work suggesting evidence for a haplotype block structure of the human genome and the body of comparative sequencing studies addressing haplotype and LD structures at the gene level. Haplotype approaches will be evaluated along the dimensions preselection of variants versus complete DNA sequence information, role of LD and stages in the process of disease gene identification. Overall, the content of haplotypes is conceived as a function of available technologies to evaluate genetic variation and general advances in human genome research.

Association of a promoter haplotype (-1542G/-525C) in the tumor necrosis factor receptor associated factor-interacting protein gene with low bone mineral density in Japanese women

Osteoporosis, a multifactorial common disease, is believed to result from the interplay of multiple environmental and genetic factors that regulate bone mineral density (BMD). Tumor necrosis factor receptor associated factor-interacting protein (I-TRAF) is an essential effecter of the tumor necrosis factor receptor-signaling cascade, one of the most potent bone-resorbing systems. In genetic studies of 382 Japanese adult women, we found that genotypes of two promoter variations of I-TRAF gene, -1542T/G and -525G/C, were similarly associated with radial BMD levels. Two variations were in almost complete linkage disequilibrium (D' = 0.978, r(2) = 0.917, chi(2) = 695, 2, P = 3.4 x 10(-153)), and there were two exclusive haplotypes (-1542T/-525C, frequency 0.74, and -1542G/-525G, frequency 0.24) among our test subjects. When BMD values were compared among the three haplotypic categories (-1542G/-525G homozygotes, heterozygotes, and -1542T/-525C homozygotes), BMD was lowest among -1542G/-525C homozygotes (mean +/- SD = 0.382 +/- 0.060 g/cm(2)), highest among -1542T/-525G homozygotes (0.405 +/- 0.051 g/cm(2)), and intermediate among heterozygotes (0.395 +/- 0.056 g/cm(2)) (r = 0.11, P = 0.030). The observed trend supported a codominant effect of the relevant haplotype of I-TRAF gene in determination of radial BMD. These results suggested that variation of I-TRAF might be an important determinant for postmenopausal osteoporosis.

Linkage disequilibrium and inference of ancestral recombination in 538 single-nucleotide polymorphism clusters across the human genome

The prospect of using linkage disequilibrium (LD) for fine-scale mapping in humans has attracted considerable attention, and, during the validation of a set of single-nucleotide polymorphisms (SNPs) for linkage analysis, a set of data for 4,833 SNPs in 538 clusters was produced that provides a rich picture of local attributes of LD across the genome. LD estimates may be biased depending on the means by which SNPs are first identified, and a particular problem of ascertainment bias arises when SNPs identified in small heterogeneous panels are subsequently typed in larger population samples. Understanding and correcting ascertainment bias is essential for a useful quantitative assessment of the landscape of LD across the human genome. Heterogeneity in the population recombination rate, rho=4Nr, along the genome reflects how variable the density of markers will have to be for optimal coverage. We find that ascertainment-corrected rho varies along the genome by more than two orders of magnitude, implying great differences in the recombinational history of different portions of our genome. The distribution of rho is unimodal, and we show that this is compatible with a wide range of mixtures of hotspots in a background of variable recombination rate. Although rho is significantly correlated across the three population samples, some regions of the genome exhibit population-specific spikes or troughs in rho that are too large to be explained by sampling. This result is consistent with differences in the genealogical depth of local genomic regions, a finding that has direct bearing on the design and utility of LD mapping and on the National Institutes of Health HapMap project.

Haplotype blocks and linkage disequilibrium in the human genome

There is great interest in the patterns and extent of linkage disequilibrium (LD) in humans and other species. Characterizing LD is of central importance for gene-mapping studies and can provide insights into the biology of recombination and human demographic history. Here, we review recent developments in this field, including the recently proposed 'haplotype-block' model of LD. We describe some of the recent data in detail and compare the observed patterns to those seen in simulations.

Randomly distributed crossovers may generate block-like patterns of linkage disequilibrium: an act of genetic drift

There is considerable interest in identifying and characterizing block-like patterns of linkage disequilibrium (LD; haplotype blocks) in the human genome as these may facilitate the identification of complex disease genes via genome-wide association studies. Although recombination hot-spots have been suggested as the primary mechanism to explain the block-like pattern of LD, other forces, such as genetic drift, may also be important. To this end, we have studied the effect of various recombination models on patterns of LD by using extensive simulations. As expected, haplotype blocks were observed under a model allowing recombination hot-spots. However, we also observed similar block-like patterns in the models where recombination crossovers are randomly and uniformly distributed, and we demonstrate that these blocks are generated by genetic drift. We caution that genetic drift may be an alternative mechanism (in addition to recombination hot-spots) that can lead to block-like patterns of LD. Our findings highlight the necessity of characterizing haplotype blocks in world-wide populations.

Common 5' beta-globin RFLP haplotypes harbour a surprising level of ancestral sequence mosaicism

Blocks of linkage disequilibrium (LD) in the human genome represent segments of ancestral chromosomes. To investigate the relationship between LD and genealogy, we analysed diversity associated with restriction fragment length polymorphism (RFLP) haplotypes of the 5' beta-globin gene complex. Genealogical analyses were based on sequence alleles that spanned a 12.2-kb interval, covering 3.1 kb around the psibeta gene and 6.2 kb of the delta-globin gene and its 5' flanking sequence known as the R/T region. Diversity was sampled from a Kenyan Luo population where recent malarial selection has contributed to substantial LD. A single common sequence allele spanning the 12.2-kb interval exclusively identified the ancestral chromosome bearing the "Bantu" beta(s) (sickle-cell) RFLP haplotype. Other common 5' RFLP haplotypes comprised interspersed segments from multiple ancestral chromosomes. Nucleotide diversity was similar between psibeta and R/T-delta-globin but was non-uniformly distributed within the R/T-delta-globin region. High diversity associated with the 5' R/T identified two ancestral lineages that probably date back more than 2 million years. Within this genealogy, variation has been introduced into the 3' R/T by gene conversion from other ancestral chromosomes. Diversity in delta-globin was found to lead through parts of the main genealogy but to coalesce in a more recent ancestor. The well-known recombination hotspot is clearly restricted to the region 3' of delta-globin. Our analyses show that, whereas one common haplotype in a block of high LD represents a long segment from a single ancestral chromosome, others are mosaics of short segments from multiple ancestors related in genealogies of unsuspected complexity.

Linkage disequilibrium in human populations

Whereas the human linkage map appears on limited evidence to be constant over populations, maps of linkage disequilibrium (LD) vary among populations that differ in gene history. The greatest difference is between populations of sub-Saharan origin and populations remotely derived from Africa after a major bottleneck that reduced their heterozygosity and altered their Malecot parameters, increasing the intercept M that reflects association in founders and decreasing the exponential decline epsilon. Variation among populations within this ethnic dichotomy is much smaller. These observations validate use of a cosmopolitan LD map based on a sizeable sample representing a large population reliably typed for markers at high density. Then an LD map for a region or isolate within an ethnic group may be created by fitting the sample LD to the cosmopolitan map, estimating Malecot parameters simultaneously. The cosmopolitan map scaled by epsilon recovers 95% of the information that a local map at the same density gives and therefore more than the information in a low-resolution local map. Relative to a Eurasian cosmopolitan map the scaling factors are estimated to be 0.82 for isolates of European descent, 1.53 for Yorubans, and 1.74 for African Americans. These observations are consistent with a common bottleneck (perhaps but not necessarily speciation) approximately 173,500 years ago, if the bottleneck associated with migration out of Africa was 100,000 years ago. Eurasian populations (especially isolates with numerous cases) are efficient for genome scans, and populations of recent African origin (such as African Americans) are efficient for identification of causal polymorphisms within a candidate sequence.

Skin pigmentation, biogeographical ancestry and admixture mapping

Ancestry informative markers (AIMs) are genetic loci showing alleles with large frequency differences between populations. AIMs can be used to estimate biogeographical ancestry at the level of the population, subgroup (e.g. cases and controls) and individual. Ancestry estimates at both the subgroup and individual level can be directly instructive regarding the genetics of the phenotypes that differ qualitatively or in frequency between populations. These estimates can provide a compelling foundation for the use of admixture mapping (AM) methods to identify the genes underlying these traits. We present details of a panel of 34 AIMs and demonstrate how such studies can proceed, by using skin pigmentation as a model phenotype. We have genotyped these markers in two population samples with primarily African ancestry, viz. African Americans from Washington D.C. and an African Caribbean sample from Britain, and in a sample of European Americans from Pennsylvania. In the two African population samples, we observed significant correlations between estimates of individual ancestry and skin pigmentation as measured by reflectometry (R(2)=0.21, P<0.0001 for the African-American sample and R(2)=0.16, P<0.0001 for the British African-Caribbean sample). These correlations confirm the validity of the ancestry estimates and also indicate the high level of population structure related to admixture, a level that characterizes these populations and that is detectable by using other tests to identify genetic structure. We have also applied two methods of admixture mapping to test for the effects of three candidate genes (TYR, OCA2, MC1R) on pigmentation. We show that TYR and OCA2 have measurable effects on skin pigmentation differences between the west African and west European parental populations. This work indicates that it is possible to estimate the individual ancestry of a person based on DNA analysis with a reasonable number of well-defined genetic markers. The implications and applications of ancestry estimates in biomedical research are discussed.

Using haplotype blocks to map human complex trait loci

Understanding of linkage disequilibrium (LD) in human populations could facilitate the discovery of genes that influence complex human diseases. The "HapMap" project is now underway to characterize patterns of LD in the human genome. A pilot study showed "haplotype blocks" in 51 regions scattered throughout the genome. These intriguing results raise important questions about the nature of recombination, and highlight practical issues of marker collection, the influence of statistical modelling on apparent block structure, and the levels of genotyping necessary for studies of common diseases. Knowledge of local disequilibrium patterns may help identify common polymorphisms involved in complex disease, but completely new analytical methods and experimental designs will be required to identify important rare variants.

Association of a haplotype (196Phe/532Ser) in the interleukin-1-receptor-associated kinase (IRAK1) gene with low radial bone mineral density in two independent populations

Osteoporosis, a multifactorial common disease, is believed to result from the interplay of multiple environmental and genetic determinants, including factors that regulate bone mineral density. Interleukin-1 (IL-1) is one of the most potent bone-resorbing factors, and interleukin-1-associated kinase 1 (IRAK1) is an essential effector of the IL-1 receptor signaling cascade. In genetic studies of two independent populations of postmenopausal women (cohort A: 220 individuals and cohort T: 126 individuals) from separated geographical regions of Japan, we found that radial bone mineral density levels had similar associations with IRAK1 genotypes in both populations. Two amino acid-substituting variations in the gene, encoding Phe196Ser and Ser532Leu, were in complete linkage disequilibrium (D' = 1.0000, r2 = 1.0000, chi2 = 192.000, p = 1.2 x 10(-43)), and we found two exclusive haplotypes (196F/532S, frequency 0.74; 196S/532L, frequency 0.26) of the IRAK1 gene among our test subjects. In both populations, a significant association with decreased radial bone mineral density was identified with haplotype 196F/532S (in cohort A: r = 0.21, p = 0.0017; in cohort T: r = 0.23, p = 0.011). Radial bone mineral density was lowest among 196F/532S homozygotes, highest among 196S/532L homozygotes, and intermediate among heterozygotes. Accelerated bone loss also correlated with the 196F/532S haplotype in a 5-year follow-up. These results suggest that variation of IRAK1 may be an important determinant of postmenopausal osteoporosis, in part through the mechanism of accelerated postmenopausal bone loss.

Linkage disequilibrium and haplotype diversity in the genes of the renin-angiotensin system: findings from the family blood pressure program

Association studies of candidate genes with complex traits have generally used one or a few single nucleotide polymorphisms (SNPs), although variation in the extent of linkage disequilibrium (LD) within genes markedly influences the sensitivity and precision of association studies. The extent of LD and the underlying haplotype structure for most candidate genes are still unavailable. We sampled 193 blacks (African-Americans) and 160 whites (European-Americans) and estimated the intragenic LD and the haplotype structure in four genes of the renin-angiotensin system. We genotyped 25 SNPs, with all but one of the pairs spaced between 1 and 20 kb, thus providing resolution at small scale. The pattern of LD within a gene was very heterogeneous. Using a robust method to define haplotype blocks, blocks of limited haplotype diversity were identified at each locus; between these blocks, LD was lost owing to the history of recombination events. As anticipated, there was less LD among blacks, the number of haplotypes was substantially larger, and shorter haplotype segments were found, compared with whites. These findings have implications for candidate-gene association studies and indicate that variation between populations of European and African origin in haplotype diversity is characteristic of most genes.

Linkage disequilibrium mapping identifies a 390 kb region associated with CYP2D6 poor drug metabolising activity

The cytochrome p450 enzyme, CYP2D6, metabolises approximately 20% of marketed drugs. CYP2D6 multiple variants are associated with altered enzyme activities. Genotyping 1018 Caucasians for CYP2D6 polymorphisms (G1846A, delT1707, delA2549 and A2935C), known to result in the recessive CYP2D6 poor drug metaboliser (PM) phenotype, identified 41 individuals with predicted PM phenotype. These 41 individuals were classified as 'cases'. Single nucleotide polymorphisms (SNPs) mapping within an 880 kb region flanking CYP2D6, were identified to evaluate potential association between genetic variation and the CYP2D6 PM phenotype. The 41 PM cases and 977 controls were genotyped and analysed for 27 SNPs. Associations were observed across a 390 kb region between 14 SNPs and the PM phenotype (P values from 6.20 x 10(-4) to 4.54 x 10(-35)). Haplotype analysis revealed more significant levels of association (P = 3.54 x 10(-56)). Strong (D' > 0.7) linkage disequilibrium (LD) between SNPs was observed across the same 390 kb region associated with the CYP2D6 phenotype. The observed phenotype:genotype association reached genome-wide levels of significance, and supports the strategy for potential application of LD mapping and whole genome association scans to pharmacogenetic studies.

Structure of linkage disequilibrium in plants

Future advances in plant genomics will make it possible to scan a genome for polymorphisms associated with qualitative and quantitative traits. Before this potential can be realized, we must understand the nature of linkage disequilibrium (LD) within a genome. LD, the nonrandom association of alleles at different loci, plays an integral role in association mapping, and determines the resolution of an association study. Recently, association mapping has been exploited to dissect quantitative trait loci (QTL). With the exception of maize and Arabidopsis, little research has been conducted on LD in plants. The mating system of the species (selfing versus outcrossing), and phenomena such as population structure and recombination hot spots, can strongly influence patterns of LD. The basic patterns of LD in plants will be better understood as more species are analyzed.

Association study of polymorphisms in the GluR6 kainate receptor gene (GRIK2) with schizophrenia

The glutamatergic dysfunction hypothesis of schizophrenia suggests genes involved in glutamatergic transmission as candidates for schizophrenia-susceptibility genes. The GluR6 kainate receptor gene GRIK2 is located on chromosome 6q16.3-q21, a schizophrenia susceptibility region, as suggested by multiple linkage studies. We examined 15 SNPs evenly distributed in the entire GRIK2 region (>700 kb) in Japanese patients with schizophrenia (n=100) and controls (n=100). Neither genotype nor allele frequency showed a significant association with the disorder. We constructed 2-SNP haplotypes from the 15 SNPs. Although we observed three long linkage disequilibrium blocks (>150 kb) within the GRIK2 region, none of the pairwise haplotypes showed a significant association with the disorder. Therefore, we conclude that GRIK2 does not play a major role in the pathogenesis of schizophrenia in the Japanese population.

Haplotype and linkage disequilibrium architecture for human cancer-associated genes

To facilitate association-based linkage studies we have studied the linkage disequilibrium (LD) and haplotype architecture around five genes of interest for cancer risk: ATM, BRCA1, BRCA2, RAD51, and TP53. Single nucleotide polymorphisms (SNPs) were identified and used to construct haplotypes that span 93-200 kb per locus with an average SNP density of 12 kb. These markers were genotyped in four ethnically defined populations that contained 48 each of African Americans, Asian Americans, Hispanic Americans, and European Americans. Haplotypes were inferred using an expectation maximization (EM) algorithm, and the data were analyzed using D', R(2), Fisher's exact P-values, and the four-gamete test for recombination. LD levels varied widely between loci from continuously high LD across 200 kb to a virtual absence of LD across a similar length of genome. LD structure also varied at each gene and between populations studied. This variation indicates that the success of linkage-based studies will require a precise description of LD at each locus and in each population to be studied. One striking consistency between genes was that at each locus a modest number of haplotypes present in each population accounted for a high fraction of the total number of chromosomes. We conclude that each locus has its own genomic profile with regard to LD, and despite this there is the widespread trend of relatively low haplotype diversity. As a result, a low marker density should be adequate to identify haplotypes that represent the common variation at a locus, thereby decreasing costs and increasing efficacy of association studies.

Haplotype block structure and its applications to association studies: power and study designs

Recent studies have shown that the human genome has a haplotype block structure, such that it can be divided into discrete blocks of limited haplotype diversity. In each block, a small fraction of single-nucleotide polymorphisms (SNPs), referred to as "tag SNPs," can be used to distinguish a large fraction of the haplotypes. These tag SNPs can potentially be extremely useful for association studies, in that it may not be necessary to genotype all SNPs; however, this depends on how much power is lost. Here we develop a simulation study to quantitatively assess the power loss for a variety of study designs, including case-control designs and case-parental control designs. First, a number of data sets containing case-parental or case-control samples are generated on the basis of a disease model. Second, a small fraction of case and control individuals in each data set are genotyped at all the loci, and a dynamic programming algorithm is used to determine the haplotype blocks and the tag SNPs based on the genotypes of the sampled individuals. Third, the statistical power of tests was evaluated on the basis of three kinds of data: (1) all of the SNPs and the corresponding haplotypes, (2) the tag SNPs and the corresponding haplotypes, and (3) the same number of randomly chosen SNPs as the number of tag SNPs and the corresponding haplotypes. We study the power of different association tests with a variety of disease models and block-partitioning criteria. Our study indicates that the genotyping efforts can be significantly reduced by the tag SNPs, without much loss of power. Depending on the specific haplotype block-partitioning algorithm and the disease model, when the identified tag SNPs are only 25% of all the SNPs, the power is reduced by only 4%, on average, compared with a power loss of approximately 12% when the same number of randomly chosen SNPs is used in a two-locus haplotype analysis. When the identified tag SNPs are approximately 14% of all the SNPs, the power is reduced by approximately 9%, compared with a power loss of approximately 21% when the same number of randomly chosen SNPs is used in a two-locus haplotype analysis. Our study also indicates that haplotype-based analysis can be much more powerful than marker-by-marker analysis.

Nucleotide variability at G6pd and the signature of malarial selection in humans

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common enzymopathy in humans. Deficiency alleles for this X-linked disorder are geographically correlated with historical patterns of malaria, and the most common deficiency allele in Africa (G6PD A-) has been shown to confer some resistance to malaria in both hemizygous males and heterozygous females. We studied DNA sequence variation in 5.1 kb of G6pd from 47 individuals representing a worldwide sample to examine the impact of selection on patterns of human nucleotide diversity and to infer the evolutionary history of the G6PD A- allele. We also sequenced 3.7 kb of a neighboring locus, L1cam, from the same set of individuals to study the effect of selection on patterns of linkage disequilibrium. Despite strong clinical evidence for malarial selection maintaining G6PD deficiency alleles in human populations, the overall level of nucleotide heterozygosity at G6pd is typical of other genes on the X chromosome. However, the signature of selection is evident in the absence of genetic variation among A- alleles from different parts of Africa and in the unusually high levels of linkage disequilibrium over a considerable distance of the X chromosome. In spite of a long-term association between Plasmodium falciparum and the ancestors of modern humans, patterns of nucleotide variability and linkage disequilibrium suggest that the A- allele arose in Africa only within the last 10,000 years and spread due to selection.

SNP genotyping on a genome-wide amplified DOP-PCR template

With the increasing demand for higher throughput single nucleotide polymorphism (SNP) genotyping, the quantity of genomic DNA often falls short of the number of assays required. We investigated the use of degenerate oligonucleotide primed polymerase chain reaction (DOP-PCR) to generate a template for our SNP genotyping methodology of fluorescence polarization template-directed dye-terminator incorporation detection. DOP-PCR employs a degenerate primer (5'-CCGACTCGAGNNNNNNATGTGG-3') to produce non-specific uniform amplification of DNA. This approach has been successfully applied to microsatellite genotyping. We compared genotyping of DOP-PCR-amplified genomic DNA to genomic DNA as a template. Results were analyzed with respect to feasibility, allele loss of alleles, genotyping accuracy and storage conditions in a high-throughput genotyping environment. DOP-PCR yielded overall satisfactory results, with a certain loss in accuracy and quality of the genotype assignments. Accuracy and quality of genotypes generated from the DOP-PCR template also depended on storage conditions. Adding carrier DNA to a final concentration of 10 ng/microl improved results. In conclusion, we have successfully used DOP-PCR to amplify our genomic DNA collection for subsequent SNP genotyping as a standard process.

The impact of genotyping error on haplotype reconstruction and frequency estimation

The choice of genotyping families vs unrelated individuals is a critical factor in any large-scale linkage disequilibrium (LD) study. The use of unrelated individuals for such studies is promising, but in contrast to family designs, unrelated samples do not facilitate detection of genotyping errors, which have been shown to be of great importance for LD and linkage studies and may be even more important in genotyping collaborations across laboratories. Here we employ some of the most commonly-used analysis methods to examine the relative accuracy of haplotype estimation using families vs unrelateds in the presence of genotyping error. The results suggest that even slight amounts of genotyping error can significantly decrease haplotype frequency and reconstruction accuracy, that the ability to detect such errors in large families is essential when the number/complexity of haplotypes is high (low LD/common alleles). In contrast, in situations of low haplotype complexity (high LD and/or many rare alleles) unrelated individuals offer such a high degree of accuracy that there is little reason for less efficient family designs. Moreover, parent-child trios, which comprise the most popular family design and the most efficient in terms of the number of founder chromosomes per genotype but which contain little information for error detection, offer little or no gain over unrelated samples in nearly all cases, and thus do not seem a useful sampling compromise between unrelated individuals and large families. The implications of these results are discussed in the context of large-scale LD mapping projects such as the proposed genome-wide haplotype map.

Human genome sequence variation and the influence of gene history, mutation and recombination

Variation in the human genome sequence is key to understanding susceptibility to disease in modern populations and the history of ancestral populations. Unlocking this information requires knowledge of the patterns and underlying causes of human sequence diversity. By applying a new population-genetic framework to two genome-wide polymorphism surveys, we find that the human genome contains sizeable regions (stretching over tens of thousands of base pairs) that have intrinsically high and low rates of sequence variation. We show that the primary determinant of these patterns is shared genealogical history. Only a fraction of the variation (at most 25%) is due to the local mutation rate. By measuring the average distance over which genealogical histories are typically preserved, these data provide the first genome-wide estimate of the average extent of correlation among variants (linkage disequilibrium). The results are best explained by extreme variability in the recombination rate at a fine scale, and provide the first empirical evidence that such recombination 'hot spots' are a general feature of the human genome and have a principal role in shaping genetic variation in the human population.

Analytic strategies for stroke genetics

Strokes occur in first-degree relatives of stroke patients at rates approximating 40%; estimated population prevalences range from 3%-12%. The recent discovery of the Strk1 gene in Icelandic families, associated with a 7-fold increased stroke risk in carriers, strongly implicates a genetic contribution to stroke independent of other genetically based risk factors such as hypertension or diabetes. In this article we review the evidence for genetic contributions to stroke. We describe roles for genetic and molecular epidemiology in evaluating this complex disease. We delineate study design issues and analytic strategies to determine genetic association.

Distribution of nonrandom associations between pairs of protein loci along the third chromosome of Drosophila melanogaster

The within-chromosome distribution of gametic disequilibrium (GD) between protein loci, and the underlying evolutionary factors of this distribution, are still largely unknown. Here, we report a detailed study of GD between a large number of protein loci (15) spanning 87% of the total length of the third chromosome of Drosophila melanogaster in a large sample of haplotypes (600) drawn from a single natural population. We used a sign-based GD estimation method recently developed for multiallelic systems, which considerably increases both the statistical power and the accuracy of estimation of the intensity of GD. We found that strong GD between pairs of protein loci was widespread throughout the chromosome. In total, 22% of both the pairs of alleles and pairs of loci were in significant GD, with mean intensities (as measured by D' coefficients) of 0.43 and 0.31, respectively. In addition, strong GD often occurs between loci that are far apart. By way of illustration, 32% of the allele pairs in significant GD occurred within pairs of loci separated by effective frequencies of recombination (EFRs) of 15-20 cM, the mean D' value being 0.49. These observations are in sharp contrast with previous studies showing that GD between protein loci is rarely found in natural populations of outcrossing species, even between very closely linked loci. Interestingly, we found that most instances of significant interallelic GD (68%) involved functionally related protein loci. Specifically, GD was markedly more frequent between protein loci related by the functions of hormonal control, molybdenum control, antioxidant defense system, and reproduction than between loci without known functional relationship, which is indicative of epistatic selection. Furthermore, long-distance GD between functionally related loci (mean EFR 9 cM) suggests that epistatic interactions must be very strong along the chromosome. This evidence is hardly compatible with the neutral theory and has far-reaching implications for understanding the multilocus architecture of the functional genome. Our findings also suggest that GD may be a useful tool for discovering networks of functionally interacting proteins.

A first-generation linkage disequilibrium map of human chromosome 22

DNA sequence variants in specific genes or regions of the human genome are responsible for a variety of phenotypes such as disease risk or variable drug response. These variants can be investigated directly, or through their non-random associations with neighbouring markers (called linkage disequilibrium (LD)). Here we report measurement of LD along the complete sequence of human chromosome 22. Duplicate genotyping and analysis of 1,504 markers in Centre d'Etude du Polymorphisme Humain (CEPH) reference families at a median spacing of 15 kilobases (kb) reveals a highly variable pattern of LD along the chromosome, in which extensive regions of nearly complete LD up to 804 kb in length are interspersed with regions of little or no detectable LD. The LD patterns are replicated in a panel of unrelated UK Caucasians. There is a strong correlation between high LD and low recombination frequency in the extant genetic map, suggesting that historical and contemporary recombination rates are similar. This study demonstrates the feasibility of developing genome-wide maps of LD.

Single nucleotide polymorphism seeking long term association with complex disease

Successful investigation of common diseases requires advances in our understanding of the organization of the genome. Linkage disequilibrium provides a theoretical basis for performing candidate gene or whole-genome association studies to analyze complex disease. However, to constructively interrogate SNPs for these studies, technologies with sufficient throughput and sensitivity are required. A plethora of suitable and reliable methods have been developed, each of which has its own unique advantage. The characteristics of the most promising genotyping and polymorphism scanning technologies are presented. These technologies are examined both in the context of complex disease investigation and in their capacity to face the unique physical and molecular challenges (allele amplification, loss of heterozygosity and stromal contamination) of solid tumor research.