The UK Biobank resource with deep phenotyping and genomic data

The molecular genetics of hand preference revisited

Hand preference is a prominent behavioural trait linked to human brain asymmetry. A handful of genetic variants have been reported to associate with hand preference or quantitative measures related to it. Most of these reports were on the basis of limited sample sizes, by current standards for genetic analysis of complex traits. Here we performed a genome-wide association analysis of hand preference in the large, population-based UK Biobank cohort (N = 331,037). We used gene-set enrichment analysis to investigate whether genes involved in visceral asymmetry are particularly relevant to hand preference, following one previous report. We found no evidence supporting any of the previously suggested variants or genes, nor that genes involved in visceral laterality have a role in hand preference. It remains possible that some of the previously reported genes or pathways are relevant to hand preference as assessed in other ways, or else are relevant within specific disorder populations. However, some or all of the earlier findings are likely to be false positives, and none of them appear relevant to hand preference as defined categorically in the general population. Our analysis did produce a small number of novel, significant associations, including one implicating the microtubule-associated gene MAP2 in handedness.

Network Medicine in the Age of Biomedical Big Data

Network medicine is an emerging area of research dealing with molecular and genetic interactions, network biomarkers of disease, and therapeutic target discovery. Large-scale biomedical data generation offers a unique opportunity to assess the effect and impact of cellular heterogeneity and environmental perturbations on the observed phenotype. Marrying the two, network medicine with biomedical data provides a framework to build meaningful models and extract impactful results at a network level. In this review, we survey existing network types and biomedical data sources. More importantly, we delve into ways in which the network medicine approach, aided by phenotype-specific biomedical data, can be gainfully applied. We provide three paradigms, mainly dealing with three major biological network archetypes: protein-protein interaction, expression-based, and gene regulatory networks. For each of these paradigms, we discuss a broad overview of philosophies under which various network methods work. We also provide a few examples in each paradigm as a test case of its successful application. Finally, we delineate several opportunities and challenges in the field of network medicine. We hope this review provides a lexicon for researchers from biological sciences and network theory to come on the same page to work on research areas that require interdisciplinary expertise. Taken together, the understanding gained from combining biomedical data with networks can be useful for characterizing disease etiologies and identifying therapeutic targets, which, in turn, will lead to better preventive medicine with translational impact on personalized healthcare.

Stroke genetics: discovery, biology, and clinical applications

Stroke, a leading cause of long-term disability and death worldwide, has a heritable component. Recent gene discovery efforts have expanded the number of known single-gene disorders associated with stroke and have linked common variants at approximately 35 genetic loci to stroke risk. These discoveries have highlighted novel mechanisms and pathways implicated in stroke related to large artery atherosclerosis, cardioembolism, and small vessel disease, and defined shared genetic influences with related vascular traits. Genetics has also successfully established causal relationships with risk factors and holds promise for prioritising targets for exploration in clinical trials. Genome-wide polygenic scores enable the identification of high-risk individuals before the emergence of vascular risk factors. Challenges ahead include a better understanding of rare variants and ancestral differences for integration of genetics into precision medicine, integration with other omics data, uncovering the genetic factors that govern stroke recurrence and stroke outcome, and the conversion of genetic discoveries to novel therapies.

Genetic studies of accelerometer-based sleep measures yield new insights into human sleep behaviour

Sleep is an essential human function but its regulation is poorly understood. Using accelerometer data from 85,670 UK Biobank participants, we perform a genome-wide association study of 8 derived sleep traits representing sleep quality, quantity and timing, and validate our findings in 5,819 individuals. We identify 47 genetic associations at P < 5 × 10-8, of which 20 reach a stricter threshold of P < 8 × 10-10. These include 26 novel associations with measures of sleep quality and 10 with nocturnal sleep duration. The majority of identified variants associate with a single sleep trait, except for variants previously associated with restless legs syndrome. For sleep duration we identify a missense variant (p.Tyr727Cys) in PDE11A as the likely causal variant. As a group, sleep quality loci are enriched for serotonin processing genes. Although accelerometer-derived measures of sleep are imperfect and may be affected by restless legs syndrome, these findings provide new biological insights into sleep compared to previous efforts based on self-report sleep measures.

Genetic Architectures of Childhood- and Adult-Onset Asthma Are Partly Distinct

The extent to which genetic risk factors are shared between childhood-onset (COA) and adult-onset (AOA) asthma has not been estimated. On the basis of data from the UK Biobank study (n = 447,628), we found that the variance in disease liability explained by common variants is higher for COA (onset at ages between 0 and 19 years; h²_g = 25.6%) than for AOA (onset at ages between 20 and 60 years; h²_g = 10.6%). The genetic correlation (r_g) between COA and AOA was 0.67. Variation in age of onset among COA-affected individuals had a low heritability (h²_g = 5%), which we confirmed in independent studies and also among AOA-affected individuals. To identify subtype-specific genetic associations, we performed a genome-wide association study (GWAS) in the UK Biobank for COA (13,962 affected individuals) and a separate GWAS for AOA (26,582 affected individuals) by using a common set of 300,671 controls for both studies. We identified 123 independent associations for COA and 56 for AOA (37 overlapped); of these, 98 and 34, respectively, were reproducible in an independent study (n = 262,767). Collectively, 28 associations were not previously reported. For 96 COA-associated variants, including five variants that represent COA-specific risk factors, the risk allele was more common in COA- than in AOA-affected individuals. Conversely, we identified three variants that are stronger risk factors for AOA. Variants associated with obesity and smoking had a stronger contribution to the risk of AOA than to the risk of COA. Lastly, we identified 109 likely target genes of the associated variants, primarily on the basis of correlated expression quantitative trait loci (up to n = 31,684). GWAS informed by age of onset can identify subtype-specific risk variants, which can help us understand differences in pathophysiology between COA and AOA and so can be informative for drug development.

Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases

Microbiome-wide association studies on large population cohorts have highlighted associations between the gut microbiome and complex traits, including type 2 diabetes (T2D) and obesity1. However, the causal relationships remain largely unresolved. We leveraged information from 952 normoglycemic individuals for whom genome-wide genotyping, gut metagenomic sequence and fecal short-chain fatty acid (SCFA) levels were available2, then combined this information with genome-wide-association summary statistics for 17 metabolic and anthropometric traits. Using bidirectional Mendelian randomization (MR) analyses to assess causality3, we found that the host-genetic-driven increase in gut production of the SCFA butyrate was associated with improved insulin response after an oral glucose-tolerance test (P = 9.8 × 10-5), whereas abnormalities in the production or absorption of another SCFA, propionate, were causally related to an increased risk of T2D (P = 0.004). These data provide evidence of a causal effect of the gut microbiome on metabolic traits and support the use of MR as a means to elucidate causal relationships from microbiome-wide association findings.

Sequence variants associating with urinary biomarkers

Urine dipstick tests are widely used in routine medical care to diagnose kidney and urinary tract and metabolic diseases. Several environmental factors are known to affect the test results, whereas the effects of genetic diversity are largely unknown. We tested 32.5 million sequence variants for association with urinary biomarkers in a set of 150 274 Icelanders with urine dipstick measurements. We detected 20 association signals, of which 14 are novel, associating with at least one of five clinical entities defined by the urine dipstick: glucosuria, ketonuria, proteinuria, hematuria and urine pH. These include three independent glucosuria variants at SLC5A2, the gene encoding the sodium-dependent glucose transporter (SGLT2), a protein targeted pharmacologically to increase urinary glucose excretion in the treatment of diabetes. Two variants associating with proteinuria are in LRP2 and CUBN, encoding the co-transporters megalin and cubilin, respectively, that mediate proximal tubule protein uptake. One of the hematuria-associated variants is a rare, previously unreported 2.5 kb exonic deletion in COL4A3. Of the four signals associated with urine pH, we note that the pH-increasing alleles of two variants (POU2AF1, WDR72) associate significantly with increased risk of kidney stones. Our results reveal that genetic factors affect variability in urinary biomarkers, in both a disease dependent and independent context.

Off the street phasing (OTSP): no hassle haplotype phasing for molecular PGD applications

PURPOSE

Pre-implantation genetic diagnosis (PGD) for molecular disorders requires the construction of parental haplotypes. Classically, haplotype resolution ("phasing") is obtained by genotyping multiple polymorphic markers in both parents and at least one additional relative. However, this process is time-consuming, and immediate family members are not always available. The recent availability of massive genomic data for many populations promises to eliminate the needs for developing family-specific assays and for recruiting additional family members. In this study, we aimed to validate population-assisted haplotype phasing for PGD.

METHODS

Targeted sequencing of CFTR gene variants and ~ 1700 flanking polymorphic SNPs (± 2 Mb) was performed on 54 individuals from 12 PGD families of (a) Full Ashkenazi (FA; n = 16), (b) mixed Ashkenazi (MA; n = 23 individuals with at least one Ashkenazi and one non-Ashkenazi grandparents), or (c) non-Ashkenazi (NA; n = 15) descent. Heterozygous genotype calls in each individual were phased using various whole genome reference panels and appropriate computational models. All computationally derived haplotype predictions were benchmarked against trio-based phasing.

RESULTS

Using the Ashkenazi reference panel, phasing of FA was highly accurate (99.4% ± 0.2% accuracy); phasing of MA was less accurate (95.4% ± 4.5% accuracy); and phasing of NA was predictably low (83.4% ± 6.6% accuracy). Strikingly, for founder mutation carriers, our haplotyping approach facilitated near perfect phasing accuracy (99.9% ± 0.1% and 98.2% ± 2.8% accuracy for W1282X and delF508 carriers, respectively).

CONCLUSIONS

Our results demonstrate the feasibility of replacing classical haplotype phasing with population-based phasing with uncompromised accuracy.

Parent of origin genetic effects on methylation in humans are common and influence complex trait variation

Parent-of-origin effects (POE) exist when there is differential expression of alleles inherited from the two parents. A genome-wide scan for POE on DNA methylation at 639,238 CpGs in 5,101 individuals identifies 733 independent methylation CpGs potentially influenced by POE at a false discovery rate ≤ 0.05 of which 331 had not previously been identified. Cis and trans methylation quantitative trait loci (mQTL) regulate methylation variation through POE at 54% (399/733) of the identified POE-influenced CpGs. The combined results provide strong evidence for previously unidentified POE-influenced CpGs at 171 independent loci. Methylation variation at 14 of the POE-influenced CpGs is associated with multiple metabolic traits. A phenome-wide association analysis using the POE mQTL SNPs identifies a previously unidentified imprinted locus associated with waist circumference. These results provide a high resolution population-level map for POE on DNA methylation sites, their local and distant regulators and potential consequences for complex traits.

Sequence variation at ANAPC1 accounts for 24% of the variability in corneal endothelial cell density

The corneal endothelium is vital for transparency and proper hydration of the cornea. Here, we conduct a genome-wide association study of corneal endothelial cell density (cells/mm2), coefficient of cell size variation (CV), percentage of hexagonal cells (HEX) and central corneal thickness (CCT) in 6,125 Icelanders and find associations at 10 loci, including 7 novel. We assess the effects of these variants on various ocular biomechanics such as corneal hysteresis (CH), as well as eye diseases such as glaucoma and corneal dystrophies. Most notably, an intergenic variant close to ANAPC1 (rs78658973[A], frequency = 28.3%) strongly associates with decreased cell density and accounts for 24% of the population variance in cell density (β = -0.77 SD, P = 1.8 × 10-314) and associates with increased CH (β = 0.19 SD, P = 2.6 × 10-19) without affecting risk of corneal diseases and glaucoma. Our findings indicate that despite correlations between cell density and eye diseases, low cell density does not increase the risk of disease.

Genomic Research Through an Indigenous Lens: Understanding the Expectations

Indigenous scholars are leading initiatives to improve access to genetic and genomic research and health care based on their unique cultural contexts and within sovereign-based governance models created and accepted by their peoples. In the past, Indigenous peoples' engagement with genomicresearch was hampered by a lack of standardized guidelines and institutional partnerships, resulting in group harms. This article provides a comparative analysis of research guidelines from Canada, New Zealand, Australia, and the United States that pertain to Indigenous peoples. The goals of the analysis are to identify areas that need attention, support Indigenous-led governance, and promote the development of a model research policy framework for genomic research and health care that has international relevance for Indigenous peoples.

Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates

Sleep is an essential state of decreased activity and alertness but molecular factors regulating sleep duration remain unknown. Through genome-wide association analysis in 446,118 adults of European ancestry from the UK Biobank, we identify 78 loci for self-reported habitual sleep duration (p < 5 × 10-8; 43 loci at p < 6 × 10-9). Replication is observed for PAX8, VRK2, and FBXL12/UBL5/PIN1 loci in the CHARGE study (n = 47,180; p < 6.3 × 10-4), and 55 signals show sign-concordant effects. The 78 loci further associate with accelerometer-derived sleep duration, daytime inactivity, sleep efficiency and number of sleep bouts in secondary analysis (n = 85,499). Loci are enriched for pathways including striatum and subpallium development, mechanosensory response, dopamine binding, synaptic neurotransmission and plasticity, among others. Genetic correlation indicates shared links with anthropometric, cognitive, metabolic, and psychiatric traits and two-sample Mendelian randomization highlights a bidirectional causal link with schizophrenia. This work provides insights into the genetic basis for inter-individual variation in sleep duration implicating multiple biological pathways.

Dysregulation of a long noncoding RNA reduces leptin leading to a leptin-responsive form of obesity

Quantitative changes in leptin concentration lead to alterations in food intake and body weight, but the regulatory mechanisms that control leptin gene expression are poorly understood. Here we report that fat-specific and quantitative leptin expression is controlled by redundant cis elements and trans factors interacting with the proximal promoter together with a long noncoding RNA (lncOb). Diet-induced obese mice lacking lncOb show increased fat mass with reduced plasma leptin levels and lose weight after leptin treatment, whereas control mice do not. Consistent with this finding, large-scale genetic studies of humans reveal a significant association of single-nucleotide polymorphisms (SNPs) in the region of human lncOb with lower plasma leptin levels and obesity. These results show that reduced leptin gene expression can lead to a hypoleptinemic, leptin-responsive form of obesity and provide a framework for elucidating the pathogenic mechanism in the subset of obese patients with low endogenous leptin levels.

Genome wide analysis for mouth ulcers identifies associations at immune regulatory loci

Mouth ulcers are the most common ulcerative condition and encompass several clinical diagnoses, including recurrent aphthous stomatitis (RAS). Despite previous evidence for heritability, it is not clear which specific genetic loci are implicated in RAS. In this genome-wide association study (n = 461,106) heritability is estimated at 8.2% (95% CI: 6.4%, 9.9%). This study finds 97 variants which alter the odds of developing non-specific mouth ulcers and replicate these in an independent cohort (n = 355,744) (lead variant after meta-analysis: rs76830965, near IL12A, OR 0.72 (95% CI: 0.71, 0.73); P = 4.4e−483). Additional effect estimates from three independent cohorts with more specific phenotyping and specific study characteristics support many of these findings. In silico functional analyses provide evidence for a role of T cell regulation in the aetiology of mouth ulcers. These results provide novel insight into the pathogenesis of a common, important condition.

Oral ulcerations are sores of the mucous membrane of the mouth and highly prevalent in the population. Here, in a genome-wide association study, the authors identify 97 loci associated with mouth ulcers highlighting genes involved in T cell-mediated immunity and T_H1 responses.

An atlas of polygenic risk score associations to highlight putative causal relationships across the human phenome

The age of large-scale genome-wide association studies (GWAS) has provided us with an unprecedented opportunity to evaluate the genetic liability of complex disease using polygenic risk scores (PRS). In this study, we have analysed 162 PRS (p<5×10-05) derived from GWAS and 551 heritable traits from the UK Biobank study (N = 334,398). Findings can be investigated using a web application (http:‌//‌mrcieu.‌mrsoftware.org/‌PRS‌_atlas/), which we envisage will help uncover both known and novel mechanisms which contribute towards disease susceptibility. To demonstrate this, we have investigated the results from a phenome-wide evaluation of schizophrenia genetic liability. Amongst findings were inverse associations with measures of cognitive function which extensive follow-up analyses using Mendelian randomization (MR) provided evidence of a causal relationship. We have also investigated the effect of multiple risk factors on disease using mediation and multivariable MR frameworks. Our atlas provides a resource for future endeavours seeking to unravel the causal determinants of complex disease.

A genome-wide association analysis identifies 16 novel susceptibility loci for carpal tunnel syndrome

Carpal tunnel syndrome (CTS) is a common and disabling condition of the hand caused by entrapment of the median nerve at the level of the wrist. It is the commonest entrapment neuropathy, with estimates of prevalence ranging between 5-10%. Here, we undertake a genome-wide association study (GWAS) of an entrapment neuropathy, using 12,312 CTS cases and 389,344 controls identified in UK Biobank. We discover 16 susceptibility loci for CTS with p < 5 × 10^-8. We identify likely causal genes in the pathogenesis of CTS, including ADAMTS17, ADAMTS10 and EFEMP1, and using RNA sequencing demonstrate expression of these genes in surgically resected tenosynovium from CTS patients. We perform Mendelian randomisation and demonstrate a causal relationship between short stature and higher risk of CTS. We suggest that variants within genes implicated in growth and extracellular matrix architecture contribute to the genetic predisposition to CTS by altering the environment through which the median nerve transits.

[Tumor banks and complex data management: Current and future challenges]

Tumor banks are asked to clinical and translationnal research project development in oncology. They strongly participate to the assessment, then to the validation of diagnostic, prognostic and predictive biomarkers. The progressive change of these structures leads to induce a professionalization of their functioning and to identify them as key actors in oncology by the stakeholders of the public and private worlds. The progresses made in biotechnologies and therapeutics are rapidly modifying the impact and the proper functioning of the biobanks. These latter are now facing different challenges, in particular for their sustainability. Among the major issues, the integration of the clinical and biological data becoming increasingly complex leads to urgently consider an optimization of the role of different biobanks in France. Their goal is to be an attractive counterpart face to the international competition. The purpose of this review is to briefly describe the current evolution of the biobanks, then their present and future challenges, and finally the role made by the pathologists in these new issues in oncology field.

Schizophrenia risk and reproductive success: a Mendelian randomization study

Schizophrenia is a debilitating and heritable mental disorder associated with lower reproductive success. However, the prevalence of schizophrenia is stable over populations and time, resulting in an evolutionary puzzle: how is schizophrenia maintained in the population, given its apparent fitness costs? One possibility is that increased genetic liability for schizophrenia, in the absence of the disorder itself, may confer some reproductive advantage. We assessed the correlation and causal effect of genetic liability for schizophrenia with number of children, age at first birth and number of sexual partners using data from the Psychiatric Genomics Consortium and UK Biobank. Linkage disequilibrium score regression showed little evidence of genetic correlation between genetic liability for schizophrenia and number of children (r g = 0.002, p = 0.84), age at first birth (r g = -0.007, p = 0.45) or number of sexual partners (r g = 0.007, p = 0.42). Mendelian randomization indicated no robust evidence of a causal effect of genetic liability for schizophrenia on number of children (mean difference: 0.003 increase in number of children per doubling in the natural log odds ratio of schizophrenia risk, 95% confidence interval (CI): -0.003 to 0.009, p = 0.39) or age at first birth (-0.004 years lower age at first birth, 95% CI: -0.043 to 0.034, p = 0.82). We find some evidence of a positive effect of genetic liability for schizophrenia on number of sexual partners (0.165 increase in the number of sexual partners, 95% CI: 0.117-0.212, p = 5.30×10-10). These results suggest that increased genetic liability for schizophrenia does not confer a fitness advantage but does increase mating success.

Biological and clinical insights from genetics of insomnia symptoms

Insomnia is a common disorder linked with adverse long-term medical and psychiatric outcomes. The underlying pathophysiological processes and causal relationships of insomnia with disease are poorly understood. Here we identified 57 loci for self-reported insomnia symptoms in the UK Biobank (n = 453,379) and confirmed their effects on self-reported insomnia symptoms in the HUNT Study (n = 14,923 cases and 47,610 controls), physician-diagnosed insomnia in the Partners Biobank (n = 2,217 cases and 14,240 controls), and accelerometer-derived measures of sleep efficiency and sleep duration in the UK Biobank (n = 83,726). Our results suggest enrichment of genes involved in ubiquitin-mediated proteolysis and of genes expressed in multiple brain regions, skeletal muscle, and adrenal glands. Evidence of shared genetic factors was found between frequent insomnia symptoms and restless legs syndrome, aging, and cardiometabolic, behavioral, psychiatric, and reproductive traits. Evidence was found for a possible causal link between insomnia symptoms and coronary artery disease, depressive symptoms, and subjective well-being.

Efficient Implementation of Penalized Regression for Genetic Risk Prediction

Polygenic risk scores (PRS) combine many single-nucleotide polymorphisms into a score reflecting the genetic risk of developing a disease. Privé, Aschard, and Blum present an efficient implementation of penalized logistic regression...

Polygenic Risk Scores (PRS) combine genotype information across many single-nucleotide polymorphisms (SNPs) to give a score reflecting the genetic risk of developing a disease. PRS might have a major impact on public health, possibly allowing for screening campaigns to identify high-genetic risk individuals for a given disease. The “Clumping+Thresholding” (C+T) approach is the most common method to derive PRS. C+T uses only univariate genome-wide association studies (GWAS) summary statistics, which makes it fast and easy to use. However, previous work showed that jointly estimating SNP effects for computing PRS has the potential to significantly improve the predictive performance of PRS as compared to C+T. In this paper, we present an efficient method for the joint estimation of SNP effects using individual-level data, allowing for practical application of penalized logistic regression (PLR) on modern datasets including hundreds of thousands of individuals. Moreover, our implementation of PLR directly includes automatic choices for hyper-parameters. We also provide an implementation of penalized linear regression for quantitative traits. We compare the performance of PLR, C+T and a derivation of random forests using both real and simulated data. Overall, we find that PLR achieves equal or higher predictive performance than C+T in most scenarios considered, while being scalable to biobank data. In particular, we find that improvement in predictive performance is more pronounced when there are few effects located in nearby genomic regions with correlated SNPs; for instance, in simulations, AUC values increase from 83% with the best prediction of C+T to 92.5% with PLR. We confirm these results in a data analysis of a case-control study for celiac disease where PLR and the standard C+T method achieve AUC values of 89% and of 82.5%. Applying penalized linear regression to 350,000 individuals of the UK Biobank, we predict height with a larger correlation than with the best prediction of C+T (∼65% instead of ∼55%), further demonstrating its scalability and strong predictive power, even for highly polygenic traits. Moreover, using 150,000 individuals of the UK Biobank, we are able to predict breast cancer better than C+T, fitting PLR in a few minutes only. In conclusion, this paper demonstrates the feasibility and relevance of using penalized regression for PRS computation when large individual-level datasets are available, thanks to the efficient implementation available in our R package bigstatsr.

Saliva molecular inflammatory profiling in female migraine patients responsive to adjunctive cervical non-invasive vagus nerve stimulation: the MOXY Study

Background

Rising evidence indicate that oxytocin and IL-1β impact trigemino-nociceptive signaling. Current perspectives on migraine physiopathology emphasize a cytokine bias towards a pro-inflammatory status. The anti-nociceptive impact of oxytocin has been reported in preclinical and human trials. Cervical non-invasive vagus nerve stimulation (nVNS) emerges as an add-on treatment for the preventive and abortive use in migraine. Less is known about its potential to modulate saliva inflammatory signaling in migraine patients. The rationale was to perform inter-ictal saliva measures of oxytocin and IL-1ß along with headache assessment in migraine patients with 10 weeks adjunctive nVNS compared to healthy controls.

Methods

12 migraineurs and 12 suitably matched healthy control were studied with inter-ictal saliva assay of pro- and anti-neuroinflammatory cytokines using enzyme-linked immuno assay techniques along with assessment of headache severity/frequency and associated functional capacity at baseline and after 10 weeks adjunctive cervical nVNS.

Results

nVNS significantly reduced headache severity (VAS), frequency (headache days and total number of attacks) and significantly improved sleep quality compared to baseline (p < 0.01). Inter-ictal saliva oxytocin and IL-1β were significantly elevated pre- as well as post-nVNS compared to healthy controls (p < 0.01) and similarly showed changes that may reflect the observed clinical effects.

Conclusions

Our results add to accumulating evidence for a therapeutic efficacy of adjunct cervical non-invasive vagus nerve stimulation in migraine patients. This study failed to provide an evidence-derived conclusion addressed to the predictive value and usefulness of saliva assays due to its uncontrolled study design. However, saliva screening of mediators associated with trigemino-nociceptive traffic represents a novel approach, thus deserve future targeted headache research.

Trial registration This study was indexed at the German Register for Clinical Trials (DRKS No. 00011089) registered on 21.09.2016

Biological insights into multiple birth: genetic findings from UK Biobank

The tendency to conceive spontaneous dizygotic (DZ) twins is a complex trait with important contributions from both environmental factors and genetic disposition. In earlier work, we identified the first two genes as maternal susceptibility loci for DZ twinning. The aim of this study was to identify genetic variants influencing multiple births and to genetically correlate the findings across a broad range of traits. We performed a genome-wide association study (GWAS) in 8962 participants with Caucasian ancestry from UK Biobank who reported being part of a multiple birth, and 409,591 singleton controls. We replicated the association between FSHB, SMAD3 and twinning in the gene-based (but not SNP-based) test, which had been established in previous genome-wide association analyses in mothers with dizygotic twin offspring. Additionally, we report a novel genetic variant associated with multiple birth, rs428022 at 15q23 (p = 2.84 × 10^-8) close to two genes: PIAS1 and SKOR1. Finally, we identified meaningful genetic correlations between being part of a multiple birth and other phenotypes (anthropometric traits, health-related traits, and fertility-related measures). The outcomes of this study provide important new insights into the genetic aetiology of multiple births and fertility, and open up novel directions for fertility and reproduction research.

Imputation of behavioral candidate gene repeat variants in 486,551 publicly-available UK Biobank individuals

Some of the most widely studied variants in psychiatric genetics include variable number tandem repeat variants (VNTRs) in SLC6A3, DRD4, SLC6A4, and MAOA. While initial findings suggested large effects, their importance with respect to psychiatric phenotypes is the subject of much debate with broadly conflicting results. Despite broad interest, these loci remain absent from the largest available samples, such as the UK Biobank, limiting researchers' ability to test these contentious hypotheses rigorously in large samples. Here, using two independent reference datasets, we report out-of-sample imputation accuracy estimates of >0.96 for all four VNTR variants and one modifying SNP, depending on the reference and target dataset. We describe the imputation procedures of these candidate variants in 486,551 UK Biobank individuals, and have made the imputed variant data available to UK Biobank researchers. This resource, provided to the scientific community, will allow the most rigorous tests to-date of the roles of these variants in behavioral and psychiatric phenotypes.

Genetic Data from Nearly 63,000 Women of European Descent Predicts DNA Methylation Biomarkers and Epithelial Ovarian Cancer Risk

DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P < 7.94 × 10-7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. SIGNIFICANCE: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.

Prioritizing putative influential genes in cardiovascular disease susceptibility by applying tissue-specific Mendelian randomization

BACKGROUND

The extent to which changes in gene expression can influence cardiovascular disease risk across different tissue types has not yet been systematically explored. We have developed an analysis pipeline that integrates tissue-specific gene expression, Mendelian randomization and multiple-trait colocalization to develop functional mechanistic insight into the causal pathway from a genetic variant to a complex trait.

METHODS

We undertook an expression quantitative trait loci-wide association study to uncover genetic variants associated with both nearby gene expression and cardiovascular traits. Fine-mapping was performed to prioritize possible causal variants for detected associations. Two-sample Mendelian randomization (MR) was then applied using findings from genome-wide association studies (GWAS) to investigate whether changes in gene expression within certain tissue types may influence cardiovascular trait variation. We subsequently used Bayesian multiple-trait colocalization to further interrogate the findings and also gain insight into whether DNA methylation, as well as gene expression, may play a role in disease susceptibility. Finally, we applied our analysis pipeline genome-wide using summary statistics from large-scale GWAS.

RESULTS

Eight genetic loci were associated with changes in gene expression and measures of cardiovascular function. Our MR analysis provided evidence of tissue-specific effects at multiple loci, of which the effects at the ADCY3 and FADS1 loci for body mass index and cholesterol, respectively, were particularly insightful. Multiple-trait colocalization uncovered evidence which suggested that changes in DNA methylation at the promoter region upstream of FADS1/TMEM258 may also affect cardiovascular trait variation along with gene expression. Furthermore, colocalization analyses uncovered evidence of tissue specificity between gene expression in liver tissue and cholesterol levels. Applying our pipeline genome-wide using summary statistics from GWAS uncovered 233 association signals at loci which represent promising candidates for further evaluation.

CONCLUSIONS

Disease susceptibility can be influenced by differential changes in tissue-specific gene expression and DNA methylation. The approach undertaken in our study can be used to elucidate mechanisms in disease, as well as helping prioritize putative causal genes at associated loci where multiple nearby genes may be co-regulated. Future studies which continue to uncover quantitative trait loci for molecular traits across various tissue and cell types will further improve our capability to understand and prevent disease.

Identification of 12 genetic loci associated with human healthspan

Aging populations face diminishing quality of life due to increased disease and morbidity. These challenges call for longevity research to focus on understanding the pathways controlling healthspan. We use the data from the UK Biobank (UKB) cohort and observe that the risks of major chronic diseases increased exponentially and double every eight years, i.e., at a rate compatible with the Gompertz mortality law. Assuming that aging drives the acceleration in morbidity rates, we build a risk model to predict the age at the end of healthspan depending on age, gender, and genetic background. Using the sub-population of 300,447 British individuals as a discovery cohort, we identify 12 loci associated with healthspan at the whole-genome significance level. We find strong genetic correlations between healthspan and all-cause mortality, life-history, and lifestyle traits. We thereby conclude that the healthspan offers a promising new way to interrogate the genetics of human longevity.

Mendelian randomization provides support for obesity as a risk factor for meningioma

Little is known about the causes of meningioma. Obesity and obesity-related traits have been reported in several epidemiological observational studies to be risk factors for meningioma. We performed an analysis of genetic variants associated with obesity-related traits to assess the relationship with meningioma risk using Mendelian randomization (MR), an approach unaffected by biases from temporal variability and reverse causation that might have affected earlier investigations. We considered 11 obesity-related traits, identified genetic instruments for these factors, and assessed their association with meningioma risk using data from a genome-wide association study comprising 1,606 meningioma patients and 9,823 controls. To evaluate the causal relationship between the obesity-related traits and meningioma risk, we consider the estimated odds ratio (OR) of meningioma for each genetic instrument. We identified positive associations between body mass index (odds ratio [OR_SD] = 1.27, 95% confidence interval [CI] = 1.03–1.56, P = 0.028) and body fat percentage (OR_SD = 1.28, 95% CI = 1.01–1.63, P = 0.042) with meningioma risk, albeit non-significant after correction for multiple testing. Associations for basal metabolic rate, diastolic blood pressure, fasting glucose, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, systolic blood pressure, total cholesterol, triglycerides and waist circumference with risk of meningioma were non-significant. Our analysis provides additional support for obesity being associated with an increased risk of meningioma.

Integrating predicted transcriptome from multiple tissues improves association detection

Integration of genome-wide association studies (GWAS) and expression quantitative trait loci (eQTL) studies is needed to improve our understanding of the biological mechanisms underlying GWAS hits, and our ability to identify therapeutic targets. Gene-level association methods such as PrediXcan can prioritize candidate targets. However, limited eQTL sample sizes and absence of relevant developmental and disease context restrict our ability to detect associations. Here we propose an efficient statistical method (MultiXcan) that leverages the substantial sharing of eQTLs across tissues and contexts to improve our ability to identify potential target genes. MultiXcan integrates evidence across multiple panels using multivariate regression, which naturally takes into account the correlation structure. We apply our method to simulated and real traits from the UK Biobank and show that, in realistic settings, we can detect a larger set of significantly associated genes than using each panel separately. To improve applicability, we developed a summary result-based extension called S-MultiXcan, which we show yields highly concordant results with the individual level version when LD is well matched. Our multivariate model-based approach allowed us to use the individual level results as a gold standard to calibrate and develop a robust implementation of the summary-based extension. Results from our analysis as well as software and necessary resources to apply our method are publicly available.

We develop a new method, MultiXcan, to test the mediating role of gene expression variation on complex traits, integrating information available across multiple tissue studies. We show this approach has higher power than traditional single-tissue methods. We extend this method to use only summary-statistics from public GWAS. We apply these methods to 222 complex traits available in the UK Biobank cohort, and 109 complex traits from public GWAS and discuss the findings.

A loss-of-function variant in ALOX15 protects against nasal polyps and chronic rhinosinusitis

Nasal polyps (NP) are lesions on the nasal and paranasal sinus mucosa and are a risk factor for chronic rhinosinusitis (CRS). We performed genome-wide association studies on NP and CRS in Iceland and the UK (using UK Biobank data) with 4,366 NP cases, 5,608 CRS cases, and >700,000 controls. We found 10 markers associated with NP and 2 with CRS. We also tested 210 markers reported to associate with eosinophil count, yielding 17 additional NP associations. Of the 27 NP signals, 7 associate with CRS and 13 with asthma. Most notably, a missense variant in ALOX15 that causes a p.Thr560Met alteration in arachidonate 15-lipoxygenase (15-LO) confers large genome-wide significant protection against NP (P = 8.0 × 10^-27, odds ratio = 0.32; 95% confidence interval = 0.26, 0.39) and CRS (P = 1.1 × 10^-8, odds ratio = 0.64; 95% confidence interval = 0.55, 0.75). p.Thr560Met, carried by around 1 in 20 Europeans, was previously shown to cause near total loss of 15-LO enzymatic activity. Our findings identify 15-LO as a potential target for therapeutic intervention in NP and CRS.

Balayage du génome et repérage des personnes à risque

The accumulation of extensive repositories linking phenotypic and genetic information, together with new computation methods, makes it possible to derive polygenic scores for susceptibility to common diseases that turn out to have strong predictive power. These will be clinically useful to identify individuals at high risk who may be eligible for protective interventions.

Reduced signal for polygenic adaptation of height in UK Biobank

Several recent papers have reported strong signals of selection on European polygenic height scores. These analyses used height effect estimates from the GIANT consortium and replication studies. Here, we describe a new analysis based on the the UK Biobank (UKB), a large, independent dataset. We find that the signals of selection using UKB effect estimates are strongly attenuated or absent. We also provide evidence that previous analyses were confounded by population stratification. Therefore, the conclusion of strong polygenic adaptation now lacks support. Moreover, these discrepancies highlight (1) that methods for correcting for population stratification in GWAS may not always be sufficient for polygenic trait analyses, and (2) that claims of differences in polygenic scores between populations should be treated with caution until these issues are better understood.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Artificial Intelligence and Personalized Medicine

The development of high-throughput, data-intensive biomedical research assays and technologies has created a need for researchers to develop strategies for analyzing, integrating, and interpreting the massive amounts of data they generate. Although a wide variety of statistical methods have been designed to accommodate 'big data,'  experiences with the use of artificial intelligence (AI) techniques suggest that they might be particularly appropriate. In addition,  the results of the application of these assays reveal a great heterogeneity in the pathophysiologic factors and processes that contribute to disease, suggesting that there is a need to tailor, or 'personalize,' medicines to the nuanced and often unique features possessed by individual patients. Given how important data-intensive assays are to revealing appropriate intervention targets and strategies for  treating an individual with a disease, AI can play an important role in the development of personalized medicines. We describe many areas where AI can play such a role and argue that AI's ability to advance personalized medicine will depend critically on not only the refinement of relevant assays, but also on ways of storing, aggregating, accessing, and ultimately integrating, the data they produce. We also point out the limitations of many AI techniques in developing personalized medicines as well as consider areas for further research.

Polygenic risk scores: a biased prediction?

A new study highlights the biases and inaccuracies of polygenic risk scores (PRS) when predicting disease risk in individuals from populations other than those used in their derivation. The design bias of workhorse tools used for research, particularly genotyping arrays, contributes to these distortions. To avoid further inequities in health outcomes, the inclusion of diverse populations in research, unbiased genotyping, and methods of bias reduction in PRS are critical.

Detecting the dominance component of heritability in isolated and outbred human populations

Inconsistencies between published estimates of dominance heritability between studies of human genetic isolates and human outbred populations incite investigation into whether such differences result from particular trait architectures or specific population structures. We analyse simulated datasets, characteristic of genetic isolates and of unrelated individuals, before analysing the isolate of Cilento for various commonly studied traits. We show the strengths of using genetic relationship matrices for variance decomposition over identity-by-descent based methods in a population isolate and that heritability estimates in isolates will avoid the downward biases that may occur in studies of samples of unrelated individuals; irrespective of the simulated distribution of causal variants. Yet, we also show that precise estimates of dominance in isolates are demonstrably problematic in the presence of shared environmental effects and such effects should be accounted for. Nevertheless, we demonstrate how studying isolates can help determine the existence or non-existence of dominance for complex traits, and we find strong indications of non-zero dominance for low-density lipoprotein level in Cilento. Finally, we recommend future study designs to analyse trait variance decomposition from ensemble data across multiple population isolates.

Dissection of genetic variation and evidence for pleiotropy in male pattern baldness

Male pattern baldness (MPB) is a sex-limited, age-related, complex trait. We study MPB genetics in 205,327 European males from the UK Biobank. Here we show that MPB is strongly heritable and polygenic, with pedigree-heritability of 0.62 (SE = 0.03) estimated from close relatives, and SNP-heritability of 0.39 (SE = 0.01) from conventionally-unrelated males. We detect 624 near-independent genome-wide loci, contributing SNP-heritability of 0.25 (SE = 0.01), of which 26 X-chromosome loci explain 11.6%. Autosomal genetic variance is enriched for common variants and regions of lower linkage disequilibrium. We identify plausible genetic correlations between MPB and multiple sex-limited markers of earlier puberty, increased bone mineral density (r_g = 0.15) and pancreatic β-cell function (r_g = 0.12). Correlations with reproductive traits imply an effect on fitness, consistent with an estimated linear selection gradient of -0.018 per MPB standard deviation. Overall, we provide genetic insights into MPB: a phenotype of interest in its own right, with value as a model sex-limited, complex trait.

Male pattern baldness (MPB) is a polygenic trait that affects the majority of European men. Here, Yap et al. estimate heritability, partitioned by autosomes and the X-chromosome, of MPB in the UK Biobank cohort, perform GWAS for MPB and find genetic correlation with other sex-specific traits.

GWAS identifies 14 loci for device-measured physical activity and sleep duration

Physical activity and sleep duration are established risk factors for many diseases, but their aetiology is poorly understood, partly due to relying on self-reported evidence. Here we report a genome-wide association study (GWAS) of device-measured physical activity and sleep duration in 91,105 UK Biobank participants, finding 14 significant loci (7 novel). These loci account for 0.06% of activity and 0.39% of sleep duration variation. Genome-wide estimates of ~ 15% phenotypic variation indicate high polygenicity. Heritability is higher in women than men for overall activity (23 vs. 20%, p = 1.5 × 10-4) and sedentary behaviours (18 vs. 15%, p = 9.7 × 10-4). Heritability partitioning, enrichment and pathway analyses indicate the central nervous system plays a role in activity behaviours. Two-sample Mendelian randomisation suggests that increased activity might causally lower diastolic blood pressure (beta mmHg/SD: -0.91, SE = 0.18, p = 8.2 × 10-7), and odds of hypertension (Odds ratio/SD: 0.84, SE = 0.03, p = 4.9 × 10-8). Our results advocate the value of physical activity for reducing blood pressure.

Genome-wide study of hair colour in UK Biobank explains most of the SNP heritability

Natural hair colour within European populations is a complex genetic trait. Previous work has established that MC1R variants are the principal genetic cause of red hair colour, but with variable penetrance. Here, we have extensively mapped the genes responsible for hair colour in the white, British ancestry, participants in UK Biobank. MC1R only explains 73% of the SNP heritability for red hair in UK Biobank, and in fact most individuals with two MC1R variants have blonde or light brown hair. We identify other genes contributing to red hair, the combined effect of which accounts for ~90% of the SNP heritability. Blonde hair is associated with over 200 genetic variants and we find a continuum from black through dark and light brown to blonde and account for 73% of the SNP heritability of blonde hair. Many of the associated genes are involved in hair growth or texture, emphasising the cellular connections between keratinocytes and melanocytes in the determination of hair colour.

Natural hair colour in Europeans is a complex genetic trait. Here, the authors carry out a genome-wide association study using UK BioBank data, suggesting that in combination with pigmentation genes, variants with roles in hair texture and growth can affect hair colouration or our perception of it.

Imprint of assortative mating on the human genome

Preference for mates with similar phenotypes; that is, assortative mating, is widely observed in humans1-5 and has evolutionary consequences6-8. Under Fisher's classical theory6, assortative mating is predicted to induce a signature in the genome at trait-associated loci that can be detected and quantified. Here, we develop and apply a method to quantify assortative mating on a specific trait by estimating the correlation (θ) between genetic predictors of the trait from single nucleotide polymorphisms on odd- versus even-numbered chromosomes. We show by theory and simulation that the effect of assortative mating can be quantified in the presence of population stratification. We applied this approach to 32 complex traits and diseases using single nucleotide polymorphism data from ~400,000 unrelated individuals of European ancestry. We found significant evidence of assortative mating for height (θ = 3.2%) and educational attainment (θ = 2.7%), both of which were consistent with theoretical predictions. Overall, our results imply that assortative mating involves multiple traits and affects the genomic architecture of loci that are associated with these traits, and that the consequence of mate choice can be detected from a random sample of genomes.

Genome-wide association meta-analysis yields 20 loci associated with gallstone disease

Gallstones are responsible for one of the most common diseases in the Western world and are commonly treated with cholecystectomy. We perform a meta-analysis of two genome-wide association studies of gallstone disease in Iceland and the UK, totaling 27,174 cases and 736,838 controls, uncovering 21 novel gallstone-associated variants at 20 loci. Two distinct low frequency missense variants in SLC10A2, encoding the apical sodium-dependent bile acid transporter (ASBT), associate with an increased risk of gallstone disease (Pro290Ser: OR = 1.36 [1.25-1.49], P = 2.1 × 10-12, MAF = 1%; Val98Ile: OR = 1.15 [1.10-1.20], P = 1.8 × 10-10, MAF = 4%). We demonstrate that lower bile acid transport by ASBT is accompanied by greater risk of gallstone disease and highlight the role of the intestinal compartment of the enterohepatic circulation of bile acids in gallstone disease susceptibility. Additionally, two low frequency missense variants in SERPINA1 and HNF4A and 17 common variants represent novel associations with gallstone disease.

Are Bigger Brains Smarter? Evidence From a Large-Scale Preregistered Study

A positive relationship between brain volume and intelligence has been suspected since the 19th century, and empirical studies seem to support this hypothesis. However, this claim is controversial because of concerns about publication bias and the lack of systematic control for critical confounding factors (e.g., height, population structure). We conducted a preregistered study of the relationship between brain volume and cognitive performance using a new sample of adults from the United Kingdom that is about 70% larger than the combined samples of all previous investigations on this subject ( N = 13,608). Our analyses systematically controlled for sex, age, height, socioeconomic status, and population structure, and our analyses were free of publication bias. We found a robust association between total brain volume and fluid intelligence ( r = .19), which is consistent with previous findings in the literature after controlling for measurement quality of intelligence in our data. We also found a positive relationship between total brain volume and educational attainment ( r = .12). These relationships were mainly driven by gray matter (rather than white matter or fluid volume), and effect sizes were similar for both sexes and across age groups.

Genome-wide association studies of brain imaging phenotypes in UK Biobank

The genetic architecture of brain structure and function is largely unknown. To investigate this, we carried out genome-wide association studies of 3,144 functional and structural brain imaging phenotypes from UK Biobank (discovery dataset 8,428 subjects). Here we show that many of these phenotypes are heritable. We identify 148 clusters of associations between single nucleotide polymorphisms and imaging phenotypes that replicate at P < 0.05, when we would expect 21 to replicate by chance. Notable significant, interpretable associations include: iron transport and storage genes, related to magnetic susceptibility of subcortical brain tissue; extracellular matrix and epidermal growth factor genes, associated with white matter micro-structure and lesions; genes that regulate mid-line axon development, associated with organization of the pontine crossing tract; and overall 17 genes involved in development, pathway signalling and plasticity. Our results provide insights into the genetic architecture of the brain that are relevant to neurological and psychiatric disorders, brain development and ageing.

Association of Whole-Genome and NETRIN1 Signaling Pathway-Derived Polygenic Risk Scores for Major Depressive Disorder and White Matter Microstructure in the UK Biobank

Background

Major depressive disorder is a clinically heterogeneous psychiatric disorder with a polygenic architecture. Genome-wide association studies have identified a number of risk-associated variants across the genome and have reported growing evidence of NETRIN1 pathway involvement. Stratifying disease risk by genetic variation within the NETRIN1 pathway may provide important routes for identification of disease mechanisms by focusing on a specific process, excluding heterogeneous risk-associated variation in other pathways. Here, we sought to investigate whether major depressive disorder polygenic risk scores derived from the NETRIN1 signaling pathway (NETRIN1-PRSs) and the whole genome, excluding NETRIN1 pathway genes (genomic-PRSs), were associated with white matter microstructure.

Methods

We used two diffusion tensor imaging measures, fractional anisotropy (FA) and mean diffusivity (MD), in the most up-to-date UK Biobank neuroimaging data release (FA: n = 6401; MD: n = 6390).

Results

We found significantly lower FA in the superior longitudinal fasciculus (β = −.035, p_corrected = .029) and significantly higher MD in a global measure of thalamic radiations (β = .029, p_corrected = .021), as well as higher MD in the superior (β = .034, p_corrected = .039) and inferior (β = .029, p_corrected = .043) longitudinal fasciculus and in the anterior (β = .025, p_corrected = .046) and superior (β = .027, p_corrected = .043) thalamic radiation associated with NETRIN1-PRS. Genomic-PRS was also associated with lower FA and higher MD in several tracts.

Conclusions

Our findings indicate that variation in the NETRIN1 signaling pathway may confer risk for major depressive disorder through effects on a number of white matter tracts.

An efficient linear mixed model framework for meta-analytic association studies across multiple contexts

Linear mixed models (LMMs) can be applied in the meta-analyses of responses from individuals across multiple contexts, increasing power to detect associations while accounting for confounding effects arising from within-individual variation. However, traditional approaches to fitting these models can be computationally intractable. Here, we describe an efficient and exact method for fitting a multiple-context linear mixed model. Whereas existing exact methods may be cubic in their time complexity with respect to the number of individuals, our approach for multiple-context LMMs (mcLMM) is linear. These improvements allow for large-scale analyses requiring computing time and memory magnitudes of order less than existing methods. As examples, we apply our approach to identify expression quantitative trait loci from large-scale gene expression data measured across multiple tissues as well as joint analyses of multiple phenotypes in genome-wide association studies at biobank scale.