Purpose

DNA repair deficiencies have been postulated to play a role in the development and progression of cardiovascular disease (CVD). The hypothesis is that DNA damage accumulating with age may induce cell death, which promotes formation of unstable plaques. Defects in DNA repair mechanisms may therefore increase the risk of CVD events. We examined whether the joints effect of common genetic variants in 5 DNA repair pathways may influence the risk of CVD events.

Methods

The PLINK set-based test was used to examine the association to myocardial infarction (MI) of the DNA repair pathway in GWAS data of 866 subjects of the GENetic DEterminants of Restenosis (GENDER) study and 5,244 subjects of the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER) study. We included the main DNA repair pathways (base excision repair, nucleotide excision repair, mismatch repair, homologous recombination and non-homologous end-joining (NHEJ)) in the analysis.

Results

The NHEJ pathway was associated with the occurrence of MI in both GENDER (P = 0.0083) and PROSPER (P = 0.014). This association was mainly driven by genetic variation in the MRE11A gene (PGENDER = 0.0001 and PPROSPER = 0.002). The homologous recombination pathway was associated with MI in GENDER only (P = 0.011), for the other pathways no associations were observed.

Conclusion

This is the first study analyzing the joint effect of common genetic variation in DNA repair pathways and the risk of CVD events, demonstrating an association between the NHEJ pathway and MI in 2 different cohorts.

Thyroid hormone is essential for normal metabolism and development, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over their life span. In addition, even mild alterations in thyroid function are associated with weight changes, atrial fibrillation, osteoporosis, and psychiatric disorders. To identify novel variants underlying thyroid function, we performed a large meta-analysis of genome-wide association studies for serum levels of the highly heritable thyroid function markers TSH and FT4, in up to 26,420 and 17,520 euthyroid subjects, respectively. Here we report 26 independent associations, including several novel loci for TSH (PDE10A, VEGFA, IGFBP5, NFIA, SOX9, PRDM11, FGF7, INSR, ABO, MIR1179, NRG1, MBIP, ITPK1, SASH1, GLIS3) and FT4 (LHX3, FOXE1, AADAT, NETO1/FBXO15, LPCAT2/CAPNS2). Notably, only limited overlap was detected between TSH and FT4 associated signals, in spite of the feedback regulation of their circulating levels by the hypothalamic-pituitary-thyroid axis. Five of the reported loci (PDE8B, PDE10A, MAF/LOC440389, NETO1/FBXO15, and LPCAT2/CAPNS2) show strong gender-specific differences, which offer clues for the known sexual dimorphism in thyroid function and related pathologies. Importantly, the TSH-associated loci contribute not only to variation within the normal range, but also to TSH values outside the reference range, suggesting that they may be involved in thyroid dysfunction. Overall, our findings explain, respectively, 5.64% and 2.30% of total TSH and FT4 trait variance, and they improve the current knowledge of the regulation of hypothalamic-pituitary-thyroid axis function and the consequences of genetic variation for hypo- or hyperthyroidism.

Author Summary

Levels of thyroid hormones are tightly regulated by TSH produced in the pituitary, and even mild alterations in their concentrations are strong indicators of thyroid pathologies, which are very common worldwide. To identify common genetic variants associated with the highly heritable markers of thyroid function, TSH and FT4, we conducted a meta-analysis of genome-wide association studies in 26,420 and 17,520 individuals, respectively, of European ancestry with normal thyroid function. Our analysis identified 26 independent genetic variants regulating these traits, several of which are new, and confirmed previously detected polymorphisms affecting TSH (within the PDE8B gene and near CAPZB, MAF/LOC440389, and NR3C2) and FT4 (within DIO1) levels. Gender-specific differences in the genetic effects of several variants for TSH and FT4 levels were identified at several loci, which offer clues to understand the known sexual dimorphism in thyroid function and pathology. Of particular clinical interest, we show that TSH-associated loci contribute not only to normal variation, but also to TSH values outside reference range, suggesting that they may be involved in thyroid dysfunction. Overall, our findings add to the developing landscape of the regulation of thyroid homeostasis and the consequences of genetic variation for thyroid related diseases.

Glycosylation of immunoglobulin G (IgG) influences IgG effector function by modulating binding to Fc receptors. To identify genetic loci associated with IgG glycosylation, we quantitated N-linked IgG glycans using two approaches. After isolating IgG from human plasma, we performed 77 quantitative measurements of N-glycosylation using ultra-performance liquid chromatography (UPLC) in 2,247 individuals from four European discovery populations. In parallel, we measured IgG N-glycans using MALDI-TOF mass spectrometry (MS) in a replication cohort of 1,848 Europeans. Meta-analysis of genome-wide association study (GWAS) results identified 9 genome-wide significant loci (P<2.27×10−9) in the discovery analysis and two of the same loci (B4GALT1 and MGAT3) in the replication cohort. Four loci contained genes encoding glycosyltransferases (ST6GAL1, B4GALT1, FUT8, and MGAT3), while the remaining 5 contained genes that have not been previously implicated in protein glycosylation (IKZF1, IL6ST-ANKRD55, ABCF2-SMARCD3, SUV420H1, and SMARCB1-DERL3). However, most of them have been strongly associated with autoimmune and inflammatory conditions (e.g., systemic lupus erythematosus, rheumatoid arthritis, ulcerative colitis, Crohn's disease, diabetes type 1, multiple sclerosis, Graves' disease, celiac disease, nodular sclerosis) and/or haematological cancers (acute lymphoblastic leukaemia, Hodgkin lymphoma, and multiple myeloma). Follow-up functional experiments in haplodeficient Ikzf1 knock-out mice showed the same general pattern of changes in IgG glycosylation as identified in the meta-analysis. As IKZF1 was associated with multiple IgG N-glycan traits, we explored biomarker potential of affected N-glycans in 101 cases with SLE and 183 matched controls and demonstrated substantial discriminative power in a ROC-curve analysis (area under the curve = 0.842). Our study shows that it is possible to identify new loci that control glycosylation of a single plasma protein using GWAS. The results may also provide an explanation for the reported pleiotropy and antagonistic effects of loci involved in autoimmune diseases and haematological cancer.

Author Summary

After analysing glycans attached to human immunoglobulin G in 4,095 individuals, we performed the first genome-wide association study (GWAS) of the glycome of an individual protein. Nine genetic loci were found to associate with glycans with genome-wide significance. Of these, four were enzymes that directly participate in IgG glycosylation, thus the observed associations were biologically founded. The remaining five genetic loci were not previously implicated in protein glycosylation, but the most of them have been reported to be relevant for autoimmune and inflammatory conditions and/or haematological cancers. A particularly interesting gene, IKZF1 was found to be associated with multiple IgG N-glycans. This gene has been implicated in numerous diseases, including systemic lupus erythematosus (SLE). We analysed N-glycans in 101 cases with SLE and 183 matched controls and demonstrated their substantial biomarker potential. Our study shows that it is possible to identify new loci that control glycosylation of a single plasma protein using GWAS. Our results may also provide an explanation for opposite effects of some genes in autoimmune diseases and haematological cancer.

Background

Coronary restenosis after percutaneous coronary intervention still remains a significant problem, despite all medical advances. Unraveling the mechanisms leading to restenosis development remains challenging. Many studies have identified genetic markers associated with restenosis, but consistent replication of the reported markers is scarce. The aim of the current study was to analyze the joined effect of previously in literature reported candidate genes for restenosis in the GENetic DEterminants of Restenosis (GENDER) databank.

Methodology/Principal Findings

Candidate genes were selected using a MEDLINE search including the terms ‘genetic polymorphism’ and ‘coronary restenosis’. The final set included 36 genes. Subsequently, all single nucleotide polymorphisms (SNPs) in the genomic region of these genes were analyzed in GENDER using set-based analysis in PLINK. The GENDER databank contains genotypic data of 2,571,586 SNPs of 295 cases with restenosis and 571 matched controls. The set, including all 36 literature reported genes, was, indeed, significantly associated with restenosis, p = 0.024 in the GENDER study. Subsequent analyses of the individual genes demonstrated that the observed association of the complete set was determined by 6 of the 36 genes.

Conclusion

Despite overt inconsistencies in literature, with regard to individual candidate gene studies, this is the first study demonstrating that the joint effect of all these genes together, indeed, is associated with restenosis.

Summary

In order to study family based association in the presence of linkage we extend a generalized linear mixed model proposed for genetic linkage analysis (Lebrec and van Houwelingen, 2007) by adding a genotypic effect to the mean. The corresponding score test is a weighted FBAT statistic, where the weight depends on the linkage effect and on other genetic and shared environmental effects. For testing of genetic association in the presence of gene covariate interaction, we propose a linear regression method where the family-specific score statistic is regressed on family-specific covariates. Both statistics are straightforward to compute. Simulation results show that adjusting the weight for the within-family variance structure may be a powerful approach in the presence of environmental effects. The test statistic for genetic association in the presence of gene-covariate interaction improved the power for detecting association. For illustration we analyze the Rheumatoid Arthritis data from GAW15. Adjusting for smoking and anti-CCP increased the significance of the association with the DR locus.

Mechanisms underlying the variation in human life expectancy are largely unknown, but lipid metabolism and especially lipoprotein size was suggested to play an important role in longevity. We have performed comprehensive lipid phenotyping in the Leiden Longevity Study (LLS). By applying multiple logistic regression analysis we tested for the first time the effects of parameters in lipid metabolism (i.e., classical serum lipids, lipoprotein particle sizes, and apolipoprotein E levels) on longevity independent of each other. Parameters in lipid metabolism were measured in offspring of nonagenarian siblings from 421 families of the LLS (n = 1,664; mean age, 59 years) and in the partners of the offspring as population controls (n = 711; mean age, 60 years). In the initial model, where lipoprotein particles sizes, classical serum lipids and apolipoprotein E were included, offspring had larger low-density lipoprotein (LDL) particle sizes (p = 0.017), and lower triglyceride levels (p = 0.026), indicating that they displayed a more beneficial lipid profile. After backwards regression only LDL size (p = 0.014) and triglyceride levels (p = 0.05) were associated with offspring from long-lived families. Sex-specific backwards regression analysis revealed that LDL particle sizes were associated with male longevity (increase in log odds ratio (OR) per unit = 0.21; p = 0.023). Triglyceride levels (decrease OR per unit = 0.22; p = 0.01), but not LDL particle size, were associated with female longevity. Due to the analysis of a comprehensive lipid profile, we confirmed an important role of lipid metabolism in human longevity, with LDL size and triglyceride levels as major predicting factors.

Genotype imputation has become an essential tool in the analysis of genome-wide association scans. This technique allows investigators to test association at ungenotyped genetic markers, and to combine results across studies that rely on different genotyping platforms. In addition, imputation is used within long-running studies to reuse genotypes produced across generations of platforms. Typically, genotypes of controls are reused and cases are genotyped on more novel platforms yielding a case–control study that is not matched for genotyping platforms. In this study, we scrutinize such a situation and validate GWAS results by actually retyping top-ranking SNPs with the Sequenom MassArray platform. We discuss the needed quality controls (QCs). In doing so, we report a considerable discrepancy between the results from imputed and retyped data when applying recommended QCs from the literature. These discrepancies appear to be caused by extrapolating differences between arrays by the process of imputation. To avoid false positive results, we recommend that more stringent QCs should be applied. We also advocate reporting the imputation quality measure (RT2) for the post-imputation QCs in publications.

The MUTYH gene is involved in base excision repair. MUTYH mutations predispose to recessively inherited colorectal polyposis and cancer. Here, we evaluate an association with breast cancer (BC), following up our previous finding of an elevated BC frequency among Dutch bi-allelic MUTYH mutation carriers. A case–control study was performed comparing 1,469 incident BC patients (ORIGO cohort), 471 individuals displaying features suggesting a genetic predisposition for BC, but without a detectable BRCA1 or BRCA2 mutation (BRCAx cohort), and 1,666 controls. First, for 303 consecutive patients diagnosed before age 55 years and/or with multiple primary breast tumors, the MUTYH coding region and flanking introns were sequenced. The remaining subjects were genotyped for five coding variants, p.Tyr179Cys, p.Arg309Cys, p.Gly396Asp, p.Pro405Leu, and p.Ser515Phe, and four tagging SNPs, c.37-2487G>T, p.Val22Met, c.504+35G>A, and p.Gln338His. No bi-allelic pathogenic MUTYH mutations were identified. The pathogenic variant p.Gly396Asp and the variant of uncertain significance p.Arg309Cys occurred twice as frequently in BRCAx subjects as compared to incident BC patients and controls (p = 0.13 and p = 0.15, respectively). The likely benign variant p.Val22Met occurred less frequently in patients from the incident BC (p = 0.03) and BRCAx groups (p = 0.11), respectively, as compared to the controls. Minor allele genotypes of several MUTYH variants showed trends towards association with lobular BC histology. This extensive case–control study could not confirm previously reported associations of MUTYH variants with BC, although it was too small to exclude subtle effects on BC susceptibility.

Electronic supplementary material

The online version of this article (doi:10.1007/s10549-012-1965-0) contains supplementary material, which is available to authorized users.

The Leiden Longevity Study consists of families that express extended survival across generations, decreased morbidity in middle-age, and beneficial metabolic profiles. To identify which pathways drive this complex phenotype of familial longevity and healthy aging, we performed a genome-wide gene expression study within this cohort to screen for mRNAs whose expression changes with age and associates with longevity. We first compared gene expression profiles from whole blood samples between 50 nonagenarians and 50 middle-aged controls, resulting in identification of 2,953 probes that associated with age. Next, we determined which of these probes associated with longevity by comparing the offspring of the nonagenarians (50 subjects) and the middle-aged controls. The expression of 360 probes was found to change differentially with age in members of the long-lived families. In a RT-qPCR replication experiment utilizing 312 controls, 332 offspring and 79 nonagenarians, we confirmed a nonagenarian specific expression profile for 21 genes out of 25 tested. Since only some of the offspring will have inherited the beneficial longevity profile from their long-lived parents, the contrast between offspring and controls is expected to be weak. Despite this dilution of the longevity effects, reduced expression levels of two genes, ASF1A and IL7R, involved in maintenance of chromatin structure and the immune system, associated with familial longevity already in middle-age. The size of this association increased when controls were compared to a subfraction of the offspring that had the highest probability to age healthily and become long-lived according to beneficial metabolic parameters. In conclusion, an “aging-signature” formed of 21 genes was identified, of which reduced expression of ASF1A and IL7R marked familial longevity already in middle-age. This indicates that expression changes of genes involved in metabolism, epigenetic control and immune function occur as a function of age, and some of these, like ASF1A and IL7R, represent early features of familial longevity and healthy ageing.

Genetic Analysis Workshop 17 (GAW17) provided a platform for evaluating existing statistical genetic methods and for developing novel methods to analyze rare variants that modulate complex traits. In this article, we present an overview of the 1000 Genomes Project exome data and simulated phenotype data that were distributed to GAW17 participants for analyses, the different issues addressed by the participants, and the process of preparation of manuscripts resulting from the discussions during the workshop.

Recently we proposed a novel two-step approach to test for pathway effects in disease progression. The goal of this approach is to study the joint effect of multiple single-nucleotide polymorphisms that belong to certain genes. By using random effects, our approach acknowledges the correlations within and between genes when testing for pathway effects. Gene-gene and gene-environment interactions can be included in the model. The method can be implemented with standard software, and the distribution of the test statistics under the null hypothesis can be approximated by using standard chi-square distributions. Hence no extensive permutations are needed for computations of the p-value. In this paper we adapt and apply the method to family data, and we study its performance for sequence data from Genetic Analysis Workshop 17. For the set of unrelated subjects, the performance of the new test was disappointing. We found a power of 6% for the binary outcome and of 18% for the quantitative trait Q1. For family data the new approach appears to perform well, especially for the quantitative outcome. We found a power of 39% for the binary outcome and a power of 89% for the quantitative trait Q1.

Analyzing sequencing data is difficult because of the low frequency of rare variants, which may result in low power to detect associations. We consider pathway analysis to detect multiple common and rare variants jointly and to investigate whether analysis at the pathway level provides an alternative strategy for identifying susceptibility genes. Available pathway analysis methods for data from genome-wide association studies might not be efficient because these methods are designed to detect common variants. Here, we investigate the performance of several existing pathway analysis methods for sequencing data. In particular, we consider the global test, which does not consider linkage disequilibrium between the variants in a gene. We improve the performance of the global test by assigning larger weights to rare variants, as proposed in the weighted-sum approach. Our conclusion is that straightforward application of pathway analysis is not satisfactory; hence, when common and rare variants are jointly analyzed, larger weights should be assigned to rare variants.

In genome-wide association studies (GWAS) of complex traits, single SNP analysis is still the most applied approach. However, the identified SNPs have small effects and provide limited biological insight. A more appropriate approach to interpret GWAS data of complex traits is to analyze the combined effect of a SNP set grouped per pathway or gene region. We used this approach to study the joint effect on human longevity of genetic variation in two candidate pathways, the insulin/insulin-like growth factor (IGF-1) signaling (IIS) pathway and the telomere maintenance (TM) pathway. For the analyses, we used genotyped GWAS data of 403 unrelated nonagenarians from long-lived sibships collected in the Leiden Longevity Study and 1,670 younger population controls. We analyzed 1,021 SNPs in 68 IIS pathway genes and 88 SNPs in 13 TM pathway genes using four self-contained pathway tests (PLINK set-based test, Global test, GRASS and SNP ratio test). Although we observed small differences between the results of the different pathway tests, they showed consistent significant association of the IIS and TM pathway SNP sets with longevity. Analysis of gene SNP sets from these pathways indicates that the association of the IIS pathway is scattered over several genes (AKT1, AKT3, FOXO4, IGF2, INS, PIK3CA, SGK, SGK2, and YWHAG), while the association of the TM pathway seems to be mainly determined by one gene (POT1). In conclusion, this study shows that genetic variation in genes involved in the IIS and TM pathways is associated with human longevity.

Electronic supplementary material

The online version of this article (doi:10.1007/s11357-011-9340-3) contains supplementary material, which is available to authorized users.

Cytokines are major immune system regulators. Previously, innate cytokine profiles determined by lipopolysaccharide stimulation were shown to be highly heritable. To identify regulating genes in innate immunity, we analyzed data from a genome-wide linkage scan using microsatellites in osteoarthritis (OA) patients (The GARP study) and their innate cytokine data on interleukin (IL)-1β, IL-1Ra, IL-10 and tumor necrosis factor (TNF)α. A confirmation cohort consisted of the Leiden 85-Plus study. In this study, a linkage analysis was followed by manual selection of candidate genes in linkage regions showing LOD scores over 2.5. An single-nucleotide polymorphism (SNP) gene tagging method was applied to select SNPs on the basis of the highest level of gene tagging and possible functional effects. QTDT was used to identify the SNPs associated with innate cytokine production. Initial association signals were modeled by a linear mixed model. Through these analyses, we identified 10 putative genes involved in the regulation of TNFα. SNP rs6679497 in gene CD53 showed significant association with TNFα levels (P=0.001). No association of this SNP was observed with OA. A novel gene involved in the innate immune response of TNFα is identified. Genetic variation in this gene may have a role in diseases and disorders in which TNFα is closely involved.

Background

Various cytokines and inflammatory mediators are known to be involved in the pathogenesis of rheumatoid arthritis (RA). We hypothesized that polymorphisms in selected inflammatory response and tissue repair genes contribute to the susceptibility to and severity of RA.

Methods

Polymorphisms in TNFA, IL1B, IL4, IL6, IL8, IL10, PAI1, NOS2a, C1INH, PARP, TLR2 and TLR4 were genotyped in 376 Caucasian RA patients and 463 healthy Caucasian controls using single base extension. Genotype distributions in patients were compared with those in controls. In addition, the association of polymorphisms with the need for anti-TNF-α treatment as a marker of RA severity was assessed.

Results

The IL8 781 CC genotype was associated with early onset of disease. The TNFA -238 G/A polymorphism was differentially distributed between RA patients and controls, but only when not corrected for age and gender. None of the polymorphisms was associated with disease severity.

Conclusions

We here report an association between IL8 781 C/T polymorphism and age of onset of RA. Our findings indicate that there might be a role for variations in genes involved in the immune response and in tissue repair in RA pathogenesis. Nevertheless, additional larger genomic and functional studies are required to further define their role in RA.

OBJECTIVE

At least 20 type 2 diabetes loci have now been identified, and several of these are associated with altered β-cell function. In this study, we have investigated the combined effects of eight known β-cell loci on insulin secretion stimulated by three different secretagogues during hyperglycemic clamps.

RESEARCH DESIGN AND METHODS

A total of 447 subjects originating from four independent studies in the Netherlands and Germany (256 with normal glucose tolerance [NGT]/191 with impaired glucose tolerance [IGT]) underwent a hyperglycemic clamp. A subset had an extended clamp with additional glucagon-like peptide (GLP)-1 and arginine (n = 224). We next genotyped single nucleotide polymorphisms in TCF7L2, KCNJ11, CDKAL1, IGF2BP2, HHEX/IDE, CDKN2A/B, SLC30A8, and MTNR1B and calculated a risk allele score by risk allele counting.

RESULTS

The risk allele score was associated with lower first-phase glucose-stimulated insulin secretion (GSIS) (P = 7.1 × 10−6). The effect size was equal in subjects with NGT and IGT. We also noted an inverse correlation with the disposition index (P = 1.6 × 10−3). When we stratified the study population according to the number of risk alleles into three groups, those with a medium- or high-risk allele score had 9 and 23% lower first-phase GSIS. Second-phase GSIS, insulin sensitivity index and GLP-1, or arginine-stimulated insulin release were not significantly different.

CONCLUSIONS

A combined risk allele score for eight known β-cell genes is associated with the rapid first-phase GSIS and the disposition index. The slower second-phase GSIS, GLP-1, and arginine-stimulated insulin secretion are not associated, suggesting that especially processes involved in rapid granule recruitment and exocytosis are affected in the majority of risk loci.

OBJECTIVE

Recently, results from a meta-analysis of genome-wide association studies have yielded a number of novel type 2 diabetes loci. However, conflicting results have been published regarding their effects on insulin secretion and insulin sensitivity. In this study we used hyperglycemic clamps with three different stimuli to test associations between these novel loci and various measures of β-cell function.

RESEARCH DESIGN AND METHODS

For this study, 336 participants, 180 normal glucose tolerant and 156 impaired glucose tolerant, underwent a 2-h hyperglycemic clamp. In a subset we also assessed the response to glucagon-like peptide (GLP)-1 and arginine during an extended clamp (n = 123). All subjects were genotyped for gene variants in JAZF1, CDC123/CAMK1D, TSPAN8/LGR5, THADA, ADAMTS9, NOTCH2/ADAMS30, DCD, VEGFA, BCL11A, HNF1B, WFS1, and MTNR1B.

RESULTS

Gene variants in CDC123/CAMK1D, ADAMTS9, BCL11A, and MTNR1B affected various aspects of the insulin response to glucose (all P < 6.9 × 10−3). The THADA gene variant was associated with lower β-cell response to GLP-1 and arginine (both P < 1.6 × 10−3), suggesting lower β-cell mass as a possible pathogenic mechanism. Remarkably, we also noted a trend toward an increased insulin response to GLP-1 in carriers of MTNR1B (P = 0.03), which may offer new therapeutic possibilities. The other seven loci were not detectably associated with β-cell function.

CONCLUSIONS

Diabetes risk alleles in CDC123/CAMK1D, THADA, ADAMTS9, BCL11A, and MTNR1B are associated with various specific aspects of β-cell function. These findings point to a clear diversity in the impact that these various gene variants may have on (dys)function of pancreatic β-cells.

Background

Fc gamma receptors (FcγRs) play a crucial role in immunity by linking IgG antibody-mediated responses with cellular effector and regulatory functions. Genetic variants in these receptors have been previously identified as risk factors for several chronic inflammatory conditions. The present study aimed to investigate the presence of copy number variations (CNVs) in the FCGR3B gene and its potential association with the autoimmune disease rheumatoid arthritis (RA).

Methodology/Principal Findings

CNV of the FCGR3B gene was studied using Multiplex Ligation Dependent Probe Amplification (MLPA) in 518 Dutch RA patients and 304 healthy controls. Surprisingly, three independent MLPA probes targeting the FCGR3B promoter measured different CNV frequencies, with probe#1 and #2 measuring 0 to 5 gene copies and probe#3 showing little evidence of CNV. Quantitative-PCR correlated with the copy number results from MLPA probe#2, which detected low copy number (1 copy) in 6.7% and high copy number (≥3 copies) in 9.4% of the control population. No significant difference was observed between RA patients and the healthy controls, neither in the low copy nor the high copy number groups (p-values = 0.36 and 0.71, respectively). Sequencing of the FCGR3B promoter region revealed an insertion/deletion (indel) that explained the disparate CNV results of MLPA probe#1. Finally, a non-significant trend was found between the novel -256A>TG indel and RA (40.7% in healthy controls versus 35.9% in RA patients; P = 0.08).

Conclusions/Significance

The current study highlights the complexity and poor characterization of the FCGR3B gene sequence, indicating that the design and interpretation of genotyping assays based on specific probe sequences must be performed with caution. Nonetheless, we confirmed the presence of CNV and identified novel polymorphisms in the FCGR3B gene in the Dutch population. Although no association was found between RA and FCGR3B CNV, the possible protective effect of the -256A>TG indel polymorphism must be addressed in larger studies.

Background

Markers for longevity that reflect the health condition and predict healthy aging are extremely scarce. Such markers are, however, valuable in aging research. It has been shown previously that the N-glycosylation pattern of human immunoglobulin G (IgG) is age-dependent. Here we investigate whether N-linked glycans reflect early features of human longevity.

Methodology/Principal Findings

The Leiden Longevity Study (LLS) consists of nonagenarian sibling pairs, their offspring, and partners of the offspring serving as control. IgG subclass specific glycosylation patterns were obtained from 1967 participants in the LLS by MALDI-TOF-MS analysis of tryptic IgG Fc glycopeptides. Several regression strategies were applied to evaluate the association of IgG glycosylation with age, sex, and longevity. The degree of galactosylation of IgG decreased with increasing age. For the galactosylated glycoforms the incidence of bisecting GlcNAc increased as a function of age. Sex-related differences were observed at ages below 60 years. Compared to males, younger females had higher galactosylation, which decreased stronger with increasing age, resulting in similar galactosylation for both sexes from 60 onwards. In younger participants (<60 years of age), but not in the older age group (>60 years), decreased levels of non-galactosylated glycoforms containing a bisecting GlcNAc reflected early features of longevity.

Conclusions/Significance

We here describe IgG glycoforms associated with calendar age at all ages and the propensity for longevity before middle age. As modulation of IgG effector functions has been described for various IgG glycosylation features, a modulatory effect may be expected for the longevity marker described in this study.

Our aim is to develop methods for mapping genes related to age at onset in general pedigrees. We propose two score tests, one derived from a gamma frailty model with pairwise likelihood and one derived from a log-normal frailty model with approximated likelihood around the null random effect. The score statistics are weighted nonparametric linkage statistics, with weights depending on the age at onset. These tests are correct under the null hypothesis irrespective of the weight used. They are simple, robust, computationally fast, and can be applied to large, complex pedigrees. We apply these methods to simulated data and to the Genetic Analysis Workshop 16 Framingham Heart Study data set. We investigate the time to the first of three events: hard coronary heart disease, diabetes, or death from any cause. We use a two-step procedure. In the first step, we estimate the population parameters under the null hypothesis of no linkage. In the second step, we apply the score tests, using the population parameters estimated in the first step.

We describe an empirical Bayesian linear model for integration of functional gene annotation data with genome-wide association data. Using case-control study data from the North American Rheumatoid Arthritis Consortium and gene annotation data from the Gene Ontology, we illustrate how the method can be used to prioritize candidate genes for further investigation.

We investigated efficient case-control association analysis using family data. The outcome of interest was coronary heart disease. We employed existing and new methods that take into account the correlations among related individuals to obtain the proper type I error rates. The methods considered for autosomal single-nucleotide polymorphisms were: 1) generalized estimating equations-based methods, 2) variance-modified Cochran-Armitage (MCA) trend test incorporating kinship coefficients, and 3) genotypic modified quasi-likelihood score test. Additionally, for X-linked single-nucleotide polymorphisms we proposed a two-degrees-of-freedom test. Performance of these methods was tested using Framingham Heart Study 500 k array data.

Background

In haplotype-based candidate gene studies a problem is that the genotype data are unphased, which results in haplotype ambiguity. The measure [1] quantifies haplotype predictability from genotype data. It is computed for each individual haplotype, and for a measure of global relative efficiency a minimum value is suggested. Alternatively, we developed methods directly based on the information content of haplotype frequency estimates to obtain global relative efficiency measures: and based on A- and D-optimality, respectively. All three methods are designed for single populations; they can be applied in cases only, controls only or the whole data. Therefore they are not necessarily optimal for haplotype testing in case-control studies.

Results

A new global relative efficiency measure was derived to maximize power of a simple test statistic that compares haplotype frequencies in cases and controls. Application to real data showed that our proposed method gave a clear and summarizing measure for the case-control study conducted. Additionally this measure might be used for selection of individuals, who have the highest potential for improving power by resolving phase ambiguity.

Conclusion

Instead of using relative efficiency measure for cases only, controls only or their combined data, we link uncertainty measure to case-control studies directly. Hence, our global efficiency measure might be useful to assess whether data are informative or have enough power for estimation of a specific haplotype risk.

Background

Elevated circulating levels of C-reactive protein (CRP), interleukin (IL)-6 and fibrinogen (FG) have been repeatedly associated with many adverse outcomes in patients with chronic obstructive pulmonary disease (COPD). To date, it remains unclear whether and to what extent systemic inflammation is primary or secondary in the pathogenesis of COPD.

The aim of this study was to examine the association between haplotypes of CRP, IL6 and FGB genes, systemic inflammation, COPD risk and COPD-related phenotypes (respiratory impairment, exercise capacity and body composition).

Methods

Eighteen SNPs in three genes, representing optimal haplotype-tagging sets, were genotyped in 355 COPD patients and 195 healthy smokers. Plasma levels of CRP, IL-6 and FG were measured in the total study group. Differences in haplotype distributions were tested using the global and haplotype-specific statistics.

Results

Raised plasma levels of CRP, IL-6 and fibrinogen were demonstrated in COPD patients. However, COPD population was very heterogeneous: about 40% of patients had no evidence of systemic inflammation (CRP < 3 mg/uL or no inflammatory markers in their top quartile). Global test for haplotype effect indicated association of CRP gene and CRP plasma levels (P = 0.0004) and IL6 gene and COPD (P = 0.003). Subsequent analysis has shown that IL6 haplotype H2, associated with an increased COPD risk (p = 0.004, OR = 4.82; 1.64 to 4.18), was also associated with very low CRP levels (p = 0.0005). None of the genes were associated with COPD-related phenotypes.

Conclusion

Our findings suggest that common genetic variation in CRP and IL6 genes may contribute to heterogeneity of COPD population associated with systemic inflammation.