Linkage analysis in the presence of errors I: complex-valued recombination fractions and complex phenotypes

Linkage is a phenomenon that correlates the genotypes of loci, rather than the phenotypes of one locus to the genotypes of another. It is therefore necessary to convert the observed trait phenotypes into trait-locus genotypes, which can then be analyzed for coinheritance with marker-locus genotypes. However, if the mode of inheritance of the trait is not known accurately, this conversion can often result in errors in the inferred trait-locus genotypes, which, in turn, can lead to the misclassification of the recombination status of meioses. As a result, the recombination fraction can be overestimated in two-point analysis, and false exclusions of the true trait locus can occur in multipoint analysis. We propose a method that increases the robustness of multipoint analysis to errors in the mode of inheritance assumptions of the trait, by explicitly allowing for misclassification of trait-locus genotypes. To this end, the definition of the recombination fraction is extended to the complex plane, as Theta=straight theta+straightepsiloni; theta is the recombination fraction between actual ("real") genotypes of marker and trait loci, and straightepsilon is the probability of apparent but false ("imaginary") recombinations between the actual and inferred trait-locus genotypes. "Complex" multipoint LOD scores are proven to be stochastically equivalent to conventional two-point LOD scores. The greater robustness to modeling errors normally associated with two-point analysis can thus be extended to multiple two-point analysis and multipoint analysis. The use of complex-valued recombination fractions also allows the stochastic equivalence of "model-based" and "model-free" methods to be extended to multipoint analysis.     

Molecular and statistical approaches to the detection and correction of errors in genotype databases.

Errors in genotyping data have been shown to have a significant effect on the estimation of recombination fractions in high-resolution genetic maps. Previous estimates of errors in existing databases have been limited to the analysis of relatively few markers and have suggested rates in the range 0.5%-1.5%. The present study capitalizes on the fact that within the Centre d'Etude du Polymorphisme Humain (CEPH) collection of reference families, 21 individuals are members of more than one family, with separate DNA samples provided by CEPH for each appearance of these individuals. By comparing the genotypes of these individuals in each of the families in which they occur, an estimated error rate of 1.4% was calculated for all loci in the version 4.0 CEPH database. Removing those individuals who were clearly identified by CEPH as appearing in more than one family resulted in a 3.0% error rate for the remaining samples, suggesting that some error checking of the identified repeated individuals may occur prior to data submission. An error rate of 3.0% for version 4.0 data was also obtained for four chromosome 5 markers that were retyped through the entire CEPH collection. The effects of these errors on a multipoint map were significant, with a total sex-averaged length of 36.09 cM with the errors, and 19.47 cM with the errors corrected. Several statistical approaches to detect and allow for errors during linkage analysis are presented. One method, which identified families containing possible errors on the basis of the impact on the maximum lod score, showed particular promise, especially when combined with the limited retyping of the identified families. The impact of the demonstrated error rate in an established genotype database on high-resolution mapping is significant, raising the question of the overall value of incorporating such existing data into new genetic maps.     

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

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

The genetics of reading disability in an often excluded sample: novel loci suggested for reading disability in rolandic epilepsy

Background

Reading disability (RD) is a common neurodevelopmental disorder with genetic basis established in families segregating “pure” dyslexia. RD commonly occurs in neurodevelopmental disorders including Rolandic Epilepsy (RE), a complex genetic disorder. We performed genomewide linkage analysis of RD in RE families, testing the hypotheses that RD in RE families is genetically heterogenenous to pure dyslexia, and shares genetic influences with other sub-phenotypes of RE.

Methods

We initially performed genome-wide linkage analysis using 1000 STR markers in 38 US families ascertained through a RE proband; most of these families were multiplex for RD. We analyzed the data by two-point and multipoint parametric LOD score methods. We then confirmed the linkage evidence in a second US dataset of 20 RE families. We also resequenced the SEMA3C gene at the 7q21 linkage locus in members of one multiplex RE/RD pedigree and the DISC1 gene in affected pedigrees at the 1q42 locus.

Results

In the discovery dataset there was suggestive evidence of linkage for RD to chromosome 7q21 (two-point LOD score 3.05, multipoint LOD 3.08) and at 1q42 (two-point LOD 2.87, multipoint LOD 3.03). Much of the linkage evidence at 7q21 derived from families of French-Canadian origin, whereas the linkage evidence at 1q42 was well distributed across all the families. There was little evidence for linkage at known dyslexia loci. Combining the discovery and confirmation datasets increased the evidence at 1q42 (two-point LOD = 3.49, multipoint HLOD = 4.70), but decreased evidence at 7q21 (two-point LOD = 2.28, multipoint HLOD  = 1.81), possibly because the replication sample did not have French Canadian representation.

Discussion

Reading disability in rolandic epilepsy has a genetic basis and may be influenced by loci at 1q42 and, in some populations, at 7q21; there is little evidence of a role for known DYX loci discovered in “pure” dyslexia pedigrees. 1q42 and 7q21 are candidate novel dyslexia loci.     

Identification of a Major Susceptibility Locus for Restless Legs Syndrome on Chromosome 12q

Restless legs syndrome (RLS) is a neurological disorder characterized by leg paresthesia associated with an irresistible urge to move that often interferes with nocturnal sleep, leading to chronic sleep deprivation. To map genes that may play a role in the vulnerability to RLS, a genomewide scan was conducted in a large French-Canadian family. Significant linkage was established on chromosome 12q, for a series of adjacent microsatellite markers with a maximum two-point LOD score of 3.42 (recombination fraction.05; P=6x10(-4); autosomal recessive mode of inheritance), whereas multipoint linkage calculations yielded a LOD score of 3.59. Haplotype analysis refined the genetic interval, positioning the RLS-predisposing gene in a 14.71-cM region between D12S1044 and D12S78. These findings represent the first mapping of a locus conferring susceptibility to RLS.     

TetraploidMap: construction of a linkage map in autotetraploid species

TetraploidMap is a suite of Fortran 90 routines run from Microsoft Windows with a text-based input and output. TetraploidMap enables the user to assemble a linkage map from dominant and codominant (multiallelic) marker loci scored for the parents and full-sib progeny of a cross in an autotetraploid species. It includes routines for the inference of the parental genotypes, identification of linkage groups, two-point analysis to estimate the recombination frequency and LOD score between all pairs of marker in a linkage group, and locus ordering by simulated annealing.     

Autosomal dominant distal myopathy: linkage to chromosome 14.

We have studied a family segregating a form of autosomal dominant distal myopathy (MIM 160500) and containing nine living affected individuals. The myopathy in this family is closest in clinical phenotype to that first described by Gowers in 1902. A search for linkage was conducted using microsatellite, VNTR, and RFLP markers. In total, 92 markers on all 22 autosomes were run. Positive linkage was obtained with 14 of 15 markers tested on chromosome 14, with little indication of linkage elsewhere in the genome. Maximum two-point LOD scores of 2.60 at recombination fraction .00 were obtained for the markers MYH7 and D14S64--the family structure precludes a two-point LOD score > or = 3. Recombinations with D14S72 and D14S49 indicate that this distal myopathy locus, MPD1, should lie between these markers. A multipoint analysis assuming 100% penetrance and using the markers D14S72, D14S50, MYH7, D14S64, D14S54, and D14S49 gave a LOD score of exactly 3 at MYH7. Analysis at a penetrance of 80% gave a LOD score of 2.8 at this marker. This probable localization of a gene for distal myopathy, MPD1, on chromosome 14 should allow other investigators studying distal myopathy families to test this region for linkage in other types of the disease, to confirm linkage or to demonstrate the likely genetic heterogeneity.     

Attention-deficit/hyperactivity disorder in a population isolate: linkage to loci at 4q13.2, 5q33.3, 11q22, and 17p11

Attention-deficit/hyperactivity disorder (ADHD [MIM 143465]) is the most common behavioral disorder of childhood. Twin, adoption, segregation, association, and linkage studies have confirmed that genetics plays a major role in conferring susceptibility to ADHD. We applied model-based and model-free linkage analyses, as well as the pedigree disequilibrium test, to the results of a genomewide scan of extended and multigenerational families with ADHD from a genetic isolate. In these families, ADHD is highly comorbid with conduct and oppositional defiant disorders, as well as with alcohol and tobacco dependence. We found evidence of linkage to markers at chromosomes 4q13.2, 5q33.3, 8q11.23, 11q22, and 17p11 in individual families. Fine mapping applied to these regions resulted in significant linkage in the combined families at chromosomes 4q13.2 (two-point allele-sharing LOD score from LODPAL = 4.44 at D4S3248), 5q33.3 (two-point allele-sharing LOD score from LODPAL = 8.22 at D5S490), 11q22 (two-point allele-sharing LOD score from LODPAL = 5.77 at D11S1998; multipoint nonparametric linkage [NPL]-log[P value] = 5.49 at approximately 128 cM), and 17p11 (multipoint NPL-log [P value] >12 at approximately 12 cM; multipoint maximum location score 2.48 [alpha = 0.10] at approximately 12 cM; two-point allele-sharing LOD score from LODPAL = 3.73 at D17S1159). Additionally, suggestive linkage was found at chromosome 8q11.23 (combined two-point NPL-log [P value] >3.0 at D8S2332). Several of these regions are novel (4q13.2, 5q33.3, and 8q11.23), whereas others replicate already-published loci (11q22 and 17p11). The concordance between results from different analytical methods of linkage and the replication of data between two independent studies suggest that these loci truly harbor ADHD susceptibility genes.     

Evidence for a prostate cancer-susceptibility locus on chromosome 20

Recent studies suggest that hereditary prostate cancer is a complex disease involving multiple susceptibility genes and variable phenotypic expression. While conducting a genomewide search on 162 North American families with > or =3 members affected with prostate cancer (PRCA), we found evidence for linkage to chromosome 20q13 with two-point parametric LOD scores >1 at multiple sites, with the highest two-point LOD score of 2.69 for marker D20S196. The maximum multipoint NPL score for the entire data set was 3.02 (P=.002) at D20S887. On the basis of findings from previous reports, families were stratified by the presence (n=116) or absence (n=46) of male-to-male transmission, average age of diagnosis (<66 years, n=73; > or =66 years, n=89), and number of affected individuals (<5, n=101; > or =5, n=61) for further analysis. The strongest evidence of linkage was evident with the pedigrees having <5 family members affected with prostate cancer (multipoint NPL 3.22, P=.00079), a later average age of diagnosis (multipoint NPL 3.40, P=.0006), and no male-to-male transmission (multipoint NPL 3.94, P=.00007). The group of patients having all three of these characteristics (n=19) had a multipoint NPL score of 3.69 (P=.0001). These results demonstrate evidence for a PRCA susceptibility locus in a subset of families that is distinct from the groups more likely to be linked to previously identified loci.     

Localization of a gene for benign adult familial myoclonic epilepsy to chromosome 8q23.3-q24.1.

Benign adult familial myoclonic epilepsy is an autosomal dominant idiopathic epileptic syndrome characterized by adult-onset tremulous finger movement, myoclonus, epileptic seizures, and nonprogressive course. It was recently recognized in Japanese families. In this study, we report that the gene locus is assigned to the distal long arm of chromosome 8, by linkage analysis in a large Japanese kindred with a maximum two-point LOD score of 4.31 for D8S555 at recombination fraction of 0 (maximum multipoint LOD score of 5.42 for the interval between D8S555 and D8S1779). Analyses of recombinations place the locus within an 8-cM interval, between D8S1784 and D8S1694, in which three markers, D8S1830, D8S555, and D8S1779, show no recombination with the phenotypes. Although three other epilepsy-related loci on chromosome 8q have been recognized-one on chromosome 8q13-21 (familial febrile convulsion) and two others on chromosome 8q24 (KCNQ3 and childhood absence epilepsy)-the locus assigned here is distinct from these three epilepsy-related loci. This study establishes the presence of a new epilepsy-related locus on 8q23.3-q24.11.     

Autosomal dominant avascular necrosis of femoral head in two Taiwanese pedigrees and linkage to chromosome 12q13

Avascular necrosis of the femoral head (ANFH) is a debilitating disease that commonly leads to destruction of the hip joint in adults. The etiology of ANFH is unknown, but previous studies have indicated that heritable thrombophilia (increased tendency to form thrombi) and hypofibrinolysis (reduced ability to lyse thrombi), alcohol intake, and steroid use are risk factors for ANFH. We recently identified two families with ANFH showing autosomal dominant inheritance. By applying linkage analysis to a four-generation pedigree, we excluded linkage between the family and three genes related to thrombophilia and hypofibrinolysis: protein C, protein S, and plasminogen activator inhibitor. Furthermore, by a genomewide scan, a significant two-point LOD score of 3.45 (recombination fraction [theta] = 0) was obtained between the family with ANFH and marker D12S85 on chromosome 12. High-resolution mapping was conducted in a second family with ANFH and replicated the linkage to D12S368 (pedigree I: LOD score 2.47, theta = 0.05; pedigree II: LOD score 2.81, theta = 0.10). When an age-dependent-penetrance model was applied, the combined multipoint LOD score was 6.43 between D12S1663 and D12S85. Thus, we mapped the candidate gene for autosomal dominant ANFH to a 15-cM region between D12S1663 and D12S1632 on chromosome 12q13.     

Linkage disequilibrium inflates type I error rates in multipoint linkage analysis when parental genotypes are missing

OBJECTIVES: Describe the inflation in nonparametric multipoint LOD scores due to inter-marker linkage disequilibrium (LD) across many markers with varied allele frequencies. METHOD: Using simulated two-generation families with and without parents, we conducted nonparametric multipoint linkage analysis with 2 to 10 markers with minor allele frequencies (MAF) of 0.5 and 0.1. RESULTS: Misspecification of population haplotype frequencies by assuming linkage equilibrium caused inflated multipoint LOD scores due to inter-marker LD when parental genotypes were not included. Inflation increased as more markers in LD were included and decreased as markers in equilibrium were added. When marker allele frequencies were unequal, the r2 measure of LD was a better predictor of inflation than D'. CONCLUSION: This observation strongly supports the evaluation of LD in multipoint linkage analyses, and further suggests that unaccounted for LD may be suspected when two-point and multipoint linkage analyses show a marked disparity in regions with elevated r2 measures of LD. Given the increasing popularity of high-density genome-wide SNP screens, inter-marker LD should be a concern in future linkage studies.     

Linkage of a Gene Causing High Bone Mass to Human Chromosome 11 (11q12-13)

The purpose of this paper is to report the linkage of a genetic locus (designated "HBM") in the human genome to a phenotype of very high spinal bone density, using a single extended pedigree. We measured spinal bone-mineral density, spinal Z(BMD), and collected blood from 22 members of this kindred. DNA was genotyped on an Applied Biosystems model 377 (ABI PRISM Linkage Mapping Sets; Perkin Elmer Applied Biosystems), by use of fluorescence-based marker sets that included 345 markers. Both two-point and multipoint linkage analyses were performed, by use of affected/unaffected and quantitative-trait models. Spinal Z(BMD) for affected individuals (N = 12) of the kindred was 5.54 +/- 1.40; and for unaffected individuals (N = 16) it was 0.41 +/- 0.81. The trait was present in affected individuals 18-86 years of age, suggesting that HBM influences peak bone mass. The only region of linkage was to a series of markers on chromosome 11 (11q12-13). The highest LOD score (5.21) obtained in two-point analysis, when a quantitative-trait model was used, was at D11S987. Multipoint analysis using a quantitative-trait model confirmed the linkage, with a LOD score of 5.74 near marker D11S987. HBM demonstrates the utility of spinal Z(BMD) as a quantitative bone phenotype that can be used for linkage analysis. Osteoporosis pseudoglioma syndrome also has been mapped to this region of chromosome 11. Identification of the causal gene for both traits will be required for determination of whether a single gene with different alleles that determine a wide range of peak bone densities exists in this region.     

Familial juvenile hyperuricemic nephropathy: localization of the gene on chromosome 16p11.2-and evidence for genetic heterogeneity

Familial juvenile hyperuricemic nephropathy (FJHN), is an autosomal dominant renal disease characterized by juvenile onset of hyperuricemia, gouty arthritis, and progressive renal failure at an early age. Using a genomewide linkage analysis in three Czech affected families, we have identified, on chromosome 16p11.2, a locus for FJHN and have found evidence for genetic heterogeneity and reduced penetrance of the disease. The maximum two-point LOD score calculated with allowance for heterogeneity (HLOD) was 4.70, obtained at recombination fraction 0, with marker D16S3036; multipoint linkage analysis yielded a maximum HLOD score of 4.76 at the same location. Haplotype analysis defined a 10-cM candidate region between flanking markers D16S501 and D16S3113, exhibiting crossover events with the disease locus. The candidate interval contains several genes expressed in the kidney, two of which-uromodulin and NADP-regulated thyroid-hormone-binding protein-represent promising candidates for further analysis.     

Shwachman-Diamond syndrome with exocrine pancreatic dysfunction and bone marrow failure maps to the centromeric region of chromosome 7

Shwachman-Diamond syndrome (SDS) is an autosomal recessive disorder characterized by exocrine pancreatic insufficiency and hematologic and skeletal abnormalities. A genomewide scan of families with SDS was terminated at approximately 50% completion, with the identification of chromosome 7 markers that showed linkage with the disease. Finer mapping revealed significant linkage across a broad interval that included the centromere. The maximum two-point LOD score was 8.7, with D7S473, at a recombination fraction of 0. The maximum multipoint LOD score was 10, in the interval between D7S499 and D7S482 (5.4 cM on the female map and 0 cM on the male map), a region delimited by recombinant events detected in affected children. Evidence from all 15 of the multiplex families analyzed provided support for the linkage, consistent with a single locus for SDS. However, the presence of several different mutations is suggested by the heterogeneity of disease-associated haplotypes in the candidate region.     

Confirmation of linkage of Van der Woude syndrome to chromosome 1q32: evidence of association with STR alleles suggests possible unique origin of the disease mutation

Van der Woude syndrome (VWS) is an autosomal dominant craniofacial disorder with high penetrance and variable expression. Its clinical features are variably expressed, but include cleft lip and/or cleft palate, lip pits and hypodontia. All VWS families studied to date map the disease gene to a < 2 cM region of chromosome 1q32, with no evidence of locus heterogeneity. The aim of this study is to refine the localization of the VWS gene and to further assess possible heterogeneity. We analyzed four multiplex VWS families. All available members were clinically assessed and genotyped for 19 short tandem repeat markers on chromosome 1 in the VWS candidate gene region. We performed two-point and multipoint limit of detection (LOD) score analyses using a high penetrance autosomal dominant model. All families showed positive LOD scores without any recombination in the candidate region. The largest two-point LOD score was 5.87. Our assay method for short tandem repeat (STR) markers provided highly accurate size estimation of marker allele fragment sizes, and therefore enabled us to determine the specific alleles segregating with the VWS gene in each of our four families. We observed a striking pattern of STR allele sharing at several closely linked loci among our four Caucasian VWS families recruited at three different locations in the US. These results suggest the possibility of a unique origin for a mutation responsible for many or most cases of VWS.     

Linkage analysis in the presence of errors III: marker loci and their map as nuisance parameters

In linkage and linkage disequilibrium (LD) analysis of complex multifactorial phenotypes, various types of errors can greatly reduce the chance of successful gene localization. The power of such studies-even in the absence of errors-is quite low, and, accordingly, their robustness to errors can be poor, especially in multipoint analysis. For this reason, it is important to deal with the ramifications of errors up front, as part of the analytical strategy. In this study, errors in the characterization of marker-locus parameters-including allele frequencies, haplotype frequencies (i.e., LD between marker loci), recombination fractions, and locus order-are dealt with through the use of profile likelihoods maximized over such nuisance parameters. It is shown that the common practice of assuming fixed, erroneous values for such parameters can reduce the power and/or increase the probability of obtaining false positive results in a study. The effects of errors in assumed parameter values are generally more severe when a larger number of less informative marker loci, like the highly-touted single nucleotide polymorphisms (SNPs), are analyzed jointly than when fewer but more informative marker loci, such as microsatellites, are used. Rather than fixing inaccurate values for these parameters a priori, we propose to treat them as nuisance parameters through the use of profile likelihoods. It is demonstrated that the power of linkage and/or LD analysis can be increased through application of this technique in situations where parameter values cannot be specified with a high degree of certainty.     

A genomewide linkage scan for quantitative-trait loci for obesity phenotypes

Obesity is an increasingly serious health problem in the world. Body mass index (BMI), percentage fat mass, and body fat mass are important indices of obesity. For a sample of pedigrees that contains >10,000 relative pairs (including 1,249 sib pairs) that are useful for linkage analyses, we performed a whole-genome linkage scan, using 380 microsatellite markers to identify genomic regions that may contain quantitative-trait loci (QTLs) for obesity. Each pedigree was ascertained through a proband who has extremely low bone mass, which translates into a low BMI. A major QTL for BMI was identified on 2q14 near the marker D2S347 with a LOD score of 4.04 in two-point analysis and a maximum LOD score (MLS) of 4.44 in multipoint analysis. The genomic region near 2q14 also achieved an MLS >2.0 for percentage of fat mass and body fat mass. For the putative QTL on 2q14, as much as 28.2% of BMI variation (after adjustment for age and sex) may be attributable to this locus. In addition, several other genomic regions that may contain obesity-related QTLs are suggested. For example, 1p36 near the marker D1S468 may contain a QTL for BMI variation, with a LOD score of 2.75 in two-point analysis and an MLS of 2.09 in multipoint analysis. The genomic regions identified in this and earlier reports are compared for further exploration in extension studies that use larger samples and/or denser markers for confirmation and fine-mapping studies, to eventually identify major functional genes involved in obesity.     

Identification of a locus on chromosome 14q for idiopathic basal ganglia calcification (Fahr disease).

Idiopathic basal ganglia calcification (IBGC) is a neurodegenerative syndrome that is associated with a variety of movement disorders and neurobehavioral and cognitive manifestations. Despite numerous clinical, pathological, and biochemical investigations, its etiology remains unknown. We have identified a multigenerational family with dominantly inherited IBGC and, in 24 members of this family, performed a whole-genome scan using polymorphic microsatellite markers to identify the first chromosomal locus for this disorder (IBGC1). A maximum two-point LOD score of 3.37 was obtained at marker D14S1014, and a maximum multipoint LOD score of 4.95 was obtained between D14S75 and D14S306. The minimal haplotype shared by affected patients extended over a 17.1-cM region bounded by D14S70 and D14S66, which is potentially further narrowed to a 13.3-cM region by a recombination observed in a patient with probable affected status. The age at onset appeared to be decreasing by an average of >20 years with each transmission, which is consistent with genetic anticipation.     

The gene for human fibronectin glomerulopathy maps to 1q32, in the region of the regulation of complement activation gene cluster.

Fibronectin glomerulopathy (GFND) is a newly recognized autosomal dominant disease of the kidney that results in albuminuria, microscopic hematuria, hypertension, renal tubular acidosis type IV, and end-stage renal disease in the 2d to 6th decade of life. The disease is characterized histologically by massive deposits of fibronectin (Fn) present in the subendothelial spaces of renal glomerular capillaries. The cause of human GFND is unknown. In order to localize a candidate gene for GFND, we performed linkage analysis of a large, 193-member pedigree containing 13 affected individuals. Since we had previously excluded the genes for Fn and uteroglobin as candidate genes for GFND, a total-genome search for linkage was performed. Examination of 306 microsatellite markers resulted in a maximum two-point LOD score of 4.17 at a recombination fraction of. 00 for marker D1S249, and a maximum multipoint LOD score of 4.41 for neighboring marker D1S2782. By detection of recombination events, a critical genetic interval of 4.1 cM was identified, which was flanked by markers D1S2872 and D1S2891. These findings confirm that GFND is a distinct disease entity among the fibrillary glomerulopathies. Gene identification will provide insights into the molecular interactions of Fn in GFND, as well as in genetically unaltered conditions.