Model-free linkage analysis with covariates confirms linkage of prostate cancer to chromosomes 1 and 4

As with many complex genetic diseases, genome scans for prostate cancer have given conflicting results, often failing to provide replication of previous findings. One factor contributing to the lack of consistency across studies is locus heterogeneity, which can weaken or even eliminate evidence for linkage that is present only in a subset of families. Currently, most analyses either fail to account for locus heterogeneity or attempt to account for it only by partitioning data sets into smaller and smaller portions. In the present study, we model locus heterogeneity among affected sib pairs with prostate cancer by including covariates in the linkage analysis that serve as surrogate measures of between-family linkage differences. The model is a modification of the Olson conditional logistic model for affected relative pairs. By including Gleason score, age at onset, male-to-male transmission, and/or number of affected first-degree family members as covariates, we detected linkage near three locations that were previously identified by linkage (1q24-25 [HPC1; LOD score 3.25, P=.00012], 1q42.2-43 [PCAP; LOD score 2.84, P=.0030], and 4q [LOD score 2.80, P=.00038]), near the androgen-receptor locus on Xq12-13 (AR; LOD score 3.06, P=.00053), and at five new locations (LOD score > 2.5). Without covariates, only a few weak-to-moderate linkage signals were found, none of which replicate findings of previous genome scans. We conclude that covariate-based linkage analysis greatly improves the likelihood that linked regions will be found by incorporation of information about heterogeneity within the sample.     

Genome scan of human systemic lupus erythematosus by regression modeling: evidence of linkage and epistasis at 4p16-15.2

Systemic lupus erythematosus (SLE) is a complex autoimmune disorder involving at least hormonal, environmental, and genetic factors. Familial aggregation, a 2%-3% sibling recurrence rate, monozygotic twin concordance >20%, association with several candidate genes, as well as the results of five genome scans support a genetic component. We present here the results of a genome scan of 126 pedigrees multiplex for SLE, including 469 sibling pairs (affected and unaffected) and 175 affected relative pairs. Using the revised multipoint Haseman-Elston regression technique for concordant and discordant sibling pairs and a conditional logistic regression technique for affected relative pairs, we identify a novel linkage to chromosome 4p16-15.2 (P=.0003 and LOD=3.84) and present evidence of an epistatic interaction between chromosome 4p16-15.2 and chromosome 5p15 in our European American families. We confirm the evidence of linkage to chromosome 4p16-15.2 in European American families using data from an independent pedigree collection. In addition, our data support the published results of three independent studies for nine purportedly linked regions and agree with the previously published results from a subset of these data for three regions. In summary, results from two new analytical techniques establish and confirm linkage with SLE at 4p16-15.2, indicate epistasis between 4p16-15.2 and 5p15, and confirm other linkage effects with SLE that have been reported elsewhere.     

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.     

Genetic linkage to chromosome 22q12 for a heavy-smoking quantitative trait in two independent samples

We conducted a genomewide linkage screen of a simple heavy-smoking quantitative trait, the maximum number of cigarettes smoked in a 24-h period, using two independent samples: 289 Australian and 155 Finnish nuclear multiplex families, all of which were of European ancestry and were targeted for DNA analysis by use of probands with a heavy-smoking phenotype. We analyzed the trait, using a regression of identity-by-descent allele sharing on the sum and difference of the trait values for relative pairs. Suggestive linkage was detected on chromosome 22 at 27-29 cM in each sample, with a LOD score of 5.98 at 26.96 cM in the combined sample. After additional markers were used to localize the signal, the LOD score was 5.21 at 25.46 cM. To assess the statistical significance of the LOD score in the combined sample, 1,000 simulated genomewide screens were conducted, resulting in an empirical P value of .006 for the LOD score of 5.21. This linkage signal is driven mainly by the microsatellite marker D22S315 (22.59 cM), which had a single-point LOD score of 5.41 in the combined sample and an empirical P value <.001 from 1,000 simulated genomewide screens. This marker is located within an intron of the gene ADRBK2, encoding the beta-adrenergic receptor kinase 2. Fine mapping of this linkage region may reveal variants contributing to heaviness of smoking, which will lead to a better understanding of the genetic mechanisms underlying nicotine dependence.     

Linkage of a gene for familial hypobetalipoproteinemia to chromosome 3p21.1-22

Familial hypobetalipoproteinemia (FHBL) is an apparently autosomal dominant disorder of lipid metabolism characterized by less than fifth percentile age- and sex-specific levels of apolipoprotein beta (apobeta) and low-density lipoprotein-cholesterol. In a minority of cases, FHBL is due to truncation-producing mutations in the apobeta gene on chromosome 2p23-24. Previously, we reported on a four-generation FHBL kindred in which we had ruled out linkage of the trait to the apobeta gene. To locate other loci containing genes for low apobeta levels in the kindred, a genomewide search was conducted. Regions on 3p21.1-22 with two-point LOD scores >1.5 were identified. Additional markers were typed in the region of these signals. Two-point LOD scores in the region of D3S2407 increased to 3.35 at O = 0. GENEHUNTER confirmed this finding with an nonparametric multipoint LOD score of 7.5 (P=.0004). Additional model-free analyses were conducted with the square root of the apobeta level as the phenotype. Results from the Loki and SOLAR programs further confirmed linkage of FHBL to 3p21.1-22. Weaker linkage to a region near D19S916 was also indicated by Loki and SOLAR. Thus, a heretofore unidentified genetic susceptibility locus for FHBL may reside on chromosome 3.     

Performance comparison of two-point linkage methods using microsatellite markers flanking known disease locations

The Genetic Analysis Workshop 14 simulated data presents an interesting, challenging, and plausible example of a complex disease interaction in a dataset. This paper summarizes the ease of detection for each of the simulated Kofendrerd Personality Disorder (KPD) genes across all of the replicates for five standard linkage statistics. Using the KPD affection status, we have analyzed the microsatellite markers flanking each of the disease genes, plus an additional 2 markers that were not linked to any of the disease loci. All markers were analyzed using the following two-point linkage methods: 1) a MMLS, which is a standard admixture LOD score maximized over theta, alpha, and mode of inheritance, 2) a MLS calculated by GENEHUNTER, 3) the Kong and Cox LOD score as computed by MERLIN, 4) a MOD score (standard heterogeneity LOD maximized over theta, alpha, and a grid of genetic model parameters), and 5) the PPL, a Bayesian statistic that directly measures the strength of evidence for linkage to a marker. All of the major loci (D1-D4) were detectable with varying probabilities in the different populations. However, the modifier genes (D5 and D6) were difficult to detect, with similar distributions under the null and alternative across populations and statistics. The pooling of the four datasets in each replicate (n = 350 pedigrees) greatly improved the chance of detecting the major genes using all five methods, but failed to increase the chance to detect D5 and D6.     

Comprehensive multi-stage linkage analyses identify a locus for adult height on chromosome 3p in a healthy Caucasian population

There have been a number of genome-wide linkage studies for adult height in recent years. These studies have yielded few well-replicated loci, and none have been further confirmed by the identification of associated gene variants. The inconsistent results may be attributable to the fact that few studies have combined accurate phenotype measures with informative statistical modelling in healthy populations. We have performed a multi-stage genome-wide linkage analysis for height in 275 adult sibling pairs drawn randomly from the Victorian Family Heart Study (VFHS), a healthy population-based Caucasian cohort. Height was carefully measured in a standardised fashion on regularly calibrated equipment. Following genome-wide identification of a peak Z-score of 3.14 on chromosome 3 at 69 cM, we performed a fine-mapping analysis of this region in an extended sample of 392 two-generation families. We used a number of variance components models that incorporated assortative mating and shared environment effects, and we observed a peak LOD score of approximately 3.5 at 78 cM in four of the five models tested. We also demonstrated that the most prevalent model in the literature gave the worst fit, and the lowest LOD score (2.9) demonstrating the importance of appropriate modelling. The region identified in this study replicates the results of other genome-wide scans of height and bone-related phenotypes, strongly suggesting the presence of a gene important in bone growth on chromosome 3p. Association analyses of relevant candidate genes should identify the genetic variants responsible for the chromosome 3p linkage signal in our population.     

Genomewide linkage scan identifies a novel susceptibility locus for restless legs syndrome on chromosome 9p

Restless legs syndrome (RLS) is a common neurological disorder that affects 5%-12% of all whites. To genetically dissect this complex disease, we characterized 15 large and extended multiplex pedigrees, consisting of 453 subjects (134 affected with RLS). A familial aggregation analysis was performed, and SAGE FCOR was used to quantify the total genetic contribution in these families. A weighted average correlation of 0.17 between first-degree relatives was obtained, and heritability was estimated to be 0.60 for all types of relative pairs, indicating that RLS is a highly heritable trait in this ascertained cohort. A genomewide linkage scan, which involved >400 10-cM-spaced markers and spanned the entire human genome, was then performed for 144 individuals in the cohort. Model-free linkage analysis identified one novel significant RLS-susceptibility locus on chromosome 9p24-22 with a multipoint nonparametric linkage (NPL) score of 3.22. Suggestive evidence of linkage was found on chromosome 3q26.31 (NPL score 2.03), chromosome 4q31.21 (NPL score 2.28), chromosome 5p13.3 (NPL score 2.68), and chromosome 6p22.3 (NPL score 2.06). Model-based linkage analysis, with the assumption of an autosomal-dominant mode of inheritance, validated the 9p24-22 linkage to RLS in two families (two-point LOD score of 3.77; multipoint LOD score of 3.91). Further fine mapping confirmed the linkage result and defined this novel RLS disease locus to a critical interval. This study establishes RLS as a highly heritable trait, identifies a novel genetic locus for RLS, and will facilitate further cloning and identification of the genes for RLS.     

A genomewide scan for attention-deficit/hyperactivity disorder in an extended sample: suggestive linkage on 17p11

Attention-deficit/hyperactivity disorder (ADHD [MIM 143465]) is a common, highly heritable neurobehavioral disorder of childhood onset, characterized by hyperactivity, impulsivity, and/or inattention. As part of an ongoing study of the genetic etiology of ADHD, we have performed a genomewide linkage scan in 204 nuclear families comprising 853 individuals and 270 affected sibling pairs (ASPs). Previously, we reported genomewide linkage analysis of a "first wave" of these families composed of 126 ASPs. A follow-up investigation of one region on 16p yielded significant linkage in an extended sample. The current study extends the original sample of 126 ASPs to 270 ASPs and provides linkage analyses of the entire sample, using polymorphic microsatellite markers that define an approximately 10-cM map across the genome. Maximum LOD score (MLS) analysis identified suggestive linkage for 17p11 (MLS=2.98) and four nominal regions with MLS values >1.0, including 5p13, 6q14, 11q25, and 20q13. These data, taken together with the fine mapping on 16p13, suggest two regions as highly likely to harbor risk genes for ADHD: 16p13 and 17p11. Interestingly, both regions, as well as 5p13, have been highlighted in genomewide scans for autism.     

Intracranial aneurysms in Finnish families: confirmation of linkage and refinement of the interval to chromosome 19q13.3

We recently reported a two-stage genomewide screen of 48 sib pairs affected with intracranial aneurysms (IAs) that revealed suggestive linkage to chromosome 19q13, with a LOD score of 2.58. The region supporting linkage spanned ~22 cM. Here, we report a follow-up study of the locus at 19q13, with a sample size expanded to 139 affected sib pairs, along with 83 other affected relative pairs (222 affected relative pairs in total). Suggestive linkage was observed in both independent sample sets, and linkage was significant in the combined set at 70 cM (LOD score 3.50; P=.00006) and at 80 cM (LOD score 3.93; P=.00002). Linkage was highly significant at 70 cM (LOD score 5.70; P=.000001) and at 80 cM (LOD score 3.99; P=.00005) when a covariate measuring the number of affected individuals in the nuclear family was included. To evaluate further the contribution to the linkage signal from families with more than two affected relatives, we performed model-based linkage analysis with a recessive model and a range of penetrances, and we obtained maximum linkage at 70 cM (LOD score 3.16; P=.00007) with a penetrance of 0.3. We then estimated location by using GENEFINDER. The most likely location for a gene predisposing to IAs in the Finnish population is in a region with a 95% confidence interval of 11.6 cM (P=.00007) centered 2.0 cM proximal to D19S246.     

Linkage of Familial Schizophrenia to Chromosome 13q32

Over the past 4 years, a number of investigators have reported findings suggestive of linkage to schizophrenia, with markers on chromosomes 13q32 and 8p21, with one recent study by Blouin et al. reporting significant linkage to these regions. As part of an ongoing genome scan, we evaluated microsatellite markers spanning chromosomes 8 and 13, for linkage to schizophrenia, in 21 extended Canadian families. Families were analyzed under autosomal dominant and recessive models, with broad and narrow definitions of schizophrenia. All models produced positive LOD scores with markers on 13q, with higher scores under the recessive models. The maximum three-point LOD scores were obtained under the recessive-broad model: 3.92 at recombination fraction (theta).1 with D13S793, under homogeneity, and 4.42 with alpha=.65 and straight theta=0 with D13S793, under heterogeneity. Positive LOD scores were also obtained, under all models, for markers on 8p. Although a maximum two-point LOD score of 3.49 was obtained under the dominant-narrow model with D8S136 at straight theta=0.1, multipoint analysis with closely flanking markers reduced the maximum LOD score in this region to 2. 13. These results provide independent significant evidence of linkage of a schizophrenia-susceptibility locus to markers on 13q32 and support the presence of a second susceptibility locus on 8p21.     

Linkage analysis of extremely discordant and concordant sibling pairs identifies quantitative trait loci influencing variation in human menopausal age

Age at natural menopause may be used as parameter for evaluating the rate of ovarian aging. Environmental factors determine only a small part of the large variation in menopausal age. Studies have shown that genetic factors are likely to be involved in variation in menopausal age. To identify quantitative-trait loci for this trait, we performed a genomewide linkage study with age at natural menopause as a continuous quantitative phenotype in Dutch sister pairs, through use of a selective sampling scheme. A total of 165 families were ascertained using extreme selected sampling and were genotyped for 417 markers. Data were analyzed by Haseman-Elston regression and by an adjusted variance-components analysis. Subgroup analyses for early and late menopausal age were conducted by Haseman-Elston regression. In the adjusted variance-components analysis, 12 chromosomes had a LOD score of > or =1.0. Two chromosomal regions showed suggestive linkage: 9q21.3 (LOD score 2.6) and Xp21.3 (LOD score 3.1). Haseman-Elston regression showed rather similar locations of the peaks but yielded lower LOD scores. A permutation test to obtain empirical P values resulted in a significant peak on the X chromosome. To our knowledge, this is the first study to attempt to identify loci responsible for variability in menopausal age and in which several chromosomal regions were identified with suggestive and significant linkage. Although the finding of the region on the X chromosome comes as no surprise, because of its widespread involvement in premature ovarian failure, the definition of which particular gene is involved is of great interest. The region on chromosome 9 deserves further consideration. Both findings require independent confirmation.     

Genome screen to identify susceptibility genes for Parkinson disease in a sample without parkin mutations

Parkinson disease (PD) is a common neurodegenerative disorder characterized by bradykinesia, resting tremor, muscular rigidity, and postural instability, as well as by a clinically significant response to treatment with levodopa. Mutations in the alpha-synuclein gene have been found to result in autosomal dominant PD, and mutations in the parkin gene produce autosomal recessive juvenile-onset PD. We have studied 203 sibling pairs with PD who were evaluated by a rigorous neurological assessment based on (a) inclusion criteria consisting of clinical features highly associated with autopsy-confirmed PD and (b) exclusion criteria highly associated with other, non-PD pathological diagnoses. Families with positive LOD scores for a marker in an intron of the parkin gene were prioritized for parkin-gene testing, and mutations in the parkin gene were identified in 22 families. To reduce genetic heterogeneity, these families were not included in subsequent genome-screen analysis. Thus, a total of 160 multiplex families without evidence of a parkin mutation were used in multipoint nonparametric linkage analysis to identify PD-susceptibility genes. Two models of PD affection status were considered: model I included only those individuals with a more stringent diagnosis of verified PD (96 sibling pairs from 90 families), whereas model II included all examined individuals as affected, regardless of their final diagnostic classification (170 sibling pairs from 160 families). Under model I, the highest LOD scores were observed on chromosome X (LOD score 2.1) and on chromosome 2 (LOD score 1.9). Analyses performed with all available sibling pairs (model II) found even greater evidence of linkage to chromosome X (LOD score 2.7) and to chromosome 2 (LOD score 2.5). Evidence of linkage was also found to chromosomes 4, 5, and 13 (LOD scores >1.5). Our findings are consistent with those of other linkage studies that have reported linkage to chromosomes 5 and X.     

Chromosome-12 Mapping of Late-Onset Alzheimer Disease among Caribbean Hispanics

Linkage to chromosome 12p for familial Alzheimer disease (AD) has been inconsistent. Using 35 markers near the centromere of chromosome 12, we investigated 79 Caribbean Hispanic families with AD. Two-point linkage analysis using affected sib pairs yielded LOD scores of 3.15 at D12S1623 and 1.43 at D12S1042. The LOD score at D12S1623 decreased to 1.62 in families with late-onset (age >65 years) AD (LOAD), but the LOD score at D12S1042 was unchanged. Among families negative for the apolipoprotein E (APOE-epsilon 4) allele, the LOD score for D12S1623 was lower (1.01), whereas that for D12S1042 increased to 1.73. Among families positive for the APOE-epsilon 4 allele, none of the LOD scores reached 1. Multipoint affected-relative-pair analysis showed peaks at D12S1623 (nonparametric linkage [NPL] score 1.52; P=.028) and near D12S1042, at D12S1057 (NPL score 1.57; P=.027). NPL scores for both D12S1623 and D12S1057 increased in families affected with LOAD, but, in APOE-epsilon 4-negative families, only scores for the region flanking D12S1623 remained elevated (NPL score 1.74; P=.013). This study of Caribbean Hispanics with familial AD extends and provides modest evidence of linkage to loci on chromosome 12p. Linkage varied by age at onset of AD and by the presence or absence of the APOE-epsilon 4 allele.     

A genomewide scan for intelligence identifies quantitative trait loci on 2q and 6p

Between 40% and 80% of the variation in human intelligence (IQ) is attributable to genetic factors. Except for many rare mutations resulting in severe cognitive dysfunction, attempts to identify these factors have not been successful. We report a genomewide linkage scan involving 634 sibling pairs designed to identify chromosomal regions that explain variation in IQ. Model-free multipoint linkage analysis revealed evidence of a significant quantitative-trait locus for performance IQ at 2q24.1-31.1 (LOD score 4.42), which overlaps the 2q21-33 region that has repeatedly shown linkage to autism. A second region revealed suggestive linkage for both full-scale and verbal IQs on 6p25.3-22.3 (LOD score 3.20 for full-scale IQ and 2.33 for verbal IQ), overlapping marginally with the 6p22.3-21.31 region implicated in reading disability and dyslexia.     

Comparison of genome screens for two independent cohorts provides replication of suggestive linkage of bone mineral density to 3p21 and 1p36

Low bone mineral density (BMD) is a major risk factor for osteoporotic fracture. Studies of BMD in families and twins have shown that this trait is under strong genetic control. To identify regions of the genome that contain quantitative trait loci (QTL) for BMD, we performed independent genomewide screens, using two complementary study designs. We analyzed unselected nonidentical twin pairs (1,094 pedigrees) and highly selected, extremely discordant or concordant (EDAC) sib pairs (254 pedigrees). Nonparametric multipoint linkage (NPL) analyses were undertaken for lumbar spine and total-hip BMD in both cohorts and for whole-body BMD in the unselected twin pairs. The maximum evidence of linkage in the unselected twins (spine BMD, LOD 2.7) and the EDAC pedigrees (spine BMD, LOD 2.1) was observed at chromosome 3p21 (76 cM and 69 cM, respectively). These combined data indicate the presence, in this region, of a gene that regulates BMD. Furthermore, evidence of linkage in the twin cohort (whole-body BMD; LOD 2.4) at chromosome 1p36 (17 cM) supports previous findings of suggestive linkage to BMD in the region. Weaker evidence of linkage (LOD 1.0-2.3) in either cohort, but not both, indicates the locality of additional QTLs. These studies validate the use, in linkage analysis, of large cohorts of unselected twins phenotyped for multiple traits, and they highlight the importance of conducting genome scans in replicate populations as a prelude to positional cloning and gene discovery.     

Fine mapping of the chromosome 12 late-onset Alzheimer disease locus: potential genetic and phenotypic heterogeneity

Apolipoprotein E (APOE) is the only confirmed susceptibility gene for late-onset Alzheimer disease (AD). In a recent genomic screen of 54 families with late-onset AD, we detected significant evidence for a second late-onset AD locus located on chromosome 12 between D12S373 and D12S390. Linkage to this region was strongest in 27 large families with at least one affected individual without an APOE-4 allele, suggesting that APOE and the chromosome 12 locus might have independent effects. We have since genotyped several additional markers across the region, to refine the linkage results. In analyzing these additional data, we have addressed the issue of heterogeneity in the data set by weighting results by clinical and neuropathologic features, sibship size, and APOE genotype. When considering all possible affected sib pairs (ASPs) per nuclear family, we obtained a peak maximum LOD score between D12S1057 and D12S1042. The magnitude and location of the maximum LOD score changed when different weighting schemes were used to control for the number of ASPs contributed by each nuclear family. Using the affected-relative-pair method implemented in GENEHUNTER-PLUS, we obtained a maximum LOD score between D12S398 and D12S1632, 25 cM from the original maximum LOD score. These results indicate that family size influences the location estimate for the chromosome 12 AD gene. The results of conditional linkage analysis by use of GENEHUNTER-PLUS indicated that evidence for linkage to chromosome 12 was stronger in families with affected individuals lacking an APOE-4 allele; much of this evidence came from families with affected individuals with neuropathologic diagnosis of dementia with Lewy bodies (DLB). Taken together, these results indicate that the chromosome 12 locus acts independently of APOE to increase the risk of late-onset familial AD and that it may be associated with the DLB variant of AD.     

Comparison of microsatellites versus single-nucleotide polymorphisms in a genome linkage screen for prostate cancer-susceptibility Loci

Prostate cancer is one of the most common cancers among men and has long been recognized to occur in familial clusters. Brothers and sons of affected men have a 2-3-fold increased risk of developing prostate cancer. However, identification of genetic susceptibility loci for prostate cancer has been extremely difficult. Although the suggestion of linkage has been reported for many chromosomes, the most promising regions have been difficult to replicate. In this study, we compare genome linkage scans using microsatellites with those using single-nucleotide polymorphisms (SNPs), performed in 467 men with prostate cancer from 167 families. For the microsatellites, the ABI Prism Linkage Mapping Set version 2, with 402 microsatellite markers, was used, and, for the SNPs, the Early Access Affymetrix Mapping 10K array was used. Our results show that the presence of linkage disequilibrium (LD) among SNPs can lead to inflated LOD scores, and this seems to be an artifact due to the assumption of linkage equilibrium that is required by the current genetic-linkage software. After excluding SNPs with high LD, we found a number of new LOD-score peaks with values of at least 2.0 that were not found by the microsatellite markers: chromosome 8, with a maximum model-free LOD score of 2.2; chromosome 2, with a LOD score of 2.1; chromosome 6, with a LOD score of 4.2; and chromosome 12, with a LOD score of 3.9. The LOD scores for chromosomes 6 and 12 are difficult to interpret, because they occurred only at the extreme ends of the chromosomes. The greatest gain provided by the SNP markers was a large increase in the linkage information content, with an average information content of 61% for the SNPs, versus an average of 41% for the microsatellite markers. The strengths and weaknesses of microsatellite versus SNP markers are illustrated by the results of our genome linkage scans.     

Direct power comparisons between simple LOD scores and NPL scores for linkage analysis in complex diseases.

Several methods have been proposed for linkage analysis of complex traits with unknown mode of inheritance. These methods include the LOD score maximized over disease models (MMLS) and the "nonparametric" linkage (NPL) statistic. In previous work, we evaluated the increase of type I error when maximizing over two or more genetic models, and we compared the power of MMLS to detect linkage, in a number of complex modes of inheritance, with analysis assuming the true model. In the present study, we compare MMLS and NPL directly. We simulated 100 data sets with 20 families each, using 26 generating models: (1) 4 intermediate models (penetrance of heterozygote between that of the two homozygotes); (2) 6 two-locus additive models; and (3) 16 two-locus heterogeneity models (admixture alpha = 1.0,.7,.5, and.3; alpha = 1.0 replicates simple Mendelian models). For LOD scores, we assumed dominant and recessive inheritance with 50% penetrance. We took the higher of the two maximum LOD scores and subtracted 0.3 to correct for multiple tests (MMLS-C). We compared expected maximum LOD scores and power, using MMLS-C and NPL as well as the true model. Since NPL uses only the affected family members, we also performed an affecteds-only analysis using MMLS-C. The MMLS-C was both uniformly more powerful than NPL for most cases we examined, except when linkage information was low, and close to the results for the true model under locus heterogeneity. We still found better power for the MMLS-C compared with NPL in affecteds-only analysis. The results show that use of two simple modes of inheritance at a fixed penetrance can have more power than NPL when the trait mode of inheritance is complex and when there is heterogeneity in the data set.     

Modeling the effect of an associated single-nucleotide polymorphism in linkage studies

For linkage analysis in affected sibling pairs, we propose a regression model to incorporate information from a disease-associated single-nucleotide polymorphism located under the linkage peak. This model can be used to study if the associated single-nucleotide polymorphism marker partly explains the original linkage peak. Two sources of information are used for performing this task, namely the genotypes of the parents and the genotypes of the siblings. We applied the methods to three significantly disease-associated single-nucleotide polymorphisms and five microsatellite markers at the end of chromosome 3 of replicate 1 of Aipotu population. Two out of five of the microsatellite markers showed a LOD score higher than 3. The question to be answered was whether one of the single-nucleotide polymorphisms partly explains these high LOD scores. We did not have the answers when we analyzed the data.