Genome-wide association studies (GWAS) have identified many loci contributing to variation in complex traits, yet the majority of loci that contribute to the heritability of complex traits remain ...elusive. Large study populations with sufficient statistical power are required to detect the small effect sizes of the yet unidentified genetic variants. However, the analysis of huge cohorts, like UK Biobank, is challenging. Here, we present an atlas of genetic associations for 118 non-binary and 660 binary traits of 452,264 UK Biobank participants of European ancestry. Results are compiled in a publicly accessible database that allows querying genome-wide association results for 9,113,133 genetic variants, as well as downloading GWAS summary statistics for over 30 million imputed genetic variants (>23 billion phenotype-genotype pairs). Our atlas of associations (GeneATLAS, http://geneatlas.roslin.ed.ac.uk ) will help researchers to query UK Biobank results in an easy and uniform way without the need to incur high computational costs.
Association between whole blood viscosity (WBV) and an increased risk of cardiovascular disease (CVD) has been reported. However, the causal relationship between WBV and CVD remains not thoroughly ...investigated. The aim of this study was to investigate the causal relation between WBV and CVD.
Two-sample Mendelian randomization (MR) was employed, with inverse variance weighting (IVW) as the primary method, to investigate the casual relationship between WBV and CVD. The calculated WBV and medical records of 378,210 individuals participating in the UK Biobank study were divided into halves and analyzed.
The means of calculated WBVs were 16.9 (standard deviation: 0.8) and 55.1 (standard deviation: 17.2) for high shear rate (HSR) and low shear rate (LSR), respectively. 37,859 (10.0%) major cardiovascular events (MACE) consisted of 23,894 (6.3%) cases of myocardial infarction (MI), 9,245 (2.4%) cases of ischemic stroke, 10,377 (2.7%) cases of revascularization, and 5,703 (1.5%) cases of coronary heart disease-related death. In the MR analysis, no evidence was found indicating a causal effect of WBV on MACE (IVW p-value for HSR = 0.81, IVW p-value for LSR = 0.47), MI (0.92, 0.83), ischemic stroke (0.52, 0.74), revascularization (0.71, 0.54), and coronary heart disease-related death (0.83, 0.70). The lack of sufficient evidence for causality persisted in other MR methods, including weighted median and MR-egger.
The Mendelian randomization analysis conducted in this study does not support a causal relationship between calculated WBV and CVD.
Relatives provide the basic material for the study of inheritance of human disease. However, the methodologies for the estimation of heritability and the interpretation of the results have been ...controversial. The debate arises from the plethora of methods used, the validity of the methodological assumptions and the inconsistent and sometimes erroneous genetic interpretations made. We will discuss how to estimate disease heritability, how to interpret it, how biases in heritability estimates arise and how heritability relates to other measures of familial disease aggregation.
Sex differences are a common feature of human traits; however, the role sex determination plays in human genetic variation remains unclear. The presence of gene-by-sex (GxS) interactions implies that ...trait genetic architecture differs between men and women. Here, we show that GxS interactions and genetic heterogeneity among sexes are small but common features of a range of high-level complex traits.
We analyzed 19 complex traits measured in 54,040 unrelated men and 59,820 unrelated women from the UK Biobank cohort to estimate autosomal genetic correlations and heritability differences between men and women. For 13 of the 19 traits examined, there is evidence that the trait measured is genetically different between males and females. We find that estimates of genetic correlations, based on ~114,000 unrelated individuals and ~19,000 related individuals from the same cohort, are largely consistent. Genetic predictors using a sex-specific model that incorporated GxS interactions led to a relative improvement of up to 4 % (mean 1.4 % across all relevant phenotypes) over those provided by a sex-agnostic model. This further supports the hypothesis of the presence of sexual genetic heterogeneity across high-level phenotypes.
The sex-specific environment seems to play a role in changing genotype expression across a range of human complex traits. Further studies of GxS interactions for high-level human traits may shed light on the molecular mechanisms that lead to biological differences between men and women. However, this may be a challenging endeavour due to the likely small effects of the interactions at individual loci.
Phenotypic correlations among partners for traits such as longevity or late-onset disease have been found to be comparable to phenotypic correlations in first-degree relatives. How these correlations ...arise in late life is poorly understood. Here we introduce a novel paradigm to establish the presence of indirect assortment on factors correlated across generations, by examining correlations between parents of couples, i.e., in-laws. Using correlations in additive genetic values we further corroborate the presence of indirect assortment on heritable factors. Specifically, using couples from the UK Biobank cohort, we show that longevity and disease history of the parents of White British couples are correlated, with correlations of up to 0.09. The correlations in parental longevity are replicated in the FamiLinx cohort, a larger and geographically more diverse historical ancestry dataset spanning a broader time frame. These correlations in parental longevity significantly (pval < 0.0093 for all pairs of parents) exceed what would be expected due to variations in lifespan based on year and location of birth. For cardiovascular diseases, in particular hypertension, we find significant correlations (r = 0.028, pval = 0.005) in genetic values among partners, supporting a model where partners assort for risk factors to some extent genetically correlated with cardiovascular disease. Partitioning the relative importance of indirect assortative mating and shared common environment will require large, well-characterized longitudinal cohorts aimed at understanding phenotypic correlations among none-blood relatives. Identifying the factors that mediate indirect assortment on longevity and human disease risk will help to unravel factors affecting human disease and ultimately longevity.
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.
We investigated whether the plasma level of 25-hydroxyvitamin D (25-OHD) after a diagnosis of colorectal cancer (CRC) influences survival outcome.
We prospectively studied 1,598 patients with stage I ...to III CRC. We sought association between plasma 25-OHD and stage-specific survival and tested for interaction between 25-OHD level and variation at the vitamin D receptor (VDR) gene locus. Blood was sampled postoperatively, and plasma was assayed for 25-OHD by liquid chromatography-tandem mass spectrometry. VDR polymorphisms (rs1544410, rs10735810, rs7975232, rs11568820) were genotyped, and haplotypes were inferred by using BEAGLE software. We tested for association between survival and 25-OHD, VDR genotype/haplotype, and after applying a VDR genotype-25-OHD interaction term. We conducted Kaplan-Meier survival analysis and used Cox proportional hazards models to estimate adjusted hazard ratios (HRs).
We found strong associations between plasma 25-OHD concentration and CRC-specific (P = .008) and all-cause mortality (P = .003). Adjusted HRs were 0.68 (95% CI, 0.50 to 0.90) and 0.70 (95% CI, 0.55 to 0.89), respectively (highest v lowest 25-OHD tertile), particularly in stage II disease (HR, 0.44; P = .004 for CRC-specific mortality). We detected gene-environment interactions between 25-OHD concentration and rs11568820 genotype for CRC-specific (P = .008) and all-cause (P = .022) mortality, number of protective alleles (P = .004 and P = .018, respectively), and GAGC haplotype at the VDR locus for all-cause mortality (P = .008).
In patients with stage I to III CRC, postoperative plasma vitamin D is associated with clinically important differences in survival outcome, higher levels being associated with better outcome. We observed interactions between 25-OHD level and VDR genotype, suggesting a causal relationship between vitamin D and survival. The influence of vitamin D supplementation on CRC outcome will require further investigation.
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.
Previous studies have shown that aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) lower colorectal cancer (CRC) risk. However, the lowest effective NSAID dose, treatment duration, and ...effects on survival are not defined. In a large population-based case-control study, we have explored the relationship between NSAID dose and duration, CRC risk and overall CRC-specific survival.
The relationship between NSAID use and CRC risk was examined in 2279 cases and 2907 controls. Subjects completed food-frequency and lifestyle questionnaires. NSAID categories were low-dose aspirin (75 mg), non-aspirin NSAIDs (NA-NSAIDs) and any NSAID. Users were defined as taking >4 tablets/week for >1 month. ORs were calculated by logistic regression models and adjusted for potential confounding factors. Effect of NSAID use on all-cause and CRC-specific mortality was estimated using Logrank tests and Cox's hazard models.
In all, 354 cases (15.5%) were taking low-dose aspirin compared to 526 controls (18.1%). Low-dose aspirin use was associated with decreased CRC risk (OR 0.78 95% CI 0.65 to 0.92, p=0.004), evident after 1 year and increasing with duration of use (p(trend)=0.004). NA-NSAID and any NSAID use were also inversely associated with CRC. There was no demonstrable effect of NSAIDS on all-cause (HR 1.11, p=0.22, 0.94-1.33) or CRC-specific survival (HR 1.01, p=0.93, 0.83-1.23).
This is the first study to demonstrate a protective effect against CRC associated with the lowest dose of aspirin (75 mg per day) after only 5 years use in the general population. NSAID use prior to CRC diagnosis does not influence survival from the disease.
The functional annotation of genomes, including chromatin accessibility and modifications, is important for understanding and effectively utilizing the increased amount of genome sequences reported. ...However, while such annotation has been well explored in a diverse set of tissues and cell types in human and model organisms, relatively little data are available for livestock genomes, hindering our understanding of complex trait variation, domestication, and adaptive evolution. Here, we present the first complete global landscape of regulatory elements in cattle and explore the dynamics of chromatin states in rumen epithelial cells induced by the rumen developmental regulator-butyrate.
We established the first global map of regulatory elements (15 chromatin states) and defined their coordinated activities in cattle, through genome-wide profiling for six histone modifications, RNA polymerase II, CTCF-binding sites, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial primary cells (REPC), rumen tissues, and Madin-Darby bovine kidney epithelial cells (MDBK). We demonstrated that each chromatin state exhibited specific enrichment for sequence ontology, transcription, methylation, trait-associated variants, gene expression-associated variants, selection signatures, and evolutionarily conserved elements, implying distinct biological functions. After butyrate treatments, we observed that the weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states, occurred concomitantly with significant alterations in gene expression and DNA methylation, which was significantly associated with heifer conception rate and stature economic traits.
Our results demonstrate the crucial role of functional genome annotation for understanding genome regulation, complex trait variation, and adaptive evolution in livestock. Using butyrate to induce the dynamics of the epigenomic landscape, we were able to establish the correlation among nutritional elements, chromatin states, gene activities, and phenotypic outcomes.