The effect of three-monthly albendazole treatment on Th2 responses: Differential effects on IgE and IL-5.

Helminth parasites induce a strong Th2 response, characterized by high levels of IgE and elevated signature cytokines such as IL-5. As many global deworming programmes are underway, there is concern that this might lead to emergence of Th1-mediated pathologies when the counterbalancing helminth-induced Th2 response is absent. Therefore, we assessed the effect of deworming on Th2-mediated responses in a household-clustered randomized controlled trial in Indonesia. Total plasma IgE and whole-blood IL-5 responses to mitogen phytohaemagglutinin (PHA) were measured in 1494 and 682 subjects, respectively, at baseline, 9 and 21 months after three-monthly single-dose treatment with albendazole or placebo. Anthelmintic treatment did not result in complete removal of helminth infections in the community. However, treatment significantly decreased IgE levels in albendazole- compared to placebo-treated subjects. IL-5 responses to PHA were not significantly affected by anthelmintic treatment and tended to increase in albendazole-treated subjects, indicating that intensive treatment of helminth parasites has different outcomes on B-cell (IgE levels) and T-cell (IL-5) responses. The data shows that 2 years of deworming can have differential effects on responses typified as Th2-mediated, which needs to be taken into account when examining the impact of helminths on noncommunicable diseases.

Maternal and child cytokine relationship in early life is not altered by cytokine gene polymorphisms.

The development of immune responses is influenced by the interaction between environmental and genetic factors. Our previous study showed a close association between maternal and young infant's cytokine responses. The question is how this association evolves over time and the contribution of genetic polymorphisms to this association. Five cytokines in mitogen-stimulated whole blood culture were measured from pregnant mothers and their children aged 2, 5, 12, 24 and 48 months. Cytokine gene polymorphisms were determined in both mothers and children. High production of maternal interleukin (IL)-10, tumour necrosis factor-α (TNF-α) and interferon-γ (IFN-γ) was significantly associated with higher levels of the corresponding cytokines in their children at 2 months (T2), but the association decreased over time. Maternal single-nucleotide polymorphism (SNP) in IFN-γ gene, rs3181032, was found to be associated with child's IFN-γ levels at T2 only, whereas maternal IL-10 rs4579758 and child's TNF-α rs13215091 were associated with child's corresponding cytokines at later ages but not at T2. In the final models including the gene polymorphisms, maternal cytokines were still the strongest determinant of child cytokines. Maternal cytokine during pregnancy, which could be a proxy for child's environmental factors, showed its highest impact at early age, with no or little influence from genetic factors.

CD226 (DNAM-1) is associated with susceptibility to juvenile idiopathic arthritis.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is considered a complex genetic autoimmune disease. We investigated the association of genetic variants previously implicated in JIA, autoimmunity and/or immunoregulation, with susceptibility to JIA. METHODS: A genetic association study was performed in 639 JIA patients and 1613 healthy controls of northwest European descent. Ninety-three single nucleotide polymorphisms (SNP) were genotyped in a candidate gene approach. Results of the entire JIA patient group (all subtypes) were compared with results obtained, alternatively, with a clinically homogeneous patient group including only oligoarticular and rheumatoid factor (RF) negative polyarticular JIA patients (n=493). Meta-analyses were performed for all SNPs that have been typed in other Caucasian JIA cohorts before. RESULTS: SNPs in or near PTPN22, VTCN1, the IL2-IL21 region, ANKRD55 and TNFA were confirmed to be associated with JIA (p<0.05), strengthening the evidence for involvement of these genes in JIA. In the majority of these replicated SNPs, effect sizes were larger when analysing a homogeneous patient cohort than when analysing all subtypes. We identified two novel associations with oligoarticular and RF-negative polyarticular JIA: CD226 rs763361 (OR 1.30, 95% CI 1.12 to 1.51, p=0.0006) and CD28 rs1980422 (OR 1.29, 95% CI 1.07 to 1.55, p=0.008). Meta-analyses including reported studies confirmed the association of both SNPs with susceptibility to JIA (OR 1.16, p=0.001 and OR 1.18, p=0.001, for rs763361 and rs1980422, respectively). CONCLUSIONS: The CD226 gene has been identified as novel association with JIA, and a SNP near CD28 as a suggestive association. Both genes are probable candidate risk factors, since they are involved in costimulation of T cells.

Mortality in osteoarthritis patients.

OBJECTIVES: To investigate whether all-cause mortality and deaths due to cardiovascular disease are increased in patients who have consulted primary or secondary health care with symptoms and signs of osteoarthritis (OA). METHOD: This study included 383 patients with symptomatic OA at multiple sites from the Genetics ARthrosis and Progression (GARP) study (mean age 60 years, 82% women, 3693 person-years of follow-up) and 459 patients with primary hand, knee, or hip OA from the Osteoarthritis Care Clinic (OCC) study (mean age 61 years, 88% women, 1890 person-years of follow-up). Standardized mortality ratios (SMRs) with 95% confidence intervals (CIs) were calculated for all-cause mortality and causes of deaths in comparison to the general population. Cox proportional hazard ratios (HRs) with 95% CIs were used to associate baseline characteristics with all-cause mortality. RESULTS: In the GARP study, 26 patients died whereas 48 deaths were expected (SMR 0.54, 95% CI 0.37-0.79). The SMR was 0.47 (95% CI 0.29-0.76) in women and 0.73 (95% CI 0.39-1.35) in men. Similar results were found in the OCC study (SMR 0.45, 95% CI 0.25-0.82). Malignancy and cardiovascular disease were the main causes of deaths in GARP. Male sex (HR 3.04, 95% CI 1.38-6.69), increasing age (HR 1.10, 95% CI 1.05-1.16), and self-reported cancer (HR 8.29, 95% CI 3.12-22.03) were associated with increased mortality in GARP. CONCLUSIONS: Patients consulting health care for their OA are not at higher risk of death than the general population. These results suggest that the management of OA patients may not need to focus specifically on the treatment of cardiovascular risk factors and comorbidities.

A genetic variant in the region of MMP-9 is associated with serum levels and progression of joint damage in rheumatoid arthritis.

OBJECTIVES: The severity of joint destruction is highly variable between rheumatoid arthritis (RA) patients. The majority of its heritability is still unexplained. Several autoimmune diseases share genetic risk variants that may also influence disease progression. We aimed to identify genetic risk factors for the severity of joint damage in RA by studying genetic susceptibility loci of several autoimmune diseases. METHODS: In phase 1, 3143 sets of x-rays of 646 Dutch RA patients taken over 7 years (Sharp van der Heijde (SHS) scored) were studied. Genotyping was done by Immunochip. Associations of single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) >0.01 and joint destruction were analysed. In phase 2, 686 North American RA patients with 926 SHS-scored x-rays over 15 years of follow-up were evaluated. In both phases multiple testing corrections were done for the number of uncorrelated SNPs; the thresholds for significance were p<1.1×10(-6) and p<0.0036. Matrix metalloproteinase 9 (MMP-9) levels were measured with ELISA in baseline serum samples. RESULTS: In phase 1, 109 SNPs associated significantly with joint destruction (p<1.1×10(-6)). Of these, 76 were located in the HLA region; the 33 non-HLA variants were studied in phase 2. Here two variants were associated with the severity of joint destruction: rs451066 on chromosome 14 (p=0.002, MAF=0.20) and rs11908352 on chromosome 20 (p=0.002, MAF=0.21). Rs11908352 is located near the gene encoding MMP-9. Serum levels of MMP-9 were significantly associated with the rs11908352 genotypes (p=0.007). CONCLUSIONS: These data indicate that two loci that confer risk to other autoimmune diseases also affect the severity of joint destruction in RA. Rs11908352 may influence joint destruction via MMP-9 production.

Genetic variation in VTCN1 (B7-H4) is associated with course of disease in juvenile idiopathic arthritis.

OBJECTIVE: The course of disease in juvenile idiopathic arthritis (JIA) is unpredictable with episodes of activity and remission. In order to identify predictive factors, 93 SNPs, JIA subtype, age at onset and ANA status were studied in relation to disease course. METHODS: Genetic and clinical parameters were analysed in a cohort of 272 Caucasian patients with persistent oligoarthritis (n=129), extended oligoarthritis (n=57) and rheumatoid factor negative polyarthritis (n=86). Categories of disease course (remitting (n=65), intermediate (n=96) and unremitting (n=111)) were designed based on the cumulative time spent in active disease in the first 2 years. RESULTS: Univariate analysis revealed association of the course of disease with JIA subtype (p=5.7*10(-5)) and three SNPs; VTCN1 rs10 923 223 (p=4.4*10(-5)), VTCN1 rs12 046 117 (p=0.017) and CDK6 rs42 041 (p=0.038). In a subsequent multivariate ordinal logistic regression analysis, VTCN1 rs10 923 223 (OR 0.41, 95%-CI 0.26 to 0.63) and JIA subtype (OR 3.8, 95%-CI 2.0 to 7.2; OR 2.5, 95%-CI 1.4 to 4.2, for extended oligoarthritis and RF-negative polyarthritis vs persistent oligoarthritis, respectively) were the strongest independent factors for course of disease. CONCLUSIONS: This study provides evidence that VTCN1, encoding B7-H4, is associated with course of disease in selected subtypes of JIA. VTCN1 might be useful in predicting the course of disease.

Sex-specific effects of naturally occurring variants in the dopamine receptor D2 locus on insulin secretion and type 2 diabetes susceptibility.

AIMS: Modulation of dopamine receptor D2 (DRD2) activity affects insulin secretion in both rodents and isolated pancreatic ß-cells. We hypothesized that single nucleotide polymorphisms in the DRD2/ANKK1 locus may affect susceptibility to type 2 diabetes in humans. METHODS: Four potentially functional variants in the coding region of the DRD2/ANKK1 locus (rs1079597, rs6275, rs6277, rs1800497) were genotyped and analysed for type 2 diabetes susceptibility in up to 25 000 people (8148 with type 2 diabetes and 17687 control subjects) from two large independent Dutch cohorts and one Danish cohort. In addition, 340 Dutch subjects underwent a 2-h hyperglycaemic clamp to investigate insulin secretion. Since sexual dimorphic associations related to DRD2 polymorphisms have been previously reported, we also performed a gender-stratified analysis. RESULTS: rs1800497 at the DRD2/ANKK1 locus was associated with a significantly increased risk for type 2 diabetes in women (odds ratio 1.14 (1.06-1.23); P = 4.1*104) but not in men (odds ratio 1.00 (95% CI 0.93-1.07); P = 0.92) or the combined group. Although rs1800497 was not associated with insulin secretion, we did find another single nucleotide polymorphism in this locus, rs6275, to be associated with increased first-phase glucose-stimulated insulin secretion in women (P = 5.5*104) but again not in men (P = 0.34). CONCLUSION: The present data identify DRD2/ANKK1 as a potential sex-specific type 2 diabetes susceptibility gene.

A genetic variant in granzyme B is associated with progression of joint destruction in rheumatoid arthritis.

OBJECTIVE: Genetic factors account for an estimated 45-58% of the variance in joint destruction in rheumatoid arthritis (RA). The serine proteinase granzyme B induces target cell apoptosis, and several in vitro studies suggest that granzyme B is involved in apoptosis of chondrocytes. Serum levels of granzyme B are increased in RA and are also associated with radiographic erosions. The aim of this study was to investigate GZMB as a candidate gene accounting for the severity of joint destruction in RA. METHODS: A total of 1,418 patients with 4,885 radiograph sets of the hands and feet from 4 independent cohorts were studied. First, explorative analyses were performed in 600 RA patients in the Leiden Early Arthritis Clinic cohort. Fifteen single-nucleotide polymorphisms (SNPs) tagging GZMB were tested. Significantly associated SNPs were genotyped in data sets representing patients from the Groningen, Sheffield, and Lund cohorts. In each data set, the relative increase in the annual rate of progression in the presence of a genotype was assessed. Data were summarized in a meta-analysis. The association of GZMB with the RNA expression level of the GZMB genomic region was tested by mapping expression quantitative trait loci (QTLs) on 1,469 whole blood samples. RESULTS: SNP rs8192916 was significantly associated with the rate of joint destruction in the first cohort and in the meta-analysis of all data sets. Patients homozygous for the minor allele of rs8192916 had a higher rate of joint destruction per year compared with other patients (P = 7.8 × 10(-4)). Expression QTL of GZMB identified higher expression in the presence of the minor allele of rs8192916 (P = 2.27 × 10(-5)). CONCLUSION: SNP rs8192916 located in GZMB is associated with the progression of joint destruction in RA as well as with RNA expression in whole blood.

Genetic variants in the region of the C1q genes are associated with rheumatoid arthritis.

Rodent models for arthritis implicate a role for complement in disease development and progression. In humans, complement deposition has been observed in inflamed synovia of rheumatoid arthritis (RA) patients. In this study we analysed whether genetic variants of complement component C1q predispose to RA. We genotyped single nucleotide polymorphisms (SNPs) in and around the C1q genes, C1qA, C1qB and C1qC, in a Dutch set of 845 RA cases and 1046 controls. Replication was sought in a sample set from North America (868 cases/1193 controls), and a meta-analysis was performed in a combined samples set of 8000 cases and 23 262 controls of European descent. We determined C1q serum levels in relation to C1q genotypes. In the discovery phase, five of the 13 SNPs tested in the C1q genes showed a significant association with RA. Additional analysis of the genomic area around the C1q genes revealed that the strongest associating SNPs were confined to the C1q locus. Within the C1q locus we observed no additional signal independent of the strongest associating SNP, rs292001 [odds ratio (OR) = 0·72 (0·58-0·88), P = 0·0006]. The variants of this SNP were associated with different C1q serum levels in healthy controls (P = 0·006). Interestingly, this SNP was also associated significantly in genome-wide association studies (GWAS) from the North American Rheumatoid Arthritis Consortium study, confirming the association with RA [OR = 0·83 (0·69-1·00), P = 0·043]. Combined analysis, including integrated data from six GWAS studies, provides support for the genetic association. Genetic variants in C1q are correlated with C1q levels and may be a risk for the development of RA.

Genetic studies on components of the Wnt signalling pathway and the severity of joint destruction in rheumatoid arthritis.

BACKGROUND: Progression of joint destruction in rheumatoid arthritis (RA) is partly heritable; knowledge of genetic factors may increase our understanding of the mechanisms underlying joint destruction. The activity of the Wnt/ß-catenin pathway influences osteoblast differentiation. Dickkopf-1 (Dkk-1) and sclerostin (Sost) are negative regulators and lipoprotein receptor-related protein-5 (LRP-5) and Kremen-1 are transmembrane receptors involved in this pathway. OBJECTIVE: To study variants in the genes encoding these proteins in relation to progression of joint destruction. METHODS: 1418 patients with RA of four cohorts with 4885 sets of hands and feet x-rays were studied. Explorative analyses were performed on 600 patients with RA from Leiden on single nucleotide polymorphisms (SNPs) tagging Dkk-1, Sost, Kremen-1 and LRP-5. SNPs significantly associating with joint damage progression were subsequently genotyped in cohorts from Groningen (NL), Sheffield (UK) and Lund (Sweden). Data were summarised in meta-analyses. Serum levels of functional Dkk-1 and sclerostin were measured and studied in relation to genotypes. RESULTS: In the first cohort, six Dkk-1, three Sost, one Kremen-1 and 10 LRP-5 SNPs were significantly associated with radiological progression of joint destruction. Three Dkk-1 SNPs were associated significantly with progression of joint damage in the meta-analysis, also after correction for multiple testing (rs1896368, rs1896367 and rs1528873). Two Sost SNPs tended to significance (rs4792909 and rs6503475, p=0.07 after false discovery rate correction). Gene-gene interactions between SNPs on Dkk-1 and Sost were seen. Serum levels of Dkk-1 were significantly correlated with the genotypes in rs1896368 (p=0.02). CONCLUSIONS: Patients with RA carrying risk alleles of genetic variants in Dkk-1 have higher serum levels of functional Dkk-1 and more progressive joint destruction over time.

Comparison of methodologies for analysing the progression of joint destruction in rheumatoid arthritis.

OBJECTIVES: Progression of joint destruction is an important phenotypic feature in rheumatoid arthritis (RA). When factors have small effect sizes, both the avoidance of phenotypic misclassification and discerning true effects from noise are challenging. Assembling radiological measurements repeatedly in time harbours a smaller risk of misclassification than single measurements. Given serial measurements, different methods of analysis can be applied. This study evaluates different statistical methods of analysing longitudinal data. METHODS: Three statistical methods were studied: linear regression (LR), generalized estimating equations (GEE), and multivariate normal regression analysis (MRA). All were applied longitudinally, testing for differences in radiological progression rates. As genetic variants are known to have small effect sizes, two genetic variants were studied as examples: rs675520 (located in the TNFAIP3-OLIG3 region) and the presence of the human leucocyte antigen (HLA) shared epitope (SE) alleles. Radiological data for 602 early RA patients with yearly radiographs and 7-years of follow-up were used. The powers obtained with the methods and the robustness against missingness were evaluated as outcome measures. RESULTS: The presence of the rs675520 polymorphism and the HLA-SE risk genotype was associated with a 0.65-0.77 and 1.17-1.51 fold increased rate of joint destruction, respectively. The analyses performed with MRA resulted in smaller 95% confidence intervals (CIs) than the analyses using LR or GEE. In addition, the 95% CIs increased with the number of radiographs per patient. The power of MRA was higher than that of GEE. MRA was more robust against selective missingness than GEE or LR with a two-step approach (LR(ts)). CONCLUSIONS: A multivariate normal regression model on subsequent radiographs is a powerful and robust method for analysing longitudinal joint destruction data.

Genetic variants in IL15 associate with progression of joint destruction in rheumatoid arthritis: a multicohort study.

BACKGROUND: Interleukin (IL)-15 levels are increased in serum, synovium and bone marrow of patients with rheumatoid arthritis (RA). IL-15 influences both the innate and the adaptive immune response; its major role is activation and proliferation of T cells. There are also emerging data that IL-15 affects osteoclastogenesis. The authors investigated the association of genetic variants in IL15 with the rate of joint destruction in RA. METHOD: 1418 patients with 4885 x-ray sets of both hands and feet of four independent data sets were studied. First, explorative analyses were performed on 600 patients with early RA enrolled in the Leiden Early Arthritis Clinic. Twenty-five single-nucleotide polymorphisms (SNPs) tagging IL-15 were tested. Second, SNPs with significant associations in the explorative phase were genotyped in data sets from Groningen, Sheffield and Lund. In each data set, the relative increase of the progression rate per year in the presence of a genotype was assessed. Subsequently, data were summarised in an inverse weighting meta-analysis. RESULTS: Five SNPs were significantly associated with rate of joint destruction in phase 1 and typed in the other data sets. Patients homozygous for rs7667746, rs7665842, rs2322182, rs6821171 and rs4371699 had respectively 0.94-, 1.04-, 1.09-, 1.09- and 1.09-fold rate of joint destruction compared to other patients (p=4.0×10(-6), p=3.8×10(-4), p=5.0×10(-3), p=5.0×10(-3) and p=9.4×10(-3)). DISCUSSION: Independent replication was not obtained, possibly due to insufficient power. Meta-analyses of all data sets combined resulted in significant results for four SNPs (rs7667746, p<0.001; rs7665842, p<0.001; rs4371699, p=0.01; rs6821171, p=0.01). These SNPs were also significant after correction for multiple testing. CONCLUSION: Genetic variants in IL-15 are associated with progression of joint destruction in RA.

Osteoarthritis susceptibility genes influence the association between hip morphology and osteoarthritis.

OBJECTIVE: The identified osteoarthritis (OA) susceptibility genes are mainly active in skeletal development and could thus affect joint geometry. Because nonoptimal joint geometry is a risk factor for the development of OA, we investigated if and how the path that leads from nonoptimal joint geometry to OA of the hip is influenced by these genes. METHODS: The shape of the hips of subjects in the Genetics, Osteoarthritis and Progression Study, consisting of sibling pairs with symptomatic OA at multiple joint locations, was quantified by applying a statistical shape model to radiographs. Shape aspects (modes) were correlated to OA characteristics. We then tested for the association of shape modes with OA susceptibility single-nucleotide polymorphisms (SNPs) of GDF5, FRZB, and DIO2. RESULTS: Four of 23 shape modes (mode 1, mode 17, mode 18, and mode 21) were strongly associated with OA characteristics. We observed a significant interaction between carrier status of DIO2 rs12885300 and hip OA characteristics for mode 1 (P = 0.005). This indicates that this specific aspect of hip shape correlates with OA characteristics only in carriers of the susceptibility allele. CONCLUSION: Our results suggest that it is more likely that the rs12885300 SNP of DIO2 increases the vulnerability of cartilage to nonoptimal bone shapes rather than directly influencing the formation of these shapes.

Non-HLA genes modulate the risk of rheumatoid arthritis associated with HLA-DRB1 in a susceptible North American Native population.

Most of the genetic risk for rheumatoid arthritis (RA) is conferred by 'shared epitope' (SE), encoding alleles of HLA-DRB1. Specific North American Native (NAN) populations have RA prevalence rates of 2-5%, representing some of the highest rates estimated worldwide. As many NAN populations also demonstrate a high background frequency of SE, we sought to determine whether other genetic factors contribute to disease risk in this predisposed population. RA patients (n=333) and controls (n=490) from the Cree/Ojibway NAN population in Central Canada were HLA-DRB1 typed and tested for 21 single-nucleotide polymorphisms (SNPs) that have previously been associated with RA, including PTPN22, TRAF1-C5, CTLA4, PADI4, STAT4, FCRL3, CCL21, MMEL1-TNFRSF14, CDK6, PRKCQ, KIF5A-PIP4K2C, IL2RB, TNFAIP3, IL10-1082G/A and REL. Our findings indicate that SE is prevalent and represents a major genetic risk factor for RA in this population (82% cases versus 68% controls, odds ratio=2.2, 95% confidence interval 1.6-3.1, P<0.001). We also demonstrate that in the presence of SE, the minor allele of MMEL1-TNFRSF14 significantly reduces RA risk in a dominant manner, whereas TRAF1-C5 increases the risk. These findings point to the importance of non-HLA genes in determining RA risk in a population with a high frequency of disease predisposing HLA-DRB1 alleles.

The major risk alleles of age-related macular degeneration (AMD) in CFH do not play a major role in rheumatoid arthritis (RA).

Because activation of the alternative pathway (AP) of the complement system is an important aspect of both age-related macular degeneration (AMD) and rheumatoid arthritis (RA), we wished to address the question whether genetic risk factors of the AP inhibitor complement factor H (CFH) for AMD would also be risk factors for RA. For this purpose we genotyped single nucleotide polymorphisms (SNPs) in a Dutch set of RA patients and controls. Similarly, a meta-analysis using a Spanish cohort of RA as well as six large genome-wide association studies (GWAS) studies was performed. For these SNPs we analysed more than 6000 patients and 20,000 controls. The CFH variants, I62V, Y402H, IVS1 and IVS10, known to associate strongly with AMD, did not show a significant association with the risk of developing RA despite a strong statistical power to detect such differences. In conclusion, the major risk alleles of AMD in CFH do not have a similar effect on developing RA.

A mixture model with random-effects components for classifying sibling pairs.

In healthy aging research, typically multiple health outcomes are measured, representing health status. The aim of this paper was to develop a model-based clustering approach to identify homogeneous sibling pairs according to their health status. Model-based clustering approaches will be considered on the basis of linear mixed effect model for the mixture components. Class memberships of siblings within pairs are allowed to be correlated, and within a class the correlation between siblings is modeled using random sibling pair effects. We propose an expectation-maximization algorithm for maximum likelihood estimation. Model performance is evaluated via simulations in terms of estimating the correct parameters, degree of agreement, and the ability to detect the correct number of clusters. The performance of our model is compared with the performance of standard model-based clustering approaches. The methods are used to classify sibling pairs from the Leiden Longevity Study according to their health status. Our results suggest that homogeneous healthy sibling pairs are associated with a longer life span. Software is available for fitting the new models.

Robust age at onset linkage analysis in nuclear families.

OBJECTIVE: Standard methods for linkage analysis ignore the phenotype of the parents when they are not genotyped. However, this information can be useful for gene mapping. In this paper we propose methods for age at onset genetic linkage analysis in sibling pairs, taking into account parental age at onset. METHODS: Two new score statistics are derived, one from an additive gamma frailty model and one from a log-normal frailty model. The score statistics are classical non-parametric linkage (NPL) statistics weighted by a function of the age at onset of the four family members. The weight depends on information from registries (age-specific incidences) and family studies (sib-sib and father-mother correlation). RESULTS: In order to investigate how age at onset of sibs and their parents affect the information for linkage analysis the weight functions were studied for rare and common disease models, realistic models for breast cancer and human lifespan. We studied the performance of the weighted NPL methods by simulations. As illustration, the score statistics were applied to the GAW12 data. The results show that it is useful to include parental age at onset information in genetic linkage analysis.

Allele-sharing statistics using information on family history.

When conducting genetic studies for complex traits, large samples are commonly required to detect new genetic factors. A possible strategy to decrease the sample size is to reduce heterogeneity using available information. In this paper we propose a new class of model-free linkage analysis statistics which takes into account the information given by the ungenotyped affected relatives (positive family history). This information is included into the scoring function of classical allele-sharing statistics. We studied pedigrees of affected sibling pairs with one ungenotyped affected relative. We show that, for rare allele common complex diseases, the proposed method increases the expected power to detect linkage. Allele-sharing methods were applied to the symptomatic osteoarthritis GARP study where taking into account the family-history increased considerably the evidence of linkage in the region of the DIO2 susceptibility locus.

Recurrence risk due to germ line mosaicism: Duchenne and Becker muscular dystrophy.

The presence of multiple affected offspring from apparently non-carrier parents is caused by germ line mosaicism. Although germ line mosaicism has been reported for many diseases, figures for recurrence risks are known for only a few of them. In X-linked Duchenne and Becker muscular dystrophies (DMD/BMD), the recurrence risk for non-carrier females due to germ line mosaicism has been estimated to be between 14% and 20% (95% confidence interval 3-30) if the risk haplotype is transmitted. In this study, we have analyzed 318 DMD/BMD cases in which the detected mutation was de novo with the aim of obtaining a better estimate of the 'true' number of germ line mosaics and a more precise recurrence risk. This knowledge is essential for genetic counseling. Our data indicate a recurrence risk of 8.6% (4.8-12.2) if the risk haplotype is transmitted, but there is a remarkable difference between proximal (15.6%) (4.1-27.0) and distal (6.4%) (2.1-10.6) deletions. Overall, most mutations originated in the female. Deletions occur more often on the X chromosome of the maternal grandmother, whereas point mutations occur on the X chromosome of the maternal grandfather. In unhaplotyped de novo DMD/BMD families, the risk of recurrence of the mutation is 4.3%.

Score test for age at onset genetic linkage analysis in selected sibling pairs.

A new score statistic is derived, which uses information from registries (age-specific incidences) and family studies (sib-sib marginal correlation) to weight affected sibling pairs according to their age at onset. Age at onset of sibling pairs is modelled by a gamma frailty model. From this model we derive a bivariate survival function, which depends on the marginal survival and on the marginal correlation. The score statistic for linkage is a classical nonparametric linkage (NPL) statistic where the identical by descent sharing is weighted by a particular function of the age at onset data. Since the statistic is based on survival models, it can also be applied to discordant and healthy sibling pairs. Simulation studies show that the proposed method is robust and more powerful than standard NPL methods. As illustration we apply the new score statistic to data from a breast cancer study.