Phenome-wide analyses identify an association between the parent-of-origin effects dependent methylome and the rate of aging in humans.

BACKGROUND: The variation in the rate at which humans age may be rooted in early events acting through the genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions enriched for genetically controlled imprinting effects (the typical type of POE) and regions influenced by environmental effects associated with parents (the atypical POE). This part of the methylome is heavily influenced by early events, making it a potential route connecting early exposures, the epigenome, and aging. We aim to test the association of POE-CpGs with early and later exposures and subsequently with health-related phenotypes and adult aging. RESULTS: We perform a phenome-wide association analysis for the POE-influenced methylome using GS:SFHS (Ndiscovery = 5087, Nreplication = 4450). We identify and replicate 92 POE-CpG-phenotype associations. Most of the associations are contributed by the POE-CpGs belonging to the atypical class where the most strongly enriched associations are with aging (DNAmTL acceleration), intelligence, and parental (maternal) smoking exposure phenotypes. A proportion of the atypical POE-CpGs form co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased within-module methylation connectivity with age. The atypical POE-CpGs also display high levels of methylation heterogeneity, fast information loss with age, and a strong correlation with CpGs contained within epigenetic clocks. CONCLUSIONS: These results identify the association between the atypical POE-influenced methylome and aging and provide new evidence for the "early development of origin" hypothesis for aging in humans.

Phenome-wide analysis identifies parent-of-origin effects on the human methylome associated with changes in the rate of aging.

Variation in the rate at which humans age may be rooted in early life events acting through genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions either enriched for genetically controlled imprinting effects (the typical type of POE) or atypical POE introduced by environmental effects associated with parents. This part of the methylome is heavily influenced by early life events, making it a potential route connecting early environmental exposures, the epigenome and the rate of aging. Here, we aim to test the association of POE-influenced methylation of CpG dinucleotides (POE-CpG sites) with early and later environmental exposures and subsequently with health-related phenotypes and adult aging phenotypes. We do this by performing phenome-wide association analyses of the POE-influenced methylome using a large family-based population cohort (GS:SFHS, Ndiscovery=5,087, Nreplication=4,450). At the single CpG level, 92 associations of POE-CpGs with phenotypic variation were identified and replicated. Most of the associations were contributed by POE-CpGs belonging to the atypical class and the most strongly enriched associations were with aging (DNAmTL acceleration), intelligence and parental (maternal) smoking exposure phenotypes. We further found that a proportion of the atypical-POE-CpGs formed co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased internal module connectivity (strength of methylation correlation across constituent CpGs) with age. Atypical POE-CpGs also displayed high levels of methylation heterogeneity and epigenetic drift (i.e. information loss with age) and a strong correlation with CpGs contained within epigenetic clocks. These results identified associations between the atypical-POE-influenced methylome and aging and provided new evidence for the "early development of origin" hypothesis for aging in humans.

Genome-by-Trauma Exposure Interactions in Adults With Depression in the UK Biobank.

Importance: Self-reported trauma exposure has consistently been found to be a risk factor for major depressive disorder (MDD), and several studies have reported interactions with genetic liability. To date, most studies have examined gene-environment interactions with trauma exposure using genome-wide variants (single-nucleotide variations [SNVs]) or polygenic scores, both typically capturing less than 3% of phenotypic risk variance. Objective: To reexamine genome-by-trauma interaction associations using genetic measures using all available genotyped data and thus, maximizing accounted variance. Design, Setting, and Participants: The UK Biobank study was conducted from April 2007 to May 1, 2016 (follow-up mental health questionnaire). The current study used available cross-sectional genomic and trauma exposure data from UK Biobank. Participants who completed the mental health questionnaire and had available genetic, trauma experience, depressive symptoms, and/or neuroticism information were included. Data were analyzed from April 1 to August 30, 2021. Exposures: Trauma and genome-by-trauma exposure interactions. Main Outcomes and Measures: Measures of self-reported depression, neuroticism, and trauma exposure with whole-genome SNV data are available from the UK Biobank study. Here, a mixed-model statistical approach using genetic, trauma exposure, and genome-by-trauma exposure interaction similarity matrices was used to explore sources of variation in depression and neuroticism. Results: Analyses were conducted on 148 129 participants (mean [SD] age, 56 [7] years) of which 76 995 were female (52.0%). The study approach estimated the heritability (SE) of MDD to be approximately 0.160 (0.016). Subtypes of self-reported trauma exposure (catastrophic, adult, childhood, and full trauma) accounted for a significant proportion of the variance of MDD, with heritability (SE) ranging from 0.056 (0.013) to 0.176 (0.025). The proportion of MDD risk variance accounted for by significant genome-by-trauma interaction revealed estimates (SD) ranging from 0.074 (0.006) to 0.201 (0.009). Results from sex-specific analyses found genome-by-trauma interaction variance estimates approximately 5-fold greater for MDD in male participants (0.441 [0.018]) than in female participants (0.086 [0.009]). Conclusions and Relevance: This cross-sectional study used an approach combining all genome-wide SNV data when exploring genome-by-trauma interactions in individuals with MDD; findings suggest that such interactions were associated with depression manifestation. Genome-by-trauma interaction accounts for greater trait variance in male individuals, which points to potential differences in depression etiology between the sexes. The methodology used in this study can be extrapolated to other environmental factors to identify modifiable risk environments and at-risk groups to target with interventions.

Effect of genotyping strategies on the sustained benefit of single-step genomic BLUP over multiple generations.

BACKGROUND: Single-step genomic best linear unbiased prediction (ssGBLUP) allows the inclusion of information from genotyped and ungenotyped individuals in a single analysis. This avoids the need to genotype all candidates with the potential benefit of reducing overall costs. The aim of this study was to assess the effect of genotyping strategies, the proportion of genotyped candidates and the genotyping criterion to rank candidates to be genotyped, when using ssGBLUP evaluation. A simulation study was carried out assuming selection over several discrete generations where a proportion of the candidates were genotyped and evaluation was done using ssGBLUP. The scenarios compared were: (i) three genotyping strategies defined by their protocol for choosing candidates to be genotyped (RANDOM: candidates were chosen at random; TOP: candidates with the best genotyping criterion were genotyped; and EXTREME: candidates with the best and worse criterion were genotyped); (ii) eight proportions of genotyped candidates (p); and (iii) two genotyping criteria to rank candidates to be genotyped (candidates' own phenotype or estimated breeding values). The criteria of the comparison were the cumulated gain and reliability of the genomic estimated breeding values (GEBV). RESULTS: The genotyping strategy with the greatest cumulated gain was TOP followed by RANDOM, with EXTREME behaving as RANDOM at low p and as TOP with high p. However, the reliability of GEBV was higher with RANDOM than with TOP. This disparity between the trend of the gain and the reliability is due to the TOP scheme genotyping the candidates with the greater chances of being selected. The extra gain obtained with TOP increases when the accuracy of the selection criterion to rank candidates to be genotyped increases. CONCLUSIONS: The best strategy to maximise genetic gain when only a proportion of the candidates are to be genotyped is TOP, since it prioritises the genotyping of candidates which are more likely to be selected. However, the strategy with the greatest GEBV reliability does not achieve the largest gain, thus reliability cannot be considered as an absolute and sufficient criterion for determining the scheme which maximises genetic gain.

Lifestyle and Genetic Factors Modify Parent-of-Origin Effects on the Human Methylome.

BACKGROUND: parent-of-origin effects (POE) play important roles in complex disease and thus understanding their regulation and associated molecular and phenotypic variation are warranted. Previous studies mainly focused on the detection of genomic regions or phenotypes regulated by POE. Understanding whether POE may be modified by environmental or genetic exposures is important for understanding of the source of POE-associated variation, but only a few case studies addressing modifiable POE exist. METHODS: in order to understand this high order of POE regulation, we screened 101 genetic and environmental factors such as 'predicted mRNA expression levels' of DNA methylation/imprinting machinery genes and environmental exposures. POE-mQTL-modifier interaction models were proposed to test the potential of these factors to modify POE at DNA methylation using data from Generation Scotland: The Scottish Family Health Study(N=2315). FINDINGS: a set of vulnerable/modifiable POE-CpGs were identified (modifiable-POE-regulated CpGs, N=3). Four factors, 'lifetime smoking status' and 'predicted mRNA expression levels' of TET2, SIRT1 and KDM1A, were found to significantly modify the POE on the three CpGs in both discovery and replication datasets. We further identified plasma protein and health-related phenotypes associated with the methylation level of one of the identified CpGs. INTERPRETATION: the modifiable POE identified here revealed an important yet indirect path through which genetic background and environmental exposures introduce their effect on DNA methylation, motivating future comprehensive evaluation of the role of these modifiers in complex diseases. FUNDING: NSFC (81971270),H2020-MSCA-ITN(721815), Wellcome (204979/Z/16/Z,104036/Z/14/Z), MRC (MC_UU_00007/10, MC_PC_U127592696), CSO (CZD/16/6,CZB/4/276, CZB/4/710), SFC (HR03006), EUROSPAN (LSHG-CT-2006-018947), BBSRC (BBS/E/D/30002276), SYSU, Arthritis Research UK, NHLBI, NIH.

Genome-wide methylation data improves dissection of the effect of smoking on body mass index.

Variation in obesity-related traits has a genetic basis with heritabilities between 40 and 70%. While the global obesity pandemic is usually associated with environmental changes related to lifestyle and socioeconomic changes, most genetic studies do not include all relevant environmental covariates, so the genetic contribution to variation in obesity-related traits cannot be accurately assessed. Some studies have described interactions between a few individual genes linked to obesity and environmental variables but there is no agreement on their total contribution to differences between individuals. Here we compared self-reported smoking data and a methylation-based proxy to explore the effect of smoking and genome-by-smoking interactions on obesity related traits from a genome-wide perspective to estimate the amount of variance they explain. Our results indicate that exploiting omic measures can improve models for complex traits such as obesity and can be used as a substitute for, or jointly with, environmental records to better understand causes of disease.

Nontrivial Replication of Loci Detected by Multi-Trait Methods.

The ever-growing genome-wide association studies (GWAS) have revealed widespread pleiotropy. To exploit this, various methods that jointly consider associations of a genetic variant with multiple traits have been developed. Most efforts have been made concerning improving GWAS discovery power. However, how to replicate these discovered pleiotropic loci has yet to be discussed thoroughly. Unlike a single-trait scenario, multi-trait replication is not trivial considering the underlying genotype-multi-phenotype map of the associations. Here, we evaluate four methods for replicating multi-trait associations, corresponding to four levels of replication strength. Weak replication cannot justify pleiotropic genetic effects, whereas strong replication using our developed correlation methods can inform consistent pleiotropic genetic effects across the discovery and replication samples. We provide a protocol for replicating multi-trait genetic associations in practice. The described methods are implemented in the free and open-source R package MultiABEL.

Identification of epigenome-wide DNA methylation differences between carriers of APOE ε4 and APOE ε2 alleles.

BACKGROUND: The apolipoprotein E (APOE) ε4 allele is the strongest genetic risk factor for late onset Alzheimer's disease, whilst the ε2 allele confers protection. Previous studies report differential DNA methylation of APOE between ε4 and ε2 carriers, but associations with epigenome-wide methylation have not previously been characterised. METHODS: Using the EPIC array, we investigated epigenome-wide differences in whole blood DNA methylation patterns between Alzheimer's disease-free APOE ε4 (n = 2469) and ε2 (n = 1118) carriers from the two largest single-cohort DNA methylation samples profiled to date. Using a discovery, replication and meta-analysis study design, methylation differences were identified using epigenome-wide association analysis and differentially methylated region (DMR) approaches. Results were explored using pathway and methylation quantitative trait loci (meQTL) analyses. RESULTS: We obtained replicated evidence for DNA methylation differences in a ~ 169 kb region, which encompasses part of APOE and several upstream genes. Meta-analytic approaches identified DNA methylation differences outside of APOE: differentially methylated positions were identified in DHCR24, LDLR and ABCG1 (2.59 × 10-100 ≤ P ≤ 2.44 × 10-8) and DMRs were identified in SREBF2 and LDLR (1.63 × 10-4 ≤ P ≤ 3.01 × 10-2). Pathway and meQTL analyses implicated lipid-related processes and high-density lipoprotein cholesterol was identified as a partial mediator of the methylation differences in ABCG1 and DHCR24. CONCLUSIONS: APOE ε4 vs. ε2 carrier status is associated with epigenome-wide methylation differences in the blood. The loci identified are located in trans as well as cis to APOE and implicate genes involved in lipid homeostasis.

SNP and Haplotype Regional Heritability Mapping (SNHap-RHM): Joint Mapping of Common and Rare Variation Affecting Complex Traits.

We describe a genome-wide analytical approach, SNP and Haplotype Regional Heritability Mapping (SNHap-RHM), that provides regional estimates of the heritability across locally defined regions in the genome. This approach utilises relationship matrices that are based on sharing of SNP and haplotype alleles at local haplotype blocks delimited by recombination boundaries in the genome. We implemented the approach on simulated data and show that the haplotype-based regional GRMs capture variation that is complementary to that captured by SNP-based regional GRMs, and thus justifying the fitting of the two GRMs jointly in a single analysis (SNHap-RHM). SNHap-RHM captures regions in the genome contributing to the phenotypic variation that existing genome-wide analysis methods may fail to capture. We further demonstrate that there are real benefits to be gained from this approach by applying it to real data from about 20,000 individuals from the Generation Scotland: Scottish Family Health Study. We analysed height and major depressive disorder (MDD). We identified seven genomic regions that are genome-wide significant for height, and three regions significant at a suggestive threshold (p-value < 1 × 10-5) for MDD. These significant regions have genes mapped to within 400 kb of them. The genes mapped for height have been reported to be associated with height in humans. Similarly, those mapped for MDD have been reported to be associated with major depressive disorder and other psychiatry phenotypes. The results show that SNHap-RHM presents an exciting new opportunity to analyse complex traits by allowing the joint mapping of novel genomic regions tagged by either SNPs or haplotypes, potentially leading to the recovery of some of the "missing" heritability.

Author Correction: Bayesian reassessment of the epigenetic architecture of complex traits.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Bayesian reassessment of the epigenetic architecture of complex traits.

Linking epigenetic marks to clinical outcomes improves insight into molecular processes, disease prediction, and therapeutic target identification. Here, a statistical approach is presented to infer the epigenetic architecture of complex disease, determine the variation captured by epigenetic effects, and estimate phenotype-epigenetic probe associations jointly. Implicitly adjusting for probe correlations, data structure (cell-count or relatedness), and single-nucleotide polymorphism (SNP) marker effects, improves association estimates and in 9,448 individuals, 75.7% (95% CI 71.70-79.3) of body mass index (BMI) variation and 45.6% (95% CI 37.3-51.9) of cigarette consumption variation was captured by whole blood methylation array data. Pathway-linked probes of blood cholesterol, lipid transport and sterol metabolism for BMI, and xenobiotic stimuli response for smoking, showed >1.5 times larger associations with >95% posterior inclusion probability. Prediction accuracy improved by 28.7% for BMI and 10.2% for smoking over a LASSO model, with age-, and tissue-specificity, implying associations are a phenotypic consequence rather than causal.

Genomic regions influencing intramuscular fat in divergently selected rabbit lines.

Intramuscular fat (IMF) is one of the main meat quality traits for breeding programmes in livestock species. The main objective of this study was to identify genomic regions associated with IMF content comparing two rabbit populations divergently selected for this trait, and to generate a list of putative candidate genes. Animals were genotyped using the Affymetrix Axiom OrcunSNP Array (200k). After quality control, the data involved 477 animals and 93 540 SNPs. Two methods were used in this research: single marker regressions with the data adjusted by genomic relatedness, and a Bayesian multiple marker regression. Associated genomic regions were located on the rabbit chromosomes (OCU) OCU1, OCU8 and OCU13. The highest value for the percentage of the genomic variance explained by a genomic region was found in two consecutive genomic windows on OCU8 (7.34%). Genes in the associated regions of OCU1 and OCU8 presented biological functions related to the control of adipose cell function, lipid binding, transportation and localisation (APOLD1, PLBD1, PDE6H, GPRC5D and GPRC5A) and lipid metabolic processes (MTMR2). The EWSR1 gene, underlying the OCU13 region, is linked to the development of brown adipocytes. The findings suggest that there is a large component of polygenic effect behind the differences in IMF content in these two lines, as the variance explained by most of the windows was low. The genomic regions of OCU1, OCU8 and OCU13 revealed novel candidate genes. Further studies would be needed to validate the associations and explore their possible application in selection programmes.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions.

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

The genetic landscape of Scotland and the Isles.

Britain and Ireland are known to show population genetic structure; however, large swathes of Scotland, in particular, have yet to be described. Delineating the structure and ancestry of these populations will allow variant discovery efforts to focus efficiently on areas not represented in existing cohorts. Thus, we assembled genotype data for 2,554 individuals from across the entire archipelago with geographically restricted ancestry, and performed population structure analyses and comparisons to ancient DNA. Extensive geographic structuring is revealed, from broad scales such as a NE to SW divide in mainland Scotland, through to the finest scale observed to date: across 3 km in the Northern Isles. Many genetic boundaries are consistent with Dark Age kingdoms of Gaels, Picts, Britons, and Norse. Populations in the Hebrides, the Highlands, Argyll, Donegal, and the Isle of Man show characteristics of isolation. We document a pole of Norwegian ancestry in the north of the archipelago (reaching 23 to 28% in Shetland) which complements previously described poles of Germanic ancestry in the east, and "Celtic" to the west. This modern genetic structure suggests a northwestern British or Irish source population for the ancient Gaels that contributed to the founding of Iceland. As rarer variants, often with larger effect sizes, become the focus of complex trait genetics, more diverse rural cohorts may be required to optimize discoveries in British and Irish populations and their considerable global diaspora.

"Arte et Labore"-A Blackburn Rovers fan's legacy in human complex trait genetics.

Through his own research contributions on the modelling and genetic analysis of quantitative traits and through his former students and postdocs, Robin Thompson has indirectly left a major legacy in human genetics. In this short note, we highlight examples of the long-lasting relevance and impact of Robin's work in human genetics. A lone early study of marker-assisted selection developed many of the tools and approaches later exploited (often after reinvention) by the human genetics community in GWAS studies and for prediction. Furthermore, a particularly clear example of the pervasive impact of Robin's work is that REML has become the default method to estimate variance components and that genetic predictions exploiting linkage disequilibrium in the population are starting to become used in precision medicine applications.

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

In the original version of this Article, the legend in the upper panel of Figure 2 incorrectly read 'paternal imprinting' and should have read 'maternal imprinting'. This has been corrected in both the PDF and HTML versions of the Article.

Insulin resistance: Genetic associations with depression and cognition in population based cohorts.

Insulin resistance, broadly defined as the reduced ability of insulin to exert its biological action, has been associated with depression and cognitive dysfunction in observational studies. However, it is unclear whether these associations are causal and whether they might be underpinned by other shared factors. To address this knowledge gap, we capitalized on the stability of genetic biomarkers through the lifetime, and on their unidirectional relationship with depression and cognition. Specifically, we determined the association between quantitative measures of cognitive function and depression and genetic instruments of insulin resistance traits in two large-scale population samples, the Generation Scotland: Scottish Family Health Study (GS: SFHS; N = 19,994) and in the UK Biobank (N = 331,374). In the GS:SFHS, the polygenic risk score (PRS) for fasting insulin was associated with verbal intelligence and depression while the PRS for the homeostasis model assessment of insulin resistance was associated with verbal intelligence. Despite this overlap in genetic architecture, Mendelian randomization analyses in the GS:SFHS and in the UK Biobank samples did not yield evidence for causal associations from insulin resistance traits to either depression or cognition. These findings may be due to weak genetic instruments, limited cognitive measures and insufficient power but they may also indicate the need to identify other biological mechanisms that may mediate the relationship from insulin resistance to depression and cognition.

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.

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, pcorrected = .029) and significantly higher MD in a global measure of thalamic radiations (ß = .029, pcorrected = .021), as well as higher MD in the superior (ß = .034, pcorrected = .039) and inferior (ß = .029, pcorrected = .043) longitudinal fasciculus and in the anterior (ß = .025, pcorrected = .046) and superior (ß = .027, pcorrected = .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.