Haplotype analysis to determine the position of a mutation among closely linked DNA markers

The origin of Eastern European Jews revealed by autosomal, sex chromosomal and mtDNA polymorphisms

Background

This study aims to establish the likely origin of EEJ (Eastern European Jews) by genetic distance analysis of autosomal markers and haplogroups on the X and Y chromosomes and mtDNA.

Results

According to the autosomal polymorphisms the investigated Jewish populations do not share a common origin, and EEJ are closer to Italians in particular and to Europeans in general than to the other Jewish populations. The similarity of EEJ to Italians and Europeans is also supported by the X chromosomal haplogroups. In contrast according to the Y-chromosomal haplogroups EEJ are closest to the non-Jewish populations of the Eastern Mediterranean. MtDNA shows a mixed pattern, but overall EEJ are more distant from most populations and hold a marginal rather than a central position. The autosomal genetic distance matrix has a very high correlation (0.789) with geography, whereas the X-chromosomal, Y-chromosomal and mtDNA matrices have a lower correlation (0.540, 0.395 and 0.641 respectively).

Conclusions

The close genetic resemblance to Italians accords with the historical presumption that Ashkenazi Jews started their migrations across Europe in Italy and with historical evidence that conversion to Judaism was common in ancient Rome. The reasons for the discrepancy between the biparental markers and the uniparental markers are discussed.

Reviewers

This article was reviewed by Damian Labuda (nominated by Jerzy Jurka), Kateryna Makova and Qasim Ayub (nominated by Dan Graur).

A chronology of fine-scale gene mapping by linkage disequilibrium

The past decade produced several proposals for fine-scale gene mapping using linkage disequilibrium data. The suggested methods fall into two main groups, those that rely on pairwise statistics and those that rely on haplotypes. This paper reviews each strategy's development from a chronological perspective.

The trimmed-haplotype test for linkage disequilibrium

Single-marker linkage-disequilibrium (LD) methods cannot fully describe disequilibrium in an entire chromosomal region surrounding a disease allele. With the advent of myriad tightly linked microsatellite markers, we have an opportunity to extend LD analysis from single markers to multiple-marker haplotypes. Haplotype analysis has increased statistical power to disclose the presence of a disease locus in situations where it correctly reflects the historical process involved. For maximum efficiency, evidence of LD ought to come not just from a single haplotype, which may well be rare, but in addition from many similar haplotypes that could have descended from the same ancestral founder but have been trimmed in succeeding generations. We present such an analysis, called the "trimmed-haplotype method." We focus on chromosomal regions that are small enough that disequilibrium in significant portions of them may have been preserved in some pedigrees and yet that contain enough markers to minimize coincidental occurrence of the haplotype in the absence of a disease allele: perhaps regions 1-2 cM in length. In general, we could have no idea what haplotype an ancestral founder carried generations ago, nor do we usually have a precise chromosomal location for the disease-susceptibility locus. Therefore, we must search through all possible haplotypes surrounding multiple locations. Since such repeated testing obliterates the sampling distribution of the test, we employ bootstrap methods to calculate significance levels. Trimmed-haplotype analysis is performed on family data in which genotypes have been assembled into haplotypes. It can be applied either to conventional parent-affected-offspring triads or to multiplex pedigrees. We present a method for summarizing the LD evidence, in any pedigree, that can be employed in trimmed-haplotype analysis as well as in other methods.

Linkage disequilibrium mapping in isolated populations: the example of Finland revisited

Linkage disequilibrium analysis can provide high resolution in the mapping of disease genes because it incorporates information on recombinations that have occurred during the entire period from the mutational event to the present. A circumstance particularly favorable for high-resolution mapping is when a single founding mutation segregates in an isolated population. We review here the population structure of Finland in which a small founder population some 100 generations ago has expanded into 5.1 million people today. Among the 30-odd autosomal recessive disorders that are more prevalent in Finland than elsewhere, several appear to have segregated for this entire period in the "panmictic" southern Finnish population. Linkage disequilibrium analysis has allowed precise mapping and determination of genetic distances at the 0.1-cM level in several of these disorders. Estimates of genetic distance have proven accurate, but previous calculations of the confidence intervals were too small because sampling variation was ignored. In the north and east of Finland the population can be viewed as having been "founded" only after 1500. Disease mutations that have undergone such a founding bottleneck only 20 or so generations ago exhibit linkage disequilibrium and haplotype sharing over long genetic distances (5-15 cM). These features have been successfully exploited in the mapping and cloning of many genes. We review the statistical issues of fine mapping by linkage disequilibrium and suggest that improved methodologies may be necessary to map diseases of complex etiology that may have arisen from multiple founding mutations.

Linkage disequilibrium and gene mapping: an empirical least-squares approach

This paper proposes a novel approach for fine-scale mapping of disease genes that is based on the well-known linkage-disequilibrium parameter delta. Using a very simple, very general model, I show how delta can be interpreted in terms of identity-by-descent probabilities. The value of delta follows a piecewise curve along the chromosome, with the maximum occurring at the disease locus where the two pieces intersect. A semiparametric, multilocus approach is used to fit this nonlinear regression curve in order to estimate the gene location. Using the bootstrap to empirically estimate much of the probability model from the data avoids the need for many detailed population assumptions. One advantage of the approach is its use of the observed covariance structure of the data, which can be highly informative as to the gene location. I illustrate the method on the cystic fibrosis data of Kerem et al.

Molecular detection of new mutations, resolution of ambiguous results and complex genetic counseling issues in Huntington disease

Huntington disease (HD) is an autosomal dominant neurodegenerative disorder caused by expansion of a variable length (CAG)n repeat in the 5' coding region of a novel gene on chromosome 4p16.3. We provide comprehensive molecular analysis of a sporadic case of HD in which a paternally derived normal length allele expanded to an affected length allele. Linkage analysis and paternity testing confirm the paternal origin of the expansion and demonstrate that unequal crossing over during meiosis is an unlikely mechanism for de novo expansion in HD. This case identifies a complex genetic counseling issue for the families of sporadic cases since calculations of recurrence risk are not possible at this time. In addition, we describe utilization of a combination of polymerase chain reaction (PCR) based assays for examination of both the CAG repeat and an adjacent variable length CCG repeat in the huntingtin gene. The combination of these assays can increase the accuracy of molecular diagnosis for HD and may clarify any ambiguous results obtained during molecular testing of HD families.

Allelic association of microsatellites of 6p in Italian hemochromatosis patients

Hemochromatosis (HC) is an inherited disorder of iron metabolism and is frequently seen in Caucasians. The biochemical defect and the responsible gene are unknown, but the HC locus is closely linked to HLA-A on human chromosome 6 in the region 6p21.3. Although extensive studies have been performed in several populations, the precise location of the gene is still undefined. Linkage disequilibrium with HC has been detected for loci that are 3 cM apart: HLA class I and D6S105, which is located on the telomeric side of HLA-A. We have analyzed the inheritance of several multi-allele polymorphisms that map to 6p (D6S265, Y52, HLA-F, D6S306, D6S105, D6S464, D6S299) in 34 Italian HC families and in 17 unrelated patients. Significant association with HC was shown for alleles of multiple markers in the HLA-A region, for the distant marker D6S105, but not for the D6S299 marker at 4 cM from HLA-A on the telomeric side. HC status was unambiguously assigned to 70 affected and 63 unaffected chromosomes from family studies. Thirty five different haplotypes were found in 70 HC chromosomes when considering four markers most tightly associated with the disease. A predominant haplotype comprising alleles 1-3-1-8 (marker order D6S265, HLA-A, Y52, D6S105) accounted for 30% of the HC chromosomes and was absent in normals. A minority of other HC haplotypes could be related to the major haplotype by assuming single crossover events. Results of haplotype studies suggest a founder effect in the Italian population, as previously shown in Australian patients, and a possible common mutation shared with affected individuals of Celtic origin.

The diastrophic dysplasia gene encodes a novel sulfate transporter: positional cloning by fine-structure linkage disequilibrium mapping

Diastrophic dysplasia (DTD) is a well-characterized autosomal recessive osteochondrodysplasia with clinical features including dwarfism, spinal deformation, and specific joint abnormalities. The disease occurs in most populations, but is particularly prevalent in Finland owing to an apparent founder effect. DTD maps to distal chromosome 5q and, based on linkage disequilibrium studies in the Finnish population, we had previously predicted that the DTD gene should lie about 64 kb away from the CSF1R locus. Here, we report the positional cloning of the DTD gene by fine-structure linkage disequilibrium mapping. The gene lies in the predicted location, approximately 70 kb proximal to CSF1R, and encodes a novel sulfate transporter. Impaired function of its product is likely to lead to undersulfation of proteoglycans in cartilage matrix and thereby to cause the clinical phenotype of the disease. These results demonstrate the power of linkage disequilibrium mapping in isolated populations for positional cloning.

The origin of the major cystic fibrosis mutation (delta F508) in European populations

delta F508 is the most frequent cystic fibrosis (CF) mutation and accounts for approximately 70% of CF chromosomes worldwide. Three highly polymorphic microsatellite markers have been used to study the origin and evolution of delta F508 chromosomes in Europe. Haplotype data demonstrate that delta F508 occurred more than 52,000 years ago, in a population genetically distinct from any present European group, and spread throughout Europe in chronologically distinct expansions, which are responsible for the different frequencies of delta F508 in Europe.