High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis

Background

The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods.

Methods

We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls.

Results

We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk.

Results

The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

Inflammatory and neurodegenerative serum protein biomarkers increase sensitivity to detect disease activity in multiple sclerosis

Background/Objectives

Serum proteomic analysis of deeply-phenotyped samples, biological pathway modeling and network analysis were performed to elucidate the inflammatory and neurodegenerative processes of multiple sclerosis (MS) and identify sensitive biomarkers of MS disease activity (DA).

Methods

Over 1100 serum proteins were evaluated in >600 samples from three MS cohorts to identify biomarkers of clinical and radiographic (gadolinium-enhancing lesions) new MS DA. Protein levels were analyzed and associated with presence of gadolinium-enhancing lesions, clinical relapse status (CRS), and annualized relapse rate (ARR) to create a custom assay panel.

Results

Twenty proteins were associated with increased clinical and radiographic MS DA. Serum neurofilament light chain (NfL) showed the strongest univariate correlation with radiographic and clinical DA measures. Multivariate modeling significantly outperformed univariate NfL to predict gadolinium lesion activity, CRS and ARR.

Discussion

These findings provide insight regarding correlations between inflammatory and neurodegenerative biomarkers and clinical and radiographic MS DA.

Funding

Octave Bioscience, Inc (Menlo Park, CA).

Gut microbiome of multiple sclerosis patients and paired household healthy controls reveal associations with disease risk and course

Changes in gut microbiota have been associated with several diseases. Here, the International Multiple Sclerosis Microbiome Study (iMSMS) studied the gut microbiome of 576 MS patients (36% untreated) and genetically unrelated household healthy controls (1,152 total subjects). We observed a significantly increased proportion of Akkermansia muciniphila, Ruthenibacterium lactatiformans, Hungatella hathewayi, and Eisenbergiella tayi and decreased Faecalibacterium prausnitzii and Blautia species. The phytate degradation pathway was over-represented in untreated MS, while pyruvate-producing carbohydrate metabolism pathways were significantly reduced. Microbiome composition, function, and derived metabolites also differed in response to disease-modifying treatments. The therapeutic activity of interferon-β may in part be associated with upregulation of short-chain fatty acid transporters. Distinct microbial networks were observed in untreated MS and healthy controls. These results strongly support specific gut microbiome associations with MS risk, course and progression, and functional changes in response to treatment.

Specific hypomethylation programs underpin B cell activation in early multiple sclerosis

Epigenetic changes have been consistently detected in different cell types in multiple sclerosis (MS). However, their contribution to MS pathogenesis remains poorly understood partly because of sample heterogeneity and limited coverage of array-based methods. To fill this gap, we conducted a comprehensive analysis of genome-wide DNA methylation patterns in four peripheral immune cell populations isolated from 29 MS patients at clinical disease onset and 24 healthy controls. We show that B cells from new-onset untreated MS cases display more significant methylation changes than other disease-implicated immune cell types, consisting of a global DNA hypomethylation signature. Importantly, 4,933 MS-associated differentially methylated regions in B cells were identified, and this epigenetic signature underlies specific genetic programs involved in B cell differentiation and activation. Integration of the methylome to changes in gene expression and susceptibility-associated regions further indicates that hypomethylated regions are significantly associated with the up-regulation of cell activation transcriptional programs. Altogether, these findings implicate aberrant B cell function in MS etiology.

Next-Generation Sequencing Identifies Extended HLA Class I and II Haplotypes Associated With Early-Onset and Late-Onset Myasthenia Gravis in Italian, Norwegian, and Swedish Populations

Genetic susceptibility to myasthenia gravis (MG) associates with specific HLA alleles and haplotypes at the class I and II regions in various populations. Previous studies have only examined alleles at a limited number of HLA loci that defined only broad serotypes or alleles defined at the protein sequence level. Consequently, genetic variants in noncoding and untranslated HLA gene segments have not been fully explored but could also be important determinants for MG. To gain further insight into the role of HLA in MG, we applied next-generation sequencing to analyze sequence variation at eleven HLA genes in early-onset (EO) and late-onset (LO) non-thymomatous MG patients positive for the acetylcholine receptor (AChR) antibodies and ethnically matched controls from Italy, Norway, and Sweden. For all three populations, alleles and haplotype blocks present on the ancestral haplotype AH8.1 were associated with risk in AChR-EOMG patients. HLA-B*08:01:01:01 was the dominant risk allele in Italians (OR = 3.28, P = 1.83E-05), Norwegians (OR = 3.52, P = 4.41E-16), and in Swedes HLA-B*08:01 was the primary risk allele (OR = 4.24, P <2.2E-16). Protective alleles and haplotype blocks were identified on the HLA-DRB7, and HLA-DRB13.1 class II haplotypes in Italians and Norwegians, whereas in Swedes HLA-DRB7 exhibited the main protective effect. For AChR-LOMG patients, the HLA-DRB15.1 haplotype and associated alleles were significantly associated with susceptibility in all groups. The HLA-DR13-HLA-DR-HLA-DQ haplotype was associated with protection in all AChR-LOMG groups. This study has confirmed and extended previous findings that the immunogenetic predisposition profiles for EOMG and LOMG are distinct. In addition, the results are consistent with a role for non-coding HLA genetic variants in the pathogenesis of MG.

High Resolution Haplotype Analyses of Classical HLA Genes in Families With Multiple Sclerosis Highlights the Role of HLA-DP Alleles in Disease Susceptibility

Multiple sclerosis (MS) susceptibility shows strong genetic associations with HLA alleles and haplotypes. We genotyped 11 HLA genes in 477 non-Hispanic European MS patients and their 954 unaffected parents using a validated next-generation sequencing (NGS) methodology. HLA haplotypes were assigned unequivocally by tracing HLA allele transmissions. We explored HLA haplotype/allele associations with MS using the genotypic transmission disequilibrium test (gTDT) and multiallelic TDT (mTDT). We also conducted a case-control (CC) study with all patients and 2029 healthy unrelated ethnically matched controls. We performed separate analyses of 54 extended multi-case families by reviewing transmission of haplotype blocks. The haplotype fragment including DRB5*01:01:01~DRB1*15:01:01:01 was significantly associated with predisposition (gTDT: p < 2.20e-16; mTDT: p =1.61e-07; CC: p < 2.22e-16) as reported previously. A second risk allele, DPB1*104:01 (gTDT: p = 3.69e-03; mTDT: p = 2.99e-03; CC: p = 1.00e-02), independent from the haplotype bearing DRB1*15:01 was newly identified. The allele DRB1*01:01:01 showed significant protection (gTDT: p = 8.68e-06; mTDT: p = 4.50e-03; CC: p = 1.96e-06). Two DQB1 alleles, DQB1*03:01 (gTDT: p = 2.86e-03; mTDT: p = 5.56e-02; CC: p = 4.08e-05) and DQB1*03:03 (gTDT: p = 1.17e-02; mTDT: p = 1.16e-02; CC: p = 1.21e-02), defined at two-field level also showed protective effects. The HLA class I block, A*02:01:01:01~C*03:04:01:01~B*40:01:02 (gTDT: p = 5.86e-03; mTDT: p = 3.65e-02; CC: p = 9.69e-03) and the alleles B*27:05 (gTDT: p = 6.28e-04; mTDT: p = 2.15e-03; CC: p = 1.47e-02) and B*38:01 (gTDT: p = 3.20e-03; mTDT: p = 6.14e-03; CC: p = 1.70e-02) showed moderately protective effects independently from each other and from the class II associated factors. By comparing statistical significance of 11 HLA loci and 19 haplotype segments with both untruncated and two-field allele names, we precisely mapped MS candidate alleles/haplotypes while eliminating false signals resulting from ‘hitchhiking’ alleles. We assessed genetic burden for the HLA allele/haplotype identified in this study. This family-based study including the highest-resolution of HLA alleles proved to be powerful and efficient for precise identification of HLA genotypes associated with both, susceptibility and protection to development of MS.

Cell type-specific transcriptomics identifies neddylation as a novel therapeutic target in multiple sclerosis

Multiple sclerosis is an autoimmune disease of the CNS in which both genetic and environmental factors are involved. Genome-wide association studies revealed more than 200 risk loci, most of which harbour genes primarily expressed in immune cells. However, whether genetic differences are translated into cell-specific gene expression profiles and to what extent these are altered in patients with multiple sclerosis are still open questions in the field. To assess cell type-specific gene expression in a large cohort of patients with multiple sclerosis, we sequenced the whole transcriptome of fluorescence-activated cell sorted T cells (CD4+ and CD8+) and CD14+ monocytes from treatment-naive patients with multiple sclerosis (n = 106) and healthy subjects (n = 22). We identified 479 differentially expressed genes in CD4+ T cells, 435 in monocytes, and 54 in CD8+ T cells. Importantly, in CD4+ T cells, we discovered upregulated transcripts from the NAE1 gene, a critical subunit of the NEDD8 activating enzyme, which activates the neddylation pathway, a post-translational modification analogous to ubiquitination. Finally, we demonstrated that inhibition of NEDD8 activating enzyme using the specific inhibitor pevonedistat (MLN4924) significantly ameliorated disease severity in murine experimental autoimmune encephalomyelitis. Our findings provide novel insights into multiple sclerosis-associated gene regulation unravelling neddylation as a crucial pathway in multiple sclerosis pathogenesis with implications for the development of tailored disease-modifying agents.

A splice acceptor variant in HLA-DRA affects the conformation and cellular localization of the class II DR alpha-chain

Class II human leucocyte antigen (HLA) proteins are involved in the immune response by presenting pathogen-derived peptides to CD4⁺ T lymphocytes. At the molecular level, they are constituted by α/β-heterodimers on the surface of professional antigen-presenting cells. Here, we report that the acceptor variant (rs8084) in the HLA-DRA gene mediates the transcription of an alternative version of the α-chain lacking 25 amino acids in its extracellular domain. Molecular dynamics simulations suggest this isoform undergoes structural refolding which in turn affects its stability and cellular trafficking. The short HLA-DRA isoform cannot reach the cell surface, although it is still able to bind the corresponding β-chain. Conversely, it remains entrapped within the endoplasmic reticulum where it is targeted for degradation. Furthermore, we demonstrate that the short isoform can be transported to the cell membrane via interactions with the peptide-binding site of canonical HLA heterodimers. Altogether, our findings indicate that short HLA-DRA functions as a novel intact antigen for class II HLA molecules.

Household paired design reduces variance and increases power in multi-city gut microbiome study in multiple sclerosis

BACKGROUND

Evidence for a role of human gut microbiota in multiple sclerosis (MS) risk is mounting, yet large variability is seen across studies. This is, in part, due to the lack of standardization of study protocols, sample collection methods, and sequencing approaches.

OBJECTIVE

This study aims to address the effect of a household experimental design, sample collection, and sequencing approaches in a gut microbiome study in MS subjects from a multi-city study population.

METHODS

We analyzed 128 MS patient and cohabiting healthy control pairs from the International MS Microbiome Study (iMSMS). A total of 1005 snap-frozen or desiccated Q-tip stool samples were collected and evaluated using 16S and shallow whole-metagenome shotgun sequencing.

RESULTS

The intra-individual variance observed by different collection strategies was dramatically lower than inter-individual variance. Shallow shotgun highly correlated with 16S sequencing. Participant house and recruitment site accounted for the two largest sources of microbial variance, while higher microbial similarity was seen in household-matched participants as hypothesized. A significant proportion of the variance in dietary intake was also dominated by geographic distance.

CONCLUSION

A household pair study largely overcomes common inherent limitations and increases statistical power in population-based microbiome studies.

Association Between Serum Neurofilament Light Chain Levels and Long-term Disease Course Among Patients With Multiple Sclerosis Followed up for 12 Years

Importance

Blood sample-based biomarkers that are associated with clinically meaningful outcomes for patients with multiple sclerosis (MS) have not been developed.

Objective

To evaluate the potential of serum neurofilament light chain (sNFL) measurements as a biomarker of disease activity and progression in a longitudinal MS data set.

Design, Setting, and Participants

Single-center, ongoing, prospective observational cohort study of 607 patients with MS from the longitudinal EPIC (Expression, Proteomics, Imaging, Clinical) study at the University of California, San Francisco from July 1, 2004, through August 31, 2017. Clinical evaluations and sample collection were performed annually for 5 years, then at different time points for up to 12 years, with a median follow-up duration of 10 (interquartile range, 7-11) years. Serum NFL levels were measured using a sensitive single molecule array platform and compared with clinical and magnetic resonance imaging variables with the use of univariable and multivariable analyses.

Main Outcomes and Measures

The main outcomes were disability progression defined as clinically significant worsening on the Expanded Disability Status Scale (EDSS) score and brain fraction atrophy.

Results

Mean (SD) age of the 607 study participants at study entry was 42.5 (9.8) years; 423 (69.7%) were women; and all participants were of non-Hispanic European descent. Of 3911 samples sequentially collected, 3904 passed quality control for quantification of sNFL. Baseline sNFL levels showed significant associations with EDSS score (β, 1.080; 95% CI, 1.047-1.114; P < .001), MS subtype (β, 1.478; 95% CI, 1.279-1.707; P < .001), and treatment status (β, 1.120; 95% CI, 1.007-1.245; P = .04). A significant interaction between EDSS worsening and change in levels of sNFL over time was found (β, 1.015; 95% CI, 1.007-1.023; P < .001). Baseline sNFL levels alone were associated with approximately 11.6% of the variance in brain fraction atrophy at year 10. In a multivariable analysis that considered sex, age, and disease duration, baseline sNFL levels were associated with 18.0% of the variance in brain fraction atrophy at year 10. After 5 years' follow-up, active treatment was associated with lower levels of sNFL, with high-potency treatments associated with the greater decreases in sNFL levels compared with platform therapies (high-potency vs untreated: β, 0.946; 95% CI, 0.915-0.976; P < .001; high-potency vs platform: β, 0.972; 95% CI, 0.948-0.998; P = .04).

Conclusions and Relevance

This study found that statistically significant associations of sNFL with relevant clinical and neuroimaging outcomes in MS were confirmed and extended, supporting the potential of sNFL as an objective surrogate of ongoing MS disease activity. In this data set of patients with MS who received early treatment, the prognostic power of sNFL for relapse activity and long-term disability progression was limited. Further prospective studies are necessary to assess the assay's utility for decision-making in individual patients.

A systems biology approach uncovers cell-specific gene regulatory effects of genetic associations in multiple sclerosis

Genome-wide association studies (GWAS) have identified more than 50,000 unique associations with common human traits. While this represents a substantial step forward, establishing the biology underlying these associations has proven extremely difficult. Even determining which cell types and which particular gene(s) are relevant continues to be a challenge. Here, we conduct a cell-specific pathway analysis of the latest GWAS in multiple sclerosis (MS), which had analyzed a total of 47,351 cases and 68,284 healthy controls and found more than 200 non-MHC genome-wide associations. Our analysis identifies pan immune cell as well as cell-specific susceptibility genes in T cells, B cells and monocytes. Finally, genotype-level data from 2,370 patients and 412 controls is used to compute intra-individual and cell-specific susceptibility pathways that offer a biological interpretation of the individual genetic risk to MS. This approach could be adopted in any other complex trait for which genome-wide data is available.

Sex-specific Tau methylation patterns and synaptic transcriptional alterations are associated with neural vulnerability during chronic neuroinflammation

The molecular events underlying the transition from initial inflammatory flares to the progressive phase of multiple sclerosis (MS) remain poorly understood. Here, we report that the microtubule-associated protein (MAP) Tau exerts a gender-specific protective function on disease progression in the MS model experimental autoimmune encephalomyelitis (EAE). A detailed investigation of the autoimmune response in Tau-deficient mice excluded a strong immunoregulatory role for Tau, suggesting that its beneficial effects are presumably exerted within the central nervous system (CNS). Spinal cord transcriptomic data show increased synaptic dysfunctions and alterations in the NF-kB activation pathway upon EAE in Tau-deficient mice as compared to wildtype animals. We also performed the first comprehensive characterization of Tau post-translational modifications (PTMs) in the nervous system upon EAE. We report that the methylation levels of the conserved lysine residue K306 are significantly decreased in the chronic phase of the disease. By combining biochemical assays and molecular dynamics (MD) simulations, we demonstrate that methylation at K306 decreases the affinity of Tau for the microtubule network. Thus, the down-regulation of this PTM might represent a homeostatic response to enhance axonal stability against an autoimmune CNS insult. The results, altogether, position Tau as key mediator between the inflammatory processes and neurodegeneration that seems to unify many CNS diseases.

Low-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis Risk

Multiple sclerosis is a complex neurological disease, with ∼20% of risk heritability attributable to common genetic variants, including >230 identified by genome-wide association studies. Multiple strands of evidence suggest that much of the remaining heritability is also due to additive effects of common variants rather than epistasis between these variants or mutations exclusive to individual families. Here, we show in 68,379 cases and controls that up to 5% of this heritability is explained by low-frequency variation in gene coding sequence. We identify four novel genes driving MS risk independently of common-variant signals, highlighting key pathogenic roles for regulatory T cell homeostasis and regulation, IFNγ biology, and NFκB signaling. As low-frequency variants do not show substantial linkage disequilibrium with other variants, and as coding variants are more interpretable and experimentally tractable than non-coding variation, our discoveries constitute a rich resource for dissecting the pathobiology of MS.

Deconstruction of HLA-DRB1*04:01:01 and HLA-DRB1*15:01:01 class II haplotypes using next-generation sequencing in European-Americans with multiple sclerosis

BACKGROUND

The association between HLA-DRB1*15:01 with multiple sclerosis (MS) susceptibility is well established, but the contribution of the tightly associated HLA-DRB5*01:01 allele has not yet been completely ascertained. Similarly, the effects of HLA-DRB1*04:01 alleles and haplotypes, defined at the full-gene resolution level with MS risk remains to be elucidated.

OBJECTIVES

To characterize the molecular architecture of class II HLA-DR15 and HLA-DR4 haplotypes associated with MS.

METHODS

Next-generation sequencing was used to determine HLA-DQB1, HLA-DQA1, and HLA-DRB1/4/5 alleles in 1403 unrelated European-American patients and 1425 healthy unrelated controls. Effect sizes of HLA alleles and haplotypes on MS risk were measured by odds ratio (OR) with 95% confidence intervals.

RESULTS

HLA-DRB1*15:01:01:01SG (OR = 3.20, p < 2.2E-16), HLA-DRB5*01:01:01 (OR = 2.96, p < 2.2E-16), and HLA-DRB5*01:01:01v1_STR1 (OR = 8.18, p = 4.3E-05) alleles all occurred at significantly higher frequencies in MS patients compared to controls. The most significant predis-posing haplotypes were HLA-DQB1*06:02:01~ HLA-DQA1*01:02:01:01SG~HLA-DRB1*15:01:01:01SG~HLA-DRB5*01:01:01 and HLA-DQB1*06:02:01~HLA-DQA1*01:02:01:01SG~HLA-DRB1*15:01:01:01SG~HLA-DRB5*01:01:01v1_STR1 (OR = 3.19, p < 2.2E-16; OR = 9.30, p = 9.7E-05, respectively). Analyses of the HLA-DRB1*04 cohort in the absence of HLA-DRB1*15:01 haplotypes revealed that the HLA-DQB1*03:01:01:01~HLA-DQA1*03:03:01:01~HLA-DRB1*04:01:01:01SG~HLA-DRB4*01:03:01:01 haplotype was protective (OR = 0.64, p = 0.028), whereas the HLA-DQB1*03:02:01~HLA-DQA1*03:01:01~HLA-DRB1*04:01:01:01SG~HLA-DRB4*01:03:01:01 haplotype was associated with MS susceptibility (OR = 1.66, p = 4.9E-03).

CONCLUSION

HLA-DR15 haplotypes, including genomic variants of HLA-DRB5, and HLA-DR4 haplotypes affect MS risk.

Association of HLA Genetic Risk Burden With Disease Phenotypes in Multiple Sclerosis

IMPORTANCE

Although multiple HLA alleles associated with multiple sclerosis (MS) risk have been identified, genotype-phenotype studies in the HLA region remain scarce and inconclusive.

OBJECTIVES

To investigate whether MS risk-associated HLA alleles also affect disease phenotypes.

DESIGN, SETTING, AND PARTICIPANTS

A cross-sectional, case-control study comprising 652 patients with MS who had comprehensive phenotypic information and 455 individuals of European origin serving as controls was conducted at a single academic research site. Patients evaluated at the Multiple Sclerosis Center at University of California, San Francisco between July 2004 and September 2005 were invited to participate. Spinal cord imaging in the data set was acquired between July 2013 and March 2014; analysis was performed between December 2014 and December 2015.

MAIN OUTCOMES AND MEASURES

Cumulative HLA genetic burden (HLAGB) calculated using the most updated MS-associated HLA alleles vs clinical and magnetic resonance imaging outcomes, including age at onset, disease severity, conversion time from clinically isolated syndrome to clinically definite MS, fractions of cortical and subcortical gray matter and cerebral white matter, brain lesion volume, spinal cord gray and white matter areas, upper cervical cord area, and the ratio of gray matter to the upper cervical cord area. Multivariate modeling was applied separately for each sex data set.

RESULTS

Of the 652 patients with MS, 586 had no missing genetic data and were included in the HLAGB analysis. In these 586 patients (404 women [68.9%]; mean [SD] age at disease onset, 33.6 [9.4] years), HLAGB was higher than in controls (median [IQR], 0.7 [0-1.4] and 0 [-0.3 to 0.5], respectively; P = 1.8 × 10-27). A total of 619 (95.8%) had relapsing-onset MS and 27 (4.2%) had progressive-onset MS. No significant difference was observed between relapsing-onset MS and primary progressive MS. A higher HLAGB was associated with younger age at onset and the atrophy of subcortical gray matter fraction in women with relapsing-onset MS (standard β = -1.20 × 10-1; P = 1.7 × 10-2 and standard β = -1.67 × 10-1; P = 2.3 × 10-4, respectively), which were driven mainly by the HLA-DRB1*15:01 haplotype. In addition, we observed the distinct role of the HLA-A*24:02-B*07:02-DRB1*15:01 haplotype among the other common DRB1*15:01 haplotypes and a nominally protective effect of HLA-B*44:02 to the subcortical gray atrophy (standard β = -1.28 × 10-1; P = 5.1 × 10-3 and standard β = 9.52 × 10-2; P = 3.6 × 10-2, respectively).

CONCLUSIONS AND RELEVANCE

We confirm and extend previous observations linking HLA MS susceptibility alleles with disease progression and specific clinical and magnetic resonance imaging phenotypic traits.

Multiple sclerosis risk loci and disease severity in 7,125 individuals from 10 studies

OBJECTIVE

We investigated the association between 52 risk variants identified through genome-wide association studies and disease severity in multiple sclerosis (MS).

METHODS

Ten unique MS case data sets were analyzed. The Multiple Sclerosis Severity Score (MSSS) was calculated using the Expanded Disability Status Scale at study entry and disease duration. MSSS was considered as a continuous variable and as 2 dichotomous variables (median and extreme ends; MSSS of ≤5 vs >5 and MSSS of <2.5 vs ≥7.5, respectively). Single nucleotide polymorphisms (SNPs) were examined individually and as both combined weighted genetic risk score (wGRS) and unweighted genetic risk score (GRS) for association with disease severity. Random-effects meta-analyses were conducted and adjusted for cohort, sex, age at onset, and HLA-DRB1*15:01.

RESULTS

A total of 7,125 MS cases were analyzed. The wGRS and GRS were not strongly associated with disease severity after accounting for cohort, sex, age at onset, and HLA-DRB1*15:01. After restricting analyses to cases with disease duration ≥10 years, associations were null (p value ≥0.05). No SNP was associated with disease severity after adjusting for multiple testing.

CONCLUSIONS

The largest meta-analysis of established MS genetic risk variants and disease severity, to date, was performed. Results suggest that the investigated MS genetic risk variants are not associated with MSSS, even after controlling for potential confounders. Further research in large cohorts is needed to identify genetic determinants of disease severity using sensitive clinical and MRI measures, which are critical to understanding disease mechanisms and guiding development of effective treatments.

A non-synonymous single-nucleotide polymorphism associated with multiple sclerosis risk affects the EVI5 interactome

Despite recent progress in the characterization of genetic loci associated with multiple sclerosis (MS) risk, the ubiquitous linkage disequilibrium operating across the genome has stalled efforts to distinguish causative variants from proxy single-nucleotide polymorphisms (SNPs). Here, we have identified through fine mapping and meta-analysis EVI5 as the most plausible disease risk gene within the 1p22.1 locus. We further show that an exonic SNP associated with risk induces changes in superficial hydrophobicity patterns of the coiled-coil domain of EVI5, which, in turns, affects the EVI5 interactome. Immunoprecipitation of wild-type and mutated EVI5 followed by mass spectrometry generated a roster of disease-specific interactors functionally linked to lipid metabolism. Among the exclusive binding partners of the risk variant, we describe the novel interaction with sphingosine 1-phosphate lyase (SGPL1)-a key enzyme for the creation of the sphingosine-1 phosphate gradient, which is relevant to the pathogenic process and therapeutic management of MS.

Mitochondrial DNA sequence variation in multiple sclerosis

OBJECTIVE

To assess the influence of common mitochondrial DNA (mtDNA) sequence variation on multiple sclerosis (MS) risk in cases and controls part of an international consortium.

METHODS

We analyzed 115 high-quality mtDNA variants and common haplogroups from a previously published genome-wide association study among 7,391 cases from the International Multiple Sclerosis Genetics Consortium and 14,568 controls from the Wellcome Trust Case Control Consortium 2 project from 7 countries. Significant single nucleotide polymorphism and haplogroup associations were replicated in 3,720 cases and 879 controls from the University of California, San Francisco.

RESULTS

An elevated risk of MS was detected among haplogroup JT carriers from 7 pooled clinic sites (odds ratio [OR] = 1.15, 95% confidence interval [CI] = 1.07-1.24, p = 0.0002) included in the discovery study. The increased risk of MS was observed for both haplogroup T (OR = 1.17, 95% CI = 1.06-1.29, p = 0.002) and haplogroup J carriers (OR = 1.11, 95% CI = 1.01-1.22, p = 0.03). These haplogroup associations with MS were not replicated in the independent sample set. An elevated risk of primary progressive (PP) MS was detected for haplogroup J participants from 3 European discovery populations (OR = 1.49, 95% CI = 1.10-2.01, p = 0.009). This elevated risk was borderline significant in the US replication population (OR = 1.43, 95% CI = 0.99-2.08, p = 0.058) and remained significant in pooled analysis of discovery and replication studies (OR = 1.43, 95% CI = 1.14-1.81, p = 0.002). No common individual mtDNA variants were associated with MS risk.

CONCLUSIONS

Identification and validation of mitochondrial genetic variants associated with MS and PPMS may lead to new targets for treatment and diagnostic tests for identifying potential responders to interventions that target mitochondria.

An ImmunoChip study of multiple sclerosis risk in African Americans

The aims of this study were: (i) to determine to what degree multiple sclerosis-associated loci discovered in European populations also influence susceptibility in African Americans; (ii) to assess the extent to which the unique linkage disequilibrium patterns in African Americans can contribute to localizing the functionally relevant regions or genes; and (iii) to search for novel African American multiple sclerosis-associated loci. Using the ImmunoChip custom array we genotyped 803 African American cases with multiple sclerosis and 1516 African American control subjects at 130 135 autosomal single nucleotide polymorphisms. We conducted association analysis with rigorous adjustments for population stratification and admixture. Of the 110 non-major histocompatibility complex multiple sclerosis-associated variants identified in Europeans, 96 passed stringent quality control in our African American data set and of these, >70% (69) showed over-representation of the same allele amongst cases, including 21 with nominally significant evidence for association (one-tailed test P < 0.05). At a further eight loci we found nominally significant association with an alternate correlated risk-tagging single nucleotide polymorphism from the same region. Outside the regions known to be associated in Europeans, we found seven potentially associated novel candidate multiple sclerosis variants (P < 10(-4)), one of which (rs2702180) also showed nominally significant evidence for association (one-tailed test P = 0.034) in an independent second cohort of 620 African American cases and 1565 control subjects. However, none of these novel associations reached genome-wide significance (combined P = 6.3 × 10(-5)). Our data demonstrate substantial overlap between African American and European multiple sclerosis variants, indicating common genetic contributions to multiple sclerosis risk.

SNP imputation bias reduces effect size determination

Imputation is a commonly used technique that exploits linkage disequilibrium to infer missing genotypes in genetic datasets, using a well-characterized reference population. While there is agreement that the reference population has to match the ethnicity of the query dataset, it is common practice to use the same reference to impute genotypes for a wide variety of phenotypes. We hypothesized that using a reference composed of samples with a different phenotype than the query dataset would introduce imputation bias. To test this hypothesis we used GWAS datasets from Amyotrophic Lateral Sclerosis (ALS), Parkinson Disease (PD), and Crohn's Disease (CD). First, we masked and then performed imputation of 100 disease-associated markers and 100 non-associated markers from each study. Two references for imputation were used in parallel: one consisting of healthy controls and another consisting of patients with the same disease. We assessed the discordance (imprecision) and bias (inaccuracy) of imputation by comparing predicted genotypes to those assayed by SNP-chip. We also assessed the bias on the observed effect size when the predicted genotypes were used in a GWAS study. When healthy controls were used as reference for imputation, a significant bias was observed, particularly in the disease-associated markers. Using cases as reference significantly attenuated this bias. For nearly all markers, the direction of the bias favored the non-risk allele. In GWAS studies of the three diseases (with healthy reference controls from the 1000 genomes as reference), the mean OR for disease-associated markers obtained by imputation was lower than that obtained using original assayed genotypes. We found that the bias is inherent to imputation as using different methods did not alter the results. In conclusion, imputation is a powerful method to predict genotypes and estimate genetic risk for GWAS. However, a careful choice of reference population is needed to minimize biases inherent to this approach.

Prognostic biomarkers of IFNb therapy in multiple sclerosis patients

BACKGROUND

Interferon beta (IFNb) reduces relapse frequency and disability progression in patients with multiple sclerosis (MS).

OBJECTIVES

Early identification of prognostic biomarkers of IFNb-treated patients will allow more effective management of MS.

METHODS

The IMPROVE study evaluated subcutaneous IFNb versus placebo in 180 patients with relapsing-remitting MS. Magnetic resonance imaging scans, clinical assessments, and blood samples were obtained at baseline and every 4 weeks from every participant. Thirty-nine biomarkers (32 transcripts; seven proteins) were studied in 155 patients from IMPROVE. Therapeutic response was defined by absence of new combined unique lesions, relapses, and sustained increase in Expanded Disability Status Scale over 1 year. A machine learning approach was used to examine the association between biomarker expression and treatment response.

RESULTS

While baseline levels of individual genes were relatively poor predictors, combinations of three genes were able to identify subjects with sub-optimal therapeutic responses. The triplet CASP2/IRF4/IRF6, previously identified in an independent dataset, was tested among other combinations. This triplet showed acceptable predictive accuracy (0.68) and specificity (0.88), but had relatively low sensitivity (0.22) resulting in an area under the curve (AUC) of 0.63. Other combinations of biomarkers resulted in AUC of up to 0.80 (e.g. CASP2/IL10/IL12Rb1).

CONCLUSIONS

Baseline expression, or induction ratios, of specific gene combinations correlate with future therapeutic response to IFNb, and have the potential to be prognostically useful.

Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis

Using the ImmunoChip custom genotyping array, we analyzed 14,498 subjects with multiple sclerosis and 24,091 healthy controls for 161,311 autosomal variants and identified 135 potentially associated regions (P < 1.0 × 10(-4)). In a replication phase, we combined these data with previous genome-wide association study (GWAS) data from an independent 14,802 subjects with multiple sclerosis and 26,703 healthy controls. In these 80,094 individuals of European ancestry, we identified 48 new susceptibility variants (P < 5.0 × 10(-8)), 3 of which we found after conditioning on previously identified variants. Thus, there are now 110 established multiple sclerosis risk variants at 103 discrete loci outside of the major histocompatibility complex. With high-resolution Bayesian fine mapping, we identified five regions where one variant accounted for more than 50% of the posterior probability of association. This study enhances the catalog of multiple sclerosis risk variants and illustrates the value of fine mapping in the resolution of GWAS signals.

Oligoclonal bands and age at onset correlate with genetic risk score in multiple sclerosis

BACKGROUND

Many genetic risk variants are now well established in multiple sclerosis (MS), but the impact on clinical phenotypes is unclear.

OBJECTIVE

To investigate the impact of established MS genetic risk variants on MS phenotypes, in well-characterized MS cohorts.

METHODS

Norwegian MS patients (n = 639) and healthy controls (n = 530) were successfully genotyped for 61 established MS-associated single nucleotide polymorphisms (SNPs). Data including and excluding Major Histocompatibility Complex (MHC) markers were summed to a MS Genetic Burden (MSGB) score. Study replication was performed in a cohort of white American MS patients (n = 1997) and controls (n = 708).

RESULTS

The total human leukocyte antigen (HLA) and the non-HLA MSGB scores were significantly higher in MS patients than in controls, in both cohorts (P << 10(-22)). MS patients, with and without cerebrospinal fluid (CSF) oligoclonal bands (OCBs), had a higher MSGB score than the controls; the OCB-positive patients had a slightly higher MSGB than the OCB-negative patients. An early age at symptom onset (AAO) also correlated with a higher MSGB score, in both cohorts.

CONCLUSION

The MSGB score was associated with specific clinical MS characteristics, such as OCBs and AAO. This study underlines the need for well-characterized, large cohorts of MS patients, and the usefulness of summarizing multiple genetic risk factors of modest effect size in genotype-phenotype analyses.

Genetic risk variants in African Americans with multiple sclerosis

OBJECTIVES

To assess the association of established multiple sclerosis (MS) risk variants in 3,254 African Americans (1,162 cases and 2,092 controls).

METHODS

Human leukocyte antigen (HLA)-DRB1, HLA-DQB1, and HLA-A alleles were typed by molecular techniques. Single nucleotide polymorphism (SNP) genotyping was conducted for 76 MS-associated SNPs and 52 ancestry informative marker SNPs selected throughout the genome. Self-declared ancestry was refined by principal component analysis of the ancestry informative marker SNPs. An ancestry-adjusted multivariate model was applied to assess genetic associations.

RESULTS

The following major histocompatibility complex risk alleles were replicated: HLA-DRB1*15:01 (odds ratio [OR] = 2.02 [95% confidence interval: 1.54-2.63], p = 2.50e-07), HLA-DRB1*03:01 (OR = 1.58 [1.29-1.94], p = 1.11e-05), as well as HLA-DRB1*04:05 (OR = 2.35 [1.26-4.37], p = 0.007) and the African-specific risk allele of HLA-DRB1*15:03 (OR = 1.26 [1.05-1.51], p = 0.012). The protective association of HLA-A*02:01 was confirmed (OR = 0.72 [0.55-0.93], p = 0.013). None of the HLA-DQB1 alleles were associated with MS. Using a significance threshold of p < 0.01, outside the major histocompatibility complex region, 8 MS SNPs were also found to be associated with MS in African Americans.

CONCLUSION

MS genetic risk in African Americans only partially overlaps with that of Europeans and could explain the difference of MS prevalence between populations.

Blood RNA profiling in a large cohort of multiple sclerosis patients and healthy controls

Multiple sclerosis (MS) is the most common autoimmune disease of the central nervous system (CNS). It is characterized by the infiltration of autoreactive immune cells into the CNS, which target the myelin sheath, leading to the loss of neuronal function. Although it is accepted that MS is a multifactorial disorder with both genetic and environmental factors influencing its development and course, the molecular pathogenesis of MS has not yet been fully elucidated. Here, we studied the longitudinal gene expression profiles of whole-blood RNA from a cohort of 195 MS patients and 66 healthy controls. We analyzed these transcriptomes at both the individual transcript and the biological pathway level. We found 62 transcripts to be significantly up-regulated in MS patients; the expression of 11 of these genes was counter-regulated by interferon treatment, suggesting partial restoration of a 'healthy' gene expression profile. Global pathway analyses linked the proteasome and Wnt signaling to MS disease processes. Since genotypes from a subset of individuals were available, we were able to identify expression quantitative trait loci (eQTL), a number of which involved two genes of the MS gene signature. However, all these eQTL were also present in healthy controls. This study highlights the challenge posed by analyzing transcripts from whole blood and how these can be mitigated by using large, well-characterized cohorts of patients with longitudinal follow-up and multi-modality measurements.

Copy number variation in pediatric multiple sclerosis

BACKGROUND

Pediatric onset multiple sclerosis (MS) accounts for 2-4% of all MS. It is unknown whether the disease shares the same underlying pathophysiology found in adult patients or an extreme early onset phenotype triggered by distinct biological mechanisms. It has been hypothesized that copy number variations (CNVs) may result in extreme early onset diseases because CNVs can have major effects on many genes in large genomic regions.

OBJECTIVES AND METHODS

The objective of the current research was to identify CNVs, with a specific focus on de novo CNVs, potentially causing early onset MS by competitively hybridizing 30 white non-Hispanic pediatric MS patients with each of their parents via comparative genomic hybridization (CGH) analysis on the Agilent 1M CGH array.

RESULTS AND DISCUSSION

We identified 10 CNVs not overlapping with any CNV regions currently reported in the Database of Genomic Variants (DGV). Fifty-five putatively de novo CNVs were also identified: all but one common in the DGV. We found the single rare CNV was a private variation harboring the SACS gene. SACS mutations cause autosomal-recessive spastic ataxia of Charlevoix-Saguenay (ARSACS) disease. Additional clinical review revealed that the patient with the SACS gene CNV shared some features of both MS and ARSACS.

CONCLUSIONS

This is the first reported study analyzing pediatric MS CNVs. While not yielding causal variation in our initial pediatric dataset, our approach confirmed diagnosis of an ARSACS-like disease in addition to MS in the affected individual, which led to a more complete understanding of the patient's disease course and prognosis.

SNP-based analysis of the HLA locus in Japanese multiple sclerosis patients

Although several major histocompatibility complex (MHC)-wide single-nucleotide polymorphism (SNP) studies have been performed in populations of European descent, none have been performed in Asian populations. The objective of this study was to identify human leukocyte antigen (HLA) loci associated with multiple sclerosis (MS) in a Japanese population genotyped for 3534 MHC region SNPs. Using a logistic regression model, two SNPs (MHC Class III SNP rs422951 in the NOTCH4 gene and MHC Class II SNP rs3997849, susceptible alleles A and G, respectively) were independently associated with MS susceptibility (204 patients; 280 controls), two (MHC Class II SNP rs660895 and MHC Class I SNP rs2269704 in the NRM gene, susceptible alleles G and G, respectively) with aquaporin-4- (AQP4-) MS susceptibility (149 patients; 280 controls) and a single SNP (MHC Class II SNP rs1694112, susceptible allele G) was significant when contrasting AQP4+ against AQP4- patients. Haplotype analysis revealed a large susceptible association, likely DRB1*04 or a locus included in the DRB1*04 haplotype, with AQP4- MS, which excluded DRB1*15:01. This study is the largest study of the HLA's contribution to MS in Japanese individuals.

Aggregation of multiple sclerosis genetic risk variants in multiple and single case families

OBJECTIVE

Multiple sclerosis (MS) is a multifactorial neurologic disease characterized by modest but tractable heritability. Genome-wide association studies have identified and/or validated multiple polymorphisms in approximately 16 genes associated with susceptibility. We aimed at investigating the aggregation of genetic MS risk markers in individuals by comparing multiple- and single-case families.

METHODS

A weighted log-additive integrative approach termed MS genetic burden (MSGB) was used to account for the well-established genetic variants from previous association studies and meta-analyses. The corresponding genetic burden and its transmission was analyzed in 1,213 independent MS families (810 sporadic and 403 multicase families).

RESULTS

MSGB analysis demonstrated a higher aggregation of susceptibility variants in multicase compared to sporadic MS families. In addition, the aggregation of non-major histocompatibility complex single nucleotide polymorphisms depended neither on gender nor on the presence or absence of HLA-DRB1*15:01 alleles. Interestingly, although a greater MSGB in siblings of MS patients was associated with an increased risk of MS (odds ratio, 2.1; p = 0.001), receiver operating characteristic curves of MSGB differences between probands and sibs (area under the receiver operator curves, 0.57 [95% confidence interval, 0.53-0.61]) show that case-control status prediction of MS cannot be achieved with the currently available genetic data.

INTERPRETATION

The primary interest in the MSGB concept resides in its capacity to integrate cumulative genetic contributions to MS risk. This analysis underlines the high variability of family load with known common variants. This novel approach can be extended to other genetically complex diseases. Despite the emphasis on assembling large case-control datasets, multigenerational, multiaffected families remain an invaluable resource for advancing the understanding of the genetic architecture of complex traits.

Genetic variation in the odorant receptors family 13 and the mhc loci influence mate selection in a multiple sclerosis dataset

Background

When selecting mates, many vertebrate species seek partners with major histocompatibility complex (MHC) genes different from their own, presumably in response to selective pressure against inbreeding and towards MHC diversity. Attempts at replication of these genetic results in human studies, however, have reached conflicting conclusions.

Results

Using a multi-analytical strategy, we report validated genome-wide relationships between genetic identity and human mate choice in 930 couples of European ancestry. We found significant similarity between spouses in the MHC at class I region in chromosome 6p21, and at the odorant receptor family 13 locus in chromosome 9. Conversely, there was significant dissimilarity in the MHC class II region, near the HLA-DQA1 and -DQB1 genes. We also found that genomic regions with significant similarity between spouses show excessive homozygosity in the general population (assessed in the HapMap CEU dataset). Conversely, loci that were significantly dissimilar among spouses were more likely to show excessive heterozygosity in the general population.

Conclusions

This study highlights complex patterns of genomic identity among partners in unrelated couples, consistent with a multi-faceted role for genetic factors in mate choice behavior in human populations.

Multiple sclerosis susceptibility alleles in African Americans

Multiple sclerosis (MS) is an autoimmune demyelinating disease characterized by complex genetics and multifaceted gene-environment interactions. Compared to whites, African Americans have a lower risk for developing MS, but African Americans with MS have a greater risk of disability. These differences between African Americans and whites may represent differences in genetic susceptibility and/or environmental factors. SNPs from 12 candidate genes have recently been identified and validated with MS risk in white populations. We performed a replication study using 918 cases and 656 unrelated controls to test whether these candidate genes are also associated with MS risk in African Americans. CD6, CLEC16a, EVI5, GPC5, and TYK2 contained SNPs that are associated with MS risk in the African American data set. EVI5 showed the strongest association outside the major histocompatibility complex (rs10735781, OR=1.233, 95% CI=1.06-1.43, P-value=0.006). In addition, RGS1 seems to affect age of onset whereas TNFRSF1A seems to be associated with disease progression. None of the tested variants showed results that were statistically inconsistent with the effects established in whites. The results are consistent with shared disease genetic mechanisms among individuals of European and African ancestry.

Refining the association of MHC with multiple sclerosis in African Americans

Multiple sclerosis (MS) is a common demyelinating disease of the central nervous system mediated by autoimmune and neurodegenerative pathogenic mechanisms. Multiple genes account for its moderate heritability, but the only genetic region shown to have a large replicable effect on MS susceptibility is the major histocompatibility complex (MHC). Strong linkage disequilibrium (LD) across the MHC has made it difficult to fully characterize individual genetic contributions of this region to MS risk in previous studies. African Americans are at a lower risk for MS when compared with northern Europeans and Americans of European descent, but greater haplotypic diversity and distinct patterns of LD suggest that this population may be particularly informative for fine-mapping efforts. To examine the role of the MHC in African American MS, a case-control association study was performed with 499 African American MS patients and 750 African American controls that were genotyped for 6040 MHC region single nucleotide polymorphisms (SNPs). A replication data set consisting of 451 African American patients and 718 African American controls was genotyped for selected SNPs. Two MHC class II SNPs, rs2647040 and rs3135021, were significant in the replication cohort and partially tagged DRB1*15 alleles. Surprisingly, in comparison to similar studies of individuals of European descent, the MHC seems to play a smaller role in MS susceptibility in African Americans, consistent with pervasive genetic heterogeneity across ancestral groups, and may explain the difference in MS susceptibility between African Americans and individuals of European descent.

Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis

Identical (or more correctly 'monozygotic') twins are widely used to study the contributions of genetics and environment to human disease. A study that focused on three pairs of monozygotic twins, in which one twin had multiple sclerosis and the other did not, has brought the latest techniques of genome sequencing and analysis to this field, and incidentally published the first female human genome sequences. Full sequences were determined for one pair of twins, and for these and the other two pairs the mRNA transcriptome and epigenome sequences of CD4⁺ lymphocytes were determined. The striking result is that no genetic, epigenetic or transcriptome differences were found that explained why one twin had the disease and the other did not. Digging deeper into the data, eQTL (expression quantitative trait locus) mapping revealed tantalizing differences within twin pairs that merit closer examination. And some possible causes can be ruled out. Future work might usefully concentrate on studies of other cell types and epigenetic modifications.

Studies of identical twins are widely used to dissect the contributions of genes and the environment to human diseases. In multiple sclerosis, an autoimmune demyelinating disease, identical twins often show differences. This might suggest that environmental effects are most significant in this case, but genetic and epigenetic differences between identical twins have been described. Here, however, studies of identical twins show no evidence for genetic, epigenetic or transcriptome differences that could explain disease discordance.

Monozygotic or ‘identical’ twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis^{1,2,3,4,5,6,7,8}. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture^9,10,11,12. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4⁺ lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among ∼3.6 million single nucleotide polymorphisms (SNPs) or ∼0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of ∼19,000 genes in CD4⁺ T cells. Only 2 to 176 differences in the methylation of ∼2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to ∼800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

Copy number variation in African Americans

BACKGROUND

Copy number variants (CNVs) have been identified in several studies to be associated with complex diseases. It is important, therefore, to understand the distribution of CNVs within and among populations. This study is the first report of a CNV map in African Americans.

RESULTS

Employing a SNP platform with greater than 500,000 SNPs, a first-generation CNV map of the African American genome was generated using DNA from 385 healthy African American individuals, and compared to a sample of 435 healthy White individuals. A total of 1362 CNVs were identified within African Americans, which included two CNV regions that were significantly different in frequency between African Americans and Whites (17q21 and 15q11). In addition, a duplication was identified in 74% of DNAs derived from cell lines that was not present in any of the whole blood derived DNAs.

CONCLUSION

The Affymetrix 500 K array provides reliable CNV mapping information. However, using cell lines as a source of DNA may introduce artifacts. The duplication identified in high frequency in Whites and low frequency in African Americans on chromosome 17q21 reflects haplotype specific frequency differences between ancestral groups. The generation of the CNV map will be a valuable tool for identifying disease associated CNVs in African Americans.

Genome-wide association analysis of susceptibility and clinical phenotype in multiple sclerosis

Multiple sclerosis (MS), a chronic disorder of the central nervous system and common cause of neurological disability in young adults, is characterized by moderate but complex risk heritability. Here we report the results of a genome-wide association study performed in a 1000 prospective case series of well-characterized individuals with MS and group-matched controls using the Sentrix HumanHap550 BeadChip platform from Illumina. After stringent quality control data filtering, we compared allele frequencies for 551 642 SNPs in 978 cases and 883 controls and assessed genotypic influences on susceptibility, age of onset, disease severity, as well as brain lesion load and normalized brain volume from magnetic resonance imaging exams. A multi-analytical strategy identified 242 susceptibility SNPs exceeding established thresholds of significance, including 65 within the MHC locus in chromosome 6p21.3. Independent replication confirms a role for GPC5, a heparan sulfate proteoglycan, in disease risk. Gene ontology-based analysis shows a functional dichotomy between genes involved in the susceptibility pathway and those affecting the clinical phenotype.

Uncoupling the roles of HLA-DRB1 and HLA-DRB5 genes in multiple sclerosis

Genetic susceptibility to multiple sclerosis (MS) is associated with the MHC located on chromosome 6p21. This signal maps primarily to a 1-Mb region encompassing the HLA class II loci, and it segregates often with the HLA-DQB1*0602, -DQA1*0102, -DRB1*1501, -DRB5*0101 haplotype. However, the identification of the true predisposing gene or genes within the susceptibility haplotype has been handicapped by the strong linkage disequilibrium across the locus. African Americans have greater MHC haplotypic diversity and distinct patterns of linkage disequilibrium, which make this population particularly informative for fine mapping efforts. The purpose of this study was to establish the telomeric boundary of the HLA class II region affecting susceptibility to MS by assessing genetic association with the neighboring HLA-DRB5 gene as well as seven telomeric single nucleotide polymorphisms in a large, well-characterized African American dataset. Rare DRB5*null individuals were previously described in African populations. Although significant associations with both HLA-DRB1 and HLA-DRB5 loci were present, HLA-DRB1*1503 was associated with MS in the absence of HLA-DRB5, providing evidence for HLA-DRB1 as the primary susceptibility gene. Interestingly, the HLA-DRB5*null subjects appear to be at increased risk for developing secondary progressive MS. Thus, HLA-DRB5 attenuates MS severity, a finding consistent with HLA-DRB5's proposed role as a modifier in experimental autoimmune encephalomyelitis. Additionally, conditional haplotype analysis revealed a susceptibility signal at the class III AGER locus independent of DRB1. The data underscore the power of the African American MS dataset to identify disease genes by association in a region of high linkage disequilibrium.

Genetic variation in nitric oxide synthase 2A (NOS2A) and risk for multiple sclerosis

Multiple sclerosis (MS) is a chronic inflammatory disorder of the central nervous system with a strong genetic component. Variation in the major histocompatibility complex on chromosome 6p21, specifically the HLA-DRB1*15 haplotype, is the strongest genetic factor for MS, yet it is estimated to account for only a portion of risk for the disease. Previous evidence has implicated the nitric oxide synthase gene (NOS2A) encoding inducible NOS on chromosome 17q11 as a potential MS susceptibility gene. To determine whether variation in the NOS2A gene contributes to MS risk, we investigated a total of 50 polymorphisms within or flanking the locus for evidence of association using a comprehensive analytical strategy. A total of 6265 members from 1858 well-characterized MS families were utilized. No evidence for overtransmission of any individual single-nucleotide polymorphism allele or haplotype to the MS-affected individuals was observed. Furthermore, different transmission rates were not observed in either DRB1*15-positive or DRB1*15-negative family subgroups, or when extreme clinical outcomes characterizing disease progression were examined. The very largest study of NOS2A variation in MS, to date, excludes even a modest role for this locus in susceptibility.

Genome-wide pharmacogenomic analysis of the response to interferon beta therapy in multiple sclerosis

OBJECTIVE

To identify promising candidate genes linked to interindividual differences in the efficacy of interferon beta therapy. Recombinant interferon beta therapy is widely used to reduce disease activity in multiple sclerosis (MS). However, up to 50% of patients continue to have relapses and worsening disability despite therapy.

DESIGN

We used a genome-wide pharmacogenomic approach to identify single-nucleotide polymorphism (SNP) allelic differences associated with interferon beta therapy response.

SETTING

Four collaborating centers in the Mediterranean Basin. Data Coordination Center at the University of California, San Francisco.

PATIENTS

A cohort of 206 patients with relapsing-remitting MS followed up prospectively for 2 years after initiation of treatment.

INTERVENTION

DNA was pooled and hybridized to Affymetrix 100K GeneChips. Pooling schemes were designed to minimize confounding batch effects and increase confidence by technical replication.

MAIN OUTCOME MEASURES

Single-nucleotide polymorphism detection. Comparison of allelic frequencies between good responders and nonresponders to interferon beta therapy.

RESULTS

A multianalytical approach detected significant associations between several SNPs and treatment response, which were validated by individual DNA genotyping on an independent platform. After the validation stage was complete, 81 additional individuals were added to the analysis to increase power. We found that responders and nonresponders had significantly different genotype frequencies for SNPs located in many genes, including glypican 5, collagen type XXV alpha1, hyaluronan proteoglycan link protein, calpastatin, and neuronal PAS domain protein 3.

CONCLUSIONS

The reported results address the question of genetic heterogeneity in MS and the response to immunotherapy by analysis of the correlation between different genotypes and clinical response to interferon beta therapy. Many of the detected differences between responders and nonresponders were genes associated with ion channels and signal transduction pathways. The study also suggests that genetic variants in heparan sulfate proteoglycan genes may be of clinical interest in MS as predictors of the response to therapy. In addition to new insights into the mechanistic biology of interferon beta, these results help define the molecular basis of interferon beta therapy response heterogeneity.

The autoimmune disease-associated IL12B and IL23R polymorphisms in multiple sclerosis

Recently published genetic studies in psoriasis, inflammatory bowel disease, and ankylosing spondylitis identified significant associations with IL12B and IL23R polymorphisms. An important role for the IL-12/IL-23 pathway in multiple sclerosis (MS) is supported by immunologic studies in patients and animal models. To determine whether IL12B/IL23R disease-associated polymorphisms play a role in susceptibility to MS, we genotyped 910 MS-nuclear families, totaling 3132 individuals. Family-based association analysis was performed. There was no evidence of transmission distortion of any of the tested alleles in this data set.

Interleukin 7 receptor alpha chain (IL7R) shows allelic and functional association with multiple sclerosis

Multiple sclerosis is a demyelinating neurodegenerative disease with a strong genetic component. Previous genetic risk studies have failed to identify consistently linked regions or genes outside of the major histocompatibility complex on chromosome 6p. We describe allelic association of a polymorphism in the gene encoding the interleukin 7 receptor alpha chain (IL7R) as a significant risk factor for multiple sclerosis in four independent family-based or case-control data sets (overall P = 2.9 x 10(-7)). Further, the likely causal SNP, rs6897932, located within the alternatively spliced exon 6 of IL7R, has a functional effect on gene expression. The SNP influences the amount of soluble and membrane-bound isoforms of the protein by putatively disrupting an exonic splicing silencer.

A second major histocompatibility complex susceptibility locus for multiple sclerosis

Objective

Variation in the major histocompatibility complex (MHC) on chromosome 6p21 is known to influence susceptibility to multiple sclerosis with the strongest effect originating from the HLA-DRB1 gene in the class II region. The possibility that other genes in the MHC independently influence susceptibility to multiple sclerosis has been suggested but remains unconfirmed.

Methods

Using a combination of microsatellite, single nucleotide polymorphism, and human leukocyte antigen (HLA) typing, we screened the MHC in trio families looking for evidence of residual association above and beyond that attributable to the established DRB1*1501 risk haplotype. We then refined this analysis by extending the genotyping of classical HLA loci into independent cases and control subjects.

Results

Screening confirmed the presence of residual association and suggested that this was maximal in the region of the HLA-C gene. Extending analysis of the classical loci confirmed that this residual association is partly due to allelic heterogeneity at the HLA-DRB1 locus, but also reflects an independent effect from the HLA-C gene. Specifically, the HLA-C*05 allele, or a variant in tight linkage disequilibrium with it, appears to exert a protective effect (p = 3.3 × 10⁻⁵).

Interpretation

Variation in the HLA-C gene influences susceptibility to multiple sclerosis independently of any effect attributable to the nearby HLA-DRB1 gene. Ann Neurol 2007

Genome-wide network analysis reveals the global properties of IFN-beta immediate transcriptional effects in humans

IFN-beta effectively controls clinical exacerbations and magnetic resonance imaging activity in most multiple sclerosis patients. However, its mechanism of action has not been yet fully elucidated. In this study we used DNA microarrays to analyze the longitudinal transcriptional profile of blood cells within a week of IFN-beta administration. Using differential expression and gene ontology analyses we found evidence of a general decrease in the cellular activity of T lymphocytes resembling the endogenous antiviral response of IFNs. In contrast, most of the differentially expressed genes (DEGs) from untreated individuals were involved in cellular physiological processes. We then used mutual information (MI) to build networks of coregulated genes in both treated and untreated individuals. Interestingly, the connectivity distribution (k) of networks generated with high MI values displayed scale-free properties. Conversely, the observed k for networks generated with suboptimal MI values approximated a Poisson distribution, suggesting that MI captures biologically relevant interactions. Gene networks from individuals treated with IFN-beta revealed a tight core of immune- and apoptosis-related genes associated with higher values of MI. In contrast, networks obtained from untreated individuals primarily reflected cellular housekeeping functions. Finally, we trained a neural network to reverse engineer the directionality of the main interactions observed at the biological process level. This is the first study that incorporates network analysis to investigate gene regulation in response to a therapeutic drug in humans. Implications of this method in the creation of personalized models of response to therapy are discussed.

Quantitative longitudinal analysis of T cell receptor repertoire expression in HIV-infected patients on antiretroviral and interleukin-2 therapy

We have developed a single-step reverse transcriptase kinetic PCR assay (kRT-PCR) to accurately determine the expression of each of the 24 TCRbetaV gene families in CD8(+) cells. We analyzed the long-term effects of highly active antiretroviral therapy (HAART) on the stability of the CD8(+) T cell receptor (TCR) repertoire in a cohort of 15 treated and 10 untreated individuals diagnosed with human immunodeficiency virus (HIV) infection. The CD4(+) TCR repertoire was studied in a second cohort receiving interleukin-2 infusions in addition to HAART. Analysis was based on kinetic (quantitative) reverse-transcription PCR (kRT-PCR) of the TCR variable B gene (TCRbetaV). Expression of each of the 24 Vbeta families was assessed at baseline immediately after infection and following initiation of HAART at 2, 4, 12, 24, and up to 192 weeks in 24-week intervals. Statistically significant family-specific expression changes were observed between treated and untreated individuals for 10 TCRbetaV families. Overall, when compared to untreated patients, a more stable expression of TCR genes was observed for HAART-treated individuals. Interestingly, this difference did not correlate with either CD4 or CD8 counts, which follow the expected curves for treated and untreated patients. When we applied our quantitative analysis to IL-2-treated patients we observed a rapid polyclonal activation of the repertoire. These results suggest that homeostasis in the T cell receptor repertoire is more robust in those patients who stay on HAART for a long time and confirm the polyclonal stimulating capacity of IL-2.

Complex gene-gene interactions in multiple sclerosis: a multifactorial approach reveals associations with inflammatory genes

The complex inheritance involved in multiple sclerosis (MS) risk has been extensively investigated, but our understanding of MS genetics remains rudimentary. In this study, we explore 51 single nucleotide polymorphisms (SNPs) in 36 candidate genes from the inflammatory pathway and test for gene-gene interactions using complementary case-control, discordant sibling pair, and trio family study designs. We used a sample of 421 carefully diagnosed MS cases and 96 unrelated, healthy controls; discordant sibling pairs from 146 multiplex families; and 275 trio families. We used multifactor dimensionality reduction to explore gene-gene interactions. Based on our analyses, we have identified several statistically significant models including both main effect models and two-locus, three-locus, and four-locus epistasis models that predict MS disease risk with between approximately 61% and 85% accuracy. These results suggest that significant epistasis, or gene-gene interactions, may exist even in the absence of statistically significant individual main effects.

Multifactor dimensionality reduction reveals gene-gene interactions associated with multiple sclerosis susceptibility in African Americans

Multiple sclerosis (MS) is a common disease of the central nervous system characterized by inflammation, myelin loss, gliosis, varying degrees of axonal pathology, and progressive neurological dysfunction. Multiple sclerosis exhibits many of the characteristics that distinguish complex genetic disorders including polygenic inheritance and environmental exposure risks. Here, we used a highly efficient multilocus genotyping assay representing variation in 34 genes associated with inflammatory pathways to explore gene-gene interactions and disease susceptibility in a well-characterized African-American case-control MS data set. We applied the multifactor dimensionality reduction (MDR) test to detect epistasis, and identified single-IL4R(Q576R)- and three-IL4R(Q576R), IL5RA(-80), CD14(-260)- locus association models that predict MS risk with 75-76% accuracy (P<0.01). These results demonstrate the importance of exploring both main effects and gene-gene interactions in the study of complex diseases.

Linkage and association with the NOS2A locus on chromosome 17q11 in multiple sclerosis

A large body of research supports a multifactorial cause in multiple sclerosis (MS), with an underlying genetic susceptibility likely acting in concert with undefined environmental exposures. Here, we used a highly efficient multilocus genotyping assay to study single nucleotide polymorphisms representing variation in 34 genes from inflammatory pathways in a well-characterized MS familial data set. Evidence of transmission distortion was present for several polymorphisms. Results for the NOS2A locus (exon 10 C/T, D346D) on chromosome 17q11 remained significant after correction for multiple testing and were reproduced in a second independent African American MS data set. In addition, linkage to a NOS2A promoter region polymorphism, (CCTTT)(n), was present in a third data set of multicase MS families. Our results provide strong evidence for linkage and association to a new candidate disease gene on chromosome 17q11 in MS and suggest that variation within NOS2A or a nearby locus contributes to disease susceptibility.