Modeling the cumulative genetic risk for multiple sclerosis from genome-wide association data

Background: Multiple sclerosis (MS) is the most common cause of chronic neurologic disability beginning in early to middle adult life. Results from recent genome-wide association studies (GWAS) have substantially lengthened the list of disease loci and provide convincing evidence supporting a multifactorial and polygenic model of inheritance. Nevertheless, the knowledge of MS genetics remains incomplete, with many risk alleles still to be revealed. Methods: We used a discovery GWAS dataset (8,844 samples, 2,124 cases and 6,720 controls) and a multi-step logistic regression protocol to identify novel genetic associations. The emerging genetic profile included 350 independent markers and was used to calculate and estimate the cumulative genetic risk in an independent validation dataset (3,606 samples). Analysis of covariance (ANCOVA) was implemented to compare clinical characteristics of individuals with various degrees of genetic risk. Gene ontology and pathway enrichment analysis was done using the DAVID functional annotation tool, the GO Tree Machine, and the Pathway-Express profiling tool. Results: In the discovery dataset, the median cumulative genetic risk (P-Hat) was 0.903 and 0.007 in the case and control groups, respectively, together with 79.9% classification sensitivity and 95.8% specificity. The identified profile shows a significant enrichment of genes involved in the immune response, cell adhesion, cell communication/ signaling, nervous system development, and neuronal signaling, including ionotropic glutamate receptors, which have been implicated in the pathological mechanism driving neurodegeneration. In the validation dataset, the median cumulative genetic risk was 0.59 and 0.32 in the case and control groups, respectively, with classification sensitivity 62.3% and specificity 75.9%. No differences in disease progression or T2-lesion volumes were observed among four levels of predicted genetic risk groups (high, medium, low, misclassified). On the other hand, a significant difference (F = 2.75, P = 0.04) was detected for age of disease onset between the affected misclassified as controls (mean = 36 years) and the other three groups (high, 33.5 years; medium, 33.4 years; low, 33.1 years). Conclusions: The results are consistent with the polygenic model of inheritance. The cumulative genetic risk established using currently available genome-wide association data provides important insights into disease heterogeneity and completeness of current knowledge in MS genetics.

A polymorphism in the HLA-DPB1 gene is associated with susceptibility to multiple sclerosis

We conducted an association study across the human leukocyte antigen (HLA) complex to identify loci associated with multiple sclerosis (MS). Comparing 1927 SNPs in 1618 MS cases and 3413 controls of European ancestry, we identified seven SNPs that were independently associated with MS conditional on the others (each ). All associations were significant in an independent replication cohort of 2212 cases and 2251 controls () and were highly significant in the combined dataset (). The associated SNPs included proxies for HLA-DRB1*15:01 and HLA-DRB1*03:01, and SNPs in moderate linkage disequilibrium (LD) with HLA-A*02:01, HLA-DRB1*04:01 and HLA-DRB1*13:03. We also found a strong association with rs9277535 in the class II gene HLA-DPB1 (discovery set , replication set , combined ). HLA-DPB1 is located centromeric of the more commonly typed class II genes HLA-DRB1, -DQA1 and -DQB1. It is separated from these genes by a recombination hotspot, and the association is not affected by conditioning on genotypes at DRB1, DQA1 and DQB1. Hence rs9277535 represents an independent MS-susceptibility locus of genome-wide significance. It is correlated with the HLA-DPB1*03:01 allele, which has been implicated previously in MS in smaller studies. Further genotyping in large datasets is required to confirm and resolve this association.

Investigation of the [−/A]8 and C1236T genetic variations within the human toll-like receptor 3 gene for association with multiple sclerosis

Multiple sclerosis (MS) is a serious cause of neurological disability among young adults. The clinical course remains difficult to predict, and the pathogenesis of the disease is still modestly understood. Autoimmunity is thought to be a key aspect of the disease, with autoreactive T cells thought to mediate central nervous system (CNS) inflammation to some extent. Toll-like receptors are known to mediate cellular recognition of pathogens by way of patterns of molecular presentation. Toll-like receptor 3 is coded by the gene TLR3 and is recognized as an important factor in virus recognition and is known to be involved in the expression of neuroprotective mediators. We set out to investigate two variations within the TLR3 gene, an 8 bp insertion-deletion \[-/A](8) and a single base-pair variation C1236T, in subjects with MS and matched healthy controls to determine whether significant differences exist in these markers in an Australian population. We used capillary gel electrophoresis and TaqMan genotyping assay techniques to resolve genotypes for each marker, respectively. Our work found no significant difference between frequencies for TLR3 \[-/A](8) by genotype (chi(2)=1.03, p=0.60) or allele (chi(2)=1.09, p=0.30). Similarly, we found no evidence for the association of TLR3 C1236T by genotype (chi(2)=0.35, p=0.84) or allele frequency (chi(2)=0.31, p=0.58). This work reveals no evidence to suggest that these markers are associated with MS in the tested population. Although the role of TLR3 and the wider toll-like receptor family remain significant in neurological and CNS inflammatory disorders, our current work does not support a role for the two tested variants in this gene with regard to MS susceptibility.

Multiple sclerosis susceptibility-associated SNPs do not influence disease severity measures in a cohort of Australian MS patients

Recent association studies in multiple sclerosis (MS) have identified and replicated several single nucleotide polymorphism (SNP) susceptibility loci including CLEC16A, IL2RA, IL7R, RPL5, CD58, CD40 and chromosome 12q13–14 in addition to the well established allele HLA-DR15. There is potential that these genetic susceptibility factors could also modulate MS disease severity, as demonstrated previously for the MS risk allele HLA-DR15. We investigated this hypothesis in a cohort of 1006 well characterised MS patients from South-Eastern Australia. We tested the MS-associated SNPs for association with five measures of disease severity incorporating disability, age of onset, cognition and brain atrophy. We observed trends towards association between the RPL5 risk SNP and time between first demyelinating event and relapse, and between the CD40 risk SNP and symbol digit test score. No associations were significant after correction for multiple testing. We found no evidence for the hypothesis that these new MS disease risk-associated SNPs influence disease severity.

An investigation of the C77G and C772T variations within the human protein tyrosine phosphatase receptor type C gene for association with multiple sclerosis in an Australian population

Multiple sclerosis (MS) is a common cause of neurological disability in young adults. The disease generally manifests in early to middle adulthood and causes various neurological deficits. Autoreactive T lymphocytes and their associated antigens have long been presumed important features of MS pathogenesis. The Protein tyrosine phosphatase receptor type C gene (PTPRC) encodes the T-cell receptor CD45. Variations within PTPRC have been previously associated with diseases of autoimmune origin such as type 1 diabetes mellitus and Graves' disease. We set out to investigate two variants within the PTPRC gene, C77G and C772T in subjects with MS and matched healthy controls to determine whether significant differences exist in these markers in an Australian population. We employed high resolution melt analysis (HRM) and restriction length polymorphism (RFLP) techniques to determine genotypic and allelic frequencies. Our study found no significant difference between frequencies for PTPRC C77G by either genotype (Χ2 = 0.65, P = 0.72) or allele (Χ2 = 0.48, P = 0.49). Similarly, we did not find evidence to suggest an association between PTPRC C772T by genotype (Χ2 = 1.06, P = 0.59) or allele (Χ2 = 0.20, P = 0.66). Linkage disequilibrium (LD) analysis showed strong linkage disequilibrium between the two tested markers (D' = 0.9970, SD = 0.0385). This study reveals no evidence to suggest that these markers are associated with MS in the tested Australian Caucasian population. Although the PTPRC gene has a significant role in regulating CD4+ and CD8+ autoreactive T-cells, interferon-beta responsiveness, and potentially other important processes, our study does not support a role for the two tested variants of this gene in MS susceptibility in the Australian population.

Genome-wide association study identifies new multiple sclerosis susceptibility loci on chromosomes 12 and 20

To identify multiple sclerosis (MS) susceptibility loci, we conducted a genome-wide association study (GWAS) in 1,618 cases and used shared data for 3,413 controls. We performed replication in an independent set of 2,256 cases and 2,310 controls, for a total of 3,874 cases and 5,723 controls. We identified risk-associated SNPs on chromosome 12q13-14 (rs703842, P = 5.4 x 10(-11); rs10876994, P = 2.7 x 10(-10); rs12368653, P = 1.0 x 10(-7)) and upstream of CD40 on chromosome 20q13 (rs6074022, P = 1.3 x 10(-7); rs1569723, P = 2.9 x 10(-7)). Both loci are also associated with other autoimmune diseases. We also replicated several known MS associations (HLA-DR15, P = 7.0 x 10(-184); CD58, P = 9.6 x 10(-8); EVI5-RPL5, P = 2.5 x 10(-6); IL2RA, P = 7.4 x 10(-6); CLEC16A, P = 1.1 x 10(-4); IL7R, P = 1.3 x 10(-3); TYK2, P = 3.5 x 10(-3)) and observed a statistical interaction between SNPs in EVI5-RPL5 and HLA-DR15 (P = 0.001).

Results of a phase IIa clinical trial of an anti-inflammatory molecule, chaperonin 10, in multiple sclerosis

Background Chaperonin 10 (Cpn10) is a mitochondrial molecule involved in protein folding. The aim of this study was to determine the safety profile of Cpn10 in patients with multiple sclerosis (MS). Methods A total of 50 patients with relapse-remitting or secondary progressive MS were intravenously administered 5 mg or 10 mg of Cpn10 weekly for 12 weeks in a double-blind, randomized, placebo controlled, phase II trial. Clinical reviews, including Expanded Disability Status Scale and magnetic resonance imaging (MRI) with Gadolinium, were undertaken every 4 weeks. Stimulation of patient peripheral blood mononuclear cells with lipopolysaccharide ex vivo was used to measure the in vivo activity of Cpn10. Results No significant differences in the frequency of adverse events were seen between treatment and placebo arms. Leukocytes from both groups of Cpn10-treated patients produced significantly lower levels of critical proinflammatory cytokines. A trend toward improvement in new Gadolinium enhancing lesions on MRI was observed, but this difference was not statistically significant. No differences in clinical outcome measures were seen. Conclusions Cpn10 is safe and well tolerated when administered to patients with MS for 3 months, however, a further extended phase II study primarily focused on efficacy is warranted.

Study of leukemia inhibitory factor polymorphism within an Australian multiple sclerosis population

OBJECTIVE: To examine a polymorphism within the 3' untranslated region of the leukemia inhibitory factor gene for an association with multiple sclerosis within an Australian case-control population. METHODS: A test group of 121 unrelated multiple sclerosis patients, of Caucasian origin, and 121 controls, matched for ethnicity, sex and age (+/-5 years) were included in the study. The LIF 3' UTR StuI polymorphism was genotyped by restriction fragment length polymorphism analysis. Statistical analysis of genotype and allele frequencies included Hardy-Weinberg law and conventional contingency table analysis incorporating the standard chi-squared test for independence. RESULTS: Allelic and genotype frequencies did not demonstrate a significant association between the case and control groups for the tested LIF 3' UTR StuI polymorphism. CONCLUSION: The results indicate that the LIF 3' UTR StuI polymorphism is not associated with multiple sclerosis, however we cannot exclude the hypothesis that other polymorphic alleles of LIF could be implicated in MS susceptibility.

Allelic variation investigation of the estrogen receptor within an Australian multiple sclerosis population

Multiple Sclerosis (MS) is a central nervous system (CNS) chronic inflammatory demyelinating disease leading to various neurological disabilities. The disorder is more prevalent for women with a ratio of 3:2 female to male. Objectives: To investigate variation within the estrogen receptor 1 (ESR1) polymorphism gene in an Australian MS case-control population using two intragenic restriction fragment length polymorphisms; the G594A located in exon 8 detected with the BtgI restriction enzyme and T938C located in intron 1, detected with PvuII. One hundred and ten Australian MS patients were studied, with patients classified clinically as Relapsing Remitting MS (RR-MS), Secondary Progressive MS (SP-MS) or Primary Progressive MS (PP-MS). Also, 110 age, sex and ethnicity matched controls were investigated as a comparative group. No significant difference in the allelic distribution frequency was found between the case and control groups for the ESR1 PvuII (P = 0.50) and Btg1 (P = 0.45) marker. Our results do not support a role for these two ESR1 markers in multiple sclerosis susceptibility, however other markers within ESR1 should not be excluded for potential involvement in the disorder.

Gene expression studies in multiple sclerosis

Multiple sclerosis (MS) is a serious neurological disorder affecting young Caucasian individuals, usually with an age of onset at 18 to 40 years old. Females account for approximately 60x of MS cases and the manifestation and course of the disease is highly variable from patient to patient. The disorder is characterised by the development of plaques within the central nervous system (CNS). Many gene expression studies have been undertaken to look at the specific patterns of gene transcript levels in MS. Human tissues and experimental mice were used in these gene-profiling studies and a very valuable and interesting set of data has resulted from these various expression studies. In general, genes showing variable expression include mainly immunological and inflammatory genes, stress and antioxidant genes, as well as metabolic and central nervous system markers. Of particular interest are a number of genes localised to susceptible loci previously shown to be in linkage with MS. However due to the clinical complexity of the disease, the heterogeneity of the tissues used in expression studies, as well as the variable DNA chips/membranes used for the gene profiling, it is difficult to interpret the available information. Although this information is essential for the understanding of the pathogenesis of MS, it is difficult to decipher and define the gene pathways involved in the disorder. Experiments in gene expression profiling in MS have been numerous and lists of candidates are now available for analysis. Researchers have investigated gene expression in peripheral mononuclear white blood cells (PBMCs), in MS animal models Experimental Allergic Encephalomyelitis (EAE) and post mortem MS brain tissues. This review will focus on the results of these studies.

Genetic investigation of methylenetetrahydrofolate reductase (MTHFR) and catechol-O-methyl transferase (COMT) in multiple sclerosis

Multiple sclerosis (MS) is a chronic neurological disease characterized by central nervous system (CNS) inflammation and demyelination. The C677T substitution variant in the methylenetetrahydrofolate reductase (MTHFR) gene has been associated with increased levels of circulating homocysteine and is a mild risk factor for vascular disease. Higher blood levels of homocysteine have also been reported in MS. Thus, the C677T mutation of the MTHFR gene may influence MS susceptibility. Noradrenaline, a neurotransmitter believed to play an immunosupressive role in neuroinflammatory disorders, is catabolized by catechol-O-methyl transferase (COMT). The COMT G158A substitution results in a three- to four-fold decreased activity of the COMT enzyme, which may influence CNS synaptic catecholamine breakdown and could also play a role in MS inflammation. We tested DNA from Australian MS patients and unaffected control subjects, matched for gender, age and ethnicity. Specifically, we genotyped the MTHFR C677T and the COMT G158A mutations. Genotype distributions showed that the homozygous mutant MTHFR genotype (T/T) and the COMT (H/H) genotype were slightly over-represented in the MS group (16% versus 11% and 24% versus 19%, respectively), but both variations failed to reach statistical significance (P=0.15 and P=0.32, respectively). Hence, results from the present study do not support a major role for either functional gene mutation in MS susceptibility.

In vitro and in vivo MMP gene expression localisation by In Situ-RT-PCR in cell culture and paraffin embedded human breast cancer cell line xenografts

Background Members of the matrix metalloproteinase (MMP) family of proteases are required for the degradation of the basement membrane and extracellular matrix in both normal and pathological conditions. In vitro, MT1-MMP (MMP-14, membrane type-1-MMP) expression is higher in more invasive human breast cancer (HBC) cell lines, whilst in vivo its expression has been associated with the stroma surrounding breast tumours. MMP-1 (interstitial collagenase) has been associated with MDA-MB-231 invasion in vitro, while MMP-3 (stromelysin-1) has been localised around invasive cells of breast tumours in vivo. As MMPs are not stored intracellularly, the ability to localise their expression to their cells of origin is difficult. Methods We utilised the unique in situ-reverse transcription-polymerase chain reaction (IS-RT-PCR) methodology to localise the in vitro and in vivo gene expression of MT1-MMP, MMP-1 and MMP-3 in human breast cancer. In vitro, MMP induction was examined in the MDA-MB-231 and MCF-7 HBC cell lines following exposure to Concanavalin A (Con A). In vivo, we examined their expression in archival paraffin embedded xenografts derived from a range of HBC cell lines of varied invasive and metastatic potential. Mouse xenografts are heterogenous, containing neoplastic human parenchyma with mouse stroma and vasculature and provide a reproducible in vivo model system correlated to the human disease state. Results In vitro, exposure to Con A increased MT1-MMP gene expression in MDA-MB-231 cells and decreased MT1-MMP gene expression in MCF-7 cells. MMP-1 and MMP-3 gene expression remained unchanged in both cell lines. In vivo, stromal cells recruited into each xenograft demonstrated differences in localised levels of MMP gene expression. Specifically, MDA-MB-231, MDA-MB-435 and Hs578T HBC cell lines are able to influence MMP gene expression in the surrounding stroma. Conclusion We have demonstrated the applicability and sensitivity of IS-RT-PCR for the examination of MMP gene expression both in vitro and in vivo. Induction of MMP gene expression in both the epithelial tumour cells and surrounding stromal cells is associated with increased metastatic potential. Our data demonstrate the contribution of the stroma to epithelial MMP gene expression, and highlight the complexity of the role of MMPs in the stromal-epithelial interactions within breast carcinoma.

An examination of MS candidate genes identified as differentially regulated in multiple sclerosis plaque tissue, using absolute and comparative real-time Q-PCR analysis

In our laboratory, we have developed methods in real-time detection and quantitative-polymerase chain reaction (Q-PCR) to analyse the relative levels of gene expression in post mortem brain tissues. We have then applied this method to examine differences in gene activity between normal white matter (NWM) and plaque tissue from multiple sclerosis (MS) patients. Genes were selected based on their association with pathology and through identification by previously conducted global gene expression analysis. Plaque tissue was obtained from secondary progressive (SP) patients displaying chronic active, as well as acute pathologies; while NWM from the same location was obtained from age- and sex-matched controls (normal patients). In this study, we used both SYBR Green I supplementation and commercially available mixes to assess both comparative and absolute levels of gene activity. The results of both methods compared favourably for four of the five genes examined (P < 0.05, Pearsons), while differences in gene expression between chronic active and acute pathologies were also identified. For example, a >50-fold increase in osteopontin (Spp1) and inositol 1-4-5 phosphate 3 kinase B (Itpkb) levels in acute plaques contrasted with the 5-fold or less increase in chronic active plaques (P < 0.05, unpaired t test). By contrast, there was no significant difference in the levels of the MS marker and calcium-dependent protease (Calpain, Capns1) in MS plaque tissue. In summary, Q-PCR analysis using SYBR Green I has allowed us to economically obtain what may be clinically significant information from small amounts of the CNS, providing an opportunity for further clinical investigations.

Variation in the vitamin D receptor gene is associated with multiple sclerosis in an Australian population

Multiple Sclerosis (MS) is a chronic inflammatory demyelinating disease of the central nervous system (CNS) resulting in accumulating neurological disability. The disorder is more prevalent at higher latitudes. To investigate VDR gene variation using three intragenic restriction fragment length polymorphisms (Apa I, Taq I and Fok I) in an Australian MS case-control population. One hundred and four Australian MS patients were studied with patients classified clinically as Relapsing Remitting MS (RR-MS), Secondary Progressive MS (SP-MS) or Primary Progressive MS (PP-MS). Also, 104 age-, sex-, and ethnicity-matched controls were investigated as a comparative group. Our results show a significant difference of genotype distribution frequency between the case and control groups for the functional exon 9 VDR marker Taq I (p(Gen) = 0.016) and interestingly, a stronger difference for the allelic frequency (p(All) = 0.0072). The Apa I alleles were also found to be associated with MS (p(All) = 0.04) but genotype frequencies were not significantly different from controls (p(Gen) = 0.1). The Taq and Apa variants are in very strong and significant linkage disequilibrium (D' = 0.96, P < 0.0001). The genotypic associations are strongest for the progressive forms of MS (SP-MS and PP-MS). Our results support a role for the VDR gene increasing the risk of developing multiple sclerosis, particularly the progressive clinical subtypes of MS.

Investigation of a neuronal nitric oxide synthase gene (NOS1) polymorphism in a multiple sclerosis population

Multiple Sclerosis (MS) is a chronic neurological disease characterized by demyelination associated with infiltrating white blood cells in the central nervous system (CNS). Nitric oxide synthases (NOS) are a family of enzymes that control the production of nitric oxide. It is possible that neuronal NOS could be involved in MS pathophysiology and hence the nNOS gene is a potential candidate for involvement in disease susceptibility. The aim of this study was to determine whether allelic variation at the nNOS gene locus is associated with MS in an Australian cohort. DNA samples obtained from a Caucasian Australian population affected with MS and an unaffected control population, matched for gender, age and ethnicity, were genotyped for a microsatellite polymorphism in the promoter region of the nNOS gene. Allele frequencies were compared using chi-squared based statistical analyses with significance tested by Monte Carlo simulation. Allelic analysis of MS cases and controls produced a chi-squared value of 5.63 with simulated P = 0.96 (OR(max) = 1.41, 95% CI: 0.926-2.15). Similarly, a Mann-Whitney U analysis gave a non-significant P-value of 0.377 for allele distribution. No differences in allele frequencies were observed for gender or clinical course subtype (P > 0.05). Statistical analysis indicated that there is no association of this nNOS variant and MS and hence the gene does not appear to play a genetically significant role in disease susceptibility.