Mendelian randomization (MR) permits causal inference between exposures and a disease. It can be compared with randomized controlled trials. Whereas in a randomized controlled trial the randomization occurs at entry into the trial, in MR the randomization occurs during gamete formation and conception. Several factors, including time since conception and sampling variation, are relevant to the interpretation of an MR test. Particularly important is consideration of the “missingness” of genotypes that can be originated by chance, genotyping errors, or clinical ascertainment. Testing for Hardy-Weinberg equilibrium (HWE) is a genetic approach that permits evaluation of missingness. In this paper, the authors demonstrate evidence of nonconformity with HWE in real data. They also perform simulations to characterize the sensitivity of HWE tests to missingness. Unresolved missingness could lead to a false rejection of causality in an MR investigation of trait-disease association. These results indicate that large-scale studies, very high quality genotyping data, and detailed knowledge of the life-course genetics of the alleles/genotypes studied will largely mitigate this risk. The authors also present a Web program (http://www.oege.org/software/hwe-mr-calc.shtml) for estimating possible missingness and an approach to evaluating missingness under different genetic models.

Objectives To determine whether more years spent in education is a causal risk factor for myopia, or whether myopia is a causal risk factor for more years in education. Design Bidirectional, two sample mendelian randomisation study. Setting Publically available genetic data from two consortiums applied to a large, independent population cohort. Genetic variants used as proxies for myopia and years of education were derived from two large genome wide association studies: 23andMe and Social Science Genetic Association Consortium (SSGAC), respectively. Participants 67 798 men and women from England, Scotland, and Wales in the UK Biobank cohort with available information for years of completed education and refractive error. Main outcome measures Mendelian randomisation analyses were performed in two directions: the first exposure was the genetic predisposition to myopia, measured with 44 genetic variants strongly associated with myopia in 23andMe, and the outcome was years in education; and the second exposure was the genetic predisposition to higher levels of education, measured with 69 genetic variants from SSGAC, and the outcome was refractive error. Results Conventional regression analyses of the observational data suggested that every additional year of education was associated with a more myopic refractive error of −0.18 dioptres/y (95% confidence interval −0.19 to −0.17; P<2e-16). Mendelian randomisation analyses suggested the true causal effect was even stronger: −0.27 dioptres/y (−0.37 to −0.17; P=4e-8). By contrast, there was little evidence to suggest myopia affected education (years in education per dioptre of refractive error −0.008 y/dioptre, 95% confidence interval −0.041 to 0.025, P=0.6). Thus, the cumulative effect of more years in education on refractive error means that a university graduate from the United Kingdom with 17 years of education would, on average, be at least −1 dioptre more myopic than someone who left school at age 16 (with 12 years of education). Myopia of this magnitude would be sufficient to necessitate the use of glasses for driving. Sensitivity analyses showed minimal evidence for genetic confounding that could have biased the causal effect estimates. Conclusions This study shows that exposure to more years in education contributes to the rising prevalence of myopia. Increasing the length of time spent in education may inadvertently increase the prevalence of myopia and potential future visual disability.

The frequency of a haplotype comprising one allele at each of two loci can be expressed as a cubic equation (the 'Hill equation'), the solution of which gives that frequency. Most haplotype and linkage disequilibrium analysis programs use iteration-based algorithms which substitute an estimate of haplotype frequency into the equation, producing a new estimate which is repeatedly fed back into the equation until the values converge to a maximum likelihood estimate (expectation-maximisation).

Multiple sclerosis (MS) is a neuro-inflammatory and neurodegenerative disease that is most prevalent in Northern Europe. Although it is known that inherited risk for MS is located within or in close proximity to immune-related genes, it is unknown when, where and how this genetic risk originated¹. Here, by using a large ancient genome dataset from the Mesolithic period to the Bronze Age², along with new Medieval and post-Medieval genomes, we show that the genetic risk for MS rose among pastoralists from the Pontic steppe and was brought into Europe by the Yamnaya-related migration approximately 5,000 years ago. We further show that these MS-associated immunogenetic variants underwent positive selection both within the steppe population and later in Europe, probably driven by pathogenic challenges coinciding with changes in diet, lifestyle and population density. This study highlights the critical importance of the Neolithic period and Bronze Age as determinants of modern immune responses and their subsequent effect on the risk of developing MS in a changing environment.

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.