Breastfeeding is associated with reduced risk of multiple sclerosis in males, predominantly among HLA-DRB1*15:01 carriers.

BACKGROUND: Breastfeeding as an infant appears protective against later development of some autoimmune diseases, but research into its influence on multiple sclerosis (MS) risk has yielded inconclusive results. OBJECTIVE: We investigated the possible impact of breastfeeding on MS risk. METHODS: We used two population-based case-control studies comprising 3670 cases and 6737 matched controls. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI) for association between MS and exposure to prolonged breastfeeding (4 months or longer) versus reduced breastfeeding (less than 4 months). A meta-analysis of case-control studies that assessed the impact of breastfeeding on MS risk among women and men was conducted. RESULTS: Prolonged breastfeeding was associated with reduced MS risk among men (OR 0.7, 95% CI 0.5-0.9) but not among women (OR 0.9, 95% CI 0.8-1.1). Among men, a synergistic effect was observed between HLA-DRB1*15:01 carrier status and reduced breastfeeding. CONCLUSIONS: Findings from the current study add to accumulating evidence that breastfeeding may be a modifiable protective factor for reducing the risk of MS in offspring. When possible, mothers should be supported to breastfeed their infants; however, the mechanism of a sex-specific biologic effect of breastfeeding on MS risk is unclear.

Genetic risk and longitudinal disease activity in systemic lupus erythematosus using targeted maximum likelihood estimation.

Systemic lupus erythematous (SLE) is a chronic autoimmune disease associated with genetic and environmental risk factors. However, the extent to which genetic risk is causally associated with disease activity is unknown. We utilized longitudinal-targeted maximum likelihood estimation to estimate the causal association between a genetic risk score (GRS) comprising 41 established SLE variants and clinically important disease activity as measured by the validated Systemic Lupus Activity Questionnaire (SLAQ) in a multiethnic cohort of 942 individuals with SLE. We did not find evidence of a clinically important SLAQ score difference (>4.0) for individuals with a high GRS compared with those with a low GRS across nine time points after controlling for sex, ancestry, renal status, dialysis, disease duration, treatment, depression, smoking and education, as well as time-dependent confounding of missing visits. Individual single-nucleotide polymorphism (SNP) analyses revealed that 12 of the 41 variants were significantly associated with clinically relevant changes in SLAQ scores across time points eight and nine after controlling for multiple testing. Results based on sophisticated causal modeling of longitudinal data in a large patient cohort suggest that individual SLE risk variants may influence disease activity over time. Our findings also emphasize a role for other biological or environmental factors.

High consumption of coffee is associated with decreased multiple sclerosis risk; results from two independent studies.

OBJECTIVE: Previous studies on consumption of caffeine and risk of multiple sclerosis (MS) have yielded inconclusive results. We aimed to investigate whether consumption of coffee is associated with risk of MS. METHODS: Using two population-representative case-control studies (a Swedish study comprising 1620 cases and 2788 controls, and a US study comprising 1159 cases and 1172 controls), participants with different habits of coffee consumption based on retrospective data collection were compared regarding risk of MS, by calculating ORs with 95% CIs. Logistic regression models were adjusted for a broad range of potential confounding factors. RESULTS: Compared with those who reported no coffee consumption, the risk of MS was substantially reduced among those who reported a high consumption of coffee exceeding 900 mL daily (OR 0.70 (95% CI 0.49 to 0.99) in the Swedish study, and OR 0.69 (95% CI 0.50 to 0.96) in the US study). Lower odds of MS with increasing consumption of coffee were observed, regardless of whether coffee consumption at disease onset or 5 or 10 years prior to disease onset was considered. CONCLUSIONS: In accordance with studies in animal models of MS, high consumption of coffee may decrease the risk of developing MS. Caffeine, one component of coffee, has neuroprotective properties, and has been shown to suppress the production of proinflammatory cytokines, which may be mechanisms underlying the observed association. However, further investigations are needed to determine whether exposure to caffeine underlies the observed association and, if so, to evaluate its mechanisms of action.

Annotation of functional variation within non-MHC MS susceptibility loci through bioinformatics analysis.

There is a strong and complex genetic component to multiple sclerosis (MS). In addition to variation in the major histocompatibility complex (MHC) region on chromosome 6p21.3, 110 non-MHC susceptibility variants have been identified in Northern Europeans, thus far. The majority of the MS-associated genes are immune related; however, similar to most other complex genetic diseases, the causal variants and biological processes underlying pathogenesis remain largely unknown. We created a comprehensive catalog of putative functional variants that reside within linkage disequilibrium regions of the MS-associated genic variants to guide future studies. Bioinformatics analyses were also conducted using publicly available resources to identify plausible pathological processes relevant to MS and functional hypotheses for established MS-associated variants.

MHC associations with clinical and autoantibody manifestations in European SLE.

Systemic lupus erythematosus (SLE) is a clinically heterogeneous disease affecting multiple organ systems and characterized by autoantibody formation to nuclear components. Although genetic variation within the major histocompatibility complex (MHC) is associated with SLE, its role in the development of clinical manifestations and autoantibody production is not well defined. We conducted a meta-analysis of four independent European SLE case collections for associations between SLE sub-phenotypes and MHC single-nucleotide polymorphism genotypes, human leukocyte antigen (HLA) alleles and variant HLA amino acids. Of the 11 American College of Rheumatology criteria and 7 autoantibody sub-phenotypes examined, anti-Ro/SSA and anti-La/SSB antibody subsets exhibited the highest number and most statistically significant associations. HLA-DRB1*03:01 was significantly associated with both sub-phenotypes. We found evidence of associations independent of MHC class II variants in the anti-Ro subset alone. Conditional analyses showed that anti-Ro and anti-La subsets are independently associated with HLA-DRB1*0301, and that the HLA-DRB1*03:01 association with SLE is largely but not completely driven by the association of this allele with these sub-phenotypes. Our results provide strong evidence for a multilevel risk model for HLA-DRB1*03:01 in SLE, where the association with anti-Ro and anti-La antibody-positive SLE is much stronger than SLE without these autoantibodies.

The rs4774 CIITA missense variant is associated with risk of systemic lupus erythematosus.

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor required for human leukocyte antigens (HLA) class II MHC-restricted antigen presentation. MHC genes, including the HLA class II DRB1*03:01 allele, are strongly associated with systemic lupus erythematosus (SLE). Recently the rs4774 CIITA missense variant (+1632G/C) was reported to be associated with susceptibility to multiple sclerosis. In the current study, we investigated CIITA, DRB1*03:01 and risk of SLE using a multi-stage analysis. In stage 1, 9 CIITA variants were tested in 658 cases and 1363 controls (N=2021). In stage 2, rs4774 was tested in 684 cases and 2938 controls (N=3622). We also performed a meta-analysis of the pooled 1342 cases and 4301 controls (N=5643). In stage 1, rs4774(*)C was associated with SLE (odds ratio (OR)=1.24, 95% confidence interval (95% CI)=1.07-1.44, P=4.2 × 10(-3)). Similar results were observed in stage 2 (OR=1.16, 95% CI=1.02-1.33, P=8.5 × 10(-3)) and the meta-analysis of the combined data set (OR=1.20, 95% CI=1.09-1.33, P(meta)=2.5 × 10(-4)). In all three analyses, the strongest evidence for association between rs4774(*)C and SLE was present in individuals who carried at least one copy of DRB1*03:01 (P(meta)=1.9 × 10(-3)). Results support a role for CIITA in SLE, which appears to be stronger in the presence of DRB1*03:01.

CIITA is not associated with risk of developing rheumatoid arthritis.

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor regulating genes required for human leukocyte antigen (HLA) class II MHC-restricted antigen presentation. MHC genes, particularly HLA class II, are strongly associated with risk of developing rheumatoid arthritis (RA). Given the strong biological relationship between CIITA and HLA class II genes, a comprehensive investigation of CIITA variation in RA was conducted. This study tested 31 CIITA single-nucleotide polymorphisms in 2542 RA cases and 3690 controls (N=6232). All individuals were of European ancestry, as determined by ancestry informative genetic markers. No evidence for association between CIITA variation and RA was observed after a correction for multiple testing was applied. This is the largest study to fully characterize common genetic variation in CIITA, including an assessment of haplotypes. Results exclude even a modest role for common CIITA polymorphisms in susceptibility to RA.

A non-synonymous SNP within membrane metalloendopeptidase-like 1 (MMEL1) is associated with multiple sclerosis.

Several single-nucleotide polymorphism (SNP) genome-wide association studies (GWASs) have been completed in multiple sclerosis (MS). Follow-up studies of the variants with the most promising rankings, especially when supplemented by informed candidate gene selection, have proven to be extremely successful. In this study we report the results of a multi-stage replication analysis of the putatively associated SNPs identified in the Wellcome Trust Case Control Consortium non-synonymous SNP (nsSNP) screen. In total, the replication sample consisted of 3444 patients and 2595 controls. A combined analysis of the nsSNP screen and replication data provides evidence implicating a novel additional locus, rs3748816 in membrane metalloendopeptidase-like 1 (MMEL1; odds ratio=1.16, P=3.54 × 10⁻6) in MS susceptibility.

A candidate gene study of CLEC16A does not provide evidence of association with risk for anti-CCP-positive rheumatoid arthritis.

CLEC16A, a putative immunoreceptor, was recently established as a susceptibility locus for type I diabetes and multiple sclerosis. Subsequently, associations between CLEC16A and rheumatoid arthritis (RA), Addison's disease and Crohn's disease have been reported. A large comprehensive and independent investigation of CLEC16A variation in RA was pursued. This study tested 251 CLEC16A single-nucleotide polymorphisms in 2542 RA cases (85% anti-cyclic citrullinated peptide (anti-CCP) positive) and 2210 controls (N=4752). All individuals were of European ancestry, as determined by ancestry informative genetic markers. No evidence for significant association between CLEC16A variation and RA was observed. This is the first study to fully characterize common genetic variation in CLEC16A including assessment of haplotypes and gender-specific effects. The previously reported association between RA and rs6498169 was not replicated. Results show that CLEC16A does not have a prominent function in susceptibility to anti-CCP-positive RA.

Genetic variation within the HLA class III influences T1D susceptibility conferred by high-risk HLA haplotypes.

Human leukocyte antigen (HLA) class II DRB1 and DQB1 represent the major type I diabetes (T1D) genetic susceptibility loci; however, other genes in the HLA region are also involved in T1D risk. We analyzed 1411 pedigrees (2865 affected individuals) from the type I diabetes genetics consortium genotyped for HLA classical loci and for 12 single-nucleotide polymorphisms (SNPs) in the class III region previously shown to be associated with T1D in a subset of 886 pedigrees. Using the transmission disequilibrium test, we compared the proportion of SNP alleles transmitted from within the high-risk DR3 and DR4 haplotypes to affected offspring. Markers rs4151659 (mapping to CFB) and rs7762619 (mapping 5' of LTA) were the most strongly associated with T1D on DR3 (P=1.2 x 10(-9) and P=2 x 10(-12), respectively) and DR4 (P=4 x 10(-15) and P=8 x 10(-8), respectively) haplotypes. They remained significantly associated after stratifying individuals in analyses for B*1801, A*0101-B*0801, DPB1*0301, DPB1*0202, DPB1*0401 or DPB1*0402. Rs7762619 and rs4151659 are in strong linkage disequilibrium (LD) (r(2)=0.82) with each other, but a joint analysis showed that the association for each SNP was not solely because of LD. Our data support a role for more than one locus in the class III region contributing to risk of T1D.

Supervised machine learning and logistic regression identifies novel epistatic risk factors with PTPN22 for rheumatoid arthritis.

Investigating genetic interactions (epistasis) has proven difficult despite the recent advances of both laboratory methods and statistical developments. With no 'best' statistical approach available, combining several analytical methods may be optimal for detecting epistatic interactions. Using a multi-stage analysis that incorporated supervised machine learning and methods of association testing, we investigated epistatic interactions with a well-established genetic factor (PTPN22 1858T) in a complex autoimmune disease (rheumatoid arthritis (RA)). Our analysis consisted of four principal stages: Stage I (data reduction)-identifying candidate chromosomal regions in 292 affected sibling pairs, by predicting PTPN22 concordance using multipoint identity-by-descent probabilities and a supervised machine learning algorithm (Random Forests); Stage II (extension analysis)-testing detailed genetic data within candidate chromosomal regions for epistasis with PTPN22 1858T in 677 cases and 750 controls using logistic regression; Stage III (replication analysis)-confirmation of epistatic interactions in 947 cases and 1756 controls; Stage IV (combined analysis)-a pooled analysis including all 1624 RA cases and 2506 control subjects for final estimates of effect size. A total of seven replicating epistatic interactions were identified. SNP variants within CDH13, MYO3A, CEP72 and near WFDC1 showed significant evidence for interaction with PTPN22, affecting susceptibility to RA.

Analysis of maternal-offspring HLA compatibility, parent-of-origin and non-inherited maternal effects for the classical HLA loci in type 1 diabetes.

AIM: Type 1 diabetes (T1D) is a complex trait for which variation in the classical human leucocyte antigen (HLA) loci within the Major Histocompatibility Complex (MHC) significantly influences disease risk. To date, HLA class II DR-DQ genes confer the strongest known genetic effect in T1D. HLA loci may also influence T1D through additional inherited or non-inherited effects. Evidence for the role of increased maternal-offspring HLA compatibility, and both parent-of-origin (POO) and non-inherited maternal HLA (NIMA) effects in autoimmune disease has been previously established. The current study tested hypotheses that classical HLA loci influence T1D through these mechanisms, in addition to genetic transmission of particular risk alleles. METHODS: The Type 1 Diabetes Genetics Consortium (T1DGC) cohort was of European descent and consisted of 2271 affected sib-pair families (total n = 11 023 individuals). Class I genes HLA-A, Cw and B, and class II genes HLA-DRB1, DQA1, DQB1, DPA1 and DPB1 were studied. The pedigree disequilibrium test was used to examine transmission of HLA alleles to individuals with T1D. Conditional logistic regression was used to model compatibility relationships between mother-offspring and father-offspring for all HLA loci. POO and NIMA effects were investigated by comparing frequencies of maternal and paternal transmitted and non-transmitted HLA alleles for each locus. Analyses were also stratified by gender of T1D-affected offspring. RESULTS: Strong associations were observed for all classical HLA loci except for DPA1, as expected. Compatibility differences between mother-offspring and father-offspring were not observed for any HLA loci. Furthermore, POO and NIMA HLA effects influencing T1D were not present. CONCLUSIONS: Maternal-offspring HLA compatibility, POO and NIMA effects for eight classical HLA loci were investigated. Results suggest that these HLA-related effects are unlikely to play a major role in the development of T1D.

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.

The MHC2TA -168A/G polymorphism and risk for rheumatoid arthritis: a meta-analysis of 6861 patients and 9270 controls reveals no evidence for association.

BACKGROUND: An association between major histocompatibility complex (MHC) genes, particularly those within the class II HLA region, and rheumatoid arthritis (RA) is well established, and accounts for an estimated 30% of the genetic component in RA. The MHC class II transactivator gene (MHC2TA) on chromosome 16p13 has recently emerged as the most important transcription factor regulating genes required for class II MHC-restricted antigen presentation. Previous studies of a promoter region polymorphism (-168A/G, rs3087456) in the MHC2TA gene and RA have yielded conflicting results. OBJECTIVE: To assess the association of the MHC2TA -168A/G polymorphism (rs3087456) and risk for RA by meta-analysis. METHODS: Meta-analysis was performed for 6861 patients with RA and 9270 controls from 10 case-control studies. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated for each study. Summary ORs and 95% CIs were calculated for random effects models. RESULTS: No effect was observed for the G risk allele (OR 1.02, 95% CI 0.93 to 1.12, p = 0.70) or the GG risk genotype (OR 1.14, 95% CI 0.95 to 1.36, p = 0.16). CONCLUSIONS: Our results indicate that the MHC2TA -168A/G polymorphism (rs3087456) is not associated with RA yet underscore the importance of including shared epitope allele carrier status, secondary phenotypes and more complete characterisation of MHC2TA variation in future studies.

Familial effects on the clinical course of multiple sclerosis.

BACKGROUND: Familial factors influence susceptibility to multiple sclerosis (MS) but it is unknown whether there are additional effects on the natural history of the disease. METHOD: We evaluated 1,083 families with > or =2 first-degree relatives with MS for concordance of age at onset, clinical course, and disease severity and investigated transmission patterns of these clinical features in affected parent-child pairs. RESULTS: There is concordance for age at onset for all families (correlation coefficient 0.14; p < 0.001), as well as for affected siblings (correlation coefficient 0.15; p < 0.001), and affected parent-child pairs (correlation coefficient 0.12; p = 0.03) when each is evaluated separately. Concordance for year of onset is present among affected siblings (correlation coefficient 0.18; p < 0.001) but not the parent-child group (correlation coefficient 0.08; p = 0.15). The clinical course is similar between siblings (kappa 0.12; p < 0.001) but not affected parents and their children (kappa -0.04; p = 0.09). This influence on the natural history is present in all clinical subgroups of relapsing-remitting, and primary and secondary progressive MS, reflecting a familial effect on episodic and progressive phases of the disease. There is no concordance for disease severity within any of the considered family groups (correlation coefficients: all families analyzed together, 0.02, p = 0.53; affected sibling group, 0.02, p = 0.61; affected parent-child group, 0.02, p = 0.69). Furthermore, there are no apparent transmission patterns of any of the investigated clinical features in affected parent-child pairs and no evidence for anticipation or effects of genetic loading. CONCLUSION: Familial factors do not significantly affect eventual disease severity. However, they increase the probability of a progressive clinical course, either from onset or after a phase of relapsing remitting disease. The familial effect is more likely to reflect genetic than environmental conditions. The results are relevant for counseling patients and have implications for the design of studies seeking to identify factors that influence the natural history of the disease.

Men transmit MS more often to their children vs women: the Carter effect.

OBJECTIVE: Multiple sclerosis (MS) is approximately twice as common among women as men. If men have greater physiologic resistance to MS, they might theoretically require stronger genetic predisposition than women to overcome this resistance. In this circumstance, men would be expected to transmit the disease more often to their children, a phenomenon known as the Carter effect. The authors evaluated whether the Carter effect is present in MS. METHODS: The authors studied 441 children (45 with definite MS) of an affected father or mother (197 families of interest) from 3598 individuals in 206 multiplex pedigrees. The authors compared transmission of MS from affected men with transmission from affected women. RESULTS: Fathers with MS transmitted the disease to their children more often (transmitted: 18, not transmitted: 99) than mothers with MS (transmitted: 27, not transmitted: 296) (p = 0.032; OR: 1.99, 95% CI: 1.05, 3.77). Adjusting for both the sex of the affected child and multiple transmissions from a single affected parent, the sex of the affected parent remained as an independent risk factor for transmission of MS to children, fathers transmitting more often than mothers (p = 0.036; OR: 2.21, 95% CI: 1.05, 4.63). CONCLUSIONS: The authors have demonstrated the Carter effect in multiple sclerosis (MS). These observations may be explained by greater genetic loading in men that leads to relative excess paternal vs maternal transmission. Linkage analysis in genetic studies of MS may be more informative if patrilineal transmission were given additional weighting.

Allelic association of sequence variants in the herpes virus entry mediator-B gene (PVRL2) with the severity of multiple sclerosis.

Discrepant findings have been reported regarding an association of the apolipoprotein E (APOE) gene with the clinical course of multiple sclerosis (MS). To resolve these discrepancies, we examined common sequence variation in six candidate genes residing in a 380-kb genomic region surrounding and including the APOE locus for an association with MS severity. We genotyped at least three polymorphisms in each of six candidate genes in 1,540 Caucasian MS families (729 single-case and multiple-case families from the United States, 811 single-case families from the UK). By applying the quantitative transmission/disequilibrium test to a recently proposed MS severity score, the only statistically significant (P=0.003) association with MS severity was found for an intronic variant in the Herpes Virus Entry Mediator-B Gene PVRL2. Additional genotyping extended the association to a 16.6 kb block spanning intron 1 to intron 2 of the gene. Sequencing of PVRL2 failed to identify variants with an obvious functional role. In conclusion, the analysis of a very large data set suggests that genetic polymorphisms in PVRL2 may influence MS severity and supports the possibility that viral factors may contribute to the clinical course of MS, consistent with previous reports.

APOE epsilon variation in multiple sclerosis susceptibility and disease severity: some answers.

BACKGROUND: Previous studies have examined the role of APOE variation in multiple sclerosis (MS), but have lacked the statistical power to detect modest genetic influences on risk and disease severity. The meta- and pooled analyses presented here utilize the largest collection, to date, of MS cases, controls, and families genotyped for the APOE epsilon polymorphism. METHODS: Studies of MS and APOE were identified by searches of PubMed, Biosis, Web of Science, Cochrane Review, and Embase. When possible, authors were contacted for individual genotype data. Meta-analyses of MS case-control data and family-based analyses were performed to assess the association of APOE epsilon genotype with disease risk. Pooled analyses of MS cases were also performed to assess the influence of APOE epsilon genotype on disease severity. RESULTS: A total of 22 studies (3,299 MS cases and 2,532 controls) were available for meta-analysis. No effect of epsilon2 or epsilon4 status on MS risk was observed (summary OR 1.14, 95% CI 0.96-1.34 and OR 0.89, 95% CI 0.78-1.01). Results obtained from analyses of APOE genotypes in 1,279 MS families were also negative (p = 0.61). Finally, results from pooled analyses of 4,048 MS cases also argue strongly that APOE epsilon status does not distinguish a relapsing-remitting from primary progressive disease course, or influence disease severity, as measured by the Expanded Disability Status Scale and disease duration. CONCLUSION: Overall, these findings do not support a role for APOE in multiple sclerosis, and underscore the importance of using large sample sizes to detect modest genetic effects, particularly in studies of genotype-phenotype relationships.

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.