Localization of psoriasis-susceptibility locus PSORS1 to a 60-kb interval telomeric to HLA-C

Recent genome scans have established the presence of a major psoriasis-susceptibility locus in the human leukocyte antigen (HLA) complex on chromosome 6p21.3. To narrow the interval for candidate gene testing, we performed a linkage-disequilibrium analysis of 339 families, with the use of 62 physically mapped microsatellite markers spanning the major histocompatibility complex (MHC). As detected by use of the transmission/disequilibrium test (TDT), individual markers yielded significant linkage disequilibrium across most of the MHC. However, the strongest evidence for marker-trait disequilibrium was found in an approximately 300-kb region extending from the MICA gene to the corneodesmosin gene. Maximum-likelihood haplotypes were constructed across the entire MHC in the original sample and across a 1.2-Mb region of the central MHC in an expanded sample containing 139 additional families. Short (two- to five-marker) haplotypes were subjected to the TDT using a "moving-window" strategy that reduced the variability of TDT P values relative to the single-locus results. Furthermore, the expanded sample yielded a sharp peak of evidence for linkage disequilibrium that spanned approximately 170 kb and that was centered 100 kb telomeric to HLA-C. The 1.2-Mb interval was further dissected by means of recombinant ancestral haplotype analysis. This analysis identified risk haplotype 1 (RH1), which is a 60-kb fragment of ancestral haplotype 57.1, on all identifiable HLA risk haplotypes. One of these haplotypes exhibits significant linkage disequilibrium with psoriasis but does not carry Cw6, which is the HLA allele most strongly associated with the disease. These results demonstrate that RH1 is highly likely to carry the disease allele at PSORS1, and they exclude HLA-C and corneodesmosin with a high degree of confidence.     

吉林人群强直性脊柱炎6号染色体短臂上的HLA区域遗传易感基因定位研究

为了寻找中国人群中与强直性脊柱炎相关的新的易感基因及其所在位置, 在与强直性脊柱炎强连锁的6 号染色体短臂上的HLA基因区域内选取11个SNPs多态位点, 通过对中国吉林地区79名AS患者和132名正常对照者进行case-control分析, 发现TNF-a -850处TT突变基因型在AS组中的分布高于正常对照组（P=0.027）, 突变型T等位基因在AS组和正常对照组中的分布差异更为显著（P=0.002）。通过多位点之间的连锁不平衡分析发现, LTA基因、TNF-a基因、LST1基因和NCR3基因中的 5个SNPs多态位点之间存在连锁不平衡, 范围是15 kb, 在这5个SNPs多态位点组成的单体型中, TCTTC单体型在AS组和正常对照组中的分布有显著差异（c2=7.406, P=0.0065）,并且该单体型中含有具有统计学意义的TNF-a –850的突变型等位基因T。提示在LTA、TNF-a、NCR3和LST1 这4个基因构成的15 kb范围内可能存在增加AS患病易感性的位点, 可能是TNF-a –850 C→T突变, 也可能是在TNF-a –850附近的其他位点。     

Linkage disequilibrium in the North American Holstein population

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium ( r ²   >   0.8) in genomic regions of ∼50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26–113 kb by various definitions, which was larger than that observed in humans (∼10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.     

Linkage disequilibria between pairs of loci within a highly polymorphic region of chromosome 2Q.

A recently described region on chromosome 2q contains seven restriction fragment length polymorphisms (RFLPs) revealed by single-copy probes isolated from a 20-kilobase (kb) segment of a single cosmid insert. Analysis of six of these loci demonstrates modest amounts of linkage disequilibrium. This reflects the presence of a substantial number of different haplotypes in this chromosome region and indicates that the region could be used as one highly polymorphic locus. No consistent relationship is found between the amount of linkage disequilibrium and the physical distance between pairs of loci. For seven of the 10 pairs of diallelic loci studied, the observed disequilibrium can be attributed primarily to the absence of the minor haplotype from the population. These results suggest that, for small regions of the genome, factors such as mutation, genetic drift, and population admixture may have effects that outweight those of recombination. In addition, results are reviewed which show that estimates of linkage disequilibrium coefficients for tightly linked loci are very imprecise. Thus, the inference of gene order from linkage disequilibrium values must be regarded with caution.     

Power and sample size calculations for genetic case/control studies using gene-centric SNP maps: application to human chromosomes 6, 21, and 22 in three populations

Power and sample size calculations are critical parts of any research design for genetic association. We present a method that utilizes haplotype frequency information and average marker-marker linkage disequilibrium on SNPs typed in and around all genes on a chromosome. The test statistic used is the classic likelihood ratio test applied to haplotypes in case/control populations. Haplotype frequencies are computed through specification of genetic model parameters. Power is determined by computation of the test's non-centrality parameter. Power per gene is computed as a weighted average of the power assuming each haplotype is associated with the trait. We apply our method to genotype data from dense SNP maps across three entire chromosomes (6, 21, and 22) for three different human populations (African-American, Caucasian, Chinese), three different models of disease (additive, dominant, and multiplicative) and two trait allele frequencies (rare, common). We perform a regression analysis using these factors, average marker-marker disequilibrium, and the haplotype diversity across the gene region to determine which factors most significantly affect average power for a gene in our data. Also, as a 'proof of principle' calculation, we perform power and sample size calculations for all genes within 100 kb of the PSORS1 locus (chromosome 6) for a previously published association study of psoriasis. Results of our regression analysis indicate that four highly significant factors that determine average power to detect association are: disease model, average marker-marker disequilibrium, haplotype diversity, and the trait allele frequency. These findings may have important implications for the design of well-powered candidate gene association studies. Our power and sample size calculations for the PSORS1 gene appear consistent with published findings, namely that there is substantial power (>0.99) for most genes within 100 kb of the PSORS1 locus at the 0.01 significance level.     

Statistical estimation and pedigree analysis of CCR2-CCR5 haplotypes

As more SNP marker data becomes available, researchers have used haplotypes of markers, rather than individual polymorphisms, for association analysis of candidate genes. In order to perform haplotype analysis in a population-based case-control study, haplotypes must be determined by estimation in the absence of family information or laboratory methods for establishing phase. Here, we test the accuracy of the Expectation-Maximization (EM) algorithm for estimating haplotype state and frequency in the CCR2-CCR5 gene region by comparison with haplotype state and frequency determined by pedigree analysis. To do this, we have characterized haplotypes comprising alleles at seven biallelic loci in the CCR2-CCR5 chemokine receptor gene region, a span of 20 kb on chromosome 3p21. Three-generation CEPH families (n=40), totaling 489 individuals, were genotyped by the 5'nuclease assay (TaqMan). Haplotype states and frequencies were compared in 103 grandparents who were assumed to have mated at random. Both pedigree analysis and the EM algorithm yielded the same small number of haplotypes for which linkage disequilibrium was nearly maximal. The haplotype frequencies generated by the two methods were nearly identical. These results suggest that the EM algorithm estimation of haplotype states, frequency, and linkage disequilibrium analysis will be an effective strategy in the CCR2-CCR5 gene region. For genetic epidemiology studies, CCR2-CCR5 allele and haplotype frequencies were determined in African-American (n=30), Hispanic (n=24) and European-American (n=34) populations.     

Linkage disequilibrium mapping of Arabidopsis CRY2 flowering time alleles

The selfing plant Arabidopsis thaliana has been proposed to be well suited for linkage disequilibrium (LD) mapping as a means of identifying genes underlying natural trait variation. Here we apply LD mapping to examine haplotype variation in the genomic region of the photoperiod receptor CRYPTOCHROME2 and associated flowering time variation. CRY2 DNA sequences reveal strong LD and the existence of two highly differentiated haplogroups (A and B) across the gene; in addition, a haplotype possessing a radical glutamine-to-serine replacement (AS) occurs within the more common haplogroup. Growth chamber and field experiments using an unstratified population of 95 ecotypes indicate that under short-day photoperiod, the AS and B haplogroups are both highly significantly associated with early flowering. Data from six genes flanking CRY2 indicate that these haplogroups are limited to an approximately 65-kb genomic region around CRY2. Whereas the B haplogroup cannot be delimited to <16 kb around CRY2, the AS haplogroup is characterized almost exclusively by the nucleotide polymorphisms directly associated with the serine replacement in CRY2; this finding strongly suggests that the serine substitution is directly responsible for the AS early flowering phenotype. This study demonstrates the utility of LD mapping for elucidating the genetic basis of natural, ecologically relevant variation in Arabidopsis.     

Haplotype Analysis of Hemochromatosis: Evaluation of Different Linkage-Disequilibrium Approaches and Evolution of Disease Chromosomes

We applied several types of linkage-disequilibrium calculations to analyze the hereditary hemochromatosis (hh) locus. Twenty-four polymorphic markers in the major histocompatibility complex (MHC) class I region were used to haplotype hh and normal chromosomes. A total of 169 hh and 161 normal chromosomes were analyzed. Disequilibrium values were found to be high over an unusually large region beginning 150 kb centro-meric of HLA-A and extending nearly 5 Mb telomeric of it. Recombination in this region was −28% of the expected value. This low level of recombination contributes to the unusually broad region of linkage disequilibrium found with hh. The strongest disequilibrium was found at locus HLA-H (d = .84) and at locus D6S2239 (d = .85), a marker −10 kb telomeric to HLA-H. All disequilibrium methods employed in this study found peak disequilibrium at HLA-H or D6S2239. The cys282tyr mutation in HLA-H, a candidate gene for hh, was found in 85% of disease chromosomes. A haplotype phylogeny for hh chromosomes was constructed and suggests that the mutation associated with the most common haplotype occurred relatively recently. The age of the hh mutation was estimated to be −60-70 generations. Disequilibrium was maintained over a greater distance for hh-carrying chromosomes, consistent with a recent mutation for hh. Our data provide a reasonable explanation for previous difficulties in localizing the hh locus and provide an evolutionary history for disease chromosomes.     

Disequilibrium Pattern Analysis. I. Theory

We have developed a method, disequilibrium pattern analysis, for examining the disequilibrium distribution of the entire array of two locus multiallelic haplotypes in a population. It is shown that a selected haplotype will produce a distinct pattern of linkage disequilibrium values for all generations while the selection is acting. This pattern will also presumably be maintained for many generations after the selection event, until the disequilibrium pattern is eventually broken down by genetic drift and recombination. Related haplotypes, sharing an allele with a selected haplotype, assume a value of linkage disequilibrium proportional to the frequency of the unshared allele and have a single negative value of the normalized linkage disequilibrium. The analysis assumes zero linkage disequilibrium for all allelic combinations initially. The same basic results continue to apply if the selection involves a new mutant, the occurrence of which creates linkage disequilibrium for some haplotypes. The disequilibrium pattern predicted under selection is robust with respect to the influence of migration and random genetic drift. This method is applicable to population data having linked polymorphic loci including that determined from protein or DNA sequencing.     

Haplotype structure and phenotypic associations in the chromosomal regions surrounding two Arabidopsis thaliana flowering time loci

The feasibility of using linkage disequilbrium (LD) to fine-map loci underlying natural variation in Arabidopsis thaliana was investigated by looking for associations between flowering time and marker polymorphism in the genomic regions containing two candidate genes, FRI and FLC, both of which are known to contribute to natural variation in flowering. A sample of 196 accessions was used, and polymorphism was assessed by sequencing a total of 17 roughly 500-bp fragments. Using a novel Bayesian algorithm based on haplotype similarity, we demonstrate that LD could have been used to fine-map the FRI gene to a roughly 30-kb region and to identify two common loss-of-function alleles. Interestingly, because of genetic heterogeneity, simple single-marker associations would not have been able to map FRI with nearly the same precision. No clear evidence for previously unknown alleles at either locus was found, but the effect of population structure in causing false positives was evident.     

Molecular,genetic and evolutionary analysis of a paracentric inversion in Arabidopsis thaliana

Chromosomal inversions can provide windows onto the cytogenetic, molecular, evolutionary and demographic histories of a species. Here we investigate a paracentric 1.17‐Mb inversion on chromosome 4 of Arabidopsis thaliana with nucleotide precision of its borders. The inversion is created by Vandal transposon activity, splitting an F‐box and relocating a pericentric heterochromatin segment in juxtaposition with euchromatin without affecting the epigenetic landscape. Examination of the RegMap panel and the 1001 Arabidopsis genomes revealed more than 170 inversion accessions in Europe and North America. The SNP patterns revealed historical recombinations from which we infer diverse haplotype patterns, ancient introgression events and phylogenetic relationships. We find a robust association between the inversion and fecundity under drought. We also find linkage disequilibrium between the inverted region and the early flowering Col‐FRIGIDA allele. Finally, SNP analysis elucidates the origin of the inversion to South‐Eastern Europe approximately 5000 years ago and the FRI‐Col allele to North‐West Europe, and reveals the spreading of a single haplotype to North America during the 17th to 19th century. The ‘American haplotype’ was identified from several European localities, potentially due to return migration.     

Modeling haplotype block variation using Markov chains

Models of background variation in genomic regions form the basis of linkage disequilibrium mapping methods. In this work we analyze a background model that groups SNPs into haplotype blocks and represents the dependencies between blocks by a Markov chain. We develop an error measure to compare the performance of this model against the common model that assumes that blocks are independent. By examining data from the International Haplotype Mapping project, we show how the Markov model over haplotype blocks is most accurate when representing blocks in strong linkage disequilibrium. This contrasts with the independent model, which is rendered less accurate by linkage disequilibrium. We provide a theoretical explanation for this surprising property of the Markov model and relate its behavior to allele diversity.     

Extent and distribution of linkage disequilibrium in three genomic regions

The positional cloning of genes underlying common complex diseases relies on the identification of linkage disequilibrium (LD) between genetic markers and disease. We have examined 127 polymorphisms in three genomic regions in a sample of 575 chromosomes from unrelated individuals of British ancestry. To establish phase, 800 individuals were genotyped in 160 families. The fine structure of LD was found to be highly irregular. Forty-five percent of the variation in disequilibrium measures could be explained by physical distance. Additional factors, such as allele frequency, type of polymorphism, and genomic location, explained <5% of the variation. Nevertheless, disequilibrium was occasionally detectable at 500 kb and was present for over one-half of marker pairs separated by <50 kb. Although these findings are encouraging for the prospects of a genomewide LD map, they suggest caution in interpreting localization due to allelic association.     

A primary assembly of a bovine haplotype block map based on a 15,036-single-nucleotide polymorphism panel genotyped in holstein-friesian cattle

Analysis of data on 1000 Holstein-Friesian bulls genotyped for 15,036 single-nucleotide polymorphisms (SNPs) has enabled genomewide identification of haplotype blocks and tag SNPs. A final subset of 9195 SNPs in Hardy-Weinberg equilibrium and mapped on autosomes on the bovine sequence assembly (release Btau 3.1) was used in this study. The average intermarker spacing was 251.8 kb. The average minor allele frequency (MAF) was 0.29 (0.05-0.5). Following recent precedents in human HapMap studies, a haplotype block was defined where 95% of combinations of SNPs within a region are in very high linkage disequilibrium. A total of 727 haplotype blocks consisting of > or =3 SNPs were identified. The average block length was 69.7 +/- 7.7 kb, which is approximately 5-10 times larger than in humans. These blocks comprised a total of 2964 SNPs and covered 50,638 kb of the sequence map, which constitutes 2.18% of the length of all autosomes. A set of tag SNPs, which will be useful for further fine-mapping studies, has been identified. Overall, the results suggest that as many as 75,000-100,000 tag SNPs would be needed to track all important haplotype blocks in the bovine genome. This would require approximately 250,000 SNPs in the discovery phase.     

Recombination and gene conversion in a 170-kb genomic region of Arabidopsis thaliana

Arabidopsis thaliana is a highly selfing plant that nevertheless appears to undergo substantial recombination. To reconcile its selfing habit with the observations of recombination, we have sampled the genetic diversity of A. thaliana at 14 loci of approximately 500 bp each, spread across 170 kb of genomic sequence centered on a QTL for resistance to herbivory. A total of 170 of the 6321 nucleotides surveyed were polymorphic, with 169 being biallelic. The mean silent genetic diversity (pi(s)) varied between 0.001 and 0.03. Pairwise linkage disequilibria between the polymorphisms were negatively correlated with distance, although this effect vanished when only pairs of polymorphisms with four haplotypes were included in the analysis. The absence of a consistent negative correlation between distance and linkage disequilibrium indicated that gene conversion might have played an important role in distributing genetic diversity throughout the region. We tested this by coalescent simulations and estimate that up to 90% of recombination is due to gene conversion.     

Linkage disequilibrium mapping in domestic dog breeds narrows the progressive rod-cone degeneration interval and identifies ancestral disease-transmitting chromosome

Canine progressive rod-cone degeneration (prcd) is a retinal disease previously mapped to a broad, gene-rich centromeric region of canine chromosome 9. As allelic disorders are present in multiple breeds, we used linkage disequilibrium (LD) to narrow the approximately 6.4-Mb interval candidate region. Multiple dog breeds, each representing genetically isolated populations, were typed for SNPs and other polymorphisms identified from BACs. The candidate region was initially localized to a 1.5-Mb zero recombination interval between growth factor receptor-bound protein 2 (GRB2) and SEC14-like 1 (SEC14L). A fine-scale haplotype of the region was developed, which reduced the LD interval to 106 kb and identified a conserved haplotype of 98 polymorphisms present in all prcd-affected chromosomes from 14 different dog breeds. The findings strongly suggest that a common ancestor transmitted the prcd disease allele to many of the modern dog breeds and demonstrate the power of the LD approach in the canine model.     

Confirmation of linkage disequilibrium between haplotype B (XV-2c,allele 1; KM-19, allele 2) and cystic fibrosis allele in the French population

Summary In 237 French families with cystic fibrosis (CF) restricted fragment length polymorphisms (RFLPs) were detected by two DNa probes, XV-2c and KM-19, which are tightly linked to the CF allele. As in other European populations linkage disequilibrium is found between the haplotype B (XV-2c, allele 1: KM-19, allele 2) and the CF allele. Linkage disequilibrium alters the probability that a person bearing a given haplotype is a carrier.     

Positional cloning of disease genes: advantages of genetic isolates.

Genetic isolates with a history of a small founder population, long-lasting isolation and population bottlenecks represent exceptional resources in the identification of disease genes. Specific rare, monogenic diseases become enriched, and families with multiple affected individuals occur frequently enough to be used in linkage analyses for locus identification. Further, the vast majority of cases are caused by the same mutation, and disease alleles reveal linkage disequilibrium (LD) with markers over significant genetic intervals; this facilitates disease locus identification by similarity search for a shared genotype or haplotype in small study samples consisting of few affected individuals. LD observed in disease alleles adds power to linkage analyses and helps to define the exact location of disease loci on the genetic map. Typically, based on the linkage disequilibrium and the ancient haplotype, the critical DNA region can be defined from the original 1- to 2-cM resolution obtained in linkage analysis to 50-200 kb, greatly facilitating the targeting of physical cloning and sequencing efforts. These advantages have been well demonstrated in the positional cloning of several rare monogenic diseases enriched in population isolates like the example of Finland used here. How useful genetic isolates will prove to be in the identification of complex disease genes is dependent on the genealogical history of the isolate, including the size of the founding population and the expansion rate during the history of the population.     

An analysis of linkage disequilibrium in the interleukin-1 gene cluster, using a novel grouping method for multiallelic markers.

In population- and family-based association studies, it is useful to have some knowledge of the patterns of linkage disequilibrium that exist between markers in candidate regions. When such studies are carried out with multiallelic markers, it is often convenient to group the alleles into a biallelic system, for analysis. In this study, we specifically examined the interleukin-1 (IL-1) gene cluster on chromosome 2, a region containing candidates for many inflammatory and autoimmune disorders. Data were collected on eight markers, four of which were multiallelic. Using these data, we investigated the effect of three allele-grouping strategies, including a novel method, on the detection of linkage disequilibrium. The novel approach, termed the "delta method," measures the deviation from the expected haplotype frequencies under linkage equilibrium, for each allelic combination. This information is then used to group the alleles, in an attempt to avoid the grouping together of alleles at one locus that are in opposite disequilibrium with the same allele at the second locus. The estimate haplotype frequencies (EH) program was used to estimate haplotype frequencies and the disequilibrium measure. In our data it was found that the delta method compared well with the other two strategies. Using this method, we found that there was a reasonable correlation between disequilibrium and physical distance in the region (r=-.540, P=.001, one-tailed). We also identified a common, eight-locus haplotype of the IL-1 gene cluster.     

Genetic polymorphism at two linked loci,Sod and Est-6, in Drosophila melanogaster

We have examined the patterns of polymorphism at two linked loci, Sod and Est-6, separated by nearly 1000 kb on the left arm of chromosome 3 of Drosophila melanogaster. The evidence suggests that natural selection has been involved in shaping the polymorphisms. At the Sod locus, a fairly strong (s>0.01) selective sweep, started ≥2600 years ago, increased the frequency of a rare haplotype, F(A), to about 50% frequency in populations of Europe, Asia, and the Americas. More recently, an F(A) allele mutated to an S allele, which has increased to frequencies 5–15% in populations of Europe, Asia and North America. All S alleles are identical (or very nearly) in sequence and differ by one nucleotide substitution (which accounts for the F→S electrophoretic difference) from F(A) alleles. At the Est-6 locus, the evidence indicates both directional and balancing selection impacting separately the promoter and the coding regions of the gene, with linkage disequilibrium occurring within each region. Some linkage disequilibrium also exists between the two genes.