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.

PTPRC (CD45) is not associated with the development of multiple sclerosis in U.S. patients.

A C-->G nucleotide transition in exon 4 of PTPRC (encoding protein-tyrosine phosphatase receptor-type C, also known as CD45) was recently reported to be genetically associated with the development of multiple sclerosis (MS). We performed an extensive evaluation of this polymorphism using large family-based and case-control comparisons. Overall, we observed no evidence of genetic association between the PTPRC polymorphism and MS susceptibility or disease course.

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.

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.

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.

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.

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.

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.

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.

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.

Multiple sclerosis: genomic rewards.

A large body of immunologic, epidemiologic, and genetic data indicate that tissue injury in multiple sclerosis (MS) results from an abnormal immune response to one or more myelin antigens that develops in genetically susceptible individuals after exposure to an as-yet undefined causal agent. The genetic component of MS etiology is believed to result from the action of several genes of moderate effect. The incomplete penetrance of MS susceptibility alleles probably reflects interactions with other genes, post transcriptional regulatory mechanisms, and significant nutritional and environmental influences. Equally significant, it is also likely that genetic heterogeneity exists, meaning that specific genes influence susceptibility and pathogenesis in some affects but not in others. Results in multiplex MS families confirm the genetic importance of the MHC region in conferring susceptibility of MS. Susceptibility may be mediated by the class II genes themselves (DR, DQ or both), related to the known function of these molecules in the normal immune response, e.g. antigen binding and presentation and T cell repertoire determination. The possibility that other genes in the MHC or the telomeric region of the MHC are responsible for the observed genetic effect cannot be excluded. The data also indicate that although the MHC region plays a significant role in MS susceptibility, much of the genetic effect in MS remains to be explained. Some loci may be involved in the initial pathogenic events, while others could influence the development and progression of the disease. The past few years have seen real progress in the development of laboratory and analytical approaches to study non-Mendelian complex genetic disorders and in defining the pathological basis of demyelination, setting the stage for the final characterization of the genes involved in MS susceptibility and pathogenesis. Their identification and characterization is likely to define the basic etiology of the disease, improve risk assessment and influence therapeutics.

Sequence variation in the transforming growth factor-beta1 (TGFB1) gene and multiple sclerosis susceptibility.

Genome screenings in multiple sclerosis (MS) have identified multiple susceptibility regions supporting a polygenic model for this disease. Evidence for linkage was consistently observed at ch.19q13 suggesting the presence of an MS gene(s) in this region. Several interesting candidate genes are encoded within this region, including transforming growth factor-beta 1 (TGFB1) and interleukin-11 (IL11). Both are multifunctional cytokines with significant and well-characterized immunomodulatory properties. We performed a comprehensive evaluation of common polymorphisms within the TGFB1 and IL11 loci and three closely flanking microsatellite markers (D19S421, CEA, D19S908) in 161 stringently ascertained and clinically characterized MS multiplex families using tests of both linkage (lod score, sib-pair analysis) and association (pedigree disequilibrium test or PDT). Patients and families were stratified by HLA-DR2 status to search for two-locus interactions. Suggestive evidence for linkage and association to CEA (lod score = 1.25, theta = 0.20, p = 0.015, respectively), located 0.4 cM from TGFB1, was observed in DR2 positive families only. Distinct clinical phenotypes were also examined and an association between a TGFB1 haplotype and a mild disease course was present (p = 0.008), raising the possibility that TGFB1 or a nearby locus may influence disease expression.

Evolution of the CCR5 Delta32 mutation based on haplotype variation in Jewish and Northern European population samples.

The chemokine receptor 5 (CCR5) serves as a fusion cofactor for macrophage-tropic strains of HIV-1. In addition, CCR5 has been shown to mediate the entry of poxviruses into target cells. Individuals homozygous for the Delta32 deletion-mutation have no surface expression of CCR5 and are highly protected against HIV-1 infection. To gain insights into the evolution of the mutation in modern populations, the relatively high frequency of the Delta32-ccr5 allele in some European and Jewish populations is explored here by examining haplotypes of 3p21.3 constructed of five polymorphic marker loci surrounding CCR5. By sampling Ashkenazi, non-Ashkenazi and non-Jewish populations, we utilize the natural experiment that occurred as a consequence of the Jewish Diaspora, and demonstrate that a single mutation was responsible for all copies of Delta32. This mutation must have moved from Northern European populations to the Ashkenazi Jews where evidence suggests that Delta32 carriers of both groups were favored by repeated occurrence of epidemic small pox beginning in the 8th century AD.

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.