Influence of Sex on the Microbiota of the Human Face.

The role of the microbiota in health and disease has long been recognized and, so far, the cutaneous microbiota in humans has been widely investigated. The research regarded mainly the microbiota variations between body districts and disease skin states (i.e., atopic dermatitis, psoriasis, acne). In fact, relatively little information is available about the composition of the healthy skin microbiota. The cosmetic industry is especially interested in developing products that maintain and/or improve a healthy skin microbiota. Therefore, in the present work, the authors chose to investigate in detail the structure and composition of the basal bacterial community of the face. Ninety-six cheek samples (48 women and 48 men) were collected in the same season and the same location in central northern Italy. Bacterial DNA was extracted, the 16S rDNA gene was amplified by PCR, the obtained amplicons were subjected to next generation sequencing. The principal members of the community were identified at the genus level, and statistical analyses showed significant variations between the two sexes. This study identified abundant members of the facial skin microbiota that were rarely reported before in the literature and demonstrated the differences between male and female microbiota in terms of both community structure and composition.

A Scoping Review on Body Fluid Biomarkers for Prognosis and Disease Activity in Patients with Multiple Sclerosis.

Multiple sclerosis (MS) is a complex demyelinating disease of the central nervous system, presenting with different clinical forms, including clinically isolated syndrome (CIS), which is a first clinical episode suggestive of demyelination. Several molecules have been proposed as prognostic biomarkers in MS. We aimed to perform a scoping review of the potential use of prognostic biomarkers in MS clinical practice. We searched MEDLINE up to 25 November 2021 for review articles assessing body fluid biomarkers for prognostic purposes, including any type of biomarkers, cell types and tissues. Original articles were obtained to confirm and detail the data reported by the review authors. We evaluated the reliability of the biomarkers based on the sample size used by various studies. Fifty-two review articles were included. We identified 110 molecules proposed as prognostic biomarkers. Only six studies had an adequate sample size to explore the risk of conversion from CIS to MS. These confirm the role of oligoclonal bands, immunoglobulin free light chain and chitinase CHI3L1 in CSF and of serum vitamin D in the prediction of conversion from CIS to clinically definite MS. Other prognostic markers are not yet explored in adequately powered samples. Serum and CSF levels of neurofilaments represent a promising biomarker.

A multi-step genomic approach prioritized TBKBP1 gene as relevant for multiple sclerosis susceptibility.

BACKGROUND: Over 200 genetic loci have been associated with multiple sclerosis (MS) explaining ~ 50% of its heritability, suggesting that additional mechanisms may account for the "missing heritability" phenomenon. OBJECTIVE: To analyze a large cohort of Italian individuals to identify markers associated with MS with potential functional impact in the disease. METHODS: We studied 2571 MS and 3234 healthy controls (HC) of continental Italian origin. Discovery phase included a genome wide association study (1727 MS, 2258 HC), with SNPs selected according to their association in the Italian cohort only or in a meta-analysis of signals with a cohort of European ancestry (4088 MS, 7144 HC). Top associated loci were then tested in two Italian cohorts through array-based genotyping (903 MS, 884 HC) and pool-based target sequencing (588 MS, 408 HC). Finally, functional prioritization through conditional eQTL and mQTL has been performed. RESULTS: Top associated signals overlap with already known MS loci on chromosomes 3 and 17. Three SNPs (rs4267364, rs8070463, rs67919208), all involved in the regulation of TBKBP1, were prioritized to be functionally relevant. CONCLUSIONS: No evidence of novel signal of association with MS specific for the Italian continental population has been found; nevertheless, two MS loci seems to play a relevant role, raising the interest to further investigations for TBKBP1 gene.

The impact of lifetime coffee and tea loads on Multiple Sclerosis severity.

BACKGROUND AND AIMS: The association between lifestyle factors and Multiple Sclerosis (MS) disease severity and progression has been investigated to a lesser extent compared with susceptibility to the disease. We aimed to assess the impact of lifetime coffee and tea consumption on MS severity. METHODS: Design: cross-sectional study. Two hundred and eight patients (139 females and 69 males) consecutively recruited at the Department of Neurology in Novara, Italy were asked about their lifetime consumption of coffee and tea. The lifetime intensity of consumption (cups/day) was estimated as the weighted sum of the mean number of standard cups drunk per day at different ages. A measure of cumulative lifetime load of the exposure was expressed in terms of cup-years. Disease severity was estimated by the Multiple Sclerosis Severity Score (MSSS). HLA-DRB1∗15 and HLA-A∗02 genotyping was performed in 167 patients. RESULTS: The MSSS was not associated with the status of coffee or tea consumer, or the amount of cups/day or cup-years. The Odds Ratios (OR) for falling in the upper tertile of the MSSS distribution was 1.30 (95% Confidence Interval (CI): 0.47-3.58) for coffee consumers of 1-3 cups/day and 1.14 (95%CI: 0.33-3.95) for 4-8 cups/day vs. non-consumers. The OR was 0.69 (95%CI: 0.35-1.34) for tea consumers vs. non-consumers. However, heavy consumers of coffee (4-8 cups/day) more frequently had a progressive form than small consumers (1-3 cups/day) and non-consumers (19% vs. 14% vs. 0%), and had a significantly higher age at MS onset (36.6 ± 10.3; 31.5 ± 9.5; 28.6 ± 8.1 years, p = 0.001). Although not reaching statistical significance, coffee consumers positive for HLA-A∗02 had a six-fold risk of being in the worst tertile compared to never consumers, whereas the risk was only 1.3 for coffee consumers negative for the same allele. CONCLUSIONS: Coffee or tea intake is not associated with different severity of MS. However, we cannot exclude a possible effect of higher doses of coffee for the subgroup of progressive patients.

Targeted Next-Generation Sequencing for the Identification of Genetic Predictors of Radiation-Induced Late Skin Toxicity in Breast Cancer Patients: A Preliminary Study.

Normal tissue radiosensitivity is thought to be influenced by an individual's genetic background. However, the specific genetic variants underlying the risk of late skin reactions following radiotherapy for breast cancer remain elusive. To unravel the genetic basis for radiation-induced late skin toxicity, we carried out targeted next-generation sequencing of germline DNA samples from 48 breast cancer patients with extreme late skin toxicity phenotypes, consisting of 24 cases with grade 2-3 subcutaneous fibrosis and/or grade 2-3 telangiectasia (LENT-SOMA scales) and 24 controls with grade 0 fibrosis and grade 0 telangiectasia. In this exploratory study, a total of five single-nucleotide variants (SNVs) located in three genes (TP53, ERCC2, and LIG1) reached nominal levels of statistical significance (p < 0.05). In the replication study, which consisted of an additional 45 cases and 192 controls, none of the SNVs identified by targeted NGS achieved nominal replication. Nevertheless, TP53 rs1042522 (G > C, Pro72Arg) in the replication cohort had an effect (OR per C allele: 1.52, 95%CI: 0.82-2.83, p = 0.186) in the same direction as in the exploratory cohort (OR per C allele: 4.70, 95%CI: 1.51-14.6, p = 0.007) and was found be nominally associated to the risk of radiation-induced late skin toxicity in the overall combined cohort (OR per C allele: 1.79, 95%CI: 1.06-3.02, p = 0.028). These results raise the possibility of an association between TP53 rs1042522 and risk of radiation-induced late skin toxicity in breast cancer patients; however, large replication studies are warranted for conclusive evidence.

An Investigation of the Role of Common and Rare Variants in a Large Italian Multiplex Family of Multiple Sclerosis Patients.

Known multiple sclerosis (MS) susceptibility variants can only explain half of the disease's estimated heritability, whereas low-frequency and rare variants may partly account for the missing heritability. Thus, here we sought to determine the occurrence of rare functional variants in a large Italian MS multiplex family with five affected members. For this purpose, we combined linkage analysis and next-generation sequencing (NGS)-based whole exome and whole genome sequencing (WES and WGS, respectively). The genetic burden attributable to known common MS variants was also assessed by weighted genetic risk score (wGRS). We found a significantly higher burden of common variants in the affected family members compared to that observed among sporadic MS patients and healthy controls (HCs). We also identified 34 genes containing at least one low-frequency functional variant shared among all affected family members, showing a significant enrichment in genes involved in specific biological processes-particularly mRNA transport-or neurodegenerative diseases. Altogether, our findings point to a possible pathogenic role of different low-frequency functional MS variants belonging to shared pathways. We propose that these rare variants, together with other known common MS variants, may account for the high number of affected family members within this MS multiplex family.

Genomic and functional evaluation of TNFSF14 in multiple sclerosis susceptibility.

Among multiple sclerosis (MS) susceptibility genes, the strongest non-human leukocyte antigen (HLA) signal in the Italian population maps to the TNFSF14 gene encoding LIGHT, a glycoprotein involved in dendritic cell (DC) maturation. Through fine-mapping in a large Italian dataset (4,198 patients with MS and 3,903 controls), we show that the TNFSF14 intronic SNP rs1077667 is the primarily MS-associated variant in the region. Expression quantitative trait locus (eQTL) analysis indicates that the MS risk allele is significantly associated with reduced TNFSF14 messenger RNA levels in blood cells, which is consistent with the allelic imbalance in RNA-Seq reads (P < 0.0001). The MS risk allele is associated with reduced levels of TNFSF14 gene expression (P < 0.01) in blood cells from 84 Italian patients with MS and 80 healthy controls (HCs). Interestingly, patients with MS are lower expressors of TNFSF14 compared to HC (P < 0.007). Individuals homozygous for the MS risk allele display an increased percentage of LIGHT-positive peripheral blood myeloid DCs (CD11c+, P = 0.035) in 37 HCs, as well as in in vitro monocyte-derived DCs from 22 HCs (P = 0.04). Our findings suggest that the intronic variant rs1077667 alters the expression of TNFSF14 in immune cells, which may play a role in MS pathogenesis.

Validation of an Algorithm to Detect Multiple Sclerosis Cases in Administrative Health Databases in Piedmont (Italy): An Application to the Estimate of Prevalence by Age and Urbanization Level.

INTRODUCTION: Italy is considered a high-risk country for multiple sclerosis (MS). Exploiting electronic health archives (EHAs) is highly useful to continuously monitoring the prevalence of the disease, as well as the care delivered to patients and its outcomes. The aim of this study was to validate an EHA-based algorithm to identify MS patients, suitable for epidemiological purposes, and to estimate MS prevalence in Piedmont (North Italy). METHODS: MS cases were identified, in the period between January 1, 2012 and December 31, 2017, linking data from 4 different sources: hospital discharges, drug prescriptions, exemptions from co-payment to health care, and long-term care facilities. Sensitivity of the algorithm was tested through record linkage with a cohort of 656 neurologist-confirmed MS cases; specificity was tested with a cohort of 2,966,293 residents presumably not affected by MS. Undercount was estimated by a capture-recapture method. We calculated crude, and age- and gender-specific prevalence. We also calculated age-adjusted prevalence by level of urbanization of the municipality of residence. RESULTS: On December 31, 2017, the algorithm identified 8,850 MS cases. Sensitivity was 95.9%, specificity was 99.97%, and the estimated completeness of ascertainment was 91.9%. The overall prevalence, adjusted for undercount, was 152 per 100,000 among men and 286 among women; it increased with increasing age and reached its peak value in the 45- to 54-year class, followed by a progressive reduction. The age-adjusted prevalence of residents in cities was 15% higher than in those living in the countryside. DISCUSSION/CONCLUSION: We validated an algorithm based on EHAs to identify cases of MS for epidemiological use. The prevalence of MS, adjusted for undercount, was among the highest in Italy. We also found that the prevalence was higher in highly urbanized areas.

Contribution of Rare and Low-Frequency Variants to Multiple Sclerosis Susceptibility in the Italian Continental Population.

Genome-wide association studies identified over 200 risk loci for multiple sclerosis (MS) focusing on common variants, which account for about 50% of disease heritability. The goal of this study was to investigate whether low-frequency and rare functional variants, located in MS-established associated loci, may contribute to disease risk in a relatively homogeneous population, testing their cumulative effect (burden) with gene-wise tests. We sequenced 98 genes in 588 Italian patients with MS and 408 matched healthy controls (HCs). Variants were selected using different filtering criteria based on allelic frequency and in silico functional impacts. Genes showing a significant burden (n = 17) were sequenced in an independent cohort of 504 MS and 504 HC. The highest signal in both cohorts was observed for the disruptive variants (stop-gain, stop-loss, or splicing variants) located in EFCAB13, a gene coding for a protein of an unknown function (p < 10-4). Among these variants, the minor allele of a stop-gain variant showed a significantly higher frequency in MS versus HC in both sequenced cohorts (p = 0.0093 and p = 0.025), confirmed by a meta-analysis on a third independent cohort of 1298 MS and 1430 HC (p = 0.001) assayed with an SNP array. Real-time PCR on 14 heterozygous individuals for this variant did not evidence the presence of the stop-gain allele, suggesting a transcript degradation by non-sense mediated decay, supported by the evidence that the carriers of the stop-gain variant had a lower expression of this gene (p = 0.0184). In conclusion, we identified a novel low-frequency functional variant associated with MS susceptibility, suggesting the possible role of rare/low-frequency variants in MS as reported for other complex diseases.

Periostin Circulating Levels and Genetic Variants in Patients with Non-Alcoholic Fatty Liver Disease.

Circulating periostin has been suggested as a possible biomarker in non-alcoholic fatty liver disease (NAFLD) in Asian studies. In the present study, we aimed to test its still controversial relevance in a Caucasian population. In patients with histologically-proven NAFLD (N. = 74; 10 with hepatocellular carcinoma, HCC) plasma periostin concentrations were analyzed. POSTN haplotype analysis was based on rs9603226, rs3829365, and rs1029728. Hepatitis C patients (N. = 81, 7 HCC) and healthy subjects (N. = 27) were used as controls. The median plasma periostin concentration was 11.6 ng/mL without differences amongst groups; it was not influenced by age, liver fibrosis or steatosis. However, possession of haplotype two (rs9603226 = G, rs3829365 = C, rs1028728 = A) was associated with lower circulating periostin compared to other haplotypes. Moreover, periostin was higher in HCC patients. At multivariate analysis, HCC remained the only predictor of high periostin. In conclusion, plasma periostin concentrations in Caucasians NAFLD patients are not influenced by the degree of liver disease, but are significantly higher in HCC. Genetically-determined differences may account for some of the variability. These data suggest extreme caution in predicting a possible future role of periostin antagonists as a rational therapeutic alternative for NAFLD, but show a potential periostin role in the management of NAFLD-associated HCC.

Author Correction: Next Generation Sequencing of Pooled Samples: Guideline for Variants' Filtering.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Analysis of the GCG repeat length in NIPA1 gene in C9orf72-mediated ALS in a large Italian ALS cohort.

Amyotrophic lateral sclerosis (ALS) is characterized by degeneration of upper and lower motor neurons. The hexanucleotide repeat expansion in C9orf72 gene (C9orf72-HRE) is the most frequent genetic cause of ALS. Since many ALS pedigrees showed incomplete penetrance, several genes have been analyzed as possible modifiers. Length of the GCG repeat tract in NIPA1 (non-imprinted in Prader-Willi/Angelman syndrome 1) gene has been recently investigated as a possible modifier factor for C9orf72-HRE patients with contrasting findings. To disclose the possible role of NIPA1 GCG repeat length as modifier of the disease risk in C9orf72-HRE carriers, we analyzed a large cohort of 532 Italian ALS cases enriched in C9orf72-HRE carriers (172 cases) and 483 Italian controls. This sample size is powered (92% power, p = 0.05) to replicate the modifier effect observed in literature. We did not observe higher frequency of NIPA1 long alleles (> 8 GCG) in C9orf72-HRE carriers (3.5%) compared with C9orf72-HRE negative patients (4.1%) and healthy controls (5%). For the latter comparison, we meta-analyzed our data with currently available literature data, and no statistically significant effect was observed (p = 0.118). In conclusion, we did not confirm a role of NIPA1 repeat length as a modifier of the C9orf72 ALS disease risk.

Characterization of the c9orf72 GC-rich low complexity sequence in two cohorts of Italian and Turkish ALS cases.

Large expansions of a noncoding GGGGCC repeat in the C9orf72 gene are the main cause of amyotrophic lateral sclerosis (ALS). The GGGGCC repeat is contiguous with another GC-rich region. Recent studies reported a significantly higher frequency of insertions/deletions within the GC-rich region in patients carrying the GGGGCC expansion. A GTGGT motif comprised within the GC-rich region, which joins two 100% GC sequences, was frequently deleted, supporting the hypothesis that these deletions could make the region more prone to slippage and pathological expansion. To confirm this hypothesis, we sequenced the GC-rich region adjacent the GGGGCC repeat in ALS patients, 116 C9orf72 expansion carriers, 219 non-carriers, and 223 healthy controls, from Italian and Turkish cohorts. Deletions were significantly more frequent in C9orf72 expansion carriers (6%) compared to non-carrier ALS patients (0.46%, OR =14.00, 95% CI =1.71-306.59, p = 0.003), to controls (0%, OR =16.29, 95% CI =2.12-725.99, p = 4.86 × 10-4) and to the whole cohort of non-carriers (0.2%, OR =28.51, 95% CI =3.47-618.91, p = 9.58 × 10-5). Among expansion carriers, deletions with or without the GTGGT motif were equally distributed (4 vs. 3). The frequency of insertions was not statistically different between C9orf72 expansion carriers and any other group including the whole cohort of non-carriers (p = 0.439, Fisher's exact test). Our data confirmed the association between deletions within GC-rich region and the GGGGCC expansion in Italian and Turkish cases, although we did not confirm a role of the GTGGT element deletion. Further studies will be therefore necessary to assess the causal relationships between contiguous deletions of the GC-rich region and the GGGGCC expansion.

HLA alleles modulate EBV viral load in multiple sclerosis.

BACKGROUND: The etiopathology of multiple sclerosis (MS) is believed to include genetic and environmental factors. Human leukocyte antigen (HLA) alleles, in particular,  are associated with disease susceptibility, whereas Epstein Barr Virus (EBV) infection has long been suspected to play a role in disease pathogenesis. The aim of the present study is to evaluate correlations between HLA alleles and EBV infection in MS. METHODS: HLA alleles, EBV viral load (VL) and serum anti-EBV antibody titers were evaluated in EBV-seropositive MS patients (N = 117) and age- and sex-matched healthy controls (HC; N = 89). RESULTS: Significantly higher DNA viral loads (p = 0.048) and EBNA-1 antibody titer (p = 0.0004) were seen in MS compared to HC. EBV VL was higher in HLA-B*07+ (p = 0.02) and HLA-DRB1*15+ (p = 0.02) MS patients, whereas it was lower in HLA-A*02+ (p = 0.04) subjects. EBV VL was highest in HLA-A*02-/B*07+/DRB1*15+ patients and lowest in HLA-A*A02+/B*07-/DRB1*15- individuals (p < 0.0001). HLA-B*07 resulted the most associated allele to EBV VL after multiple regression analysis considering altogether the three alleles, (p = 0.0001). No differences were observed in anti-EBV antibody titers in relationship with HLA distribution. CONCLUSIONS: Host HLA-B*07 allele influence EBV VL in MS. As HLA-class I molecules present antigens to T lymphocytes and initiate immune response against viruses, these results could support a role for EBV in MS.

Overexpression of the Cytokine BAFF and Autoimmunity Risk.

BACKGROUND: Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug-targetable pathways. METHODS: Using case-control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus-specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence-based fine mapping, cross-population and cross-phenotype analyses, and gene-expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. RESULTS: A variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease-risk allele was also associated with up-regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion-deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up-regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. CONCLUSIONS: A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.).

Next Generation Sequencing of Pooled Samples: Guideline for Variants' Filtering.

Sequencing large number of individuals, which is often needed for population genetics studies, is still economically challenging despite falling costs of Next Generation Sequencing (NGS). Pool-seq is an alternative cost- and time-effective option in which DNA from several individuals is pooled for sequencing. However, pooling of DNA creates new problems and challenges for accurate variant call and allele frequency (AF) estimation. In particular, sequencing errors confound with the alleles present at low frequency in the pools possibly giving rise to false positive variants. We sequenced 996 individuals in 83 pools (12 individuals/pool) in a targeted re-sequencing experiment. We show that Pool-seq AFs are robust and reliable by comparing them with public variant databases and in-house SNP-genotyping data of individual subjects of pools. Furthermore, we propose a simple filtering guideline for the removal of spurious variants based on the Kolmogorov-Smirnov statistical test. We experimentally validated our filters by comparing Pool-seq to individual sequencing data showing that the filters remove most of the false variants while retaining majority of true variants. The proposed guideline is fairly generic in nature and could be easily applied in other Pool-seq experiments.

Variations in the high-mobility group-A2 gene (HMGA2) are associated with idiopathic short stature.

BACKGROUND: Several association studies confirmed high-mobility group-A2 gene (HMGA2) polymorphisms as the most relevant variants contributing to height variability. Animal models and deletions in humans suggest that alterations of HMGA2 might be relevant in causing short stature. Together, these observations led us to investigate the involvement of HMGA2 in idiopathic short stature (ISS) through an association study and a mutation screening. METHODS: We conducted an association study (155 ISS patients and 318 normal stature controls) with three HMGA2 single-nucleotide polymorphisms (SNPs) (SNPs rs1042725, rs7968682, and rs7968902) using a TaqMan-based assay. The patients were then analyzed by direct sequencing and multiplex ligation-dependent probe amplification (MLPA) to detect point mutations and genomic micro-rearrangements. RESULTS: Considering a recessive model, an OR value >1 was observed for genotypes rs7968682 TT (Odds ratio (OR) = 1.72, confidence interval (CI): 1.14-2.58) and rs1042725 TT (OR = 1.51, CI: 1.00-2.28) in accordance to the effect exhibited by the single alleles in the general population. None of the patients carried possibly causative HMGA2 mutations. CONCLUSION: Besides the already known role in determining variability in human height, HMGA2 polymorphisms also contribute to susceptibility to ISS. Moreover, we here report the first mutation screening performed in ISS concluding that HMGA2 has not a significant impact on the monogenic form of ISS.

The burden of multiple sclerosis variants in continental Italians and Sardinians.

BACKGROUND: Recent studies identified > 100 non-HLA (human leukocyte antigen) multiple sclerosis (MS) susceptibility variants in Northern European populations, but their role in Southern Europeans is largely unexplored. OBJECTIVE: We aimed to investigate the cumulative impact of those variants in two Mediterranean populations: Continental Italians and Sardinians. METHODS: We calculated four weighted Genetic Risk Scores (wGRS), using up to 102 non-HLA MS risk variants and 5 HLA MS susceptibility markers in 1691 patients and 2194 controls from continental Italy; and 2861 patients and 3034 controls from Sardinia. We then assessed the differences between populations using Nagelkerke's R(2) and the area under the Receiver Operating Characteristic (ROC) curves. RESULTS: As expected, the genetic burden (mean wGRS value) was significantly higher in MS patients than in controls, in both populations. Of note, the burden was significantly higher in Sardinians. Conversely, the proportion of variability explained and the predictive power were significantly higher in continental Italians. Notably, within the Sardinian patients, we also observed a significantly higher burden of non-HLA variants in individuals who do not carry HLA risk alleles. CONCLUSIONS: The observed differences in MS genetic burden between the two Mediterranean populations highlight the need for more genetic studies in South Europeans, to further expand the knowledge of MS genetics.

Genetic variants are major determinants of CSF antibody levels in multiple sclerosis.

Immunological hallmarks of multiple sclerosis include the production of antibodies in the central nervous system, expressed as presence of oligoclonal bands and/or an increased immunoglobulin G index-the level of immunoglobulin G in the cerebrospinal fluid compared to serum. However, the underlying differences between oligoclonal band-positive and -negative patients with multiple sclerosis and reasons for variability in immunoglobulin G index are not known. To identify genetic factors influencing the variation in the antibody levels in the cerebrospinal fluid in multiple sclerosis, we have performed a genome-wide association screen in patients collected from nine countries for two traits, presence or absence of oligoclonal bands (n = 3026) and immunoglobulin G index levels (n = 938), followed by a replication in 3891 additional patients. We replicate previously suggested association signals for oligoclonal band status in the major histocompatibility complex region for the rs9271640*A-rs6457617*G haplotype, correlated with HLA-DRB1*1501, and rs34083746*G, correlated with HLA-DQA1*0301 (P comparing two haplotypes = 8.88 × 10(-16)). Furthermore, we identify a novel association signal of rs9807334, near the ELAC1/SMAD4 genes, for oligoclonal band status (P = 8.45 × 10(-7)). The previously reported association of the immunoglobulin heavy chain locus with immunoglobulin G index reaches strong evidence for association in this data set (P = 3.79 × 10(-37)). We identify two novel associations in the major histocompatibility complex region with immunoglobulin G index: the rs9271640*A-rs6457617*G haplotype (P = 1.59 × 10(-22)), shared with oligoclonal band status, and an additional independent effect of rs6457617*G (P = 3.68 × 10(-6)). Variants identified in this study account for up to 2-fold differences in the odds of being oligoclonal band positive and 7.75% of the variation in immunoglobulin G index. Both traits are associated with clinical features of disease such as female gender, age at onset and severity. This is the largest study population so far investigated for the genetic influence on antibody levels in the cerebrospinal fluid in multiple sclerosis, including 6950 patients. We confirm that genetic factors underlie these antibody levels and identify both the major histocompatibility complex and immunoglobulin heavy chain region as major determinants.