Search

Testing the effect of rare variants on phenotypic variation is difficult due to the need for extremely large cohorts to identify associated variants given expected effect sizes. An alternative approach is to investigate the effect of rare genetic variants on low-level genomic traits, such as gene expression or DNA methylation (DNAm), as effect sizes are expected to be larger for low-level compared to higher-order complex traits. Here, we investigate DNAm in healthy ageing populations - the Lothian Birth cohorts of 1921 and 1936 and identify both transient and stable outlying DNAm levels across the genome. We find an enrichment of rare genetic variants within 1kb of DNAm sites in individuals with stable outlying DNAm, implying genetic control of this extreme variation. Using a family-based cohort, the Brisbane Systems Genetics Study, we observed increased sharing of DNAm outliers among more closely related individuals, consistent with these outliers being driven by rare genetic variation. We demonstrated that outlying DNAm levels have a functional consequence on gene expression levels, with extreme levels of DNAm being associated with gene expression levels towards the tails of the population distribution. Overall, this study demonstrates the role of rare variants in the phenotypic variation of low-level genomic traits, and the effect of extreme levels of DNAm on gene expression.

BackgroundAmyotrophic lateral sclerosis (ALS), the most predominant form of Motor Neuron Disease (MND), is a progressive and fatal neurodegenerative condition that spreads throughout the neuromotor system by afflicting upper and lower motor neurons. Lower motor neurons project from the central nervous system and innervate muscle fibres at motor endplates, which degrade over the course of the disease leading to muscle weakness. The direction of neurodegeration from or to the point of neuromuscular junctions and the role of muscle itself in pathogenesis has continued to be a topic of debate in ALS research.MethodsTo assess the variation in gene expression between affected and nonaffected muscle tissue that might lead to this local degeneration of motor units, we generated RNA-seq skeletal muscle transcriptomes from 28 MND cases and 18 healthy controls and conducted differential expression analyses on gene-level counts, as well as an isoform switching analysis on isoform-level counts.ResultsWe identified 52 differentially-expressed genes (Benjamini-Hochberg-adjustedp< 0.05) within this comparison, including 38 protein coding, 9 long non-coding RNA, and 5 pseudogenes. Of protein-coding genes, 31 were upregulated in cases including with notable genes including the collagenicCOL25A1(p= 3.1 × 10−10),SAA1which is released in response to tissue injury (p= 3.6 × 10−5) as well as others of the SAA family, and the actin-encodingACTC1(p= 2.3 × 10−5). Additionally, we identified 17 genes which exhibited a functional isoform switch with likely functional consequences between cases and controls.ConclusionsOur analyses provide evidence of increased tissue generation in MND cases, which likely serve to compensate for the degeneration of motor units and skeletal muscle.

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value

The recent increase in obesity levels across many countries is likely to be driven by nongenetic factors. The epigenetic modification DNA methylation (DNAm) may help to explore this as it is sensitive to both genetic and environmental exposures. While the relationship between DNAm and body fat traits has been extensively studied [1–9], there is limited literature on the shared associations of DNAm variation across such traits. Akin to genetic correlation estimates, which measure the degree of common genetic control between two traits, here we introduce an approach to evaluate the similarities in DNAm associations between traits, DNAm correlations. As DNAm can be both a cause and consequence of complex traits [5, 10, 11], DNAm correlations have the potential to provide novel insights into trait relationships above that currently obtained from genetic and phenotypic correlations. Utilising 7,519 unrelated individuals from Generation Scotland (GS), we calculated DNAm correlations using the bivariate OREML framework in the OSCA software [12] to investigate the shared associations of DNAm variation between traits. For each trait we also estimated the shared contribution of DNAm between sexes. We identified strong, positive DNAm correlations between each of the body fat traits (BMI, body fat % and waist to hip ratio; ranging from 0.96 to 1.00), finding larger associations than those identified by genetic and phenotypic correlations. We identified a significant deviation from 1 in the rDNAmfor BMI between males and females, with sex-specific DNAm changes associated with BMI identified at eight DNAm probes. Employing genome-wide DNAm correlations to evaluate the similarities in the associations of DNAm with complex traits has provided novel insight into obesity related traits beyond that provided by genetic correlations.

Microarray technology has been used to measure genome-wide DNA methylation in thousands of individuals. These studies typically test the associations between individual DNA methylation sites (“probes”) and complex traits or diseases. The results can be used to generate methylation profile scores (MPS) to predict outcomes in independent data sets. Although there are many parallels between MPS and polygenic (risk) scores (PGS), there are key differences. Here, we review motivations, methods, and applications of DNA methylation-based trait prediction, with a focus on common diseases. We contrast MPS with PGS, highlighting where assumptions made in genetic modeling may not hold in epigenetic data.

Endometriosis is a leading cause of pain and infertility affecting millions of women globally. Identifying biologic and genetic effects on DNA methylation (DNAm) in endometrium increases understanding of mechanisms that influence gene regulation predisposing to endometriosis and offers an opportunity for novel therapeutic target discovery. Herein, we characterize variation in endometrial DNAm and its association with menstrual cycle phase, endometriosis, and genetic variants through analysis of genome-wide genotype data and methylation at 759,345 DNAm sites in endometrial samples from 984 deeply-phenotyped participants. We identify significant differences in DNAm profiles between menstrual cycle phases and at four DNAm sites between stage III/IV endometriosis and controls. We estimate that 15.4% of the variation in endometriosis is captured by DNAm, and identify DNAm networks associated with endometriosis. DNAm quantitative trait locus (mQTL) analysis identified 118,185 independentcis-mQTL including some tissue-specific effects. We find significant differences in DNAm profiles between endometriosis sub- phenotypes and a significant association between genetic regulation of methylation in endometrium and disease risk, providing functional evidence for genomic targets contributing to endometriosis risk and pathogenesis.

The complement system, including complement components 3 and 4 (C3, C4), traditionally has been linked to innate immunity. More recently, complement components have also been implicated in brain development and the risk of schizophrenia. Based on a large, population-based case-cohort study, we measured the blood concentrations of C3 and C4 in 68,768 neonates. We found a strong correlation between the concentrations of C3 and C4 (phenotypic correlation = 0.65,P-value < 1.0×10−100, genetic correlation = 0.38,P-value = 1.9×10−35). A genome-wide association study (GWAS) for C4 protein concentration identified 36 independent loci, 30 of which were in or near the major histocompatibility complex on chromosome 6 (which includes theC4gene), while six loci were found on six other chromosomes. A GWAS for C3 identified 15 independent loci, seven of which were located in theC3gene on chromosome 19, and eight loci on five other chromosomes. We found no association between (a) measured neonatal C3 and C4 concentrations, imputed C4 haplotypes, or predictedC4gene expression, with (b) schizophrenia (SCZ), bipolar disorder (BIP), depression (DEP), autism spectrum disorder, attention deficit hyperactivity disorder or anorexia nervosa diagnosed in later life. Mendelian randomisation (MR) suggested a small positive association between higher C4 protein concentration and an increased risk of SCZ, BIP, and DEP, but these findings did not persist in more stringent analyses. Evidence from MR supported causal relationships between C4 concentration and several autoimmune disorders: systemic lupus erythematosus (SLE, OR and 95% confidence interval, 0.37, 0.34 – 0.42); type-1 diabetes (T1D, 0.54, 0.50 - 0.58); multiple sclerosis (MS, 0.68, 0.63 - 0.74); rheumatoid arthritis (0.85, 0.80 - 0.91); and Crohn’s disease (1.26, 1.19 - 1.34). A phenome-wide association study (PheWAS) in UK Biobank confirmed that the genetic correlates of C4 concentration were associated a range of autoimmune disorders including coeliac disease, thyrotoxicosis, hypothyroidism, T1D, sarcoidosis, psoriasis, SLE and ankylosing spondylitis. We found no evidence of associations between C3 versus mental or autoimmune disorders based on either MR or PheWAS. In general, our results do not support the hypothesis that C4 is causally associated with the risk of SCZ (nor several other mental disorders). We provide new evidence to support the hypothesis that higher C4 concentration is associated with lower risks of autoimmune disorders.

BackgroundCpG methylation levels can help to explain inter-individual differences in phenotypic traits. Few studies have explored whether identifying CpG subsets based on biological and statistical properties can maximise predictions while minimising array content.MethodsVariance component analyses and penalised regression (epigenetic predictors) were used to test the influence of (i) the number of CpGs considered, (ii) mean CpG methylation variability and (iii) methylation QTL status on the variance captured in eighteen traits by blood DNA methylation. Training and test sets comprised ≤4,450 and ≤2,578 unrelated individuals from Generation Scotland, respectively.ResultsAs the number of CpG sites under consideration decreased, so too did the estimates from the variance components and prediction analyses. Methylation QTL status and mean CpG variability did not influence variance components. However, relative effect sizes were 15% larger for epigenetic predictors based on CpGs with methylation QTLs compared to sites without methylation QTLs. Relative effect sizes were 45% larger for predictors based on CpGs with mean beta-values between 10%-90% compared to those using hypo- or hypermethylated CpGs (beta-value ≤10% or ≥90%).ConclusionArrays with fewer CpGs could reduce costs, leading to increased sample sizes for analyses. Our results show that reducing array content can restrict prediction metrics and careful attention must be given to the biological and distribution properties of CpGs in array content selection.

Background Amyotrophic lateral sclerosis (ALS) is a complex, late-onset, neurodegenerative disease with a genetic contribution to disease liability. Genome-wide association studies (GWAS) have identified ten risk loci to date, including the TNIP1/GPX3 locus on chromosome five. Given association analysis data alone cannot determine the most plausible risk gene for this locus, we undertook a comprehensive suite of in silico, in vivo and in vitro studies to address this. Methods The Functional Mapping and Annotation (FUMA) pipeline and five tools (conditional and joint analysis (GCTA-COJO), Stratified Linkage Disequilibrium Score Regression (S-LDSC), Polygenic Priority Scoring (PoPS), Summary-based Mendelian Randomisation (SMR-HEIDI) and transcriptome-wide association study (TWAS) analyses) were used to perform bioinformatic integration of GWAS data (Ncases = 20,806, Ncontrols = 59,804) with ‘omics reference datasets including the blood (eQTLgen consortium N = 31,684) and brain (N = 2581). This was followed up by specific expression studies in ALS case-control cohorts (microarray Ntotal = 942, protein Ntotal = 300) and gene knockdown (KD) studies of human neuronal iPSC cells and zebrafish-morpholinos (MO). Results SMR analyses implicated both TNIP1 and GPX3 (p < 1.15 × 10−6), but there was no simple SNP/expression relationship. Integrating multiple datasets using PoPS supported GPX3 but not TNIP1. In vivo expression analyses from blood in ALS cases identified that lower GPX3 expression correlated with a more progressed disease (ALS functional rating score, p = 5.5 × 10−3, adjusted R2 = 0.042, Beffect = 27.4 ± 13.3 ng/ml/ALSFRS unit) with microarray and protein data suggesting lower expression with risk allele (recessive model p = 0.06, p = 0.02 respectively). Validation in vivo indicated gpx3 KD caused significant motor deficits in zebrafish-MO (mean difference vs. control ± 95% CI, vs. control, swim distance = 112 ± 28 mm, time = 1.29 ± 0.59 s, speed = 32.0 ± 2.53 mm/s, respectively, p for all gpx3 expression, with no phenotype identified with tnip1 KD or gpx3 overexpression. Conclusions These results support GPX3 as a lead ALS risk gene in this locus, with more data needed to confirm/reject a role for TNIP1. This has implications for understanding disease mechanisms (GPX3 acts in the same pathway as SOD1, a well-established ALS-associated gene) and identifying new therapeutic approaches. Few previous examples of in-depth investigations of risk loci in ALS exist and a similar approach could be applied to investigate future expected GWAS findings.

Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNA methylation (DNAm) signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample (Generation Scotland; n = 9537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore-disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors, and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification.Although our genetic code does not change throughout our lives, our genes can be turned on and off as a result of epigenetics. Epigenetics can track how the environment and even certain behaviors add or remove small chemical markers to the DNA that makes up the genome. The type and location of these markers may affect whether genes are active or silent, this is, whether the protein coded for by that gene is being produced or not. One common epigenetic marker is known as DNA methylation. DNA methylation has been linked to the levels of a range of proteins in our cells and the risk people have of developing chronic diseases. Blood samples can be used to determine the epigenetic markers a person has on their genome and to study the abundance of many proteins. Gadd, Hillary, McCartney, Zaghlool et al. studied the relationships between DNA methylation and the abundance of 953 different proteins in blood samples from individuals in the German KORA cohort and the Scottish Lothian Birth Cohort 1936. They then used machine learning to analyze the relationship between epigenetic markers found in people’s blood and the abundance of proteins, obtaining epigenetic scores or ‘EpiScores’ for each protein. They found 109 proteins for which DNA methylation patterns explained between at least 1% and up to 58% of the variation in protein levels. Integrating the ‘EpiScores’ with 14 years of medical records for more than 9000 individuals from the Generation Scotland study revealed 137 connections between EpiScores for proteins and a future diagnosis of common adverse health outcomes. These included diabetes, stroke, depression, Alzheimer’s dementia, various cancers, and inflammatory conditions such as rheumatoid arthritis and inflammatory bowel disease. Age-related chronic diseases are a growing issue worldwide and place pressure on healthcare systems. They also severely reduce quality of life for individuals over many years. This work shows how epigenetic scores based on protein levels in the blood could predict a person’s risk of several of these diseases. In the case of type 2 diabetes, the EpiScore results replicated previous research linking protein levels in the blood to future diagnosis of diabetes. Protein EpiScores could therefore allow researchers to identify people with the highest risk of disease, making it possible to intervene early and prevent these people from developing chronic conditions as they age.

Ankylosing Spondylitis (AS) is a highly heritable inflammatory arthritis which occurs more frequently in men than women. In their recent publication examining sex differences in the genetic aetiology of common complex traits and diseases, Bernabeu et al. (2021) observe differences in heritability of AS between sexes, and a genome-wide significant genotype by sex interaction in risk of AS at the major histocompatability (MHC) locus1. The authors then present evidence suggesting that this genotype by sex interaction arises primarily as a result of differential expression of the gene MICA across the sexes in skeletal muscle tissue. Through a series of conditional association analyses in the UK Biobank, reanalysis of the GTEx gene expression resource and RNASeq experiments on peripheral blood cells from AS cases and controls, we show that the genotype by sex interaction the authors’ report is unlikely to be a result of variation in MICA, but probably reflects a known interaction between the HLA-B gene, sex and risk of AS. We demonstrate that the diagnostic accuracy of AS in the UK Biobank is low, particularly amongst women, likely explaining some of the observed differences in heritability across the sexes and the difficulty in precisely locating association signals in the cohort.

The mechanisms underlying how monozygotic (or identical) twins arise are yet to be determined. Here, the authors investigate this in an epigenome-wide association study, showing that monozygotic twinning has a characteristic DNA methylation signature in adult somatic tissues. Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin., Netherlands Organization for Scientific Research (NWO); Biobanking and Biomolecular Research Infrastructure (BBMRI–NL); NWO Large Scale infrastructures; X-Omics

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