Knee flexor strength and biceps femoris activation deficits following hamstring strain injury

Background: Hamstring strain injuries (HSI) are prevalent in sport and re-injury rates have been high for many years. Maladaptation following HSI are implicated in injury recurrence however nervous system function following HSI has received little attention. Aim: To determine if recreational athletes with a history of unilateral HSI, who have returned to training and competition, will exhibit lower levels of voluntary activation (VA) and median power frequency (MPF) in the previously injured limb compared to the uninjured limb at long muscle lengths. Methods: Twenty-eight recreational athletes were recruited. Of these, 13 athletes had a history of unilateral HSI and 15 had no history of HSI. Following familiarisation, all athletes undertook isokinetic dynamometry testing and surface electromyography assessment of the biceps femoris long head and medial hamstrings during concentric and eccentric contractions at ± 180 and ± 60deg/s. Results: The previously injured limb was weaker at all contraction speeds compared to the uninjured limb (+180deg/s mean difference(MD) = 9.3Nm, p = 0.0036; +60deg/s MD = 14.0Nm, p = 0.0013; -60deg/s MD = 18.3Nm, p = 0.0007; -180deg/s MD = 20.5Nm, p = 0.0007) whilst VA was only lower in the biceps femoris long head during eccentric contractions (-60deg/s MD = 0.13, p = 0.0025; -180deg/s MD = 0.13, p = 0.0003). There were no between limb differences in medial hamstring VA or MPF from either biceps femoris long head or medial hamstrings in the injured group. The uninjured group showed no between limb differences with any of the tested variables. Conclusion: Previously injured hamstrings were weaker than the contralateral uninjured hamstring at all tested speeds and contraction modes. During eccentric contractions biceps femoris long head VA was lower in the previously injured limb suggesting neural control of biceps femoris long head may be altered following HSI. Current rehabilitation practices have been unsuccessful in restoring strength and VA following HSI. Restoration of these markers should be considered when determining the success of rehabilitation from HSI. Further investigations are required to elucidate the full impact of lower levels of biceps femoris long head VA following HSI on rehabilitation outcomes and re-injury risk.

Identification of GABA receptors in chick cornea

Purpose: The cornea has an important role in vision, is highly innervated and many neurotransmitter receptors are present, e.g., muscarine, melatonin, and dopamine receptors. γ-aminobutyric acid (GABA) is the most important inhibitory neurotransmitter in the retina and central nervous system, but it is unknown whether GABA receptors are present in cornea. The aim of this study was to determine if GABA receptors are located in chick cornea. Methods: Corneal tissues were collected from 25, 12-day-old chicks. Real time PCR, western blot, and immunohistochemistry were used to determine whether alpha1 GABAA, GABAB, and rho1 GABAC receptors were expressed and located in chick cornea. Results: Corneal tissue was positive for alpha1 GABAA and rho1 GABAC receptor mRNA (PCR) and protein (western blot) expression but was negative for GABAB receptor mRNA and protein. Alpha1 GABAA and rho1 GABAC receptor protein labeling was observed in the corneal epithelium using immunohistochemistry. Conclusions: These investigations clearly show that chick cornea possesses alpha1 GABAA, and rho1 GABAC receptors, but not GABAB receptors. The purpose of the alpha1 GABAA and rho1 GABAC receptors in cornea is a fascinating unexplored question.

The hyperthermic patient

Concepts used in this chapter include: Thermoregulation:- Thermoregulation refers to the body’s sophisticated, multi-system regulation of core body temperature. This hierarchical system extends from highly thermo-sensitive neurons in the preoptic region of the brain proximate to the rostral hypothalamus, down to the brain stem and spinal cord. Coupled with receptors in the skin and spine, both central and peripheral information on body temperature is integrated to inform and activate the homeostatic mechanisms which maintain our core temperature at 37oC1. Hyperthermia:- An imbalance between the metabolic and external heat accumulated in the body and the loss of heat from the body2. Exertional heat stroke:- A disorder of excessive heat production coupled with insufficient heat dissipation which occurs in un-acclimated individuals who are engaging in over-exertion in hot and humid conditions. This phenomenon includes central nervous system dysfunction and critical dysfunction to all organ systems including renal, cardiovascular, musculoskeletal and hepatic functions. Non-exertional heat stroke:- In contrast to exertional heatstroke as a consequence of high heat production during strenuous exercise, non-exertional heatstroke results from prolonged exposure to high ambient temperature. The elderly, those with chronic health conditions and children are particularly susceptible.3 Rhabdomylosis:- An acute, sometimes fatal disease characterised by destruction of skeletal muscle. In exertional heat stroke, rhabdomylosis occurs in the context of strenuous exercise when mechanical and/or metabolic stress damages the skeletal muscle, causing elevated serum creatine kinease. Associated with this is the potential development of hyperkalemia, myoglobinuria and renal failure. Malignant hyperthermia:- Malignant hyperthermia is “an inherited subclinical myopathy characterised by a hypermetabolic reaction during anaesthesia. The reaction is related to skeletal muscle calcium dysregulation triggered by volatile inhaled anaesthetics and/or succinylcholine.”4 Presentation includes skeletal muscle rigidity, mixed metabolic and respiratory acidosis, tachycardia, hyperpyrexia, rhabdomylosis, hyperkalaemia, elevated serum creatine kinease, multi-organ failure, disseminated intravascular coagulation and death.5

Cerebral cortex activation mapping upon electrical muscle stimulation by 32-channel time-domain functional near-infrared spectroscopy

The application of different EMS current thresholds on muscle activates not only the muscle but also peripheral sensory axons that send proprioceptive and pain signals to the cerebral cortex. A 32-channel time-domain fNIRS instrument was employed to map regional cortical activities under varied EMS current intensities applied on the right wrist extensor muscle. Eight healthy volunteers underwent four EMS at different current thresholds based on their individual maximal tolerated intensity (MTI), i.e., 10 % < 50 % < 100 % < over 100 % MTI. Time courses of the absolute oxygenated and deoxygenated hemoglobin concentrations primarily over the bilateral sensorimotor cortical (SMC) regions were extrapolated, and cortical activation maps were determined by general linear model using the NIRS-SPM software. The stimulation-induced wrist extension paradigm significantly increased activation of the contralateral SMC region according to the EMS intensities, while the ipsilateral SMC region showed no significant changes. This could be due in part to a nociceptive response to the higher EMS current intensities and result also from increased sensorimotor integration in these cortical regions.

Resequencing and fine-mapping of the chromosome 12q13-14 locus associated with multiple sclerosis refines the number of implicated genes

Multiple sclerosis (MS) is a common chronic inflammatory disease of the central nervous system. Susceptibility to the disease is affected by both environmental and genetic factors. Genetic factors include haplotypes in the histocompatibility complex (MHC) and over 50 non-MHC loci reported by genome-wide association studies. Amongst these, we previously reported polymorphisms in chromosome 12q13-14 with a protective effect in individuals of European descent. This locus spans 288 kb and contains 17 genes, including several candidate genes which have potentially significant pathogenic and therapeutic implications. In this study, we aimed to fine-map this locus. We have implemented a two-phase study: a variant discovery phase where we have used next-generation sequencing and two target-enrichment strategies [long-range polymerase chain reaction (PCR) and Nimblegen's solution phase hybridization capture] in pools of 25 samples; and a genotyping phase where we genotyped 712 variants in 3577 healthy controls and 3269 MS patients. This study confirmed the association (rs2069502, P = 9.9 × 10−11, OR = 0.787) and narrowed down the locus of association to an 86.5 kb region. Although the study was unable to pinpoint the key-associated variant, we have identified a 42 (genotyped and imputed) single-nucleotide polymorphism haplotype block likely to harbour the causal variant. No evidence of association at previously reported low-frequency variants in CYP27B1 was observed. As part of the study we compared variant discovery performance using two target-enrichment strategies. We concluded that our pools enriched with Nimblegen's solution phase hybridization capture had better sensitivity to detect true variants than the pools enriched with long-range PCR, whilst specificity was better in the long-range PCR-enriched pools compared with solution phase hybridization capture enriched pools; this result has important implications for the design of future fine-mapping studies.

Studies on the pathophysiology and genetic basis of migraine

Migraine is a neurological disorder that affects the central nervous system causing painful attacks of headache. A genetic vulnerability and exposure to environmental triggers can influence the migraine phenotype. Migraine interferes in many facets of people’s daily life including employment commitments and their ability to look after their families resulting in a reduced quality of life. Identification of the biological processes that underlie this relatively common affliction has been difficult because migraine does not have any clearly identifiable pathology or structural lesion detectable by current medical technology. Theories to explain the symptoms of migraine have focused on the physiological mechanisms involved in the various phases of headache and include the vascular and neurogenic theories. In relation to migraine pathophysiology the trigeminovascular system and cortical spreading depression have also been implicated with supporting evidence from imaging studies and animal models. The objective of current research is to better understand the pathways and mechanisms involved in causing pain and headache to be able to target interventions. The genetic component of migraine has been teased apart using linkage studies and both candidate gene and genome-wide association studies, in family and case-control cohorts. Genomic regions that increase individual risk to migraine have been identified in neurological, vascular and hormonal pathways. This review discusses knowledge of the pathophysiology and genetic basis of migraine with the latest scientific evidence from genetic studies.

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.

Semaphorin–plexin signalling genes associated with human breast tumourigenesis

Introduction Gene expression profiling has enabled us to demonstrate the heterogeneity of breast cancers. The potential of a tumour to grow and metastasise is partly dependant on its ability to initiate angiogenesis or growth and remodelling of new blood vessels, usually from a pre-existing vascular network, to ensure delivery of oxygen, nutrients, and growth factors to rapidly dividing transformed cells along with access to the systemic circulation. Cell–cell signalling of semaphorin ligands through interaction with their plexin receptors is important for the homeostasis and morphogenesis of many tissues and has been widely studied for a role in neural connectivity, cancer, cell migration and immune responses. This study investigated the role of four semaphorin/plexin signalling genes in human breast cancers in vivo and in vitro. Materials and methods mRNA was extracted from formalin fixed paraffin embedded archival breast invasive ductal carcinoma tissue samples of progressive grades (grades I–III) and compared to tissue from benign tumours. Gene expression profiles were determined by microarray using the Affymetrix GeneChip® Human Genome U133 Plus 2.0 Arrays and validated by Q-PCR using a Corbett RotorGene 6000. Following validation, the gene expression profile of the identified targets was correlated with those of the human breast cancer cell lines MCF-7 and MDA-MD-231. Results The array data revealed that 888 genes were found to be significantly (p ≤ 0.05) differentially expressed between grades I and II tumours and 563 genes between grade III and benign tumours. From these genes, we identified four genes involved in semaphorin–plexin signalling including SEMA4D which has previously been identified as being involved in increased angiogenesis in breast cancers, and three other genes, SEMA4F, PLXNA2 and PLXNA3, which in the literature were associated with tumourigenesis, but not directly in breast tumourigenesis. The microarray analysis revealed that SEMA4D was significantly (P = 0.0347) down-regulated in the grade III tumours compared to benign tumours; SEMA4F, was significantly (P = 0.0159) down-regulated between grades I and II tumours; PLXNA2 was significantly (P = 0.036) down-regulated between grade III and benign tumours and PLXNA3 significantly (P = 0.042) up-regulated between grades I and II tumours. Gene expression of SEMA4D was validated using Q-PCR, demonstrating the same expression profile in both data sets. When the sample set was increased to incorporate more cases, SEMA4D continued to follow the same expression profile, including statistical significance for the differences observed and small standard deviations. In vitro the same pattern was present where expression for SEMA4D was significantly higher in MDA-MB-231 cells when compared to MCF-7 cells. The expression of SEMA4F, PLXNA2 and PLXNA3 could not be validated using Q-PCR, however in vitro analysis of these three genes revealed that both SEMA4F and PLXNA3 followed the microarray trend in expression, although they did not reach significance. In contrast, PLXNA2 demonstrated statistical significance and was in concordance with the literature. Discussion We, and others, have proposed SEMA4D to be a gene with a potentially protective effect in benign tumours that contributes to tumour growth and metastatic suppression. Previous data supports a role for SEMA4F as a tumour suppressor in the peripheral nervous system but our data seems to indicate that the gene is involved in tumour progression in breast cancer. Our in vitro analysis of PLXNA2 revealed that the gene has higher expression in more aggressive breast cancer cell types. Finally, our in vitro analysis on PLXNA3 also suggest that this gene may have some form of growth suppressive role in breast cancer, in addition to a similar role for the gene previously reported in ovarian cancer. From the data obtained in this study, SEMA4D may have a role in more aggressive and potentially metastatic breast tumours. Conclusions Semaphorins and their receptors, the plexins, have been implicated in numerous aspects of neural development, however their expression in many other epithelial tissues suggests that the semaphorin–plexin signalling system also contributes to blood vessel growth and development. These findings warrant further investigation of the role of semaphorins and plexins and their role in normal and tumour-induced angiogenesis in vivo and in vitro. This may represent a new front of attack in anti-angiogenic therapies of breast and other cancers.

Variants in the human potassium channel gene (KCNN3) are associated with migraine in a high risk genetic isolate

The calcium-activated potassium ion channel gene (KCNN3) is located in the vicinity of the familial hemiplegic migraine type 2 locus on chromosome 1q21.3. This gene is expressed in the central nervous system and plays a role in neural excitability. Previous association studies have provided some, although not conclusive, evidence for involvement of this gene in migraine susceptibility. To elucidate KCNN3 involvement in migraine, we performed gene-wide SNP genotyping in a high-risk genetic isolate from Norfolk Island, a population descended from a small number of eighteenth century Isle of Man ‘Bounty Mutineer’ and Tahitian founders. Phenotype information was available for 377 individuals who are related through the single, well-defined Norfolk pedigree (96 were affected: 64 MA, 32 MO). A total of 85 SNPs spanning the KCNN3 gene were genotyped in a sub-sample of 285 related individuals (76 affected), all core members of the extensive Norfolk Island ‘Bounty Mutineer’ genealogy. All genotyping was performed using the Illumina BeadArray platform. The analysis was performed using the statistical program SOLAR v4.0.6 assuming an additive model of allelic effect adjusted for the effects of age and sex. Haplotype analysis was undertaken using the program HAPLOVIEW v4.0. A total of four intronic SNPs in the KCNN3 gene displayed significant association (P < 0.05) with migraine. Two SNPs, rs73532286 and rs6426929, separated by approximately 0.1 kb, displayed complete LD (r 2 = 1.00, D′ = 1.00, D′ 95% CI = 0.96–1.00). In all cases, the minor allele led to a decrease in migraine risk (beta coefficient = 0.286–0.315), suggesting that common gene variants confer an increased risk of migraine in the Norfolk pedigree. This effect may be explained by founder effect in this genetic isolate. This study provides evidence for association of variants in the KCNN3 ion channel gene with migraine susceptibility in the Norfolk genetic isolate with the rarer allelic variants conferring a possible protective role. This the first comprehensive analysis of this potential candidate gene in migraine and also the first study that has utilised the unique Norfolk Island large pedigree isolate to implicate a specific migraine gene. Studies of additional variants in KCNN3 in the Norfolk pedigree are now required (e.g. polyglutamine variants) and further analyses in other population data sets are required to clarify the association of the KCNN3 gene and migraine risk in the general outbred population.

Common variants in the regulative regions of GRIA1 and GRIA3 receptor genes are associated with migraine susceptibility

Background Glutamate is the principal excitatory neurotransmitter in the central nervous system which acts by the activation of either ionotropic (AMPA, NMDA and kainate receptors) or G-protein coupled metabotropic receptors. Glutamate is widely accepted to play a major role in the path physiology of migraine as implicated by data from animal and human studies. Genes involved in synthesis, metabolism and regulation of both glutamate and its receptors could be, therefore, considered as potential candidates for causing/predisposing to migraine when mutated. Methods The association of polymorphic variants of GRIA1-GRIA4 genes which encode for the four subunits (GluR1-GluR4) of the alpha-amino-3- hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) receptor for glutamate was tested in migraineurs with and without aura (MA and MO) and healthy controls. Results Two variants in the regulative regions of GRIA1 (rs2195450) and GRIA3 (rs3761555) genes resulted strongly associated with MA (P = 0.00002 and P = 0.0001, respectively), but not associated with MO, suggesting their role in cortical spreading depression. Whereas the rs548294 variant in GRIA1 gene showed association primarily with MO phenotype, supporting the hypothesis that MA and MO phenotypes could be genetically related. These variants modify binding sites for transcription factors altering the expression of GRIA1 and GRIA3 genes in different conditions. Conclusions This study represents the first genetic evidence of a link between glutamate receptors and migraine.

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.

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.

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.