The -308 tumour necrosis factor-alpha gene polymorphism predicts therapeutic response to TNFalpha-blockers in rheumatoid arthritis and spondyloarthritis patients

OBJECTIVE: To examine whether the G-to-A polymorphism at position -308 in the promoter of the tumour necrosis factor-alpha (TNFalpha) gene influences the therapeutic response to TNFalpha-blockers in patients with rheumatoid arthritis (RA), psoriatic arthritis (PsA) and ankylosing spondylitis (AS). METHODS: A total of 54 patients with RA, 10 with PsA and 22 with AS were genotyped by polymerase chain reaction for the -308 TNFalpha promoter polymorphism. They were treated with infliximab (n = 63), adalimumab (n = 10) or etanercept (n = 13). Clinical response was assessed after 24 weeks by the Disease Activity Score in 28 joints (DAS28) for RA and PsA, and the Bath Ankylosing Spondylitis Activity Index (BASDAI) for AS patients. RESULTS: All patients with the A/A genotype (n = 3, all RA) and two patients with the A/G genotype (AS) failed to respond to anti-TNF treatment. Irrespective of the underlying disease, moderate response (n = 44) was predominantly associated with the A/G genotype (A/G 18/22, G/G 4/22), whereas good response (n = 59) was exclusively seen in patients with the G/G genotype. The average improvement in the DAS28 score was 0.83 in the A/A, 1.50 in the A/G and 2.64 in the G/G group of RA and PsA patients (P < 0.0001). The BASDAI score in AS improved on average by 1.21 in the A/G and by 3.30 in the G/G group (P < 0.005). CONCLUSIONS: The data suggest that humans with a TNFalpha -308 G/G genotype are better responders to anti-TNFalpha treatment than those with A/A or A/G genotypes independent of the treated rheumatic disease (RA, PsA or AS).

The common G-allele of interleukin-18 single-nucleotide polymorphism is a genetic risk factor for atopic asthma. The SAPALDIA Cohort Study

BACKGROUND: IL-18 is a pleiotrophic cytokine involved in both, T-helper type 1 (Th1) and Th2 differentiation. Recently genetic variants in the IL-18 gene have been associated with increased risk of atopy and asthma. OBJECTIVE: To examine the relationship of a genetic, haplotype-tagging promotor variant -137G/C in the IL-18 gene with atopic asthma in a large, well-characterized and population-based study of adults. METHODS: Prospective cohort study design was used to collect interview and biological measurement data at two examination time-points 11 years apart. Multivariate logistic regression analysis was used to assess the association of genotype with asthma and atopy. RESULTS: The G-allele of the IL-18 promotor variant (-137G/C) was associated with a markedly increased risk for the prevalence of physician-diagnosed asthma with concomitant skin reactivity to common allergens. Stratification of the asthma cases by skin reactivity to common allergens revealed an exclusive association of IL-18 -137 G-allele with an increased prevalence of atopic asthma (adjusted odds ratio (OR): 3.63; 95% confidence interval: (1.64-8.02) for GC or GG carriers vs. CC carriers), and no according association with asthma and concomitant negative skin reactivity (adjusted OR: 1.13; 0.66-1.94). The interaction between IL-18 -137G/C genotype and positive skin prick test was statistically significant (P=0.029). None of 74 incident asthma cases with atopy at baseline exhibited the CC genotype. CONCLUSION: Our results strongly suggest that this variant of the IL-18 gene is an important genetic determinant involved in the development of atopic asthma.

TRIO LOGIC REGRESSION - DETECTION OF SNP - SNP INTERACTIONS IN CASE-PARENT TRIOS

Statistical approaches to evaluate higher order SNP-SNP and SNP-environment interactions are critical in genetic association studies, as susceptibility to complex disease is likely to be related to the interaction of multiple SNPs and environmental factors. Logic regression (Kooperberg et al., 2001; Ruczinski et al., 2003) is one such approach, where interactions between SNPs and environmental variables are assessed in a regression framework, and interactions become part of the model search space. In this manuscript we extend the logic regression methodology, originally developed for cohort and case-control studies, for studies of trios with affected probands. Trio logic regression accounts for the linkage disequilibrium (LD) structure in the genotype data, and accommodates missing genotypes via haplotype-based imputation. We also derive an efficient algorithm to simulate case-parent trios where genetic risk is determined via epistatic interactions.

Role of the NFKB1 -94ins/delATTG promoter polymorphism in IBD and potential interactions with polymorphisms in the CARD15/NOD2, IKBL, and IL-1RN genes

BACKGROUND: Recently, an association of the NFKB1 polymorphism -94ins/delATTG with ulcerative colitis (UC) has been reported. This 4-bp insertion/deletion polymorphism is localized in the promoter region of the NFKB1 gene and appears to be functionally relevant. The aim of the present study was to confirm the association of the -94ins/delATTG (W/D) NFKB1 promoter polymorphism with UC in a population of German origin and to test for a potential association with Crohn's disease (CD). Furthermore, potential interactions of the -94ins/delATTG polymorphism with the IKBL and the IL-1RN genes should be determined. MATERIALS AND METHODS: The study population comprised 630 patients with CD, 365 patients with UC, and 974 healthy controls. Genotyping was performed using polymerase chain reaction and restriction fragment length polymorphism analysis. For statistical evaluation, the chi-square test and the Fisher exact test were used. RESULTS: No significant association of the W/D NFKB1 polymorphism with CD or UC was detected. In addition, no significant interactions between the -94ins/delATTG NFKB1 polymorphism and polymorphisms within the IKBL and the IL-1RN genes, respectively, were found in CD or UC. Also, no significant interactions of the NFKB1 polymorphism with mutations of the CARD15/NOD2 gene and with clinical phenotypes were detected in CD. Moreover, no associations of the NFKB1 polymorphism were found in UC depending on disease localization. CONCLUSIONS: The present study could not confirm the reported association of the -94ins/delATTG NFKB1 polymorphism with UC and also found no evidence for a role of this polymorphism in CD. The results do not give evidence for a role of this NFKB1 polymorphism in the pathogenesis of UC and CD.

Microbiological composition associated with interleukin-1 gene polymorphism in subjects undergoing supportive periodontal therapy

BACKGROUND: Interleukin-1 gene polymorphism (IL-1 gene) has been associated with periodontitis. The present study examined the subgingival microbiota by IL-1 gene status in subjects undergoing supportive periodontal therapy (SPT). METHODS: A total of 151 subjects with known IL-1 gene status (IL-1A +4845/IL-1B -3954) (IL-1 gene) were included in this study. Clinical data and subgingival plaque samples (40 taxa) were collected. These taxa were determined by the checkerboard DNA-DNA hybridization method. RESULTS: Gender, smoking habits (n-par tests), age, and clinical periodontal conditions did not differ by IL-1 gene status. IL-1 gene-negative subjects had a higher total bacterial load (mean difference, 480.4 x 10(5); 95% confidence interval [CI], 77 to 884 x 10(5); P <0.02). The levels of Actinobacillus actinomycetemcomitans (mean difference, 30.7 x 10(5); 95% CI, 2.2 to 59.5 x 10(5); P <0.05), Eubacterium nodatum (mean difference, 4.2 x 10(5); 95% CI, 0.6 to 7.8 x 10(5); P <0.02), Porphyromonas gingivalis (mean difference, 17.9 x 10(5); 95% CI, 1.2 to 34.5 x 10(5); P <0.05), and Streptococcus anginosus (mean difference, 4.0 x 10(5); 95% CI, 0.2 to 7.2 x 10(5); P <0.05) were higher in IL-1 gene-negative subjects, an observation specifically found at sites with probing depths <5.0 mm. CONCLUSIONS: Bleeding on probing did not differ by IL gene status, reflecting clinical SPT efficacy. IL-1 gene-negative subjects had higher levels of periodontal pathogens. This may suggest that among subjects undergoing SPT, a lower bacterial load is required in IL-1 gene-positive subjects to develop the same level of periodontitis as in IL-1 gene-negative subjects.

A MULTILEVEL MODEL TO ADDRESS BATCH EFFECTS IN COPY NUMBER ESTIMATION USING SNP ARRAYS

Submicroscopic changes in chromosomal DNA copy number dosage are common and have been implicated in many heritable diseases and cancers. Recent high-throughput technologies have a resolution that permits the detection of segmental changes in DNA copy number that span thousands of basepairs across the genome. Genome-wide association studies (GWAS) may simultaneously screen for copy number-phenotype and SNP-phenotype associations as part of the analytic strategy. However, genome-wide array analyses are particularly susceptible to batch effects as the logistics of preparing DNA and processing thousands of arrays often involves multiple laboratories and technicians, or changes over calendar time to the reagents and laboratory equipment. Failure to adjust for batch effects can lead to incorrect inference and requires inefficient post-hoc quality control procedures that exclude regions that are associated with batch. Our work extends previous model-based approaches for copy number estimation by explicitly modeling batch effects and using shrinkage to improve locus-specific estimates of copy number uncertainty. Key features of this approach include the use of diallelic genotype calls from experimental data to estimate batch- and locus-specific parameters of background and signal without the requirement of training data. We illustrate these ideas using a study of bipolar disease and a study of chromosome 21 trisomy. The former has batch effects that dominate much of the observed variation in quantile-normalized intensities, while the latter illustrates the robustness of our approach to datasets where as many as 25% of the samples have altered copy number. Locus-specific estimates of copy number can be plotted on the copy-number scale to investigate mosaicism and guide the choice of appropriate downstream approaches for smoothing the copy number as a function of physical position. The software is open source and implemented in the R package CRLMM available at Bioconductor (http:www.bioconductor.org).

EXPLORATION, NORMALIZATION, AND GENOTYPE CALLS OF HIGH DENSITY OLIGONUCLEOTIDE SNP ARRAY DATA

In most microarray technologies, a number of critical steps are required to convert raw intensity measurements into the data relied upon by data analysts, biologists and clinicians. These data manipulations, referred to as preprocessing, can influence the quality of the ultimate measurements. In the last few years, the high-throughput measurement of gene expression is the most popular application of microarray technology. For this application, various groups have demonstrated that the use of modern statistical methodology can substantially improve accuracy and precision of gene expression measurements, relative to ad-hoc procedures introduced by designers and manufacturers of the technology. Currently, other applications of microarrays are becoming more and more popular. In this paper we describe a preprocessing methodology for a technology designed for the identification of DNA sequence variants in specific genes or regions of the human genome that are associated with phenotypes of interest such as disease. In particular we describe methodology useful for preprocessing Affymetrix SNP chips and obtaining genotype calls with the preprocessed data. We demonstrate how our procedure improves existing approaches using data from three relatively large studies including one in which large number independent calls are available. Software implementing these ideas are avialble from the Bioconductor oligo package.

EFFICIENT EVALUATION OF RANKING PROCEDURES WHEN THE NUMBER OF UNITS IS LARGE WITH APPLICATION TO SNP IDENTIFICATION

Simulation-based assessment is a popular and frequently necessary approach to evaluation of statistical procedures. Sometimes overlooked is the ability to take advantage of underlying mathematical relations and we focus on this aspect. We show how to take advantage of large-sample theory when conducting a simulation using the analysis of genomic data as a motivating example. The approach uses convergence results to provide an approximation to smaller-sample results, results that are available only by simulation. We consider evaluating and comparing a variety of ranking-based methods for identifying the most highly associated SNPs in a genome-wide association study, derive integral equation representations of the pre-posterior distribution of percentiles produced by three ranking methods, and provide examples comparing performance. These results are of interest in their own right and set the framework for a more extensive set of comparisons.

Assessment of use of microsatellite polymorphism analysis for improving spatial distribution tracking of echinococcus multilocularis

Alveolar echinococcosis (AE)--caused by the cestode Echinococcus multilocularis--is a severe zoonotic disease found in temperate and arctic regions of the northern hemisphere. Even though the transmission patterns observed in different geographical areas are heterogeneous, the nuclear and mitochondrial targets usually used for the genotyping of E. multilocularis have shown only a marked genetic homogeneity in this species. We used microsatellite sequences, because of their high typing resolution, to explore the genetic diversity of E. multilocularis. Four microsatellite targets (EmsJ, EmsK, and EmsB, which were designed in our laboratory, and NAK1, selected from the literature) were tested on a panel of 76 E. multilocularis samples (larval and adult stages) obtained from Alaska, Canada, Europe, and Asia. Genetic diversity for each target was assessed by size polymorphism analysis. With the EmsJ and EmsK targets, two alleles were found for each locus, yielding two and three genotypes, respectively, discriminating European isolates from the other groups. With NAK1, five alleles were found, yielding seven genotypes, including those specific to Tibetan and Alaskan isolates. The EmsB target, a tandem repeated multilocus microsatellite, found 17 alleles showing a complex pattern. Hierarchical clustering analyses were performed with the EmsB findings, and 29 genotypes were identified. Due to its higher genetic polymorphism, EmsB exhibited a higher discriminatory power than the other targets. The complex EmsB pattern was able to discriminate isolates on a regional and sectoral level, while avoiding overdistinction. EmsB will be used to assess the putative emergence of E. multilocularis in Europe.

Beta-D-glucoside utilization by Mycoplasma mycoides subsp. mycoides SC: possible involvement in the control of cytotoxicity towards bovine lung cells

BACKGROUND: Contagious bovine pleuropneumonia (CBPP) caused by Mycoplasma mycoides subsp. mycoides small-colony type (SC) is among the most serious threats for livestock producers in Africa. Glycerol metabolism-associated H2O2 production seems to play a crucial role in virulence of this mycoplasma. A wide number of attenuated strains of M. mycoides subsp. mycoides SC are currently used in Africa as live vaccines. Glycerol metabolism is not affected in these vaccine strains and therefore it does not seem to be the determinant of their attenuation. A non-synonymous single nucleotide polymorphism (SNP) in the bgl gene coding for the 6-phospho-beta-glucosidase (Bgl) has been described recently. The SNP differentiates virulent African strains isolated from outbreaks with severe CBPP, which express the Bgl isoform Val204, from strains to be considered less virulent isolated from CBPP outbreaks with low mortality and vaccine strains, which express the Bgl isoform Ala204. RESULTS: Strains of M. mycoides subsp. mycoides SC considered virulent and possessing the Bgl isoform Val204, but not strains with the Bgl isoform Ala204, do trigger elevated levels of damage to embryonic bovine lung (EBL) cells upon incubation with the disaccharides (i.e., beta-D-glucosides) sucrose and lactose. However, strains expressing the Bgl isoform Val204 show a lower hydrolysing activity on the chromogenic substrate p-nitrophenyl-beta-D-glucopyranoside (pNPbG) when compared to strains that possess the Bgl isoform Ala204. Defective activity of Bgl in M. mycoides subsp. mycoides SC does not lead to H2O2 production. Rather, the viability during addition of beta-D-glucosides in medium-free buffers is higher for strains harbouring the Bgl isoform Val204 than for those with the isoform Ala204. CONCLUSION: Our results indicate that the studied SNP in the bgl gene is one possible cause of the difference in bacterial virulence among strains of M. mycoides subsp. mycoides SC. Bgl does not act as a direct virulence factor, but strains possessing the Bgl isoform Val204 with low hydrolysing activity are more prone to survive in environments that contain high levels of beta-D-glucosides, thus contributing in some extent to mycoplasmaemia.

Functional polymorphism in ABCA1 influences age of symptom onset in coronary artery disease patients

ATP-binding-cassette-transporter-A1 (ABCA1) plays a pivotal role in intracellular cholesterol removal, exerting a protective effect against atherosclerosis. ABCA1 gene severe mutations underlie Tangier disease, a rare Mendelian disorder that can lead to premature coronary artery disease (CAD), with age of CAD onset being two decades earlier in mutant homozygotes and one decade earlier in heterozygotes than in mutation non-carriers. It is unknown whether common polymorphisms in ABCA1 could influence age of symptom onset of CAD in the general population. We examined common promoter and non-synonymous coding polymorphisms in relation to age of symptom onset in a group of CAD patients (n = 1164), and also carried out in vitro assays to test effects of the promoter variations on ABCA1 promoter transcriptional activity and effects of the coding variations on ABCA1 function in mediating cellular cholesterol efflux. Age of symptom onset was found to be associated with the promoter - 407G > C polymorphism, being 2.82 years higher in C allele homozygotes than in G allele homozygotes and intermediate in heterozygotes (61.54, 59.79 and 58.72 years, respectively; P = 0.002). In agreement, patients carrying ABCA1 haplotypes containing the -407C allele had higher age of symptom onset. Patients of the G/G or G/C genotype of the -407G > C polymorphism had significant coronary artery stenosis (>75%) at a younger age than those of the C/C genotype (P = 0.003). Reporter gene assays showed that ABCA1 haplotypes bearing the -407C allele had higher promoter activity than haplotypes with the -407G allele. Functional analyses of the coding polymorphisms showed an effect of the V825I substitution on ABCA1 function, with the 825I variant having higher activity in mediating cholesterol efflux than the wild-type (825V). A trend towards higher symptom onset age in 825I allele carriers was observed. The data indicate an influence of common ABCA1 functional polymorphisms on age of symptom onset in CAD patients.

Heterosis for biomass-related traits in Arabidopsis investigated by quantitative trait loci analysis of the triple testcross design with recombinant inbred lines

Arabidopsis thaliana has emerged as a leading model species in plant genetics and functional genomics including research on the genetic causes of heterosis. We applied a triple testcross (TTC) design and a novel biometrical approach to identify and characterize quantitative trait loci (QTL) for heterosis of five biomass-related traits by (i) estimating the number, genomic positions, and genetic effects of heterotic QTL, (ii) characterizing their mode of gene action, and (iii) testing for presence of epistatic effects by a genomewide scan and marker x marker interactions. In total, 234 recombinant inbred lines (RILs) of Arabidopsis hybrid C24 x Col-0 were crossed to both parental lines and their F1 and analyzed with 110 single-nucleotide polymorphism (SNP) markers. QTL analyses were conducted using linear transformations Z1, Z2, and Z3 calculated from the adjusted entry means of TTC progenies. With Z1, we detected 12 QTL displaying augmented additive effects. With Z2, we mapped six QTL for augmented dominance effects. A one-dimensional genome scan with Z3 revealed two genomic regions with significantly negative dominance x additive epistatic effects. Two-way analyses of variance between marker pairs revealed nine digenic epistatic interactions: six reflecting dominance x dominance effects with variable sign and three reflecting additive x additive effects with positive sign. We conclude that heterosis for biomass-related traits in Arabidopsis has a polygenic basis with overdominance and/or epistasis being presumably the main types of gene action.