Multilocus lod scores in large pedigrees: combination of exact and approximate calculations

To detect the positions of disease loci, lod scores are calculated at multiple chromosomal positions given trait and marker data on members of pedigrees. Exact lod score calculations are often impossible when the size of the pedigree and the number of markers are both large. In this case, a Markov Chain Monte Carlo (MCMC) approach provides an approximation. However, to provide accurate results, mixing performance is always a key issue in these MCMC methods. In this paper, we propose two methods to improve MCMC sampling and hence obtain more accurate lod score estimates in shorter computation time. The first improvement generalizes the block-Gibbs meiosis (M) sampler to multiple meiosis (MM) sampler in which multiple meioses are updated jointly, across all loci. The second one divides the computations on a large pedigree into several parts by conditioning on the haplotypes of some 'key' individuals. We perform exact calculations for the descendant parts where more data are often available, and combine this information with sampling of the hidden variables in the ancestral parts. Our approaches are expected to be most useful for data on a large pedigree with a lot of missing data.     

Detection of linkage for heterogeneous disorders by using multipoint linkage analysis.

We have compared the efficiency of the lod score test which assumes heterogeneity (lod2) to the standard lod score test which assumes homogeneity (lod1) when three-point linkage analysis is used in successive map intervals. If it is assumed that a gene located midway between two linked marker loci is responsible for a proportion of disease cases, then the lod1 test loses power relative to the lod2 test, as the proportion of linked families decreases, as the flanking markers are more closely linked, and as more map intervals are tested. Moreover, when multipoint analysis is used, linkage for a disease gene is more likely to be incorrectly excluded from a complete and dense linkage map if true genetic heterogeneity is ignored. We thus conclude that, in general, the lod2 linkage test is more efficient for detecting a true linkage when a complete genetic marker map is screened for a heterogeneous disorder.     

Closing in on a breast cancer gene on chromosome 17q.

Linkage of early-onset familial breast and ovarian cancer to 11 markers on chromosome 17q12-q21 defines an 8-cM region which is very likely to include the disease gene BRCA 1. The most closely linked marker is D17S579, a highly informative CA repeat polymorphism. D17S579 has no recombinants with inherited breast or ovarian cancer in 79 informative meioses in the seven families with early-onset disease (lod score 9.12 at zero recombination). There is no evidence for linkage heterogeneity in the families with early-onset disease. The proportion of older-onset breast cancer attributable to BRCA 1 is not yet determinable, because both inherited and sporadic cases occur in older-onset families.     

D19S51 is closely linked with and maps distal to the myotonic dystrophy locus on 19q.

Recent genetic linkage studies have mapped the myotonic dystrophy (DM) locus to 19q13.3. All closely linked DM markers identified to date have been located on the centromeric side of the disease locus, with a relatively large genetic interval (9 cM) observed between the nearest distal marker and DM. We show here that the recently described marker p134C is tightly linked to DM (peak lod score 35.8 at peak recombination fraction .006) and confirm the previous suggestion that the p134C locus, D19S51 maps distal to the disease locus. D19S51 and the closest proximal flanking loci, ERCC1 and D19S115 (pE0.8), define a small genetic interval of less than 2 cM that contains the DM locus.     

Linkage analysis of families with fragile-X mental retardation, using a novel RFLP marker (DXS 304).

A new polymorphic DNA marker U6.2, defining the locus DXS304, was recently isolated and mapped to the Xq27 region of the X chromosome. In the previous communication we describe a linkage study encompassing 16 fragile-X families and using U6.2 and five previously described polymorphic markers at Xq26-q28. One recombination event was observed between DXS304 and the fragile-X locus in 36 informative meioses. Combined with information from other reports, our results suggest the following order of the examined loci on Xq: cen-F9-DXS105-DXS98-FRAXA-DXS304-(DXS52-F8 -DXS15). The locus DXS304 is closely linked to FRAXA, giving a peak lod score of 5.86 at a corresponding recombination fraction of .00. On the basis of the present results, it is apparent that U6.2 is a useful probe for carrier and prenatal diagnosis in fragile-X families.     

Hereditary Hyperparathyroidism–Jaw Tumor Syndrome: The Endocrine Tumor Gene HRPT2 Maps to Chromosome 1q21-q31

The syndrome of hereditary hyperparathyroidism and jaw tumors (HPT-JT) is characterized by inheritance, in an autosomal dominant pattern, of recurrent parathyroid adenomas, fibro-osseous tumors of the mandible and/or maxilla, Wilms tumor, and parathyroid carcinoma. This syndrome is clinically and genetically distinct from other endocrine neoplasia syndromes and appears to result from mutation of an endocrine tumor gene designated “HRPT2.” We studied five HPT-JT families (59 persons, 20 affected); using PCR-based markers, we instituted a genomewide linkage search after excluding several candidate genes. Lod scores were calculated at various recombination fractions (θ), penetrance 90%. We mapped HRPT2 to the long arm of chromosome 1 (1q21-q31). The maximal lod score was 6.10 at θ = .0 with marker D1S212, or >10⁶ odds in favor of linkage. In six hereditary Wilms tumor families (96 persons, 29 affected), we found no linkage to 1q markers closely linked with HRPT2 (lod scores ?15.6 [D1S191] and ?17.8 [D1S196], θ = .001). Nine parathyroid adenomas and one Wilms tumor from nine members of three HPT-JT families were examined for loss of heterozygosity at linked loci. The parathyroid adenomas and Wilms tumor showed no loss of heterozygosity for these DNA markers. Our data establish that HRPT2, an endocrine tumor gene on the long arm of chromosome 1, is responsible for the HPT-JT syndrome but not for the classical hereditary Wilms tumor syndrome.     

Fragile X syndrome: diagnosis using highly polymorphic microsatellite markers.

We describe two highly polymorphic microsatellite AC repeat sequences, VK23AC and VK14AC, which are closely linked to the fragile X at Xq27.3. Both VK23AC (DXS297) and VK14AC (DXS292) are proximal to the fragile site. Two-point linkage analysis in 31 fragile X families gave (a) a recombination frequency of 1% (range 0.00%-4%) with a maximum lod score of 32.04 for DXS297 and (b) a recombination frequency of 7% (range of 3%-15%) with a maximum lod score of 12.87 for DXS292. Both of these polymorphisms are applicable to diagnosis by linkage in families with fragile X syndrome. A multipoint linkage map of genetic markers at Xq27.3 was constructed from genotyping these polymorphisms in the CEPH pedigrees. The DXS292 marker is in the DXS98-DXS297 interval and in 3 cM proximal to DXS297.     

New distal marker closely linked to the fragile X locus

Summary We have isolated II-10, a new X-chromosomal probe that identifies a highly informative two-allele TaqI restriction fragment length polymorphism at locus DXS466. Using somatic cell hybrids containing distinct portions of the long arm of the X chromosome, we could localize DXS466 between DXS296 and DXS304, both of which are closely linked distal markers for fragile X. This regional localization was supported by the analysis, in fragile X families, of recombination events between these three loci, the fragile X locus and locus DXS52, the latter being located at a more distal position. DXS466 is closely linked to the fragile X locus with a peak lod score of 7.79 at a recombination fraction of 0.02. Heterozygosity of DXS466 is approximately 50%. Its close proximity and relatively high informativity make DXS466 a valuable new diagnostic DNA marker for fragile X.     

A refined genetic map of the region of chromosome 17 surrounding the von recklinghausen neurofibromatosis (NF1) gene

The von Recklinghausen neurofibromatosis (NF1) gene has been mapped to the pericentromeric region of chromosome 17. We conducted linkage analyses of NF1 by using 10 polymorphic DNA markers from this chromosomal region. We ascertained 20 American Caucasian NF1 families (163 individuals, 98 NF1 affected) in Michigan and Ohio and also studied a large family ascertained primarily in North Carolina. The following markers were used in this study: HHH202, TH17.19, D17Z1, ERBA1, EW203, EW206, EW207, EW301, CRI-L581, and CRI-L946. NF1 did not recombine with either TH17.19 or HHH202 in any of the informative meioses surveyed (maximum lod scores of 17.04 and 7.21, respectively, at a recombination fraction of .00), indicating that these markers map very close to the NF1 gene. We also report evidence of three instances of recombination between NF1 and the centromeric marker D17Z1 (maximum lod score of 13.43 at a recombination fraction of .04), as well as two crossovers between pairs of marker loci. We find no evidence of locus heterogeneity, and our results support the localization of the NF1 gene to proximal chromosome 17q.     

A new polymorphic marker very closely linked to DXS52 in the q28 region of the human X chromosome

Summary We have isolated an X chromosome probe, St35.691 (DXS305), which detects two RFLPs with TaqI and PstI, whose combined heterozygosity is about 60%. This probe has been assigned to Xq28 by physical and genetic mapping and is very closely linked to DXS52, DXS15, and the coagulation factor VIII gene (F8C). The best estimate of the recombination fraction for the DXS52-DXS305 interval is 0.014, with a lod score of 50.1. Multipoint analysis places DXS305 on the same side of F8C as DXS52, but complete ordering of the three loci was not possible with our present data. This highly informative marker should be useful in the precise mapping of the many disease genes that have been assigned to the Xq28 band.     

Genetic mapping of autosomal dominant Charcot-Marie-Tooth disease in a large French-Acadian kindred: identification of new linked markers on chromosome 17.

We have performed linkage analysis in a large French-Acadian kindred segregating one form of autosomal dominant Charcot-Marie-Tooth disease (CMTD) (type IA) using 17 polymorphic DNA markers spanning human chromosome 17 and demonstrate linkage to several markers in the pericentromeric region, including DNA probes pA10-41, EW301, S12-30, pTH17.19, c11-2B, and p11-2c11.5. Linkage of markers pA10-41 and EW301 to CMTD type IA has been reported elsewhere. Four new markers, 1516, 1517, 1541, and LL101, which map to chromosome 17 have been identified. The marker 1516 appears to be closely linked to the CMTD locus on chromosome 17 as demonstrated by a maximum lod score of 3.42 at theta (recombination fraction) = 0. This marker has been mapped to 17p11.2 using a somatic cell hybrid constructed from a patient with Smith-Magenis syndrome [46,XY, del(17)(p11.2p11.2)]. A lod score of 6.16 has been obtained by multipoint linkage analysis with 1516 and two markers from 17q11.2, pTH17.19, and c11-2B. The markers 1517 and 1541 have been mapped to 17p12-17q11.2 and demonstrate maximum lod scores of 2.35 and 0.63 at recombination values of .1 and .2, respectively. The marker LL101 has been mapped to 17p13.105-17p13.100 and demonstrates a maximum lod score of 1.56 at a recombination value of .1. Our study confirms the localization of CMTD type IA to the pericentromeric region of chromosome 17.     

The comparison of RFLP,RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis

The utility of RFLP (restriction fragment length polymorphism), RAPD (random-amplified polymorphic DNA), AFLP (amplified fragment length polymorphism) and SSR (simple sequence repeat, microsatellite) markers in soybean germplasm analysis was determined by evaluating information content (expected heterozygosity), number of loci simultaneously analyzed per experiment (multiplex ratio) and effectiveness in assessing relationships between accessions. SSR markers have the highest expected heterozygosity (0.60), while AFLP markers have the highest effective multiplex ratio (19). A single parameter, defined as the marker index, which is the product of expected heterozygosity and multiplex ratio, may be used to evaluate overall utility of a marker system. A comparison of genetic similarity matrices revealed that, if the comparison involved both cultivated (Glycine max) and wild soybean (Glycine soja) accessions, estimates based on RFLPs, AFLPs and SSRs are highly correlated, indicating congruence between these assays. However, correlations of RAPD marker data with those obtained using other marker systems were lower. This is because RAPDs produce higher estimates of interspecific similarities. If the comparisons involvedG. max only, then overall correlations between marker systems are significantly lower. WithinG. max, RAPD and AFLP similarity estimates are more closely correlated than those involving other marker systems.Abbreviations RFLP restriction fragment length plymorphism - RAPD random-amplified polymorphic DNA - AFLP amplified fragment length polymorphism - SSR simple sequence repeat - PCR polymerase chain reaction - TBE Tris-borate-EDTA buffer - MI marker index - SENA sum of effective numbers of alleles     

A new deletion polymorphism at D5S71 raises the linkage information on adenomatous polyposis coli: implications for presymptomatic diagnosis

Summary Two independent study-groups, one in Britain and the other in the United States, were the first to report linkage between APC and a TaqI restriction fragment length polymorphism (RFLP) at D5S71 (probe C11p11) on chromosome 5q. They found no recombinants in about 50 informative meioses. The same TaqI RFLP was found to be uninformative for linkage in 15 Dutch polyposis families. The recently reported four base-pair deletion polymorphism (DEL1) at D5S71 has raised the polymorphism information content of this marker from 0.17 to 0.40 in the Dutch population. Seven of 20 polyposis families screened for the DEL1 as well as the TaqI polymorphism gave a combined peak lod score of 5.68 with no recombinants in 37 informative meioses. These data, together with those so far reported in the literature, raise the peak lod score to 17.09 at a recombination fraction of 0.05, the 95% upper confidence limit being 0.09. In combination with the use of another informative marker, D5S81 (probe YN5.48) closely mapping on the other side of APC, the presymptomatic diagnosis of the disease can be made with more than 99.9% certainty. It has to be stressed, however, that the the possible existence of more than one polyposis locus cannot, as yet, be excluded.     

DXS539, a polymorphic DNA marker proximal of the fragile-X gene

We report a new polymorphic DNA marker (pJH89, DXS539) proximal to the fragile-X site. The pJH89 probe identifies a TaqI and a NcoI restriction fragment length polymorphism (combined heterozygosity of 42%) and is linked to the fragile-X locus with a maximal LOD score of 12 at 4 cM. Multipoint linkage analysis and physical mapping studies indicate that the pJH89 probe is located within the interval defined by the markers DXS369 and DXS548.     

Hereditary neuralgic amyotrophy: evidence for genetic homogeneity and mapping to chromosome 17q25

Hereditary neuralgic amyotrophy (HNA) is a rare autosomal dominant disorder on chromosome 17q, associated with recurrent, episodic, painful brachial plexus neuropathy. Dysmorphic features, including hypotelorism, long nasal bridge and facial asymmetry, are frequently associated with HNA. To assess genetic homogeneity, determine the cytogenetic location, and identify flanking markers for the HNA locus, six pedigrees were studied with multiple DNA markers from distal chromosome 17q. The results in all pedigrees supported linkage of the HNA locus to chromosome 17. A maximum combined lod score (Ζ = 10.94, ￡ = 0.05) was obtained with marker D17S939 and the maximum multipoint lod score was 22.768 in the interval defined by D17S802– D17S939. An analysis of crossovers placed the HNA locus within an approximate 4.0-cM interval flanked by D17S1603 and D17S802. Analysis of DNA from a human/mouse somatic cell hybrid with linked markers suggests that band 17q25 harbors the HNA locus. These results support genetic homogeneity within HNA and define a specific interval and a precise cytogenetic location in chromosome 17q25 for this disorder. Received: 24 June 1997 / Accepted: 21 August 1997     

The genetic locus for free sialic acid storage disease maps to the long arm of chromosome 6.

Salla disease (SD), or adult-type free sialic acid storage disease, is an autosomal recessive lysosomal storage disorder characterized by impaired transport of free sialic acid across the lysosomal membrane and severe psychomotor retardation. Random linkage analysis of a sample of 27 Finnish families allowed us to localize the SD locus to the long arm of chromosome 6. The highest lod score of 8.95 was obtained with a microsatellite marker of locus D6S286 at theta = .00. Evidence for linkage disequilibrium was observed between the SD locus and the alleles of three closely linked markers, suggesting that the length of the critical region for the SD locus is in the order of 190 kb.     

A gene for cleidocranial dysplasia maps to the short arm of chromosome 6.

Cleidocranial dysplasia (CCD) is an autosomal dominant generalized bone dysplasia characterized by mild-to-moderate short stature, clavicular aplasia or hypoplasia, supernumerary and ectopic teeth, delayed eruption of secondary teeth, a characteristic craniofacial appearance, and a variety of other skeletal anomalies. We have performed linkage studies in five families with CCD, with 24 affected and 20 unaffected individuals, using microsatellite markers spanning two candidate regions on chromosomes 8q and 6. The strongest support for linkage was with chromosome 6p microsatellite marker D6S282 with a two-point lod score of 4.84 (theta = .03). Furthermore, the multipoint lod score was 5.70 in the interval between D6S282 and D6S291. These data show that the gene for autosomal dominant CCD is located within a 19-cM interval on the short arm of chromosome 6, between D6S282 and D6S291.     

Evaluation of designs for reference families for livestock linkage mapping experiments

Abstract

The development of dense linkage maps consisting of highly polymorphic loci for livestock species is technically feasible. However, linkage mapping experiments are expensive as they involve many animals and marker typings per animal. To minimize costs of developing linkage maps for livestock species, optimizing designs for mapping studies is necessary. This study provides a general framework for evaluating the efficiency of designs for reference families consisting of two‐ or three‐ generation full‐sib or half‐sib families selected from a segregating population. The influence of number of families, number of offspring per family, family structure (either half‐sib or full‐sib) and marker polymorphism is determined. Evaluation is done for two markers with a recombination rate of .20 and for a marker and a dominant single gene with a recombination rate of .20. Two evaluation criteria are used: expected maximum lod score for detection of linkage and accuracy of an estimated recombination rate defined as probability that the true recombination rate is in an interval around the estimated recombination rate. First, for several designs the contribution of reference families to expected maximum lod score and accuracy is given. Second, the required number of families in a design to obtain a certain value for the evaluation criteria is calculated when number of offspring per family, family structure and marker polymorphism are specified. The required numbers increase when designs are optimized not only for expected maximum lod score but also for accuracy. The required number of animals to map a dominant single gene is very large. Therefore, a set of reference families should be designed for strictly mapping marker loci. Examples illustrate how tabulated results can be generalized to determine the values for a wide range of designs containing two‐ or three‐generation full‐sib or half‐sib families.     

The genes for apolipoprotein all (APOA2) and the Duffy blood group (FY) are linked on chromosome 1 in man

We have cloned a cDNA probe for human apolipoprotein AII and used it to analyze linkage relationships on chromosome 1. We found no recombinations between APOA2 and the gene coding for the Duffy blood group antigens (FY) in the 19 meioses examined. Our maximal lod score is 4.2 at zero recombination rate. K. Berg (1987, Cytogenet. Cell Genet. 46:579) found a maximal score of 2.5 at recombination fraction 0.14 in 54 meioses. When results from both studies are combined, the most likely distance between FY and APOA2 is about 10% recombination with a combined lod score of 5.6 for both sexes.     

A mapped set of genetic markers for human chromosome 9

A genetic map of markers for human chromosome 9, spanning a genetic distance of 147 cM in males and 231 cM in females, has been constructed from linkage studies with 19 loci in a large panel of reference families. The markers included four classical systems previously assigned to chromosome 9, and restriction fragment length polymorphisms of two cloned genes, ABL oncogene and argininosuccinase synthetase pseudogene 3 (ASSP3). The remaining 13 marker loci, with an average heterozygosity of 42%, were defined by arbitrary DNA probes newly ascertained from genomic libraries; seven of them were variable number of tandem repeat (VNTR) loci. A subset of 7 of the 19 linked markers is proposed for a primary map that could detect linkage with a genetic defect within the covered region of chromosome 9, provided that at least 45 phase-known meioses were available for study in an affected family.