Autism-related dietary preferences mediate autism-gut microbiome associations.

There is increasing interest in the potential contribution of the gut microbiome to autism spectrum disorder (ASD). However, previous studies have been underpowered and have not been designed to address potential confounding factors in a comprehensive way. We performed a large autism stool metagenomics study (n = 247) based on participants from the Australian Autism Biobank and the Queensland Twin Adolescent Brain project. We found negligible direct associations between ASD diagnosis and the gut microbiome. Instead, our data support a model whereby ASD-related restricted interests are associated with less-diverse diet, and in turn reduced microbial taxonomic diversity and looser stool consistency. In contrast to ASD diagnosis, our dataset was well powered to detect microbiome associations with traits such as age, dietary intake, and stool consistency. Overall, microbiome differences in ASD may reflect dietary preferences that relate to diagnostic features, and we caution against claims that the microbiome has a driving role in ASD.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis.

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma. Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness. To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness. The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry. We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10-8) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

Detection and replication of epistasis influencing transcription in humans.

Epistasis is the phenomenon whereby one polymorphism's effect on a trait depends on other polymorphisms present in the genome. The extent to which epistasis influences complex traits and contributes to their variation is a fundamental question in evolution and human genetics. Although often demonstrated in artificial gene manipulation studies in model organisms, and some examples have been reported in other species, few examples exist for epistasis among natural polymorphisms in human traits. Its absence from empirical findings may simply be due to low incidence in the genetic control of complex traits, but an alternative view is that it has previously been too technically challenging to detect owing to statistical and computational issues. Here we show, using advanced computation and a gene expression study design, that many instances of epistasis are found between common single nucleotide polymorphisms (SNPs). In a cohort of 846 individuals with 7,339 gene expression levels measured in peripheral blood, we found 501 significant pairwise interactions between common SNPs influencing the expression of 238 genes (P < 2.91 × 10(-16)). Replication of these interactions in two independent data sets showed both concordance of direction of epistatic effects (P = 5.56 × 10(-31)) and enrichment of interaction P values, with 30 being significant at a conservative threshold of P < 9.98 × 10(-5). Forty-four of the genetic interactions are located within 5 megabases of regions of known physical chromosome interactions (P = 1.8 × 10(-10)). Epistatic networks of three SNPs or more influence the expression levels of 129 genes, whereby one cis-acting SNP is modulated by several trans-acting SNPs. For example, MBNL1 is influenced by an additive effect at rs13069559, which itself is masked by trans-SNPs on 14 different chromosomes, with nearly identical genotype-phenotype maps for each cis-trans interaction. This study presents the first evidence, to our knowledge, for many instances of segregating common polymorphisms interacting to influence human traits.

Nick Martin and the Genetics of Depression: Sample Size, Sample Size, Sample Size.

Nick Martin is a pioneer in recognizing the need for large sample size to study the complex, heterogeneous and polygenic disorders of common mental disorders. In the predigital era, questionnaires were mailed to thousands of twin pairs around Australia. Always quick to adopt new technology, Nick's studies progressed to phone interviews and then online. Moreover, Nick was early to recognize the value of collecting DNA samples. As genotyping technologies improved over the years, these twin and family cohorts were used for linkage, candidate gene and genome-wide association studies. These cohorts have underpinned many analyses to disentangle the complex web of genetic and lifestyle factors associated with mental health. With characteristic foresight, Nick is chief investigator of our Australian Genetics of Depression Study, which has recruited 16,000 people with self-reported depression (plus DNA samples) over a time frame of a few months - analyses are currently ongoing. The mantra of sample size, sample size, sample size has guided Nick's research over the last 30 years and continues to do so.

Harmonised collection of data in youth mental health: Towards large datasets.

OBJECTIVE: The current international trend is to create large datasets with existing data and/or deposit newly collected data into repositories accessible to the scientific community. These practices lead to more efficient data sharing, better detection of small effects, modelling of confounders, establishment of sample generalizability and identification of differences between any given disorders. In Australia, to facilitate such data-sharing and collaborative opportunities, the Neurobiology in Youth Mental Health Partnership was created. This initiative brings together specialised researchers from around Australia to work towards a better understanding of the cross-diagnostic neurobiology of youth mental health and the translation of this knowledge into clinical practice. One of the mandates of the partnership was to develop a protocol for harmonised prospective collection of data across research centres in the field of youth mental health in order to create large datasets. METHODS: Four key research modalities were identified: clinical assessments, brain imaging, neurocognitive assessment and collection of blood samples. This paper presents the consensus set of assessments/data collection that has been selected by experts in each domain. CONCLUSION: The use of this core set of data will facilitate the pooling of psychopathological and neurobiological data into large datasets allowing researchers to tackle important questions requiring very large numbers. The aspiration of this transdiagnostic approach is a better understanding of the mechanisms underlying mental illnesses.

Evidence for mitochondrial genetic control of autosomal gene expression.

The mitochondrial and nuclear genomes coordinate and co-evolve in eukaryotes in order to adapt to environmental changes. Variation in the mitochondrial genome is capable of affecting expression of genes on the nuclear genome. Sex-specific mitochondrial genetic control of gene expression has been demonstrated in Drosophila melanogaster, where males were found to drive most of the total variation in gene expression. This has potential implications for male-related health and disease resulting from variation in mtDNA solely inherited from the mother. We used a family-based study comprised of 47,323 gene expression probes and 78 mitochondrial SNPs (mtSNPs) from n = 846 individuals to examine the extent of mitochondrial genetic control of gene expression in humans. This identified 15 significant probe-mtSNP associations (P<10-8) corresponding to 5 unique genes on the mitochondrial and nuclear genomes, with three of these genes corresponding to mitochondrial genetic control of gene expression in the nuclear genome. The associated mtSNPs for three genes (one cis and two trans associations) were replicated (P < 0.05) in an independent dataset of n = 452 unrelated individuals. There was no evidence for sexual dimorphic gene expression in any of these five probes. Sex-specific effects were examined by applying our analysis to males and females separately and testing for differences in effect size. The MEST gene was identified as having the most significantly different effect sizes across the sexes (P≈10-7). MEST was similarly expressed in males and females with the G allele; however, males with the C allele are highly expressed for MEST, while females show no expression of the gene. This study provides evidence for the mitochondrial genetic control of expression of several genes in humans, with little evidence found for sex-specific effects.

Improving Phenotypic Prediction by Combining Genetic and Epigenetic Associations.

We tested whether DNA-methylation profiles account for inter-individual variation in body mass index (BMI) and height and whether they predict these phenotypes over and above genetic factors. Genetic predictors were derived from published summary results from the largest genome-wide association studies on BMI (n ∼ 350,000) and height (n ∼ 250,000) to date. We derived methylation predictors by estimating probe-trait effects in discovery samples and tested them in external samples. Methylation profiles associated with BMI in older individuals from the Lothian Birth Cohorts (LBCs, n = 1,366) explained 4.9% of the variation in BMI in Dutch adults from the LifeLines DEEP study (n = 750) but did not account for any BMI variation in adolescents from the Brisbane Systems Genetic Study (BSGS, n = 403). Methylation profiles based on the Dutch sample explained 4.9% and 3.6% of the variation in BMI in the LBCs and BSGS, respectively. Methylation profiles predicted BMI independently of genetic profiles in an additive manner: 7%, 8%, and 14% of variance of BMI in the LBCs were explained by the methylation predictor, the genetic predictor, and a model containing both, respectively. The corresponding percentages for LifeLines DEEP were 5%, 9%, and 13%, respectively, suggesting that the methylation profiles represent environmental effects. The differential effects of the BMI methylation profiles by age support previous observations of age modulation of genetic contributions. In contrast, methylation profiles accounted for almost no variation in height, consistent with a mainly genetic contribution to inter-individual variation. The BMI results suggest that combining genetic and epigenetic information might have greater utility for complex-trait prediction.

Hypermetabolism in ALS is associated with greater functional decline and shorter survival.

OBJECTIVE: To determine the prevalence of hypermetabolism, relative to body composition, in amyotrophic lateral sclerosis (ALS) and its relationship with clinical features of disease and survival. METHODS: Fifty-eight patients with clinically definite or probable ALS as defined by El Escorial criteria, and 58 age and sex-matched control participants underwent assessment of energy expenditure. Our primary outcome was the prevalence of hypermetabolism in cases and controls. Longitudinal changes in clinical parameters between hypermetabolic and normometabolic patients with ALS were determined for up to 12 months following metabolic assessment. Survival was monitored over a 30-month period following metabolic assessment. RESULTS: Hypermetabolism was more prevalent in patients with ALS than controls (41% vs 12%, adjusted OR=5.4; p<0.01). Change in body weight, body mass index and fat mass (%) was similar between normometabolic and hypermetabolic patients with ALS. Mean lower motor neuron score (SD) was greater in hypermetabolic patients when compared with normometabolic patients (4 (0.3) vs 3 (0.7); p=0.04). In the 12 months following metabolic assessment, there was a greater change in Revised ALS Functional Rating Scale score in hypermetabolic patients when compared with normometabolic patients (-0.68 points/month vs -0.39 points/month; p=0.01). Hypermetabolism was inversely associated with survival. Overall, hypermetabolism increased the risk of death during follow-up to 220% (HR 3.2, 95% CI 1.1 to 9.4, p=0.03). CONCLUSIONS AND RELEVANCE: Hypermetabolic patients with ALS have a greater level of lower motor neuron involvement, faster rate of functional decline and shorter survival. The metabolic index could be important for informing prognosis in ALS.

Study protocol for the Australian autism biobank: an international resource to advance autism discovery research.

BACKGROUND: The phenotypic and genetic heterogeneity of autism spectrum disorder (ASD) presents considerable challenges in understanding etiological pathways, selecting effective therapies, providing genetic counselling, and predicting clinical outcomes. With advances in genetic and biological research alongside rapid-pace technological innovations, there is an increasing imperative to access large, representative, and diverse cohorts to advance knowledge of ASD. To date, there has not been any single collective effort towards a similar resource in Australia, which has its own unique ethnic and cultural diversity. The Australian Autism Biobank was initiated by the Cooperative Research Centre for Living with Autism (Autism CRC) to establish a large-scale repository of biological samples and detailed clinical information about children diagnosed with ASD to facilitate future discovery research. METHODS: The primary group of participants were children with a confirmed diagnosis of ASD, aged between 2 and 17 years, recruited through four sites in Australia. No exclusion criteria regarding language level, cognitive ability, or comorbid conditions were applied to ensure a representative cohort was recruited. Both biological parents and siblings were invited to participate, along with children without a diagnosis of ASD, and children who had been queried for an ASD diagnosis but did not meet diagnostic criteria. All children completed cognitive assessments, with probands and parents completing additional assessments measuring ASD symptomatology. Parents completed questionnaires about their child's medical history and early development. Physical measurements and biological samples (blood, stool, urine, and hair) were collected from children, and physical measurements and blood samples were collected from parents. Samples were sent to a central processing site and placed into long-term storage. DISCUSSION: The establishment of this biobank is a valuable international resource incorporating detailed clinical and biological information that will help accelerate the pace of ASD discovery research. Recruitment into this study has also supported the feasibility of large-scale biological sample collection in children diagnosed with ASD with comprehensive phenotyping across a wide range of ages, intellectual abilities, and levels of adaptive functioning. This biological and clinical resource will be open to data access requests from national and international researchers to support future discovery research that will benefit the autistic community.

Fine-mapping of the HNF1B multicancer locus identifies candidate variants that mediate endometrial cancer risk.

Common variants in the hepatocyte nuclear factor 1 homeobox B (HNF1B) gene are associated with the risk of Type II diabetes and multiple cancers. Evidence to date indicates that cancer risk may be mediated via genetic or epigenetic effects on HNF1B gene expression. We previously found single-nucleotide polymorphisms (SNPs) at the HNF1B locus to be associated with endometrial cancer, and now report extensive fine-mapping and in silico and laboratory analyses of this locus. Analysis of 1184 genotyped and imputed SNPs in 6608 Caucasian cases and 37 925 controls, and 895 Asian cases and 1968 controls, revealed the best signal of association for SNP rs11263763 (P = 8.4 × 10(-14), odds ratio = 0.86, 95% confidence interval = 0.82-0.89), located within HNF1B intron 1. Haplotype analysis and conditional analyses provide no evidence of further independent endometrial cancer risk variants at this locus. SNP rs11263763 genotype was associated with HNF1B mRNA expression but not with HNF1B methylation in endometrial tumor samples from The Cancer Genome Atlas. Genetic analyses prioritized rs11263763 and four other SNPs in high-to-moderate linkage disequilibrium as the most likely causal SNPs. Three of these SNPs map to the extended HNF1B promoter based on chromatin marks extending from the minimal promoter region. Reporter assays demonstrated that this extended region reduces activity in combination with the minimal HNF1B promoter, and that the minor alleles of rs11263763 or rs8064454 are associated with decreased HNF1B promoter activity. Our findings provide evidence for a single signal associated with endometrial cancer risk at the HNF1B locus, and that risk is likely mediated via altered HNF1B gene expression.

Genetic and environmental exposures constrain epigenetic drift over the human life course.

Epigenetic mechanisms such as DNA methylation (DNAm) are essential for regulation of gene expression. DNAm is dynamic, influenced by both environmental and genetic factors. Epigenetic drift is the divergence of the epigenome as a function of age due to stochastic changes in methylation. Here we show that epigenetic drift may be constrained at many CpGs across the human genome by DNA sequence variation and by lifetime environmental exposures. We estimate repeatability of DNAm at 234,811 autosomal CpGs in whole blood using longitudinal data (2-3 repeated measurements) on 478 older people from two Scottish birth cohorts--the Lothian Birth Cohorts of 1921 and 1936. Median age was 79 yr and 70 yr, and the follow-up period was ∼10 yr and ∼6 yr, respectively. We compare this to methylation heritability estimated in the Brisbane Systems Genomics Study, a cross-sectional study of 117 families (offspring median age 13 yr; parent median age 46 yr). CpG repeatability in older people was highly correlated (0.68) with heritability estimated in younger people. Highly heritable sites had strong underlying cis-genetic effects. Thirty-seven and 1687 autosomal CpGs were associated with smoking and sex, respectively. Both sets were strongly enriched for high repeatability. Sex-associated CpGs were also strongly enriched for high heritability. Our results show that a large number of CpGs across the genome, as a result of environmental and/or genetic constraints, have stable DNAm variation over the human lifetime. Moreover, at a number of CpGs, most variation in the population is due to genetic factors, despite some sites being highly modifiable by the environment.

The genetic regulation of transcription in human endometrial tissue.

Study question: Do genetic effects regulate gene expression in human endometrium? Summary answer: This study demonstrated strong genetic effects on endometrial gene expression and some evidence for genetic regulation of gene expression in a menstrual cycle stage-specific manner. What is known already: Genetic effects on expression levels for many genes are tissue specific. Endometrial gene expression varies across menstrual cycle stages and between individuals, but there are limited data on genetic control of expression in endometrium. Study design, size, duration: We analysed genome-wide genotype and gene expression data to map cis expression quantitative trait loci (eQTL) in endometrium. Participants/materials, setting, methods: We recruited 123 women of European ancestry. DNA samples from blood were genotyped on Illumina HumanCoreExome chips. Total RNA was extracted from endometrial tissues. Whole-transcriptome profiles were characterized using Illumina Human HT-12 v4.0 Expression Beadchips. We performed eQTL mapping with ~8 000 000 genotyped and imputed single nucleotide polymorphisms (SNPs) and 12 329 genes. Main results and the role of chance: We identified a total of 18 595 cis SNP-probe associations at a study-wide level of significance (P < 1 × 10-7), which correspond to independent eQTLs for 198 unique genes. The eQTLs with the largest effect in endometrial tissue were rs4902335 for CHURC1 (P = 1.05 × 10-32) and rs147253019 for ZP3 (P = 8.22 × 10-30). We further performed a context-specific eQTL analysis to investigate if genetic effects on gene expression regulation act in a menstrual cycle-specific manner. Interestingly, five cis-eQTLs were identified with a significant stage-by-genotype interaction. The strongest stage interaction was the eQTL for C10ORF33 (PYROXD2) with SNP rs2296438 (P = 2.0 × 10-4), where we observe a 2-fold difference in the average expression levels of heterozygous samples depending on the stage of the menstrual cycle. Large scale data: The summary eQTL results are publicly available to browse or download. Limitations, reasons for caution: A limitation of the present study was the relatively modest sample size. It was not powered to identify trans-eQTLs and larger sample sizes will also be needed to provide better power to detect cis-eQTLs and cycle stage-specific effects, given the substantial changes in expression across the menstrual cycle for many genes. Wider implications of the findings: Identification of endometrial eQTLs provides a platform for better understanding genetic effects on endometriosis risk and other endometrial-related pathologies. Study funding/competing interest(s): Funding for this work was provided by NHMRC Project Grants GNT1026033, GNT1049472, GNT1046880, GNT1050208, GNT1105321 and APP1083405. There are no competing interests.

Genomewide Association Study of Alcohol Dependence Identifies Risk Loci Altering Ethanol-Response Behaviors in Model Organisms.

BACKGROUND: Alcohol dependence (AD) shows evidence for genetic liability, but genes influencing risk remain largely unidentified. METHODS: We conducted a genomewide association study in 706 related AD cases and 1,748 unscreened population controls from Ireland. We sought replication in 15,496 samples of European descent. We used model organisms (MOs) to assess the role of orthologous genes in ethanol (EtOH)-response behaviors. We tested 1 primate-specific gene for expression differences in case/control postmortem brain tissue. RESULTS: We detected significant association in COL6A3 and suggestive association in 2 previously implicated loci, KLF12 and RYR3. None of these signals are significant in replication. A suggestive signal in the long noncoding RNA LOC339975 is significant in case:control meta-analysis, but not in a population sample. Knockdown of a COL6A3 ortholog in Caenorhabditis elegans reduced EtOH sensitivity. Col6a3 expression correlated with handling-induced convulsions in mice. Loss of function of the KLF12 ortholog in C. elegans impaired development of acute functional tolerance (AFT). Klf12 expression correlated with locomotor activation following EtOH injection in mice. Loss of function of the RYR3 ortholog reduced EtOH sensitivity in C. elegans and rapid tolerance in Drosophila. The ryanodine receptor antagonist dantrolene reduced motivation to self-administer EtOH in rats. Expression of LOC339975 does not differ between cases and controls but is reduced in carriers of the associated rs11726136 allele in nucleus accumbens (NAc). CONCLUSIONS: We detect association between AD and COL6A3, KLF12, RYR3, and LOC339975. Despite nonreplication of COL6A3, KLF12, and RYR3 signals, orthologs of these genes influence behavioral response to EtOH in MOs, suggesting potential involvement in human EtOH response and AD liability. The associated LOC339975 allele may influence gene expression in human NAc. Although the functions of long noncoding RNAs are poorly understood, there is mounting evidence implicating these genes in multiple brain functions and disorders.

Shared genetic control of expression and methylation in peripheral blood.

BACKGROUND: Expression QTLs and epigenetic marks are often employed to provide an insight into the possible biological mechanisms behind GWAS hits. A substantial proportion of the variation in gene expression and DNA methylation is known to be under genetic control. We address the proportion of genetic control that is shared between these two genomic features. RESULTS: An exhaustive search for pairwise phenotypic correlations between gene expression and DNA methylation in samples from human blood (n = 610) was performed. Of the 5 × 10(9) possible pairwise tests, 0.36 % passed Bonferroni corrected p-value cutoff of 9.9 × 10(-12). We determined that the correlation structure between probe pairs was largely due to blood cell type specificity of the expression and methylation probes. Upon adjustment of the expression and methylation values for observed blood cellular composition (n = 422), the number of probe pairs which survived Bonferroni correction reduced by more than 5400 fold. Of the 614 correlated probe pairs located on the same chromosome, 75 % share at least one methylation and expression QTL at nominal 10(-5) p-value cutoff. Those probe pairs are located within 1Mbp window from each other and have a mean of absolute value of genetic correlation equal to 0.69, further demonstrating the high degree of shared genetic control. CONCLUSIONS: Overall, this study demonstrates notable genetic covariance between DNA methylation and gene expression and reaffirms the importance of correcting for cell-counts in studies on non-homogeneous tissues.

Inference of the genetic architecture underlying BMI and height with the use of 20,240 sibling pairs.

Evidence that complex traits are highly polygenic has been presented by population-based genome-wide association studies (GWASs) through the identification of many significant variants, as well as by family-based de novo sequencing studies indicating that several traits have a large mutational target size. Here, using a third study design, we show results consistent with extreme polygenicity for body mass index (BMI) and height. On a sample of 20,240 siblings (from 9,570 nuclear families), we used a within-family method to obtain narrow-sense heritability estimates of 0.42 (SE = 0.17, p = 0.01) and 0.69 (SE = 0.14, p = 6 × 10(-)(7)) for BMI and height, respectively, after adjusting for covariates. The genomic inflation factors from locus-specific linkage analysis were 1.69 (SE = 0.21, p = 0.04) for BMI and 2.18 (SE = 0.21, p = 2 × 10(-10)) for height. This inflation is free of confounding and congruent with polygenicity, consistent with observations of ever-increasing genomic-inflation factors from GWASs with large sample sizes, implying that those signals are due to true genetic signals across the genome rather than population stratification. We also demonstrate that the distribution of the observed test statistics is consistent with both rare and common variants underlying a polygenic architecture and that previous reports of linkage signals in complex traits are probably a consequence of polygenic architecture rather than the segregation of variants with large effects. The convergent empirical evidence from GWASs, de novo studies, and within-family segregation implies that family-based sequencing studies for complex traits require very large sample sizes because the effects of causal variants are small on average.

Meta-analysis of Genome-Wide Association Studies for Extraversion: Findings from the Genetics of Personality Consortium.

Extraversion is a relatively stable and heritable personality trait associated with numerous psychosocial, lifestyle and health outcomes. Despite its substantial heritability, no genetic variants have been detected in previous genome-wide association (GWA) studies, which may be due to relatively small sample sizes of those studies. Here, we report on a large meta-analysis of GWA studies for extraversion in 63,030 subjects in 29 cohorts. Extraversion item data from multiple personality inventories were harmonized across inventories and cohorts. No genome-wide significant associations were found at the single nucleotide polymorphism (SNP) level but there was one significant hit at the gene level for a long non-coding RNA site (LOC101928162). Genome-wide complex trait analysis in two large cohorts showed that the additive variance explained by common SNPs was not significantly different from zero, but polygenic risk scores, weighted using linkage information, significantly predicted extraversion scores in an independent cohort. These results show that extraversion is a highly polygenic personality trait, with an architecture possibly different from other complex human traits, including other personality traits. Future studies are required to further determine which genetic variants, by what modes of gene action, constitute the heritable nature of extraversion.

Genetic loci for Epstein-Barr virus nuclear antigen-1 are associated with risk of multiple sclerosis.

BACKGROUND: Infection with the Epstein-Barr virus (EBV) is associated with an increased risk of multiple sclerosis (MS). OBJECTIVE: We sought genetic loci influencing EBV nuclear antigen-1 (EBNA-1) IgG titers and hypothesized that they may play a role in MS risk. METHODS: We performed a genome-wide association study (GWAS) of anti-EBNA-1 IgG titers in 3599 individuals from an unselected twin family cohort, followed by a meta-analysis with data from an independent EBNA-1 GWAS. We then examined the shared polygenic risk between the EBNA-1 GWAS (effective sample size (Neff) = 5555) and a large MS GWAS (Neff = 15,231). RESULTS: We identified one locus of strong association within the human leukocyte antigen (HLA) region, of which the most significantly associated genotyped single nucleotide polymorphism (SNP) was rs2516049 (p = 4.11 × 10-9). A meta-analysis including data from another EBNA-1 GWAS in a cohort of Mexican-American families confirmed that rs2516049 remained the most significantly associated SNP (p = 3.32 × 10-20). By examining the shared polygenic risk, we show that the genetic risk for elevated anti-EBNA-1 titers is positively correlated with the development of MS, and that elevated EBNA-1 titers are not an epiphenomena secondary to MS. In the joint meta-analysis of EBNA-1 titers and MS, loci at 1p22.1, 3p24.1, 3q13.33, and 10p15.1 reached genome-wide significance (p < 5 × 10-8). CONCLUSIONS: Our results suggest that apart from the confirmed HLA region, the association of anti-EBNA-1 IgG titer with MS risk is also mediated through non-HLA genes, and that studies aimed at identifying genetic loci influencing EBNA immune response provides a novel opportunity to identify new and characterize existing genetic risk factors for MS.