A web application to perform linkage disequilibrium and linkage analyses on a computational grid

Motivation: Unravelling the genetic architecture of complex traits requires large amounts of data, sophisticated models and large computational resources. The lack of user-friendly software incorporating all these requisites is delaying progress in the analysis of complex traits. Methods: Linkage disequilibrium and linkage analysis (LDLA) is a high-resolution gene mapping approach based on sophisticated mixed linear models, applicable to any population structure. LDLA can use population history information in addition to pedigree and molecular markers to decompose traits into genetic components. Analyses are distributed in parallel over a large public grid of computers in the UK. Results: We have proven the performance of LDLA with analyses of simulated data. There are real gains in statistical power to detect quantitative trait loci when using historical information compared with traditional linkage analysis. Moreover, the use of a grid of computers significantly increases computational speed, hence allowing analyses that would have been prohibitive on a single computer. © The Author 2009. Published by Oxford University Press. All rights reserved.

Migraine association and linkage studies of an endothelial nitric oxide synthase (NOS3) gene polymorphism

Migraine shows strong familial aggregation. However, the number of genes involved in the disorder is unknown and not identified. Nitric oxide is involved in the central processing of pain stimuli and plays an important role in the regulation of basal or stimulated vasodilation. Nitric oxide synthase, which controls the synthesis of nitric oxide, could possibly be a cause, or candidate gene, in migraine etiology. In this study, we detected a polymorphism for endothelial nitric oxide synthase by polymerase chain reaction and tested this for association and linkage to migraine. Results from the study did not show an association of the nitric oxide synthase microsatellite when tested in 91 affected and 85 unaffected individuals. Using the FASTLINK program for parametric linkage analysis, the polymorphism did not show significant linkage to migraine when tested in four migraine pedigrees composed of 116 individuals, 52 affected. Total LOD scores excluded linkage up to 8.5 cM between the nitric oxide synthase polymorphism and migraine. Results using the nonparametric affected pedigree member form of analysis also did not support a role for this gene in migraine etiology.

Regional chromosomal assignment of human renin gene to 1q12→qter and use in linkage studies in Charcot-Marie-Tooth disease

The gene for renin, previously mapped to human chromosome 1, was further localized to 1q12 → qter using human-mouse somatic cell hybrid DNAs. The renin DNA probe used (λ HR5) could detect a HindIII restriction fragment length polymorphism. When used in studies of 12 informative families, no linkage could be found between the renin and Charcot-Marie-Tooth disease. Furthermore, an association of any renin allele with hypertension was not apparent.

Chromosome I linkage studies in Charcot-Marie-Tooth neuropathy type I

Charcot-Marie-Tooth neuropathy type 1 (CMT1) is an autosomal dominant disorder of peripheral nerve. The gene for CMT1 was originally localized to chromosome 1 by linkage to the Duffy blood group, but it has since been shown that not all CMT1 pedigrees show this linkage. We report here the results of linkage studies using five chromosome 1 markers - Duffy (Fy), antithrombin III (AT3), renin (REN), β-nerve growth factor (NGFB), and salivary amylase (AMY1) - in 16 CMT1 pedigrees. The total lod scores exclude close linkage of CMT1 to any of these markers. However, individual families show probable linkage of CMT1 to Duffy, AT3, and/or AMY1. No linkage was indicated with REN or NGFB. These results indicate that possible location of a CMT1 gene between the AMY1 and AT3 loci at p21 and q23, respectively, on chromosome 1 and support the theory that there is at least one other CMT1 gene.

Heterogeneity evidence and linkage studies on Charcot-Marie-Tooth disease

Charcot-Marie-Tooth neuropathy type 1 (CMT1) is an autosomal dominant disorder originally localized to chromosome 1 by linkage to the Duffy blood group. Studies have since shown that the disorder may be heterogeneous, as not all families show this linkage. We tested genetic heterogeneity by the HOMOG computer program in 15 CMT1 pedigrees informative for Duffy. We detected no evidence for heterogeneity in this sample, but when we combined results with previously published lod scores, heterogeneity was statistically significant. Twelve of the 15 families studied did not show linkage to Duffy. We found six of these families to be informative for a chromosome 19 marker, apolipoprotein CII(ApoC2). Despite a previous report showing probable linkage of a non-Duffy-linked CMT1 pedigree to two chromosome 19 markers, we did not detect significant linkage of ApoC2 to CMT1 in these families.

Isolation and use of chromosome 1 probes for linkage studies on Charcot-Marie-Tooth disease

Nine probes were isolated from a human chromosome 1 enriched library and mapped to regions of chromosome 1 using somatic cell hybrid lines. One clone, LR67, which mapped 1q12→q23 detected a BglI RFLP. This probe, as well as 4 other known chromosome 1 markers, α-spectrin, Factor XIIIB, DR10 and DR78, were used for linkage studies in 15 Charcot-Marie-Tooth disease (CMT1) families. Close linking of CMT1 to any of the 5 markers was not indicated. Total lod scores excluded linkage of CMT1 to LR67 and to DR10 at 5 cM or less, to DR78 and 10 cM or less, α-spectrin at 15 cM or less and Factor XIIIB at 20 cM or less. Possible linkage, however, was shown between LR67 and CMT1 at a distance of 30 cM. Also linkage at a distance of 5 cM was detected between this probe and α-spectrin.

Familial typical migraine: Linkage to chromosome 19p13 and evidence for genetic heterogeneity

Migraine is a frequent familial disorder that, in common with most multifactorial disorders, has an unknown etiology. The authors identified several families with multiple individuals affected by typical migraine using a single set of diagnostic criteria and studied these families for cosegregation between the disorder and markers on chromosome 19, the location of a mutation that causes a rare form of familial hemiplegic migraine (FHM). One large tested family showed both cosegregation and significant allele sharing for markers situated within or adjacent to the FHM locus. Multipoint GENEHUNTER results indicated significant excess allele sharing across a 12.6- cM region containing the FHM Ca2+ channel gene, CACNL1A4 (maximum nonparametric linkage Z score = 6.64, p = 0.0026), with a maximum parametric lod score of 1.92 obtained for a (CAG)(n) triplet repeat polymorphism situated in exon 47 of this gene. The CAG expansion did not, however, appear to be the cause of migraine in this pedigree. Other tested families showed neither cosegregation nor excess allele sharing to chromosome 19 markers. HOMOG analysis indicated heterogeneity, generating a maximum HLOD score of 3.6. It was concluded that Chr19 mutations either in the CACNL1A4 gene or a closely linked gene are implicated in some pedigrees with familial typical migraine, and that the disorder is genetically heterogeneous.

A linkage and cross-sectional study of hypertension and obesity using a poly(A) Alu-repeat polymorphism at the glucagon receptor gene locus (17q25)

1. Previous glucagon receptor gene (GCGR) studies have shown a Gly40Ser mutation to be more prevalent in essential hypertension and to affect glucagon binding affinity to its receptor. An Alu-repeat poly(A) polymorphism colocalized to GCGR was used in the present study to test for association and linkage in hypertension as well as association in obesity development. 2. Using a cross-sectional approach, 85 hypertensives and 95 normotensives were genotyped using polymerase chain reaction primers flanking the Alu-repeat. Both hypertensive and normotensive populations were subdivided into lean and obese categories based on body mass index (BMI) to determine involvement of this variant in obesity. For the linkage study, 89 Australian Caucasian hypertension affected sibships (174 sibpairs) were genotyped and the results were analysed using GENE-HUNTER, Mapmaker Sibs, ERPA and SPLINK (all freely available from http://linlkage.rockefeller. edu/soft/list.html). 3. Cross-sectional results for both hypertension and obesity were analysed using Chi-squared and Monte Carlo analyses. Results did not show an association of this variant with either hypertension (χ2 = 6.9, P = 0.14; Monte Carlo χ2 = 7.0, P = 0.11; n = 5000) or obesity (χ2 = 3.3, P = 0.35; Monte Carlo χ2 = 3.26, P = 0.34; n = 5000). In addition, results from the linkage study using hypertensive sib-pairs did not indicate linkage of the poly(A) repent with hypertension. Hence, results did not indicate a role far the Alu-repeat in either hypertension or obesity. However, as the heterozygosity of this poly(A) repeat is low (35%), a larger number of hypertensive sib-pairs may be required to draw definitive conclusions.

Association and linkage analyses of restriction fragment length polymorphisms for the human renin and antithrombin III genes in essential hypertension

Essential hypertension is a highly hereditable disorder in which genetic influences predominate over environmental factors. The molecular genetic profiles which predispose to essential hypertension are not known. In rats with genetic hypertension, there is some recent evidence pointing to linkage of renin gene alleles with blood pressure. The genes for renin and antithrombin III belong to a conserved synteny group which, in humans, spans the q21.3-32.3 region of chromosome I and, in rats, is linkage group X on chromosome 13. The present study examined the association of particular human renin gene (REN) and antithrombin III gene (AT3) polymorphisms with essential hypertension by comparing the frequency of specific alleles for each of these genes in 50 hypertensive offspring of hypertensive parents and 91 normotensive offspring of normotensive parents. In addition, linkage relationships were examined in hypertensive pedigrees with multiple affected individuals. Alleles of a REN HindIII restriction fragment length polymorphism (RFLP) were detected using a genomic clone, λHR5, to probe Southern blots of HindIII-cut leucocyte DNA, and those for an AT3 Pstl RFLP were detected by phATIII 113 complementary DNA probe. The frequencies of each REN allele in the hypertensive group were 0.76 and 0.24 compared with 0.74 and 0.26 in the normotensive group. For AT3, hypertensive allele frequencies were 0.49 and 0.51 compared with normotensive values of 0.54 and 0.46. These differences were not significant by χ2 analysis (P > 0.2). Linkage analysis of a family (data from 16 family members, 10 of whom were hypertensive), informative for both markers, without an age-of-onset correction, and assuming dominant inheritance of hypertension, complete penetrance and a disease frequency of 20%, did not indicate linkage of REN with hypertension, but gave a positive, although not significant, logarithm of the odds for linkage score of 0.784 at a recombination fraction of 0 for AT3 linkage to hypertension. In conclusion, the present study could find no evidence for an association of a REN HindIII RFLP with essential hypertension or for a linkage of the locus defined by this RFLP in a family segregating for hypertension. In the case of an AT3 Pstl RFLP, although association analysis was negative, linkage analysis suggested possible involvement (odds of 6:1 in favour) of a gene located near the 1q23 locus with hypertension in one informative family.

A linkage study across the T cell receptor A and T cell receptor B loci in families with rheumatoid arthritis

Objective. To examine whether the T cell receptor (TCR) A or TCRB loci exhibit linkage with disease in multiplex rheumatoid arthritis (RA) families. Methods. A linkage study was performed in 184 RA families from the UK Arthritis and Rheumatism Council Repository, each containing at least 1 affected sibpair. The microsatellites D14S50, TCRA, and D14S64 spanning the TCRA locus and D7S509, Vβ6.7, and D7S688 spanning the TCRB locus were used as DNA markers. The subjects were genotyped using a semiautomated polymerase chain reaction-based method. Two-point and multipoint linkage analyses were performed. Results. Nonparametric single-marker likelihood odds (LOD) scores were 0.49 (P = 0.07) for D14S50, 0.65 (P = 0.04) for TCRA, 0.07 (P = 0.29) for D14S64, 0.01 (P = 0.43) for D7S509, 0.0 (P = 0.50) for Vβ6.7, and 0.0 (P = 0.50) for D7S688. By multipoint analysis, there was no evidence of linkage at TCRB (LOD score 0), and the maximum LOD score at the TCRA locus was 0.37 (at D14S50). The presence of a susceptibility locus (LOD score < -2.0) was excluded, with lambda ≤ 1.8 at TCRA and ≤1.4 at TCRB. Conclusion. These linkage studies provide no significant evidence of a major germline-encoded TCRA or TCRB component of susceptibility to RA.

Whole-genome linkage analysis of rheumatoid arthritis susceptibility loci in 252 affected sibling pairs in the United Kingdom

Objective. To undertake a systematic wholegenome screen to identify regions exhibiting genetic linkage to rheumatoid arthritis (RA). Methods. Two hundred fifty-two RA-affected sibling pairs from 182 UK families were genotyped using 365 highly informative microsatellite markers. Microsatellite genotyping was performed using fluorescent polymerase chain reaction primers and semiautomated DNA sequencing technology. Linkage analysis was undertaken using MAPMAKER/SIBS for single-point and multipoint analysis. Results. Significant linkage (maximum logarithm of odds score 4.7 [P = 0.000003] at marker D6S276, 1 cM from HLA-DRB1) was identified around the major histocompatibility complex (MHC) region on chromosome 6. Suggestive linkage (P < 7.4 × 10-4) was identified on chromosome 6q by single- and multipoint analysis. Ten other sites of nominal linkage (P < 0.05) were identified on chromosomes 3p, 4q, 7p, 2 regions of 10q, 2 regions of 14q, 16p, 21q, and Xq by single-point analysis and on 3 sites (1q, 14q, and 14q) by multipoint analysis. Conclusion. Linkage to the MHC region was confirmed. Eleven non-HLA regions demonstrated evidence of suggestive or nominal linkage, but none reached the genome-wide threshold for significant linkage (P = 2.2 × 10-5). Results of previous genome screens have suggested that 6 of these regions may be involved in RA susceptibility.

Non-major-histocompatibility-complex genetics of ankylosing spondylitis

There is strong evidence from twin and family studies indicating that a substantial proportion of the heritability of susceptibility to ankylosing spondylitis (AS) and its clinical manifestations is encoded by non-major-histocompatibility-complex genes. Efforts to identify these genes have included genomewide linkage studies and candidate gene association studies. One region, the interleukin (IL)-1 gene complex on chromosome 2, has been repeatedly associated with AS in both Caucasians and Asians. It is likely that more than one gene in this complex is involved in AS, with the strongest evidence to date implicating IL-1A. Identifying the genes underlying other linkage regions has been difficult due to the lack of obvious candidates and the low power of most studies to date to identify genes of the small to moderate magnitude that are likely to be involved. The field is moving towards genomewide association analysis, involving much larger datasets of unrelated cases and controls. Early successes using this approach in other diseases indicates that it is likely to identify genes in common diseases like AS, but there remains the risk that the common-variant, common-disease hypothesis will not hold true in AS. Nonetheless, it is appropriate for the field to be cautiously optimistic that the next few years will bring great advances in our understanding of the genetics of this condition.

The genetics of osteoporosis: Future diagnostic possibilities

The risk of developing osteoporosis is determined by the interaction of several mostly unknown genes and environmental factors. Genetic studies in osteoporosis have largely focussed on association studies of a small number of candidate genes, with few linkage studies performed, and large areas of the genome remaining unexplored. Identifying the genes involved in osteoporosis would be a major advance in our understanding of the causation of the disease, and lead to advances in diagnosis, risk prediction, and potentially preventive and therapeutic measures.

Combined analysis of three whole genome linkage scans for Ankylosing Spondylitis

Objective. Ankylosing spondylitis (AS) is a debilitating chronic inflammatory condition with a high degree of familiality (λs=82) and heritability (>90%) that primarily affects spinal and sacroiliac joints. Whole genome scans for linkage to AS phenotypes have been conducted, although results have been inconsistent between studies and all have had modest sample sizes. One potential solution to these issues is to combine data from multiple studies in a retrospective meta-analysis. Methods: The International Genetics of Ankylosing Spondylitis Consortium combined data from three whole genome linkage scans for AS (n=3744 subjects) to determine chromosomal markers that show evidence of linkage with disease. Linkage markers typed in different centres were integrated into a consensus map to facilitate effective data pooling. We performed a weighted meta-analysis to combine the linkage results, and compared them with the three individual scans and a combined pooled scan. Results: In addition to the expected region surrounding the HLA-B27 gene on chromosome 6, we determined that several marker regions showed significant evidence of linkage with disease status. Regions on chromosome 10q and 16q achieved 'suggestive' evidence of linkage, and regions on chromosomes 1q, 3q, 5q, 6q, 9q, 17q and 19q showed at least nominal linkage in two or more scans and in the weighted meta-analysis. Regions previously associated with AS on chromosome 2q (the IL-1 gene cluster) and 22q (CYP2D6) exhibited nominal linkage in the meta-analysis, providing further statistical support for their involvement in susceptibility to AS. Conclusion: These findings provide a useful guide for future studies aiming to identify the genes involved in this highly heritable condition. . Published by on behalf of the British Society for Rheumatology.

Suggestive linkage of the parathyroid receptor type 1 to osteoporosis

We have investigated the role of 23 candidate genes in the control of bone mineral density (BMD) by linkage studies in families of probands with osteoporosis (lumbar spine [LS] or femoral neck [FN] BMD T score < -2.5) and low BMD relative to an age- and gender-matched cohort (Z score < -2.0). One hundred and fifteen probands (35 male, 80 female) and 499 of their first- or second-degree relatives (223 males and 276 females) were recruited for the study. BMD was measured at the LS and FN using dual-energy X-ray absorptiometry and expressed as age- and gender-matched Z scores corrected for body mass index. The candidate genes studied were the androgen receptor, type I collagen A1 (COLIA1), COLIA2, COLIIA1, vitamin D receptor (VDR), colony-stimulating factor 1, calcium-sensing receptor, epidermal growth factor (EGF), estrogen receptor 1 (ESR1), fibrillin type 1, insulin-like growth factor 1, interleukin-1 alpha (IL-1α), interleukin-4 (IL-4), interleukin-6 (IL-6), interleukin-11 (IL-11), osteopontin, parathyroid hormone (PTH), PTH-related peptide, PTH receptor type 1 (PTHR1), transforming growth factor-beta 1, and tumor necrosis factors alpha and beta. Sixty-four microsatellites lying close to or within these genes were investigated for linkage with BMD. Using the program MapMaker/Sibs there was suggestive evidence of linkage between BMD and PTHR1 (maximum LOD score obtained [MLS] 2.7-3.5). Moderate evidence of linkage was also observed with EGF (MLS 1.8), COLIA1 (MLS 1.7), COLIIA1/VDR (MLS 1.7), ESR1 (MLS 1.4), IL-1α (MLS 1.4), IL-4 (MLS 1.2), and IL-6 (MLS 1.2). Variance components analysis using the program ACT, correcting for proband-wise ascertainment, also showed evidence of linkage (p ≤0.05) at markers close to or within the candidate genes IL- 1α, PTHR1, IL-6, and COLIIA1/VDR. Further studies will be required to confirm these findings, to refine the location of gene responsible for the observed linkage, and to screen the candidate genes targeted at these loci for mutations.