Program report: GENECOUNTING support programs

Evaluation of Common Variants in Matrix Metalloproteinase-9 Gene with Lumbar Disc Herniation in Han Chinese Population

OBJECTIVE

Lumbar disc herniation (LDH) is a common and frequent orthopedic disease with strong genetic determinants. The disruption of the intervertebral disc extracellular matrix has been found to play a key role in the development of LDH, suggesting that abnormal matrix metalloproteinases (MMPs) may promote the degradation of the disc matrix. MMP-9, an important member of the MMP family, is a good candidate for the LDH susceptibility gene. The present study aimed to investigate the association of common variants in the MMP-9 gene with the risk, severity, and clinical characteristic variables of LDH.

MATERIALS AND METHODS

Fourteen tag single nucleotide polymorphisms (SNPs) entirely covering the region of the MMP-9 gene were analyzed in a sample of 845 patients and 1751 healthy controls.

RESULTS

The SNP rs17576 was found to be significantly associated with susceptibility to LDH (OR = 0.77, p = 0.0002), which was also confirmed by haplotype-based analyses (rs79845319-rs17576-rs45437897, global p < 0.001). Our results indicated that the A allele of rs17576 reduced the risk of LDH by ∼23% on average. Furthermore, the G allele of rs17576 was found to correlate with more severe grades of disc degeneration.

CONCLUSION

Our results provide additional evidence supporting an important role of the MMP-9 gene in the pathogenesis of LDH.

Genetic polymorphism of matrix metalloproteinase-1 and coronary artery disease susceptibility: a case-control study in a Han Chinese population

Coronary artery disease (CAD) receives intensive research due to its high incidence and severe impact on the quality of life. One member of the matrix metalloproteinases (MMPs), MMP-1, has been reported to be associated with CAD. To identify the markers contributing to the genetic susceptibility to CAD, nine single-nucleotide polymorphisms (rs1799750, rs498186, rs475007, rs514921, rs494379, rs996999, rs2071232, rs1938901, and rs2239008) throughout the MMP-1 gene were genotyped using MALDI-TOF within the MassARRAY system, and the allele and genotype distributions were compared between 438 healthy controls and 411 patients with CAD from a Chinese Han population. The analysis revealed a weak association between the rs1799750 (in the promoter region) genotype distribution and CAD (p=0.022). An increased risk of CAD was significantly associated with the 2G allele of rs1799750 (p=0.005, odds ratio=1.329, 95% confidence interval=1.090-1.620, after Bonferroni corrections). Strong linkage disequilibrium was observed in three blocks (D'>0.9). Significantly more C-2G (rs498186-rs1799750) haplotypes (p=0.001 after Bonferroni corrections) were found in CAD subjects. These findings point to a role for the polymorphism in the MMP-1 promoter in CAD among a Han Chinese population and may be informative for future genetic or biological studies on CAD.

Analysis of interleukin-8 gene variants reveals their relative importance as genetic susceptibility factors for chronic periodontitis in the Han population

Interleukin (IL)-8, an important chemokine that regulates the inflammatory response, plays an important role in periodontitis. Previous studies indicate that certain IL-8 gene polymorphisms are associated with periodontitis susceptibility in some populations. However, the literature is somewhat contradictory, and not all IL-8 polymorphisms have been examined, particularly in Han Chinese individuals. The aim of this study was to investigate the association of every IL-8 SNP with chronic periodontitis in Han Chinese individuals. We analyzed 23 SNPs with minor allele frequency (MAF)≥0.01, which were selected from 219 SNPs in the NCBI dbSNP and preliminary HapMap data analyses from a cohort of 400 cases and 750 controls from genetically independent Han Chinese individuals. Single SNP, haplotype and gender-specific associations were performed. We found that rs4073 and rs2227307 were significantly associated with chronic periodontitis. Further haplotype analysis indicated that a haplotype block (rs4073-rs2227307-rs2227306) that spans the promoter and exon1 of IL-8 was highly associated with chronic periodontitis. Additionally, the ATC haplotype in this block was increased 1.5-fold in these cases. However, when analyzing the samples by gender, no significant gender-specific associations in IL-8 were observed, similar to the results of haplotype association analyses in female and male subgroups. Our results provide further evidence that IL-8 is associated with chronic periodontitis in Han Chinese individuals. Furthermore, our results confirm previous reports suggesting the intriguing possibilities that IL-8 plays a role in the pathogenesis of chronic periodontitis and that this gene may be involved in the etiology of this condition.

Annexin A11 (ANXA11) gene polymorphisms are associated with sarcoidosis in a Han Chinese population: a case-control study

OBJECTIVES

To further identify the single-nucleotide polymorphisms (SNPs) that contribute to the genetic susceptibility to sarcoidosis, we examined the potential association between sarcoidosis and 15 SNPs of the ANXA11 gene.

DESIGN

A case-control study.

SETTING

A tuberculosis unit in a hospital of the university in China.

PARTICIPANTS

Participants included 412 patients with sarcoidosis and 418 healthy controls.

METHODS

The selected SNPs were genotyped using the MALDI-TOF in the MassARRAY system.

RESULTS

Statistically significant differences were found in the allelic or genotypic frequencies of the rs2789679, rs1049550 and rs2819941 in the ANXA11 gene between patients with sarcoidosis and controls. The rs2789679 A allele (p=0.00004, OR=1.42, 95% CI 1.17 to 1.73) and rs2819941 T allele (p=0.0006, OR=1.41, 95% CI 1.16 to 1.71) were significantly more frequent in patients with sarcoidosis compared with controls. The frequency of the rs1049550 T allele (p=0.000002, OR=0.61, 95% CI 0.49 to 0.74) in patients with sarcoidosis was significantly lower than that in controls. The multi-SNP model reveals that rs1049550 is the only independent SNP association effect after accounting for the other two marginally associated SNPs. In block 2 (rs1049550-rs2573351), the T-C haplotype occurred significantly less frequently (p=0.001), whereas the C-C haplotypes occurred more frequently (p=0.0001) in patients with sarcoidosis than controls. Furthermore, genotype frequency distribution revealed that, in rs1049550, the CC genotype was significantly more in patients with chest X-ray (CXR) stage I sarcoidosis than in patients with CXR stage II-IV sarcoidosis (p=0.012).

CONCLUSIONS

These findings point to a role for the polymorphisms of ANXA11 in sarcoidosis in a Chinese Han population, and may be informative for future genetic studies on sarcoidosis.

4q22.1 contributes to bone mineral density and osteoporosis susceptibility in postmenopausal women of Chinese Han population

Osteoporosis is a multifactorial disease in which genetic determinants are modulated by hormonal, environmental and nutritional factors. An important clinical risk factor in the pathogenesis of osteoporosis is the presence of genetics polymorphism in/around susceptibility genes/regions. This study explored whether the region of 4q22.1, which confers risk of developing osteoporosis in some populations, associated with bone mineral density and osteoporosis susceptibility in postmenopausal women of Han Chinese. We investigated 32 SNPs with minor allele frequencies ≥0.05 between 20 kb upstream and 20 kb downstream (40 kb window) of rs6532023, mapping in the 4q22.1 region, which was reported to be significantly associated with osteoporosis in previous studies. We found that rs6532023 was significantly associated with bone mineral density and osteoporosis (corrected p = 0.015) in our sample, including 440 cases and 640 controls, and allele G was supposed as a risk factor while T worked as a protective factor. Further genotype association analyses suggested a similar pattern (corrected p = 0.040). Additionally, analyses by haplotypes indicated that a haplotype block rs7683315-rs6532023-rs1471400-rs1471403 in the region associated with bone mineral density and osteoporosis (global p = 0.032), and risk haplotype A-G-G-C had almost 1.5-fold increased in the cases. To our knowledge, this is the first report to examine 4q22.1 region polymorphisms and osteoporosis in Han Chinese. Our results provide further evidence for an effect of the region of 4q22.1 on the etiology of osteoporosis and suggest that 4q22.1 may be a genetic risk factor for bone mineral density and osteoporosis.

A population-based association study of 2q32.3 and 8q21.3 loci with schizophrenia in Han Chinese

It has been reported that two new schizophrenia susceptibility loci (2q32.3 and 8q21.3) containing two single-nucleotide polymorphisms (SNPs; rs17662626 and rs7004633) have been identified in Europeans by a genome-wide association study. To determine if the two regions are associated with schizophrenia in Han Chinese, which are distinct from Europeans, we analyzed 9 SNPs, including rs17662626 and rs7004633, within 2 regions involving 1430 cases and 1570 controls from the Han population. Single SNP association, haplotype association, and gender-specific association analyses were performed. Single SNP analyses revealed that there was no association with schizophrenia in the region of 2q32.3, but rs7004633 mapping to the region of 8q21.3 was significantly associated with schizophrenia (p = 5.1 × 10(-5); OR = 1.274; 95% CI, 1.134-1.429). Further genotype and haplotype association analyses for the region of 8q21.3 suggested a similar pattern. Additionally, analyses by haplotypes indicated that a haplotype block in the region of 8q21.3 highly associated with schizophrenia and one haplotype in this haploblock had a 1.5-fold increase in the cases. Our results provide further evidence regarding the association of the region of 8q21.3 with schizophrenia in Han Chinese, as well as Europeans, which confirmed the previous report that 8q21.3 may play important roles in the etiology of schizophrenia.

Association study of a new schizophrenia susceptibility locus of 10q24.32-33 in a Han Chinese population

Recently, a new schizophrenia susceptibility locus 10q24.32-q24.33, containing two single-nucleotide polymorphisms (SNPs: rs7914558 and rs11191580), was identified in a genome-wide association study. To examine if this locus is associated with schizophrenia in the Han Chinese population, we analyzed six SNPs, including two SNPs within the region of interest. The sample consisted of 1430 schizophrenia cases and 1570 controls from genetically independent members of the Han population. Single-SNP association, haplotype association and sex-specific association analyses were performed. Three SNPs, rs7914558 (p=1.41×10(-4); OR=1.11; 95% CI 1.05-1.17), rs12220375 (p=1.18×10(-4); OR=1.06; 95% CI 1.03-1.09) and rs11191580 (p=3.03×10(-4); OR=1.05; 95% CI 1.02-1.10), mapped to the locus and were significantly associated in our sample set. Further genotype and haplotype association analyses suggested a similar pattern. Similar to results from European populations, our results provide further evidence that this region associated with schizophrenia in Han Chinese. Results also confirm previous reports suggesting that 10q24.32-q24.33 offers an intriguing new insight into the pathogenesis of schizophrenia and may play an important role in its etiology.

Tests of linkage and allelic association between markers in the 1p36 PRKCZ (protein kinase C zeta) gene region and bipolar affective disorder

Three linkage studies of families with multiple cases of bipolar disorder and/or unipolar affective disorder have confirmed the involvement of the chromosome 1p36 region in the etiology of affective disorders with LOD scores of 2.7, 3.6, and 3.97. We investigated the protein kinase C zeta gene (PRKCZ) as a susceptibility locus for bipolar disorder because it is highly brain expressed and is localized close to the marker D1S243 which was linked to affective disorder in a single large UCL bipolar disorder family with a LOD of 3.1. PRKCZ encodes an unusual type of protein kinase which affects axonal differentiation through Wnt-signaling. We genotyped four microsatellite markers and nine single nucleotide polymorphism (SNP) markers within or near the PRKCZ gene in the UCL case-control sample of 600 bipolar disorder patients and up to 605 supernormal controls. Markers D1S243 and rs3128396 were significantly associated with bipolar disorder (empirical P = 0.037 and P = 0.040, respectively). We also included data from eight SNPs which were genotyped as part of our GWA study on bipolar disorder for association analysis. Tests of haplotypic association found that a haplotype block comprising markers rs3128296, rs2503706, and rs3128309 was associated with bipolar disorder (empirical P = 0.004). A previous linkage study had shown greater evidence for linkage within female cases compared to males. Therefore, to assess if the association was sex-specific, we performed a female-only allelic-association analysis, which resulted in SNPs rs3128296 and rs3128309 becoming associated with bipolar disorder (P = 0.004 and P = 0.016, respectively). PRKCZ may play a role in susceptibility to bipolar affective disorder.

Association of PDE4B polymorphisms and schizophrenia in Northwestern Han Chinese

The phosphodiesterase 4B (PDE4B) is a candidate susceptibility gene for schizophrenia (SCZ), interacting with DISC1, a known genetic risk factor for SCZ. To examine if variants within PDE4B gene are associated with SCZ in Northwestern Han Chinese, and if these effects vary in gender-specific subgroup, we analyzed 20 SNPs, selected from previous studies and preliminary HapMap data analyses with minor allele frequency (MAF) ≥ 20%, in a cohort of 428 cases and 572 controls from genetically independent Northwestern Han Chinese. Single SNP association, haplotype association and sex-specific association analysis were performed. We found that rs472952 is significantly associated with SCZ and rs7537440 is associated with SCZ in females. Further analysis indicated that a haplotype block spanning PDE4B2 splice site is highly associated with SCZ and several haplotypes in this block have about twofold to threefold increase in cases. Our results provide further evidence that PDE4B may play important roles in the etiology of SCZ.

Analysis of new lactotransferrin gene variants in a case-control study related to periodontal disease in dog

The molecular and genetic research has contributed to a better understanding of the periodontal disease (PD) in humans and has shown that many genes play a role in the predisposition and progression of this complex disease. Variations in human lactotransferrin (LTF) gene appear to affect anti-microbial functions of this molecule, influencing the PD susceptibility. PD is also a major health problem in small animal practice, being the most common inflammatory disease found in dogs. Nevertheless, the research in genetic predisposition to PD is an unexplored subject in this species. This work aims to contribute to the characterization of the genetic basis of canine PD. In order to identify genetic variations and verify its association with PD, was performed a molecular analysis of LTF gene in a case-control approach, including 40 dogs in the PD cases group and 50 dogs in the control group. In this study were detected and characterized eight new single nucleotide variations in the dog LTF gene. Genotype and allele frequencies of these variations showed no statistically significant differences between the control and PD cases groups. Our data do not give evidence for the contribution of these LTF variations to the genetic background of canine PD. Nevertheless, the sequence variant L/15_g.411C > T leads to an aminoacid change (Proline to Leucine) and was predicted to be possibly damaging to the LTF protein. Further investigations would be of extreme value to clarify the biological importance of these new findings.

Assessing the contribution family data can make to case-control studies of rare variants

When pathogenic variants are rare then even among cases the proportion of subjects possessing a variant might be low, meaning that very large samples might be required to conclusively demonstrate evidence of an effect. Relatives of subjects within a case-control sample might provide useful additional information. The method of model-free linkage analysis implemented in MFLINK was adapted to incorporate linkage disequilibrium (LD) parameters in order to test for an effect of a putative pathogenic variant in complete LD with a disease locus. The effect of adding in to the analysis relatives of cases and controls found to carry the variant was investigated. When affected siblings or cousins of cases possessing the variant were incorporated they had a large effect on the results obtained. The evidence for involvement increased or reduced as expected, depending on whether or not the relatives themselves were found to possess the variant. The size of the effect was large relative to that expected from just increasing the size of a standard case-control sample. Affected relatives offer a valuable resource to assist the interpretation of case-control studies of rare variants. The method is capable of including other relative types and can deal with complex pedigrees.

Genetic association and sequencing of the insulin-like growth factor 1 gene in bipolar affective disorder

Insulin-like growth factor 1 (IGF1) has been shown to have an important role in brain development and function. Studies of IGF1 administration in rodents have shown that it has an anxiolytic and antidepressant effect. A genome-wide association study (GWAS) of the first University College London (UCL) cohort of 506 bipolar affective disorder subjects and 510 controls was carried out. The exons and flanking regions of IGF1 were resequenced, any new polymorphisms found were genotyped in an enlarged UCL sample of 937 cases and 941 controls. GWAS data gave good evidence of allelic and haplotypic association between multiple IGF1 SNP's and bipolar disorder (BD). New polymorphisms were found by resequencing IGF1 region. Data from GWAS and the new markers showed that twelve out of 43 SNPs showed association with BD with the four most significant SNPs having values of 3.7 × 10(-5) , 8.4 × 10(-4) , 2.6 × 10(-4) , and 2.5 × 10(-4) . A 5' promoter microsatellite polymorphism previously correlated with plasma lipoprotein concentration was also associated with BD (P = 0.013). Haplotypic association confirmed association with BD with significance values similar to the single marker SNP values. The marker rs12426318 has also been found to be associated with BD in a second sample. A test of gene wide significance with permutation testing for all markers genotyped at IGF1 was also significant. These data implicate IGF1 as a candidate gene to cause genetic susceptibility to BD.

A threonine to isoleucine missense mutation in the pericentriolar material 1 gene is strongly associated with schizophrenia

Markers at the pericentriolar material 1 gene (PCM1) have shown genetic association with schizophrenia in both a University College London (UCL) and a USA-based case-control sample. In this paper we report a statistically significant replication of the PCM1 association in a large Scottish case-control sample from Aberdeen. Resequencing of the genomic DNA from research volunteers who had inherited haplotypes associated with schizophrenia showed a threonine to isoleucine missense mutation in exon 24 which was likely to change the structure and function of PCM1 (rs370429). This mutation was found only as a heterozygote in 98 schizophrenic research subjects and controls out of 2246 case and control research subjects. Among the 98 carriers of rs370429, 67 were affected with schizophrenia. The same alleles and haplotypes were associated with schizophrenia in both the London and Aberdeen samples. Another potential aetiological base pair change in PCM1 was rs445422, which altered a splice site signal. A further mutation, rs208747, was shown by electrophoretic mobility shift assays to create or destroy a promoter transcription factor site. Five further non-synonymous changes in exons were also found. Genotyping of the new variants discovered in the UCL case-control sample strengthened the evidence for allelic and haplotypic association (P=0.02-0.0002). Given the number and identity of the haplotypes associated with schizophrenia, further aetiological base pair changes must exist within and around the PCM1 gene. PCM1 protein has been shown to interact directly with the disrupted-in-schizophrenia 1 (DISC1) protein, Bardet-Biedl syndrome 4, and Huntingtin-associated protein 1, and is important in neuronal cell growth. In a separate study we found that clozapine but not haloperidol downregulated PCM1 expression in the mouse brain. We hypothesize that mutant PCM1 may be responsible for causing a subtype of schizophrenia through abnormal cell division and abnormal regeneration in dividing cells in the central nervous system. This is supported by our previous finding of orbitofrontal volumetric deficits in PCM1-associated schizophrenia patients as opposed to temporal pole deficits in non-PCM1-associated schizophrenia patients. Caution needs to be exercised in interpreting the actual biological effects of the mutations we have found without further cell biology. However, the DNA changes we have found deserve widespread genotyping in multiple case-control populations.

High-density SNP screen of sodium channel genes by haplotype tagging and DNA pooling for association with idiopathic generalized epilepsy

We have investigated seven voltage-gated sodium channel genes for association with idiopathic generalized epilepsy (IGE). Probands and control DNA were grouped into pools and used to screen 85 single-nucleotide polymorphisms (SNPs), mostly HapMap SNPs tagging the common variation in these genes. Twelve SNPs exhibiting an allele frequency difference between pools were genotyped individually in our sample of 232 probands, 313 controls, and 95 parent-proband trios. Two SNPs, in SCN1A and SCN8A, were associated by allele and genotype at nominal level of significance, but were not significant after Bonferroni correction. Two SCN2A SNPs (rs3943809 and rs16850331) were associated by case-control with a subgroup with IGE and history of febrile seizures and also by transmission disequilibrium test (TDT) in parent-proband trios. Both SNPs are part of a linkage disequilibrium (LD) cluster of 38 SNPs, but none are obvious functional variants. The association of rs3943809 with the febrile seizure subgroup (p = 0.0004) remains significant after the conservative Bonferroni correction for multiple testing.

Evidence for the association of the DAOA (G72) gene with schizophrenia and bipolar disorder but not for the association of the DAO gene with schizophrenia

Background

Previous linkage and association studies have implicated the D-amino acid oxidase activator gene (DAOA)/G30 locus or neighbouring region of chromosome 13q33.2 in the genetic susceptibility to both schizophrenia and bipolar disorder. Four single nucleotide polymorphisms (SNPs) within the D-amino acid oxidase (DAO) gene located at 12q24.11 have also been found to show allelic association with schizophrenia.

Methods

We used the case control method to test for genetic association with variants at these loci in a sample of 431 patients with schizophrenia, 303 patients with bipolar disorder and 442 ancestrally matched supernormal controls all selected from the UK population.

Results

Ten SNPs spanning the DAOA locus were genotyped in these samples. In addition three SNPs were genotyped at the DAO locus in the schizophrenia sample. Allelic association was detected between the marker rs3918342 (M23), 3' to the DAOA gene and both schizophrenia (χ² = 5.824 p = 0.016) and bipolar disorder (χ² = 4.293 p = 0.038). A trend towards association with schizophrenia was observed for two other DAOA markers rs3916967 (M14, χ² = 3.675 p = 0.055) and rs1421292 (M24; χ² = 3.499 p = 0.062). A test of association between a three marker haplotype comprising of the SNPs rs778293 (M22), rs3918342 (M23) and rs1421292 (M24) and schizophrenia gave a global empirical significance of p = 0.015. No evidence was found to confirm the association of genetic markers at the DAO gene with schizophrenia.

Conclusion

Our results provide some support for a role for DAOA in susceptibility to schizophrenia and bipolar disorder.

Case-control studies show that a non-conservative amino-acid change from a glutamine to arginine in the P2RX7 purinergic receptor protein is associated with both bipolar- and unipolar-affective disorders

Three linkage studies of bipolar disorder have implicated chromosome 12q24.3 with lod scores of over 3.0 and several other linkage studies have found lods between 2 and 3. Fine mapping within the original chromosomal linkage regions has identified several loci that show association with bipolar disorder. One of these is the P2RX7 gene encoding a central nervous system-expressed purinergic receptor. A non-synonymous single nucleotide polymorphism, rs2230912 (P2RX7-E13A, G allele) and a microsatellite marker NBG6 were both previously found to be associated with bipolar disorder (P=0.00071 and 0.008, respectively). rs2230912 has also been found to show association with unipolar depression. The effect of the polymorphism is non-conservative and results in a glutamine to arginine change (Gln460Arg), which is likely to affect P2RX7 dimerization and protein-protein interactions. We have confirmed the allelic associations between bipolar disorder and the markers rs2230912 (P2RX7-E13A, G allele, P=0.043) and NBG6 (P=0.010) in a London-based sample of 604 bipolar cases and 560 controls. When we combined these data with the published case-control studies of P2RX7 and mood disorder (3586 individuals) the association between rs2230912 (Gln460Arg) and affective disorders became more robust (P=0.002). The increase in Gln460Arg was confined to heterozygotes rather than homozygotes suggesting a dominant effect (odds ratio 1.302, CI=1.129-1.503). Although further research is needed to prove that the Gln460Arg change has an aetiological role, it is so far the most convincing mutation to have been found with a role for increasing susceptibility to bipolar and genetically related unipolar disorders.

A simple method for assessing the strength of evidence for association at the level of the whole gene

Introduction

It is expected that different markers may show different patterns of association with different pathogenic variants within a given gene. It would be helpful to combine the evidence implicating association at the level of the whole gene rather than just for individual markers or haplotypes. Doing this is complicated by the fact that different markers do not represent independent sources of information.

Method

We propose combining the p values from all single locus and/or multilocus analyses of different markers according to the formula of Fisher, X = ∑(−2ln(p_i)), and then assessing the empirical significance of this statistic using permutation testing. We present an example application to 19 markers around the HTRA2 gene in a case-control study of Parkinson’s disease.

Results

Applying our approach shows that, although some individual tests produce low p values, overall association at the level of the gene is not supported.

Discussion

Approaches such as this should be more widely used in assimilating the overall evidence supporting involvement of a gene in a particular disease. Information can be combined from biallelic and multiallelic markers and from single markers along with multimarker analyses. Single genes can be tested or results from groups of genes involved in the same pathway could be combined in order to test biologically relevant hypotheses. The approach has been implemented in a computer program called COMBASSOC which is made available for downloading.

CLUMPHAP: a simple tool for performing haplotype-based association analysis

The completion of the HapMap Project and the development of high-throughput single nucleotide polymorphism genotyping technologies have greatly enhanced the prospects of identifying and characterizing the genetic variants that influence complex traits. In principle, association analysis of haplotypes rather than single nucleotide polymorphisms may better capture an underlying causal variant, but the multiple haplotypes can lead to reduced statistical power due to the testing of (and need to correct for) a large number of haplotypes. This paper presents a novel method based on clustering similar haplotypes to address this issue. The method, implemented in the CLUMPHAP program, is an extension of the CLUMP program designed for the analysis of multi-allelic markers (Sham and Curtis [1995] Ann. Hum. Genet. 59(Pt1):97-105). CLUMPHAP performs a hierarchical clustering of the haplotypes and then computes the chi(2) statistic between each haplotype cluster and disease; the statistical significance of the largest of the chi(2) statistics is obtained by permutation testing. A significant result suggests that the presence of a disease-causing variant in the haplotype cluster is over-represented in cases. Using simulation studies, we have compared CLUMPHAP and more widely used approaches in terms of their statistical power to identify an untyped susceptibility locus. Our results show that CLUMPHAP tends to have greater power than the omnibus haplotype test and is comparable in power to multiple regression locus-coding approaches.

Interleukin-6 polymorphisms are associated with pathogenic bacteria in subjects with periodontitis

BACKGROUND

Growing evidence suggests that individual genetic susceptibility may influence the host's response to infections. Previously, we showed that a common variation in the interleukin (IL)-6 gene was associated with increased odds of detection of common periodontal pathogens from individuals with aggressive periodontitis. The aim of this study was to investigate the association between IL-6 polymorphisms and periodontopathogenic bacteria in a larger, ethnically mixed population of subjects with periodontitis.

METHODS

Genomic DNA was extracted from 107 subjects diagnosed with severe forms of periodontitis to study a cluster of polymorphisms in inflammatory genes, including IL-6. The presence of Aggregatibacter actinomycetemcomitans (previously Actinobacillus actinomycetemcomitans), Porphyromonas gingivalis, and Tannerella forsythia (previously T. forsythensis) in their subgingival biofilm was determined by polymerase chain reaction analysis. Serum IL-6 and C-reactive protein (CRP) concentrations were analyzed by enzyme-linked immunosorbent assay.

RESULTS

Multiple logistic regression analysis revealed that the IL-6 -6106 polymorphism was associated with the detection of A. actinomycetemcomitans (P = 0.009; odds ratio [OR] = 3.5; 95% confidence interval [CI]: 1.38 to 9.16) and the concomitant detection of A. actinomycetemcomitans and P. gingivalis (P = 0.015; OR = 3.6; 95% CI: 1.28 to 10.04). The IL-6 -174 polymorphism was associated with increased odds of the concomitant detection of A. actinomycetemcomitans and P. gingivalis (P = 0.042; OR = 2.8; 95% CI: 1.04 to 7.75). Haplotype analysis of all five IL-6 polymorphisms confirmed an association with the detection of A. actinomycetemcomitans (P = 0.046). The IL-6 -6106 polymorphism was also associated with CRP serum levels at multivariate analysis (P = 0.024).

CONCLUSIONS

These findings confirm the hypothesis that complex interactions between the microbiota and host genome are at the basis of susceptibility to periodontitis. Periodontal disease may represent a useful model to study the pathways and mechanisms of host-pathogen interactions in inflammatory diseases.

Investigation into the ability of SNP chipsets and microsatellites to detect association with a disease locus

We wished to investigate the ability of different SNP chipsets to detect association with a disease and to investigate the linkage disequilibrium (LD) relationships between microsatellites and nearby SNPs in order to assess their potential usefulness to detect association.

SNP genotypes were obtained from HapMap and microsatellite genotypes from CEPH. 5000 SNPs were simulated as disease genes which increased penetrance from 0.01 to 0.02 in a sample of 400 cases and 400 controls. The power of flanking SNPs to detect association was tested using sets of 1, 2, 3 or 4 markers analysed with haplotype analysis or logistic regression and using either all HapMap markers or those from the Affymetrix 500K, Illumina 300K or Illumina 550K chipsets. Additionally, LD relationships between 10 microsatellites and SNPs within 2Mb of each other were studied.

The power for one of the markers to detect association at p = 0.001 was around 0.4. Power was slightly better for logistic regression than haplotype analysis and for two-marker as opposed to single marker analysis but analysing with larger numbers markers had little benefit. The Illumina 550K marker set was better able to detect association than the other two and was almost as powerful as using all HapMap markers. Microsatellites had detectable LD with only a small number of nearby SNPs and the pattern of LD was very variable.

Available chipsets have quite good ability to detect association although obviously results will be critically dependent on the nature of the genetic effect on risk, sample size and the actual LD relationships of the susceptibility polymorphisms involved. Microsatellites seem ill-suited for systematic studies to detect association.

Association between interleukin-6 promoter haplotypes and aggressive periodontitis

BACKGROUND

Interleukin-6 (IL-6) polymorphisms have been shown to affect IL-6 promoter activity. This study investigated the possible role of IL-6 genetic polymorphisms and haplotypes in the predisposition to aggressive periodontitis (AgP).

MATERIAL AND METHODS

A case-control association study on 224 AgP patients and 231 healthy controls was performed in order to detect differences in genotype distributions of five single nucleotide polymorphisms (SNPs) located in the promoter region of the IL-6 gene.

RESULTS

The IL-6 -1363 polymorphism was associated with a diagnosis of AgP in subjects of all ethnicities (p=0.006, adjusted logistic regression). The -1480 SNP was associated with LAgP in subjects of all ethnicities (p=0.003). The -1480 and -6106 polymorphisms were associated with Localized AgP in Caucasians (n=24) (p=0.007 and 0.010, respectively). Haplotypes determined by the -1363 and -1480 polymorphisms were also associated with LAgP (p=0.001) in Caucasians.

CONCLUSIONS

This study supports the hypothesis of a link between IL-6 genetic factors and AgP and highlights the importance of two IL-6 polymorphisms (-1363 and -1480) in modulating disease phenotype and susceptibility.

An investigation of the alpha1A-adrenergic receptor gene and antipsychotic-induced side-effects

Antipsychotic treatment is hampered by the induction of side-effects such as tardive dyskinesia (TD), weight gain, sedation and extrapyramidal side-effects (EPS). Identification of the factors related to their development would facilitate their avoidance and the improvement of antipsychotic treatment. It has been hypothesised that genetic variants in drug targeted receptors may contribute to the development of side-effects. In this study, we have investigated the possible influence of genetic variants (-563-C/T, -4155-G/C and -4884-A/G) of the alpha(1A)-adrenergic receptor, an important target of atypical antipsychotic drugs, and development of side-effects after antipsychotic medication in a sample of N = 427 US Caucasian patients. We found several marginal associations (p < 0.05) between alpha(1A)-adrenergic genetic variants and antipsychotic-induced side-effects which did not reach statistical significance after corrections for multiple analyses. These results do not support a major role of alpha(1A)-adrenergic genetic variants in obesity and other side-effects observed after prolonged treatment with antipsychotic medications.

Comparison of artificial neural network analysis with other multimarker methods for detecting genetic association

BACKGROUND

Debate remains as to the optimal method for utilising genotype data obtained from multiple markers in case-control association studies. I and colleagues have previously described a method of association analysis using artificial neural networks (ANNs), whose performance compared favourably to single-marker methods. Here, the performance of ANN analysis is compared with other multi-marker methods, comprising different haplotype-based analyses and locus-based analyses.

RESULTS

Of several methods studied and applied to simulated SNP datasets, heterogeneity testing of estimated haplotype frequencies using asymptotic p values rather than permutation testing had the lowest power of the methods studied and ANN analysis had the highest power. The difference in power to detect association between these two methods was statistically significant (p = 0.001) but other comparisons between methods were not significant. The raw t statistic obtained from ANN analysis correlated highly with the empirical statistical significance obtained from permutation testing of the ANN results and with the p value obtained from the heterogeneity test.

CONCLUSION

Although ANN analysis was more powerful than the standard haplotype-based test it is unlikely to be taken up widely. The permutation testing necessary to obtain a valid p value makes it slow to perform and it is not underpinned by a theoretical model relating marker genotypes to disease phenotype. Nevertheless, the superior performance of this method does imply that the widely-used haplotype-based methods for detecting association with multiple markers are not optimal and efforts could be made to improve upon them. The fact that the t statistic obtained from ANN analysis is highly correlated with the statistical significance does suggest a possibility to use ANN analysis in situations where large numbers of markers have been genotyped, since the t value could be used as a proxy for the p value in preliminary analyses.

Association study of dysbindin gene with clinical and outcome measures in a representative cohort of Italian schizophrenic patients

There is evidence suggesting that Dysbindin (DTNBP1) is a susceptibility gene for schizophrenia in Caucasian, Chinese, and Japanese populations. We sought to determine if dysbindin was associated with schizophrenia and its symptoms in a representative group of schizophrenic patients from a Community-Based Mental Health Service (CMHS) in Verona, Italy. A prevalence cohort of schizophrenic patients (n = 141) was assessed at baseline and then 3 and 6 years later. Eighty patients and 106 healthy controls were genotyped for polymorphisms in dysbindin. We tested if diagnosis, clinical symptoms as measured by the Brief Psychiatric Rating Scale (BPRS), and functioning as measured by the Global Assessment of Functioning Scale (GAF), were associated with the presence of certain dysbindin polymorphisms. Finally, using the longitudinal clinical data, we tested if patients carrying dysbindin high-risk haplotypes had a more unfavorable longitudinal clinical outcome. A trend towards statistical association (P = 0.058) between schizophrenia and rs2619538 was found. Using GENECOUNTING software, we found that rs2619538-P1583 (P = 0.048), P1320-P1757 (P = 0.034), and rs2619538-P1583-P1578 (P = 0.040) haplotypes occurred more often in cases compared to controls before correction for multiple testing. The rs2619538-P1583 haplotype was more likely to be transmitted to subjects with more severe and persistent psychopathology. These preliminary results are compatible with the view that DTNBP1 is a susceptibility factor for schizophrenia, and is associated with worse psychopathology.

Minor differences in haplotype frequency estimates can produce very large differences in heterogeneity test statistics

BACKGROUND

Tests for association between a haplotype and disease are commonly performed using a likelihood ratio test for heterogeneity between case and control haplotype frequencies. Using data from a study of association between heroin dependence and the DRD2 gene, we obtained estimated haplotype frequencies and the associated likelihood ratio statistic using two different computer programs, MLOCUS and GENECOUNTING. We also carried out permutation testing to assess the empirical significance of the results obtained.

RESULTS

Both programs yielded similar, though not identical, estimates for the haplotype frequencies. MLOCUS produced a p value of 1.8*10-15 and GENECOUNTING produced a p value of 5.4*10-4. Permutation testing produced a p value 2.8*10-4.

CONCLUSION

The fact that very large differences occur between the likelihood ratio statistics from the two programs may reflect the fact that the haplotype frequencies for the combined group are not constrained to be equal to the weighted averages of the frequencies for the cases and controls, as they would be if they were directly observed rather than being estimated. Minor differences in haplotype frequency estimates can result in very large differences in the likelihood ratio statistic and associated p value.

Allelic association studies of genome wide association data can reveal errors in marker position assignments

BACKGROUND

Genome wide association (GWA) studies provide the opportunity to develop new kinds of analysis. Analysing pairs of markers from separate regions might lead to the detection of allelic association which might indicate an interaction between nearby genes.

METHODS

396,591 markers typed in 541 subjects were studied. 7.8*1010 pairs of markers were screened and those showing initial evidence for allelic association were subjected to more thorough investigation along with 10 flanking markers on either side.

RESULTS

No evidence was detected for interaction. However 6 markers appeared to have an incorrect map position according to NCBI Build 35. One of these was corrected in Build 36 and 2 were dropped. The remaining 3 were left with map positions inconsistent with their allelic association relationships.

DISCUSSION

Although no interaction effects were detected the method was successful in identifying markers with probably incorrect map positions.

CONCLUSION

The study of allelic association can supplement other methods for assigning markers to particular map positions. Analyses of this type may usefully be applied to data from future GWA studies.

Fine mapping by genetic association implicates the chromosome 1q23.3 gene UHMK1, encoding a serine/threonine protein kinase, as a novel schizophrenia susceptibility gene

BACKGROUND

Linkage studies by us and others have confirmed that chromosome 1q23.3 is a susceptibility locus for schizophrenia. Based on this information, several research groups have published evidence that markers within both the RGS4 and CAPON genes, which are 700 kb apart, independently showed allelic association with schizophrenia. Tests of allelic association with both of these genes in our case control sample were negative. Therefore, we carried out further fine mapping between the RGS4 and CAPON genes.

METHODS

Twenty-nine SNP and microsatellite markers in the 1q23.3 region were genotyped in the United Kingdom based sample of 450 cases and 450 supernormal control subjects.

RESULTS

We detected positive allelic association after the eighth marker was genotyped and found that three microsatellite markers (p = .011, p = .014, p = .049) and two SNPs (p = .004, p = .043) localized in the 700 kb region between the RGS4 and CAPON genes, within the UHMK1 gene, were associated with schizophrenia. Tests of significance for marker rs10494370 remained significant following Bonferroni correction (alpha = .006) for multiple tests. Tests of haplotypic association were also significant for UHMK1 (p = .009) using empirical permutation tests, which make it unnecessary to further correct for both multiple alleles and multiple markers.

CONCLUSIONS

These results provide preliminary evidence that the UHMK1 gene increases susceptibility to schizophrenia. Further confirmation in adequately powered samples is needed. UHMK1 is a serine threonine kinase nuclear protein and is highly expressed in regions of the brain implicated in schizophrenia.

A genetic association study of chromosome 11q22-24 in two different samples implicates the FXYD6 gene, encoding phosphohippolin, in susceptibility to schizophrenia

Previous linkage analyses of families with multiple cases of schizophrenia by us and others have confirmed the involvement of the chromosome 11q22-24 region in the etiology of schizophrenia, with LOD scores of 3.4 and 3.1. We now report fine mapping of a susceptibility gene in the 11q22-24 region, determined on the basis of a University College London (UCL) sample of 496 cases and 488 supernormal controls. Confirmation was then performed by the study of an Aberdeen sample consisting of 858 cases and 591 controls (for a total of 2,433 individuals: 1,354 with schizophrenia and 1,079 controls). Seven microsatellite or single-nucleotide polymorphism (SNP) markers localized within or near the FXYD6 gene showed empirically significant allelic associations with schizophrenia in the UCL sample (for D11S1998, P=.021; for rs3168238, P=.009; for TTTC20.2, P=.048; for rs1815774, P=.049; for rs4938445, P=.010; for rs4938446, P=.025; for rs497768, P=.023). Several haplotypes were also found to be associated with schizophrenia; for example, haplotype Hap-F21 comprising markers rs10790212-rs4938445-rs497768 was found to be associated with schizophrenia, by a global permutation test (P=.002). Positive markers in the UCL sample were then genotyped in the Aberdeen sample. Two of these SNPs were found to be associated with schizophrenia in the Scottish sample (for rs4938445, P=.044; for rs497768, P=.037). The Hap-F21 haplotype also showed significant association with schizophrenia in the Aberdeen sample, with the same alleles being associated (P=.013). The FXYD6 gene encodes a protein called "phosphohippolin" that is highly expressed in regions of the brain thought to be involved in schizophrenia. The protein functions by modulating the kinetic properties of Na,K-ATPase to the specific physiological requirements of the tissue. Etiological base-pair changes in FXYD6 or in associated promoter/control regions are likely to cause abnormal function or expression of phosphohippolin and to increase genetic susceptibility to schizophrenia.

NADPH oxidase (CYBA) and FcgammaR polymorphisms as risk factors for aggressive periodontitis: a case-control association study

INTRODUCTION

Neutrophils (PMN) in aggressive periodontitis (AgP) patients have been reported to be hyperactive especially with regards to superoxide production. Polymorphisms in genes influencing PMN function have been proposed as candidate risk factors for AgP. The aim of this study was to test the association of specific gene polymorphisms affecting PMN functions with AgP.

MATERIALS AND METHODS

Two hundred and twenty-four patients with confirmed diagnosis of AgP and 231 subjects with healthy periodontium took part in the study. A blood sample was collected from subjects and genotypes for p22phox (CYBA) NADPH oxidase, FP, Fcalpha and Fcgamma receptors were analysed in a blind fashion.

RESULTS

The C242T p22phox NADPH oxidase T allele was significantly associated with AgP in a multiple logistic regression model adjusting for confounders, and this was observed for all subjects [p = 0.002, odds ratio (OR) = 1.87, 95% confidence interval (CI) = 1.27-2.83] and Caucasians (p = 0.009, OR=2.07, 95% CI = 1.20-3.59). Concomitant presence of C242T p22phox NADPH oxidase T allele and FcgammaRIIIb NA1 homozygosity was associated with the generalized AgP phenotype in Caucasians (p = 0.001, OR = 30.35, 95% CI = 3.81-241.97).

CONCLUSIONS

C242T p22phox NADPH oxidase and FcgammaR polymorphisms may predispose to AgP through a modulation of neutrophil superoxide production.

Fine mapping of a susceptibility locus for bipolar and genetically related unipolar affective disorders, to a region containing the C21ORF29 and TRPM2 genes on chromosome 21q22.3

Linkage analyses of bipolar families have confirmed that there is a susceptibility locus near the telomere on chromosome 21q. To fine map this locus we carried out tests of allelic association using 30 genetic markers near the telomere at 21q22.3 in 600 bipolar research subjects and 450 ancestrally matched supernormal control subjects. We found significant allelic association with the microsatellite markers D21S171 (P=0.016) and two closely linked single-nucleotide polymorphisms, rs1556314 (P=0.008) and rs1785467 (P=0.025). A test of association with a three locus haplotype across the susceptibility region was significant with a permutation test of P=0.011. A two SNP haplotype was also significantly associated with bipolar disorder (P=0.01). Only two brain expressed genes, TRPM2 and C21ORF29 (TSPEAR), are present in the associated region. TRPM2 encodes a calcium channel receptor and TSPEAR encodes a peptide with repeats associated with epilepsy in the mouse. DNA from subjects who had inherited the associated marker alleles was sequenced. A base pair change (rs1556314) in exon 11 of TRPM2, which caused a change from an aspartic acid to a glutamic acid at peptide position 543 was found. This SNP showed the strongest association with bipolar disorder (P=0.008). Deletion of exon 11 of TRPM2 is known to cause dysregulation of cellular calcium homeostasis in response to oxidative stress. A second nonconservative change from arginine to cysteine at position 755 in TRPM2 (ss48297761) was also detected. A third nonconservative change from histidine to glutamic acid was found in exon 8 of TSPEAR. These changes need further investigation to establish any aetiological role in bipolar disorder.