A simple note on how to save money in linkage analysis

Genetic linkage maps are often based on maximum-likelihood estimates of recombination fractions which are converted into map units by mapping functions. This paper presents a cost analysis of linkage analysis for a segregating F₂␣population with codominant or dominant molecular markers and a qualitative monogenic dominant–recessive trait. For illustration, a disease-resistance trait is considered, where the susceptible allele is recessive. Three sub-populations of the F₂ can be used for linkage analysis [susceptible (= recessive) individuals, resistant (= dominant) individuals, complete F₂]. While it is well-known that recessive individuals are more informative than dominant individuals, it is not obvious a priori, which of the three sub-populations should be preferred, when costs of phenotyping and genotyping are taken into consideration. A comparative economic analysis of alternative procedures of linkage detection based on these three sub-populations does exhibit a clear economic superiority of the sub-population of susceptible (= recessive) individuals, when costs of genotyping are high. This cost-effectiveness is due to the higher information content of this sub-population compared to the sub-population of dominant (= resistant) individuals and also compared to the complete F₂. Our final conclusion/recommendation is as follows: If the cost to genotype an individual is sufficiently large compared with the cost to phenotype an individual, then linkage analysis and genetic mapping should be only based on susceptible (= recessive) individuals. Conversely, if the cost of phenotyping exceeds that for genotyping, it may be preferable to genotype all plants. The exact conditions under which a strategy is preferable are described in the paper.     

The autosomal dominant familial exudative vitreoretinopathy locus maps on 11q and is closely linked to D11S533.

Autosomal dominant familial exudative vitreoretinopathy (adFEVR) is a hereditary disorder characterized by the incomplete vascularization of the peripheral retina. The primary biochemical defect in adFEVR is unknown. The adFEVR locus has tentatively been assigned to 11q by linkage studies. We report the results of an extended multipoint linkage analysis of two families with adFEVR by using five markers (INT2, D11S533, D11S527, D11S35, and CD3D) from 11q13-q23. Pairwise linkage data obtained in the two families were rather similar and hence have not provided evidence for genetic heterogeneity. The highest complied two-point lod score (3.67, at a recombination fraction of .07) was obtained for the disease locus versus D11S533. Multipoint analyses showed that the adFEVR locus maps most likely, with a maximum location score of over 20, between D11S533/D11S527 and D11S35, at recombination rates of .147 and .104, respectively. Close linkage without recombination (maximum lod score 11.26) has been found between D11S533 and D11S527.     

Determining the sample size for co-dominant molecular marker-assisted linkage detection for a monogenic qualitative trait by controlling the type-I and type-II errors in a segregating F2 population

Tests for linkage are usually performed using the lod score method. A critical question in linkage analyses is the choice of sample size. The appropriate sample size depends on the desired type-I error and power of the test. This paper investigates the exact type-I error and power of the lod score method in a segregating F(2) population with co-dominant markers and a qualitative monogenic dominant-recessive trait. For illustration, a disease-resistance trait is considered, where the susceptible allele is recessive. A procedure is suggested for finding the appropriate sample size. It is shown that recessive plants have about twice the information content of dominant plants, so the former should be preferred for linkage detection. In some cases the exact alpha-values for a given nominal alpha may be rather small due to the discrete nature of the sampling distribution in small samples. We show that a gain in power is possible by using exact methods.     

Evidence for heterogeneity in facioscapulohumeral muscular dystrophy (FSHD)

Facioscapulohumeral muscular dystrophy (FSHD) is a slowly progressive primary disease of muscle which is usually inherited as an autosomal dominant disorder. FSHD has been localized to the long arm of chromosome 4, specifically to the 4q3.5-qter region. Initially published linkage studies showed no evidence for heterogeneity in FSHD. In the present study we have examined individuals in seven FSHD families. Two-point lod scores show significant evidence for linkage for D4S163 (lod score 3.04 at recombination fraction .21) and D4S139 (lod score 3.84 at recombination fraction .20). D4S171 also gave a positive score (lod score 2.56 at recombination fraction .24). Significant evidence for heterogeneity was found for each of the three markers. Multipoint linkage analysis in this region resulted in a peak multipoint lod score of 6.47. The multipoint analysis supported the two-point studies with odds of 20:1 showing linkage and heterogeneity over linkage and homogeneity. Five of the seven families gave a posterior probability of >95% of being of the linked type, while two families appeared unlinked to this region of 4q (P < .01%). Individuals in the two unlinked families met the clinical criteria for the diagnosis of FSHD, including facial weakness, clavicular flattening, scapula winging, proximal muscle weakness, and myopathic changes on muscle biopsies without inflammatory or mitochondrial pathology. This study demonstrates genetic heterogeneity in FSHD and has important implications for both genetic counseling and the elucidation of the etiology of FSHD.     

Biased estimation of the recombination fraction using half-sib families and informative offspring

A maximum-likelihood method to estimate the recombination fraction and its sampling variance using informative and noninformative half-sib offspring is derived. Estimates of the recombination fraction are biased up to 20 cM when noninformative offspring are discarded. In certain scenarios, the sampling variance can be increased or reduced up to fivefold due to the bias in estimating the recombination fraction and the LOD score can be reduced up to 5 units when discarding noninformative offspring. Comparison of the estimates of recombination fraction, map distance, and LOD score when constructing a genetic map with 251 two-point linkage analyses and six families of Norwegian cattle was carried out to evaluate the implications of discarding noninformative offspring in practical situations. The average discrepancies in absolute value (average difference when using and neglecting noninformative offspring) were 0.0146, 1.64 cM, and 2.61 for the recombination fraction, map distance, and the LOD score, respectively. A method for simultaneous estimation of allele frequencies in the dam population and a transmission disequilibrium parameter is proposed. This method might account for the bias in estimating allele frequencies in the dam population when the half-sib offspring is selected for production traits.     

Further evidence for localization of the gene of erythrokeratodermia variabilis

Summary Close genetic linkage between erythrokeratodermia variabilis (EKV) and the Rh blood group system has been reported by our group. Here we describe the results of a linkage analysis in another EKV kindred, in which the disease segregated with the CDe genotype. Among 18 informative individuals, 1 recombinant was found. A maximum lod score of 4.21 was calculated at a recombination fraction of 0.03–0.04. Addition of this lod score to the earlier reported results gives a maximum lod score of 9.93 for linkage between EKV and Rh at a recombination fraction of 0.03 (95% confidence limits 0.008–0.11).     

Localisation of a gene for Papillon-Lefèvre syndrome to chromosome 11q14–q21 by homozygosity mapping

Papillon-Lefèvre syndrome is an autosomal recessively inherited palmoplantar keratoderma of unknown aetiology associated with severe periodontitis leading to premature loss of dentition. Three consanguineous families, two of Turkish and one of German origin, and three multiplex families, one of Ethiopian and two of German origin, with 11 affected and 6 unaffected siblings in all were studied. A targeted genome search was initially attempted to several candidate gene regions but failed to demonstrate linkage. Therefore a genome-wide linkage scan using a combination of homozygosity mapping and traditional linkage analysis was undertaken. Linkage was obtained with marker D11S937 with a maximum two-point lod score of Z _max = 6.1 at recombination fraction θ = 0.00 on chromosome 11q14–q21 near the metalloproteinase gene cluster. Multipoint likelihood calculations gave a maximum lod score of 7.35 between D11S901 and D11S1358. A 9.2-cM region homozygous by descent in the affected members of the three consanguineous families lies between markers D11S1989 and D11S4176 harbouring the as yet unknown Papillon-Lefèvre syndrome gene. Haplotype analyses in all the families studied support this localisation. This study has identified a further locus harbouring a gene for palmoplantar keratoderma and one possibly involved in periodontitis. Received: 19 July 1997 / Accepted: 22 August 1997     

Mapping of the locus for X-linked cardioskeletal myopathy with neutropenia and abnormal mitochondria (Barth syndrome) to Xq28.

X-linked cardioskeletal myopathy with neutropenia and abnormal mitochondria is clinically characterized by congenital dilated cardiomyopathy, skeletal myopathy, recurrent bacterial infections, and growth retardation. We analyzed linkage between the disease locus and X-chromosomal markers in a family with seven carriers, four patients, and eight unaffected sons of carriers. Highest lod scores obtained by two-point linkage analysis were 2.70 for St14.1 (DXS52, TaqI) at a recombination fraction of zero and 2.53 for cpX67 (DXS134) at a recombination fraction of zero. Multipoint linkage analysis resulted in a maximum lod score of 5.24 at the position of St35.691 (DXS305). The most distal recombination detected in this family was located between the markers II-10 (DXS466) and DX13 (DXS15). These data indicate the location of the mutated gene at Xq28.     

Mapping of the gene for X-linked amelogenesis imperfecta by linkage analysis.

X-linked Amelogenesis imperfecta (AI) is a genetic disorder affecting the formation of enamel. In the present study two families, one with X-linked dominant and one with X-linked recessive AI, were studied by linkage analysis. Eleven cloned RFLP markers of known regional location were used. Evidence was obtained for linkage between the AI locus and the marker p782, defining the locus DXS85 at Xp22, by using two-point analysis. No recombination was scored between these two loci in 15 informative meioses, and a peak lod score (Zmax) of 4.45 was calculated at zero recombination fraction. Recombination was observed between the more distal locus DXS89 and AI, giving a peak lod score of 3.41 at a recombination fraction of .09. Recombination was also observed between the AI locus and the more proximal loci DXS43 and DXS41 (Zmax = 0.09 at theta max = 0.31 and Zmax = 0.61 at theta max = 0.28, respectively). Absence of linkage was observed between the AI locus and seven other loci, located proximal to DXS41 or on the long arm of the X chromosome. On the basis of two-point linkage analysis and analysis of crossover events, we propose the following order of loci at Xp22: DXS89-(AI, DXS85)-DXS43-DXS41-Xcen.     

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.     

Two families with familial amyotrophic lateral sclerosis are linked to a novel locus on chromosome 16q

Amyotrophic lateral sclerosis (ALS) is a fatal adult-onset disease in which motor neurons in the brain and spinal cord degenerate by largely unknown mechanisms. ALS is familial (FALS) in 10% of cases, and the inheritance is usually dominant, with variable penetrance. Mutations in copper/zinc super oxide dismutase (SOD1) are found in 20% of familial and 3% of sporadic ALS cases. Five families with ALS and frontotemporal dementia (ALS-FTD) are linked to 9q21, whereas one family with pure ALS is linked to 18q21. We identified two large European families with ALS without SOD1 mutations or linkage to known FALS loci and conducted a genomewide linkage screen using 400 microsatellite markers. In both families, two-point LOD scores >1 and a haplotype segregating with disease were demonstrated only across regions of chromosome 16. Subsequent fine mapping in family 1 gave a maximum two-point LOD score of 3.62 at D16S3137 and a three-point LOD score of 3.85 for markers D16S415 and D16S3137. Haplotype analysis revealed no recombination > approximately 30 cM, (flanking markers at D16S3075 and D16S3112). The maximum two-point LOD score for family 2 was 1.84 at D16S415, and the three-point LOD score was 2.10 for markers D16S419 and D16S415. Definite recombination occurred in several individuals, which narrowed the shared haplotype in affected individuals to a 10.1-cM region (flanking markers: D16S3396 and D16S3112). The region shared by both families on chromosome 16q12 corresponds to approximately 4.5 Mb on the Marshfield map. Bioinformatic analysis of the region has identified 18 known genes and 70 predicted genes in this region, and sequencing of candidate genes has now begun.     

Autosomal dominant macrothrombocytopenia with leukocyte inclusions (May-Hegglin anomaly) is linked to chromosome 22q12-13

Macrothrombocytopenia with leukocyte inclusions (May-Hegglin anomaly) is a rare autosomal dominant disorder characterized by thrombocytopenia, giant platelets, and D?hle body-like inclusions in leukocytes. To determine the genetic basis of this disorder, we performed a genome-wide screen for linkage in three families with May-Hegglin anomaly. For the pooled analysis of the three families, three markers on chromosome 22 had two-point logarithm-of-difference (lod) scores greater than 3, with a maximum lod score of 3.91 at a recombination fraction (theta) of 0.076 for marker D22S683. Within the largest family (MHA-1), the maximum lod score was 5.36 at theta=0 at marker D22S445. Fine mapping of recombination events using eight adjacent markers indicated that the minimal disease region of family MHA-1 alone is in the approximately 26 cM region from D22S683 to the telomere. The maximum lod score for the three families combined was 5.84 at theta=0 for marker IL2RB. With the assumption of locus homogeneity, haplotype analysis of family MHA-4 indicated the disease region is centromeric to marker D22S1045. These data best support a minimal disease region from D22S683 to D22S1045, a span of about 1 Mb of DNA that contains 17 known genes and 4 predicted genes. Further analysis of this region will identify the genetic basis of May-Hegglin anomaly, facilitating subsequent characterization of the biochemical role of the disease gene in platelet formation.     

Genetic heterogeneity in X-linked amelogenesis imperfecta.

The AMELX gene located at Xp22.1-p22.3 encodes for the enamel protein amelogenin and has been implicated as the gene responsible for the inherited dental abnormality X-linked amelogenesis imperfecta (XAI). Three families with XAI have been investigated using polymorphic DNA markers flanking the position of AMELX. Using two-point linkage analysis, linkage was established between XAI and several of these markers in two families, with a combined lod score of 6.05 for DXS16 at theta = 0.04. This supports the involvement of AMELX, located close to DXS16, in the XAI disease process (AIH1) in those families. Using multipoint linkage analysis, the combined maximum lod score for these two families was 7.30 for a location of AIH1 at 2 cM distal to DXS16. The support interval around this location extended about 8 cM proximal to DXS92, and the AIH1 location could not be precisely defined by multipoint mapping. Study of recombination events indicated that AIH1 lies in the interval between DXS143 and DXS85. There was significant evidence against linkage to this region in the third family, indicating locus heterogeneity in XAI. Further analysis with markers on the long arm of the X chromosome showed evidence of linkage to DXS144E and F9 with no recombination with either of these markers. Two-point analysis gave a peak lod score at DXS144E with a maximum lod score of 2.83 at theta = 0, with a peak lod score in multipoint linkage analysis of 2.84 at theta = 0. The support interval extended 9 cM proximal to DXS144E and 14 cM distal to F9.(ABSTRACT TRUNCATED AT 250 WORDS)     

结合SADF与RAPD标记构建家蚕连锁图

用40个选择性扩增多态性引物和137个随机性扩增多态性引物对家蚕Bombyx mori F₂群体进行分子标记的筛选。获得544个符合遗传分离比的多态性位点,并用Mapmaker软件进行连锁分析。在最小LOD值为4,最大重组值为0.2条件下拆分连锁群,并对处于同一连锁群的位点进行合理排序,计算出位点间的重组值后转换为图谱距离,构建了一个家蚕的分子连锁图。     

Suggestive evidence for linkage of schizophrenia to markers on chromosome 13 in Caucasian but not Oriental populations

Previously we reported suggestive evidence for linkage of schizophrenia to markers on chromosome 13q14.1–q32. We have now studied an additional independent sample of 44 pedigrees consisting of 34 Taiwanese, 9 English and 1 Welsh family in an attempt to replicate this finding. Narrow and broad models based on Research Diagnostic Criteria or the Diagnostic and Statistical Manual of Mental Disorders, third edition, revised, were used to define the schizophrenia phenotype. Under a dominant genetic model, two-point lod scores obtained for most of the markers were negative except that marker D13S122 gave a total lod score of 1.06 (θ = 0.2, broad model). As combining pedigrees from different ethnic origins may be inappropriate, we combined this replication sample and our original sample, and then divided the total sample into Caucasian (English and Welsh pedigrees) and Oriental (Taiwanese and Japanese pedigrees) groups. The Caucasian pedigrees produced maximized admixture two-point lod scores (A-lod) of 1.41 for the marker D13S119 (θ = 0.2, α = 1.0) and 1.54 for D13S128 (θ = 0, α = 0.3) with nearby markers also producing positive A-lod scores. When five-point model-free linkage analysis was applied to the Caucasian sample, a maximum lod score of 2.58 was obtained around the markers D13S122 and D13S128, which are located on chromosome 13q32. The linkage results for the Oriental group were less positive than the Caucasian group. Our results again suggest that there is a potential susceptibility locus for schizophrenia on chromosome 13q14.1–q32, especially in the Caucasian population. Received: 13 September 1996     

Additional glomangioma families link to chromosome 1p: no evidence for genetic heterogeneity

Venous malformations are a common abnormality of the vasculature that may occur sporadically or, more rarely, as an autosomal dominant trait. One familial form of venous malformations has previously been linked to chromosome 9p. Mutations in the gene encoding Tie2, an endothelial specific receptor tyrosine kinase, have been identified in four different families. Glomangiomas are a subtype of venous malformations with glomus cell involvement. These cutaneous lesions can be inherited as an autosomal dominant disease with reduced penetrance and variable expressivity. We present evidence of linkage to chromosome 1p21-1p22 using four new glomangioma families, with a combined maximum two-point lod score of 7.32 at marker D1S2804. Markers D1S2129 and D1S2881 define the 24-cM linkage interval determined by recombination within affected individuals. A recent report also showed linkage of the glomangioma locus to chromosome 1p. A total of 9 families now map to this region, suggesting a decreased likelihood of locus heterogenity in familial glomangiomas. Investigation of candidate genes within the interval should provide new insights into lesion formation in inherited venous malformations.     

Confirmation of linkage to 1q21-31 in a Danish autosomal dominant juvenile-onset glaucoma family and evidence of genetic heterogeneity

Autosomal dominant juvenile-onset open-angle glaucoma has been mapped to 1q21-31 in a number of American families. Our study confirms linkage in a Danish five-generation dominant juvenile-onset glaucoma family with a maximum two-point lod score of 6.67 at the D1S210 locus. Multipoint linkage analysis in a nine-generation Swedish family with dominant juvenile-onset glaucoma and iris hypoplasia excludes linkage to the region of approximately 18 cM between loci D1S104 and D1S218, shown to contain the previously mapped glaucoma gene. This study thus provides support for genetic heterogeneity with respect to dominant juvenile-onset glaucoma.     

Mapping of the gene for G blood group antigens to chromosome 15 in swine*

Linkage analysis between size variants of the centromeric region of chromosome 15 and G blood group alleles produced a lod score of 5.03. The maximum likelihood estimate of the recombination fraction is θ= 0.24 with a 95 % confidence interval of 0.08 < θ < 0.40. Since this is the second time that linkage between a chromosome 15 marker and the G blood group locus has been shown, the assignment of the G blood group locus to chromosome 15 is confirmed.     

Criteria to optimize designs for detection and estimation of linkage between marker loci from segregating populations containing several families

Summary Construction of a genome map of highly polymorphic markers has become possible in the past decade. Establishing a complete marker map is an enormous task. Therefore, designs to map molecular markers should be optimal. Designs to detect and estimate linkage between markers from segregating populations were studied. Two measures of design quality were used. The expectation of the maximum lod score indicates the possibility of designs to detect linkage. The accuracy of estimating recombination rate was measured as the probability that the true recombination rate is in a specified internal given the estimate. Accurate approximate methods were developed for rapid evaluation of designs. Seven family types (e.g., double backcross) can be distinguished that describe all families in a segregating population. The family type influences the expected maximum lod score and the accuracy of estimation. The frequency of favorable family types increased with increasing marker polymorphism. At a true recombination rate of 0.20,27 observations on offspring when five alleles were segregating, and 55 observations on offspring when two alleles were segregating, were necessary to obtain an expected maximum lod score of 3. The probability that the true recombination rate was between 0.15 and 0.25, given an estimate of 0.20, was about 0.85 for a design with 40 families with ten offspring and two alleles segregating and for a design with ten families with ten offspring and six alleles segregating. For smaller designs, accuracies were less, approximate evaluation of accuracy was not justified and, on average, true recombination rates were much greater than estimated given a specified value for the estimated recombination rate.     

Lack of association and linkage between HLA and familial polyposis coli

Summary We investigated possible association of and linkage between HLA and familial polyposis coli (FPC). In 182 individuals from 66 pedigrees of FPC and 108 individuals from a normal population, HLA-A,-B, and-C antigens were determined. When the frequencies of HLA antigens in 66 unrelated patients and in normal controls were compared, no association of FPC with HLA was observed. For the linkage analysis, HLA haplotypes of 17 affected sib pairs were investigated by the affected sib pair method. The number of pairs which shared two, one, and no haplotypes identical by descent was not significantly different from the number expected with random occurrence (P>0.95). Finally, seven families were analyzed using Morton's sequential test. A maximum lod score of-0.056 at a recombination fraction of 0.4, and a lod of-3.089 at a recombination fraction of 0.05 were obtained. Therefore, there is neither an association of nor linkage between FPC and HLA.