Cross-Trait Mendelian Randomization Study to Investigate Whether Migraine Is a Risk Factor for Multiple Sclerosis.

BACKGROUND AND OBJECTIVES: Migraine is common among people with multiple sclerosis (MS), but the reasons for this are unknown. We tested 3 hypothesized mechanisms for this observed comorbidity, including migraine is a risk factor of MS, genetic variants are shared between the conditions, and migraine is because of MS. METHODS: Data were from 2 sources: publicly available summary statistics from genome-wide association studies of MS (N = 115,748) and migraine (N = 375,752 and N = 361,141) and a case-control study of MS recruited from the Kaiser Permanente Northern California Health Plan (N = 1,991). For the latter participants, migraine status was ascertained using a validated electronic health record migraine probability algorithm or self-report. Using the public summary statistics, we used 2-sample Mendelian randomization to test whether a migraine genetic instrumental variable was associated with MS. We used linkage disequilibrium score regression and LOGODetect to ascertain whether MS and migraine shared genetic variants across the genome and regionally. Using the Northern California MS cohort, we used logistic regression to identify whether people with both MS and migraine had different odds of clinical characteristics (e.g., age at MS onset, Perceived Deficits Questionnaire, and depression) or MS-specific risk factors (e.g., body mass index, smoking status, and infectious mononucleosis status) compared with people with MS without migraine. RESULTS: We did not find evidence supporting migraine as a causal risk factor of MS (p = 0.29). We did, however, identify 4 major histocompatibility complex (MHC) loci shared between MS and migraine. Among the Northern California MS cohort, 774 (39%) experienced migraine. People with both MS and migraine from this cohort were more likely to ever smoke (odds ratio [OR] = 1.30, 95% CI: 1.08-1.57), have worse self-reported cognitive deficits (OR = 1.04, 95% CI: 1.02-1.06), and ever experience depression (OR = 1.48, 95% CI: 1.22-1.80). DISCUSSION: Our findings do not support migraine as a causal risk factor of MS. Several genetic variants, particularly in the MHC, may account for some of the overlap. It seems likely that migraine within the context of MS is because of MS. Identifying what increases the risk of migraine within MS might lead to an improved treatment and quality of life.

Global expression and CpG methylation analysis of primary endothelial cells before and after TNFa stimulation reveals gene modules enriched in inflammatory and infectious diseases and associated DMRs.

Endothelial cells are a primary site of leukocyte recruitment during inflammation. An increase in tumor necrosis factor-alpha (TNFa) levels as a result of infection or some autoimmune diseases can trigger this process. Several autoimmune diseases are now treated with TNFa inhibitors. However, genomic alterations that occur as a result of TNF-mediated inflammation are not well understood. To investigate molecular targets and networks resulting from increased TNFa, we measured DNA methylation and gene expression in 40 human umbilical vein endothelial cell primary cell lines before and 24 hours after stimulation with TNFa via microarray. Weighted gene co-expression network analysis identified 15 gene groups (modules) with similar expression correlation patterns; four modules showed a strong association with TNFa treatment. Genes in the top TNFa-associated module were all up-regulated, had the highest proportion of hypomethylated regions, and were associated with 136 Disease Ontology terms, including autoimmune/inflammatory, infectious and cardiovascular diseases, and cancers. They included chemokines CXCL1, CXCL10 and CXCL8, and genes associated with autoimmune diseases including HLA-C, DDX58, IL4, NFKBIA and TNFAIP3. Cardiovascular and metabolic disease genes, including APOC1, ACLY, ELOVL6, FASN and SCD, were overrepresented in a module that was not associated with TNFa treatment. Of 223 hypomethylated regions identified, several were in promoters of autoimmune disease GWAS loci (ARID5B, CD69, HDAC9, IL7R, TNIP1 and TRAF1). Results reveal specific gene groups acting in concert in endothelial cells, delineate those driven by TNFa, and establish their relationship to DNA methylation changes, which has strong implications for understanding disease etiology and precision medicine approaches.

miRNA contributions to pediatric-onset multiple sclerosis inferred from GWAS.

OBJECTIVE: Onset of multiple sclerosis (MS) occurs in childhood for approximately 5% of cases (pediatric MS, or ped-MS). Epigenetic influences are strongly implicated in MS pathogenesis in adults, including the contribution from microRNAs (miRNAs), small noncoding RNAs that affect gene expression by binding target gene mRNAs. Few studies have specifically examined miRNAs in ped-MS, but individuals developing MS at an early age may carry a relatively high burden of genetic risk factors, and miRNA dysregulation may therefore play a larger role in the development of ped-MS than in adult-onset MS. This study aimed to look for evidence of miRNA involvement in ped-MS pathogenesis. METHODS: GWAS results from 486 ped-MS cases and 1362 controls from the U.S. Pediatric MS Network and Kaiser Permanente Northern California membership were investigated for miRNA-specific signals. First, enrichment of miRNA-target gene network signals was evaluated using MIGWAS software. Second, SNPs in miRNA genes and in target gene binding sites (miR-SNPs) were tested for association with ped-MS, and pathway analysis was performed on associated target genes. RESULTS: MIGWAS analysis showed that miRNA-target gene signals were enriched in GWAS (P = 0.038) and identified 39 candidate biomarker miRNA-target gene pairs, including immune and neuronal signaling genes. The miR-SNP analysis implicated dysregulation of miRNA binding to target genes in five pathways, mainly involved in immune signaling. INTERPRETATION: Evidence from GWAS suggests that miRNAs play a role in ped-MS pathogenesis by affecting immune signaling and other pathways. Candidate biomarker miRNA-target gene pairs should be further studied for diagnostic, prognostic, and/or therapeutic utility.

BMI1 enhancer polymorphism underlies chromosome 10p12.31 association with childhood acute lymphoblastic leukemia.

Genome-wide association studies of childhood acute lymphoblastic leukemia (ALL) have identified regions of association at PIP4K2A and upstream of BMI1 at chromosome 10p12.31-12.2. The contribution of both loci to ALL risk and underlying functional variants remain to be elucidated. We carried out single nucleotide polymorphism (SNP) imputation across chromosome 10p12.31-12.2 in Latino and non-Latino white ALL cases and controls from two independent California childhood leukemia studies, and additional Genetic Epidemiology Research on Aging study controls. Ethnicity-stratified association analyses were performed using logistic regression, with meta-analysis including 3,133 cases (1,949 Latino, 1,184 non-Latino white) and 12,135 controls (8,584 Latino, 3,551 non-Latino white). SNP associations were identified at both BMI1 and PIP4K2A. After adjusting for the lead PIP4K2A SNP, genome-wide significant associations remained at BMI1, and vice-versa (pmeta < 10-10 ), supporting independent effects. Lead SNPs differed by ethnicity at both peaks. We sought functional variants in tight linkage disequilibrium with both the lead Latino SNP among Admixed Americans and lead non-Latino white SNP among Europeans. This pinpointed rs11591377 (pmeta = 2.1 x 10-10 ) upstream of BMI1, residing within a hematopoietic stem cell enhancer of BMI1, and which showed significant preferential binding of the risk allele to MYBL2 (p = 1.73 x 10-5 ) and p300 (p = 1.55 x 10-3 ) transcription factors using binomial tests on ChIP-Seq data from a SNP heterozygote. At PIP4K2A, we identified rs4748812 (pmeta = 1.3 x 10-15 ), which alters a RUNX1 binding motif and demonstrated chromosomal looping to the PIP4K2A promoter. Fine-mapping chromosome 10p12 in a multi-ethnic ALL GWAS confirmed independent associations and identified putative functional variants upstream of BMI1 and at PIP4K2A.

Hypomethylation of CYP2E1 and DUSP22 Promoters Associated With Disease Activity and Erosive Disease Among Rheumatoid Arthritis Patients.

OBJECTIVE: Epigenetic modifications have previously been associated with rheumatoid arthritis (RA). In this study, we aimed to determine whether differential DNA methylation in peripheral blood cell subpopulations is associated with any of 4 clinical outcomes among RA patients. METHODS: Peripheral blood samples were obtained from 63 patients in the University of California, San Francisco RA cohort (all satisfied the American College of Rheumatology classification criteria; 57 were seropositive for rheumatoid factor and/or anti-cyclic citrullinated protein). Fluorescence-activated cell sorting was used to separate the cells into 4 immune cell subpopulations (CD14+ monocytes, CD19+ B cells, CD4+ naive T cells, and CD4+ memory T cells) per individual, and 229 epigenome-wide DNA methylation profiles were generated using Illumina HumanMethylation450 BeadChips. Differentially methylated positions and regions associated with the Clinical Disease Activity Index score, erosive disease, RA Articular Damage score, Sharp score, medication at time of blood draw, smoking status, and disease duration were identified using robust regression models and empirical Bayes variance estimators. RESULTS: Differential methylation of CpG sites associated with clinical outcomes was observed in all 4 cell types. Hypomethylated regions in the CYP2E1 and DUSP22 gene promoters were associated with active and erosive disease, respectively. Pathway analyses suggested that the biologic mechanisms underlying each clinical outcome are cell type-specific. Evidence of independent effects on DNA methylation from smoking, medication use, and disease duration were also identified. CONCLUSION: Methylation signatures specific to RA clinical outcomes may have utility as biomarkers or predictors of exposure, disease progression, and disease severity.

Correlates of Prenatal and Early-Life Tobacco Smoke Exposure and Frequency of Common Gene Deletions in Childhood Acute Lymphoblastic Leukemia.

Tobacco smoke exposure has been associated with risk of childhood acute lymphoblastic leukemia (ALL). Understanding the relationship between tobacco exposures and specific mutations may yield etiologic insights. We carried out a case-only analysis to explore whether prenatal and early-life tobacco smoke exposure influences the formation of leukemogenic genomic deletions. Somatic copy number of 8 genes frequently deleted in ALL (CDKN2A, ETV6, IKZF1, PAX5, RB1, BTG1, PAR1 region, and EBF1) was assessed in 559 pretreatment tumor samples from the California Childhood Leukemia Study. Parent and child's passive tobacco exposure was assessed using interview-assisted questionnaires as well as DNA methylation in aryl-hydrocarbon receptor repressor (AHRR), a sentinel epigenetic biomarker of exposure to maternal smoking during pregnancy. Multivariable Poisson regressions were used to test the association between the smoking exposures and total number of deletions. Deletion burden varied by subtype, with a lower frequency in high-hyperdiploid and higher frequency in ETV6-RUNX1 fusion ALL. The total number of deletions per case was positively associated with tobacco smoke exposure, in particular for maternal ever-smoking (ratio of means, RM, 1.31; 95% CI, 1.08-1.59), maternal smoking during pregnancy (RM, 1.48; 95% CI, 1.12-1.94), and during breastfeeding (RM, 2.11; 95% CI, 1.48-3.02). The magnitude of association with maternal ever-smoking was stronger in male children compared with females (Pinteraction = 0.04). The total number of deletions was also associated with DNA methylation at the AHRR epigenetic biomarker (RM, 1.32; 95% CI, 1.02-1.69). Our results suggest that prenatal and early-life tobacco smoke exposure increase the frequency of somatic deletions in children who develop ALL. Cancer Res; 77(7); 1674-83. ©2017 AACR.

Causal Effect of Genetic Variants Associated With Body Mass Index on Multiple Sclerosis Susceptibility.

Multiple sclerosis (MS) is an autoimmune disease with both genetic and environmental risk factors. Recent studies indicate that childhood and adolescent obesity double the risk of MS, but this association may reflect unmeasured confounders rather than causal effects of obesity. We used separate-sample Mendelian randomization to estimate the causal effect of body mass index (BMI) on susceptibility to MS. Using data from non-Hispanic white members of the Kaiser Permanente Medical Care Plan of Northern California (KPNC) (2006-2014; 1,104 cases of MS and 10,536 controls) and a replication data set from Sweden (the Epidemiological Investigation of MS (EIMS) and the Genes and Environment in MS (GEMS) studies, 2005-2013; 5,133 MS cases and 4,718 controls), we constructed a weighted genetic risk score using 97 variants previously established to predict BMI. Results were adjusted for birth year, sex, education, smoking status, ancestry, and genetic predictors of MS. Estimates in KPNC and Swedish data sets suggested that higher genetically induced BMI predicted greater susceptibility to MS (odds ratio = 1.13, 95% confidence interval: 1.04, 1.22 for the KPNC sample; odds ratio = 1.09, 95% confidence interval: 1.03, 1.15 for the Swedish sample). Although the mechanism remains unclear, to our knowledge, these findings support a causal effect of increased BMI on susceptibility to MS for the first time, and they suggest a role for inflammatory pathways that characterize both obesity and the MS disease process.

Mendelian randomization shows a causal effect of low vitamin D on multiple sclerosis risk.

OBJECTIVE: We sought to estimate the causal effect of low serum 25(OH)D on multiple sclerosis (MS) susceptibility that is not confounded by environmental or lifestyle factors or subject to reverse causality. METHODS: We conducted mendelian randomization (MR) analyses using an instrumental variable (IV) comprising 3 single nucleotide polymorphisms found to be associated with serum 25(OH)D levels at genome-wide significance. We analyzed the effect of the IV on MS risk and both age at onset and disease severity in 2 separate populations using logistic regression models that controlled for sex, year of birth, smoking, education, genetic ancestry, body mass index at age 18-20 years or in 20s, a weighted genetic risk score for 110 known MS-associated variants, and the presence of one or more HLA-DRB1*15:01 alleles. RESULTS: Findings from MR analyses using the IV showed increasing levels of 25(OH)D are associated with a decreased risk of MS in both populations. In white, non-Hispanic members of Kaiser Permanente Northern California (1,056 MS cases and 9,015 controls), the odds ratio (OR) was 0.79 (p = 0.04, 95% confidence interval (CI): 0.64-0.99). In members of a Swedish population from the Epidemiological Investigation of Multiple Sclerosis and Genes and Environment in Multiple Sclerosis MS case-control studies (6,335 cases and 5,762 controls), the OR was 0.86 (p = 0.03, 95% CI: 0.76-0.98). A meta-analysis of the 2 populations gave a combined OR of 0.85 (p = 0.003, 95% CI: 0.76-0.94). No association was observed for age at onset or disease severity. CONCLUSIONS: These results provide strong evidence that low serum 25(OH)D concentration is a cause of MS, independent of established risk factors.

Genetic contribution to variation in DNA methylation at maternal smoking-sensitive loci in exposed neonates.

Epigenome-wide DNA methylation association studies have identified highly replicable genomic loci sensitive to maternal smoking during gestation. The role of inter-individual genetic variation in influencing DNA methylation, leading to the possibility of confounding or bias of such associations, has not been assessed. We investigated whether the DNA methylation levels at the top 10 CpG sites previously associated with exposure to maternal smoking during gestation were associated with individual genetic variation at the genome-wide level. Genome-wide association tests between DNA methylation at the top 10 candidate CpG and genome-wide SNPs were performed in 736 case and control participants of the California Childhood Leukemia Study. Three of the strongest maternal-smoking sensitive CpG sites in newborns were significantly associated with SNPs located proximal to each gene: cg18146737 in the GFI1 gene with rs141819830 (P = 8.2×10-44), cg05575921 in the AHRR gene with rs148405299 (P = 5.3×10-10), and cg12803068 in the MYO1G gene with rs61087368 (P = 1.3×10-18). For the GFI1 CpG cg18146737, the underlying genetic variation at rs141819830 confounded the association between maternal smoking and DNA methylation in our data (the regression coefficient changed from -0.02 [P = 0.139] to -0.03 [P = 0.015] after including the genotype). Our results suggest that further studies using DNA methylation at cg18146737, cg05575921, or cg12803068 that aim to assess exposure to maternal smoking during gestation should include genotype at the corresponding SNP. New methods are required for adequate and routine inclusion of genotypic influence on DNA methylation in epigenome-wide association studies to control for potential confounding.

Pathway Analysis of Genome-wide Association Study in Childhood Leukemia among Hispanics.

BACKGROUND: The incidence of acute lymphoblastic leukemia (ALL) is nearly 20% higher among Hispanics than non-Hispanic Whites. Previous studies have shown evidence for association between risk of ALL and variation within IKZF1, ARID5B, CEBPE, CDKN2A, GATA3, and BM1-PIP4K2A genes. However, variants identified only account for <10% of the genetic risk of ALL. METHODS: We applied pathway-based analyses to genome-wide association study (GWAS) data from the California Childhood Leukemia Study to determine whether different biologic pathways were overrepresented in childhood ALL and major ALL subtypes. Furthermore, we applied causal inference and data reduction methods to prioritize candidate genes within each identified overrepresented pathway, while accounting for correlation among SNPs. RESULTS: Pathway analysis results indicate that different ALL subtypes may involve distinct biologic mechanisms. Focal adhesion is a shared mechanism across the different disease subtypes. For ALL, the top five overrepresented Kyoto Encyclopedia of Genes and Genomes pathways include axon guidance, protein digestion and absorption, melanogenesis, leukocyte transendothelial migration, and focal adhesion (PFDR < 0.05). Notably, these pathways are connected to downstream MAPK or Wnt signaling pathways which have been linked to B-cell malignancies. Several candidate genes for ALL, such as COL6A6 and COL5A1, were identified through targeted maximum likelihood estimation. CONCLUSIONS: This is the first study to show distinct biologic pathways are overrepresented in different ALL subtypes using pathway-based approaches, and identified potential gene candidates using causal inference methods. IMPACT: The findings demonstrate that newly developed bioinformatics tools and causal inference methods can provide insights to furthering our understanding of the pathogenesis of leukemia. Cancer Epidemiol Biomarkers Prev; 25(5); 815-22. ©2016 AACR.

Feasibility study for remote assessment of cognitive function in multiple sclerosis.

BACKGROUND: Cognitive impairment is common in multiple sclerosis (MS), and affects employment and quality of life. Large studies are needed to identify risk factors for cognitive decline. Currently, a MS-validated remote assessment for cognitive function does not exist. Studies to determine feasibility of large remote cognitive function investigations in MS have not been published. OBJECTIVE: To determine whether MS patients would participate in remote cognitive studies. We utilized the Modified Telephone Interview for Cognitive Status (TICS-M), a previously validated phone assessment for cognitive function in healthy elderly populations to detect mild cognitive impairment. We identified factors that influenced participation rates. We investigated the relationship between MS risk factors and TICS-M score in cases, and score differences between cases and control individuals. METHODS: The TICS-M was administered to MS cases and controls. Linear and logistic regression models were utilized. RESULTS: 11.5% of eligible study participants did not participate in cognitive testing. MS cases, females and individuals with lower educational status were more likely to refuse (p<0.001). Cases who did complete testing did not differ in terms of perceived cognitive deficit compared to cases that did participate. More severe disease, smoking, and being male were associated with a lower TICS-M score among cases (p<0.001). The TICS-M score was significantly lower in cases compared to controls (p=0.007). CONCLUSIONS: Our results demonstrate convincingly that a remotely administered cognitive assessment is quite feasible for conducting large epidemiologic studies in MS, and lay the much needed foundation for future work that will utilize MS-validated cognitive measures.